Possible pentaquarks with heavy quarks

Inspired by the discovery of two pentaquarks Pc(4380)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{c}(4380)$$\end{document} and Pc(4450)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{c}(4450)$$\end{document} at the LHCb detector, we study possible hidden-charm molecular pentaquarks in the framework of quark delocalization color screening model. Our results suggest that both Nηc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\eta _{c}$$\end{document} with IJP=1212-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{IJ}^{P}=\frac{1}{2}\frac{1}{2}^{-}$$\end{document} and NJ/ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{NJ}/\psi $$\end{document} with IJP=1232-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{IJ}^{P}=\frac{1}{2}\frac{3}{2}^{-}$$\end{document} are bounded by channels coupling. However, NJ/ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{NJ}/\psi $$\end{document} with IJP=1232-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{IJ}^{P}=\frac{1}{2}\frac{3}{2}^{-}$$\end{document} may be a resonance state in the D-wave Nηc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\eta _{c}$$\end{document} scattering process. Moreover, Pc(4380)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{c}(4380)$$\end{document} can be explained as the molecular pentaquark of Σc∗D\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma ^{*}_{c}D$$\end{document} with quantum numbers IJP=1232-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{IJ}^{P}=\frac{1}{2}\frac{3}{2}^{-}$$\end{document}. The state Σc∗D∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma ^{*}_{c}D^{*}$$\end{document} with IJP=1252-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{IJ}^{P}=\frac{1}{2}\frac{5}{2}^{-}$$\end{document} is a resonance, it may not be a good candidate of the observed Pc(4450)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{c}(4450)$$\end{document} because of the opposite parity of the state to Pc(4380)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{c}(4380)$$\end{document}, although the mass of the state is not far from the experimental value. In addition, the calculation is extended to the hidden-bottom pentaquarks, and similar properties to that of hidden-charm pentaquarks system are obtained.


