Pentaquarks with hidden charm as hadroquarkonia

We consider hidden charm pentaquarks as hadroquarkonium states in a QCD inspired approach. Pentaquarks arise naturally as bound states of quarkonia excitations and ordinary baryons. The LHCb 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} pentaquark is interpreted as a ψ′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\psi '$$\end{document}-nucleon bound state with spin-parity JP=3/2-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^P=3/2^-$$\end{document}. The partial decay width Γ(Pc(4450)→J/ψ+N)≈11\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varGamma (P_c(4450)\rightarrow J/\psi +N)\approx 11$$\end{document} MeV is calculated and turned out to be in agreement with the experimental data for 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}. The 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} pentaquark is predicted to be a member of one of the two almost degenerate hidden-charm baryon octets with spin-parities JP=1/2-,3/2-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^{P}=1/2^-,3/2^-$$\end{document}. The masses and decay widths of the octet pentaquarks are calculated. The widths are small and comparable with the width of the 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} pentaquark, and the masses of the octet pentaquarks satisfy the Gell-Mann–Okubo relation. Interpretation of pentaquarks as loosely bound ΣcD¯∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varSigma _c\bar{D}^*$$\end{document} and Σc∗D¯∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varSigma _c^*\bar{D}^*$$\end{document} deuteronlike states is also considered. We determine quantum numbers of these bound states and calculate their masses in the one-pion exchange scenario. The hadroquarkonium and molecular approaches to exotic hadrons are compared and the relative advantages and drawbacks of each approach are discussed.


Introduction
Long anticipated heavy hadron states with hidden charm (and/or beauty) finally arrived in the recent years (see, e.g., the review in [1] and references therein) and are here to stay. Four-quark states with hidden charm were discovered first, and the heavy pentaquarks followed [2].
There are at least four possible scenarios for the dynamics of the LHCb pentaquarks. In the QCD inspired scenario one assumes that pentaquarks arise as a result of chromoelectric dipole interaction between a small quarkonium (charmoa e-mail: eides@pa.uky.edu b e-mail: Victor.Petrov@thd.pnpi.spb.ru c e-mail: maxim.polyakov@tp2.ruhr-uni-bochum.de nium) and a large baryon [3][4][5][6][7] (heavy quarkonium interaction with nuclei was considered in [8,9], see also references in [10]). Smallness of the quarkonium state is due to large masses of the heavy quarks. The strength of the quarkoniumproton interaction in this case is determined by in principle calculable quarkonium chromoelectric polarizability and by the proton energy-momentum density. The last one is normalized to the proton mass and is to a large extent modelindependent. Pentaquarks in this scenario look like atomlike systems with a small nucleus whose role plays the quarkonium state and light nucleon quarks that play the role of the atomic electrons. The characteristic feature of this scenario is that the pentaquark decay into a charmonium state with hidden charm and an ordinary baryon is by far the dominant mode of decay. Decays into states with open charm are strongly suppressed because they can go only via exchange by a heavy open charm meson. We will discuss the hadrocharmonium scenario in more detail below.
Molecular-like scenarios initiated in [11] rely on an analogy between heavy exotic hadrons with hidden charm and molecules. In this scenario charmed constituents of hiddencharm hadrons preserve their individuality and form bound states. There are two kinds of the molecular-like scenarios. In the first one charmed hadrons interact via exchange of light mesons and form hidden charmed pentaquarks with the binding energy at the level of hundreds of MeV (see review in [12] and references therein). In the second kind of molecular scenario initiated in [13][14][15] the binding energy is at the level of tens of MeV, and the heavy exotic hadrons are bound due to the one-pion exchange. This approach mimics the loosely bound deuteron. In the deuteron the S-wave one-pion central potential is not strong enough to bind the proton and neutron, and the much stronger noncentral tensor potential does not contribute to the S-wave. Binding in the deuteron arises because the tensor potential supports coupling between the Sand D-states. This mechanism was generalized for the case of the tetraquark mesons with hidden charm in [13,15]. Below we will develop this approach and apply it to the dynamical interpretation of the newly discovered pentaquarks. The generic feature of both molecular scenarios is the necessity to introduce a small distance repulsive cutoff. Technically this cutoff is needed to avoid collapse of the would be bound state. Physically, cutoff arises because due to the finite size both of the constituent hadrons and the exchanged mesons, the potential picture does not work at small distances. Similar cutoff is routinely introduced in nuclear physics, see, e.g., [16]. The distances between the open charm constituents in the molecular scenario are relatively large, what strongly impedes possible decays into final states with hidden-charm mesons like J/ψ. This seems to be a problem for this scenario, since the LHCb pentaquarks were discovered as bumps in the invariant mass distributions of J/ψ N . We will discuss the role of the cutoff and other features of the molecular scenario in more detail below.
One more popular idea is to treat heavy exotic hadrons with hidden charm as "true" tetra-and pentaquarks. This approach to heavy pentaquarks with hidden charm was initiated in [17] and developed further in numerous later publications, see, e.g., [18][19][20][21][22] and references in [12,23]. The idea is that the LHCb pentaquarks are diquark-diquark-antiquark bound states. The characteristic feature of this approach is that the hidden-charm pentaquarks arise as compact structures, more or less on par with the ordinary hadrons. The mere assumption about the diquark-diquark-antiquark structure of pentaquarks allows one to develop a rather reach phenomenology. Consideration of the color and flavor assignments for diquarks leads to the prediction of the flavor pentaquark multiplets [17]. The multiplet pattern in the diquarkdiquark-antiquark is qualitatively different than the one in the hadrocharmonium and molecular approaches and can serve as an experimental signature that allows to choose between different models. The SU (3) flavor symmetry was used to predict ratios of partial weak decay widths of bottom baryons to a pseudoscalar meson and pentaquark [20,21]. These predictions were further developed in the framework of the effective Hamiltonian approach to pentaquarks in [22] where the flavor SU (3) symmetry was amended by the heavy quark symmetry. A whole spectrum of new pentaquark states with definite properties arises in this approach. Also a very interesting set selection rules for weak decays of bottom baryons to pentaquarks is predicted. All these results could be used as a guideline in experimental searches for pentaquarks in bottom baryon decays.
In one more scenario pentaquarks are considered as molecular-like bound states of a "baryon" and a "meson" with an open color [24]. There are also suggestions in the literature that the LHCb results could be explained by some kinematical effects [25][26][27][28] without need for pentaquarks.
Below we will concentrate on the hadrocharmonium and molecular approaches to pentaquarks. We will describe in detail the interpretation of the LHCb pentaquarks as bound states of charmonium and the nucleon suggested in our previous paper [6], calculate the pentaquark masses and widths, and predict new pentaquark states. We will also present some new results in molecular approach and compare our results in hadrocharmonium and molecular approaches with the predictions of other authors.

