Pentaquark and Tetraquark states
Yan-Rui Liu, Hua-Xing Chen, Wei Chen, Xiang Liu, Shi-Lin Zhu

TL;DR
This review updates recent experimental and theoretical progress on multiquark states, especially hidden heavy flavor systems, highlighting new discoveries, theoretical models, and lattice QCD insights into exotic tetraquark and pentaquark states.
Contribution
It provides a comprehensive update on recent experimental findings, theoretical models, and lattice QCD simulations of multiquark states, emphasizing the role of QCD-based predictions and upcoming searches.
Findings
Confirmation of hidden-charm pentaquarks by LHCb
Evidence supporting molecular interpretations of multiquark states
Predictions of fully heavy tetraquarks and their production at colliders
Abstract
The past seventeen years have witnessed tremendous progress on the experimental and theoretical explorations of the multiquark states. The hidden-charm and hidden-bottom multiquark systems were reviewed extensively in [Phys. Rept. 639 (2016) 1-121]. In this article, we shall update the experimental and theoretical efforts on the hidden heavy flavor multiquark systems in the past three years. Especially the LHCb collaboration not only confirmed the existence of the hidden-charm pentaquarks but also provided strong evidence of the molecular picture. Besides the well-known and states, we shall discuss more interesting tetraquark and pentaquark systems either with one, two, three or even four heavy quarks. Some very intriguing states include the fully heavy exotic tetraquark states and doubly heavy tetraquark states , where is a heavy quark.…
| Hadron | CMI | Th. | Ex. | (Th.-Ex.) | Hadron | CMI | Th. | Ex. | (Th.-Ex.) |
|---|---|---|---|---|---|---|---|---|---|
| 246.6 | 139.6 | 107 | 882.3 | 775.3 | 107 | ||||
| 602.8 | 493.7 | 109 | 1001.7 | 891.8 | 110 | ||||
| 882.3 | 782.7 | 100 | |||||||
| 1136.7 | 1019.5 | 117 | |||||||
| 1980.7 | 1869.7 | 111 | 2121.5 | 2010.3 | 111 | ||||
| 2157.7 | 1968.3 | 189 | 2300.6 | 2112.2 | 188 | ||||
| 5380.9 | 5279.5 | 102 | 5425.7 | 5324.7 | 101 | ||||
| 5556.3 | 5366.9 | 189 | 5605.4 | 5415.4 | 190 | ||||
| 3364.4 | 2983.9 | 381 | 3477.5 | 3096.9 | 381 | ||||
| 10059.2 | 9399.0 | 660 | 10121.1 | 9460.3 | 661 | ||||
| 6724.6 | 6274.9 | 450 | [33] | 6795.0 | |||||
| 937.9 | 938.3 | 0 | 1232.3 | 1232.0 | 0 | ||||
| 1183.7 | 1189.4 | -6 | 1377.3 | 1382.8 | -6 | ||||
| 1116.5 | 1115.7 | 1 | |||||||
| 1330.6 | 1314.9 | 16 | 1524.2 | 1531.8 | -8 | ||||
| 1672.9 | 1672.5 | 0 | |||||||
| 2300.8 | 2286.5 | 14 | |||||||
| 2454.4 | 2454.0 | 0 | 2518.4 | 2518.4 | 0 | ||||
| 2529.8 | 2467.9 | 62 | |||||||
| 2613.5 | 2577.4 | 36 | 2681.5 | 2645.5 | 36 | ||||
| 2774.5 | 2695.2 | 79 | 2846.5 | 2765.9 | 81 | ||||
| 5629.0 | 5619.6 | 9 | |||||||
| 5811.4 | 5811.3 | 0 | 5832.2 | 5832.1 | 0 | ||||
| 5858.0 | 5794.5 | 64 | |||||||
| 5973.7 | 5935.0 | 39 | 5993.7 | 5955.3 | 39 | ||||
| 6137.9 | 6046.1 | 92 | 6157.1 | ||||||
| 3777.0 | 3621.4 | 156 | 3841.0 | ||||||
| 3950.2 | 4022.2 | ||||||||
| 10458.2 | 10479.0 | ||||||||
| 10637.9 | 10657.1 | ||||||||
| 7106.6 | |||||||||
| 7132.4 | 7158.7 | ||||||||
| 7281.2 | |||||||||
| 7312.9 | 7338.3 | ||||||||
| 5200.0 | |||||||||
| 15172.7 | |||||||||
| 8489.0 | 8521.6 | ||||||||
| 11813.2 | 11845.8 |
| State | Upper limit | Estimation |
|---|---|---|
| 4977 | ||
| 14937 | ||
| 11610 | ||
| 8285 | ||
| 11577 | ||
| 8253 | ||
| 3669 | ||
| 10347 | ||
| 3799 | ||
| 10474 | ||
| 7048 | ||
| 7174 | ||
| 3733 | ||
| 10368 | ||
| 3871 | ||
| 10493 | ||
| 6996 | ||
| 7022 | ||
| 7117 | ||
| 7149 |
| Tetraquarks | Pentaquarks |
|---|---|
| , | |
| , | , |
| , | |
| , | |
| Current | Structure | [GeV2] | Borel Mass [GeV2] | Mass [GeV] | (, ) |
|---|---|---|---|---|---|
| () | |||||
| () | |||||
| () | |||||
| () | |||||
| () | |||||
| () |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Pentaquark and Tetraquark states
Yan-Rui Liu,1111These authors equally contribute to this work. Hua-Xing Chen,211footnotemark: 1 Wei Chen,311footnotemark: 1 Xiang Liu,4,5222Corresponding author: [email protected] Shi-Lin Zhu6,7,8333Corresponding author: [email protected]
1School of Physics, Shandong University, Jinan 250100, China
2School of Physics, Beihang University, Beijing 100191, China
3School of Physics, Sun Yat-Sen University, Guangzhou 510275, China
4School of Physical Science and Technology,
Lanzhou University, Lanzhou 730000, China
5Research Center for Hadron and CSR Physics, Lanzhou University and
Institute of Modern Physics of CAS, Lanzhou 730000, China
6School of Physics and State Key Laboratory of Nuclear Physics and Technology,
Peking University, Beijing 100871, China
7Collaborative Innovation Center of Quantum Matter, Beijing 100871, China
8Center of High Energy Physics, Peking University, Beijing 100871, China
Abstract
The past seventeen years have witnessed tremendous progress on the experimental and theoretical explorations of the multiquark states. The hidden-charm and hidden-bottom multiquark systems were reviewed extensively in Ref. [1]. In this article, we shall update the experimental and theoretical efforts on the hidden heavy flavor multiquark systems in the past three years. Especially the LHCb collaboration not only confirmed the existence of the hidden-charm pentaquarks but also provided strong evidence of the molecular picture. Besides the well-known and states, we shall discuss more interesting tetraquark and pentaquark systems either with one, two, three or even four heavy quarks. Some very intriguing states include the fully heavy exotic tetraquark states and doubly heavy tetraquark states , where is a heavy quark. The states may be produced at LHC while the system may be searched for at BelleII and LHCb. Moreover, we shall pay special attention to various theoretical schemes such as the chromomagnetic interaction (CMI), constituent quark model, meson exchange model, heavy quark and heavy diquark symmetry, QCD sum rules, Faddeev equation for the three body systems, Skyrme model and the chiral quark-soliton model, and the lattice QCD simulations. We shall emphasize the model-independent predictions of various models which are truly/closely related to Quantum Chromodynamics (QCD). For example, the chromomagnetic interaction arises from the gluon exchange which is fundamental and universal in QCD and responsible for the mixing and mass splitting of the conventional mesons and baryons within the same multiplet. The same CMI mechanism shall also be responsible for the mixing of the different color configurations and mass splittings within the multiplets in the multiquark sector. There have also accumulated many lattice QCD simulations through multiple channel scattering on the lattice in recent years, which provide deep insights into the underlying structure/dynamics of the states. In terms of the recent states, the lattice simulations of the charmed baryon and anti-charmed meson scattering are badly needed. We shall also discuss some important states which may be searched for at BESIII, BelleII and LHCb in the coming years.
Contents
1 Introduction
In the past decades, hadron physicists show great interest in hunting for evidences of the multiquark states. As a new form of matter beyond conventional mesons () and baryons (), multiquark states containing more than three quarks are of special importance in the hadron family. Especially, with the observations of various charmonium-like states, the study of multiquark states has entered upon a new era.
In this article, we will give a concise review on the research progress of the tetraquark () and pentaquark () states, which are typical multiquark matters. Before doing that, let us first look back on the history of particle/hadron physics. It is well known that the development of quantum mechanics is closely related to the study of atomic spectroscopy, which reveals the mysterious veil of the atom’s microstructure. Here, we find a similar situation that the study of hadron spectroscopy is bringing us new insights into the internal structure and dynamics of hadrons.
In the early 1960s, many strongly interacting particles were observed in particle/nucleon experiments, which were named as “hadron” by L. B. Okun later [2]. Based on these observations, M. Gell-Mann and G. Zweig independently developed the classification scheme for hadrons—the quark model [3, 4, 5]. Especially, the name “quark” was given by M. Gell-Mann to denote the substructure of hadrons. The quark model achieves a great success, and it is a milestone in the development of particle physics. A well-known example is the prediction of the baryon with three quarks, that is the . It was discovered in 1964 [6] after its existence, mass, and decay products had been predicted in 1961 independently by M. Gell-Mann [7] and Y. Ne’eman [8].
Then we would like to mention the establishment of the Cornell model. In 1974, as the first state in the charmonium family, the was discovered by the E598 Collaboration [9] and the SLAC-SP-017 Collaboration [10], which confirmed the existence of the quark [11]. After that, a series of charmonia were discovered, such as the [12], [13], [14], [15], and [16], etc. Based on these experimental observations, the Cornell model was proposed [17, 18], which uses the Cornell potential [19] to depict the interaction between the and - quarks. This potential consists of both the Coulomb-type and linear potentials, which can well reproduce masses of the above observed charmonia at that time. Inspired by the Cornell model, various potential models were developed [20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33], among which the Godfrey-Isgur model [33] is quite popular. This model contains the semi-relativistic kinetic and potential energy terms, which can be applied to quantitatively describe not only meson spectra [33] but also baryon spectra [34]. These phenomenological models, inspired by the charmonium family, greatly improve our understanding of the internal structure and dynamics of conventional mesons and baryons.
Nowadays several hundreds of hadrons have been observed in particle experiments, most of which can be categorized into two families: baryons made of three quarks and mesons made of one quark and one antiquark [35]. They are formed by the , , , , and quarks/antiquarks (the quark has a very short lifetime and therefore does not form hadrons), which are governed by the strong interaction. Although the mechanism of the color confinement remains one of the most difficult problems in particle physics, various phenomenological models were proposed to quantitatively describe the hadron spectroscopy. Due to the joint efforts from the particle/hadron theorists and experimentalists, the meson and baryon families have been nearly complete, and almost all the ground-state hadrons have been discovered. Especially, the doubly charmed baryon was recently discovered by LHCb Collaboration [36].
Other than hundreds of conventional mesons and baryons, there were only tens of multiquark candidates observed in particle experiments, although the concept of the multiquark states appeared together with the birth of the quark model. For example, M. Gell-Mann indicated [3]: ”Baryons can now be constructed from quarks by using the combinations , , etc., while mesons are made out of , , etc”, and G. Zweig [5] also wrote: ”In general, we would expect that baryons are built not only from the product of these aces, , but also from , , etc., where denotes an anti-ace. Similarly, mesons could be formed from , , etc”. Due to the constraint of reliable experimental and theoretical techniques, the study of the multiquark states was far from satisfactory before 2003. At that time, theorists mainly focused on a) the light scalar mesons , , , and , b) the scalar mesons , , and , and c) the , and discussed their possible interpretations as multiquark states and glueballs (see review articles [37, 38] for details).
2003 is an important year in the history of the multiquark states. Since 2003, there has been continuous progress in this field. With the accumulation of experimental data, a series of charmonium-like states were reported (see Fig. 1), which stimulated us to reveal their exotic inner structures. We shall try our best to convey the progress and excitement to the readers in the present review. Here are several typical examples:
- •
The BES and BESIII Collaborations reported a series of light-flavor multiquark candidates, including the observed in the decay [39], the in [40, 41], the in [42], the in [43, 44], and so on. These observations inspired extensive discussions, and various exotic interpretations such as multiquark states were proposed. Since the present review mainly focuses on the heavy-flavor multiquark states, interested readers may consult Ref. [38] for more information.
- •
In 1997, Diakonov et al. predicted the existence of a light-flavor pentaquark with the chiral soliton model [45], which has a mass around 1530 MeV, spin , isospin [math] and strangeness . In 2003, LEPS Collaboration claimed that a narrow particle consistent with the above prediction was observed in the reaction [46]. Very quickly, the , as a pentaquark candidate, became a super star at that time [47]. Unfortunately, it was not confirmed by the subsequent series of high precise experiments [48]. The is probably not a genuine resonance [49]. The rise and fall of the unveiled our ignorance on the non-perturbative behavior of quantum chromodynamics. The lessons and experience with the have been helpful in the search of the pentaquark and tetraquark states.
- •
In 2003, BaBar Collaboration observed a narrow heavy-light state in the invariant mass spectrum [50]. Since the mass of the is about MeV lower than the quark model prediction of the -wave charmed-strange meson with [33], its tetraquark explanation was proposed in Refs. [51, 52, 53]. A similar situation happened to the observed by CLEO Collaboration [54]. The strong channel coupling between the S-wave DK() scattering states and bare quark model states plays a pivotal role in lowering the masses of the and [55, 56].
- •
Still in 2003, , as the first particle in charmonium-like family, was discovered by Belle Collaboration [57]. Since the mass of the is close to the threshold of the pair, its assignment as the molecular state is very popular. We note that its nature is still under heated debates today. Similar to the and , there also exists a low mass puzzle related to the . Its mass is much lower than the quark model prediction of the charmonium state [33]. Again, the coupled channel effect may mediate this difference [58].
In the past 17 years, more and more candidates of the exotic hadrons were observed in the Belle, BaBar, BESIII, D0, CDF, CMS, and LHCb experiments, such as a) the charmonium-like states , , and , b) the bottomonium-like states and , and c) the hidden-charm pentaquark states and , etc. In Figs. 2-3, we summarize these observations concisely. An extensive review of the states observed before 2016 can be found in Ref. [1]. Here we briefly introduce the experimental observations after 2016:
- •
In the invariant mass spectrum of the decay, LHCb established the existence of the structure with significance [59, 60]. This state has a mass MeV, width MeV, and quantum numbers . Before this observation, CDF and CMS had reported an evidence of a structure around 4274 MeV in the invariant mass spectrum [61, 62]. Besides the structure, two extra structures with higher masses were found in the invariant mass spectrum, which are the and with . Their resonance parameters in units of MeV are [59, 60]
[TABLE]
- •
In 2017, BESIII Collaboration presented more precise measurements of the cross section at the center-of-mass energy from 3.77 to 4.60 GeV [63] and the cross section at the center-of-mass energy from 3.896 to 4.600 GeV [64]. Three vector structures , and were reported with the resonance parameters in units of MeV:
[TABLE]
In fact this updated analysis of the cross section [64] shows that the contains two substructures, the and .
- •
The Belle Collaboration carried out an analysis of the cross section, and concluded that a broad structure exists in the invariant mass spectrum, which has a mass MeV and width MeV [65]. Their result shows that is favored over at the level of [65].
- •
A new charged structure was reported by BESIII Collaboration in the invariant mass spectrum of the process [66]. BESIII confirmed the vector charmonium-like structure in the invariant mass spectrum.
- •
In 2018, LHCb found the evidence of a structure in the invariant mass spectrum of the decay, which was named as the with a mass MeV and width MeV [67]. The LHCb measurement also indicated that the structure may be described under the spin-parity or assignments [67].
- •
Very recently, the LHCb Collaboration announced the observation of three new pentaquarks [68] at the Rencontres de Moriond QCD conference, as shown in Fig. 4. In the measured invariant mass spectrum, a new pentaquark was discovered with a significance. The new LHCb analysis further found that the is composed of two substructures and with 5.4 significance. Their resonance parameters are collected as following
[TABLE]
The isospin of , and is since these three pentaquarks were discovered in the channel. The lies below the threshold, while the masses of the and are slightly lower than the threshold. The observation of , and clearly confirms the hidden-charm molecular pentaquarks.
The present situation of the charmonium-like and states is a bit similar to a) that of the ground-state mesons and baryons in 1960s, and b) that of the charmonia in 1980s. With such abundant novel phenomenon, the most important task is to identify the genuine tetraquark and pentaquark signals. There have accumulated hundreds of investigations of the tetraquark and pentaquark systems with various phenomenological models. In this review, we will introduce several typical phenomenological methods/models and their applications to multiquark states.
Our previous review focused on the hidden-charm and hidden-bottom multiquark systems [1]. In the present article, we will summarize the experimental and theoretical progress on the hidden-charm tetraquark and pentaquark states in the past three years. Besides the hidden heavy flavor systems, we shall discuss more interesting tetraquark and pentaquark systems either with open heavy flavor or with three/four heavy quarks. These new states include the exotic tetraquarks and systems where is a heavy quark. The states may be produced at LHC while the system may be searched for at BelleII. Moreover, we will summarize theoretical predictions of the multiquark states from various formalisms such as the chromomagnetic interaction (CMI), constituent quark model, meson exchange model, heavy quark and heavy diquark symmetry, QCD sum rules, Faddeev equation for the three body systems, Skyrme model and the chiral quark-soliton model, and the lattice QCD simulations. We shall pay special attention to those theoretical schemes which are not covered or not addressed in depth in our previous review [1], such as the chromomagnetic interaction, heavy quark symmetry etc. We shall emphasize the model-independent predictions of these models which are truly/closely related to Quantum Chromodynamics (QCD).
Besides the present article and our previous review [1], there exist many nice reviews of the XYZ states in literature [69, 70, 71, 72, 73, 74, 75, 76]. Interested readers are encouraged to consult them for a glimpse of this extremely active and vast field.
2 symmetry and Chromomagnetic interaction (CMI)
The hyperfine structure for atoms is induced by the spin-related interaction between electrons and nuclei. In a similar situation, the hyperfine structure in hadron spectroscopy is from the spin-related interaction between quarks or between quarks and antiquarks. In the Hamiltonian or Lagrangian formalism, such a term contains a color factor since quarks interact in the color space. The simple color-magnetic interaction arises from the one-gluon-exchange potential and causes the mass splittings of the conventional hadrons, whose color configuration is unique. After the inclusion of the quark mass, one obtains an effective description method for hadron masses, which is the Hamiltonian of the color-magnetic interaction or chromomagnetic interaction (CMI) model.
Although the CMI model is simple and the results may sometimes deviate from reality, it is helpful to understand the basic features of hadron spectra through the estimated masses, see Ref. [77] for detailed arguments. This type of simple models together with symmetry analysis can also be helpful to confirm dynamical models to some extent. For example, from the analysis of the mass splittings of the conventional mesons and baryons, the author of Ref. [78] concluded that the Cornell potential (Coulomb+linear) is the preferred model.
Within the multiquark systems, there may exist two or more color configurations. For example, one may decompose the fully heavy tetraquark into two clusters. There are two possible color configurations: and . The type color-electric interactions do not induce mixing between different color configurations. In contrast, the chromomagnetic interaction will not only mix two different color configurations, but also cause the mass splitting within the multiplets.
Recent developments in spectroscopy motivated the proposal of improved CMI models or methods which are becoming powerful tools in studying interesting multiquark states. We shall discuss recent developments about the model and application methods in this section.
2.1 Hamiltonian for the CMI model
A realistic quark potential model generally includes the quark kinetic term, color Coulomb term, color-spin interaction term, spin-orbital term, tensor term, and color confinement term. Considering the Hamiltonian in Ref. [79] for the -wave hadrons (where the sin-orbit and tensor interactions vanish) and ignoring the electromagnetic part, one has
[TABLE]
Here, the Gell-Mann matrix should be replaced by for an antiquark, is responsible for the quark binding, , , and are the position, momentum, and constituent mass of the -th (anti)quark, respectively, and . The average on the orbital wave function () gives
[TABLE]
where containing various effects is the effective mass for the th (anti)quark and . The coupling parameter is also written as in the literature. Note that or for mesons differs from that for baryons. This mass expression is also called Sakharov-Zeldovich formula [80, 78].
Now, the model contains only the color-spin operator and this Hamiltonian should have an symmetry. The eigenstate of is also that of . One may express the eigenvalue on a state with Casimir operators. For a state containing only identical quarks whose flavor-spin-color wave function is , one has [81, 82]
[TABLE]
where is the number of quarks, is the total spin, and the definitions for the quadratic Casimir operators are
[TABLE]
If a representation is specified by the Young diagram , one has
[TABLE]
For a quark-antiquark pair or two quarks without constraint from the Pauli principle, the color part and spin part can be calculated separately and the formula in this special case is
[TABLE]
If and one acts on a multiquark system whose quark content is (), the average value can be expressed as [83, 84]
[TABLE]
Now, is the total number of quarks and antiquarks. One may also consult Ref. [85] for more information.
In principle, the symmetry breaking should be considered. Then the coupling parameters can be different for different flavors and more complicated CMI expressions may be obtained (see Eq. (5) of Ref. [84]). In addition, mixing between different color-spin structures may be induced by the chromomagnetic interaction. The number of base states for a given is generally larger than one and we need to know all the relevant CMI expressions which can be calculated with the knowledge of group theory [86, 87, 88]. After diagonalizing the obtained matrices ’s, one obtains eigenvalues of . Alternatively, one may directly calculate the matrix elements with the constructed wave functions in color and spin spaces separately and then combine the results together [89]. If the results in the two spaces have the form
[TABLE]
where or is the label of a base state, one gets (See also Ref. [77])
[TABLE]
2.2 CMI models
In principle, the values of and in the CMI model (2.1) should be different for various systems, which can be understood from the reduction procedure in getting the Hamiltonian. However, it is difficult to exactly calculate these parameters for a given system without knowing the spatial wave function. Practically, they are extracted from the masses of conventional hadrons by assuming that quark-(anti)quark interactions are the same for all the hadron systems. This assumption certainly leads to uncertainties on mass estimations. The uncertainty caused by does not allow us to give accurate multiquark masses while the uncertainty in coupling parameters has smaller effects and the mass splittings should be more reliable. To give reasonable description for hadron spectra in the CMI model, one usually uses modified mass equation based on the adopted assumptions or chooses refined parameters in calculations. In the literature, different forms of CMI Hamiltonians can be found. Here we summarize several versions of the model in studying heavy quark multiquark states.
- •
Models with diquark assumption. The -wave correlated quark-quark state (diquark) with color= is assumed [90, 91], where the spin of the diquark can be 1 and 0 (see Ref. [92] for a review). For the light-light diquark, the spin-0 state is more tightly bound than the spin-1 state. For the heavy-light diquark, the difference between the spin-1 and spin-0 states is suppressed by the heavy quark mass. The formation mechanism of tetraquark states in the diquark configuration was pictured in Ref. [93] and the hypothesis that diquarks and antidiquarks in tetraquarks are separated by a potential barrier was introduced in Ref. [94]. In this diquark approximation, the Hamiltonian for tetraquark states has the form [95, 96]
[TABLE]
or the form [97]
[TABLE]
The former (latter) Hamiltonian is for the “type-I” (“type-II”) diquark-antidiquark model where orbital excitation between the diquark and the antidiquark is involved. In the type-II model, the spin-spin interactions between quark components in different diquarks are ignored and the coupling coefficients have different values from those in the type-I model.
- •
Model with triquark assumption. To understand the inner structure of the LHCb pentaquark states, Lebed proposed the diquark-triquark picture in Ref. [98] where the triquark is composed of a diquark and an antiquark. The diquarks from their generation are always in the color- representation [93]. This picture is different from the one proposed in Ref. [99] where the diquark in the triquark is in the color- representation. Another difference lies in the stability. The diquark-triquark state in Ref. [99] looks like a static molecule stabilized by a -wave barrier, while the state in Ref. [98] can exist in each partial wave and can last as long as the components continue to separate. Compared with the diquark-diquark-antiquark configuration, the antiquark belongs to a compact component of the pentaquark, although the two configurations have the same color structure. This diquark-triquark picture allows the assignment consistent with LHCb data and can qualitatively explain the measured widths. The effective Hamiltonian in this diquark-triquark picture can be written in a form similar to the model (13) [100],
[TABLE]
where () indicates the diquark (triquark) and the subscripts of mean the spin-spin, spin-orbital, and orbital-orbital interactions within or between different clusters.
- •
Model with the chromoelectric term. In Ref. [101], Hogaasen et al generalized the CMI model by including the chromoelectric term,
[TABLE]
where the effective quark masses have different values from those in Eq. (2.1) because of the non-vanishing . The contributions from the chromoelectric term were implicitly included in the model (2.1).
Ref. [102] adopted a modified form of (16) to discuss the masses of the doubly and triply heavy conventional baryons. The model Hamiltonian reads
[TABLE]
where , , and is the to-be-fitted mass parameter of the quark pair labeled with and .
- •
Model motivated with a QCD-string-junction picture [103]. Noting the fact that the effective quark masses and coupling parameters determined from mesons differ from those from baryons, Karliner, Rosner, and Nussinov proposed a method to estimate hadron masses in Refs. [104, 105]. Since there is no junction (but one QCD string) for conventional mesons and one junction (and two additional QCD strings) for conventional baryons, a universal constant is added to the baryon mass expression. For tetraquark states, the number of added to mass expressions is 2 because there are two junctions (and two more QCD strings). Besides, a correction term representing additional binding is usually needed. Therefore, the mass formula in this model generally has the form
[TABLE]
where is the effective quark mass term, is an integer, and is the hyperfine splitting term.
- •
Method with a reference mass scale. From the explicit calculations, one finds that the model (2.1) usually gives higher masses for conventional hadrons (see Table 2 for a comparison). According to our studies for various multiquark systems [106, 107, 108, 109, 110, 111, 112, 113], the obtained masses with this model are also the largest values we can obtain. They should be higher than the realistic masses. Here, “realistic” means “measured” by assuming that the states do exist. These observations indicate that the attractions between quark components in this simple model are not sufficiently considered. In order to reduce the uncertainties and obtain more appropriate estimations, one may adopt an alternative form of the mass formula.
Introducing a reference system and modifying the mass formula, one gets
[TABLE]
where and are the physical mass of the reference system and the corresponding CMI eigenvalue, respectively. For , one may use the mass of a reference multiquark state or use the threshold of a reference hadron-hadron system whose quark content is the same as the considered multiquark states. With this method, the problem of using extracted quark masses from conventional hadrons in multiquark systems [114] is evaded and part of missed attraction between quark components is phenomenologically compensated. Whether this manipulation gives results close to realistic masses or not can be tested in more studies. In fact, such an estimation method had been applied to the multiquark states many years ago [77, 115, 116].
In calculating the CMI matrix elements for the multiquark states, we consider all the color-spin configurations and the important color mixing effects [117]. When one adopts this method, more than one thresholds or scales may exist. Since the size of a multiquark state is expected to be larger than that of a conventional hadron and the distance between quark components may be larger, the attraction between quark components should be weaker. The highest masses seem to be more reasonable in this method and we use this assumption in our studies, although we cannot give a definite answer for the problem which threshold should be adopted. From the theoretical studies in recent years, it seems that the predicted masses with such a strategy are still lower than the realistic values. That is,
[TABLE]
where and are obtained with Eqs. (19) and (2.1), respectively.
2.3 CMI and Conventional hadrons
The parameters in CMI models are effective masses and coupling coefficients. In the simplest CMI model (2.1), one may extract effective quark masses and most of the coupling parameters with the known hadrons. There are different approaches to determine them, by fitting the masses of all known hadrons, by fitting the masses of mesons and baryons differently, or by just extracting their values with parts of hadron masses. A set of the extracted effective quark masses is MeV (), MeV, MeV, and MeV while the obtained coupling parameters are listed in Table 1. For the unknown , its value was extracted with the meson masses predicted in the Godfrey-Isgur (GI) model [33]. For the unknown , , , and , the approximation [78, 118] has been used. This approximation seems to be applicable for all the coupling parameters. The parameters presented here are consistent with those in the literature, e.g. [95, 104].
When substituting these parameters into Eq. (2.1) for various conventional hadrons and comparing the obtained masses with experiments (Table 2), it is obvious that most theoretical masses are higher than the measured masses. This observation illustrates that the effective attraction is not appropriately considered in the model (2.1). To reproduce the masses of the conventional hadrons in acceptable uncertainties in the CMI model, it is necessary to adopt modified Hamiltonian by including more effects. In the mentioned CMI models, only two models, (16) or (17) and (18), are appropriate to achieve the task. Recently, with the model (18), Karliner and Rosner obtained reasonable masses for conventional mesons and baryons in Ref. [104]. The predicted mass of , MeV, is very close to that measured by LHCb, MeV [36]. With the extended CMI model whose Hamiltonian is given in Eq. (17), Weng, Chen, and Deng explored in Ref. [102] the masses of doubly and triply heavy baryons. The parameters were obtained by fitting masses of known hadrons and they can reproduce the baryon masses well. The calculated mass ( MeV) is also close to that measured by LHCb. For other predicted masses, the values in these two methods are also consistent with each other. All the predicted and masses are not far from the lattice results [119, 120], either. The studies indicate that the CMI model is still powerful once chromoelectric interactions are appropriately considered, although the model is very simple.
2.4 CMI and
Most of the observed states belong to this category. Decades ago, such type of tetraquark states had been discussed after the excited charmonium state was observed. Since the observation of the in 2003, more and more heavy hadrons with exotic properties were observed and partly confirmed in the following years. If the states have explicit exotic quantum numbers, they are good tetraquark candidates. Otherwise, some of these states may be conventional heavy quarkonia which are affected heavily by the mixing effects with four-quark components. It’s possible some states are molecules or compact tetraquarks with non-exotic quantum numbers. In some cases, there exist non-resonant explanations. We are still far from understanding the nature of the states. Here, we first concentrate on recent theoretical studies and developments within the CMI models, especially in recent three years. One may consult studies about the states before 2016 in our previous review [1].
In Ref. [95], Maiani et al discussed the nature of in the “type-I” diquark-antidiquark model (13). As the tetraquark state, the properties of such as its narrow width and isospin breaking in decay can be understood. Within this model, they also discussed the properties of [now called ], , , and . Latter in Ref. [121], they interpreted as the first orbital excitation of a diquark-antidiquark state (). In Refs. [122, 123], the was identified as the first radial excitation of the tetraquark basic supermultiplet to which belongs. In Ref. [96], the authors studied the states with different and predicted a state at 4277 MeV decaying into .
Hogaasen, Richard, and Sorba discussed the nature of in the tetraquark configuration with the mass equation (2.1) in Ref. [124]. By introducing the annihilation term and different masses for the and quarks, they found that the may be interpreted as a tetraquark state where the pair is mostly a color-octet state. This configuration can also explain the rough ratio between the and decay modes. However, their parameters could not be applied to the hidden-bottom tetraquark states.
Cui et al. considered the partner states of the with the mass equation (2.1) [88]. They determined the coupling parameter by assuming the as a tetraquark. With the obtained parameters, the masses of all the ground states with the configurations , , and were estimated, where and .
Stancu explored the spectrum of tetraquarks with the CMI model (2.1) after the observation of the [114]. The results suggest that the could be the strange partner of with in a tetraquark interpretation. The later LHCb measurement for the [60] indicates that the mass splitting between the and is consistent with the CMI model prediction.
Tetraquark states were also studied with the CMI model (2.1) together with the possible hexaquark states [125]. The and can be assigned as and states, respectively.
In a schematic study of the S-wave () tetraquarks with the CMI model (2.1), the authors of Ref. [126] obtained the or states with masses MeV and MeV, which are compatible with the and .
In Refs. [127, 128], Ali et al considered the tetraquark interpretation for the bottomonium-like states and with . They assumed that these mesons are diquark-antidiquark states produced in the decay of the which was also assumed to be a tetraquark state [129, 130]. Using the Hamiltonian (13) and including the meson loop effects, they were able to reproduce their measured masses [35].
After the observation of the , Maiani et al interpreted it as the tetraquark state around 3880 MeV [131], which was predicted in their previous study within the diquark-antidiquark model (13) [95]. They also discussed the possibility of another peak around 3760 MeV.
In Ref. [101], Hogaasen et al used the Hamiltonian (16) to discuss possible partner states of the . Their fitted parameters can reproduce the mass spectrum of the states obtained with Eq. (2.1) in Ref. [124]. The states ( and ) may lie around 10.62 GeV, and the or is such a tetraquark candidate, whose component is almost a color-octet state. They also identified the implied in Ref. [124] as the corresponding tetraquark.
After the development of meson spectroscopy for about ten years, many candidates of exotic hadrons were reported. To understand their possible natures, the “type-II” model (14) was proposed in Ref. [97]. According to this model, the is still the state with the flavor wave function ; the other combination and can form the and with ; the and are identified to be with and with , respectively. Within these assignments, the quark-quark interactions are assumed to be stronger while the quark-antiquark interactions are negligible, which implies that the assumed diquarks in tetraquarks are more compact than those in conventional baryons. In the case, six exotic vector states can be assigned: the has the same spin structure and spin interactions as the but has an additional -wave excitation between the diquark and the antidiquark, the is the orbitally excited , the is the -wave state, the or is identified as the orbitally excited state with , and the and are the first radial excitations of and , respectively. The observed decay was argued to be a natural consequence of the diquark-antidiquark description for these two states [132].
The authors of Ref. [133] assumed the diquark-antidiquark configuration. Noticing the similar color structure between the assumed diquark-antidiquark system and the quark-antiquark system, the authors estimated the masses in a similar way,
[TABLE]
where can also be a scalar or axialvector diquark. After considering the mass shifts due to the recombined meson-meson rescattering effects, they presented the masses of , , and tetraquark states.
After assuming the as a diquark-antidiquark type tetraquark and using it as an input, the authors of Ref. [134] investigated possible assignments for its partner states through the mass splittings calculated with the CMI method (19). The masses of two states with and one state with fit very nicely to the , , and .
In Ref. [135], Lebed and Polosa proposed that the [originally called , now called ] is the lightest diquark-antidiquark state in the “type-II” model (14). This identification may explain its observed width around 20 MeV. After treating the as a state and the and as states, they further interpreted the and as states. Moreover, the is interpreted as a or tetraquark state, while the may be an orbitally excited state. In 2006 the LHCb reported other two exotic structures in the channel [60], the and . In Ref. [136] Maiani et al assigned the as the ground tetraquark, and the and as two radially excited diquark-antidiquark states. They noticed that the corresponds to two almost degenerate states with and . Predictions of the other states within this picture will be tested by upcoming experiments.
Later Zhu et al presented a study on the hidden charm tetraquarks in the diquark-antidiquark model (13) [137]. The authors proposed the following assignments: the with can be explained as , the with is a companion of , the is the radial excitation of , the and may be the tetraquarks with the flavor wave function , the and may be the radial excitations of the tetraquarks with and respectively, and the may be an orbitally excited diquark-antidiquark tetraquark.
In Ref. [138], Ali et al studied the tetraquark states with an extended Hamiltonian of the model (14) by including the tensor coupling contributions. Only two scenarios are consistent with the Belle, BaBar, and BESIII data in the diquark-antidiquark picture. Each scenario contains four states. Scenario I contains , , , and , while Scenario II contains , , , and .
In Ref. [106], we updated estimations for the mass spectrum of the tetraquark states with the CMI method (19). After assuming the to be the lowest state, we confirmed the mass splittings obtained in Ref. [114]. We found that the should be the other tetraquark, and the [139] is probably a , while the interpretation of the and needs orbital or radial excitations. In this study we also confirmed the mass inequality (20), and found that is closer to than to . Similar features may be found in other systems.
In Ref. [113], we have systematically evaluated the mass spectra of the ( and ) states. Besides the estimation method (19), we also roughly estimated the tetraquark mass spectra by treating the mass of the as an input. The quark mass differences and need to be replaced by other values, since quarks get different masses in different systems. With this replacement method, higher masses than can be obtained. Such masses satisfy the inequality (20) and we treat them as realistic values. We show the results for states in Fig. 5. Results for other systems can also be found in this study [113]. It is possible to assign both the and to be the tetraquark states. With the defined measure reflecting effective quark interactions (see Sec. 2.12), we argued that this assignment is acceptable. Another possible tetraquark candidate is the , but we need further studies to answer whether the assigned structures are correct or not. At present, our results with this updated method are higher than the masses obtained in an early paper [77] where the , , and states were all found to be stable. The and are also weakly bound, although the binding energy in the sector is very weak.
2.5 CMI and
The study of the doubly heavy tetraquark states (, ) has a long history in the literature. Various approaches have been applied to such systems, and the most interesting state is the lowest isoscalar (or its antiparticle ), named . The and its partner tetraquark states are explicitly exotic. They do not suffer from possible annihilation effects and their inner quark interactions may provide enough attractions to stabilize the four-quark system. One may consult Refs. [140, 141, 142, 143, 144, 145, 77, 146, 147, 148, 149, 150, 151, 125, 152, 153, 154, 155, 156, 157, 158] for relevant investigations before 2003. Whether they are bound or not depends on models and adopted approximations. Generally speaking, the lowest states are narrow tetraquarks. Up to now, neither the nor any of its partner states have been observed in experiments. Recently, the confirmation of the state reignited interests in these exotic states since the and systems may be related to each other with the heavy quark symmetry (see the following Sec. 5). In this subsection, we mainly concentrate on recent developments in studies with CMI models since 2003.
In Ref. [88], Cui et al investigated the mass spectrum of the tetraquark states with Eq. (2.1) using the mass of as input. The obtained mass of the most interesting state with is 3786 MeV. This value is about 100 MeV below the threshold and thus this state should be narrow. The is also a stable state. The results in Ref. [125] obtained with Eq. (2.1) indicate that the lowest state is also below the threshold.
In Refs. [159, 160], Lee et al argued that the , , and states may be stable against rearrangement decays by using a diquark model with Eq. (19). In Refs. [161, 162], the masses of two states were estimated using the model (19), one with the color structure and the other with . The lower state is 71 MeV below the threshold and the higher is 54 MeV above the threshold. Their productions in collisions were also investigated.
We have performed a systematic investigation of the (, ) states with the CMI model (19) in Ref. [109], where the mixing effects of different color-spin structures were included. Considering the uncertainties in the estimation method, the mass of the state should be larger than 3780 MeV but not exceed 4007 MeV, that of should be between 10483 MeV and 10686 MeV, that of the lowest should be between 7041 MeV and 7256 MeV, and that of the lowest should be between 7106 MeV and 7321 MeV. Most theoretical predictions fell into these ranges, see comparison tables in Refs. [163, 109]. Our results indicate that the mixing effects for the systems are not so significant. If our lower bounds on the masses were reasonable, the and states with as well as and with are probably stable. Our results for () are consistent with an early study performed in Ref. [77]. Namely, such states with are bound. The binding results in these references [77, 109] for the with and with are also consistent. If the mass of the realistic satisfies the inequality (20) and is closer to 3780 MeV than to 4007 MeV, this tetraquark seems to be around the threshold, a situation similar to the .
Two investigations based on heavy quark symmetry analysis also show the possible existence of the doubly heavy tetraquark states [164, 165], which are partly motivated by LHCb observation of the double-charm baryon [36] since the role of the diquark within is similar to that in . Later, we will further introduce this issue in Sec. 5.
Shortly after the LHCb observation of the [36], also motivated by the success of the model (18) in predicting the mass [104], Karliner and Rosner investigated the masses of the , , and in Ref. [166]. They noticed that the is just about 7 MeV above the threshold, the is about 215 MeV below the threshold, and the lowest is 11 MeV below the threshold. Therefore, they concluded that there are no stable . But there are definitely stable states. The lifetime of the was estimated to be 367 fs.
In the diquark-antidiquark model (13), the authors of Ref. [167] studied the spectra and decay properties of the () states. For the lowest , its quantum numbers are and its mass is 3.6 GeV. This state is 140 (270) MeV below the () threshold. For the other states, the masses are 3.94 GeV, 3.97 GeV, and 4.04 GeV for the , , and , respectively. The masses of the states range between 3.82 GeV and 4.14 GeV. The masses of their partner states, and , can also be found in this work. In addition to the study of the mass spectrum, the authors discussed their - meson-meson and - antibaryon-baryon decay modes, too. However, their subsequent calculation with Eq. (2.1) gives an above-threshold and a below-threshold by adopting the diquark-antidiquark configuration [168].
2.6 CMI and
Generally speaking, many exotic states have been observed in particle experiments, some of which are good tetraquark candidates, as illustrated in the hidden-heavy case. However, it is still difficult to distinguish the compact multiquark picture from the molecular picture once two or more light quarks are involved in the states. For the full-heavy () states, there does not exist appropriate binding mechanism if they are treated as loosely bound molecules. If such a state is observed, it can be identified as a compact tetraquark state bound by short-range gluon-exchange interactions. In the literature, one may also find lots of theoretical studies on the heavy-full tetraquark states with various methods, such as MIT bag model and quark potential model [169, 170, 171, 172, 173, 140, 174, 175, 141, 142, 143, 176, 77, 148]. Similar to the studies, controversial conclusions on their existence were drawn when different models were applied. In recent years, experimental developments in search of exotic states reignited our interests in the fully heavy tetraquark states. Here, we focus on recent theoretical studies on the states with CMI models.
In Refs. [177, 178], the mass spectra of the diquark-antidiquark , , and states were studied with the Hamiltonian (13). The authors determined the coupling parameters by solving the nonrelativistic Schrödinger equation [179]. The lowest and with were found to be above the and thresholds, respectively. But the lowest states with other quantum numbers are all below relevant meson-meson thresholds. These states may be stable.
In Ref. [110], we have systematically investigated the mass spectrum of the () compact tetraquark states with the method (19). Our estimated values are lower bounds on realistic tetraquark masses, and it seems that no stable states exist. From our study, the only possible states, which are stable or relatively narrow, are and . The unbinding result for the states is consistent with an early CMI calculation in Ref. [77]. Our conclusions are consistent with those drawn from dynamical studies within constituent quark models (see Sec. 3). For example, the authors of Ref. [180] concluded that no narrow bound states are expected in experiments. For comparison, their mass values also satisfy the mass inequality (20). Since their masses are closer to our than to , which feature is similar to the case, the results in this potential calculation are probably closer to the realistic case.
In Ref. [105], Karliner et al studied the () states with the model (18) motivated by the QCD-string-junction picture in the diquark-antidiquark configuration. The predicted lowest-lying state with has a mass MeV, which is just below the threshold but can decay to . The predicted lowest has a mass MeV, which is just 28 MeV above the threshold and may have a width narrow enough to be observed.
Last year, LHCb searched for a possible exotic meson in the invariant mass distribution in Ref. [181]. Such a study was motivated by theoretical predictions and by the observation of production by CMS [182]. However, no significant excess was found in the mass range between 17.5 GeV and 20 GeV. The authors of Ref. [183] studied a possible di-bottomonium with the diquark-antidiquark picture (14). The obtained mass is approximately 100 MeV below the threshold. They argued that the decay of this state into is unlikely to be observed at LHC, which is consistent with the LHCb experiment [181].
2.7 CMI and
If the states do exist, the gluon-exchange interactions would provide enough binding forces for heavy quarks, and such interactions could also bind three heavy quarks and one light quark into compact (, ) states. In such systems, no long-range particle-exchange force exists, neither. Among these states, the and tetraquarks are explicitly exotic. Although the or tetraquarks look like the excited , , , or mesons, their excitation energy may not be explained by simple orbital or radial excitations if such states were observed in the future experiments. Similar arguments were used in the prediction of the hidden-charm pentaquarks [184]. Therefore, the exotic tetraquark nature of the would be easily identified once observed. Up to now, such systems are poorly studied.
More than twenty years ago, Silvestre-Brac investigated various tetraquark systems with the CMI method (19) [77]. For the systems, no bound states with but weakly bound , , and states with were found.
In addition to the and states, the authors of Ref. [88] also presented results for the states using the mass equation (2.1). From the obtained masses, many states, e.g. the , , and , were found to be below their respective lowest threshold of rearrangement decays and should be narrow.
We have performed a systematic study of the systems in Ref. [107] with the CMI method (19). Our results indicate that the lowest tetraquarks are only near-threshold bound or resonant states. Since our estimated values are just lower bounds on their realistic masses, all such tetraquarks are probably not stable. However, there may exist relatively narrow and thus detectable states such as and if the interactions inside the diquarks are attractive, although they are above the thresholds for the rearrangement decay channels.
2.8 CMI and
In 2015, LHCb observed two pentaquark-like resonances and with opposite -parities in the decay [185]. After that, in Ref. [186], the LHCb collaboration further applied a model-independent method to the same sample and the result supports their previous model-dependent observation. Later in Ref. [187], a full amplitude analysis of decays also supports the evidence for exotic hadron contributions from the , , and . Besides the LHCb, the CLAS12 at JLab may be a possible place to observe the states in photoproduction [188, 189]. Its spin and photocouplings may be measured with the future data, too. Now, the search for the LHCb pentaquark in the photo-production process in Hall C at Jefferson Lab has been approved [190, 191]. In future, these two states and their partner states can also be searched for by PANDA/FAIR [192], EIC, and JPARC, etc.
Before the observation of the LHCb states, many studies had been performed on the existence of possible hidden-charm pentaquarks [184, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204]. Here we show several examples. In Refs. [184, 193], Wu et al predicted the hidden-charm and resonances above 4 GeV in the molecule picture. In studying the predicted and in Ref. [196], the authors adopted three kinds of quark-quark hyperfine interactions, the flavor-spin interaction based on the meson exchange, the CMI based on one-gluon exchange, and the instanton-induced interaction based on the non-perturbative QCD vacuum structure. In Ref. [200], the discovery potential of the hidden-charm baryon resonances via photoproduction was investigated. In Ref. [204], the effects of the hidden charm in the reaction near threshold were discussed. In Ref. [205], there were some discussions about meson-baryon type hidden-charm pentaquarks with the help of flavor symmetry. After the LHCb observation, more and more studies appear on the nature of the states and relevant predictions. In Ref. [206], the and were interpreted as the loosely bound and molecules, respectively, in the one-boson-exchange model. In Ref. [207], the quantum numbers of these two states were suggested to be and , respectively. In Ref. [208], the masses of the pentaquarks with hidden beauty and strangeness were simply estimated with the heavy quark symmetry. For other studies of the hidden-charm pentaquark states before 2016, one may consult Refs. [1, 209]. Here we concentrate on recent developments with CMI models.
In Ref. [210], Maiani et al argued that the observed pentaquarks by LHCb are naturally expected diquark-diquark-antiquark states in an extended picture of Refs. [95, 129]. They assigned the compositions for the two states as
[TABLE]
In Ref. [211], the pentaquarks were assumed as diquark-diquark-antiquark states where all possible spins of diquarks were considered. Their masses were estimated with the formula taking into account the spin splitting,
[TABLE]
where and are the mass and angular momentum of the diquark-diquark system. According to the obtained pentaquark masses in this analysis, the is a good candidate of diquark-diquark-antidiquark state, but the was suggested to be a broad bump in the -wave resulting from rescatterings in the channel. This scheme differs from the classification in Refs. [210, 212]. In Ref. [213], the contributions from the possible diquark-antidiquark states [133] with masses around 41894300 MeV in the channel of the decay were discussed. The authors emphasised that such exotic mesons can imitate broad bumps in the channel. The study of the pentaquark spectrum was extended to the strange case in Ref. [214]. Further discussions in this scheme can be found in Refs. [215, 216].
In Ref. [100], Zhu and Qiao studied the pentaquark states in the diquark-triquark picture where both diquark and triquark are not compact objects. Their investigated pentaquarks have the structure or (). Their masses were calculated with the model Hamiltonian (15), where the diquark masses were extracted from the and the by assuming that they are diquark-antidiquark states, and the triquark masses were estimated by summing the quark masses. Their results suggest that the is a mixed state of three diquark-triquark structures with , and the can be interpreted as the diquark-triquark pentaquark state with . The authors also predicted the spectra of , , and , and discussed possible production and decay processes for these pentaquark states.
In the () diquark-diquark-antiquark picture, Ali et al investigated the pentaquark spectrum using the extended model of (14) in Ref. [212],
[TABLE]
where and are the orbital and spin angular momenta of the pentaquark state, respectively. They also adopted the extended model (13) by including additional spin-spin terms. Identifying the observed as the state , the masses of the other - and -wave pentaquark states were predicted. According to their results, identifying the state to be the observed is problematic because its production from the decay is suppressed by the heavy quark symmetry. They suggest the LHCb Collaboration reanalyze their data to search for the lower pentaquark in the range 4140-4130 MeV. In Ref. [217], the same authors extended their studies of the and pentaquark states to the cases. The anticipated discovery modes in -baryons decays were also discussed.
In Ref. [108], we have systematically studied the (, ) pentaquark states with the method (19), by assuming to be in the color representation. The LHCb states fall in the mass region of the studied system. Most such pentaquark states were found to have -wave open-heavy decays and their widths should be broad. In contrast, the states do not decay through the -wave, so their widths should not be very broad. We also found low-mass -type states with in this model. One may understand the existence of such a state from the CMI matrix elements. In the flavor symmetric limit, two states in the representation and three states in the representations have negative CMI matrix elements. Their mixing due to the symmetry breaking effectively provides additional attraction, which leads to a low-mass -type pentaquark finally. This result agrees with the study in Ref. [218] where detailed investigations with one-gluon exchange and instanton-induced interactions were performed.
In the study of hidden-charm pentaquark states, Buccella discussed the masses of the compact with the method similar to (19) in Ref. [219], where both and are color octet states. For the case, he predicted four pentaquarks around 4360 MeV, 4409 MeV, 4491 MeV, and 4560 MeV. The first two states couple strongly to , and they appear as a resonance with a mass 4380 MeV in this channel, while the third and forth states couple strongly to and , respectively. For the , the two color octets were assumed to be in -wave.
In Ref. [220], the authors assigned the observed as a compact pentaquark with in the lowest multiplet and considered its flavor-octet partner states using the extended Gürsey-Radicati (GR) formula [221],
[TABLE]
where is a scale parameter, is the isospin, is the hypercharge, is the eigenvalue of the Casimir operator, and is a counter of or quark. With the obtained parameters which were determined by fitting the masses of the ground state baryons, the authors found a state with the mass 4377 MeV and the same quantum numbers as the (charge, spin, parity). Possible discovery channels in bottom baryon decays and the partial decay widths for all the pentaquarks were discussed. In Ref. [222], an extended quark model was developed. The mass formula for baryons has the form
[TABLE]
where the spin-flavor operator has a generalized form of the GR formula and the orbital contribution term has the form with being the linear trajectory slope. Applying the formula to the () states, the authors found that all the hidden-charm pentaquarks are below the relevant -baryon thresholds. Therefore, the LHCb states were proposed to be the excited pentaquarks within this scheme.
2.9 CMI and
Up to now, the baryon is the unique one of the experimentally confirmed doubly heavy hadrons. Its observation motivated heated discussions on the possible (, ) tetraquarks. If one adds one more light quark, the existence of the pentaquark baryons might also be possible, because the light quark interactions inside the cluster may provide stronger attractions and play an important role in binding the five quarks, especially when the color-spin mixing effects are considered. Compared with the conventional baryons, the masses are probably not very large if the complicated interactions among quark components significantly lower their masses. As a result, the properties of the excited baryons may be affected by the states through coupled channel effects. Knowledge about the basic features for the spectra may be helpful for us to understand possible structures of heavy hadrons. In the literature, there exist some explorations of such systems which focus on the - molecules. However, such systems have not been paid enough attention to.
In Ref. [111], we have systematically estimated the masses of the possible pentaquark masses in the CMI method (19). Our results suggest that there exist many states below relevant open-heavy baryon-meson thresholds, and thus stable pentaquarks are expected. Since the realistic pentaquarks probably have higher masses than the values in Ref. [111]. It is not very clear where the - type baryon-meson thresholds locate since most of the baryons have not been observed yet. Further study is needed to verify whether such states are stable or not. Our results can be a helpful guide for deeper investigations.
To search for stable pentaquark states, the authors of Ref. [223] presented a systematic analysis on the flavor-color-spin structures of the heavy quark pentaquark systems with the method (19). The authors noticed that the state is perhaps the most stable one.
2.10 CMI and
The system is the mirror-type structure of the by exchanging the heavy and light flavors. Such states are related with the triply heavy conventional baryons, which have not been observed so far. This situation is similar to the relation between mesons and tetraquark states. With the potential observations of baryons in future measurements, it is also instructive to study the feature of the states. Besides, the study of the bound state and scattering problems induced by the observed baryons can help experiments to look for exotic phenomena.
In Ref. [112], we have discussed whether the compact pentaquark states are possible with the CMI method (19) by assuming that the is always a color-octet state, similar to Ref. [108]. Our results indicate that the compact pentaquarks should not be stable against their rearrangement decays. On the other hand, the near-threshold baryon-meson molecules are possible, see Sec. 4. Therefore, the gap between the lowest threshold and the lowest compact pentaquark leaves room for the identification of hadronic molecules. If experiments could observe an exotic resonance around some low-lying threshold, its nature as a molecule will be preferred over a compact pentaquark. As a by-product, we conjectured the maximum values for the realistic masses of the conventional triply heavy baryons, which are consistent with most theoretical predictions in the literature.
2.11 Mass constraints for and baryons
With the simplest CMI model (2.1), one cannot get accurate masses for conventional hadrons because of the adopted assumption, which has been illustrated in Sec. 2.3. The estimated multiquark masses with the model or its variants might also be far from the realistic values. Irrespective of the results, however, we find that some mass constraints for conventional hadrons and may be conjectured when estimating multiquark masses with the CMI method (19).
In some cases, there may be two or more reference thresholds when one discusses a system with the help of Eq. (19). Obviously, they lead to different multiquark masses. Then one may establish an inequality for the thresholds of these reference systems. If one reference system contains a hadron without measured mass, the inequality can be used to constrain the mass of that hadron, e.g. or . For convenience, we use and ( and ) to denote the two hadrons in the reference system (). What we observed from the results for various systems is that the multiquark masses estimated with are usually higher than those with , if both and contain heavy quarks while or does not contain heavy quarks, i.e.
[TABLE]
If does not contain heavy quarks, that means
[TABLE]
For example, in estimating the masses of the states, the reference systems may be or , and one expects that the latter system leads to lower pentaquark masses. Then we get the upper bound for the mass of (see Table 3). If another system is considered, e.g. , a different value of mass can be obtained. The minimum mass should be the upper bound. In Ref. [112], we have obtained mass inequalities between the triply heavy baryons and the doubly heavy baryons with this feature, which are also collected in Table 3.
With the help of Ref. [109], we may similarly constrain the masses of doubly heavy baryons with baryons and mesons. Such inequalities are also shown in Table 3. Numerically, the mass within this approach is smaller than the upper limit with Eq. (2.1), which is illustrated by comparing masses in Tables 2 and 3. Up to now, only the doubly charmed baryon was discovered by LHCb and no other baryons containing two or more heavy quarks are observed. It is obvious that the measured mass of is consistent with the upper bound for . For other baryon states, one may compare the constraints given here with the masses in the literature [104, 224, 225, 102, 112]. We find most of the theoretical predictions satisfy such constraints.
Since the masses in the second column of Table 3 are physical values, these constraints can be tested only after all the relevant ground-state baryons are observed. In this sense these inequalities are just conjectures. With these constraints, we may judge whether an observed heavy baryon is a ground state or not, i.e., once experiments observed states higher than such bounds, they should not be ground states.
2.12 Effective CMI for multiquark states
When studying the properties of the tetraquark states [110], we checked the mass shifts due to variations of coupling parameters. What the mass shift reflects is the effective color-spin interaction between a pair of quark components. For a multiquark state, its mass in the CMI model can be expressed as a function of the effective quark masses and coupling constants
[TABLE]
In the case that mixing among different color-spin structures is considered, this mass function cannot be expressed as a simple formula like Eq. (10). However, we may approximately use a linear form to show the relation between and the parameters. To do that, we consider effective chromomagnetic interactions by making a change on the coupling constants. When one reduces a parameter , the mass of the state becomes larger or smaller. If the effective CMI between the th quark component and the th quark component is attractive (repulsive), the reduction of will lead to a larger (smaller) . This effect can be reflected by a dimensionless measure, which was defined in Ref. [112],
[TABLE]
Here is the variation of the coupling parameter , and is the corresponding variation of the multiquark mass. When , becomes a constant, then we have
[TABLE]
where of Eq. (2.1) or of Eq. (19). The validity of this mass formula in the mixing case is easy to be checked with numerical results. Obviously, negative (positive) indicates attractive (repulsive) effective CMI.
In practice, the value of can be extracted by reducing slightly, while it can also be derived with the matrix and its eigenvalues. Consider a mixed state
[TABLE]
which is a mixture of different bases ’s with corresponding coefficients ’s. If represents an eigenvalue of , one has
[TABLE]
The Hamiltonian acting on bases can be expressed as (see Eq. (12))
[TABLE]
where can be read out from the obtained . Then
[TABLE]
i.e., . With this expression, one can understand how contributions from different structures affect the multiquark mass.
Although Eq. (31) has a linear form, it does not mean that the multiquark mass is linearly related with the coupling parameters, because ’s also depend on ’s. The effects on the multiquark mass due to uncertainties of coupling parameters can only be seen with the values of roughly. If the absolute value of a is large, it is necessary to reduce the uncertainty of the corresponding as much as possible. This possibility of large can appear in some multiquark states. An application of the measure is to qualitatively guess the stability of tetraquark states. For states with the configuration , their dominant decays are through the fall-apart or rearrangement mechanism. However, there still exist some special configurations which lead to relatively stable states although their masses may be high. For example, if both the and interactions in a state are effectively attractive while other quark-antiquark interactions are repulsive, such a state should be relatively difficult to fall apart into meson-meson decay channels. From the systematic study of the (, ) systems [113], the second highest states usually satisfy this condition, which is helpful for us to identify the possible nature of the exotic mesons. However, it is not easy to qualitatively discuss the relative stability of pentaquark states with ’s.
2.13 A short summary for CMI
The CMI models play an important role in understanding the multiquark systems. Although the Hamiltonian is simple, the models do catch the basic features of spectra since the mass splittings of hadrons reply on the basic controlling symmetries of the quark world. In recent years, developments about the CMI model investigations lead to series of important results. All types of heavy quark tetraquark and pentaquark systems that could exist in nature are listed in Table 4. We have reviewed the CMI studies on four types of tetraquark systems and three types of pentaquark systems in this section. Some interesting observations are:
- •
The states. The , , and observed in , the observed in , and the observed in can be consistently assigned as the compact tetraquark states with , , , , and , respectively within the framework of the chromomagnetic interactions. Additional high-lying and relatively narrow tetraquark states are also expected whose effective chromomagnetic interactions inside the diquarks are attractive while those between quark and antiquarks are repulsive. Uncovering the inner structures of the exotic states needs more experimental investigations and further theoretical studies from various aspects.
- •
The states. From the quark mass ratio dependence of the masses, the existence of such doubly heavy tetraquarks becomes inevitable when the ratio is small enough. From the CMI model calculations, it is still unclear whether the lowest state is bound or not depending on the details of the quark interactions. An interesting possibility would be that the mass of the is around the threshold, similar to the masses of and . On the other hand, the existence of the bound seems without doubt. The and states have only repulsive color-spin interactions inside the heavy diquarks. In contrast, not only the lowest but also heavier states may exist with relatively narrow widths because the effective color-spin interaction inside the and diquarks can be both attractive. The existence of the doubly heavy tetraquarks is certainly one of the most important issues in hadron physics.
- •
The and states. The lowest compact and states with from the CMI models seem unstable if all the color-spin structures are considered. Although the compact and states may also have rearrangement decay patterns, the effectively attractive color-spin interactions inside the diquarks probably lead to relatively narrow tetraquarks. Searching for states in the channel seems to be feasible at LHC. Moreover, the relatively narrow and tetraquarks probably exist.
- •
The states. Although the and states can be assigned as compact pentaquark states in the CMI models, their inner structures still remain unresolved. An interesting observation is that a low mass isoscalar pentaquark is possible which may be searched for in the and channels. Similarly, the isoscalar , , and states may be searched for in their rearrangement strong decay channels.
- •
The and states. All the , , and (denoted as ) states contain the diquark. If the were bound, the lowest pentaquarks might also be bound states. Since the other baryons except have not been observed, it is not clear whether low mass structures such as isoscalar , , and have rearrangement decay channels or not yet. The isoscalar , , and states are probably the most stable doubly charmed pentaquark. The triply heavy () structures seem unstable.
3 Constituent quark models
Various versions of nonrelativistic and relativistic constituent quark models can be found in the literature, which were proposed to understand hadron properties. Almost all of them incorporate both the short-range one-gluon-exchange (OGE) force and the term representing the color confinement, in either the coordinate or momentum space. Some of them include the additional flavor-dependent Goldstone-boson-exchange (GBE) force from the spontaneously broken chiral symmetry and/or the screening effects from the quark-antiquark pair creation. There are also models containing only GBE and confinement potentials. In some approaches, one has to adopt the diquark and triquark approximations in order to reduce the few-body problem to a two-body problem. In this section, we focus on recent studies of compact multiquark states within constituent quark models. We will also discuss some states which are candidates of the conventional or hybrid charmonium.
One of the widely adopted potentials was proposed by Godfrey and Isgur [33, 34], which was very successful in explaining conventional meson and baryon spectra. Now we denote it as the GI model. In studying multiquark states, there also exist many types of potentials, two of which were very popular. They were proposed by Bhaduri, Cohler, and Nogami in Ref. [226], and by Semay and Silvestre-Brac in Refs. [149, 227]. The latter references provided four potentials corresponding to different fitted parameters, which we shall call as BCN and Grenoble (AL1, AP1, AL2, and AP2) in the present review. The potential between th and th quark components in the BCN model reads
[TABLE]
while that in the Grenoble model has the form
[TABLE]
with .
In the following subsections, we shall separately review the theoretical progress of various constituent quark models on the , , , , and multiquark states.
3.1
Some neutral states are probably the conventional quarkonium. In some cases, the bare quarkonium states are strongly affected by the couple-channel effects. A typical example is the [228] which is affected strongly by the scattering states. In fact, the mixing between the and may be accounted for through the quark pair creation in the vacuum, which could be equivalently incorporated by the screened potential. The core contributions could be calculated even at the hadron level through Weinberg’s compositeness theorem [229]. On the other hand, the states do not have a quarkonium core.
3.1.1 Possible quarkonium or hybrid assignments
The authors of Ref. [230] studied the charmonium spectrum and electromagnetic transitions between charmonium states in a nonrelativistic quark potential model where both the linear and screened confinement potentials are considered. Their results indicate that the possibility of the being a dominant state cannot be excluded from its radiative decay properties. See also discussions in Ref. [231]. The as the is supported by its radiative decay properties. The mass splitting between the state and the state does not support the assignment of the as the state. In the screened potential model, the and may be good candidates of the and , respectively. The or might be identified as the state. From the subsequent study on the open-charm strong decays of the higher charmonium states in Ref. [232], the authors found that it is possible to assign the , , , and as , , , and , respectively. The assignment for the and in the same model as the excited states is difficult. The adopted model does not lead to a self-consistent description for the vector and the scalar .
In Ref. [233], the authors investigated possible assignments for the four structures in a coupled channel scheme by using a nonrelativistic constituent quark model [234, 235]. They found that the seems to be a cusp because of the near coincidence of the and thresholds, while the , , and appear as conventional , , and charmonia, respectively.
In Ref. [236], a relativistic Dirac potential model inspired by the OGE interaction was developed to study the charmonium spectrum in the momentum space. Due to the relativistic nature of the model, the spin-orbit, spin-spin, and tensor effects were all automatically included in the calculation. Two types of the scalar potentials with screening factors were considered. The authors determined the model parameters by fitting the masses of the eight resonances below the threshold and reproduced the overall structure of the charmonium spectrum. Analysing the predicted masses for high-lying resonances, they identified the , , and as states with =0, 1, and 2, respectively, while the is not well described and the is slightly higher than the measured mass. The was well reproduced, a result consistent with Ref. [230], but it is difficult to accommodate the , , , , , and from their results. The authors also suggested to perform deeper dynamical studies by considering effects such as molecule/tetraquark admixture and threshold effects.
In studying the properties of the exotic states, the authors of Ref. [237] discussed their possible heavy quarkonium assignments by solving the relativistic Dirac equation with a linear confinement potential. Their results indicate that the is a mixture of two -wave charmonia, and the and are admixed states of the - waves. In Ref. [238], the spectroscopy of the charmonium states was investigated with a Coulomb plus linear potential. From their calculations, the authors assigned the , , , and to be the , , , and , respectively.
In Ref. [239], the masses and decay widths of the vector charmonium states from to were calculated by using the instantaneous Bethe-Salpeter equation with a screened Cornell potential and the model by considering the mixing effects between the and states. Most of the results were consistent with experimental data. The authors further investigated the mixing between and states, and found that it is still possible to assign the and to be the - mixed charmonium states. The branching ratios of their decays into the were predicted to be small.
In a study of the heavy quarkonium hybrids [240] based on the strong coupling regime of pNRQCD (potential nonrelativistic QCD), the authors found that most of the isospin zero states fit well either as the hybrid or standard quarkonium candidates. The is compatible with the charmonium . The is compatible with the , but mixing with may have large contributions. The and are compatible with the charmonium states and , respectively, but may contribute here. The is compatible with a hybrid state, which mixes with a spin-1 charmonium. The and are compatible with the hybrid states. These states may be affected by the threshold. The and are compatible with the charmonium state having a dominant spin-0 hybrid component. The is compatible with the charmonium , which may be affected by the threshold. The is a hybrid or a conventional charmonium. The , , and are compatible with the spin-0 hybrid. The is compatible with a hybrid state, but its mixing with the spin-1 charmonium is little and it is difficult to understand its observation in the channel. The is compatible with the charmonium . The is compatible with a spin-0 hybrid state. The is compatible with the charmonium . In the bottom case, the is compatible with the bottomonium , the [241] is compatible with the spin-0 hybrid, and the with the bottomonium .
In Ref. [242], a nonrelativistic effective field theory describing heavy quarkonium hybrids was constructed. The authors discussed possible hybrid assignments for the exotic and states, such as the , , and . In Ref. [243], Miyamoto and Yasui studied the spectra and decay widths of the hybrid quarkonia in a hyperspherical coordinate approach, and discussed the possibility of the , , , , and being the hybrid states. According to their results, in the charmonium sector, the ground state of a magnetic gluon hybrid and the first excited state of an electric gluon hybrid lie close to the . In the bottomonium sector, the first excited electric hybrid and the ground magnetic hybrid appear a few hundred MeV above the . If one of the exotic meson candidates is a hybrid, its constituent gluon was expected to be magnetic. In a nonrelativistic quark model with the Cornell like potential, the authors of Ref. [244] found that assigning the as a mixture of the conventional and the lowest hybrid state may give a plausible explanation for its () production rates.
3.1.2 Schemes of the core admixed with coupled channels
Very recently, BESIII observed a new decay mode for the , , in the processes () and found [245]. This new measurement disfavors the interpretation for the presented in Ref. [246]. An early measured ratio by BaBar [247] did not support a pure molecule interpretation, while Belle’s result of [248] indicates that the contribution in is probably not large. These measurements strongly indicate that the probably contains a larger component than the . Inspired by these experimental measurements, there has been continuous progress of theoretical studies on the coupled nature of the and other and states in recent years.
In Ref. [249], the authors investigated the quarkonium contributions to the meson molecules like the and its heavy-quark spin-flavor partners. They found that the core produces an extra attraction, and thus less attractive interaction is needed to form a bound . But such an attraction does not appear in the sector. The content in the was found to be 1030%.
With the extended Friedrichs scheme, Zhou and Xiao of Ref. [250] studied the radially excited -wave charmonium spectrum. According to their analyses, the is dynamically generated by the coupling of the bare and continuum states, with large molecule components (64-85%). The observed has a much larger width than the charmonium although their masses are close. The state was predicted to be around 3890 MeV with a pole width about 44 MeV. Later in Ref. [251], they considered the isospin breaking effects of the with the same scheme, and obtained the reasonable ratio 0.580.92 between the and channels. In Ref. [252], the bottomonium counterpart of the around 10615 MeV and the around 10771 MeV were predicted. The reason for the non-observation of the signal in the [253] and [254] channels was also discussed.
In Ref. [255], the coupled channel effects on the heavy quarkonium states like and were studied by considering the meson loop corrections. The authors found that the is a superposition of a core and meson-meson continuum components, while the should be a pure bottomonium. The authors of Ref. [256] presented an unquenched method to consider the coupling between the quark-antiquark states and meson-meson states for a given . After the investigation on channels related with the , the authors got three states. The state that can be realized as the has a dominant component and a posterior . The state identified as the has a predominant component. The third state is an almost pure state. The coupled channel effects for the bottomonium with realistic wave were also studied under the framework of model in Ref. [257].
By assuming the as a mixing state of the 2P charmonium and molecule, the authors of Ref. [258] analyzed its radiative decays into and in an effective field theory. They found that a wide range of charmonium probability in the is consistent with the experimental ratio . In the case of the destructive interferences between the long-range meson loops and the short-range counter-term, a strong constraint on the admixture is found ( is expected).
To understand the nature of , the spectral function of the charmonium state coupled to mesons was studied in Ref. [259]. The authors found two poles in the complex plane, one pole corresponding to the object and the other virtual pole just below the threshold. The as a mixed object emerges from an interplay between a state and the thresholds. With their approach, the existence of is not possible without the seed charmonium state.
In Ref. [260], the was treated as an -wave molecule containing a small charmonium component. Extracting parameters from experimental results, the authors found that the physical wave function of the is . In Ref. [261], Lu, Anwar, and Zou studied the coupling between the charmonium structures and relevant meson-antimeson structures with the pair creation model. Their results support the picture that the is a molecule having a non-negligible core.
3.1.3 Tetraquark states
In the chiral SU(3) quark model [262], the masses of the (, ) states with and were calculated. It was impossible to assign the and [now called ] as pure tetraquark states with and , respectively. From this study, narrow tetraquarks in the hidden-bottom sector are possible. The study of the spectrum in the framework of the constituent quark model in Ref. [263] did not show any bound four-quark states for low-lying , ruling out the possibility that the is a compact tetraquark system, if no additional correlations are considered in simple quark models.
In Ref. [264], the , , , , and were studied within a coupled-channel Schrödinger model by coupling the and excited quark-antiquark pairs and their OZI-allowed decay channels or . Poles matching experimental data for these mesons were all found. In Ref. [265], the nature of the and other exotic states in the bottom sector was discussed in the diquark-antidiquark configuration in a quark potential model. The authors found that the is probably a radially excited diquark-antidiquark state, and that the might be the tetraquark rather than a conventional state.
In Ref. [266], Lu and Dong presented a study for the diquark-antidiquark states in the GI model by including the color screening effects. In this model, the can be regarded as the type tetraquark state, and the as the first radially excited or state, where () denotes the scalar (axial-vector) diquark. They also studied the case that the internal orbital or radial excitation is allowed for diquarks [267, 268]. Their results indicate that the can be explained as the tetraquark composed of one diquark and one antidiquark, while the is a tetraquark composed of one and one . However, there is no assignment for the in this model, which is proposed to be the candidate of the charmonium , supported by the decay width calculation in Ref. [269] and the arguments in Ref. [270]. Their results also indicate that it is possible to assign the as the lightest state, which supports the proposal given in Ref. [135]. Besides, some tetraquark masses are close to those of the and .
In Ref. [271], possible molecular states were studied in the BCN model and a chiral constituent quark model by including the -channel one gluon exchange. The authors found the bound states in the charm sector and bound states in the bottom sector. Good candidates for the and were then obtained. In another work [272], the same authors considered the -wave (, ) systems with two chiral constituent quark models in the meson-antimeson picture. From their results, several bound states were found. If the hidden-color channels are also included, a bound state with was also found. The obtained and states with can be related to the and , respectively.
In Ref. [273], the investigation for the spectroscopy of the hidden-charm and tetraquarks was performed in a relativized diquark model. The authors found possible assignments for the , , , , , , , , and in the sector, and for the , , and in the sector, but did not for the and .
In the dynamical diquark model, the authors of Ref. [274] presented a study on the spectrum of charmoniumlike tetraquarks by solving the Schrödinger equation with the Born-Oppenheimer potentials calculated numerically on the lattice. Choosing the or as a reference diquark-antidiquark state, the authors obtained a spectrum that agrees well with the observed charmoniumlike states.
Yang and Ping adopted a chiral quark model with the exchange of , , and and considered the diquark-antidiquark and meson-meson configurations in Ref. [275]. They investigated the spectrum of states and found that the and are good candidates of the compact tetraquark states with and , respectively. The can be explained as the radial excited tetraquark state with . However, no appropriate matching states for and were found.
3.2
The doubly heavy tetraquark systems are particularly interesting depending on the spatial configurations of two heavy quarks. In the extreme case that the pair is in the color anti-triplet and stays very close to each other, the compact heavy quark pair acts like one static color source. Two light anti-quarks circle around this point-like color source. This configuration is very similar to the Helium atom in QED. We denote it as the “QCD Helium atom”. On the other hand, if the heavy quarks are well separated from each other, the pair is shared by the two heavy quarks. This is the QCD valence bond. Such a configuration is the QCD analogue of the hydrogen molecule in QED. We denote it as the ”QCD Hydrogen molecule”. In general, the tetraquark system may be the superposition of the atomic and molecular structures.
In Ref. [276], using the BCN potential and the Grenoble AL1 potential, the authors studied the binding problem of the (the lowest isoscalar , see Sec. 2.5). According to their calculation, the with the molecular structure is weakly bound against the threshold, but it can become atomic with the inclusion of the three-body force.
In Ref. [163], the masses of the tetraquark states (, ) were calculated in the diquark-antidiquark picture with a relativistic quark model. All the states are not bound, and only the state lies below the threshold. The authors of Ref. [277] studied the four-quark bound states in chiral quark model. Their calculation indicates that a state with is bound, but there is no bound state in the system. With three different quark models one of which is chiral quark model, the authors of Ref. [278] studied the (, ) spectrum. Their results show that only the state with is bound within these models.
The authors of Ref. [279] studied possible compact four-quark states with two different constituent quark models, one of which contains Goldstone boson exchanges between quarks. They noticed that the and are bound and should be narrow, which was consistent with their previous studies [280, 281, 282]. Further study of the exotic four-quark states in Ref. [283] showed two isoscalar bound states with and . In Ref. [284], the Bethe-Salpeter equations for the ground-state were established in the diquark-antidiquark picture. According to their numerical results, the and bound states with and the bound states with and should all be stable.
In Ref. [285], Park and Lee investigated the and states with the BCN potential by including the contribution from the - diquark-antidiquark color structure, which was neglected in Ref. [152]. The was found to be 100 MeV above the threshold, while the was about 100 MeV below the threshold. The - contribution was found to be negligible. Later in Ref. [286], the masses of the doubly heavy tetraquarks were updated after the observation of the . Now the tetraquark is bound by 121 MeV and the is bound by 7 MeV. In an investigation with one gluon exchange potential, the authors of Ref. [287] studied the existence of the bound states. They found that there are no stable and states, but stable states are possible.
The authors of Ref. [288] discussed approximations used in the multiquark studies, i.e., diquark, Born-Oppenheimer, Hall-Post inequalities, color-mixing, and spin-dependent corrections. They pointed out that the is at the edge of binding, and more delicate studies are still needed. For the , the spin effects or color mixing effects are needed to achieve a binding state.
3.3 and
The study in Ref. [289] with a nonrelativistic potential model indicates that the lowest all-charm tetraquark state is below the threshold. A study for the fully-heavy tetraquarks was performed in Ref. [290] with a nonrelativistic effective field theory (NREFT) at the leading order and a relativized diquark-antidiquark model. Both approaches give a ground state around 18.72 GeV with , which is below the threshold. Mass inequalities for all type () were also investigated.
The authors of Ref. [291] presented a calculation for the mass spectrum of the ground-state tetraquark states using a diffusion Monte-Carlo method to solve the non-relativistic many-body problem, whose potential is based on the flux-tube model. A state around 18.69 GeV was found to be about 100 MeV below the threshold.
In Ref. [292], the masses of the tetraquark states were calculated in a constituent quark model with the hyperspherical formalism. The authors noted that the state is 58 MeV above the threshold, state is 14 MeV above the threshold, and state is 22 MeV below the threshold. In Ref. [293], the authors investigated the spectroscopy of in a diquark-antidiquark configuration using a non-relativistic model, whose potential is a Cornell-inspired type. Their results show very compact tetraquarks, which are below their thresholds of spontaneous dissociation into low-lying charmonium pairs.
In Ref. [294], Richard et al found that full-charm and full-beauty tetraquarks are unbound in the chromoelectric model with additive potentials. If the naive color-additive model of confinement is replaced by a string-inspired interaction, however, the case might be favorable. In additional to studies about the states, the authors also discussed the existence of , , and tetraquarks in Ref. [287]. Because of lack of strongly separated mass scales, they concluded that no such bound states exist.
In a dynamical study with a nonrelativistic potential model in Ref. [180], the mass spectra of all types of full-heavy () tetraquark states were obtained. No bound states can be formed below the thresholds of any meson pairs or , and thus such states with narrow widths are not expected in experiments. This conclusion about spectrum confirmed the arguments from the CMI estimations given in Ref. [110]. In a study with a nonrelativistic chiral quark model [295], the analysis of states indicates that the states are all higher than the corresponding thresholds, in both meson-meson and diquark-antidiquark configurations. The results also confirmed the conclusions from the CMI calculation with the method (19) [110].
In Ref. [148], Silvestre-Brac and Semay did not find any bound states with a nonrelativistic quark model, whose potential was proposed by Bhaduri et al [226]. This result is roughly consistent with the estimations with the CMI model (19) [77].
3.4 and
Considering the quark delocalization color screening effects in a chiral quark model, Huang et al investigated the -like hidden-charm pentaquarks in the molecule configuration in Ref. [296]. According to their results, the can be explained as the molecule with . Other hidden-charm and hidden-bottom bound states were also investigated. But all the positive-parity states were found to be unbound. They studied the hidden-strangeness pentaquarks in Ref. [297]. In Ref. [298], the same group performed a dynamical calculation of the five-quark systems in the framework of a chiral quark model. The was again suggested to be a molecule. All the obtained positive-parity states are unbound in their calculation, unless the effective meson exchange is employed. The authors did not adopt the assignment for the as the molecule because of the inconsistent parity, although their masses are close to each other. They presented a similar study of the hidden-bottom pentaquarks in Ref. [299].
In Ref. [300], the solution for the relativistic five-quark equations favors the assignment for the quantum numbers of the pentaquarks and to be and , respectively.
In the molecule picture, the authors of Ref. [301] studied the LHCb states with a chiral constituent quark model [234, 235] and confirmed the existence of several structures near the masses. In Ref. [302], the isospin-half pentaquarks with the color-octet component were investigated with the quark cluster model. Three structures were found to be around the thresholds: one bound state, two resonances, and one large cusp. The authors argued that the may arise from these structures.
In Ref. [303], the authors studied the possible pentaquark nature of the with using a hyperfine+confinement potential. With the variational method, they found that the ground state is the isolated and state, and the excited state lies far above the observed . This observation led to the conclusion that the observed cannot be a compact pentaquark generated by the conventional two-body quark interactions. In Ref. [304], Stancu studied the effects of three-body chromoelectric interaction in pentaquark states and found that the three-body confining interaction can stabilize the pentaquark states with an appropriate sign. She suggested that perhaps the fall apart decay mode of the ground-state pentaquarks in Ref. [303] may be avoided if such an interaction is included.
In Ref. [305], Richard, Valcarce, and Vijande studied the hidden-charm pentaquarks in a constituent model with the Grenoble AL1 potential by solving the five-body problem. They found that the lowest states with and are below their lowest - and -wave thresholds, and , and thus are expected to be stable. In a recent work [306], Hiyama, Hosaka, Oka, and Richard studied the five-body scattering problem with a nonrelativistic quark model by including explicitly open channels such as , , and , etc. They also adopted the Grenoble potentials (AP1 and AL1). Their results are compatible with those obtained in Ref. [305]. According to their analysis, there does not exist any resonance at the energies of the LHCb states, while two narrow states are possible lying at 4690 MeV and 4920 MeV with and , respectively.
In Ref. [307], Stancu investigated the spectrum of the hidden-charm pentaquark states with an Goldstone-boson-exchange model, in which the mass splittings are accounted for by the flavor-spin hyperfine interaction. She found that the lowest positive-parity states with and are below the lowest negative-parity state with . The results accommodate a and a .
With the Bethe-Salpeter equation, the spectrum of the heavy pentaquark states (, ) were studied in the diquark-diquark-antiquark configuration in Ref. [308]. The masses of the like hidden-charm pentaquarks and states are MeV higher than the experimental values. The authors also calculated the masses of states. The mass of the lowest state is around 4.54 GeV which is larger than our result obtained in the CMI model [111].
In Ref. [274], the charmoniumlike pentaquark states were investigated in the dynamical diquark model by assuming that they are diquark-triquark states. Using the observed and as inputs, the authors predicted a number of other unobserved states. Especially the ground hidden-charm pentaquarks were below the even threshold.
The authors of Ref. [309] presented a study of the magnetic moments of the hidden-charm pentaquark states with the isospin and in the molecular, diquark-triquark, and diquark-diquark-antiquark models. The obtained magnetic moments are different, although the description for pentaquark masses and decay patterns can be achieved in all these models. Their results may be used to understand the inner structures and distinguish phenomenological models. In Ref. [310], the axial charges of the hidden-charm pentaquarks were studied in the chiral quark model. Their results can also be used to distinguish the underlying structures of the states.
3.5 Color flux-tube model
The flux-tube model [311, 312] was first extracted from the strong-coupling Hamiltonian lattice formulation [313] and the early descriptions of flux tubes as cylindrical bags of colored fields [314]. In the flux-tube model, the gluonic excitation was considered as the transverse vibration of the string-like flux tube between a quark-antiquark pair. The ordinary quark model was contained under the zero angular momentum approximation, while exotica can also be achieved for the non-zero angular momentum [312, 315].
The string potential in Hamiltonian was usually in a standard linear form to investigate the confinement phenomenon [315]. Such a flux-tube model has been extensively used to study the hadron masses, decays and production characteristics [315, 316, 317]. Recently, a naive color flux-tube model was developed by considering a harmonic-like confinement potential rather than a linear one to study the tetraquark states [318, 319].
In Ref. [320], Deng et. al. studied the lowest charged tetraquark states ( and ) in the framework of the color flux-tube model with a four-body confinement potential. Within the color flux-tube model, the total hamiltonian of the diquark-antidiquark can be expressed as,
[TABLE]
in which or , is the center-of-mass kinetic energy, is the momentum of the -th quark, is the quadratic confinement potential. The interaction between two different quarks is described by , which contains the one-boson-exchange potential , the -meson exchange potential , the one-gluon-exchange potential . In the diquark-antidiquark configuration, the energy of the tetraquark states was then calculated by solving the four-body Schrödinger equation. The numerical results indicated that some compact resonance states were formed with the four-body quadratic potential. These resonances can not decay into and but into and through the breakdown and recombination of the flux tubes. The authors interpreted the charged states and as the -wave tetraquarks with quantum numbers and and 2, respectively.
Later in Ref. [321], they systematically studied the charged family in the color flux-tube model and considered the spin-orbit interactions between quarks. Their investigation indicated that the , , , , , and can be described as the tetraquark states with and of and , and , and , and , and , and , respectively. However, the two heavier charged states and can not be explained as tetraquark states in this model. In Ref. [322], Zhou, Deng and Ping also employed the flux-tube model to investigate the tetraquark states with . The vector states , , , and were described as the , , and hidden-charm tetraquarks, respectively.
In Ref. [323], the hidden-charm pentaquark states were investigated within the framework of the color flux-tube model including a five-body confinement potential. The results showed that the main component of the can be described as a compact pentaquark state with the pentagonal structure and . However, the state can not be interpreted as a compact pentaquark, since the masses of the positive parity states were extracted much higher than those of the and in the color flux-tube model.
Recently, they also studied the stability of the doubly heavy tetraquark states in the color flux-tube model [324]. They solved the four-body Schrödinger equation to calculate the energies of all doubly heavy tetraquark states with both the color configurations and . Their numerical results showed that the states with and with were the most promising stable doubly heavy tetraquark states against strong interaction. The doubly heavy tetraquark states have also been addressed by a fully unitary and microscopic quark model with a triple string flip-flop potential [325]. They found several tetraquark bound states and resonances and calculated their masses and decay widths by computing the matrix and finding the pole positions in the complex energy plane.
3.6 A short summary
We summarize some interesting observations for the multiquark states from quark model investigations in recent years:
- •
The physical should be a mixed state due to the strong channel coupling between the bare quark model state and the scattering state. Assigning the as the charmonium is not supported by all quark model calculations. The can be assigned as the or (affected by ).
- •
In the charmonium or hybrid picture, the is a mixed , a hybrid affected by , or a cusp. The is the affected by . The is the or a hybrid. The is the or a hybrid. The is the or .
- •
Perhaps both the mixed - states and the states contribute to the observed and , while the and are dominantly mixed - even charmonia. Their masses can also be reproduced in tetraquark models.
- •
The masses of and may be obtained in both molecule picture and diquark-antidiquark picture. The and should be molecules, although the diquark-antidiquark interpretation for is also possible. The nature of is still unclear.
- •
Whether the lowest tetraquark state is above or below the threshold is still an open question, but the lowest state should lie below the threshold.
- •
Few models favor the bound and states. A bound state in the case might be favorable. Bound states should not exist.
- •
The quark model calculations favor the molecule interpretations for the and .
4 Meson exchange and scattering methods
In the study of the hadronic molecules, almost all the meson-exchange models are constructed at the hadron level. One first derives the effective boson-exchange potentials in coordinate or momentum space, and then solves the bound state problem or scattering problem of two hadrons. From the obtained binding energy or scattering phase shifts, one extracts the resonance information. Such a formalism is a straightforward extension of the traditional meson exchange models in nuclear force.
There are meson-exchange models based on the quark-level interactions, e.g. based on the chiral quark model. After integrating out the inner degrees of freedom inside the hadrons, hadron-level effective meson-exchange potentials can be obtained. An example is illustrated in Ref. [326]. Of course, the molecule problems can also be dealt with at the quark level directly.
The color force of the van de Waals type is not described by the meson-exchange but by two-gluon-exchange interactions. This type interaction was applied to understand the nature of the exotic states in the hadroquarkonium picture where a heavy quarkonium is embedded in a light hadronic matter. This hypothesis dated back to early nineties [327]. Although the interaction between ground-state hadrons may be weak and no bound state can be formed, probably such a force may result in bound states or resonances for radially and orbitally excited hadrons.
In the chiral limit, the eight pseudoscalar meson masses vanish as the up, down and strangeness current quark masses go to zero. In contrast, the masses of the nucleon, rho meson or heavy mesons/baryons etc remain finite or even large in the chiral limit. They are sometimes denoted as the matter fields. If they contain the light quark degree of freedom, the matter fields will interact with the pseudoscalar mesons. The rigorous framework is the well-known chiral perturbation theory.
In order to ensure the chiral invariance, one should replace the simple derivative operator by the chirally covariant derivative operator in the construction of the kinetic energy of the matter fields at the leading order Lagrangian. The chirally covariant derivative operator contains the chiral connection. In other words, the kinetic term of the matter fields always induce the contact seagull term in the form of , where the is the matter field, is the pseudoscalar meson, and is the pion decay constant. The model-independent seagull term is a direct manifestation of the chiral symmetry, which does not introduce new coupling constants.
Now we move on to the scattering process of the pseudoscalar meson and the matter fields. The physical scattering amplitude is composed of the seagull amplitude, the corrections from higher-order tree-level diagrams and loop corrections. In other words, the seagull amplitude does not saturate the whole amplitude. Sometimes the seagull amplitude may be not dominant.
In the chiral unitary approach, only the seagull amplitude is iterated (unitarized) through the Bethe-Salpeter or Lippmann-Schwinger equation because of its easy-be-dealt-with and model-independent expression. From the pole positions of the unitarized amplitudes, one can extract information of bound states or resonances. Such states are called dynamically generated states in the literature. Since the seagull amplitude may miss some important physics as illustrated above, the chiral unitary approach sometimes generates spurious signals.
In the scattering of the light vector meson and the matter fields, one may treat the vector meson as the “gauge” field and introduce the hidden local symmetry into the Lagrangian. Then the chiral unitary formalism can be extended to the processes involving light vector mesons in a similar way.
If only the pions are exchanged between two heavy hadrons, the potentials can be organized order by order in the framework of the chiral perturbation theory (PT). One first constructs the tree-level Lagrangians order by order. Then one computes the higher order chiral corrections from the loop diagrams such as the triangle diagram, the box diagrams etc where only the two-particle irreducible contributions are kept. From the above scattering amplitudes, one derives the effective potentials at the given chiral order which enters the Schrödinger equation or Lippmann-Schwinger equation in the iteration. Again, one can extract the binding energy after solving the Schrödinger equation. One may also extract the resonance parameters either from the phase shifts or from the poles of the scattering amplitudes. This approach mimics the modern nuclear force in terms of PT which is pioneered by Weinberg [328, 329].
4.1
In Ref. [1], we performed a comprehensive review of the hidden-charm tetraquark states. In this section, we mainly introduce the theoretical progress on the systems with the configuration after 2016. Different theoretical groups still focused on the interaction between a charmed meson and an anticharmed meson. The central issue remains to answer whether these observed charmoniumlike states can be understood in the hidden-charm molecular scenario.
In Ref. [330], the mass spectra for the heavy meson-antimeson bound states were discussed with the Hellmann plus one-pion-exchange potential. The results indicate that the is not a pure molecule. The authors of Ref. [331] discussed possible hadronic bound states formed by heavy-light mesons and heavy-light antimesons in the one-pion-exchange approach, where a heavy-light meson belongs to the multiplet or its chiral partner multiplet . The coupled channel results show that only the charm isosinglet molecules with are bound, which can be identified as the . The analysis indicates that the reported isotriplet is at best a near threshold resonance, and the at best also a threshold effect. The study for the hidden-bottom meson-antimeson systems leads to the identification of the isotriplet exotic molecule with as a mixture of the and . A bound isosinglet hidden-bottom molecule with analogous to the was proposed.
A recent systematic study of the possible hadronic molecules composed of the -wave heavy-light mesons can be found in Ref. [332]. The assignment of the as a shallow molecule is supported, while assigning the , , , and as molecules is not favorable. Moreover, the cannot be assigned as a molecule composed of , neither.
With the Bethe-Salpeter equation, the authors of Ref. [333] calculated the mass of the in the molecule picture. They assigned this meson to be a mixed state of three pure molecule states , , and . In Ref. [334], the molecule in the Bethe-Salpeter equation approach was discussed with the ladder and instantaneous approximations by considering , , , , and exchanges. The can be a molecular bound state.
The possibility of the as a -wave state in a quasipotential Bethe-Salpeter equation approach was discussed in Ref. [335]. A pole at MeV was produced through the -wave interaction with the channel. If this state can be interpreted as a hadronic molecule, an -wave bound state below the threshold should also exist.
In Ref. [336], He and Chen investigated the interaction in a quasipotential Bethe-Salpeter equation approach by considering one-pion-exchange interaction. They obtained a bound state at 4384 MeV with and a resonance at 4.461+39 MeV with . The former (latter) state was assigned to be the (Z(4430)). Later in [337], Chen, Xiao, and He found that the molecule assignment for the can naturally explain its observation in and its absence in .
In a study with chiral unitary method in Ref. [338], Sakai, Roca, and Oset investigated the - and - interactions. They found isoscalar bound states above 7 GeV for the , and , and systems. In Ref. [339], the interaction was investigated with an extended hidden gauge symmetry in the chiral unitary model. A bound state slightly lower than the threshold was dynamically generated in the isoscalar channel, which was related to the . An isoscalar bound state was also found.
These authors in Ref. [340] studied the and structures in a constituent quark model by considering the coupled-channels , , and . They found that these two structures are virtual states, whose effects as shown in the production line shapes can be seen as the threshold cusps.
In the hadrocharmonium picture, the authors of Ref. [341] studied possible bound states. They authors found that the obtained -wave vector-vector bound state corresponds to a mass-degenerate multiples, and the can be such a state. The degeneracy indicates that two more structures around the should exist if this picture is correct. In Ref. [342], the bound states formed by or and isoscalar mesons in the hadrocharmonium picture were investigated. From their results, the , , , , and are , , , , and states, respectively. A different assignment was proposed in Ref. [343] to interpret the , , and as dynamically generated resonances in the coupled channels of and , with some relevant ones.
Voloshin interpreted the and as two states of hadrocharmonium related by the charm quark spin symmetry [344]. Later in Ref. [345], he predicted the existence of the strange hadrocharmonium resonances and which decay dominantly into and , respectively.
In a study of the scattering problem related with near-threshold heavy-flavor resonances [346], the authors found that both the and other hadronic degrees of freedom are equally important inside the . The , , and components inside the , , and , respectively, are not so important. For the , this study favors its interpretation as the molecule. A previous study with the same approach [347] indicates that both the and are dominated by the component, 76% and 68 %, respectively. See also discussions in Ref. [348].
4.2
In Ref. [349], besides the , the authors also studied the doubly charmed vector-vector meson states within the hidden gauge formalism in a coupled channel unitary approach. They obtained a pole around 3970 MeV in the channel, which is 100 MeV above the threshold. Later in an extended work [338], the - and - interactions were systematically investigated in order to see whether there are dynamically generated states. The authors noted that the bound states in , , and systems are all possible. In Ref. [350], the doubly charmed exotic states as meson-meson molecules were investigated. The authors studied the scattering problem by solving the Lippmann-Schwinger equation, where the effective potentials were derived from the chiral constituent quark model [351]. Their results suggest the existence of a stable state with .
The authors of Ref. [352] presented a study of possible doubly charm and doubly bottom molecular states composed of a pair of heavy mesons, by solving the coupled Schrödinger equations through the one-boson-exchange (OBEP) method. The authors got bound and/or resonant states of various quantum numbers up to . In the case, they obtained the and bound states with binding energies about tens of MeV. In another coupled channel calculation for the molecule problem of - (, ) systems in the OBEP framework [353], the and bound states were also obtained. In addition, the molecule were found to be bound together with several other molecules of various configurations. Studies of their partner molecules can be found in Ref. [354].
Using the lattice potentials from Refs. [355, 356, 357] and the Born-Oppenheimer approximation, Bicudo et al studied the tetraquark states and found an state about 90 MeV below the threshold. Later in Ref. [358], an resonance with MeV and MeV was found. The authors of Ref. [359] discussed the possible states with the exotic doubly charmed configurations and based on the kaon-exchange interactions. The authors obtained bound states in the and channels.
In the framework of the chiral effective field theory, the authors of Ref. [360] derived the effective potentials up to the order , where can be the momentum of the pion, residual momentum of heavy mesons, or the - mass splitting. As shown in Fig. 6, the two pion exchange contribution which is equivalent to phenomenological and exchange contribution, was considered in their calculation. With the obtained potentials in the coordinate space, the authors found an bound state by solving the Schrödinger equation. Further study in the same framework [361] showed that two bound states, and , are possible.
4.3
In 2011, the authors of Ref. [195] predicted the existence of the hidden-charm molecular pentaquarks in the framework of one boson exchange model. Later, the hidden-charm and states were announced by LHCb in the process [185], which inspired extensive discussions of the molecular pentaquarks with the content. The research status of the hidden-charm molecular pentaquarks composed of an S-wave charmed baryon and an S-wave anti-charmed meson was covered in great details in Ref. [1]. Since the investigation along the same line continues, we discuss the progress in the past three years in this section. The coupled channel bound states problem of was studied with the one-pion-exchange model in Ref. [362]. The authors found the existence of one or two bound states with the binding energy of several MeV below the threshold. This result indicates that the can be understood as a loosely bound molecule. The extended analysis to the , , and cases also shows pentaquarks lying about 10 MeV below the corresponding thresholds. For the , it was proposed to be the Feshbach resonance. In Ref. [363], an extended study for the coupled system of -- with the complex scaling method was performed. The potential included both the pion exchange and meson exchange force. The contribution of the channel was very small, and a pentaquark candidate with a value close to the was obtained.
In Ref. [364], the coupled - meson-baryon systems were investigated in a one-meson-exchange potential model. From the obtained bound and resonant states with , the authors concluded that the assignments for the and are and , respectively, in agreement with the LHCb results. In Ref. [365], the study was extended to include coupling with compact five-quark channels and to the hidden-bottom case. The short-range effective potential from the compact five-quark channel was found to be attractive, which plays an important role in producing the states.
The author of Ref. [335] explored the interactions with a Bethe-Saltpeter equation method in the hadronic molecule picture by considering the channel coupling to . A pole near the threshold with can be found. When this -wave state is produced near threshold, the -wave pole corresponding to locates around 4390 MeV. These two poles are suggested to be related with the and , respectively. Later in Ref. [366], by replacing the quark with the quark, the author interpreted the and as the strange molecular partners of the and , respectively.
In Ref. [367], possible strange hidden-charm pentaquarks in the and systems were investigated in a one-boson-exchange model. Promising candidates including and states were predicted. This hidden-charm pentaquark with strangeness is accessible at LHCb.
The authors of Ref. [368] discussed the scale invariance in hadron molecules like the - ( denotes ) coupled state, which is related to the . The nearly on-shell pion exchange transition can generate a long-range potential which can lead to approximate scale invariance for the equation describing the molecule if the attraction is strong enough. As a result, an Efimov-like geometrical spectrum in two-hadron systems is possible. The molecules, - of , - of , - of , and - of , seem to satisfy the condition for the approximate scale invariance.
In Ref. [369], Chen, Hosaka, and Liu explored the intermediate- and short-range forces in the framework of one-boson-exchange model and studied various hadron-hadron bound state problems. They found -wave , , , and molecules. In the study of bound state problems in Ref. [370], the authors proposed a one-boson-exchange potential model by adding a screen Yukawa-like potential. From their results, the was interpreted as an molecule.
The authors of Ref. [371] presented an exploratory study of possible - and - resonances by extending the Jülich-Bonn dynamical coupled-channel framework, which is a unitary meson-baryon exchange model. The authors found one pole in each partial wave up to in the hidden-charm case and several poles in partial waves up to in the hidden-bottom case. The very narrow pole with in the hidden-charm case was interpreted as a bound state, which was predicted in Refs. [184, 193].
In Ref. [372], Huang and Ping investigated the hidden-charm and hidden-bottom pentaquark resonances in hadron-hadron scattering processes in the framework of the quark delocalization color screening model. They found a few narrow hidden-charm resonances above 4.2 GeV and some narrow hidden-bottom states above 11 GeV. Besides, they also noticed bound - states from the behavior of the low-energy phase shifts obtained with coupled channel calculations.
In Ref. [373], the authors interpreted the as a bound state generated by the charmonium-nucleon interaction in terms of charmonium chromoelectric polarizabilities and the nucleon energy-momentum distribution. They obtained two almost degenerate narrow states at the position of the with and . The authors of Ref. [374] confirmed the results and predicted the isospin-3/2 narrow bound states around 4.5 GeV and broader resonances around 4.9 GeV in the framework of the Skyrme model. Later in Ref. [375], with a QCD inspired approach, the authors interpreted the as a bound state with , which is a member of one of two almost degenerate hidden-charm baryon octets. According to their study, one has to assign the as a state but no natural interpretation was found in the hadroquarkonium picture. The authors also studied the other hadroquarkonia and compared their results with those in the one-pion exchange approach. They did not find any bound state but found an inconclusive state.
A study of the LHCb states in the hadroquarkonium picture was carried out in Ref. [327]. The authors solved the Schrödinger equation with the potential from QCD multiple expansion and got spectrum. The and can be interpreted as the and bound states, respectively. They also predicted the hadroquarkonium states in the hidden-bottom sector.
Unfortunately, no evidence of the exotic state was observed by LHCb in the decay into [376]. The hadrocharmonium picture is being challenged.
After the observation of and , their possible inner structures were proposed either as tightly bound pentaquark states [210, 211, 377, 377, 378, 98, 100] or loosely bound molecular states [206, 207]. However, the experimental data in Ref. [185] was unable to distinguish these configurations. In 26 March 2019, at the Rencontres de Moriond QCD conference, the LHCb Collaboration reported the observation of three new pentaquarks [68]. The observed may correspond to the molecule with , while the and can be identified as the molecular states with and . The current measurement strongly supports the molecular hidden-charm pentaquarks predicted in Refs. [195, 184, 193], which shall become a milestone in the exploration of the multiquark matter. After their announcement, these three pentaquark states were immediately studied in the QCD sum rules [379], the one-boson-exchange(OBE) model [380], a contact-range effective field theory and a contact-range effective field theory [381]. The isospin breaking decay pattern was studied in [382]. In Ref. [380], the authors studied these new pentaquarks in a direct calculation with the OBE model. Their result supported that the , and may correspond to the loosely bound molecule with (), with () and with (), respectively.
4.4
In Refs. [205, 383], various baryon-meson molecular states were studied by considering vector-meson exchange forces as well as the coupled channel effects with the broken flavor symmetry. In the charm number sector, the authors predicted more than 10 states with various quantum numbers. The authors of Ref. [384] predicted a possible molecule, and the authors of Ref. [385] systematically studied dynamically generated double-charm baryons in a coupled channel unitary model. In a chiral unitary approach, the authors of Ref. [386] studied various dynamically generated double-charm meson-baryon states through vector meson exchanges interactions. They found ten molecular pentaquark states.
In Ref. [387], the and interactions and related bound state problems were studied. The authors derived two-pion-exchange potentials and regularized the divergence with a phenomenological cutoff. With solutions from the Schrödinger equation, they concluded that the bound state is possible, but the bound state can only be prudently expected. In Ref. [369], Chen, Hosaka, and Liu explored the intermediate- and short-range forces in the framework of one-boson-exchange model and studied various hadron-hadron bound state problems. They found -wave , , , and molecules.
In addition to the study of the coupled system -- in Ref. [363], Shimizu and Harada also considered the coupled system - in the channel. The authors used the same one-pion-exchange potential after the antimesons are replaced by mesons. They noted that there exists a doubly charmed baryon of the type as a hadronic molecule, named , whose mass and width are close to those of .
The observation of the double-charm baryon by LHCb stimulated the investigations of the molecule-type states in the literature. In Ref. [388], Guo considered the scattering of the baryons and Goldstone bosons in a chiral effective theory. According to their calculation, there are a pair of bound and virtual states near the threshold, indicating that there exists an interesting exotic state with the quark component . A resonance pole around the threshold, a bound state below the threshold, and two resonances around the and thresholds were also found. The authors of Ref. [370] presented a study of various meson-baryon and baryon-baryon/antibaryon states in a one-boson-exchange potential model compensated by a screened Yukawa-like potential. A series of double-heavy molecules were obtained.
4.5
In addition to the molecular bound or resonant states with charm number , , , and , the states with were also studied in Refs. [205, 383]. A flavor singlet bound state with () with mass around 4.3 (4.3) GeV is possible, which results from the scattering of a baryon and a meson as well as a baryon and a meson. In the systematic study of triply charmed dynamically generated baryons in Ref. [385], the authors found one bound state around 4.4 GeV and one bound state around 4.5 GeV by coupling - and - channels.
The discovery of the charmonium-like states and and pushed the exploration of the hadronic molecular states in the past decade, which was based on the interaction between charmed hadron and anti-charmed hadron. The observation of the double-charm baryon drives us to further explore the interaction between the double-charm baryon and charmed meson, which is a natural extension of the one boson exchange model as shown in Fig. 7. The same idea was further extended to the investigation of the interaction between the double-charm baryon and charm baryon in Ref. [389].
Due to the above motivation, in a one-boson-exchange potential model, possible molecular states formed by the and a charmed state and were discussed in Ref. [390]. Two isoscalar bound states, of and of , were found. Similar bound states with the same quantum numbers were found by replacing with . In Ref. [391], possible triple-charm molecular pentaquarks and were further discussed. The and loosely bound molecules are possible.
4.6 A short summary
Since 2003, the observations of a series of charmonium-like near-threshold states have stimulated strong interest in the interactions between two heavy mesons. Later, the coupled-channel effect, various hyperfine interactions and recoil corrections were introduced into the one-boson-exchange model step by step. Within this simple framework, the hidden-charm molecular type pentaquarks were predicted. In 2015, LHCb collaboration did surprise us with two states. Very recently, LHCb collaboration updated their analysis with a ten times larger data sample, which strongly supports the molecular pentaquarks. At present, the one boson exchange model remains an effective and popular tool to deal with the hadronic interactions.
However, this framework sometimes lacks the definite predictive power because of too many unknown parameters such as the coupling constants and cutoff parameter. The original one-boson-exchange model was proposed for the nuclear force where there exists plenty of experimental data such as the deuteron binding energy and enough nucleon nucleon scattering data, which can be used to fix all the unknown model parameters. In contrast, except the pionic couplings, most of the light meson and heavy hadron interaction vertices remain unknown. Especially, the bound state or resonance solution is very sensitive to the cutoff parameter in the form factor which is introduced to suppress the ultraviolet contribution.
One may borrow the same formalism from the modern nuclear force and construct the interaction potential between two heavy hadrons in the framework of chiral perturbation theory. The pion exchange is still responsible for the long-range force. The two-pion exchange mimics the scalar meson and vector meson exchange to some extent. The contact heavy hadron interaction contributes to the short-range interaction in terms of the low energy constants which may be extracted through fitting to the lattice simulation data on the heavy hadron scattering if possible. One may expect significant progress along this direction in the future.
In short summary, many loosely bound molecular states or resonances such as the hidden-charm and open-strange pentaquark states and triple-charm pentaquark states have been proposed either within the one-boson-exchange model or with various unitary schemes, which may be searched for at LHCb and BelleII.
5 Heavy quark symmetry and multiquark states
QCD exhibits the chiral symmetry in the limit when quarks are massless, and the heavy quark symmetry [392] in the limit when quarks have infinitely large masses, both of which play important roles in understanding properties of hadrons and their interactions. The latter symmetry indicates both the heavy quark flavor symmetry (HQSS) which means that the dynamics is not affected by the exchange of heavy quark flavors, and the heavy quark spin symmetry (HQSS) which means that the dynamics is irrelevant of the heavy quark spin. For hadrons containing a single heavy quark, there exist degenerate multiplets which can be classified by the angular momentum of the light degree of freedom inside the hadrons. Recently, discussions about the application of the heavy quark symmetry to hadronic systems containing one heavy quark were performed in Ref. [393].
In the heavy quark limit, the heavy quark core of a hadron looks like a static color source, and the hadron resembles an atomic system. For conventional mesons with one heavy quark and conventional baryons with two heavy quarks, the heavy color source is of course in the representation . They can be related with the heavy antiquark-diquark symmetry (HADS) [394, 395, 396, 397] since the motion and spin interactions of the color source are negligible in this limit. An example of such relations is given in Ref. [397]
[TABLE]
For multiquark hadrons with two heavy quarks, the heavy quark core can also be in the representation . One may naively expect that the lowest states contain the color- heavy diquarks only, and then the doubly heavy tetraquarks can be related to the doubly heavy baryons. In fact, mass relations between , , , and states [164],
[TABLE]
can be found according to the heavy quark symmetry. In Ref. [397], other mass relations were given. Similarly, one may understand that relations exist between pentaquarks and tetraquarks. With these heavy quark relations and heavy quark mass corrections, one may predict relevant multiquark masses.
The heavy quark symmetry may also be used to discuss hadronic molecule problem, where no additional light particles are created and the heavy hadron components remain distinct and move non-relativistically [398]. From the existence of one molecule which consists of two or more heavy hadrons, one can expect possible existence of its partner molecules according to the heavy quark symmetry. The respective binding energies can be obtained by solving the Schrödinger equation.
5.1
In Ref. [399], the exotic charmonium spectrum in hadron-charmonium, tetraquark, and hadronic molecule scenarios were investigated by employing heavy quark spin symmetry. The spectrum patterns from different models turn out to be quite distinct. For the , a lighter partner state around 4140 MeV should exist only within the hadrocharmonium picture.
Three hadronic partner states of the with , , and should exist in the strict heavy-quark limit, if it is a molecular state [400]. Once the one-pion-exchange interaction was included, this result was found to be robust only if all relevant channels and partial waves are considered. The heavy-quark spin symmetry implied spin partners of the and , which were studied by assuming they are and molecules, respectively in Ref. [401]. With the force from short-range contact terms and the one-pion-exchange potential, the authors predicted the existence of an isovector tensor state lying a few MeV below the threshold.
In Ref. [402], the heavy quark symmetry partners of the were discussed by assuming that this exotic meson is a molecule. With the one-boson-exchange () potential model in which the cutoff parameter was extracted from , the bound state problems of these meson-antimeson systems were discussed. The location and quantum numbers of the and were correctly reproduced. Other bound states such as , , and were predicted. We refer interested readers to Ref. [403], which introduced an effective Lagrangian implementing the heavy quark symmetry to describe those molecular states. See also Refs. [404, 405] where the non-relativistic effective theory was applied to describe heavy quarkonium hybrids .
In Ref. [406], Voloshin proposed the existence of light quark spin symmetry (LQSS) in molecular resonances, and , by noticing the weak coupling of the to the channel. In the case of free heavy meson pairs, the spin structure of the two states reads
[TABLE]
In the case of interacting heavy meson pairs, if the interaction in the state differed from that in , the two states would not be the eigenstates for interaction mesons. Then, one has to assume that the interaction between the mesons does not depend on the total spin of the light quark-antiquark pair, which is the proposed LQSS. Unlike the HQSS, this approximate symmetry is not a symmetry from the underlying QCD, and is manifestly broken by the pion exchange force. Although the LQSS is unexpected, it was used to predict the existence of four new isovector negative -parity resonances, , , , and , together with the HQSS. The application of this symmetry to the charm quark case would be complicated, because of the enhanced violation of HQSS and the breaking of isospin symmetry.
5.2
In Ref. [164], Eichten and Quigg calculated the masses of the (, ) tetraquark states based on mass relations obtained from the heavy quark symmetry and finite-mass corrections. According to their results, the state is 102 MeV above the threshold, while the state is 121 MeV below the threshold. The finding is consistent with that by Karliner and Rosner [166]. Besides, this calculation indicated a stable bound state with but no stable .
The chiral Lagrangian that uses HADS to relate doubly heavy tetraquarks to singly heavy baryons was given in Ref. [165], which is an extension of Ref. [407]. The author constructed the chiral Lagrangian for multiplets containing singly heavy mesons and doubly heavy baryons. With these Lagrangians and the measured mass of the , Mehen studied properties of the five lowest-lying excited states together with doubly heavy tetraquarks. He noted that the lowest state should be stable. If its mass is below 10405 MeV, he argued that the state should also be stable.
In Ref. [398], the authors presented a model-independent argument about the existence of near threshold exotic mesons containing two heavy quarks. Based on Born-Oppenheimer and semi-classical considerations, they found the existence of parametrically narrow tetraquarks, which are close to the threshold of two heavy mesons.
5.3 and
In Ref. [408], the structure of heavy quark spin multiplets for the -wave molecular states was studied. The authors introduced the light cloud spin (LCS, the allowed angular momentum of the light degree of freedom) basis, and constructed the unitary transformation matrices in relating the LCS basis and the molecule basis. They found four types of multiplets classified by the structure of heavy quark spin and LCS, and discussed their restrictions on decay widths of the molecules. In Ref. [409], they extended the formalism to study the heavy quark spin multiplet structures of the -like pentaquarks as -wave hadronic molecules. Solving the coupled-channel Schrödinger equation in the OPEP model and comparing results with the previous study [408], they pointed out that it seems difficult to explain the masses and decay widths of the and simultaneously.
In Ref. [410], the authors discussed the possible partner state of the according to the heavy quark spin symmetry assuming it to be a bound state. They predicted the existence of a molecule with a binding energy similar to the . In fact, before the observation of states by LHCb, the authors of Ref. [198] had discussed the heavy-quark spin multiplet structures in hidden-charm molecules.
Besides the hidden-charm partner state of the , the authors of Ref. [410] also discussed the possible open-charm molecule partners of the using HADS. They predicted the molecule with , molecule with , molecule with , and molecule with . The binding energies of all these states are in the 2030 MeV range.
6 QCD sum rules
The formalism of QCD sum rules is a powerful and successful non-perturbative method [411, 412], which has been widely applied to study the mass spectra and decay properties of various exotic hadrons. Since we have thoroughly reviewed its applications to the hidden-charm pentaquark and tetraquark states in Ref. [1] and its applications to the open-charm tetraquark states in Ref. [413], we shall only briefly introduce the recent progress of this method in the present view, and we refer interested readers to Refs. [1, 413, 411, 412, 414, 415, 416, 417] for more discussions.
6.1 A short induction to QCD sum rules
A key idea of QCD sum rules is the quark-hadron duality, i.e., the equivalence of the (integrated) correlation functions at both the hadron level and the quark-gluon level.
When studying the mass spectrum of the hadron , one calculates the two-point correlation function
[TABLE]
at both the hadron and quark-gluon levels. Here is an interpolating current which has the same quantum numbers as the hadron , and their coupling is defined to be
[TABLE]
The interpolating current can partly reflect the internal structure of the hadron , but we still do not fully understand this relation [418, 419]. We shall discuss this in details in the next subsection.
At the hadron level one expresses in the form of the dispersion relation:
[TABLE]
where is the spectral density, for which one usually adopts a parametrization of one pole dominance for the ground state and a continuum contribution:
[TABLE]
At the quark-gluon level, one calculates using the method of operator product expansion (OPE), and calculates the spectral density up to certain order in the expansion. After performing the Borel transformation at both the hadron and quark-gluon levels, and approximating the continuum using above a threshold value , one obtains the sum rule equation
[TABLE]
which can be used to calculate through
[TABLE]
The above two-point correlation function is investigated when extracting hadron masses, while one can also consider the three-point correlation function to study their decay properties:
[TABLE]
Here , , and are the momenta of , , and , respectively. The above correlation function can be used to study the decay. One still calculates it at both the hadron and quark-gluon levels.
At the hadron level one expresses as:
[TABLE]
where is the coupling constant, and , , and are the relevant decay constants.
At the quark-gluon level one calculates using the method of operator product expansion (OPE). Again, by using the quark-hadron duality to relate the expressions of at the hadron and quark-gluon levels, one can calculate the coupling constant , and further evaluate the decay width of the process.
6.2 Interpolating currents and their relations to physical states
When applying the method of QCD sum rules to investigate a physical state, one always needs to construct the relevant interpolating current. However, we still do not fully understand their relations: the interpolating current sees only the quantum numbers of the physical state, so it can also couple to other physical states and thresholds having the same quantum numbers; while one can sometimes construct more than one interpolating currents, all of which couple to the same physical state.
In 2006 we first systematically constructed all the interpolating currents of in a local product form, and found that there are five independent currents for this channel [420]. They can be either in the diquark-antidiquark form ()
[TABLE]
or in the meson-meson form ()
[TABLE]
and can be related to each other through the Fierz transformation [421, 422, 423, 420]. In the above expressions, and are color indices, is the charge conjugation operator, and the superscript represents the transpose of Dirac indices. Note that only five of the ten currents in Eqs. (LABEL:define_meson_current) are independent, and all of them can be written as combinations of Eqs. (LABEL:define_diquark_current). We refer to Ref. [420] for detailed discussion.
Later we applied the same method to systematically study various tetraquark, pentaquark, dibaryon and baryonium states [420, 424, 425, 426, 427, 428, 429, 430, 431, 432, 433, 434, 435, 436, 437, 438, 439, 440, 267]. Since the internal structure of multiquark states is quite complicated, we also applied this method to systemically conventional mesons and baryons [441, 442, 443, 444, 445, 446, 447, 448, 449, 450, 451]. We clearly verified in many cases that the local diquark-antidiquark () and meson-meson () currents can be related to each other, so a straightforward conclusion is that the method of QCD sum rules can not actually differentiate the compact tetraquark and the hadronic molecule, when local four-quark currents are taken into account.
After so many QCD sum rule studies, we have obtained some experiences on hadronic interpolating currents and their relations to hadronic states, for examples,
The singly heavy baryons were systematically studied in Refs. [452, 453, 454, 455, 456, 457, 458, 459, 460], where we found that they have rich internal structures (see references of the review [413] for more relevant discussions). There can be as many as three -wave excited states of , three of , and one of . For each state we constructed one relevant interpolating current. Maybe not all of them exist in nature, but the LHCb experiment [461] did observe as many as five excited states at the same time, all of which are candidates of -wave charmed baryons. 2. 2.
The mass spectra of vector and axial-vector hidden-charm tetraquark states were systematically investigated in Ref. [462], where we found as many as eight () interpolating currents of . Comparably, there were also many vector charmonium-like states observed in particle experiments [35], including the [463], [464], [63], [465], [466], and [467], etc. 3. 3.
The hidden-charm pentaquark states having spin were systematically studied in Refs. [207, 468, 469], where we constructed hundreds of hidden-charm pentaquark interpolating currents. In 2015 the first two hidden-charm pentaquark states, and , were discovered by LHCb [185]. In this year LHCb discovered a new candidate , and at the same time they separated the into two structures, and [68]. One may expect there would exist more hidden-charm pentaquark states to be discovered in future experiments.
The internal structures of (exotic) hadrons are not so simple. In QCD sum rule studies, we can construct relevant interpolating currents to partly reflect such internal structures, but their relations are even more complicated.
To clarify this problem, we found a good subject in Ref. [418] that there are only two independent tetraquark currents of [470]:
[TABLE]
We can also construct four non-vanishing meson-meson interpolating currents with , but they all depend on the above two diquark-antidiquark ones. Their relations can be derived by using the Fierz transformation [421, 422, 423, 420].
We have separately used and to perform QCD sum rule analyses in Ref. [470], by calculating the diagonal terms:
[TABLE]
The masses extracted are both around 2.3 GeV, which can be used to explain the of .
To study the relations between interpolating currents and their relevant physical states, given that we do not know how to diagonalize physical states, in Ref. [418] we further calculated the following off-diagonal term
[TABLE]
which were found to be non-zero, suggesting that and can couple to the same physical state. Based on Eqs. (54) and (55), two new currents were obtained in Ref. [418]
[TABLE]
satisfying
[TABLE]
Hence, these two new currents are non-correlated, and should not strongly couple to the same physical state. In Ref. [418] we assumed they couple to two different states, whose masses were extracted to be
[TABLE]
with the mass splitting
[TABLE]
The mass extracted using is consistent with the experimental mass of the , suggesting that may couple to the ; while the mass extracted using is a bit larger, suggesting that the may have a partner state, whose mass is around MeV larger. The latter may be used to explain the structure in the invariant mass spectrum at around 2.4 GeV [471, 43, 472, 44, 473].
Later this method was applied in Ref. [419] to study the tetraquark states with . Again two independent interpolating currents were found, but any one of them was found to give reliable QCD sum rule results. The mass was extracted to be GeV, which can be used to explain the observed by BESIII recently [474].
In Ref. [475], the authors applied the same method to investigate four charged charmonium-like states at the same time, including the , , , and [476]. Their results suggest that the and are two -wave tetraquark states with , and the other two higher states can be established as their first radial excitations; the and both contains one “good” diquark with and one “bad” diquark with , while the and both contains two “bad” diquarks with ; an illustration is given in Fig. 8.
6.3 and
Recently, the QCD sum rule method has been extensively used to study the tetraquark systems [477, 478, 479, 480, 481, 482, 483, 484, 485, 486, 487, 488, 489, 490, 491, 492, 493, 494, 495, 496, 497, 498, 499, 500, 501, 502, 503]. Most of these works have been already reviewed in Refs. [416, 1]. In this paper, we shall only introduce some very recent studies on the hidden-charm tetraquarks and pentaquarks. We will also introduce some investigations on the doubly, triply and fully charmed and bottom tetraquark states. These tetraquark systems are probably very stable and narrow against the strong decays, and thus motivated lots of theoretical interests. In this subsection we review the QCD sum rule studies on hidden-charm tetraquarks states in recent years.
To study the structures of the and in the diquark-antidiquark model, the authors of Ref. [504] investigated the hidden-charm and hidden-bottom , , , tetraquark states with and in the framework of QCD sum rules. They constructed ten scalar and four tensor interpolating currents to calculate the two-point correlation functions. After performing the numerical analyses, they obtained the mass spectra for the hidden-charm and tetraquark. As shown in Fig. 9, the mass of the hidden-charm tetraquark extracted from the currents , and was about GeV, which was consistent with the mass of the state. However, the tensor tetraquark state was calculated to be around , which was a bit higher than the mass of . This implies that the assignment is favored for than , if it is a tetraquark. The same conclusion can be obtained for . For , their results cannot distinguish the scalar and tensor configurations, due to the nearly degenerated masses for them. It is also shown that the cannot be interpreted as a tetraquark with either or . The was also studied as a tetraquark state in Ref. [505].
In Ref. [267], the newly observed and were studied together with the and states, based on the diquark-antidiquark configuration in the framework of QCD sum rules. The authors found that the and can be both interpreted as -wave tetraquark states with the quark content and quantum numbers : the consists of one -wave “bad” diquark and one -wave “bad” antidiquark, with the antisymmetric color structure ; the consists of similar diquarks, but with the symmetric color structure . These two interpretations were remarkably similar to those obtained in Ref. [462], in which the and can be both interpreted as -wave tetraquark states of , but with distinct color structures. They also calculated the masses of the hidden-bottom partner states of the and and suggested to search for them in the invariant mass distribution. The masses and decay widths for the , , and mesons were also studied as the tetraquark states in Refs. [506, 507, 508] and molecule state in Ref. [509] by using QCD sum rules.
The QCD sum rule has also been performed to study the doubly hidden-charm/bottom tetraquark states in Ref. [510], in which the moment method was adopted instead of the Borel transformation. The moment was defined as the -th derivative of the correlation function in the Euclidean region
[TABLE]
where is the spectral function and the narrow resonance approximation was adopted in the last step. The in Eq. (61) contains the contributions of higher states and the continuum. It tends to zero as goes to infinity for a certain value of . Then the hadron mass of the lowest lying resonance can be obtained as
[TABLE]
In Refs. [510, 511], the authors studied both the and tetraquark systems with various quantum numbers. As shown in Fig. 10, all the tetraquark states were predicted to be slightly below than the mass thresholds of and , while the tetraquarks lied above the two-charmonium thresholds and . Moreover, the positive parity states with were lighter than the negative parity states.
In Refs. [512, 513], the and tetraquark states with were also studied in the Borel sum rule method. The extracted masses for these tetraquarks were slightly different with those obtained in Ref. [510].
6.4 and
The triply heavy tetraquark states with and were also studied in Ref. [514], including the explicitly exotic states and the hidden exotic states , , . The masses are lower than those estimated with the CMI model [107]. Especially, the states were predicted to lie below the bottomonium plus thresholds and their OZI(Okubo-Zweig-Iizuka)-allowed strong decays are kinematically forbidden. These tetraquarks were expected to be very narrow and the authors suggested to search for them in the final states with a meson plus a light meson or photon processes. More studies are needed to understand such triple-heavy tetraquark states.
In Ref. [515], the mass spectra for the doubly charmed/bottom , , tetraquark states with quantum numbers and have been calculated. Their numerical results indicated that the masses of the doubly bottom tetraquarks , , were below the corresponding thresholds of the two-bottom mesons, two-bottom baryons, and one doubly bottom baryon plus one anti-nucleon. These doubly bottom tetraquarks are thus predicted to be stable against the strong decays. These results were consistent with the earlier study in Ref. [516]. Further in Ref. [517], the open-flavor scalar and axial-vector and tetraquark states were also studied and their masses were below the and thresholds respectively, suggesting dominantly weak decay mechanisms.
Later, the doubly heavy tetraquark states have attracted lots of attention [518, 519, 520, 521], motivated by the discovery of doubly charmed baryon . In Ref. [519], the authors calculated both the masses and two-body strong decay width for the and tetraquarks with . They obtained MeV and MeV for the tetraquark, while MeV and MeV for the tetraquark. Since the masses of the doubly bottomed tetraquark were predicted below the two-meson thresholds, the weak decays of the axial-vector state to the scalar state were investigated using the three-point sum rule approach [520]. They recalculated the spectroscopic parameters of the and states in the framework of two-point QCD sum rules. The semileptonic decay channels were then investigated with the width and lifetime MeV and fs. Later in Ref. [521], they also studied the semileptonic decays of the scalar tetraquark to final states and . After reanalysing the spectroscopic parameters of , they used the three-point sum rules to evaluate the weak form factors and extrapolated them to the whole momentum region. The total semileptonic decay width and lifetime of the scalar tetraquark were finally obtained as MeV and ps.
6.5
Since the discovery of the and by LHCb in 2015 [185], there have been lots of theoretical studies using various methods, and the QCD sum rule approach is one of the most popular one.
In Ref. [207], we applied the method of QCD sum rules and studied the and as exotic hidden-charm pentaquarks composed of an anti-charmed meson and a charmed baryon. Our results suggest that the and have quantum numbers and , respectively. We also predicted the masses of their hidden-bottom partners to be GeV and GeV. This study was later expanded in Refs. [468, 469], where we found two mixing currents, and again our results suggest that they can be identified as hidden-charm pentaquark states having and , while there still exist other possible spin-parity assignments, such as and , which still needs to be clarified in further theoretical and experimental studies.
Especially, in Ref. [468] we systematically constructed all the possible local hidden-charm pentaquark currents with spin and quark contents , through which we found that the internal structure of hidden-charm pentaquark states is quite complicated. We derived some mass predictions, and we summarized in Table 5 all the mass predictions that are extracted from single currents and less than 4.5 GeV. The coincidence of these mass predictions with the masses of , , and measured in the recent LHCb experiment [68] drives us to the “molecular” picture [379] that the can be well explained as the bound state with ; the and can be explained as the and bound states with , respectively, while one of them may also be explained as the bound state with or the bound state with ; there is still a place for the , i.e., to be explained as the bound state with . The bound state with was not investigated in Ref. [468]. See Ref. [379] for detailed discussions of these possible interpretations.
In Ref. [522], Azizi et al performed QCD sum rules analyses on the hidden-charm pentaquark states with and . They adopted a molecular picture for states and a mixed current in a molecular form for , and their analyses suggested that the and can be considered as hidden-charm pentaquark states with and , respectively. Later in Ref. [523], they applied the same method to study the hidden-bottom pentaquark states with spin and , and predicted their masses. The strong decays and electromagnetic multipole moments of the were later studied in Refs. [524, 525, 526].
In Ref. [527], Wang constructed the diquark-diquark-antiquark type interpolating currents, and his QCD sum rule studies also supported to assign the and as and pentaquark states, respectively. This study was later expanded in Refs. [528, 529, 530].
7 Three-body system
There are various two-body interpretations of the exotic hadrons, as we have reviewed in previous sections. There are also some three-body interpretations for exotic hadrons. For example, combining the Faddeev equations with the chiral unitary model [531, 532], the authors interpreted the as a dynamically generated state in the system. In this section, we review some of these studies. We note that the three-body system have also been extensively studied within the chiral effective field theory [533, 534, 535, 536, 537] and by using Lattice QCD [538, 539, 540, 541, 542, 543, 544, 545, 546, 547, 548, 549, 550, 551, 552, 540, 553, 554] in recent years, but we shall not review them in the present paper, and just recommend interested readers to the above references for detailed discussions.
7.1 A short introduction to Faddeev equations
When studying a three-body system, usually one does not consider the “direct” three-body interactions, as shown in Fig. 11. Instead, one needs to consider all the possible two-body exchanges/interactions in a system of three particles. One can use the Faddeev equations, named after their inventor Ludvig Faddeev [555]. We refer to Ref. [556] for detailed instructions.
In a three-body system there are three different two-body subsystems. The idea of the Faddeev equations is to sum up all the pair forces in each two-body subsystem to infinite order among all three particles. First we write the Schrödinger equation for a three-body system as
[TABLE]
where is the kinetic energy of the relative motion for three particles; , , and are the three interactions in the two-body subsystems. Then we transfer Eq. (63) into its integral form:
[TABLE]
and further iterate it many times to be:
[TABLE]
This can be shown graphically in Fig. 12, where all the pair forces in each two-body subsystem have been summed up iteratedly to infinite order.
To solve Eq. (65), we decompose into three Faddeev components:
[TABLE]
where is that part of having as the last interaction to the left. As an example, we show graphically in Fig. 13.
We can further iterate to be:
[TABLE]
so that
[TABLE]
The kernel can be simplified to be
[TABLE]
where sums up to infinite order.
Finally, we arrive at three coupled equations
[TABLE]
which are called as the Faddeev equations. For identical particles, these equations can be simplified. Again, we refer Refs. [555, 556] for detailed discussions.
7.2 Applications of Faddeev equations to exotic hadrons
The Faddeev equations have been widely applied to study the quantum-mechanical three-body problem, such as three-nucleon bound states ( [557], [558] and [559]) and quasi-bound state [560] in nuclear physics. Here we review its applications to exotic hadrons in hadron physics.
In Refs. [531, 532] the authors developed a method to study the meson-meson-baryon system by solving the Faddeev equations together with the chiral unitary model. Here we briefly review its application to the [532]. There the was interpreted as a resonant state, with the dynamically generated in the channel. Hence, to study the , the authors solved the Faddeev equations in the three-body system.
The full three-body scattering matrix is
[TABLE]
where is the Faddeev partition. Note that and are used in the previous subsection. They further write as
[TABLE]
where is the two-body scattering matrix where the particle is a spectator; the term corresponds to a “disconnected” diagram and should be removed; sums up all the diagrams with the last two matrices being and :
[TABLE]
In the above expression is the three-body propagator:
[TABLE]
All the terms are written as t^{i}~{}G^{ijk}\left(t^{j}g^{jk}t^{k}\right)\Big{|}_{\rm on-shell}:
[TABLE]
where
[TABLE]
The detailed explanations of the above equations can be found in Refs. [531, 532]. As examples, we show and graphically in Fig. 14.
The two-body scattering matrices for the interactions of pseudoscalar and pseudoscalar mesons can be found in Ref. [561], while those for the interactions of pseudoscalar and vector mesons can be found in Ref. [562]. Using them as inputs, the authors solved the Faddeev equations, Eqs. (7.2), in the chiral unitary approach. They found a clear sharp peak of around 2150 MeV with a narrow width about tens of MeV, as shown in Fig. 15. This structure can be associated to the , so their results indicate that the is a dynamically generated resonant state in the system. We refer to Ref. [532] for more discussions, where the coupled channel has also been taken into account.
The above method has been systematically applied to study the meson-meson-baryon (two pseudoscalar mesons and one octet baryon) system, and many dynamically generated resonant states are proposed. For examples, the is interpreted as a dynamically generated resonance in the system with and [563]. The three mesons system has also been systematically studied within the same approach in Refs. [564, 565]. For examples, the is interpreted as a dynamically generated resonance in the system [564]; an exotic state with mass around 4140 MeV and was predicted in the system [565].
A similar approach was applied to study the system by solving the Faddeev equations with relativistic kinematics in Ref. [566], and a resonance with and was predicted to be around MeV. The -wave three-body system is investigated in Ref. [567], and the results suggest that there exist -wave bound states with quantum numbers , all of which locate at 30-40 MeV below the threshold. The systems are investigated in Ref. [568], and a unique bound state of three B mesons with and was predicted to be around 90 MeV below the threshold of three -mesons.
7.3 Fixed center approximation
The Faddeev equations are sometimes not easy to be solved, so the fixed center approximation was proposed to deal with them [569, 570, 571, 572]. In Ref. [573] this approach within the chiral unitary model was applied to study the , , , and resonances as multi- states, and later on to many other exotic structures. In this paper we briefly review its application to the three-body system [574].
In Ref. [574] the authors investigated the three-body system by considering that two of the mesons are clusterized forming an strongly binding resonance, i.e., -. In such a situation one can use the fixed center approximation to the Faddeev equations. We show this diagrammatically in Fig. 16, where the external particle interacts successively with the other two mesons forming the . The total scattering matrix is
[TABLE]
where () are the two partition functions, accounting for all the diagrams starting with the interaction of and the -th meson of the compound system:
[TABLE]
Here is the scattering amplitude, which has been systematically investigated in Ref. [575] and the results therein can be used here as inputs; is the loop function for the particle propagating inside the compound system
[TABLE]
where the form factor is
[TABLE]
and the normalization factor is
[TABLE]
We recommend interested readers to Refs. [573, 574] for detailed explanations of the above equations, where the wave function of the as well as its width are discussed.
Because there are two identical mesons, one has , and Eqs. (77), (78) and (79) can be simplified to be:
[TABLE]
The results are shown in the left panel of Fig. 17, where there is a peak at around MeV with a width of MeV. This peak does not have a standard Breit-Wigner form, so it could be due to the cusp effect.
In Ref. [574] the authors also applied the same approach to investigate the three-body system. The results are shown in the right panel of Fig. 17, where there is a peak around 2372 MeV with an approximate Breit-Wigner shape. This peak can be associated with the of [576].
The fixed center approximation to the Faddeev equations has been systematically applied to study the -multi- systems [577, 578]. For examples, in Ref. [577] the , , and are interpreted as molecules made of an increasing number of and one mesons; in Ref. [578] the -multi- system was investigated and several charmed resonances, , , and , were predicted to be around 2800-2850 MeV, 3075-3200 MeV, 3360-3375 MeV and 3775 MeV, respectively. The above method has also been applied to study the system [579], the system [580], the system [581], the system [582], the system [583], and the system [584], etc.
8 The Skyrme model and the chiral quark-soliton model
The Skyrme model was proposed by T. H. R. Skyrme in 1961 by introducing the Skyrme term to the nonlinear sigma model [585]. Within this model the baryons appear as collective excitations of the meson fields. In fact, the non-trivial topological field configuration, which is called as soliton, is identified with the baryon. Later in 1979, E. Witten developed this model by arguing that baryons indeed emerge as solitons in the large generalization of QCD [586]. He also studied static properties of nucleons in the Skyrme model with two light flavors [587]. The Skyrme model was widely applied in hadron physics. Based on the Skyrme model one might qualitatively understand the EMC observation that the quark spin contribution to the total nucleon spin was unexpectedly small [588]. This feature is intensively studied in the Skyrme model with three light flavors in order to understand the nucleon’s strangeness content [589].
In the large -limit the chiral quark-soliton model (QSM) [590, 591, 592, 593]. is quite similar to the Skyrme model, and sometimes they are taken as the same method/approach. In the Skyrme model one has integrated out the quark fields, and so directly operates with the action involving the pseudoscalar mesons without the quark fields. Actually, one can obtain the effective action of the Skyrme type from the one of the chiral quark-soliton model by a gradient expansion.
The Skyrme model and the chiral quark-soliton model with three light flavors can be naturally applied to study exotic baryons, which emerge as rotational excitations of nucleons [45]. In this paper we briefly review some of these applications, separately in the following subsections. We note that there are lots of investigations within this approach, and we refer interested readers to the lecture note by Diakonov [594], the one by Weigel [595], and the one by Ma and Harada [596] for detailed instructions.
8.1 A short introduction to the Skyrme model
The Lagrangian of the nonlinear sigma model is
[TABLE]
Here is the decay constant of the Nambu-Goldstone boson , and the field is defined as
[TABLE]
where is the meson field, and the field variables are related to through
[TABLE]
Under the chiral transformation, the field transforms as , so the Lagrangian (84) is chirally invariant. With the nonlinear realization of the chiral symmetry, the Hamiltonian density of the nonlinear sigma model can be obtained as
[TABLE]
and its energy can be expressed as
[TABLE]
Here the first term is the rotational energy, and the second term is the static energy. The static energy is not stable: under the rescaling of the space coordinates , one has .
To avoid the above stability problem of the static energy, Skyrme introduced the so-called Skyrme term, and extended the nonlinear sigma model Lagrangian into the Skyrme model:
[TABLE]
where is a dimensionless parameter, indicating the magnitude of the soliton. Its energy can be obtained as
[TABLE]
One can solve the following extremum stable condition
[TABLE]
and the obtained stabilized solutions are called as Skyrme solitons or skyrmions.
Based on Eq. (89), the Witten-Wess-Zumino term was added at the level of the action to reproduce QCD’s low energy anomalous structure [586, 597, 598]:
[TABLE]
where is the number of colors. Sometimes the symmetry breaking terms are also necessary:
[TABLE]
The Skyrme model can be used to derive a series of states [599, 600, 601, 602]. In the three-flavor case, the lowest states are the and , which can be used to describe the ground-state octet baryons of and decuplet baryons of . The remaining states are all exotic baryons, and the states at the next level are the , and .
In the chiral symmetry limit one can write the effective Hamiltonian as:
[TABLE]
where denotes the baryon spin, denotes the Casimir operator for the representation , and the other three parameters and are treated as free parameters. In the Skyrme model the symmetry breaking Hamiltonian is
[TABLE]
which leads to both the first order correction to the baryon mass
[TABLE]
and the second order correction
[TABLE]
We refer interested readers to Refs. [603, 604] for detailed discussions on the above equations. Especially, the mass of the pentaquark belonging to was predicted to be [603, 604]
[TABLE]
However, this mass value is quite sensitive to the choice of parameters. For example, the authors of Ref. [604] plot in Fig. 18 the results of the constrained fits, and one can find a rather steep rise of with . The above mass value 1530 MeV is very close to the mass of the observed by LEPS [46], whose observation quickly led to renewed interest in the Skyrme model descriptions of baryons. Unfortunately, the was not confirmed in a series of high precise particle experiments [48].
The Skyrme model has many advantages in describing hadron physics, and we refer to Refs. [605, 595, 596] for detailed discussions. Especially, the Skyrme model can be applied to study both baryon and meson physics in free space [587], while it can also be applied to study nuclear matter and the medium modified hadron properties [606, 607, 608]. In recent years, the Skyrme model is extended so that a) the effects of the higher resonances and higher chiral order terms can be self-consistently analysed [609, 610, 611, 612], and b) properties of heavy baryons containing a heavy quark can be systematically investigated [613, 614, 615, 616, 617].
8.2 A short introduction to the chiral quark-soliton model
In the previous subsection we have briefly introduced the Skyrme model. In this subsection we shall briefly introduce the chiral quark-soliton model [618, 619]. Both of them have similar group structure, and can be applied to study the pentaquark belonging to the exotic flavor representation. In fact, the prediction of the narrow width of the in the chiral quark-soliton model [45] stimulated the experimental search by LEPS [46]. We refer to Refs. [594, 620] for detailed discussions.
According to the chiral symmetry, the interaction of pseudoscalar mesons with constituent quarks can be written as:
[TABLE]
where the fields are the ordinary pseudoscalar mesons , , and . In the large limit baryons emerge as solitons of this chiral action [586]. The model of baryons based on the large limit is called the chiral quark-soliton model [620].
Actually, we can also use the word “soliton” for the self-consistent pion mean field in the nucleon. Because these Hartree-Fock states composed of valence- and sea-quarks are degenerate with respect to rotations in both space and flavor-space, we can use the following collective Hamiltonian to describe them in the rigid rotor approximation [621, 622]:
[TABLE]
where is the classical mean-field energy of the quark system, is the baryon spin, and is the Casimir operator for the representation . This Hamiltonian is very similar to the one of the Skyrme model given in Eq. (94) except a different symmetry breaking term, i.e.,
[TABLE]
where , , and are parameters of order , are Wigner rotation matrices, and are collective spin operators.
The symmetry-breaking term mixes different representations, so the collective wave functions can be written as the following linear combinations:
[TABLE]
where are the states which reduce to the representation when and the spin index is suppressed. , , and are the -dependent coefficients. For example, we show in Fig. 19 the weight diagrams for the lowest-lying baryon multiplets, which contains the baryons with the hypercharge .
Assuming the state to be a member of the flavor representation, Diakonov, Petrov and Polyakov investigated an exotic baryon with spin 1/2, isospin 0 and strangeness +1, and predicted its mass to be about 1530 MeV and its total width to be less than 15 MeV [45]. Now this state was known as the , which was first observed by LEPS [46]. However, we note that this small width actually resulted from an arithmetic error, as explained by Jaffe in Ref. [624], and correcting this error gives a width twice as large at least. See Ref. [589] for detailed discussions.
Later in Refs. [625, 626] the masses of the lowest-lying baryons were calculated by including the effects of isospin symmetry breaking arising from both the quark masses and electromagnetic self-energies. The authors used the following collective Hamiltonian in the chiral soliton model [621, 622]:
[TABLE]
where is the rotational energy, and are the symmetry-breaking terms including both isospin and flavor symmetry breaking effects:
[TABLE]
We refer interested readers to Refs. [625, 626] for detailed discussions on the above equations.
The baryon mass splittings due to the isospin and flavor symmetry breaking up to the first order were investigated in Ref. [625], and those due to the second order flavor symmetry breaking were investigated in Ref. [626]. The authors estimated the mass of the as a member of the baryon antidecuplet to be MeV, which can be used to explain the , a narrow bump-like structure observed by GRAAL in the quasi-free cross section in 2006 [627]. The authors also estimated the pion-nucleon sigma term to be MeV.
In 2015, another narrow resonant structure was observed by the GRAAL Collaboration in real Compton scattering off the proton. Its mass and width were determined to be GeV and MeV, respectively. In the chiral quark-soliton model, this structure can be identified as a member of the eikosiheptaplet () with spin 3/2 [628]. Later in Ref. [623] the authors investigated the strong and radiative decay widths of the narrow nucleon resonances and together within the chiral quark-soliton model [623]. Especially, they found , indicating that the production of is more likely to be observed in the proton channel.
9 Progresses from Lattice QCD
Lattice QCD is the unique non-perturbative theoretical framework to study the hadron spectroscopy starting from the first principle QCD Lagrangian. In a lattice QCD simulation, the hadron mass can be extracted from the time-dependence correlation function
[TABLE]
where the interpolating operator creates the state of interest from the vacuum at time and annihilates the state at a later Euclidean time . is the eigenvalue of the Hamiltonian for the system. These creation and annihilation operators should have the same quantum numbers, which are called interpolators. In principle they can couple to all physical eigenstates with the same given quantum numbers but different magnitudes of overlaps
[TABLE]
In lattice QCD, the correlation functions are usually evaluated by using a large basis of interpolators constructing of quark and gluon fields, the Dirac gamma matrices and the gauge-covariant derivative operator. To search for the tetraquark and pentaquark signatures, the two-particle interpolators are usually adopted to deduce the two-point correlation functions. If the interpolator basis is complete enough, one can extract and identify all two-particle discrete energy levels with the given quantum numbers from the full correlation functions. The exotic multiquark states are related to the extra energy levels in addition to the expected non-interacting two-particle states. The two-particle energies are obtained from the correlation function by solving the generalized eigenvalue problem (GEVP) [629, 630]:
[TABLE]
where is an eigenvector and a suitably chosen reference time . The eigenvalues feature the relevant state energies [629]
[TABLE]
Fitting this exponential decay behavior of the time dependence for each of these quantities, the exact discrete spectrum can be extracted from the effective mass plateau of the eigenvalue [631, 632, 633, 634, 635]
[TABLE]
In the infinite volume, a bound state of two particles can be defined as a discrete energy eigenstate of the Hamiltonian with energy level below the two-particle threshold within the Lüscher’s method [636, 637, 629, 638, 639]. In a two-particle scattering process, the resonant character on the lattice can be observed as the negative center of mass (c.m.) momentum , which corresponds to the attractive interaction between the two particles and renders the energy level lower than the non-interacting threshold. The Lüscher’s formula provides a direct relation of and the elastic scattering phase shift in the infinite volume, e. g. in the case of the -wave elastic scattering
[TABLE]
where is the zeta-function. In the limit of , the phase shift will approximate if there exists a true bound state at the particular energy in the infinite volume. Therefore, one usually needs to study the quantities and to really identify a bound state in a lattice simulation. We refer to a recent review article [640] for detailed discussion, and refer to Refs. [641, 642] from the viewpoints of the chiral perturbation theory.
In the last decade, there are some remarkable progresses in the study of new hadron states by lattice QCD. In this section, we try to introduce these investigations, such as the studies on the , , the charged states, the doubly heavy tetraquark states and the hidden-charm pentaquark states.
9.1
The underlying structure of the meson was studied by TWQCD Collaboration [643]. They calculated the mass spectra of the hidden-charm hybrid mesons (), molecules and diquark-antidiquark tetraquarks with , in quenched lattice QCD with exact chiral symmetry. The authors computed the time correlation functions for all operators in the above various schemes and detected all these exotic resonances. However, the lowest-lying hybrid charmonium lies about 250 MeV above the mass of . In contrast, the hidden-charm molecule () and tetraquark () states were detected with masses 4238(31)(57) MeV and 4267(68)(83) MeV respectively, which are in good agreement with the mass of . Finally, the authors suggested the to be a molecule since it has a better overlap with the molecular operator than any other ones. Besides, they also detected resonances of molecular and tetraquark states with quark contents around MeV, and with quark contents around MeV.
Before the observation of , the mass spectrum of the charmonium hybrids were predicted by the UKQCD quenched lattice NRQCD [644, 645, 646] and anisotropic lattice simulations [647]. In Ref. [648], the Hadron Spectrum Collaboration performed dynamical lattice QCD calculations to study the mass spectra of highly excited charmonium and charmonium hybrid mesons. To provide a large basis of operators with a range of spatial structures in each channel, the authors used the derivative-based construction for operators of the general form to calculate the two-point correlation functions. Their computations were presented for up and down quarks corresponding to a pion mass of 400 MeV at a single lattice spacing. They calculated the dynamical spectrum of charmonium hybrids with exhaustive quantum numbers including the exotic ones () up to 4.5 GeV. Their results identified the scheme that the hybrid meson can be interpreted as a colour-octet quark-antiquark pair coupled to a chromomagnetic gluonic excitation. In this scheme, the lightest hybrid supermultiplet consists of negative-parity states , in which the colour-octet quark-antiquark pair is in S-wave. Such a pattern suggested that the lightest gluonic excitation has an energy scale of 1.2-1.3 GeV, compared with the conventional charmonium bound states. The first excited hybrid supermultiplet composed of a P-wave colour-octet quark-antiquark pair coupled to an excited gluonic field would contain positive-parity states with . The vector hybrid charmonium with was extracted at GeV, which agrees well with the mass of and supports the hybrid interpretation of this state. The charmonium hybrid interpretation for was also performed in Ref. [649], in which the authors calculated the masses of the ground and first excited states for the hybrid mesons with and . They finally concluded that the meson can not be identified as the ground state of the hybrid, but most probably the first excited state of the hybrid charmonium state. However, the QCD sum rule calculation gave much lower mass prediction for the vector hybrid charmonium, which was around 3.4 GeV [650]. The heavy quarkonium hybrid states with various flavor structures and quantum numbers have been extensively studied in Refs. [650, 651, 495].
Since the lattice investigations are always performed with unphysical-heavy light quarks, it is significant to explore the changes in the pattern of states with light-quark masses closer to their physical values [652]. For the charmonium systems, the light-quark mass dependence enters through the sea quark content in the dynamical gauge field ensembles. In the case of the open-charm charmed mesons, the existence of a valence light quark affects the spectrum directly. To investigate the light quark dependence effect on lattice QCD simulation, the Hadron Spectrum Collaboration presented spectra of the excited charmonium, open-charm and mesons from dynamical lattice QCD calculations with a pion mass MeV in Ref. [653]. The charmonium spectra labelled by were shown in Fig. 20 together with their previous results with MeV in Ref. [648]. The masses of the low-lying states are generally consistent between these two ensembles within statistical uncertainties. Only the hyperfine splitting has a small but statistically significant increase when the pion mass decreased. For the higher states in the spectrum, the masses are generally larger with the MeV ensemble. Especially for the hybrids, there is a small but significant increase in their masses as the pion mass is reduced, which will increase the mass splitting between the hybrids and low-lying conventional charmonium states. However, the statistical uncertainties at higher energies are larger. The unstable nature of these states above threshold may be important [654, 655, 656]. The authors thus concluded that the hybrid mesons appear to show a mild increase in mass as the light quark mass is decreased, but the pattern of states and supermultiplet structure are unchanged.
In Ref. [657], Chen et al. investigated the existence of the exotic vector charmonia by constructing the interpolating operators , where is the spatial index and the chromomagnetic field tensor. Compared to the usual hybrid operator, this is a new type of the hybrid-like operator that the charm quark-antiquark pair component and the gluon field are split into two parts by an explicit spatial displacement . Such a configuration can resemble the center-of-mass motion of the recoiling against an additional degree of freedom, and thus is expected to suppress the coupling to the conventional charmonia. They accordingly calculated the two-point correlation functions by using these operators. To eliminate the contribution from the conventional charmonium states, they combined linearly the correlation functions at two specific different as , where is a tunable parameter. For the numerical analyses, they used the tadpole-improved gauge action [658, 659, 660] to generate gauge configurations on anisotropic lattices. Two lattices and with different lattice spacings were used to check the discretization artifacts. After fitting the effective mass plateaus as shown in Fig. 21, they observed a vector charmonium-like state with a mass of GeV with exotic nature. In addition, they calculated the leptonic decay constant of this exotic state to be MeV, which accordingly led to a very small leptonic decay width eV. The mass and leptonic decay width of this signal are consistent with the production and decay properties of .
9.2
The TWQCD Collaboration also studied the state as a hidden-charm four-quark meson with in quenched lattice QCD with exact chiral symmetry [661]. They used both molecular operator and tetraquark operator to compute the time-correlation functions for two lattice volumes and . For both the molecular and tetraquark operators, they detected a resonance with mass around MeV in the limit , which is identified as . Their results showed that has good overlap with the molecular operator as well as the tetraquark operator in the quenched approximation. Comparing with their previous study for , they concluded that may be more tightly bound than , which had better overlap with the molecular operator than any tetraquark ones [643]. Since only one energy level was extracted near the threshold, the result could not support or disfavor the existence of . Letting , they predicted a heavier exotic meson around MeV. By using the same method and interpolating currents as above, they investigated the mass spectrum of the tetraquark states with . They also detected an axial-vector resonance around MeV [662]. In these investigations of TWQCD Collaboration, the authors neglected the disconnected quark loop diagrams for annihilation channels. However, these diagrams provide the mixing effects between states created by and operators, which have been proven to be important for the lattice studies of light scalar tetraquarks [663, 664], and the string breaking in the static limit [665].
As indicated above, the mass spectra of the highly excited charmonium states have been calculated up to around 4.5 GeV in dynamical lattice QCD by the Hadron Spectrum Collaboration in Ref. [648]. The first radial excitation of the P-wave charmonium state is extracted about 110 MeV above the mass of , where the mass of is subtracted from the calculated mass in order to reduce the systematic error from the tuning of the bare charm quark mass. However, the results presented in Ref. [648] are for unphysically-heavy up and down quarks corresponding to MeV. The dynamical lattice simulation with lighter up and down quarks MeV showed that only a very small mass increase for the low-lying and even higher charmonium states [653]. Moreover, several other lattice simulations found one excited charmonium state near the mass of [666, 667, 668]. However, none of these simulations can unambiguously determine whether this state is or the scattering state.
In Ref. [669], the authors calculated the charmonium spectrum including higher spin and gluonic excitations. In particular, they studied the mixing of charmonia with states created by the hidden-charm molecular operators by using the variational generalized eigenvalue method as well as improved stochastic all-to-all propagator methods. Their simulations were employed with light pseudoscalar mass down to MeV and sea quarks. They investigated the binding between pairs of and anti- mesons in the pseudoscalar, vector and axialvector sectors. Only the axialvector channel was clearly attractive. Moreover, the mixing effects between the isoscalar radial charmonium states and molecular states were very small for the pseudoscalar and vector channels. But for the channel, they found a significant binding MeV. This binding is much bigger than the mass difference .
For the molecular interpretation of , many dynamical studies supported the attractive channel of and mesons. One can consult Refs. [1, 72] and references therein for these discussions. Due to the coupled channel effects, the -wave scattering state with can easily couple to the state [670, 671] and significantly lower the mass of the pure state predicted in the GI model [33].
In lattice QCD, this picture was supported by the dynamical lattice simulation in Ref. [672], which was based on one ensemble of Clover-Wilson dynamical configurations with MeV. The authors identified the low-lying and as well as the nearby discrete scattering levels and for both and channels. They counted the number of lattice states near the threshold in order to establish the existence of the . In their simulation, they chose the interpolating fields that couple to as well as the scattering states, i.e., , , (for ), and (for ). Finally, they found a candidate for the charmoniumlike state MeV below the with and , in addition to the nearby and discrete scattering states, as shown in Fig. 22. They extracted large and negative scattering length, fm, and the effective range, fm. Moreover, they found the state but no candidate for the in the channel if they used only five scattering interpolators but no interpolating fields. In the channel, they did not find a candidate for , which is consistent with the assignment of the state in PDG [35].
Later, the isoscalar charmonium state was also investigated by the Fermilab Lattice and MILC Collaborations on gauge field configurations with flavors of highly improved staggered sea quarks (HISQ) with clover (Fermilab interpretation) charm quarks and HISQ light valence quarks [673]. They used a combination of the and interpolating operators. They found an isosinglet candidate of with an energy level 13(6) MeV below the threshold, which agrees with the result in Ref. [672].
These studies were then extended by Padmanath, Lang, and Prelovsek by utilizing a large basis of the , two-meson and diquark-antidiquark interpolating fields (altogether 22 interpolators) [674]. The Wick contractions were considered including both the connected contraction diagrams and the diagrams for annihilation channels in which the light quarks do not propagate from source to sink. Their simulation was performed with and MeV, which aims at the possible charmonium-like signatures for both and channels. They found a lattice candidate for with close to the experimental data only if both the and interpolators are included. However, they found no candidate if the diquark-antidiquark and interpolators were used in the absence of , which indicates that the Fock component is crucial for the while the four-quark component alone does not produce the signal. Their simulation also implies a combined dominance of the and operators in determining the position of the energy levels, while the tetraquark operators do not significantly affect the results. Moreover, their results do not support the neutral or charged in the channel. And no signature of the exotic state (candidate for the ) with are found below 4.2 GeV. We refer readers to Refs. [675, 676, 677, 678] for the analytic investigations considering the quark mass dependence, the volume dependence and the effect from the isospin breaking.
The unphysical pion mass MeV on the lattice is still much larger than its physical mass 140 MeV. Within the molecular scheme, the long range one-pion-exchange force plays a dominant role in the formation of the loosely bound molecular state [679, 680], which decays exponentially as the pion mass increases. The present lattice simulation with the pion mass MeV was obtained from a rather small lattice volume, which has already indicated the important role of both the and components in the formation of the signal on the lattice. As the pion mass approaches the physical value 140 MeV, one shall expect a larger component within the physical state.
9.3 The charged states
The observations of the charged charmonium-like states and trigged lots of investigations of the molecular interpretations for their inner structures [680, 681, 682, 683]. In Ref. [684], Prelovsek and Leskovec searched for the in the and channel in their lattice simulation with degenerate dynamical quarks and MeV. They used six meson-meson type of interpolators and . They calculated the time-dependence of the correlation matrix by considering the connected diagrams of Wick contraction. They omitted the diagrams involving charm annihilation, which effect was proven to be suppressed due to the OZI rule for the conventional charmonium [685]. They observed four energy levels which almost exactly coincide with the energies of four non-interacting and scattering states. They did not find any additional energy level relating to the state. They suggested to perform simulations including meson-meson interpolators as well as some different type of interpolators, such as diquark-antidiquark interpolators with the same quantum numbers.
Such a simulation was performed in Ref. [686] by using 18 meson-meson interpolating operators and 4 diquark-antidiquark interpolators with structure in the channel. The two-meson operators included interpolators such as , , and . They aimed to extract and identify all 13 two-meson energy levels from the corresponding correlation functions and examine the possible additional states related to the exotic hadron. As a result, they found all expected lowest thirteen levels, appearing near the corresponding non-interacting energies of the two-particle states. A two-particle level will disappear from the spectrum or become very noisy when the corresponding operator is absent from the correlator matrix. However, the energy spectrum was unaffected after the operators were omitted. Besides these two-particle levels, they did not find any additional state below 4.2 GeV relating to an exotic candidate. They listed some possible reasons for the absence of the candidate: a) the experimental resonance might not be of dynamical origin since it has not been seen in the meson decay process; b) it might be a couple-channel threshold effect [687, 688]; c) their simulation was performed at a large unphysical MeV; d) their interpolator basis may not be complete enough to render a state in addition to thirteen two-meson states; e) the S-D mixing effect may be important for creating the experimental states. In Ref. [673], the authors used the Highly Improved Staggered Quark action to search for the with a combination of the and channels. Again, they found no evidence for the existence of . A similar analysis was performed in Ref. [689] by considering a basis of five interpolating operators and three operators with the same . Their results didn’t show any unknown energy level.
These negative results of searching for the exotic charged in lattice QCD simulations were also supported by the investigations of the resonance-like structures in the low-energy scattering of the systems by the CLQCD Collaboration [690, 691]. In Ref. [690], Chen et al. presented an exploratory lattice study of the low-energy scattering of the two-meson system by using single-channel Lüscher’s finite-size technique. The calculation was based on twisted mass fermion configuration of size with a lattice spacing of about fm with three pion mass values MeV, 420 MeV and 485 MeV. They constructed the -wave two-meson operators with to calculate the correlation functions. Since the energy being considered is very close to the threshold, they computed the scattering length and effective range to search for the evidence of the existence of . They found negative values of the scattering length for all three pion masses, which indicated a weak repulsive interaction between the two mesons and . They also checked the possibility of the bound state for the negative energy shifts. None of them was consistent with the signal of a bound state. Based on their results, they concluded that their lattice QCD simulation did not support a bound state in the channel corresponding to , at least for the pion mass values being studied. Later in Ref. [691], they performed a similar investigation for the low-energy scattering of the two-meson system in the same channel with . Their results gave negative values of the scattering lengths for all three pion masses, which again indicated a weak repulsive interaction between the two vector charmed mesons, and did not support a bound state of the two mesons in the channel corresponding to the state.
The lattice QCD studies of the spectrum and the low-energy scattering parameters with the standard Lüscher’s method gave no candidates for the and in all , , and two-meson channels. These results indicated that the and may not be the conventional resonance states. In Refs. [692, 693], the HAL QCD Collaboration studied the coupled-channel interactions from -flavor full QCD simulations at three pion masses MeV, 570 MeV, 701 MeV in order to explore the structure of . To consider the interactions among the , and channels, the authors calculated the matrix from the equal-time Nambu-Bethe-Salpeter (NBS) wave functions by using the coupled-channel HAL QCD method [694, 695, 696]. They extracted the -wave coupled-channel potential and found that all diagonal potentials , and were very weak, which indicated that the was neither a simple hadrocharmonium nor molecule state. Instead, they found the coupling and the coupling from the off-diagonal potentials were both strong, which implied that the structure of can be explained as a threshold cusp. They further calculated the invariant mass spectra and pole positions associated with the coupled-channel two-body matrix on the basis of . Their results supported the resonance as a threshold cusp induced by the strong coupling. They further made a semiphenomenological analysis of the three-body decays processes. As shown in Fig. 23, they found that the coupled-channel potential can well reproduce the experimental peak structures in the and decays. When the off-diagonal components of were turned off, the peak structures at 3.9 GeV disappeared in both cases. Accordingly, the authors concluded that the was not a conventional resonance but a threshold cusp [692, 693].
As indicated in Refs. [684, 686], including a large basis of four-quark interpolating operators in the calculations could lead to more reliable determination of finite-volume mass spectra and scattering amplitudes. In the very recent lattice QCD study [697], the Hadron Spectrum Collaboration constructed 29 hidden-charm interpolating currents and 48 doubly-charmed interpolating currents with different color-flavor-spatial-spin structures. Using these diverse bases of meson-meson operators as well as compact tetraquark operators, they computed the two-point correlation functions in the isospin-1 hidden-charm and doubly-charmed sectors using the distillation framework [698]. Their calculations were performed on an anisotropic lattice volume using 478 configurations and a Clover fermion action with flavors of dynamical quarks. The mass of the two degenerate light quarks corresponds to MeV while the strange quark was tuned so that its mass approximated the physical value [699]. The multiple energy levels associated with the meson-meson levels which could be degenerate in the non-interacting limit were extracted reliably. In Fig. 24, the mass spectra in the isospin-1 hidden-charm sector were shown for the lattice irreps corresponding to the and for the charge neutral channels, respectively. One found that every energy level had a dominant overlap onto one meson-meson operator, while the addition of tetraquark operators to a basis of meson-meson operators did not significantly affect the finite-volume spectrum. For all these channels, the number of energy levels in the computed spectra was equal to the number of non-interacting two-meson levels expected in the considered energy region and they all lie close to the non-interacting levels. The authors concluded that there were no strong indications for any bound state or narrow resonance in these channels [697]. Especially for the isovector irreps with , they didn’t find any candidate for the and states. Since the has been admitted as an isosinglet [35], this study gave no hint for its existence. However, if the exists as a tetraquark [672, 673, 674], its isospin-1 partner would appear in this channel with . However, no clear signal for such an isovector hidden-charm tetraquark state was found in Ref. [697].
Besides the studies of and in low-energy scattering of the systems, the CLQCD Collaboration also investigated the low-energy near threshold scattering of the system in both the -wave ( with ) and -wave ( with ) channels [700, 701]. In Ref. [701], the authors performed their calculations by using lattice QCD with twisted mass fermion configurations of size with lattice spacing fm at three pion masses MeV, 423.6 MeV, 488.4 MeV. For the scattering in the -wave channel, they obtained the values of the lowest to be in the range , which increased by an order of magnitude compared with their previous quenched lattice QCD result in Ref. [700]. This resulted in stronger attractive interaction between the two charmed mesons. The phase shift was also checked to satisfy , which confirmed the existence of a shallow bound state in this channel. In the -wave channel, similar conclusions were reached by inspecting the lowest values of and the quantity . Based on these results, the author concluded that the interactions between a system are attractive and a possible bound state below the threshold may exist for both pseudoscalar and axial-vector channels. Especially, the resonance candidate in the channel may provide some hints on the nature of the state [701].
9.4 The doubly-charmed/bottom tetraquark states
The idea of searching for stable tetraquark states from the first principle lattice QCD has a long history and is still attractive now. In 1990’s, there were some efforts to calculate the interactions and potential between two heavy-light mesons in lattice QCD [702, 703, 155, 156, 704, 705], inspired by many phenomenological studies of the stability for these systems [173, 142, 141, 144, 706, 145, 147, 707, 150, 151, 152].
In Ref. [703], the authors extracted an effective potential between the heavy-light mesons (HLHL systems) from the quark correlation functions in the framework of quenched lattice QCD with Kogut-Susskind fermions. The correlation functions were calculated by including the two kinds of diagrams corresponding to the pure gluon-exchange part and the quark-exchange part of the mesonic interactions. The gauge field configurations were generated on a periodic lattice with inverse gauge coupling corresponding to a lattice spacing fm. Finally, they were able to extract the potential from the Euclidean time behavior of these correlators. The result showed that the resulting potential was attractive at short heavy-quark distances and the interaction was stronger for smaller light quark masses. They also found that the quark-exchange diagram played a significant role only for distance much less than . This result was supported by the investigation in Ref. [155], where an adiabatic approximation was used to derive the binding energy potential between two heavy-light mesons in quenched SU(2)-color lattice QCD. The derived binding potential was attractive at short and medium range. The UKQCD Collaboration studied the system at fixed heavy quark separation in quenched [156] and unquenched [704] SU(3) lattice QCD. They found evidence for deep binding at small in the and (where and are total isospin and spin for the two light quarks, respectively) cases, for both quenched and unquenched results. The binding energy was about MeV at and very short-ranged. It was essentially contributed by the gluon-exchange diagram and insensitive to the light quark mass. The channel at was attractive for unquenched and repulsive for quenched, which was the only difference between the quenched and unquenched results. The quark-exchange diagram contributed at larger separations fm, where the evidence for weak binding was found. The authors concluded that the exotic -mesons exist as states stable under strong interactions.
However, the extracted potentials remain largely unexplored in the above lattice calculations due to the large statistical uncertainties. In Ref. [708], the NPLQCD Collaboration studied the potentials between two mesons in the heavy-quark limit by choosing a relatively small lattice volume in order to explore the intermediate and short-distance component of the potential. Their calculations were performed on quenched lattices with a spatial length of fm at a pion mass MeV. They found nonzero central potentials in all four spin-isospin channels, where is the total spin of the light quarks. At short distance, they found clear evidence of repulsion between the mesons in the channels and attraction in the channels. Later, the possible tetraquark bound states in HLHL system were suggested in Refs. [709, 710], based on their lattice QCD calculations of the interaction potentials for various channels in the heavy-quark limit.
Besides, the potential was also computed by using the Wilson twisted mass lattice QCD with two flavors of degenerate dynamical quarks ( MeV) [711, 712]. The authors calculated 36 independent trial state potentials using gauge field configurations. They established a simple rule to judge whether a potential is attractive or repulsive: a potential is attractive if the trial state is symmetric under the meson exchange, while the potential is repulsive if the trial state is antisymmetric under the meson exchange, where the meson exchange means the combined exchange of flavor, spin and parity. Using the dynamical results obtained in Refs. [711, 712], Bicudo and Wagner found strong indication for the existence of a tetraquark bound state [355]. Assuming that the pair of antibottom quarks was immersed in a cloud of two light quarks whose screening has little effect on the interaction, the authors suggested the following ansatz to model the heavy antiquark-antiquark potential:
[TABLE]
where is the separation of two antibottom quarks, characterizes the size of the meson and is an exponent. They focused on a scalar isosinglet and a vector isotriplet tetraquark channels, which were found to have attractive interactions between two static-light mesons in the dynamical lattice simulations [711, 712]. Fitting the ansatz Eq. (114) to the lattice results for the heavy antiquark-antiquark interactions, the authors were able to determine all three parameters and for the scalar isosinglet while only two parameters and for the vector isotriplet. They derived an analytical rule for the existence or non-existence of a bound tetraquark state for the radial equation in three dimensions. Applying this rule to angular momentum , they found the condition for having at least one bound state
[TABLE]
where is the reduced antibottom quark mass. According to this condition, the authors only found a strong indication for the existence of a bound state in the isoscalar channel. To investigate the existence of this bound state rigorously, they numerically solved the Schrödinger equation to calculate the binding for heavy tetraquarks. They confirmed the existence of this tetraquark state with a confidence level of around … and binding energy of approximately 30 MeV…57 MeV. For the vector isotriplet channel, however, both the analytical rule and binding energy obtained from solving Schrödinger equation indicated that the existence of a bound state in this channel is rather questionable. Nevertheless, possible sources of systematic errors should be considered, including the choice of the heavy quark mass and the lattice spacing, the lattice QCD finite volume effects, the unphysically heavy light quark masses.
In Ref. [357], they extended their computations for the light quark combinations [355] to similar systems with heavier strange and charm quarks, i.e. and , by investigating and quantifying systematic uncertainties in detail. They confirmed the existence of the isoscalar tetraquark state with quantum numbers and no bound tetraquark state in the isotriplet channel with . They found no and tetraquarks exist in any channel. Later they considered the light-quark mass dependence and computed the potential as a function of their spatial separation , by using the twisted mass lattice QCD (tmQCD) with a lattice spacing around 0.079 fm and three infinite-volume pion mass MeV, 480 MeV, 650 MeV [356]. This allowed for an extrapolation of the potentials to the physical pion mass, showing the tendency that the binding in the scalar isosinglet channel became stronger towards the physical pion mass, while the pion mass dependence was relatively mild close to the physical point. Accordingly, the authors concluded that the two mesons could form a tetraquark state in the scalar isosinglet channel with the binding energy MeV. For the vector isotriplet channel, there was no binding and the results were essentially independent on the pion mass.
Including the effects of the heavy quark spins, Bicudo et. al. reanalyzed the tetraquark binding energy by solving a coupled-channel Schrödinger equation with the Born-Oppenheimer approximation [713]. They found that the spin of heavy quarks decreased the binding energy. However, the attraction was sufficiently strong that the previously predicted bound tetraquark state persisted in channel. Later in Ref. [358], a new tetraquark resonance with quantum numbers was predicted.
However, the above lattice simulations for the tetraquarks were all performed with MeV in the static heavy quark limit. Recently, Francis et. al. investigated the possibility of the and () tetraquark bound states using lattice QCD ensembles with sufficiently low pion masses and 415 MeV and the near-physical [714]. The authors considered one diquark-antidiquark interpolating operator and a meson-meson operator to calculate the correlation functions. They used the NRQCD lattice action [715, 716, 717] to calculate the bottom quark propagators, which avoided the static approximation for the heavy quarks. Finally, they found unambiguous signals for the tetraquark states in the and channels, which lied 189(10) and 98(7) MeV below the corresponding free two-meson thresholds. These and tetraquarks can only decay via the weak interaction, implying that these states are rather stable. Later in Ref. [718], they extended their investigation to study the tetraquark states with the same method. They found the evidence of a stable tetraquark, which lies about MeV below the corresponding threshold.
There were also some studies on the doubly-charmed tetraquark states in lattice QCD. In Ref. [719], Ikeda et. al. investigated the S-wave meson-meson interactions in the and systems with both and , by using -flavor full QCD gauge configurations generated at MeV. They employed the relativistic heavy-quark action to treat the charm quark dynamics on the lattice. Using the HAL QCD method, the authors extracted the S-wave potentials in lattice QCD simulations and then calculated the meson-meson scattering phase shifts and scattering lengths. Their results showed that the interactions were repulsive and insensitive to the pion mass for the channels while they were attractive and growing with for the channels. However, the S-wave scattering phase shifts in these attractive channels indicated that no bound states or resonances were formed at the pion mass MeV, particularly in the channel corresponding to the doubly-charmed tetraquarks with quantum numbers .
The mass spectra for the doubly-charmed tetraquark states were also studied by the Hadron Spectrum Collaboration, using 48 interpolating currents [697]. Using these diverse bases of the meson-meson operators as well as compact tetraquark operators, they calculated the mass spectra for the flavor content with isospin-0 in the lattice irreps corresponding to respectively as shown in Fig. 26, and mass spectra for the flavor with isospin- in the irreps corresponding to respectively as shown in Fig. 27. Similar to the results for the isospin-1 hidden-charm sector, the number of energy levels in the computed spectra for the doubly-charmed sector was equal to the number of non-interacting two-meson levels and the majority of energies were at most slightly shifted from the non-interacting levels. There were no strong indications that there was any bound state or narrow resonance for the tetraquarks in these channels [697].
A comprehensive lattice QCD calculation of the tetraquark states with quark contents and in both spin zero () and spin one () sectors was presented on three dynamical highly improved staggered quark ensembles at lattice spacings of about , and fm [720]. The authors employed the overlap fermion action for the light, strange, and charm quarks [721, 722] while a NRQCD formulation for the bottom quark [716]. They performed the standard method for GEVP to obtain the effective masses of the ground state energy levels. The results of the spin one tetraquark states showed the presence of the ground state energy levels which were below their respective thresholds for all , and light flavor combinations. The mass spectra for all spin zero tetraquark ground state were found to lie above their respective thresholds, strongly disfavor the existence of any heavy bound tetraquark state with spin zero.
9.5 Other tetraquark configurations
After the was observed by the D0 Collaboration, Lang et. al. investigated the -wave scattering using lattice QCD to search for its existence with the flavor and spin-parity [723]. They employed gauge configurations with dynamical quarks and lattice spacing fm at rather low pion mass MeV. For the strange quark, they used a partially quenched setup with the valence mass close to the physical point leading to MeV [724]. The bottom quark was treated as a valence quark using the Fermilab method [725, 726]. They took into account the channel for completeness. They didn’t find a candidate for with .
In Ref. [727], the authors studied the low-lying spectrum of the system using the lattice nonrelativistic QCD methodology to search for a stable tetraquark state below the corresponding lowest noninteracting bottomonium-pair threshold. They constructed a full -wave basis for the meson-meson type of interpolating operators and the diquark-antidiquark type of operators with quantum numbers coupling to and , the coupling to the and coupling to . They employed four gluon field ensembles at lattice spacings ranging from fm, and one ensemble which had physical light-quark masses. All ensembles had and quarks in the sea. The results showed no evidence of a stable tetraquark candidate below the noninteracting bottomonium-pair thresholds in any channel by studying a full -wave color-spin basis of QCD operators. To ensure the robustness of this conclusion, they added an auxiliary scalar potential into the QCD interactions with the objective of pushing a near threshold tetraquark increasingly lower than the threshold. As a result, they found no indication of any state below the noninteracting threshold in the or channel.
9.6 Pentaquarks
The charmonium-nucleon scattering at low energy has been studied in lattice QCD before the LHCb’s observation of the hidden-charm pentaquark states. In Ref. [728], Kawanai and Sasaki studied the charmonium-nucleon potentials for both the and systems in quenched lattice QCD. The central and spin-independent potentials were calculated from the equal-time BS amplitude through the effective Schrödinger equation. They performed the quenched lattice QCD simulations on two different lattice sizes with a lattice cutoff of GeV. They used nonperturbatively improved Wilson fermions for light quarks [729] and a relativistic heavy quark (RHQ) action for the charm quark [730]. As a result, the authors found that the all charmonium-nucleon potentials were weakly attractive at short distances at MeV, while the attractive interaction in the systems were rather stronger than that of the system. However, these potentials were still not strong enough to form a bound state in the system. Moreover, the attractive interaction in the system tends to get slightly weaker as the light quark mass decreases, which may imply that the strength of the potential at the physical point becomes much weaker than that measured at MeV. Similar conclusions were also obtained previously in the preliminary quenched lattice QCD study of the and scattering for MeV [731], and in the computation of the scattering lengths of the charmonium scattering with light hadrons in full QCD at MeV [732]. However, the NPL Collaboration reported a contradictory result in Ref. [733] that the system with has a deeply bound state with binding energy of 19.8(2.6) MeV, by using lattice QCD calculations at very heavy pion mass MeV and a single lattice spacing fm.
To further discriminate the conflicting results between Refs. [728] and [733], Sugiura et. al. performed the time-dependent HAL QCD method to study the charmonium-nucleon effective central interactions in Ref. [734]. The authors focused on three -wave channels , , and , which were a part of the coupled channels that can couple to the hidden-charm pentaquark states. The calculations were employed with flavor full QCD gauge configurations at the lattice spacing fm. The clover action was used for all quarks including charm quarks. They calculated the effective central potentials in the two and one channels for different time slices. The scattering phase shift for was also obtained by solving the -wave radial Schrödinger equation. They found the attractive interactions for all channels, however, not strong enough to form bound states. This result confirmed the time-independent method observation in Ref. [728].
The first lattice simulation reaching the energy region of the hidden-charm pentaquark states and was performed in Ref. [735], in which both the -wave and -wave, both and systems were studied including all possible assignments. The authors considered the nucleon-charmonium systems with total momentum zero of the form , where is nucleon and denotes or . They performed the simulations on the ensemble with a lattice spacing fm at the pion mass MeV. The Wilson Clover action was adopted for all light quarks and charm quarks. The energy spectra were calculated in the one-channel approximation for all six irreps of the discrete lattice group : irrep with , irrep with , irrep with . The final eigenenergies were presented in Fig. 28 for the system and in Fig. 29 for the system. As seen in these plots, the extracted lattice energy spectra were exactly consistent with the prediction of the non-interacting nucleon-charmonium systems. No additional eigenstates were found, which implied that there was no strong indication for a hidden-charm pentaquark resonance in the one-channel approximation for and scatterings.
The hadroquarkonium picture [736] was proposed for the modification of the potential between a static pair (in the heavy quark limit) induced by the presence of octet and decuplet light baryons. In Ref. [737], the authors performed lattice QCD simulations on a CLS ensemble with flavors of nonperturbatively improved Wilson quarks at a pion mass MeV at a lattice spacing fm. The shift of the binding energies due to the presence of the nucleon was extracted to be down only by a few MeV, similar in strength to the deuterium binding. It was dubious whether such a small attraction survives in the infinite volume limit and supports the existences of bound states or resonances.
9.7 A short summary for lattice QCD
We summarize some important conclusions from the present lattice QCD investigations as follows:
- •
Lattice QCD simulations found an isosinglet candidate for the state close to the experimental data, only if both the conventional operators and the two-meson operators were used. However, they found no candidate if the diquark-antidiquark and interpolators were used in the absence of , implying the Fock component is crucial for the . The inclusion or exclusion of the diquark-antidiquark operators does not affect the lattice energy levels significantly.
- •
The lattice QCD simulations for the mass spectra of the hidden-charm four-quark systems didn’t find any additional energy level in the channel associated with the and states. The lattice investigations of the multiple channel scattering indicated that the is not a conventional resonance but a threshold cusp.
- •
A quenched lattice QCD calculation showed that the interactions between the system are attractive for both the pseudoscalar (-wave) and axial-vector (-wave) channels. The might be a candidate of the -wave molecular state.
- •
The hybrid charmonium mass from Lattice calculations is consistent with the experimental mass of the , which supported the hybrid interpretation for this state.
- •
Lots of lattice simulations showed that the tetraquark state with quantum numbers lies below the thresholds, implying that this state is stable against the strong and electromagnetic interaction. This observation has been confirmed at the sufficiently low pion masses MeV.
- •
There is no evidence of a tetraquark state with a mass below the lowest noninteracting threshold in the or channel in the lattice nonrelativistic QCD studies.
- •
To date, there are few lattice studies for many other structures reported in the last 16 years. Future investigations are badly needed to study these claimed signals in various channels.
- •
Most of the present lattice simulations were performed with a rather large pion mass. The pion exchange force plays a pivotal role in the formation of the loosely bound molecular states, which decays exponentially with and is very sensitive to the pion mass. Full QCD simulations near the physical point are crucial to explore the underlying structures of the near-threshold states.
- •
The lattice studies of the and scatterings in the one-channel approximation found no strong indication for the existence of and pentaquark states. The strong coupling between the with other two-hadron channels might be important for the formation of the states. The multiple channel scattering need to be investigated in the future lattice simulations to confirm or refute the existence of the states.
10 Production and decay properties
In this subsection, we shall discuss recent theoretical progress on the productions and decay properties of the multiquark states, which are crucial to understand the nature of exotic states. In particular, the non-resonant explanations for some states depend on their production processes.
The production of the multiquark systems at the quark level remains very challenging since nonperturbative effects are usually involved. One has to turn to effective formalisms such as the nonrelativistic QCD (NRQCD) framework [738, 162, 161, 739, 740] and quark combination models [741]. At the hadron level, one may check whether it is necessary to introduce new resonances through refitting the production cross sections and line shapes. There are investigations of the interference effects on line shapes, values, and resonance properties due to the nearby hadrons. There are also explorations of amplitude pole structures in understanding the nature of resonances by considering the rescattering or meson loop mechanism. An interesting observation is that the anomalous triangle singularity (ATS) in final state interactions through the triangle diagram may mimic resonance-like structures in the mass invariant distributions of measured channels, when all the three particles in the loop approach their on-shell conditions simultaneously. If this mechanism were confirmed, the number of genuine resonances might be reduced. In the literature, there are also discussions about productions of exotic hadrons in heavy ion collisions [742, 743, 744].
One may adopt the widely used Cornell model or pair creation model to study the decays of higher charmonium(like) states [745, 746]. The quark-interchange model was proposed for the rearrangement decays of the multiquark states. [747, 748, 749, 750]. The rescattering mechanism may be used to study the decays at the hadron level [751].
10.1
10.1.1 The states
In Ref. [752], Zhou, Xiao, and Zhou reanalyzed the and processes where the [now called ] and the were observed, respectively. Their results indicate that these two mesons are the same state. In Ref. [753], Baru, Hanhart, and Nefediev investigated whether the is indeed the tensor partner of the in the molecule picture [754] by analysing the production amplitudes. They found that current data favor a scalar assignment, if the is a molecule. Its nature would require some exotic interpretation (i.e., neither a regular quarkonium state nor the spin partner of the ), if it is indeed dominated by the helicity-0 contribution of the nearby tensor state.
The hadronic effects on the meson abundance in heavy ion collisions were studied in Ref. [755], where the authors found that the absorption cross sections and the time evolution of the meson abundance were strongly dependent on its structure and quantum numbers. Combining the double parton scattering and basic ideas of the color evaporation model, the authors of Ref. [756] developed a model for the tetraquark production. They predicted the production cross section at LHC TeV to be in the order of tens of nb. In Ref. [757], Wang and Zhao discussed dependence of the ratio between the semileptonic and nonleptonic decays into the . They demonstrated that the ratios would be universal and predictable if the production mechanism is through the component, while significant deviations from the predicted results would be a clear signal for the non-standard charmonium structure at the short distance.
In Ref. [758], the authors emphasized the connection between the decays and , with () since these processes together with the case arise from the transition at the quark level. Assuming the and to be tetraquark states, Hsiao and Geng predicted the branching ratios of and are and , respectively, which were expected to be accessible at LHCb. In Ref. [759], they predicted the branching fractions for , , , , , , and .
In Ref. [760], productions of the and in exclusive weak decays of were studied with the improved Bethe-Salpeter method by assigning these two mesons as and , respectively. The predicted branching ratios are Br and Br. Nonleptonic decays to these states were also investigated. In Ref. [761], the same group investigated the two-body open-charm OZI-allowed strong decays of these two mesons with the model. The decay width of is around 34 MeV and close to the width of . Therefore, the is a good candidate of . Although the width of is around 70 MeV and close to the lower limit of the width of , the ratio of is larger than the experimental data of . The authors concluded that is not the candidate of .
In Ref. [270], Liu tried to understand the nature of , , , and in the process with the rescattering mechanism. The and may be simulated by the rescattering effects, but the and could not if the quantum numbers of and are and , respectively, which indicates that the and could be genuine resonances. The arguments from the mass, decay, and production [33, 269, 762, 60, 59] lead to the suggestion that the may be the conventional . In Ref. [763], Liu and Oka investigated the radiative transition processes , , and to search for the states, such as the , , , and , by considering the rescattering processes via the charmed-strange meson loops. They found that the anomalous triangle singularity (ATS) may appear as narrow peaks, when some special conditions are satisfied, and thus nonresonance explanation is possible for their structures. Contrary to genuine resonances, the location of the ATS peaks depends on kinematic configurations, and their movements may be used to distinguish the ATS peaks from genuine resonances, once high resolution experiments are available.
In Ref. [764], the production of as a molecule was considered in the process by using the heavy meson effective theory. In Ref. [765], the production and absorption in the hot hadron gas produced in heavy ion collisions were studied. The authors found that the average number as a molecule, , is about 80 times larger than that as a tetraquark state, since the production mechanisms are quite different. Similar investigations about the and cases can be found in Refs. [766, 767, 768].
The authors of Ref. [769] studied the production and decay of the exotic hadrons in multiproduction processes at high energy hadron colliders in the non-relativistic wave function framework by treating the exotic states as hadronic molecules. The rapidity and transverse momentum distributions of the , , and at TeV in collisions were investigated.
The authors of Ref. [770] discussed the possibility of the being the state through the invariant mass spectrum in the process. The analysis for the invariant mass in the decay results in a candidate with a mass around 4080 MeV. From the decay behaviors of the , the assignments of the and as the and are reinforced respectively.
In Ref. [771], Kang and Oller introduced a near-threshold parameterization to analyze the production processes , , and . The data can be reproduced with similar quality for being a bound and/or a virtual state. The authors of Ref. [772] studied the mass distribution of the reaction, and found that the contributions of the narrow and the are needed, which provided an explanation for the large width of observed by LHCb [59, 60]. The state was found to be strongly tied to the channel.
In Ref. [773], Dai, Dias, and Oset studied the and reactions. They found that the molecular states and couple mostly to in the and channels respectively, and have big influence on the mass distribution. Later in Ref. [774], the authors studied the semileptonic decay of the meson into the , , and resonances by treating them as dynamically generated states.
In Ref. [775], the productions of the and in interactions in collisions at the LHC were investigated. The results with both and indicate that the experimental study is feasible and can be used to check the existence of these two states. In Refs. [776, 777], Braaten, He, and Ingles pointed out that the production rate of accompanied by a pion at a high energy hadron collider should be larger than that of without a pion, if the is a weakly bound charm-meson molecule. This type of production may happen through the creation of the pair at short distances followed by a rescattering process. In Ref. [778], the production of accompanied by a pion in meson decays was discussed. Such a production should be observable, which can be used to support for the interpretation of as a molecule.
In Ref. [779], the transverse momentum spectra and yields of the , , and for heavy-ion collisions at the LHC energies are predicted within the framework of the statistical hadronization model. The production yield of was found to be about 1% relative to that for . In Ref. [780], Guo proposed a method to study whether the is above or below the threshold. He found that the line shape of is very sensitive to the deviation of the mass from the threshold due to a triangle singularity. This indirect method can be applied to experiments producing copious pairs.
In Ref. [781], the production of molecule in the process was investigated by using the rescattering mechanism. The production ratios are orders of . In Ref. [782], a comparative study for the strong decays of as a charmonium state into and channels were performed with the Cornell potential model, and a generalized screen potential model. The suppression of the mode can be understood with both models, while the significant mode favors the latter model.
In Ref. [334], the decay was investigated in the Bethe-Salpeter equation approach if it is a molecule. The decay width was found to be in the range 8.749.5 keV. The authors of Ref. [783] studied the two-photon decay of the state if it as a molecule. The obtained decay width is 36 keV. The decay width into is 66 MeV in the same molecular picture [784].
10.1.2 The states
In Ref. [785], in the Fano-like interference scenario, Chen, Liu, Li, and Ke studied the asymmetric line shapes of the and structures in their discovery processes and by considering the interference between continuum and resonance contributions. They found that the broad structure can be induced by the Fano-like interference. Their numerical results indicate that the , , and may be not genuine resonances, which may explain the absence of these states in the value scan and the non-observation of their open-charm decay modes. The Fano-like phenomena were also expected in processes such as and . In Ref. [786], Chen, Liu, and Matsuki studied the interference effects in understanding the and . They found that the signals of these two structures can be reproduced by introducing the interference effects between , , and background contributions. Therefore, the and may be not genuine resonances. A better description of experimental results needs the introduction of , which was proposed to be the charmonium . In Ref. [787], the updated data of were used to categorize this state into the family in a - mixing scheme. The present experimental data seem to support the charmonium assignment. The assignments of higher charmonia states and properties of newly predicted states were also discussed.
In Ref. [745], the open-charm decay mode of was explored in the pair creation model by assuming this exotic meson to be a -wave diquark-antidiquark state [97]. The decay mechanisms of in both tetraquark and molecule pictures were studied in Ref. [788]. The , , and decay modes of were discussed in the rescattering mechanism in Ref. [751].
In Ref. [789], the role of the s-channel in the reaction near threshold was investigated. The authors found that this state gives a clear bump structure with the magnitude of 10 b, which can be tested by the PANDA experiment. In Ref. [790], productions of in the processes and were studied simultaneously, where the nucleon and its excitations play an important role in the reactions. Therefore, the PANDA experiment can also be used to study the properties of .
In Ref. [260], Qin, Xue, and Zhao investigated the production mechanism of in annihilation by assuming it to be a molecule admixed with a compact component. They found that the heavy quark spin symmetry breaking plays a crucial role for the production, and the molecule picture can describe the available observables simultaneously. They also found that is one of the important decay channels of . In Ref. [791], Cleven and Zhao demonstrated that the molecule picture of can explain the cross section line shape of around the mass region. Further in Ref. [792], Xue, Jing, Guo, and Zhao showed that the partial widths of the molecular to and are much smaller than the width to , and the contributions to the cross section line shapes of the and are rather small in most of the energy region from thresholds to about 4.6 GeV. They also found that the interferences of the with nearby charmonium states , , and produce a dip around 4.22 GeV in the cross section line shape.
In Ref. [793], Dai, Haidenbauer, and Meissner investigated the reaction at energies close to the threshold after considering the effects from the resonance and the final state interactions. The mass and width of this state were found to be around 4652 MeV and 63 MeV, respectively, and thus the and could be the same state. In Ref. [794], a combined fit was performed to the cross sections of , , , , and with three resonances , , and . The authors found that the first two resonances are sufficient to explain these cross sections below 4.6 GeV, so the , , and should be one state.
In Ref. [795], the loop contributions to the reaction around 4.0 GeV were considered with effective Lagrangians. The authors noticed a second pole dynamically generated by the meson-meson quantum fluctuations in addition to the charmonium . This pole leads to the controversial state observed in . The authors also expect that the pole could be visible in the process . The same mechanism may be applied to the and .
In Ref. [796], the relation between the and within a unitarized effective Lagrangian approach was studied. The authors found that the is not an independent resonance but a manifestation of the due to the final state interactions between the and the . In Ref. [797], Chen, Dai, Guo, and Kubis studied the processes in the dispersion theory, and analyzed the roles of the light-quark singlet and octet states in the transitions. According to their investigations, the contains a large light-quark component, and it is neither a hybrid nor a conventional charmonium state. The analysis of the ratio of the octet and singlet components indicates that the has a sizeable component, although it is not a pure molecule.
Very recently, BESIII studied the reaction at center-of-mass energies above 4.08 GeV, and observed and in this process. Neither fast rise of the cross section above the threshold nor obvious structure is visible from the measurements, which is different from the expectation of Ref. [798].
In Ref. [799], Chen et al studied the decay of into by assuming that they are and molecules, respectively. The partial width was found to be around 1.23.0 MeV. In Ref. [800], the open-charm decays of the charmonium were studied in the pair creation model. The total width is consistent with those of both and . One of the dominant decays of is .
In Ref. [746], Xiao et al. studied the decay mode of the higher vector charmonium states around 4.6 GeV in the pair creation model. They found that the partial decay widths for -wave charmonia are about a few MeV, and those for -wave charmonia are less than one MeV. The is very likely to be an -wave charmonium, if the and are the same state.
10.1.3 The states
In Ref. [801], the final state interaction effects on the was studied. In the analysis, the authors calculated the amplitudes of the , , and processes, and searched poles in the complex energy plane. They found only one pole in the nearby energy region in different Riemann sheets and concluded that the is a molecular state according to the pole counting rule [802, 803]. In Ref. [804], the same group noted that the triangle singularity mechanism did not produce the observed peak. Their study again favors the molecular assignment for the .
In Ref. [805], to understand the nature of , the authors considered the and coupled-channel T-matrix analysis in a finite volume. Comparing their calculations with the lattice results [686], they concluded that it is difficult to distinguish whether the is a resonance originating from a pole above the threshold or a virtual state below the threshold. In Ref. [806], the authors carried out a thorough amplitude analysis of the data and considered four different scenarios corresponding to the pure QCD states, virtual states, or purely kinematical enhancements. They concluded that it is not possible to distinguish between these scenarios at that time, which contradicted the conclusion of Ref. [804].
Both the and were observed in the annihilations (into and , respectively), but only the was observed in decay. In order to understand this puzzle, Yang, Wang, and Meissner studied the isospin amplitudes in the exclusive decay into in Ref. [807]. Their analysis indicates that the production of the isovector state is highly suppressed compared to the isospin singlet one, if both states are molecules. Their findings seem to support the molecular nature of .
In Ref. [808], He and Chen investigated the - interaction in the decay in a coupled-channel quasipotential Bethe-Salpeter equation approach in order to understand the origin of the and . They found a virtual state around 3870 MeV, and concluded that the two states are from the same virtual state. In Ref. [809], Zhao analyzed possible quantum numbers and various explanations of the state via its production in the decay process. He found that the state can be either induced by final state interaction effects arising from an -wave rescattering with , or it is a -wave resonance with arising from the interaction. The latter structure can be regarded as an excited state of .
In Ref. [810], Cao and Dai discussed possible assignments for the , [called in PDG [35]], and and proposed that the observed in by LHCb [67] and the observed in by Belle [811] correspond to the same state. They suggested that the should be a or state, while the could be a or state.
In Ref. [812], the authors demonstrated that the and can be consistently interpreted as kinematical singularities from the triangle diagrams. They explained why the contribution was observed in but was not. In Ref. [813], the and structures observed in were also explained with the triangle singularities, and the was predicted to be ( or ) for ().
In Ref. [814], the production of and in decays was studied with an effective Lagrangian approach by considering meson loop contributions. The branching ratios of decays to and are around and , respectively. In a light-front model, Ke and Li studied the decays of and [called in [35]] into by assuming them as and molecules, respectively [815]. They found that the seems to be a molecular state. If the is also a molecule, it would be observable in , since the partial width of is only three times smaller than that of . The nonobservation of in this channel with precise measurements would be helpful to rule out its molecule interpretation.
The role of the charged exotic states in the reaction was considered with a dispersive approach in Ref. [816]. The authors found that the state plays an important role in explaining the BESIII’s invariant mass distribution at both and 4.258 GeV [66]. The sharp narrow structure at GeV can be explained with a heavier charged state whose mass and width are about GeV and 52 MeV, respectively. At GeV, no intermediate state is necessary. The authors also concluded that the final state interaction leads to the invariant mass distribution for these four energies.
In Ref. [817], the strong decay widths of into and were calculated with the coupling constants in the QCD sum rule method with the assumption that the is a diquark-antidiquark state with . The resulting decay width of these two modes is around 66 MeV, larger than the PDG value MeV [35]. In Ref. [818], the charmless decays of the and [called in [35]] into vector-pseudoscalar (VP) type light mesons were investigated with the rescattering mechanism by assuming these mesons to be and molecules, respectively. The branching ratios of () are typically of the order ().
In Ref. [819], the strong decay widths (, , ) and () were calculated within a covariant quark model. The authors found that the tetraquark-type (molecule-type) current is in disaccord (accordance) with the experimental observation that has a much stronger coupling to than to . They found that the state is a good candidate for the compact tetraquark state, and predicted the partial decay width to be around 24 MeV.
In Ref. [750], Wang et al investigated the strong decays of , , or , , and into a heavy quarkonium plus a pion in a relativized quark model in the molecule picture. The decay proceeds through the interchange of a heavy quark and a light quark as well as the exchange of a gluon. Therefore, the exchanged or is a color-octet state. According to their results, the and have a larger coupling with than , but the partial width is much larger than . This feature is consistent with the experimental observations. The obtained decay ratios also favor the molecule assignments for the states. However, the calculation does not favor the pure molecule assignment for the , which is different from the study with the naive nonrelativistic quark model [820].
In Ref. [821], Voloshin considered the decays and . He assumed that is dominantly an -wave molecule and that the is an -wave resonance, where the transitions occur through . The rates were found to be at most several tenths percent in terms of the branching fraction for the .
10.1.4 The states
In Ref. [822], the role of the states in the () transitions was analyzed by studying the final rescattering effects with the dispersion theory. The effects in and are very small, but those in are significant. The inclusion of the states and interactions can explain the observed anomaly of . Further in Ref. [823], the effects of the two intermediate states and bottom meson loops on the dipion transition processes were studied. The authors found that the contribution from meson loops is comparable to those from the chiral contact terms and the -exchange terms. With the inclusion of ’s and meson loops, the authors can reproduce the experimental two-hump behavior of the spectra in . In the process , they expect a narrow dip around 1 GeV in the invariant mass distribution.
In Ref. [824], Guo et al studied the line shapes of the near-threshold states and with using their parametrization method. They found that the is a virtual state on the second Riemann sheet near the threshold, while the is a resonance on the third or fourth Riemann sheet near the threshold.
In Ref. [825], Bondar and Voloshin discussed the possibility that the production channels in the mass region are dominated by the (not ) due to the triangle singularity in . If this mechanism is dominant, any nonresonant background not associated with the should be suppressed. A similar structure near 11.06 GeV in the annihilations into and was also expected.
In Ref. [826], the hidden-bottom decays of the and into (), (), and were analyzed with the effective Lagrangian approach, via final state interactions or meson loop contributions. With this mechanism, the final state interactions were found to be important, and the branching ratios into and were predicted. In the same framework, the hidden-charm decays of the and states in Ref. [827] were consistent with the recent BESIII measurement [828].
In Ref. [829], the bottomoniumlike and resonances were proposed to be a mixture of the molecular and states with one mixing angle . With , the ratio and a definite interference pattern in the process were predicted.
In Ref. [830], the experimental data for the two states were analyzed simultaneously in the molecule picture using the Lippmann-Schwinger equations. The long-range pion interaction does not affect the line shapes if only -waves are considered, but the situation is slightly different once -waves are included. The study also indicates that the two states can be described by poles on the unphysical Riemann sheets, where the is associated with a virtual state just below the threshold, while the is likely an above-threshold shallow state. In Ref. [831], the same method was applied to the spin partner states ’s with . The pionful (pionless) approach leads to the threshold cusp (above-threshold hump) in the line shapes.
In Ref. [832], the production of and from the decays through bottom-meson loops was investigated. The contributions from the loops dominate those from the loops due to the larger pionic couplings in the triangle diagrams. The branching ratios of and were predicted to be around 4%.
In Ref. [833], the strong decays of and were studied with molecular type currents using the same method in studying the decay widths of and [819]. The and decay modes are suppressed but not as much as the Belle data.
Inspired by BESIII’s new measurement [828], Voloshin estimated the cross section of the processes at the maximum of the resonance to be about 0.1 pb [834], where are the predicted isovector meson-antimeson states from the considerations of the heavy quark spin symmetry and light quark spin symmetry [406].
10.2 and
In Ref. [835], the feasibility of the production and detection of at various facilities was discussed. Its production is comparable to the doubly charmed baryon production and prompt production. Thus, its observations at SELEX, LHC, and Belle should be possible.
The production of the states in the electron-positron process has been discussed in Refs. [161, 162, 739]. To search for these exotic states at factories, the authors of Ref. [739] proposed a more efficient analysis method. If a three-jet event may be identified by their energy and angular distributions, the doubly charmed states can be searched for in the most energetic jet, which may suppress the combination background fluctuations of the reconstructed hadron mass spectra.
In [836], the authors argued that the doubly charmed tetraquarks are a straightforward consequence of the diquark-antidiquark model [95, 97]. Assuming that these tetraquarks are above-threshold particles, the authors explored their decay modes and rates as well as their production mechanisms. They showed that the doubly charmed tetraquarks could be produced at LHC from the decays of and .
In Ref. [837], the authors noted that the hadronic effects on the tetraquark state are insignificant in relativistic heavy ion collisions. The authors also noticed that the and decay modes of should be the most probable channels to reconstruct it from heavy ion collisions. Besides, the order of yields at LHC and RHIC reads, (molecule)(compact).
In Ref. [838], Ali et al explored the production of , , , and (), in -boson decays. The lifetimes of these tetraquarks were estimated to be around 0.5 ps. Besides, their signature decay modes were studied. In Ref. [839], the discovery potential of these tetraquarks at LHC was discussed. According to the estimated cross section , the prospects of discovering these states are excellent.
In Ref. [168], the semileptonic 2-, 3-, and 4-body weak decays of (, ) were systematically investigated. The nonleptonic decays and possible golden decay channels at LHCb and Belle II experiments were also discussed. The weak decays of the tetraquark state to the scalar diquark-antidiquark were investigated in [520] using the QCD sum rule approach, and the lifetime of the was found to be around 9.2 fs, which is shorter than that obtained in Ref. [838].
In Refs. [840, 841], the inclusive decay mode was proposed as a potential discovery channel for the doubly bottom baryon at LHC, which is based on the fact that can only arise from the weak decay of a hadron with two bottom quarks. If the states could be observed, such a decay mode may also be used to search for the bound states.
The possible observation of the and states with at LHC was discussed in Refs. [177, 178]. The color connections of states and the induced hadronization effects on and in annihilation processes were studied in [842]. A study of production in double parton scattering was given in Ref. [756]. Possible channels to search for such a state were discussed in Refs. [843, 844].
In Ref. [845], the authors investigated the lifetime and weak decays of the fully heavy tetraquark state with . They obtained a lifetime value around 0.10.3 ps and discussed the golden channels to search for this state.
10.3
The analysis of the process in Ref. [846] by Roca, Nieves, and Oset supports the as a molecule of mostly and . In Ref. [847], an improved analysis by Roca and Oset indicates that the existence of the state cannot be undoubtedly claimed with only the fit to the and mass distributions. The investigation of and in Ref. [848] indicates that the s-/u-channel contributions are not negligible. The obtained bound states and line shape of the mass distribution favor the assignment of the as a molecule.
In Ref. [849], the process was studied with the contributions of a hidden-charm pentaquark-like molecule state [184, 193]. A clear peak in the mass distribution was found. In Ref. [850], the contributions of the hidden-charm states with and to the reaction were discussed. These states lead to clear peak structures in the total cross sections. The contributions of a hidden-beauty to the scattering were discussed similarly in Ref. [851]. The authors of Ref. [852] suggested that a clean peak should be seen for a hidden-charm resonance (mass4265 MeV), which couples to the channel in the decay .
In Ref. [853], the production of the states and their yields in ultrarelativistic heavy ion collisions were discussed in both the compact pentaquark picture and the molecule picture. In Ref. [854], the production cross section for a hidden-charm pentaquark in proton-nucleus collisions was found to be sizable by assuming if it is a molecule of a charmonium and a light baryon.
In Ref. [855], Liu and Oka discussed the rescattering effects in the reaction due to the intermediate open-charm hadrons. They showed that the triangle singularity peaks can simulate the pentaquark-like resonances in the invariant mass distributions. In Ref. [856], Guo et al proposed that the observed in the process might be due to a triangle singularity around the threshold. A discussion on triangle singularities in Ref. [857] indicates that the narrow would have an origin other than the triangle singularity from the -charmonium-proton intermediate states, if this baryon has quantum numbers .
In studying the nature of with the pole counting rule in Ref. [801], the authors suggested that the pentaquark can be related to by replacing a pair by a . They speculated that the should be a molecule with rather than .
In Ref. [858], Kim et al considered the reaction based on the hybridized Regge model and included the contributions of the LHCb states. The total cross sections were found to be about 1 nb in the vicinity of the heavy pentaquark masses. In Ref. [859], Karliner and Rosner discussed the photoproduction on deuterium, through which the simultaneous investigation of the and isospin partner states can be achieved. The photoproduction of on a deuteron target may be used to check whether it is a genuine isospin-half resonance.
The authors of Ref. [860] studied the pentaquark production in the semileptonic decays where . The branching ratios are about two orders of magnitude smaller than the decay , which is accessible to the ongoing experiments at the LHCb.
In Ref. [740], the hidden-charm pentaquark production at the colliders was discussed in the NRQCD factorization framework, by treating the pentaquark as a state with a color-octet configuration. It is possible to search for the direct pentaquark production signal at colliders. The authors of Ref. [219] discussed the production mechanism of the pentaquarks in decays and assumed that they are states. The decay happens through followed by a gluon emission () from . Then and form the kaon meson while the other quarks form the pentaquarks. The two color octet and components in -wave (-wave) lead to the () state. This mechanism can explain why only a few states with nonminimal constituents were discovered so far.
In Ref. [861], the difference between partial-wave analysis formalisms used in the construction of three-body amplitudes involving fermions was discussed, particularly in the decay . The authors found meaningful effects on the resonance pole position extraction, which is particularly relevant when several resonances overlap and the quantum number assignment is not stable. In Ref. [862], the authors discussed the role of the hidden-charm pentaquark resonance in photoproduction off and target nuclei near threshold. The presence of produces additional enhancements above threshold in the total creation cross section on nuclei, which can be tested in future JLab experiments.
In Ref. [863], Voloshin discussed the hidden-charm pentaquark production in collisions, which is possible due to the coupling of charmonia to the channel. He found that the pentaquark formation may happen with the nucleons moving slowly inside the deuteron due to the masses of the pentaquark, charmonium, and the nucleon being close to a special kinetic relation. The formation cross section of hypothetical pentaquark states decaying to was much larger than that of pentaquarks into , because of a much larger decay width.
The authors of Ref. [864] studied the strong decay mode of the states in the molecule scenario. The partial decay widths are significantly different for various assignments. The -wave pictures for the and and the assignment for the with are all allowed by the present experimental data.
In Ref. [865], the different decay properties of the and molecules were discussed. The was proposed as a molecule and the decay channel was suggested to be used to disentangle its nature. In Ref. [866], the decay properties of the and were studied with the effective Lagrangian framework by treating them as meson-baryon molecules. The and molecules have different decay branching ratios. The () was proposed to be a ( ) molecule. The cross sections for the processes and through the -channel states were also calculated. The same framework was used to discuss the decay behaviors of the strange and beauty partners of the hadronic molecules in Ref. [867].
Besides the spectrum, the authors of Ref. [375] also calculated the partial decay width of the into to be around 11 MeV by treating it as a hadroquarkonium state. Later in Ref. [868], the decays of the pentaquark were further investigated in both hadrocharmonium and molecular pictures. The authors found that the decay patterns are vastly different. The decays of the molecule into are strongly suppressed, while the opposite happens in the hadrochamonium case.
In Ref. [524], the authors studied the strong decay width of into by treating it as a molecule with or . The coupling constants were calculated with QCD sum rules, and the obtained decay width is around 187 (213) MeV for the case ().
After the new pentaquark states , , and were announced [68], Guo, Jing, Meissner and Sakai suggested to search for the in the mode in order to test its nature [382]. If the state is a molecule, the authors pointed out that the isospin breaking decay ratio would be at the level ranging from a few percent to about .
10.4 A short summary
We give a short summary for the production and decay properties of the heavy tetraquark and pentaquark states based on the recent studies in the literature.
- •
The states. In the conventional charmonium picture, the could be the while the or could be the . In the molecule picture, the could be a molecule.
- •
The states. From the theoretical studies of production and decay, it seems that not all states (, , , , , , , and ) shall exist. The interference effects, coupled channel effects, or final state interactions may make several genuine states behave like more resonance signals. Up to now, we have much deeper understandings about the nature of than 15 years ago, although more experimental measurements are still needed. Different from the , the structures in the mass region GeV seem more complicated and challenging. The conventional charmonium states, theoretically expected but elusive hybrid charmonium mesons, and molecular states coincide here. Do they mix? The nature of is still far from being understood, although it was observed almost 15 years ago.
- •
The and states. The studies of productions and decays favor the molecule assignments for the and . Whether the state is a molecule, a compact tetraquark, or a structure induced by triangle singularity is still unclear. One cannot distinguish the nature of with the present experimental decay data, either. Is it a P-wave resonance or a state caused by final state interaction effects? The and should be the and molecular states, respectively.
- •
The states. The and molecule assignments for and , respectively, can explain experimental measurements for their production and decay properties, although their nonresonance interpretations cannot be excluded.
- •
So far, a large number of investigations on exotic productions have been presented in the literature while the investigations on decays are still lacking. Further studies on decay properties of the exotic hadrons, especially at the quark level, are necessary, which play a key role in understanding the inner structures of the exotic states. On the experimental side, searching for expected multiquarks, checking predicted features, and measuring more physical quantities are definitely needed.
- •
At present, all the multiquark candidates are produced at , , , , and colliders. Future , , and colliders are still under construction. Besides these collision processes, the interactions, charmonia decays, bottomonia decays, decays, and decays are also used to study the multiquark states. From the experimental fact that all the confirmed exotic candidates were observed in the (), , or hadron decay processes, it seems that producing multiquark states in multiproduction processes is very difficult and undetectable [869, 741]. If this is true, searching for the conventional baryons should be more feasible and important than searching for compact states. The exotic hadrons we may observe in the near future would be dominantly states in the charm sector. They are indeed the structures being studied at LHC (see Table 26 of Ref. [870]).
11 Summary and perspective
At the end of this review, we would give a brief summary and perspective for the progress on the multiquark physics. It is useful and instructive to recall the establishment of the quark model before we discuss the multiquark systems. The quark model was first proposed in 1964 to deal with the classification of hadrons. In 1974, the discovery of the meson gave the direct evidence for the existence of the charm quark, which was known as “November Revolution”. Actually, the existence of the charm quark was already predicted in the GIM mechanism before its observation [11]. The bottom quark was discovered by E288 experiment in 1977 at Fermilab. After that, plenty of bound states and excitations of the heavy mesons and baryons have been observed in the following decades [35]. We show the roadmap of the conventional hadrons in Fig. 30.
Although we still can’t understand QCD confinement from the first principle, the quark model has achieved indubitable success in systematizing the properties of mesons and baryons. For the exotic hadron states such as the multiquark systems, one may turn to those phenomenological models for help which are closely related to QCD since the rigorous solutions of soft QCD are still unavailable. The interplay between the phenomenological models and experimental observations will also improve their prediction capabilities. The discovery of the in 2017 was a successful example for the collaboration between experimenters and theorists. Inspired by this observation, one will further expect the future discovery of the baryon.
Searching for the multiquark states is another central issue of the hadron physics. In the past several decades, there accumulated abundant investigations on tetraquarks and pentaquarks based on various phenomenological models. However, the progress was still rudimentary and far away from our expectation due to the lack of the sufficiently accurate experimental results. The year of 2003 is very important in the history of the exotic hadrons. Since 2003, lots of quarkonium-like states and pentaquark states have been observed in the modern facilities, such as BaBar, Belle, CLEO-c, D0, CDF, BESIII, and LHCb etc.
The desires for the nature of these new hadron states have promoted the model building. Some model predictions agree with the experimental data. For example, the masses of some states can be reproduced in the tetraquark configurations. However, there are still many observables which can not be computed accurately enough to be compared with the experimental measurements, e. g., the branching fractions and decay widths for some important channels. Moreover, some new resonances can usually be explained in different models, implying the present phenomenological investigations are unable to distinguish their underlying structures. More efforts are needed to bring down the error bars in future.
Nevertheless, we can still find some critical breakthroughs in the phenomenological aspect. The investigations for the doubly heavy tetraquarks with various methods such as the QCD sum rules [515], CMI model [109, 166], LQCD [720], and the heavy-quark symmetry prediction [164] lead to consistent predictions that the doubly bottom tetraquark states would probably be stable against the strong interaction. Such an agreement may imply critical interaction mechanisms in these tetraquark systems, which can be a breaking point for our understanding of the multiquark states. To date, most of the multiquark candidates were observed in the charmonium sector, such as the hidden-charm XYZ states and pentaquark states. From the theoretical point of view, more tetraquark and pentaquark states with different flavor and color structures and different quantum numbers should also be possible. In the past several years, the authors have made some efforts along this direction. We predicted the pentaquark states with hidden-charm and open-strange flavors. Based on the discovery of , we studied the interactions between the baryon and -wave charmed mesons and predicted the existence of the triply charmed pentaquark states.
Last but not the least, we shall emphasize the unique and pivotal role of the lattice QCD in the search of the multiquark states. The lattice QCD simulations start from the first principle and can investigate the few-hadron reaction processes and extract important information such as the resonance properties, the hadronic interactions, scattering amplitudes, phase shifts, binding energies, mass spectra, decay widths and so on. In the past few years, tremendous lattice results from quenched to unquenched calculations become accurate enough to be compared with both the experimental data and phenomenological results. Especially the simulations using a large basis of operators including the multi-hadron operators have become a unique, vital and powerful tool to explore the near-threshold exotic states. Some very successful examples include the lattice QCD investigations of the and states, the doubly bottom tetraquark states, which provide very clear and intuitive picture of the underlying structures of the exotic states. However, only the and scattering processes were studied in lattice QCD for the pentaquark systems, which is far away from pinning down the existence of the hidden-charm pentaquark states. In terms of LHCb’s recent discovery of the three hidden-charm pentaquark states, the lattice QCD simulations of the charmed baryon and anti-charmed meson are badly needed. At the present stage, many lattice simulations were performed at the unphysically large pion mass due to the expensive computational cost. However, the existence of the molecular type of exotic hadrons is very sensitive to the pion mass since the pion exchange force decays exponentially with the pion mass. Further dynamical lattice QCD simulations using ensembles with near-physical or even the physical pion mass will be extremely valuable for the systems. As shown in Fig. 31, the interplay among experiments, phenomenological models and lattice QCD and their joint progress will definitely sharpen our understanding of the new hadron spectroscopy and QCD itself.
When we were preparing this review, the LHCb Collaboration reported a new charmonium state with via collision in the PHI-PSI 2019 workshop. The BESIII Collaboration has observed the processes and for the first time [871] and measured the cross section of the process [872], which provided new useful information to study the states. Very recently, LHCb has observed CP violation in the charmed meson decays [873], which is definitely an important milestone in the history of particle physics. These excellent observations at LHCb demonstrated its distinguished performances and capabilities to study the exotic hadron states. The Belle II Collaboration released the Belle II Physics Book [874] in 2018 and started to take data with a fully instrumented detector in March, 2019. Although the main aim of Belle II experiment is searching for CP-violation and revealing violations of the symmetry between particles and anti-particles, it will also search for the exotic hadrons and make precision measurements of their properties. To some extent, BelleII shall become a factory of the charmonium-like states. The BESIII experiment at BEPCII will continue to contribute to the field of the states in the near future. They will perform unique investigations on the states around GeV, 4.38 GeV and 4.6 GeV. Moreover, the upcoming PANDA experiment will also study the multiquark and other exotic hadrons in future. The coming decade shall witness a new landscape of the hadron spectroscopy.
Acknowledgments
We would like to express our gratitude to all the collaborators and colleagues who contributed to the investigations presented here, in particular to Dian-Yong Chen, Rui Chen, Xiao-Lin Chen, Er-Liang Cui, Wei-Zhen Deng, Jun He, Ning Li, Shi-Yuan Li, Xiao-Hai Liu, Yu-Nan Liu, Zhi-Gang Luo, Li Ma, Takayuki Matsuki, Zong-Guo Si, T. G. Steele, Zhi-Feng Sun, Guan-Juan Wang, Jing Wu, Lu Zhao, Zhi-Yong Zhou. We thank Er-Liang Cui for helping prepare some relevant documents. This project is supported by the National Natural Science Foundation of China under Grants No. 11722540, No. 11775132, No. 11575008, No. 11621131001 and National Key Basic Research Program of China (2015CB856700), the China National Funds for Distinguished Young Scientists under Grant No. 11825503, the Fundamental Research Funds for the Central Universities, and the National Program for Support of Top-notch Young Professionals.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] H.-X. Chen, W. Chen, X. Liu, S.-L. Zhu, The hidden-charm pentaquark and tetraquark states, Phys. Rept. 639 (2016) 1–121. ar Xiv:1601.02092 , doi:10.1016/j.physrep.2016.05.004 . · doi ↗
- 2[2] L. B. Okun, The theory of weak interaction, in: High-energy physics. Proceedings, 11th International Conference, ICHEP’62, Geneva, Switzerland, Jul 4-11, 1962, 1962, pp. 845–866.
- 3[3] M. Gell-Mann, A Schematic Model of Baryons and Mesons, Phys. Lett. 8 (1964) 214–215. doi:10.1016/S 0031-9163(64)92001-3 . · doi ↗
- 4[4] G. Zweig, An SU(3) model for strong interaction symmetry and its breaking. Version 2, in: D. Lichtenberg, S. P. Rosen (Eds.), DEVELOPMENTS IN THE QUARK THEORY OF HADRONS. VOL. 1. 1964 - 1978, 1964, pp. 22–101.
- 5[5] G. Zweig, An S U ( 3 ) 𝑆 𝑈 3 SU(3) model for strong interaction symmetry and its breaking. Version 1, 1964.
- 6[6] V. E. Barnes, et al., Observation of a Hyperon with Strangeness -3, Phys. Rev. Lett. 12 (1964) 204–206. doi:10.1103/Phys Rev Lett.12.204 . · doi ↗
- 7[7] M. Gell-Mann, The Eightfold Way: A Theory of strong interaction symmetry, 1961.
- 8[8] Y. Ne’eman, Derivation of strong interactions from a gauge invariance, Nucl. Phys. 26 (1961) 222–229. doi:10.1016/0029-5582(61)90134-1 . · doi ↗