Introduction
Multi-quark states were studied even before the advent of quantum chromodynamics (QCD). The development of QCD accelerated multi-quark study because it is natural in QCD that there should be multi-quark states, glueballs, and quarkgluon hybrids. After more than 40 years of quark-model study, the ideas of a baryon and a meson is about to go beyond the naive picture: q 3 baryon and qq meson. The proton spin puzzle could be explained by introducing the q 3 qq component in the quark model [1,2]. In order to understand the baryon spectroscopy better, the five-quark component of the proton was proposed [3]. The baryon resonance is certainly coupled to the meson-baryon scattering state and should be studied by coupling q 3 with q 3 -qq scattering channel in a e-mail: jlping@njnu.edu.cn a quark-model approach. Although the strange pentaquark state + claimed by experimental groups 13 years ago might be questionable (the LEPS Collaboration insists on the existence of the pentaquark + [4]) and the multi-quark states might be hard to identify, the multi-quark study is indispensable for understanding the low energy QCD, because the multi-quark states can provide information unavailable in the qq meson and the q 3 baryon, especially the property of a hidden color structure.
In the past decade, many near-threshold charmoniumlike states have been observed at Belle, BaBar, BESIII, and LHCb, triggering lots of studies on the molecule-like hadrons containing heavy quarks. In the heavy quark sector, the large masses of the heavy quarks reduce the kinetic of the system, which makes it easier to form bound states or resonances. So the heavy quarks play an important role to stabilize the multi-quark systems. There were many theoretical studies of hidden-charm pentaquarks [5][6][7][8][9], especially the prediction of narrow N * and * resonances with hidden charm above 4 GeV by using the coupled-channel unitary approach [5,6], and the systematical investigation of possible hidden-charm molecular baryons with components of an anti-charmed meson and a charmed baryon within the one boson exchange model [7].
Very recently, the LHCb Collaboration observed two pentaquark-charmonium states in the J/ψ p invariant mass spectrum of 0 b → J/ψ K − p [10]. One is P c (4380) with a mass of 4380±8±29 MeV and a width of 205±18±86 MeV, and another is P c (4450) with a mass of 4449.8 ± 1.7 ± 2.5 MeV and a width of 39 ± 5 ± 19 MeV. The preferred J P assignments are of opposite parity, with one state having spin 3 2 and the other 5 2 . A lot of theoretical work has been done to explain these two states. In Ref. [11], the current experimental progress and theoretical interpretations of the states were reviewed. Chen et al. [12] interpreted these two hidden-charm states as the loosely bound c (2455)D * and * c (2520)D * molecular states by using the boson exchange model, and gave the spin parity J P = 3 2 − and 5 2 − , respectively. In Ref. [13], a Bethe-Salpeter equation approach was used to studied theD * c andD * c interactions, and then P c (4380) and P c (4450) were identified asD * c andD * c molecular states with the spin parity J P = 3 2 − and 5 2 + , respectively. A QCD sum rule investigation was performed, by which P c (4380) was suggested to be aD * c hiddencharm pentaquark with J P = 3 2 − and the P c (4450) was proposed as a mixed hidden-charm pentaquark ofD * c and D * c with J P = 5 2 + [14]. Also a coupled-channel calculation was performed to analyze the 0 b → J/ψ K − p reaction and gave support to a J P = 3 2 − assignment to P c (4450) and to its nature as a molecular state mostly made ofD * c andD * * c [15]. In Ref. [16], Meißner and Oller suggested that P c (4450) was almost entirely a χ c1 p resonance, coupling much more strongly to this channel than to J/ψ p. Kubarovsky and Voloshin [17] showed that the observed P c resonances are composites of J/ψ and excited nucleon states with the quantum numbers of N (1440) and N (1520) within a simple "baryocharmonium" model. Moreover, some people proposed various rescattering mechanisms to show that the P c (4450) state might arise from the kinematical effect [18,19]. Besides, Burns [20] explored the phenomenology of the P c (4380) and P c (4450) states, and their possible partners. Several intriguing similarities were also discussed in Ref. [20], which suggested that P c (4450) was related to the X (3872) meson. Thus, different models may give different descriptions for the resonance structures. Clearly the quark level study of these two pentaquark-charmonium states is interesting and necessary.
It is well known that the nuclear force (the interaction between nucleons) is qualitatively similar to the molecular force (the interaction between atoms). This molecular model of nuclear forces, the quark delocalization color screening model (QDCSM) [21,22], has been developed and extensively studied. In this model, quarks confined in one nucleon are allowed to delocalize to a nearby nucleon and the confinement interaction between quarks in different baryon orbits is modified to include a color screening factor. The latter is a model description of the hidden color channel-coupling effect [23]. The delocalization parameter is determined by the dynamics of the interacting quark system, and it thus allows the quark system to choose the most favorable configuration through its own dynamics in a larger Hilbert space. The model gives a good description of nucleon-nucleon and hyperonnucleon interactions and the properties of the deuteron [24][25][26][27][28]. It is also employed to calculated the baryon-baryon scattering phase shifts and the dibaryon candidates in the framework of the resonating group method (RGM) [29][30][31].
In this work, the resonating-group method (RGM) is employed to study the possible hidden-charm molecular pen-taquarks in QDCSM, and the channel-coupling effect are considered. Extension to the bottom case is straightforward and is also included in the present work. The structure of this paper is as follows. After the introduction, we present a brief introduction of the quark model used in Sect. 2. Section 3 is devoted to the numerical results and discussions. A summary is presented in the last section.