Quarkonium-nucleon interaction
Different models were suggested for description of hadrons with hidden charm. Especially appealing is the hadroquarkonium approach to tetraquarks put forward in [3] (see also [4,5]). In [6] we applied this idea to the LHCb pentaquarks. The hadroquarkonium approach is based on the simple observation that interaction of a small size heavy quarkonium with other hadrons can be considered in the framework of the QCD multipole expansion [29][30][31][32], the role of the small parameter plays the ratio of the quarkonium size and the gluon wave length. For the Quarkonium-Nucleon interaction this ratio is just the ratio of the quarkonium and nucleon sizes. In the leading order approximation we have to consider emission (or absorption) of a chromoelectric dipole gluon by a heavy quark-antiquark pair. The color singlet pair goes into a color octet state after interaction with a dipole gluon and a second dipole interaction is needed to return it to the color singlet state. As a result in the leading approximation interaction of a heavy singlet quark-antiquark pair with other hadrons is described by the effective Hamiltonian (see, e.g., [10]) where E a i is the chromoelectric field with the absorbed strong coupling constant 1 α s , and α i j is the quarkonium chromoelectric polarizability where t a i are the SU (3) c color generators in the fundamental representation, r = r 1 − r 2 describes the relative positions of the quark and antiquark, and G is the quark-antiquark Green function in the color octet channel. In the nonrelativistic approximation the heavy quark-antiquark interaction is described by an effective Coulomb potential. Spins and coordinates decouple in a nonrelativistic bound state and the effective dipole Hamiltonian for quarkonium S-states reduces to where The perturbative quarkonium chromoelectric polarizability for an arbitrary nl states was calculated long time ago [33,34] (see also calculations for nS states in the large N c limit [31,32] and the recent calculation of the 1S polarizability for an arbitrary N c in [35]). Below we will need polarizabilities for the two lowest energy levels of nonrelativistic quarkonium [33,34] where m Q is the heavy quark mass, and a 0 = 3α s /4m q is quarkonium Bohr radius. We will also need the transitional 1S-2S polarizability that can be easily calculated in the large N c limit along the lines described in [31,32] Fitting the J/ψ and ψ masses, we obtain numerical values for the Coulombic polarizabilities Charmonium is not a Coulombic system and therefore one cannot expect quantitative agreement between the charmonium polarizabilities and their perturbative values. Transitional polarizabilities |α(J/ψ − ψ )| ≈ 2 GeV −3 and |α(Υ − Υ )| ≈ 0.66 GeV −3 for charmonium and bottomonium were extracted from the phenomenological analysis of the pionic decays ψ → J/ψπ π and Υ → Υ ππ [36]. Comparing these values with the perturbative results above we see that the perturbative calculations provide at best an order-of-magnitude estimates of the true polarizabilities. Below we will use perturbative polarizabilities for such estimates and for rough comparison of the relative magnitudes of the polarizabilities.
To describe interaction of a heavy quarkonium with a nucleon we need to calculate the expectation value of the chromoelectric field squared in Eq. (3) in a nucleon state. To facilitate this calculation we represent the chromoelectric field squared as a linear combination of the covariant gluon field strength G 2 μν and the energy density of the gluon field (zero component of the gluon energy-momentum tensor T G 00 ) where g 2 is the QCD coupling constant normalized at the scale of the quarkonium radius. It arises here because the QCD coupling constant describing interaction of a small chromoelectric quark-antiquark dipole with an external field is normalized at the quarkonium size (in our notation this coupling constant is swallowed by the field).
Exploiting the QCD scale anomaly we obtain where T μ μ is the trace of the QCD energy-momentum ten- is the leading coefficient of the β-function, and g s is the running strong coupling constant at the scale of the nucleon radius. Scale dependence of the coupling constant can be safely ignored for charmonium, but could become important for bottomonium. We temporarily omit the light quark masses in the trace of the energymomentum tensor, they will be accounted for later.
With the help of the representation in Eq. (9) the effective Hamiltonian in Eq. (3) reduces to the static nonrelativistic potential that describes interaction of a heavy quarkonium with the nucleon where T μ μ (x) and T G 00 (x) are the respective tensor densities inside the proton. The energy density T G 00 (x) carried by the gluons inside the proton cannot be determined unambiguously and is model-dependent. We make a natural assumption that it is proportional to the total proton energy density T G 00 (x) = ξ T 00 (x) [37]. Such assumption worked pretty well in the case of the pion. The factor ξ depends on the normalization point and is about 1/2 at Q 2 ∼ 1 GeV 2 [37]. We will assume this value in further calculations. It is convenient to represent the interaction potential in terms of the energy density and pressure, diag This effective potential has a simple interpretation. A pointlike quarkonium serves as a tool that scans the local energy density and local pressure inside the nucleon. It could happen that the size of quarkonium is not small enough in comparison with the size of the nucleon. In such case we will need to consider higher order terms in the QCD multipole expansion in order to improve description of the Quarkonium-Nucleon interaction.
The potential in Eq. (11) is to a large extent modelindependent, its normalization is determined by the total energy of the nucleon d 3 xρ E (x) = M N and the stability condition d 3 x p(x) = 0. Only the factor ν = 1 + ξ(bg 2 s /8π 2 ) in Eq. (12) cannot be determined from the first principles, and we use the phenomenological value ξ ∼ 1/2 to obtain ν ∼ 1.5.
We are going to use the interaction potential in Eq. (11) to explore possible bound states formed by heavy quarkonia states and the nucleon. With the known normalization of the potential and the nucleon radius we could proceed in an almost model-independent way, 2 considering the potential as a potential well with the size of the nucleon. Instead we will use the local energy density ρ E (x) and pressure p(x) that were computed in the χ QSM in [38]. The potential constructed in this way automatically satisfies the normalization condition in Eq. (12) since the normalization condition for the energy density and the stability condition for the pressure hold due to equations of motion.