The quark delocalization color screening model (QDCSM)
The details of the QDCSM used in the present work can be found in Refs. [21][22][23][24][25][26][27][28][29][30][31]. Here, we just present the salient features of the model. The model Hamiltonian is where S i j is quark tensor operator; Y (x) and H (x) are standard Yukawa functions ( [32] and reference there in); T c is the kinetic energy of the center of mass; α ch is the chiral coupling constant; determined as usual from the π -nucleon coupling constant; α s is the quark-gluon coupling constant. In order to cover the wide energy range from light to heavy quarks one introduces an effective scale-dependent quarkgluon coupling α s (μ) [33]: where μ is the reduced mass of two interacting quarks. All other symbols have their usual meanings. Here, a phenomenological color screening confinement potential is used, and μ i j is the color screening parameter. For the light-flavor quark system, it is determined by fitting the deuteron properties, N N scattering phase shifts, N and N scattering phase shifts, respectively, with μ uu = 0.45, μ us = 0.19 and μ ss = 0.08, satisfying the relation, μ 2 us = μ uu * μ ss . When extending to the heavy quark case, there is no experimental data available, so we take it as a adjustable parameter. In the present work, we take μ cc = 0.01 ∼ 0.0001 fm −2 and μ uc is obtained by the relation μ 2 uc = μ uu * μ cc . All other parameters are also taken from our previous work [31], except for the charm and bottom quark masses m c and m b , which are fixed by a fitting to the masses of the charmed and bottom baryons and mesons. The values of those parameters are listed in Table 1. The calculated masses of the charmed and bottom baryons and mesons are shown in Table 2.
The quark delocalization in QDCSM is realized by specifying the single particle orbital wave function of QDCSM as a linear combination of left and right Gaussians, the single particle orbital wave functions used in the ordinary quark cluster model, Here s i , i = 1, 2, . . . , n are the generating coordinates, which are introduced to expand the relative motion wavefunction [23]. The mixing parameter (s i ) is not an adjusted one but determined variationally by the dynamics of the multi-quark system itself. This assumption allows the multiquark system to choose its favorable configuration in the interacting process. It has been used to explain the cross-over transition between hadron phase and quark-gluon plasma phase [35].