Pentaquarks as hadrocharmonium states
The LHCb pentaquarks were discovered in the analysis of the invariant mass distributions of J/ψ plus nucleon. A natural idea is that the pentaquarks arise as bound states of charmonium excitations and the nucleon. We use the nonrelativistic Quarkonium-Nucleon potential in Eq. (11) to explore this hypothesis. The nonrelativistic Schrödinger equation in the channels J/ψ + N and ψ + N has the form where μ, ψ b and E b are the reduced mass, wave function, and the binding energy, respectively. The chromoelectric polarizabilities for each channel are collected in Eq. (7). Due to the poor knowledge of polarizabilities we will vary them in a relatively wide region. Solving the Schrödinger equation Eq. (13) numerically we find that a bound state of J/ψ and the nucleon arises only when the polarizability reaches the critical value α = 5.6 GeV −3 . This value of polarizability is more than an order of magnitude larger then the perturbative α(1S) in Eq. (7), and we conclude that J/ψ does not form a bound state with the nucleon. Critical values of polarizabilities for the excited S states of charmonia are far below their perturbative values in Eq. (7) (see [6,[31][32][33][34] for the perturbative polarizabilities of the higher excited states of charmonia). Hence, such states do form bound states with the nucleon. 3 We concentrate on the lowest ψ N bound states in this paper.
A bound ψ N state with the mass of the P c (4450) pentaquark, the binding energy E b = −176 MeV, and the orbital momentum l = 0 is formed at α(2S) = 17.2 GeV −3 . This value of polarizability is well inside the error bars of the perturbative calculation of the α(2S) polarizability in Eq. (7). There are no other bound ψ N states at α(2S) = 17.2 GeV −3 .
Exploiting the uncertainty in our knowledge of the α(2S) polarizability we can also adjust it in such way as to match the light LHCb pentaquark. A bound ψ N state with the mass of the P c (4380) pentaquark, the binding energy E b = −246 MeV, and the orbital momentum l = 0 is formed at α(2S) = 20.2 GeV −3 . Again, there are no other bound ψ N states at this value of polarizability. An identification of this bound state with the P c (4380) pentaquark would mean that there are no heavier pentaquarks formed by ψ N .
Taking into account opposite parities of the observed LHCb pentaquarks [2] it is interesting to explore possible ψ N bound states with l = 1. Such state arises for the first time when polarizability reaches the value α ≈ 22.4 GeV −3 . One could try to identify this state together with a more tightly bound l = 0 bound state with the pair of the LHCb pentaquarks. The spin-parities 3/2 − for the lighter state and 5/2 + for the heavier one fit nicely the experimental data. However, the mass splitting between these states is about 300 MeV instead of the experimentally observed 70 MeV. The large mass difference between the rotational excitation and the ground state indicates that the moment of inertia of the bound state is small. This bound state moment of inertia is determined by the size of the binding potential. The binding potential is proportional to the nucleon energy density, and hence the same binding potential determines the nucleon moment of inertia. In the mean field picture of the nucleon its moment of inertia determines the energy of its rotational excitations which is about a few hundred MeV as can be seen from the N −Δ mass splitting. Due to the connection between the nucleon moment of inertia and the bound state moment of inertia we are compelled to conclude that the moment of inertia of the bound state is small. This explains large splitting between the bound states with different angular momenta. Another drawback of the scenario with two pentaquarks as l = 0 and l = 1 bound states is that it predicts that the heavier pentaquark with l = 1 has a larger decay width, what squarely contradicts the experimental data in [2]. Both due to the prediction of a too large mass splitting and an unre-alistic hierarchy of decay widths we reject the interpretation of the LHCb pentaquarks as l = 0 and l = 1 bound states of ψ N .
The hadrocharmonium interpretation of pentaquarks was tested in [7] in the framework of the Skyrme model. It turned out that the Skyrme model energy-momentum tensor densities lead to the same conclusions as the considerations in [6] which were loosely based on the chiral quark soliton model. This demonstrates that the hadrocharmonium interpretation of pentaquarks is robust and does not depend on the details of a particular nucleon model. New hadrocharmonium ψ Δ bound states with hidden charm, isospin 3/2 and masses 4.5 and 4.9 GeV were predicted in [7].
In summary, solving the Schrödinger equation we have found two theoretically acceptable values of polarizability that admit interpretation of either of the two LHCb pentaquarks as a ψ N bound state. Only one bound state exists at each value of polarizability, and, respectively, only one of the observed pentaquarks can be interpreted as a ψ N bound state. Experimentally the P c (4380) peak has a rather large width 205 ± 18 ± 86 MeV, whereas the P c (4450) peak is narrow with the width 39 ± 5 ± 19 MeV. To make a choice between the two possible hadrocharmonium interpretations of the LHCb pentaquarks we need to calculate the theoretical decay widths of the bound state solutions found above.

Partial width of the ψ N bound state
Interaction of heavy charmonia states with the nucleon is described by the nonrelativistic potential in Eq. (11). This potential is universal, only its overall strength, determined by the polarizability of the respective charmonia excitation, changes when we go from one charmonia state to another. The nonzero transition polarizability in Eq. (7) shows that there exists a similar nondiagonal potential that describes the transition J/ψ → ψ off the nucleon. Due to the coupling between the J/ψ N and ψ N channels the pentaquark that we found in the ψ N channel should arise as a resonance in the J/ψ N scattering channel. We are going to solve the twochannel scattering problem, find the resonance J/ψ N → J/ψ N scattering amplitude, and determine the width of the resonance by comparing this scattering amplitude with the standard Breit-Wigner expression.
The Hamiltonian for the two-channel nonrelativistic scattering problem has the form where The potentials V i j are obtained from the potential in Eq. (11) by substituting the respective polarizabilities from Eq. (7) instead of α.
Next we solve the scattering problem for the Schrödnger equation where E is the nonrelativistic J/ψ N energy in the center of mass frame (E = q 2 /2μ 1 , q is the relative momentum) with only the incoming plane wave ψ 1 = e i q·x in the J/ψ N channel different from zero. The transition potential V 12 is small and the perturbation theory treatment is sufficient for the scattering problem in Eq. (15). Due to coupling between the channels the incoming plane wave ψ 1 (x) leaks in the ψ N channel Here is the Green function in the ψ N channel (see Eqs. (14) and (15)). Near the resonance where E R is the resonance energy. The wave function ψ 2 (x) in Eq. (16) in its turn generates correction to the incoming plane wave ψ 1 (x), that near the resonance has the form where is the free Green function in the J/ψ N channel.
Calculating δψ 1 (x) near the resonance with the orbital momentum l at large r = |x| we obtain the wave function in the J/ψ N channel as a superposition of the incoming plane wave and the outgoing spherical wave where the resonance radial wave function R l (r ) is normalized by the condition ∞ 0 drr 2 R 2 l (r ) = 1. This wave function can be written in terms of the scattering amplitude f (θ ) (θ is the scattering angle) The scattering amplitude near the resonance has the standard Breit-Wigner form where Γ is the resonance partial decay width in the J/ψ N channel.
Comparing Eqs. (20) and (21) we obtain where q = √ 2μ 1 E R and j l (z) is the spherical Bessel function.

Phenomenology of charmonium-baryon bound states
We obtained above a candidate for the heavy LHCb pentaquark as a bound ψ N state that arises at α(2S) = 17.2 GeV −3 , and a candidate for the light LHCb pentaquark as a bound ψ N state that arises at α(2S) = 20.2 GeV −3 . Now we are in a position to calculate these bound state partial decay widths into the ψ N channel. Using the phenomenological value of the transitional polarizability α(2S → 1S) = 2 GeV −3 [36] we obtain partial widths at the level of tens of MeV for both bound states. We also made a rough estimate of the partial width for the decay of either of the ψ N bound states into J/ψ + N + π , and it turned out to be even smaller than the partial width into the J/ψ + N channel. The decays of these bound states into (anti)charmed meson + charmed baryon are strongly suppressed in this scenario, since such decays into open charm channels can go only via t-channel exchange by a heavy D-meson. Therefore the total width of each of the ψ N bound states is small, in the range of tens of MeV.
The potentials in Eq. (13) are spin-independent, so there are two degenerate bound states with J P = 1/2 − and J P = 3/2 − . The hyperfine splitting between these degenerate color-singlet bound states arises due to interference of the chromoelectric dipole E1 and the chromomagnetic quadrupole M2 transitions in charmonium. It can be described by the effective Hamiltonian where S j is the quarkonium spin, α and m Q are the same chromoelectric polarizability and the heavy quark mass as above, and only the nucleon matrix element of the product of chromoelectric and chromomagnetic fields requires calculation.
The strength of this interaction is determined by the chromoelectric polarizability and it is additionally suppressed by the heavy quark mass ∼ 1/m Q . A semiquantitative estimate of the hyperfine splitting produces a small value in the range of 5−10 MeV. Therefore we expect to find two almost degenerate pentaquark states with J P = 1/2 − and J P = 3/2 − and with the mass of the observed pentaquark 4450 MeV. It would be very interesting if the LHCb collaboration could check this hypothesis in their partial wave analysis.
Thus far we ignored flavor symmetry of ordinary baryons. Recall that the nucleon is a member of the baryon octet. The interaction potential in Eq. (11) is proportional to the matrix element of E 2 , and in the linear approximation in the quark mass it is one and the same for all members of the baryon octet. Therefore we should expect that all members of the baryon octet bind with ψ , and the respective pentaquarks also form an octet. Masses of these octet pentaquarks are just the sums of the constituent masses and the binding energies. The binding energy depends on the mass of the ordinary octet baryon B only through the reduced mass in the kinetic energy in the respective Eq. (13). Then the pentaquark octet mass splittings in the leading order in the ordinary octet mass splitting ΔM are (see the definition of the reduced mass μ 1 after Eq. (14)) We checked this result by solving Eq. (13) for each ordinary octet baryon and calculating the pentaquark binding energies (and their changes) directly. Both approaches lead to the same results.
We have also solved the two-channel scattering problems for J/ψ scattering off all members B of the ordinary baryon octet, found the respective ψ B resonances, and calculated their partial decay widths into J/ψ B. The results for the octet pentaquarks mass splittings and widths are collected in Table 1. All octet pentaquarks P B have very small decay widths into J/ψ + B. We expect that like P c (4450) they also have small total widths. The mass splittings between the octet pentaquarks in Table 1 are somewhat smaller then the mass splitting in the ordinary baryon octet as predicted by Eq. (25). Due to hyperfine splitting there are two almost degenerate pentaquark octets with J P = 1/2 − and J P = 3/2 − . With very good accuracy the states in the pentaquark octets satisfy the Gell-Mann-Okubo mass formula or, numerically, 4613 MeV≈ 4615 MeV.
Only the heavy narrow pentaquark with spin-parity 3/2 − finds a natural interpretation as a ψ N bound state in the scenario described above. Another explanation should be found for the wide P c (4380) pentaquark 4 with spin-parity 5/2 + . The hadrocharmonium approach predicts also bound states formed by the nucleon and other excited states of charmonium, besides ψ . Scanning the charmonium spectrum in search of a state that could bind with the nucleon to form a pentaquark with mass 4380 MeV and spin-parity 5/2 + we observe that χ c2 (3556) with J P = 2 + has necessary quantum numbers and mass.
To find out if χ c2 (3556) really binds with the nucleon we have to solve a dynamical problem. As a nonrelativistic heavy quark-antiquark excitation χ c2 (3556) is a P-wave state. The chromoelectric polarizability of this P-state is a two-index symmetric tensor α i j that can be calculated using Eq. (2). Then one can calculate the χ c2 (3556)-nucleon interaction potential starting with Eq. (1), like it was done above in the case of the S-wave charmonium states, and solve the respective bound state Schrödinger equation. An estimate of the perturbative polarizability tensor for the P-state shows that it has roughly the same magnitude as for the S-state. Thus we have every reason to expect that χ c2 (3556) forms a bound state with the nucleon with spin-parity 5/2 + and mass about 4380 MeV. This state could be a candidate for the observed P c (4380) pentaquark. Moreover, due to the smallness of the spin-spin interaction there should be also almost degenerate states with spin-parities 1/2 + , 3/2 + . These particles do not exhaust the reach spectrum of pentaquarks formed by the P-states of charmonia and the nucleon. We expect to find χ c0nucleon bound state with spin-parity 1/2 + , almost degenerate χ c1 -nucleon bound states with spin-parities 1/2 + , 3/2 + , and almost degenerate h c -nucleon bound states with spinparities 1/2 + , 3/2 + . All these states should be very narrow because there are no open channels for decays except decays into particles without hidden charm that are strongly suppressed in accordance with the Okubo-Zweig-Iizuka rule. Suppression of strong decays of these particles puts a question mark over the possibility to identify the χ c2 (3556)nucleon bound state with the LHCb P c (4380) pentaquark.

Are there Bottomonium-nucleon hadrocharmonium bound states?
The binding mechanism developed above could generate bound states of bottomonium and the nucleon. The perturbative polarizability in Eq. (5) (see also [31][32][33][34]) depends on the heavy quark mass, the running strong coupling constant, and the Bohr radius. Both the effective coupling constant and the Bohr radius for bottomonium are smaller than for charmonium, while the mass of the bottom quark is larger than the mass of the charmed quark. Taking into account the interplay of these effects we calculated perturbative polarizabilities for the lowest states of bottomonia where we used the bottom quark mass m b = 5105 MeV, the bottomonium Bohr radius a 0 = 3α s /4m b ≈ 0.1 fm, and the strong coupling α s (a 0 ) ≈ 0.6 in this calculations. We searched for bottomonia-nucleon bound states, solving the Schrödinger equation Eq. (13) for bottomonium. No Υ (1S)-nucleon bound state was found for a reasonable value of polarizability. The results for the Υ (2S)-nucleon system are inconclusive due to poor knowledge of polarizability. A bound state could exist but a much better handle on polarizability is needed to make a definite statement. The sizes of higher bottomonia excitations are comparable to the size of the nucleon and the dipole approximation at the root of our approach is not valid any more.

Meson exchanges, nuclei, and the deuteron
Molecular bound states of two charmed mesons as an explanation of certain states in the charmonium spectrum were suggested long time ago [11]. Nowadays molecular models of mesons and baryons with hidden charm are very popular, and there are numerous papers discussing this scenario, see, e.g., review in [12]. We would like to compare the characteristic properties of pentaquarks arising in the molecular approach with the properties of the hadrocharmonium pentaquarks. Binding potential in the molecular approach is due to the exchange of light mesons. It is modeled after the similar approach widely accepted in nuclear physics. Let us recall the basics of the meson picture of nuclear forces [16].
The forces between nucleons in nuclei are due to exchanges of light mesons: pions, η-, σ -, ρand ω-mesons. Together they generate the potential resembling the Van der Waals molecular potential. Attraction at large distances (≥ 1 fm) is due to the light pion, attraction at intermediate distances is described by σ -meson (or two-pion) exchange and the short distance repulsion is usually ascribed to the ρand ω-meson exchanges. The meson exchanges do not make sense when distances between the nucleons in nuclei become comparable to the sizes of the exchanged bosons and/or sizes of the constituents. A strong repulsion core at small distances is usually introduced in the potential. Its position is determined by the particle sizes and should be chosen in the range of 0.3-0.5 fm. This approach provides at least a qualitative description of nuclei. The nucleons in a typical nucleus are separated by the distances about 0.7 fm. At these distances the contribution of the light pion exchange to the interaction potential is strongly suppressed and is almost irrelevant in comparison with the σ, ρ, ω contributions. This is not the case for the loosely bound deuteron. The binding energy in the deuteron is very small, about 2.2 MeV, and the nucleons in the deuteron are separated by a relative distance about 2 fm. At such distances only the light pion contribution to the potential survives. Calculating the nucleonnucleon scattering amplitude in the nonrelativistic approximation we obtain the momentum space nucleon-nucleon potential where g π N N = 13.7 is the pseudoscalar nucleon-pion coupling constant, m π and M N are the pion and nucleon masses, and S i and T a are the nucleon spin and isospin operators, respectively. The potential in coordinate space is a sum of a central spin-spin potential and a tensor potential where S 12 (S 1 , S 2 , n) = 3(S 1 · n)(S 2 · n) − (S 1 · S 2 ), and formally Naively, one could hope that this one-pion exchange potential would be sufficient to describe the deuteron. This does not happen due to the problems at small distances. Both the δ-function contribution to the spin-spin potential and the singular 1/r 3 contribution to the tensor potential are unphysical, they arise from distances where the one-pion exchange makes no sense due to finite sizes of all particles. To get rid of unphysical short distance contributions one could introduce soft or hard core at small distances. Instead it is routine in nuclear physics to regularize the potential at small distances by inserting the dipole form factor [(Λ 2 − m 2 π )/(Λ 2 + q 2 )] 2 in Eq. (28), see, e.g., [16]. Naive insertion of this form factor smears the δ-function contribution in Eq. (30). Notice that the δ-function contribution to the spin-spin potential has the sign opposite to the sign of the Yukawa type contribution in Eq. (30). After such formal regularization a repulsive at large distances spin-spin potential turns into strong attraction at distances shorter than than the inverse regularization parameter. We consider this regularized δ-function contribution to the potential unphysical, and subtract it from the regularized potential. Then the regularized potentials in Eq. (30) have the form (compare [15,39]) where The functions Y (Λ, m π , r ) and Z (Λ, m π , r ) are nonsingular and positive (or negative) definite at all distances. As a result the regularized potentials are finite at zero and a repulsive (attractive) at large distances potential remains repulsive (attractive) at all distances. It turns out that the one-pion exchange allows quantitative description of the principal deuteron characteristics with any short distance modifications we just described [13,15]. There is a nontrivial mechanism at work. The attractive spin-spin potential vanishes in the chiral limit and is therefore suppressed by the factor m 2 π /M 2 N . It is not strong enough to bind the proton and neutron if the tensor potential is turned off. The tensor potential is nonzero even in the chiral limit and is thus much stronger. It couples S-and D-waves in the Schrödinger equation and due to this coupling a loosely bound deuteron arises. We obtain the experimental binding energy 2.2 MeV if we place an infinite wall at r 0 = 0.485 fm. The fraction of the D-wave squared in this case is about 7% and the deuteron root mean square (rms) radius is 1.98 fm. The regularized potentials in Eq. (31) reproduce the same binding energy at Λ = 800 MeV, see also [15,39,40]. With the regularized potentials the fraction of the D-wave squared is about 5% and the deuteron rms is 1.92 fm. In both cases the cutoff parameters have a reasonable magnitude, confining the one-pion potential to distances larger than 0.25-0.4 fm.
The one-pion exchange mechanism modified at small distances completely describes all possible nucleon-nucleon bound states. There are four different spin-isospin states of two nucleons with S = 0, 1, T = 0, 1. There is no bound state with S = T = 0 since the spin-spin potential in this case is repulsive and the tensor potential turns into zero. The spin-spin potential is attractive and the tensor one is nonzero in the state with the quantum numbers of the deuteron, S = 1, T = 0. The spin-spin potential for S = 0, T = 1 is attractive and coincides with the deuteron spin-spin potential, while the tensor potential in this case is zero. Finally, in the state with S = T = 1 the signs of both potentials are opposite to the deuteron potentials and are suppressed by the factor 3. The Schrödinger equation with the potential in Eq. (29) and any short range regularization discussed above has a bound state solution with the deuteron quantum numbers and the binding energy 2.2 MeV. There are no bound states with any other spin-isospin quantum numbers. There are also no loosely bound nucleon-antinucleon states. Onepion exchange potential changes the overall sign when one replaces one of the nucleons by an antinucleon. Then absence of the nucleon-antinucleon states follows from the previous analysis.
Modern nuclear-type potentials (see, e.g., [16]) are much more sophisticated than the primitive one-pion exchange and include exchanges by other light bosons. We do not need them for our discussion of the deuteron. As we have seen the main features of a loosely bound deuteron are due to the long distance part of the one-pion exchange potential. Quantitative description of the deuteron requires some kind of short distance cutoff whether we include exchanges by other mesons besides pion or not, and the nature of this cutoff is only obscured by other mesons.
Below we will consider applications of the one-pion exchange and light boson exchange mechanisms to the LHCb pentaquarks.

One-pion exchange and pentaquarks
The one-pion exchange mechanism was generalized for description of tetraquarks with hidden charm in [13][14][15]. The essence of this approach to the loosely bound tetraquarks is the interplay between the channels with different orbital momenta. With reasonable assumptions about the magnitude of the short distance cutoff new tetraquark states were predicted in this framework [15]. The state X (3782) that was discovered many years later [41] turned out to be one of these predicted states [42]. It is worthwhile to figure out if this nice mechanism could be applied for the description of pentaquarks, can we construct the LHCb pentaquarks from a charmed baryon and an anticharmed boson as a deuteronlike loosely bound states. Inspection of the charmed hadron spectrum shows that the sum of masses of Σ * c (2520) andD(1870) exceeds the mass of P c (4380) only by 10 MeV, and the sum of masses of Σ c (2455) andD * (2010) exceeds the mass of P c (4450) only by 15 MeV. This immediately suggests that the respective LHCb pentaquarks are deuteronlike loosely bound states. At first glance the one-pion exchange mechanism has a fair chance to support the necessary binding in both cases. However, one-pion exchange cannot bind Σ * cD because the ππ D vertex is banned by parity conservation. Therefore P c (4380) cannot be considered as a loosely bound deuteronlike Σ * cD state. 5 It remains to figure out if P c (4450) could be interpreted as a deuteronlike loosely bound Σ cD * state. Let us first obtain the momentum space one-pion potential for interaction of two arbitrary hadrons. The Goldberger-Treiman relationship g N A /F π = g π N N /M N allows to replace the coupling constant g π N N in Eq. (28) by the nucleon axial charge g N A and the pion decay constant F π = 92 MeV. The coordinate space nucleon-nucleon potential in terms of the axial charge is obtained from the expressions in Eqs. (30), (31) by the substitution The respective hadron-hadron potential is obtained from Eqs. (30), (31) by the substitution where the constants g H i A are axial charges of the respective hadrons. The axial charges of heavy charmed hadrons almost never can be derived from the experimental data. An estimate of these charges can be obtained with the help of the naive constituent quark model for light quarks as suggested in [44]. We do not expect these axial charges to be particularly accurate, but they would have at least correct signs and order of magnitude. In the leading nonrelativistic approximation the time component of the quark axial current j 5 μa =ψ q t a γ μ γ 5 ψ turns into zero, and only the spatial components of the axial current proportional to σ i survive. Then in the framework of the nonrelativistic constituent quark model the axial charge g H A of the hadron H is given by the relationship where |Ψ H is the heavy hadron state vector in terms of the light quark states, g q A is the light quark axial charge, s i and t a are the quark spin and isospin operators, and S i and T a are the respective heavy hadron operators. Summation goes over all light constituent quarks. It is easy to generalize this expression for transitional axial charges g H 1 H 2 A . Below we will assume that g q A = 1 (and gq A = −1 for antiquarks). We can make a more accurate estimate of the quark axial charge using the Goldberger-Treiman relationship and Eq. (33) for the nucleon axial charge. In this way we would obtain Taking into account inaccuracy of the nonrelativistic constituent quark model itself it does not make much sense to make a distinction between g q A = 1 and g q A = 0.81. Axial charges of the nucleon and some charmed hadrons are collected in Table 2.
Let us return to the one-pion exchange interaction between Σ c andD * . Unlike the case of the deuteron the tensor potential in this case does not commute with the total spin S = S Σ c + SD * , [S 12 (S Σ c , SD * , n), S] = 0. Only the total angular momentum J = L + S is conserved. The lowest Σ cD * bound state should be dominated by the S-wave, and seeking an interpretation for the P c (4450) pentaquark we start with the sector with J = 3/2. The tensor potential has nonzero matrix elements between S-and D-waves and the Σ cD * state with J = 3/2 is a superposition of three states |L = 0, S = 3/2 , |L = 2, S = 1/2 , and |L = 2, S = 3/2 . The Hamiltonian in the subspace with J = 3/2 and T = 1/2 has the form (it coincides with the pion contribution to the respective Hamiltonian in [45]) where μ is the Σ cD * reduced mass. The potentials (nonregularized and regularized, with the subtracted δ-function contribution) according to Eq. (30) and Eq. (31) are e −m π r 9πr 3 , In this calculation we used g Σ c A = 2/3, gD * A = 1/2 from Table 2 and the substitution The spin-spin S-wave potential is attractive and is suppressed by the factor 8/25 in comparison with the respective deuteron potential. Again, like in the deuteron case it is not strong enough to bind Σ c andD * if the tensor potential is turned off.
We looked for the bound state solutions with the Hamiltonian in Eq. (35). To cut off the singular behavior of the potential at short distances we one time amended the potential in Eq. (35) by an infinitely hard wall at small distances, and another time we used the regularized potential in Eq. (36). In a model with the wall at r 0 = 0.33 fm we find a Σ cD * bound state with J P = 3/2 − , T = 1/2 and the binding energy 14.7 MeV, exactly with the mass of the P c (4450) pentaquark. This is a deuteronlike state, the binding arises due to the nondiagonal tensor potential. The hard core radius r 0 = 0.33 fm is somewhat smaller than in the case of the deuteron, but is still not too small. The rms of the bound state is about 1.6 fm.  Table 2. We also repeated these calculations with the phenomenological axial charges, see, e.g., [45]. We again can obtain a bound state at the position of P c (4450) but now Λ = 2000 MeV, rms is 1.13 fm, and the fraction of the D-wave squared is about 10%. Dependence of the energy level on Λ with any choice of the set of coupling constants is rather steep, the binding energy changes by about 4 MeV when Λ changes by 100 MeV. We see that P c (4450) can be interpreted as a deuteronlike bound state, but requires fine tuning of the short distance regularization parameter Λ. Now, that we fixed Λ, it is natural to look for the Σ cD * bound states in the channels with J = 3/2, T = 3/2 and J = 1/2, T = 1/2, 3/2. It turns out that there are no Σ cD * bound states with other quantum numbers besides J P = 3/2 − , T = 1/2 for the values of Λ and/or the position of the hard wall determined above. A Σ cD * bound state with J P = 3/2 − , T = 1/2 and the binding energy 85 MeV was obtained in [46,47] on the basis of the one-pion exchange. It was identified with the P c (4380) LHCb pentaquark. The binding in [46,47] occurred due to the spin-spin part of the one-pion exchange potential with account only for the S-wave wave function. The one-pion potential in [46] was regularized by the dipole form factor with Λ = 2.35 GeV (Λ = 1.78 GeV in [47]). The sign of the unregularized spin-spin potential in [46,47] is opposite to the sign in Eq. (35) and in [45], and corresponds to a long-distance repulsion. The binding is due to the regularized δ-function contribution to the spin-spin potential, see Eq. (30) (compare with the regularized potential in Eq. (36), where the δ-function contribution is subtracted). The repulsive long-distance one-pion spin-spin contribution to the potential in [46,47] can be omitted without changing the results. The rms of the bound state in [46,47] is 0.42fm, what is by far too small to justify validity of the one-pion exchange approximation. 6 This rms value can be roughly estimated almost without calculations, simply from the well known asymptotic formula for the bound state wave func- , where E is the binding energy. We disagree with the sign of the one-pion potential in [46,47], do not accept the idea of binding due to the attractive smeared δ-function, cannot justify the dominant role of the one-pion exchange at the distances about 0.4 fm, and therefore cannot accept the interpretation of P c (4380) as a Σ cD * bound state due to the one-pion exchange.

Are there other deuteronlike pentaquarks?
In search for other loosely bound pentaquark states with a deuteronlike binding we considered possible Σ * cD * loosely bound states. One could expect to find bound states with the binding energy about 10 MeV and mass about 4520 MeV, slightly below M Σ * c + MD * = 4530 MeV. Like in the case of Σ cD * interaction only the total angular momentum is conserved. We consider first the sector with J = 5/2. The lowest Σ * cD * bound state should be dominated by the S-wave, and the tensor potential has nonzero matrix elements between S-and D-waves. Hence, the Σ * cD * state with J = 5/2 − is a superposition of four partial waves, |L = 0, S = 5/2 , |L = 2, S = 1/2 , |L = 2, S = 3/2 , and |L = 2, S = 5/2 . The Hamiltonian in the subspace with J = 5/2, T = 1/2 has the form The potentials (nonregularized and regularized, with the subtracted δ-function contribution) according to Eqs. (30) and (31) are Y (Λ, m π , r ), In this calculation we used g Table 2 and the substitution The spin-spin S-wave potential is attractive and is 1.5 times stronger than the respective potential in the case of Σ cD * . We looked for the eigenstates of the Hamiltonian in Eq. (37) with the hard wall at r 0 = 0.33 fm, exactly at the same position as in the case of Σ cD * above. Surprisingly, there are no shallow loosely bound states with the binding energy about dozens of MeV. Only a bound state with the binding energy 82 MeV and the mass close to the mass of P c (4450) exists. There is a steep dependence of the binding energy on the radius of the hard core. Change of this radius from r 0 = 0.33 fm to r 0 = 0.38 fm reduces the binding energy to 10 MeV and leads to prediction of a pentaquark state with J P = 5/2 − , T = 1/2 and mass 4520 MeV. We repeated these calculations with the phenomenological coupling constants (see, e.g., [45]) and the regularized potentials in Eq. (38). Again we observe a steep dependence of the binding energy on the regularization parameter Λ, see Table 3. At Λ = 1400 we obtain a state with the mass of P c (4450), J P = 5/2 − and T = 1/2. Recall that we have already found a Σ cD * bound state with the mass of P c (4450), J P = 3/2 − and T = 1/2 at Λ = 2000 MeV using the phenomenological coupling constants. Interpretation of the Σ * cD * bound state as P c (4450) could be preferable, since Λ is smaller than in the case of Σ cD * . On the other hand rms of the Σ * cD * bound state at Λ = 1400 MeV is only 0.78 fm, and the fraction of the D-wave squared is about 25%, see Fig. 2. The one-pion exchange mechanism is probably not dominant for the bound state with such rms. The steep dependence of the bound state mass on the magnitude of Λ makes prediction of the mass of this state not too reliable. A fair conclusion could be that the consideration above does not have too much predictive power.
We also looked for possible Σ * MeV. This is an almost ideal deuteronlike bound state, but again a steep dependence on Λ does not allow to make a reliable prediction.
The Σ * cD * one-pion exchange interaction was also considered in [46,47], where a bound state with the binding energy 70 MeV, J P = 5/2 − , T = 1/2, and mass 4450 MeV was obtained and identified with P c (4450). Like in the case of Σ cD * we disagree with the sign of the potential in [46,47], see discussion above.

Nuclear type potential and pentaquarks
We considered above the deuteronlike model of loosely bound pentaquarks where the long distance one-pion exchange amended by a short distance repulsive core plays the defining role. We have already mentioned that more sophisticated descriptions of nucleon-nucleon interactions include exchanges by other light mesons, see, e.g., [16]. We do not need to go into detail of these more sophisticated models, but let us try to include exchanges by σ , ρ, ω, and η in the most straightforward way. One can use the constituent quark model to estimate the meson-nucleon interaction constants, see, e.g., [48]. In the framework of this model, σ and ω inter-act with the light quark baryon number and ρ interacts with isospin. Exchange by a scalar σ is always attraction, by a vector ω it is always repulsion and by a vector ρ it depends on isospin. In the case of the deuteron the total effect of the ρ and ω exchanges is repulsion, that naturally explains the origin of the repulsive core at small distances. The regularization at small distances is still needed due to the singular behavior of the tensor potential at small distances. The deuteron problem with the nuclear type potential and phenomenological coupling constants was considered in [48]. It turned out that inclusion of the meson exchanges effectively did not change values of the deuteron parameters obtained in the one-pion exchange scenario and the value of the regularization parameter Λ also did not change. This justifies the one-pion binding scenario for the deuteron. Almost all papers on the molecular model of tetra-and pentaquarks adopt the nuclear point of view and construct nuclear type potentials that include exchanges by all light bosons, see, e.g., review in [12]. Like in nuclear physics each meson exchange is regularized at short distances by the dipole potential with a phenomenological cutoff parameter Λ. This potential is nonuniversal and the cutoff parameter is adjusted for each experimentally known tetra-or pentaquark.
Like in the case of the deuteron one way to test the reliability of the one-pion exchange scenario for pentaquarks considered above is to include in the potential exchanges by other mesons and to see if the characteristics of the bound states (rms radius, fraction of the D-wave squared, parameter Λ) would remain stable. We obtained above a Σ cD * (J P = 3/2 − , T = 1) bound state at Λ = 2000 MeV with mass 4450 MeV, rms 1.13 fm, the fraction of the D-wave squared about 10%. Now we amend the one-pion exchange potential by σ , ρ, ω, and η exchanges, and use the phenomenological coupling constants (but not the signs of the individual contributions to the potential and not the regularization without subtraction of the δ-function contribution) from [45]. Then the bound state at the position of P c (4450) arises at Λ = 1300 MeV with rms 1.46 fm, and D-wave squared fraction about 4%. The value of Λ changed significantly thus hinting that the one-pion mechanism is not too reliable for this bound state.
We have also reconsidered Σ * cD * bound states in a model with the nuclear type potential. Like in the one-pion scenario we obtain a J P = 5/2 − , T = 1/2 bound state with E b = −74 MeV, rms 0.71 fm, and D-wave squared fraction 15% at Λ = 1200 MeV. Again, like in the one-pion exchange scenario we have also found a second Σ * cD * bound state at the same Λ. This is J P = 1/2 − , T = 3/2 state with E b = −37 MeV, rms 0.83 fm, and D-wave squared fraction 1.2%. We observe a rather substantial change of parameters, especially of the second state, in comparison with the onepion exchange scenario.
We see that the one-pion exchange scenario is not stable with respect to inclusion of other light meson exchanges, and it is hard to insist that this is the dominant binding mechanism even for the loosely bound pentaquarks. Both the onepion exchange and nuclear type scenarios for pentaquarks suffer from steep dependence on the short distance regularization parameter Λ (or position of the hard wall at small distances). One can describe existing experimental data on pentaquarks with the help of nuclear type and/or one-pion exchange potentials choosing different values of Λ for different states.
There are also apparent problems with the pentaquark decays in the molecular approach. The pentaquarks were discovered in the invariant mass spectrum of J/ψ and the proton. So decay to J/ψ + N is the only reliably established pentaquark decay mode. It is hard to understand how this decay to the states without open charm can can give a substantial contribution to the total width in the molecular picture. The constituents with open charm preserve their individuality in the molecular picture and we expect the decays into states without open charm to be strongly suppressed. This presents a qualitative difference with the hadrocharmonium picture, where the charmed quarks are close to each other and decay into J/ψ + N should should give a significant contribution.

Discussion of results
We have developed a QCD based approach to dynamical interpretation of pentaquarks. In this approach pentaquarks arise as bound states of ordinary baryons and excited states of quarkonia. The binding is due to the chromoelectric dipole interaction between the quarkonia states and ordinary baryons. The strength of the quarkonium-baryon interaction is determined by the quarkonium state chromoelectric polarizability and the baryon mass. The interaction potential is proportional to the density of the baryon energymomentum distributions. In this approach we interpret the LHCb P c (4450) pentaquark as a bound ψ N state with spinparity J P = 3/2 − . We calculated its decay width into J/ψ N , Γ (P c (4450) → J/ψ + N ) ≈ 11 MeV, what is in rough agreement with the experimental data. The P c (4450) pentaquark in this approach turns out to be a member of a pentaquark flavor octet, similar to the octet of ordinary baryons. The interaction between the quarkonia states and ordinary baryons is spin-independent in the leading approximation, so there are two degenerate pentaquark octets with spin-parities J P = 1/2 − and J P = 3/2 − . This degeneracy is lifted by a small color-singlet spin-spin interaction due to the interference of the chromoelectric E1 and chromomagnetic M2 transitions in charmonium.
Experimentally acceptable spin-parities for the LHCb pentaquarks are (3/2 − , 5/2 + ), (3/2 + , 5/2 − ), and (5/2 + , 3/2 − ) [2]. With our assignment of spin-parity J P = 3/2 − to P c (4450) we have to assign J P = 5/2 + to the wide P c (4380). We cannot find a natural interpretation for such state as a ψ N bound state. We have discovered χ c0 N , χ c1 N , χ c2 N , and h c N bound states with spin-parity J P = 5/2 + in the same mass region as P c (4380). Unfortunately, these states have no open channels for decays into an ordinary baryon and a charmonium state with hidden charm. Hence, all strong decays should go via heavy quark-antiquark annihilation, and are thus strongly suppressed in accordance with the Okubo-Zweig-Iizuka rule.
From the phenomenological perspective the main problem with the hadrocharmonium approach is its apparent inability to describe the P c (4380) pentaquark. On the theoretical side, development of a reliable method to calculate charmonium polarizability is urgently required.
We also considered the one-pion exchange model for loosely bound pentaquarks. In the framework of this model the P c (4450) pentaquark could be interpreted as a Σ cD * bound state with J P = 3/2 − and T = 1/2. The P c (4380) pentaquark that did not find a satisfactory explanation in the hadrocharmonium approach, does not admit a description based on the one-pion exchange either. Another difficulty of the one-pion exchange model is connected with decays. Phenomenologically we should expect that the partial width Γ (P c (4450) → J/ψ +N ) is not small. This follows from the simple observation that the P c (4450) pentaquark was discovered in the invariant mass distribution of J/ψ N . The amplitude of this decay is proportional to the bound state wave function at zero, since to create J/ψ the charmed quark and antiquark should come closely together. On the other hand the pentaquark in the one-pion molecular model is loosely bound, the wave function is smeared over a relatively large region and as a result the bound state wave function at zero is relatively small. Suppression of the pentaquark J/ψ + N decay mode is an apparent difficulty of the one-pion exchange model.
A nuclear type potential, that includes exchanges of different light mesons allows to consider not necessarily loosely bound pentaquarks. In this approach the binding energy could be relatively large and the constituents could be at much shorter distances than in the one-pion exchange model. This approach is much more flexible than the deuteronlike mechanism. The problem of the nuclear type potential which it shares with the one-pion potential is the apparent lack of predictive power. One cannot describe even the already observed pentaquarks with one and the same short distance regularization. An adjustment of this parameter is required for each particular state.
Either of the scenarios above predicts a number of new pentaquark states. A hadrocharmonium LHCb pentaquark is a member of an SU (3) flavor octet with small mass splittings (see the discussion above). The situation with a molec-ular type pentaquark is less clear. The SU (3) partners of such pentaquark can fail to form a bound state. This is what happens in the case of the deuteron, where partners do not form due to the difference of the pion and kaon masses. The pentaquark binding energy is larger, constituents in the pentaquarks are more densely packed than in the deuteron, so the mass differences inside the meson octets are less important. Whether the molecular pentaquarks have the SU (3) partners remains an open question, that probably cannot be answered in our rough approximation.
It is interesting to explore if there exist molecular bound states of other heavy mesons and baryons due to the one-pion exchange mechanism. As a simplest possibility it was suggested some time ago [49] that Λ c (2940) could be a molecular pentaquark made ofD * and nucleon. We tested this suggestion quantitatively in the same approach as in the discussion of the P c (4450) pentaquark. It turned out that there is no bound state in the ND * channel. This happens because due to the isospin factors the attraction induced by the pions in the ND * channel is only 3/4 of the attraction between Σ andD * , while the nucleon is lighter than Σ c (the coupling constants in both channels are approximately equal). Still, molecular pentaquarks formed by heavier baryons with different quantum numbers and by other heavy mesons could exist.
Hadrocharmonium scenario also admits existence of pentaquarks formed by other baryons (for example, by the Roper resonance as suggested in [50]) and different cc states. The chromoelectric interaction is spin-independent, so the bound states in this case (if any) will come in multiplets splitted by hyperfine interaction. They also should form the SU (3) flavor multiplets.
Both the hadrocharmonium interpretation of pentaquarks and the molecular-like approach have their own drawbacks and advantages, and need further development. Experimental and theoretical research on pentaquark decay rates and branching ratios could help to discriminate between different models. We hope to address decays in the future.