The results and discussions
Here, we investigate the possible hidden-charm molecular pentaquarks with Y = 1, For the negative parity states, we calculate the S-wave channels with spin S = 1 2 , 3 2 , and 5 2 , respectively; and for the positive parity states, we calculate the P-wave channels with spin S = 1 2 , 3 2 , and 5 2 , respectively. All the channels involved are listed in Table 3. In the present calculation, we only consider the hidden-charm molecular pentaquarks which consist of two S-wave hadrons. The channel-coupling effects are also taken into account. However, we find there is not any bound state with the positive parity within our calculations. There may exist other molecular structures, which contain excited hadrons, such as χ c1 p resonance [16], J/ψ N (1440), J/ψ N (1520) [17] and so on, which are out of range of present calculation. In the following we only show the results of the negative parity states.
First, the effective potentials between two hadrons are calculated and shown in Figs. 1, 2 and 3, because an attractive potential is necessary for forming a bound state or resonance. The effective potential between two colorless clusters is defined by is the energy of the system at the separation s of two clus- Table 3 The channels involved in the calculation ters, which is obtained by the adiabatic approximation. As mentioned in Sect. 2, a phenomenological color screening confinement potential is introduced in our model. For the multi-quark systems with heavy quark, because no experimental data is available, so we take three different values of μ cc (μ cc = 0.01, 0.001, 0.0001), to check the dependence of our results on this parameter. For the J P = 1 2 − system ( Fig. 1), one sees that the potentials are all attractive for the channels N η c , N J/ψ, c D, c D * and * c D * . For the channels c D and c D * , the potentials are repulsive, so no bound states or resonances can be formed in these two channels. However, the bound states or resonances are possible for other channels due to the attractive nature of the interaction between two hadrons. The attraction between * c and D * is the largest one, followed by that of the c D * channel, which is a little larger than that of the c D channel. In addition, the attraction of N η c is almost the same as that of N J/ψ, which is the smallest one during these five attractive channels. Comparing figures (a), (b), and (c) in Fig. 1, we also find that larger values of μ cc give rise lower energy, although the variation is not very significant.
For the J P = 3 2 − system (Fig. 2), similar results to that of I J P = 1   In order to see whether or not there are any bound states or resonances, a dynamic calculation is needed. The resonating group method (RGM), described in more detail in Ref. [36], is used here. Expanding the relative motion wavefunction between two clusters in the RGM by gaussians, the integrodifferential equation of the RGM can be reduced to an algebraic equation, the generalized eigen-equation. The energy of the system can be obtained by solving the eigen-equation. In the calculation, the baryon-meson separation (|s n |) is taken to be less than 6 fm (to keep the matrix dimension sufficiently small, so that we can still manage this case).
For the J P = 1 2 − system, the single-channel calculation shows that both c D and c D * are unbound, which agree with the repulsive nature of the interaction of these two channels. For the N η c and the N J/ψ channels, the attractions are too weak to tie the two particles together, the calculation shows that they are also unbound. While, due to the stronger attractions, the obtained lowest energies of c D, c D * and * c D * are below their corresponding thresholds. The binding energies of these three states are listed in Table   4, in which 'ub' means unbound. Here we should mention how we obtain the mass of a hidden-charm molecular pentaquark. Generally, the mass of a molecular pentaquark can be written as M the. = M the. 1 + M the.   Table 6 The masses (in MeV) of the hidden-charm molecular pentaquarks with all channels coupling and the percentages of each channel in the eigen-states are qualitatively similar to the conclusion of Ref. [5,6], in which they predicted two new N * states (the c D molecular state N * (4265) and the c D * molecular state N * (4415)) in the coupled-channel unitary approach. Meanwhile, the chiral quark-model calculation also supported the existence of the S-wave c D bound state [37]. At the same time, we also do a channel-coupling calculation. In this work, two kinds of channel coupling are performed. The first one is the coupling of three closed channels ( c D, c D * and * c D * ). The results, the lowest three eigenenergies and the percentages of coupling channels for the three eigen-states, are shown in Table 5. Taking the results of μ cc = 0.01 as an example, we can see that the main component of the lowest eigen-states is c D, ∼95.5%, and the energy is pushed down a little, compared with the singlechannel calculation, 4300-4296. The main component of the second lowest state is c D * with the percentage of 95.1%; and the main component of the third lowest state is * c D * with the percentage of 94.4%. The three eigen-energies are all smaller than the thresholds of the corresponding main channels, and they are stable against the change of the baryonmeson separations. The large percentage of the main component and the small change of energy lead us to infer that the channel coupling is very weak. However, these three closed channels can be coupled to the other four open channels, N η c , N J/ψ, c D and c D * . The results of this channelcoupling calculation are shown in Table 6. We obtain a stable state, the mass of which is lower than the threshold of N η c , and the main component of this state is N η c , with the percentage of 41.7%. This shows that N η c of J P = 1 2 − is bounded by channel coupling in our quark-model calculation, the energy is 3881-3884 (MeV). In addition, we also obtain several quasi-stable states, the masses of which are smaller than the thresholds of the corresponding main channels, but they fluctuate around the eigen-energies obtained in the three closed channel-coupling calculation.  Table 4. These three states also exist when they are coupled together, the masses and the percentages of each channel of the lowest three eigen-states are shown in Table 7. We can see that the mass of the first eigen-state is about 4362-4368 MeV and the main channel is * c D with the percentage of 91.0-95.5%; the mass of the second eigen-state is about 4445-4451 MeV and the main channel is c D * with the percentage of 96.2-98.5%; the mass of the third eigen-state is about 4551-4555 MeV and the main channel is * c D * with the percentage of 94.6-96.2%. From the above results, we find that the mass of the first eigen-state is close to the mass of the observed P c (4380), a pentaquark reported by the LHCb Collaboration. Therefore, in our quark-model calculation the main component of the P c (4380) is * c D with the quantum number − . In addition, the mass of the second eigen-state is close to the mass of another reported pentaquark P c (4450). Nevertheless, the opposite parity of the state to P c (4380) may prevent one from making this assignment. Moreover, all these closed channels can be coupled to the open channels N J/ψ and c D * . The results of the couplings of these five channels are shown in Table 6. There is a stable state, the mass of which is lower than the threshold of N J/ψ, and the main channel of this state is N J/ψ, with the ratio of 80.8-62.1%. This shows that the N J/ψ of J P = 3 scattering process. In fact, the possible existence of a nuclear bound quarkonium state was proposed more than 20 years ago by Brodsky et al. [38]; Gao et al. [39] also predicted the existence of the N φ bound state, which is very similar to the N J/ψ state; and the recent lattice QCD calculation also supported the existence of the strangenium-nucleus and the charmonium-nucleus bound states [40]. Therefore, searching for the N J/ψ resonance state is interesting work for the future. In addition, there are also several quasi-stable states in the full channel-coupling calculation, the masses of which are lower than the thresholds of the corresponding main channels and they fluctuate several MeVs around their central values. It is just the behavior of a resonance. The amplitude of the fluctuation can be taken as the decay width of the quasi-states. In a quark-model calculation, the decay width of the P c (4380) candidate is too small to match the experimental value. Further study is needed to check whether these eigen-states are resonance states in the N J/ψ and c D * scattering process and to calculate the widths of other decay modes. As mentioned above, by taking into account the channelcoupling effect, a bound state N η c is obtained for the J P = 1 2 − system; and another bound state N J/ψ is obtained for the J P = 3 2 − system. In these two systems, the coupling between the calculated S-wave channels is through the central force. In order to see the strength of this channel coupling, we calculate the transition potentials of these channels. Here, we take the result of the J P = 3  Fig. 4. Obviously, it is a strong coupling among these channels that makes Fig. 4 The transition potentials of channels N J/ψ, c D * , c D * , * c D, and * c D * , which are labeled channels 1, 2, 3, 4, and 5 respectively N J/ψ the bound state. The mechanism to form a bound state has been proposed before. Swanson proposed that the admixtures of ρ J/ψ and ω J/ψ states were important for forming X (3872) state [41], which was also demonstrated in Ref. [42] by Fernández-Caramés et al. The mechanism also applied to the study of H -dibaryon [43], in which the single channel is unbound, but when coupled to the channels N and , it becomes a bound state. The effect of channel coupling of the J P = 3 2 − system is the same as the one of the For the J P = 5 2 − system, it includes only one channel * c D * , and it is a bound state with the mass of 4512-4517 (MeV), which is a little higher than that of P c (4450). Although the width of the S-wave * c D * decaying to the D-wave N η c and J/ψ p (tensor interaction induced decay) is generally small, which can be used to explain why the width of P c (4450) is much narrower than that of P c (4380), J P = 5 2 − may not be a good candidate of P c (4450) because of the opposite parity of the state to P c (4380).
In the previous discussion, the hidden-charm molecular pentaquarks were investigated. We also extend the study to the hidden-bottom pentaquarks because of the heavy flavor symmetry. Here we take the value of μ bb = 0.0001. The numerical results are listed in Tables 8, 9, and 10. The results are similar to the hidden-charm molecular pentaquarks. For Table 8 The binding energies (in MeV) of the hidden-bottom molecular pentaquarks of I = 1    . However, it may not a good candidate of the observed P c (4450) because of the opposite parity of the state to P c (4380). Besides, the calculation is also extended to the hidden-bottom pentaquarks. The results are similar to the case of the hidden-charm molecular pentaquarks.
QDCSM, which was developed to study the multi-quark states, is an extension of the naïve quark model. As is well known, the quark model plays an important role in the development of hadron physics. The discovery of − is based on the prediction of the quark concept of Gell-Mann-Zweig. The naïve quark model of Isgur et al. gave a remarkable description of the properties of ground-state hadrons. On application to the excited states of hadron, hadron-hadron interaction, and multi-quark systems, extensions to the naïve quark model have to be made. Based on the different extension of the naive quark model, a proliferation of bound states or resonances are predicted. The recent progress of experiments on "XY Z" particles, P + c pentaquarks, and dibaryons such as d * [44][45][46] is encouraging. However, some precaution as regards the proliferation of quark-model bound states has to be posed. So far, there is no multi-quark state unambiguously identified by the experiments. For particular multi-quark state, there exist different points of view. For example, Vijande et al. studied the four-quark system ccnn in the constituent quark model by using different types of quark-quark potentials, and no four-quark bound states have been found [47], whereas the diquark-antidiquark picture was used by Maiani et al. to explain the state X (3872) [48]. Vijande et al. also searched for the doubly heavy dibaryons in a simple quark model, but no bound or metastable state was found [49], whereas H -like dibaryons with heavy quarks were proposed in Ref. [50]. More theoretical and experimental work are needed to distinguish the different extension of the quark model. The critical development of the quark model may be the unquenching quark model, where the valence quarks and real/virtual quark pair are treated equally. Caramés and Valcarce have studied the possible multi-quark contributions to the charm baryon spectrum by considering higher order Fork space components [51]. By incorporating new ingredients, the phenomenological quark model is expected to describe ordinary and exotic hadrons well.