Mass spectrum of the hidden-charm pentaquarks in the compact diquark model

The LHCb collaboration have recently updated their analysis of the resonant $J/\psi p$ mass spectrum in the decay $\Lambda_b^0 \to J/\psi p K^-$, making use of their combined Run 1 and Run 2 data. In the updated analysis, three narrow states, $P_c (4312)^+$, $P_c (4440)^+$,and $P_c (4457)^+$, are observed. The spin-parity assignments of these states are not yet known. We interpret these narrow resonances as compact hidden-charm diquark-diquark-antiquark pentaquarks. Using an effective Hamiltonian, based on constituent quarks and diquarks, we calculate the pentaquark mass spectrum for the complete $SU (3)_F$ lowest $S$- and $P$-wave multiplets, taking into account dominant spin-spin, spin-orbit, orbital and tensor interactions. The resulting spectrum is very rich and we work out the quark flavor compositions, masses, and $J^P$ quantum numbers of the pentaquarks. However, heavy quark symmetry restricts the observable states in $\Lambda_b$-baryon, as well as in the decays of the other weakly-decaying $b$-baryons, $\Xi_b$ and $\Omega_b$. In addition, some of the pentaquark states are estimated to lie below the $J/\psi p$ threshold in $\Lambda_b$-decays (and corresponding thresholds in $\Xi_b$- and $\Omega_b$-decays). They decay via $c \bar c$ annihilation into light hadrons or a dilepton pair, and are expected to be narrower than the $P_c$-states observed. We anticipate their discovery, as well as of the other pentaquark states present in the spectrum at the LHC, and in the long-term future at a Tera-$Z$ factory.


Introduction
Recently, the LHCb collaboration has presented an updated account of the resonant J/ψ p mass spectrum in the decay Λ 0 b → J/ψ p K − , based on the combined Run 1 and Run 2 data, adding up to 9 fb −1 [1]. In this analysis, which supersedes their earlier findings from 2015 [2], nominal fits of the data have been performed with an incoherent sum of Breit-Wigner amplitudes, which have resulted in the observation of three peaks, whose masses, decay widths (with 95% C.L. upper limits), and the ratio R, defined as JHEP10(2019)256 are given in table 1. The state P c (4450) + in the 2015 data [2], is now replaced by two narrow states, P c (4440) + and P c (4457) + . In addition, a third narrow peak, P c (4312) + , having the mass M = (4311.9 ± 0.7 +6. 8 −0.6 ) MeV, is also observed. The spin-parity, J P , assignments of the three narrow states, which are crucial to decipher the underlying dynamics of the pentaquark states, are not yet determined. The broad peak P c (4380) + from the earlier data [2] is neither confirmed nor refuted, as the current LHCb analysis is not sensitive to broad resonances. Hence, it is entirely conceivable that more P + c -like structures are present in the J/ψ p invariant-mass spectrum, anticipated in the compact pentaquark interpretation, similar to the excited Λ * -baryon spectrum in the K p channel, expected in the quark model and confirmed in data [3]. For the discussion of the 2015 LHCb data on pentaquarks and other multiquark hadrons and references to the earlier works see the reviews [3][4][5][6][7][8][9][10][11][12][13][14] and [15].
The new pentaquarks reported by LHCb [1] have triggered a surge in theoretical papers interpreting the three narrow resonances as loosely-bound hadronic-molecule states in various incarnations [16][17][18][19][20][21][22][23][24][25]. A more up to date list of references can be seen in the proceedings of a recent Workshop on Exotic Hadrons [26]. Indeed, in the mass region of interest, a number of kinematical thresholds are present, such as Σ + (E thr = 4459.9 ± 0.9 MeV) [3]. The masses of the observed resonances in the hadron molecule interpretation are essentially the sums of their respective hadronic constituents, which is striking. Thus, not only the three observed narrow pentaquarks can be accommodated, having the spin-parity and isospin (J P = 1/2 − , I = 1/2) and (J P = 3/2 − , I = 1/2), but several other states, also with a negative parity due to their assumed S-wave character, are predicted. These emphatic claims are at best tentative and await a confirmation of the spin-parity assignments of the three observed states. The nature of the P c (4312) + as a candidate for the Σ + cD molecular state was also analyzed based on the S-matrix principles [27] and it was found that the attraction in this system is not enough for a bound state.
We emphasize that the decays Λ b → Σ ( * ) c + X are anticipated to be suppressed by Heavy Quark Symmetry (HQS), due to the mismatch of the spectator diquark spin (a conserved quantity in the HQS limit) in Λ b -and Σ ( * ) c -baryons. This suppression is wellknown, but is often ignored in the literature on the hadron molecular approach for the new P c -states, though there are some rare exceptions, see, for example, [28]. There is an ample evidence of this suppression in the PDG tables, in which the non-leptonic transitions Λ b → Λ ( * ) c + X dominate and the Λ b → Σ ( * ) c + X transitions are rather sparse [3]. In the exact HQS limit, such decays are forbidden, but as the HQS is not an exact symmetry and is broken by power (1/m b ) and QCD (O(α s (m b ))) corrections, they are allowed but have reduced decay rates. It has important implications for the analysis of the P c -states in the hadron molecule approach as well. In particular, it implies that the branching ratios B(Λ b → P + c K − ), which in the hadron-molecule interpretations of the P c -states are supposed to be induced by the intermediate (Σ cD ( * ) ) K − states, followed by rescattering, are subject to the HQS-implied suppression. Since the ratio R in (1.1), which is the product B(Λ b → P + c K − ) × B(P + c → J/ψ p) in units of B(Λ b → J/ψ p K − ), is well-measured for all three states by LHCb (see  Table 1. Masses, decay widths (with 95% C.L. upper limits), and the ratio R, of the three narrow J/ψ p resonances observed by the LHCb collaboration in the Λ b → J/ψ p K − decay [1].
expected to be small in the hadron-molecule interpretation. It is relevant to point out that model-dependent 90% C.L. upper limits on B(P + c → J/ψ p) of 4.6%, 2.3%, and 3.8% for P + c (4312), P + c (4440), and P + c (4457), respectively, have been posted by the GlueX collaboration [29], assuming they have J P = 3/2 − (i.e., S-wave J/ψ p states). While this is currently not robust enough an argument against the P c 's being hadron molecules, due to inherent assumptions about the theoretical estimates of the photoproduction cross sections, but improved measurements will soon test this interpretation in the P c -photoproduction experiments [29,30]. The role of the photoproduction process γ p → J/ψ p in constraining the models of the new P c -states is noted in [31,32]. Another process to measure the P c production is the antiproton-deuterium collisions, in which case the charm quark and antiquark entering the charmonium state J/ψ or η c are the result of the proton-antiproton annihilation [33]. Thus, the internal consistency of the hadron molecule approach for the new P c states remains to be checked.
In the preceding Letter [34], two of us have argued that there also exists a prima facie case for the three narrow P + c -resonances to be considered as candidates for compact pentaquarks. The basic idea of this approach is that highly correlated colored diquarks play a key role in the physics of multiquark states [35][36][37], and they are at work in the underlying dynamics of the P + c -states. A brief account of the new pentaquarks was presented in [34] in the effective Hamiltonian framework based on constituent quarks and diquarks, using isospin and heavy quark symmetry. In this paper, we extend the compact diquark template to cover the complete SU(3) F -multiplets of compact hidden-charm pentaquarks in the lowest S-and P -wave states. Some of them can be searched for in the decays of the Λ bbaryon, but some others can only be reached in the decays of the other weakly-decaying b-baryons, Ξ b and Ω b .
The pentaquark dynamics depends upon how the five constituents, i.e., the 4 quarks and an antiquark, are structured. Since quarks transform as a triplet 3 of the color SU(3)group, the diquarks resulting from the direct product 3 ⊗ 3 =3 ⊕ 6, are thus either a color anti-triplet3 or a color sextet 6. Of these only the color3 configuration is kept, as suggested by perturbative arguments. Another justification is that the color sextet diquarks are heavier than the corresponding color anti-triplet ones, and can be integrated out for spectroscopic considerations. This remains to be quantified, which can be eventually done in lattice QCD. Both spin-0 and spin-1 diquarks are, however, allowed. In the case of a diquark [q q ] consisting of two light quarks (q , q = u, d, s), the spin-0 diquarks are believed to be more tightly bound than the spin-1, and this hyperfine splitting has JHEP10(2019)256 implications for the spectroscopy. For the heavy-light diquarks, such as [cq] or [bq], this splitting is suppressed by 1/m c for the [cq]-diquark or by 1/m b for the [bq]-diquark, and hence both spin configurations are treated at par. Thus, the constituents of the hiddencharm pentaquarks in the compact diquark model are [cq]3,c3, and [q q ]3, which make up a color singlet. However, it is still a three-body problem to solve and there are several dynamical possibilities to model their interconnections. We follow here the intuitive picture in which the heavier components form a nucleus and the lighter one is in an orbit around this nucleus, as it is energetically easier to excite light degrees of freedom. Among the three constituents, the light diquark [q q ]3 is the lightest. Moreover, its spin-parity quantum numbers are fixed due to heavy-quark symmetry constraints. So, we keep the light diquark as emerging intact in the b-baryon decays, and put it in the orbit for the P -wave states, with the heavier components, carrying a charm quark or a charm antiquark, acting as a nucleus of (an almost) static color source. This is shown in figure 1 and the color and flavor structure in this approach can be written as [[c3(cq)3] 3 (q q )3]. Such a description is closer to the doubly-charm tetraquarks, with the quark content [(cc) 3 (q q )3], and doublycharm antibaryons, having [(cc) 3q3 ], as all three systems have two charm (anti)quarks at the center [38][39][40][41][42][43][44][45][46][47][48][49]. We note that a colored diquark-triquark template for pentaquarks θ δ ≡ Q (q 1 q 2 )3 3 (Qq 3 )3, through the color-triplet binding mechanism, was first used by Lebed [50]. There is also the third possibility to introduce a model for the internal structure of a pentaquark [51,52] in which the light and heavy diquarks bind first into a tetraquark in a color-triplet 3 state which further interacts with the charm antiquark, [[(cq)3(q q )3] 3c3 ]. Which of the three internal structures discussed provides a realistic template for the bound charm-anticharm pentaquarks remains to be seen.
In this paper we use the L -S scheme for describing the angular structure of pentaquarks. The spin of the heavy [cq]-diquark is S hd = S c + S q with the allowed values S hd = 0, 1, where S c and S q are the spins of the charm-and light u-, d-or s-quark, respectively. Combining the heavy diquark spin with the spin Sc of the charm antiquark, we get the spin of the doubly-heavy triquark, S t = S hd + Sc, having the values S t = 1/2, 3/2. The spin of the light [q q ]-diquark is S ld = S q + S q , where S ld = 0, 1. The total spin of the pentaquark is the sum of the triquark and the light-diquark spins, S = S t + S ld with the set of values S t = 1/2, 3/2, 5/2. In the case of the orbitally excited pentaquarks, the orbital angular momentum L of the pentaquark is the sum of two terms: L t , arising from the triquark system, and L ld , which determines the relative motion of the light diquark around the doubly heavy triquark. The total orbital angular momentum L of the pentaquark is then obtained with the help of the momentum sum rules from quantum mechanics, i.e., it takes a value in the range L = |L t − L ld |, . . . , L t + L ld . After L and S are constructed, the total angular momentum J of the pentaquark follows from J = L + S. The parity of the states depends on the orbital angular momentum as P = (−1) L+1 . It is negative for the S-wave states and positive for the states with one unit of the orbital angular momentum.
Mass estimates are worked out in an effective Hamiltonian approach, which apart from the constituent quark and diquark masses, includes dominant spin-spin, spin-orbit, orbital, and tensor interactions. The effective Hamiltonian for the S-wave states in this formulation is the same as the one used earlier in the analysis of the 2015 LHCb data JHEP10(2019)256 on pentaquarks [51,52], but the P -wave (and higher orbital angular) states are packaged differently. In addition to this, we incorporate the tensor interaction in the effective Hamiltonian, which affect the P -states, and work out the detailed mass spectrum in this framework, using the parameters fixed from earlier studies of baryons and hidden-charm tetraquarks. The resulting pentaquark spectrum is very rich. However, imposing the spin conservation in the heavy-quark symmetry limit, we argue that only that part of the pentaquark spectrum is expected to be observed in Λ b -decays, in which the pentaquarks have a "good" light diquark, i.e., having spin S ld = 0. This reduces the number of anticipated pentaquark states in Λ b -decays greatly. The same is also true for the decays of the Ξ 0,− b -baryons which are from the same SU(3) F -triplet as the Λ b -baryon. The other weakly decaying b-baryon, Ω b , in which the light diquark has the spin-1, belongs to the SU(3) F -sextet. Its decays will yield pentaquarks from the SU(3) F -decuplets, i.e., have a spin-1 light diquark (S ld = 1) as a constituent [51,53]. Some of the unflavored pentaquark states are estimated to lie below the J/ψ p (even below η c p) threshold and they will decay via the cc-annihilation into light hadrons or a lepton pair (e + e − and µ + µ − ), giving rise to states narrower than the observed peaks, which can also be searched in the P + c → p + − modes in Λ b → P c K − decays. The same holds for the decays of some other pentaquark states, produced in the decays of the Ξ b -and Ω b -baryons, in which case the corresponding pentaquarks with one, two, or three strange quarks will give rise to narrow Λ 0 + − , Ξ 0 + − , and Ω − + − final states, respectively.
We also note that a detailed spectrum of the tetraquark and pentaquark states has also been worked out in the dynamical diquark model [54], using Born-Oppenheimer potentials calculated numerically on the lattice. Likewise, the P + c -states have been studied in different color-bindings, such as {(cc) ((ud) u) ; L = 0} [28] and as (cc) (uud) hadrocharmonia [55].
This paper is organized as follows. In section 2, we introduce the doubly-heavy triquark -light diquark model of pentaquark, and define the state vectors having a total angular momentum quantum number J by |S hd , S t , L t ; S ld , L ld ; S, L J , where S hd , S t , and S ld are the spins of the heavy diquark, doubly-heavy triquark and light diquark, respectively. The corresponding sets of the S-wave state vectors (with L t = L ld = L = 0) with the "good" (S ld = 0) and "bad" (S ld = 1) light diquarks are presented in tables 2 and 3, respectively. State vectors of the P -wave pentaquarks with the ground-state triquark (L t = 0, L ld = L = 1) and "good" light diquark with the spin S ld = 0 ("bad" light diquark with the spin S ld = 1), are given in table 4 (table 5). In section 3, we give the effective Hamiltonian used to work out the pentaquark mass spectrum. In section 4, analytical expressions for effective Hamiltonian matrix elements are presented taking into account the dominant spin-spin interactions in the heavy and light diquarks and in the hidden-charm triquark. For the P -wave states, additional contributions and mixings due to the orbital, spin-orbit and tensor interactions are included. In all the cases, we diagonalize the mass matrices and the analytical equations for all the S-and P -wave pentaquark masses are presented. Section 5 contains the values of various input parameters, with the constituent diquark masses given in table 6 and the spin-spin couplings, KQ Q and (K QQ )3, extracted from spectra of mesons and baryons, are given in table 7. Our predictions for the unflavored pentaquark masses having the quark flavors (c[cq][q q ]), with q, q , and q being u-and JHEP10(2019)256 Figure 1. A picture of pentaquarks in the diquark model involving the heavy diquark [cq]3 and charm antiquark (c)3, which form the triquark system, which combines with the light diquark [q q ]3 to make a color singlet pentaquark. The subscripts indicate that all three constituents are color anti-triplets; q, q , and q are light quarks each of which can be u-, d-, or s-quark. d-quarks, assuming the isospin symmetry, are given in table 8. Masses of the hiddencharm strange pentaquarks are shown in tables 10-14. In this section, a comparison of our approach with the dynamical diquark model [54] is also presented. We continue with a discussion of dominant decay channels of ground-state pentaquarks containing a scalar light diquark in section 6 and conclude in section 7. In appendix A, corrections due to the spin-spin interactions between the constituents of the doubly-heavy triquark and the light diquark are detailed; in appendix B, cumbersome mass derivations for the P -wave pentaquark states are presented, and in appendix C the details and results of the χ 2 -analysis of the orbitally-excited Ω * c -baryons are shown.
2 Doubly-heavy triquark -light diquark model of pentaquark In the pentaquark picture shown in figure 1, there are two color flux tubes, with the first stretched between the charm diquark and the charm antiquark from which the diquark is in the color-antitriplet state3. With the antiquark being also a color antitriplet state3, their product is decomposed into two irreducible representations,3 ×3 = 3 +6, from which the color triplet 3 is kept. This color-triplet triquark makes a color-singlet bound state with the light diquark through the second flux tube in much the same way as the quark and diquark bind in an ordinary baryon. This approach differs from the one [51,52] used previously to get the mass spectrum of the S-and P -wave pentaquark states with the spin-parities J P = 1/2 ± , 3/2 − , and 5/2 ± , where the two diquarks are assumed to be organized first into a color-triplet tetraquark state which afterwards interacts with the charm antiquark. The double-charm triquark system [[cq]c] 3 is expected to be dynamically similar to the double-charm antidiquark, {cc} 3 , as far as their color and masses are concerned. The role of the latter, forming an almost static source, has been discussed at great length in a number of papers for the doubly-heavy anti-baryons, such asΞ cc (=ccq), and tetraquarks T (ccqq ) [38,[41][42][43][44][45][46]. The doubly-heavy triquark and the doubly-heavy diquark also differ in that the former has a light quark and the latter none. However, under the assumption JHEP10(2019)256  that the diquark [cq]3 is bound, and the dynamics is essentially determined by the color configuration and masses, the two are expected not to deviate from each other, at least as a first approximation. With this in mind, we first detail the doubly-heavy triquark -light diquark picture, modify the effective Hamiltonian for the S-wave pentaquarks [51,52] by keeping the most relevant terms for the mass determination, and then extend it for the Pstates by including the orbital, spin-orbit and tensor interactions between the hidden-charm triquark and light diquark.
The effective Hamiltonian for the ground-state pentaquarks is described in terms of two constituent diquark masses of the heavy diquark m [cq] ≡ m hd and of the light one m [q q ] ≡ m ld , the spin-spin interactions between the quarks in each diquark shell, and spin-spin interactions between the diquarks. To these are added the constituent mass m c of the charm antiquark and its spin-spin interactions with each of the diquarks.
Using the quantum numbers introduced above, we specify the complete orthogonal set of basis vectors for the hidden-charm pentaquark states. As already stated, we define the state vectors having a total angular momentum quantum number J by The corresponding sets of the S-wave state vectors (with L t = L ld = L = 0) with the "good" (S ld = 0) and "bad" (S ld = 1) light diquark are presented in tables 2 and 3, respectively. Among the overall ten states with a fixed light-quark flavor content, there are five states with J P = 1/2 − , four states with J P = 3/2 − and one with J P = 5/2 − . In view of the propensity of these states, we note that in the heavy-quark symmetry limit, the production of the states with the quantum numbers shown in table 3 is forbidden in Λ b -or Ξ b -decays. Nevertheless, they are present in the pentaquark mass spectrum, and they are reachable in other processes.  Table 4. Spin-parity J P and state vectors of the P -wave pentaquarks with the groundstate triquark (L t = 0) and "good" light diquark with the spin S ld = 0. The horizontal line demarcates the spin S hd of the heavy diquark.

JHEP10(2019)256
J P |S hd , S t , L t ; S ld , L ld ; S, L J  Table 5. Spin-parity J P and state vectors of the P -wave pentaquarks with the groundstate triquark (L t = 0) and "bad" light diquark with the spin S ld = 1. The horizontal line demarcates the spin S hd of the heavy diquark.
For the orbitally excited states, one needs to specify which part of the pentaquark wave function is excited. In the triquark-diquark template used here, and shown in figure 1, the heavy triquark state consists of the charm diquark and charm antiquark. Since both are heavy, the triquark is an (almost) static system and, hence, L t = 0 is its most probable orbital quantum number. Thus, the orbital excitation is generated by the light diquark (i.e., L = L ld ). With this, we list the lowest-lying orbitally excited states (L = 1) with the "good" (S ld = 0) and the "bad" (S ld = 1) light diquark in tables 4 and 5, respectively. In total, there are 25 pentaquark states with a fixed light-quark flavor content which are divided into nine states with J P = 1/2 + , ten with J P = 3/2 + , five with J P = 5/2 + and one with J P = 7/2 + .
The pentaquark spectrum emerging from the underlying diquark picture is very rich.
Contrasting it with the current experimental situation, with only three known hidden-JHEP10(2019)256 charm pentaquark states P c (4312) ± , P c (4440) ± , and P c (4457) ± observed in Λ b -baryon decays, some selection rules have to be applied to restrict the number of observable pentaquarks. Since in the heavy-quark symmetry limit, the spin of the light quarks is conserved in heavy-baryon decays, and S ld = 0 in Λ b -baryon, the states with the light-diquark spin S ld = 1 are suppressed. This means that we need to consider only the states in tables 2 and 4, for the S-and P -wave pentaquarks, respectively. This reduces the number of states to three (S-wave) and seven (P -wave). The states with S ld = 1 can, however, be produced in Ω b -baryon decays or in prompt production processes.
A related issue is the presence of two identical light quarks in the Fock space of a hidden-charm pentaquark. Among the light quark q present in the triquark, and the light quarks q and q in the diquark, there are two identical quarks. Treated as free quarks, they are subject to Pauli's exclusion principle, affecting the resulting spectrum of hadrons. This, for example, would be the case in the hadroquarkonium picture of the hidden-charm pentaquarks. However, if the q q -pair is bound in a diquark, which is a boson, then this exclusion does not apply so long as the diquark is not broken into its fermionic constituents. We illustrate this situation by focusing on the discovery mode Λ b → P + c K − . This decay is induced at the quark level by the b [ud] → cc s [ud] + (uū) vac transition. In this, the [ud]-diquark, which has well-defined quantum numbers (spin S [ud] = 0 and isospin I [ud] = 0), is a boson. In the heavy-quark symmetry limit, assumed in our paper, the light degrees of freedom remain conserved, and the light [ud]-diquark emerges in tact from the decay. Indeed, this imposes restrictions on the quantum numbers of final states, as we have emphasized earlier. If this symmetry is broken, say, by 1/m b terms, then the u-and d-quarks in the [ud]-diquark are no longer bound, and in the example here, one has two identical quarks, the u-quark produced from the vacuum excitation and the u-quark from the break-up of the spectator diquark, and the Pauli's exclusion principle would apply. In the heavy-quark symmetry limit, we have the spin-1/2 u-quark from the vacuum excitation, and the spin-0 [ud]-diquark, and they are distinctly different building blocks. In this case, no exclusion restrictions from identical quarks apply, which is what we have used in building the pentaquark spectrum.

Effective Hamiltonian for pentaquark spectrum
We calculate the mass spectrum of pentaquark under the assumption that their underlying structure is given byc, [cq], and [q q ], with q, q , and q being any of the light u-, d-, and s-quarks. We also assume the isospin symmetry among the states. For this, we extend the effective Hamiltonian proposed for the tetraquark spectroscopy by Maiani et al. [56]. The effective Hamiltonian for the S-wave pentaquark mass spectrum can be written as follows: (3.1) The first term in the Hamiltonian (3.1) is related with the color triquark:

JHEP10(2019)256
where m c and m hd are the constituent masses of the charm antiquark and charm diquark, respectively. The last three terms describe the spin-spin interactions in the charm diquark and between the diquark constituents and the charm antiquark. Among the three spin-spin couplings (K cq )3, Kc q , and Kc c , the spin-spin interaction inside the diquark (K cq )3 is argued to be the dominant one [56].
The second term H ld in the Hamiltonian (3.1) contains the operators responsible for the spin-spin interaction in the light diquark and its interaction with the triquark: where m ld is the constituent mass of the light diquark, consisting of the light quarks q and q , and the contribution of the spin-spin interaction in this diquark to the pentaquark mass is determined by the coupling (K q q )3. The last term in H ld is responsible for all the possible spin-spin interactions between the constituents of the light diquark and doublyheavy triquark: Here, the coefficients with the tilde differ from the ones introduced above as the former show the strengths of the spin-spin interactions inside the compact objects like diquarks and triquark while the later ones determine the strengths between the constituents of two objects, the heavy triquark and light diquark, and are strongly suppressed. The rationale of this is specific to the diquark model, in which the hadronic size of the diquarks is much smaller or compared to the overall size of the multiquark hadrons. Thus, only local spinspin interactions (within a diquark or triquark) are allowed. This then accounts for all possible spin-spin interactions and completes the content of the effective Hamiltonian (3.1) for the masses of the ground-state pentaquarks. A comment is in order: Despite the fact that this picture, which is physically more motivated, differs form the one used in [51,52], for the analysis and predictions of the ground-state pentaquarks, actually all the terms in the effective Hamiltonian (3.1) are the same. The general form of the effective Hamiltonian for the orbitally-excited pentaquark mass spectrum can be written as follows: In addition to the spin-spin interactions introduced for the ground-state pentaquarks described above, the extended effective Hamiltonian (3.5) includes terms explicitly dependent on the internal orbital angular momentum L t of the hidden-charm triquark and orbital momentum L ld of the light diquark relative to the triquark system. The corresponding Hamiltonian is called H L in (3.5). The terms relevant for the tensor interactions in each subsystem specified above are subsumed in H T in (3.5). As already stated, we assume that the triquark state is an S-wave (i.e., L t = 0). Thus, L = L ld . Inclusion of the triquark orbital angular momentum L t will be commented later on. The terms in H L which contain JHEP10(2019)256 the orbital angular momentum operator L are as follows: where quantities A t , A ld , and B parametrize the strengths of the triquark spin-orbit, lightdiquark spin-orbit and orbital momentum couplings, respectively. The last term in (3.5) represents the tensor interaction among the heavy triquark and light diquark: where R determines the position of the light diquark relative to the heavy triquark and b is the strength of the tensor interaction. The inclusion of a possible excitation in the triquark system (L t = 0) results into additional terms in the effective Hamiltonian: (3.8) The structure of ∆H L which includes the angular momentum operator L t in the triquark system is similar to (3.6): where quantities Ac, A hd , and B t parametrize the strengths of the charm antiquark and diquark spin-orbit and triquark orbital momentum couplings, respectively. As the triquark contains charm quark and antiquark, one expects a suppression of Ac, A hd , and B t by the c-quark mass in comparison with the parameters A t , A ld , and B entering (3.6). In the numerical calculations, we set them to zero. The second term in (3.8) represents the tensor interaction among the charm diquark and antiquark in the triquark: where r determines the position of the charm diquark relative to the charm antiquark and b t is the strength of this interaction. In the phenomenological analysis of the hiddencharm orbitally excited Y -tetraquarks, it was shown [57] that a similar strength entering the tensor interaction between the charm diquark and antidiquark is consistent with zero. Following [57], it is reasonable to assume that b t in (3.10) is also small, and hence we neglect this interaction in the analysis reported below by setting ∆H T = 0.
Returning to the discussion of the tensor contribution, the expression in eq. (3.7) is of the following general form: where S 1 and S 2 are the spins of two particles and n = R/R is the unit vector in the direction of R. This notation was used by us previously [57] in the analysis of the orbitally JHEP10(2019)256 excited (L = 1) Ω c -baryons and Y -tetraquarks. In both cases S 1 was identified with the spin of the diquark, S Q , and S 2 is either the spin of the charm quark, S c , in the Ω c -baryon, or the spin of the antidiquark, SQ, in the Y -tetraquark. The scalar operator Q(S 1 , S 2 ) in eq. (3.11) is expressed as the convolution of spins with the tensor operator: For further applications, we need the matrix elements of the operator N ij between the states with the same fixed value L of the orbital angular momentum operator L which can be obtained with the help of the identity from Landau and Lifshitz [58]: It is obvious that this matrix element is trivial when the two-particle quantum state is in the S-wave, L = 0. In terms of Q(S 1 , S 2 ), eq. (3.7) for L = 1 can be expressed as follows and similarly ∆H T = b t Q(S hd , Sc), if b t = 0. In eq. (3.14), the mixed-spin Q-operator can be reexpressed in terms of the single spin ones: where S = S t + S ld . To calculate the single-spin Q-operators inside the brackets, it is necessary to account for appropriate commutation relations of the components of L and S X , where S X = S t , S ld , and S are the spins of the triquark, light-diquark, and the total spin of the pentaquark, respectively. Setting L = 1, one finds easily [59]: The matrix elements can be computed directly by applying the operators (L · S X ) to the products of states corresponding to the individual spins and angular momenta. More effectively, one can use Wigner 6j-symbols, as is customary for analogous cases in atomic and nuclear physics. An explicit example of their application to the analysis of the orbitally-excited Ω c -baryons and Y -tetraquarks is presented in the appendix of ref. [57]. So, to get the mass predictions one needs to specify a concrete basis of the orbitally-excited pentaquark states in which the mass predictions can be done. Hamiltonian presented above in section 3. Note that this is simpler for the Model II by Maiani et al. [56], but becomes more involved in the Model I [35], where additional mixings between the spins of (anti)quarks in compact shells are included. As the later couplings are suppressed in comparison with the spin-spin interactions inside the shells, we neglect the later mixings and restrict ourselves with the Model II.
Restricting ourselves to the light u-and d-quarks, and assuming isospin symmetry, the universal contribution entering all the pentaquark states is defined as M 0 , which is the sum of the constituent masses of the heavy and light diquarks and charm antiquark: The two terms left in eq. (3.2) are responsible for the contributions of the spin-spin interactions between the charm antiquark and the quarks inside the heavy diquark, which together form the triquark state. To calculate their impact, it is convenient to use the Wigner 6jsymbols, which allows us to describe a recoupling between three angular momenta, say j 1 , j 2 and j 3 , which are combined into the state with the total momentum J. Following the angular momentum sum rules in quantum mechanics, one combines two of the three momenta first and then adds the third momentum. Thus, there are three possibilities: in the the first, the momenta j 1 and j 2 couple together in the state with the angular momentum j 12 , which after combining with j 3 results in the state with total angular momentum J; in the second, the momenta j 1 and j 3 produce the state with the momentum j 13 and it is combined with j 2 to get the state with the momentum J; and in the third, the momenta j 2 and j 3 produce the state with the momentum j 23 which then combines with j 1 to obtain the state with the momentum J. Each way of summing up the momenta results in its own set of basis vectors. If, for example, the momentum j 2 is coupled initially to j 1 and one wants to recouple it to j 3 , there are two extra momenta j 12 and j 23 , each of which has its own set of basis vectors, |(j 1 , j 2 ) j 12 , j 3 ; J and |j 1 , (j 2 , j 3 ) j 23 ; J . The Wigner 6j-symbols describe the transformation from one basis to the other as follows [60]: (4.4) Here, Wigner 6j-symbols are represented by the curly brackets. To find the matrix elements of the spin-spin operators 2 Kc q (Sc · S q ) and 2 Kc c (Sc · S c ) in eq. (3.2) with JHEP10(2019)256 the help of the Wigner 6j-symbols, it is convenient to rewrite the pentaquark vector state |S hd , S t , L t ; S ld , L ld ; S, L J as |Sc, (S q , S c ) S hd ; S t and |Sc, (S c , S q ) S hd ; S t , respectively, (here, all the unnecessary fixed quantum numbers are omitted from the vector and the charm antiquark spin, Sc, is shown explicitly) and make the angular momenta recoupling as in eq. (4.4): The required matrix elements can be written in the form: where spins of all the quarks are replaced by their numerical values, Sc = S c = S q = 1/2. Up to the factors Kc q and Kc c , both equations coincide, so the contribution of these two terms entering the effective Hamiltonian to the mass formulae is proportional to the sum of the coupling strengths, Kc q + Kc c . The matrix elements of the operators in eq. (3.4) have in general six different coupling strengths. For the pentaquark having the light u-and d-quarks only, isospin symmetry results into relations among the couplings. This simplifies the spin-spin interaction (3.4) substantially. For the pentaquarks having one or more strange quarks, some of the relations no longer hold, as SU(3) F -symmetry is broken. Details of the matrix element calculation in this case are presented in appendix A. This contribution will be neglected in the numerical analysis as we are working in the Model II [56].
We apply the above formalism for the calculation of the mass spectrum beginning from the pentaquark states from table 2. From three states presented there, two states JHEP10(2019)256 with J P = 1/2 − mix due to the spin-spin interaction of the charm antiquark and the heavy diquark, and the third one with J P = 3/2 − remains unmixed. For the later state, the mass m S0 3 (here, the superscript denotes the S-wave pentaquark with the "good" light diquark, S ld = 0) is the average of the effective Hamiltonian over this state: where M 0 is defined in (4.1). The former two states with J P = 1/2 − , after sandwiching the effective Hamiltonian, yield the following (2 × 2) mass matrix: Diagonalizing this matrix yields the masses of the other two states: where two ratios r hd and r ld of the couplings are defined as: The averaged mass of the states with masses given in (4.11) and (4.12) has the following value:m 14) and the mass of the third state lies above by an amount We continue calculations of the mass spectrum taking the S-wave pentaquark states with the "bad" light diquark from table 3. From the seven states presented there, two states with J P = 1/2 − and the other two with J P = 3/2 − , both pairs having the triquark spin S t = 1/2, mix due to the spin-spin interaction of the charm antiquark and the heavy diquark. The other three states with the triquark spin S t = 3/2 remain unmixed due to this interaction but can mix through the spin-spin interactions between (anti)quarks in the two separated shells -the heavy triquark and light diquark. As mentioned earlier, such types of spin-spin interactions are suppressed and neglected in this analysis. For the later three states, enumerated according to their entries in table 3, the masses m S1 5,6,7 are the averages of the effective Hamiltonian over the states considered (here, the superscript denotes the S-wave pentaquark with the "bad" light diquark, S ld = 1):

JHEP10(2019)256
So, one can see that taking into account only the spin-spin interactions inside the light diquark and heavy triquark, these states are degenerated in mass. The pairs of states with J P = 1/2 − and J P = 3/2 − , after sandwiching the effective Hamiltonian, yield identical (2 × 2) mass matrices: (4.17) Diagonalizing the matrix, the states, being twice degenerate, have the masses: where r hd and r ld are defined in eq. (4.13). The averaged masses of the states with masses in (4.18) and (4.19) have the following value: and the mass gap between these average values and the massm S0 12 (4.14) is completely determined by the strength of the spin-spin interaction in the light diquark: The degeneracy is lifted if the spin-spin interactions between the triquark and the lightdiquark constituents are taken into account as suggested by Maiani et al. in the Model I [35].

P -wave pentaquarks
As it was discussed earlier in section 3, one needs to include the terms dependent on the internal angular momenta of the pentaquark system which can be written as additional terms (3.6) and (3.7) in the pentaquark effective Hamiltonian. Let the vector states |S hd , S t , L t ; S ld , L ld ; S, L J , appearing in tables 4 and 5, be denoted as |L, (S t , S ld ) S ; J and transformed to the state |(L, S t ) Jt , S ld ; J . After the transformation (4.4) is applied, we get: (4.23) The required matrix elements can be written in the form: where the value of the orbital angular momentum, L = 1, is substituted.

JHEP10(2019)256
Next, we denote the state |S hd , S t , L t ; S ld , L ld ; S, L J as |L, (S ld , S t ) S ; J and transform it to the state |(L, S ld ) J ld , S t ; J . After the transformation (4.4) is applied, we get: (4.25) With this, the matrix elements can be written as: where again the value of the orbital angular momentum, L = 1, is substituted. We calculate the mass spectrum of the P -states, starting with the pentaquark states with the "good" light diquark from table 4. From the seven states presented there, two states with J P = 1/2 + and the other two with J P = 3/2 + , both pairs having the triquark spin S t = 1/2, mix due to the spin-spin interaction of the charm antiquark and heavy diquark. The other three states with the triquark spin S t = 3/2 remain unmixed due to this interaction but can mix through the spin-spin interactions between the (anti)quarks entering two separated shells -the heavy triquark and light diquark, as also discussed earlier for the S-wave states. As mentioned earlier, such types of spin-spin interactions are suppressed and neglected in this analysis. From the terms entering the angular-momentum dependent parts (3.6) and (3.7) of the effective Hamiltonian, one obtains the spin-orbit, orbital, and tensor contributions to matrix elements. For the three states, enumerated as the fifth, sixth and seventh according to their rows in table 4, the masses m P 0 5,6,7 are the matrix elements of the effective Hamiltonian over these states (here, the superscript denotes the P -wave pentaquark with the "good" light diquark, S ld = 0): These states are non-degenerate due to the triquark spin-orbit interaction, i.e., A t = 0. The pair of states with J P = 1/2 − , after sandwiching the effective Hamiltonian, yields the following (2 × 2) mass matrix: (4.28) while the masses of the pair with J P = 3/2 − are determined by the matrix: (4.29)

JHEP10(2019)256
Diagonalizing these matrices, we get the masses: where r hd and r ld are defined in eq. (4.13).
We continue the calculations of the mass spectrum for the pentaquark states with the "bad" light diquark from table 5. From the eighteen states presented, there are six states with J P = 1/2 + , seven with J P = 3/2 + , four with J P = 5/2 + , and the last one with J P = 7/2 + . Except for the J P = 7/2 + state, all the others mix due to the spinspin interactions, but also because of the spin-orbit and tensor terms (3.6) and (3.7) in the effective Hamiltonian. For the last state, enumerated as the eighteenth according to table 5, the mass m P 1 18 is the average of the effective Hamiltonian over this state (here, the superscript denotes the P -wave pentaquark with the "bad" light diquark, S ld = 1): As mentioned earlier, there are four states with spin-parity J P = 5/2 + which can be divided into two pairs according to their mixing mechanism: the fifth and the tenth states in table 5 mix due to the spin-spin interaction while for the fifteenth and the seventeenth states, their mixing is determined by the spin-orbit and tensor interactions. We start from the first of the mentioned two states. Their mass matrix is as follows: Diagonalizing the matrix, we obtain the masses: For the remaining two states, the mass matrix has the form:

JHEP10(2019)256
So, the masses of the pentaquarks are two eigenvalues of this matrix: Further, there are seven states with the spin-parity J P = 3/2 + which can be separated according to the triquark spin S t = 1/2 (four states) and S t = 3/2 (three states) as shown in table 5. The former four states mix under the spin-spin, spin-orbit and tensor interactions while three states with S t = 3/2 mix due to the spin-orbit and tensor interactions only. The details of the mass derivation can be found in appendix B; here we present explicit equations for the masses only: where λ 1,2,3 are determined in eq. (B.10).
To derive the expressions for the masses of the four states with the spin-parity J P = 3/2 + and the triquark spin S t = 1/2, one needs to find the eigenvalues of a non-diagonal symmetric (4 × 4) matrix. The explicit form of the mass matrix and the details of calculations can be found in appendix B. The set of masses is as follows: Finally, there are six states with the spin-parity J P = 1/2 + which can be divided according to the triquark spin S t = 1/2 (four states) and S t = 3/2 (two states) as shown in table 5. The four states mix under the spin-spin, spin-orbit and tensor interactions while the two states with S t = 3/2 mix through the spin-orbit and tensor interactions only. The mass matrix of the later two states is as follows:

JHEP10(2019)256
The pentaquark masses are the two eigenvalues of this matrix: For the last four states with the spin-parity J P = 1/2 + and triquark spin S t = 1/2, the mass matrix is again a non-diagonal symmetric (4 × 4) matrix. It is written explicitly in appendix B and details of calculations of its eigenvalues can be also found there. The resulting set of masses is presented below: where With these expressions, the analytical calculations of the pentaquark mass spectrum in Model II is done and we present the numerical estimates of the pentaquark masses in the next section.

Input parameters
Working within the Constituent Quark-Diquark Model [56], input parameters are the masses of the constituents, charm quark and two diquarks, spin-spin couplings, and other parameters related with the orbital or radial excitations. In the present paper we analyse the ground-state pentaquarks and their first orbital excitations only. To estimate the constituent quark masses, there are two possibilities: either extract them from the masses of  Table 6. Constituent diquark masses (in MeV) used for calculating the hidden-charm pentaquark mass spectrum. Here, q denotes u-and d-quarks and isospin symmetry is assumed. The brackets [. . .] and {. . .} indicate the antisymmetric (S = 0) and symmetric (S = 1) spin configurations of the quarks in the diquark, respectively. (Taken from [62] for light diquarks and from [15] for the heavy-light ones).
known mesons or from baryons [15]. Quark masses obtained from the baryon spectrum are larger by typically 50 MeV [15]. This can be exemplified by the charm quark mass, which is estimated as m m c = 1667 MeV from the D-meson spectrum, as opposed to m b c = 1710 MeV from the charm baryon masses, yielding a mass difference of 43 MeV [15]. In particular, with the m b c value as an input, predictions for the charm baryon masses were obtained [61] which differ from the experimentally observed masses [3] by about 10 MeV, which may be viewed as an error on the input charm-quark mass. We use m b c in the numerical analysis. Based on the same arguments, we accept m b q = (362 ± 10) MeV and m b s = (540 ± 10) MeV as the light and strange quark masses, respectively [15].
The diquark masses are presented in table 6. In contrast to quarks, which have spin S q = 1/2, diquarks, being composite objects, have two possible spin configurations from which the antisymmetric one corresponding to the diquark spin S = 0 is energetically more favorable. Both configurations are allowed if the flavors of the quarks are different, indicated by right brackets, [. . .], in table 6 (for S = 0), and by curly brackets, {. . .}, (for S = 1). For the diquarks, having two quarks of the same flavor, only the spin-symmetric configuration is allowed by Bose statistics, as indicated for the case of the diquark with two s-quarks in table 6. In the limit of the exact isospin symmetry, the light u-and d-quarks are massdegenerate and denoted by q in  [15]. To get the diquark-mass error estimates, one should compare baryon masses predicted within this model and in experiment. Based on the analysis presented in [61], in which the predicted and measured mass differences for unflavored (i.e., having u-or d-quarks), strange and charm baryons do not exceed 15 MeV, we take it as a measure of the uncertainty in diquark masses.
The spin-spin couplings, KQ Q and (K QQ )3, extracted from the spectra of mesons and baryons, respectively, are presented in table 7. These values are taken from [15] except for (K cc )3, which is from [61]. The factor two difference in Kc c and (K cc )3 reflects the fact that the one-gluon exchange assumption is used to relate these couplings. Concerning the (K cq )3-coupling, it follows from table 7 that the value (K cq )3 = 15 MeV extracted from the charm baryon spectroscopy is numerically substantially smaller than all others. In fact, this is the spin-spin interaction between the charm quark and the two quarks from a light diquark, but not between the quarks within the charm diquark [ Table 7.
Spin-spin couplings, KQ Q and (K QQ )3, extracted from the spectra of mesons and baryons, respectively (borrowed from [15] and for cc from [61]). Here, q ( ) denotes u-and d-quarks and isospin symmetry is assumed. estimated as (K cq )3 = 67 MeV [15,56]. The later value was obtained from the Constituent Quark-Diquark Model analysis of the hidden-charm exotic mesons in which the underlying structure is assumed to be a charm diquark and charm antidiquark. In our numerical estimates, we use the later value. The spin-spin coupling obtained from the charm baryons can be used to determine (K cq )3 in eq. (3.4), associated with the spin-spin interactions between the charm quark in the doubly-heavy triquark and quarks from the light diquark. The same is also true for the spin-spin couplings between the charm and strange quarks, for which the value (K cs )3 = 2 MeV [57, 59] 1 follows from the analysis of the narrow orbitallyexcited Ω * c -baryons recently observed by the LHCb [63] and Belle [64] collaborations. To understand the uncertainty of (K cs )3, we performed the χ 2 -analysis based on the masses of these baryons [63,64], which yields (K cs )3 = (2.01 ± 0.20) MeV. The details of this analysis are in appendix C. This yields 10% error on (K cs )3, which we assume as an uncertainty on all the spin-spin couplings. The pentaquark masses involve the ratios of the couplings (4.13) which for the input values are evaluated as: r hd = 0.73 and r ld = 1.07. According to the Model II, these ratios should be of order one which explicitly demonstrates that the strengths of all possible spin-spin interactions inside the diquarks and triquark are approximately the same.
In estimating the P -wave pentaquark mass spectrum, values of the couplings in the orbital angular momentum term and the spin-orbit ones are required. In the previous analysis [51,52], the spin-orbit coupling was taken to be A P = 52 MeV, which follows from the analysis of the vector hidden-charm Y -tetraquarks. Updated analysis of these states in [57], which takes into account in addition the tensor interaction between the diquark and antidiquark, results in the smaller value, a Y = (22 ± 3) MeV. A similar analysis [57,59], performed for the orbitally-excited Ω c -baryons, yields even smaller values, 2 a 1 = (13.45 ± 0.13) MeV and a 2 = (12.94 ± 0.36) MeV, for the spin-orbit interactions of the {ss}-diquark and charm antiquark, respectively. The uncertainties in these values result from the χ 2 -analysis of the Ω * c -baryon masses [63,64] and the details of this analysis are also presented in appendix C. The errors in these couplings are typically a few percent and are substantially smaller than the relative error from the Y -tetraquark spectrum which is approximately 14%.

JHEP10(2019)256
For P -wave pentaquarks, the spin-orbit couplings are responsible for the interactions of the triquark and the light diquark, denoted by A t and A ld , respectively, which are, in general, different. However, based on the observation that the couplings of the spin-orbit interaction in Ω * c -baryons are close to each other, a 1 a 2 13 MeV, one can use the same approximate relation A t A ld = 13 MeV in the pentaquark system but, as we shall show later, this is not supported by the data on the newly observed pentaquark states [1] for the plausible spin-parity assignment, as discussed in [34].
For the mass spectrum of orbitally-excited pentaquarks with the "bad" light diquark, the tensor coupling b is required. It was extracted from the spectrum of the orbitallyexcited Ω c -baryons and Y -tetraquarks, yielding b Ωc = 13.5 MeV [57,59] 3 and b Y = (−136 ± 6) MeV [57], respectively. Thus, this coupling is quite different in the charm baryons and in the hidden-charm tetraquarks and, moreover, the central values are of opposite signs in these hadrons. Here, more experimental input is needed. To be definite, we use b = (13.3 ± 0.5) MeV, which is obtained from the orbitally-excited Ω c -baryons rather than from the Y -tetraquarks.
The last parameter to be specified is the orbital coupling B, several estimates of which are available from different hadrons. Some representative values are: B(cc) = 457 MeV [57], obtained from charmonia, B(Ω c ) = 325 MeV, from Ω c -baryons, B(Y ) = 362 MeV and B(Y ) = 505 MeV, from the hidden-charm Y -tetraquarks. 4 In the analysis [51,52], the coupling B was derived from the experimentally observed mass of the P c (4450) + having the (preferred) spin-parity assignment J P = 5/2 + . The orbital coupling was determined from the P c (4450) + mass to be B(P c ) = 220 MeV for the Model II [56]. With the observation of three narrow pentaquarks in the updated LHCb data, the experimental picture is changed drastically. Following the interpretation of these peaks as suggested in [34], which assumes that two resonances are orbitally-excited states and the lowest mass peak is the S-wave state, it is possible to determine the orbital coupling B phenomenologically, discussed in the next subsection.

Predictions for the hidden-charm unflavored pentaquarks
To get an estimate of the mass spectrum of the ground-state hidden-charm unflavored pentaquarks, it is enough to know the diquark and charm quark masses and spin-spin couplings, discussed in detail in the previous subsection. For the P -wave pentaquark masses, the values of the orbital angular momentum coupling B, the spin-orbit couplings, A t and A ld , and the tensor coupling b are required. Some of these couplings can be determined from the measured masses of the observed resonances [1]. In this subsection we consider two possible assignments, both assume that the heavier states, P c (4440) + and P c (4457) + , are the P -wave pentaquarks, while the lowest mass P c (4312) + can be either S-or P -wave state, and show that the second scenario is physically unacceptable. Assuming that the 3 The χ 2 -analysis of the Ω * c -baryons presented in appendix C gives slightly smaller value bΩ c = (13.30 ± 0.48) MeV. 4 The values presented are determined by two possible assignments of the vector tetraquarks discussed in [57].

JHEP10(2019)256
P c (4312) + pentaquark is an S-wave state, we predict the masses and J P quantum numbers of yet unmeasured resonances. We discuss first the pentaquarks P c (4440) + and P c (4457) + . As they replace the former narrow P c (4450) + state with the preferred spin-parity J P = 5/2 + , we tentatively assign this spin-parity to one of the observed two states, P c (4457) + . In this case, the lighter partner, P c (4440) + , most probably has the spin-parity J P = 3/2 + . In the Model II by Maiani et al. [56], the mass splitting of the two positive-parity states is related to the spin-orbit coupling A t of the triquark: where the error is obtained by adding the experimental errors on the masses in quadrature. From (5.1), the value of the coefficient A t follows immediately: It is not surprising that A t is found numerically small, as the doubly-heavy triquark is almost static. The third narrow state P c (4312) + can have several J P assignments. Identifying it with the J P = 3/2 − state, one can work out the mass difference between this state and the heavier pentaquarks, P c (4440) + and P c (4457) + . This is determined by the orbital B and the triquark spin-orbit A t couplings. The strength of the latter is already known from the P c (4440) + and P c (4457) + mass splitting (5.2), and the mass difference, say, between P c (4312) + and P c (4457) + , allows us to read off the strength of the orbital interaction: Alternatively, assuming that the spin-spin couplings and the constituent (quark and diquarks) masses are known, the strength of B can also be determined from the masses of P c (4440) + and P c (4457) + only, as follows:  [57]. With these parameters, the mass of the third pentaquark M = 4240 MeV with J P = 3/2 − , also shown in the solid box, is somewhat lower than the mass of the observed P c (4312) + peak, but is still in the right ball-park. The second possibility is to assign the lowest mass state P c (4312) + with the one having the spin-parity J P = 3/2 + and mass M = 4361 MeV, or with J P = 1/2 + and M = 4351 MeV. Both predictions are rather close to the observed pentaquark mass. In the case of J P = 1/2 + assignment for P c (4312) + , Model II predicts for the mass difference: which does not depend on B and A t but is determined by the spin-spin couplings from the doubly-heavy triquark. If we fix the spin-spin coupling (K cq )3 = 67 MeV, eq. (5.6) can be solved analytically with the roots: The first solution is unacceptable, as it implies Kc q + Kc c → ∞. The second solution in (5.7) yields a negative and, hence, unphysical value for r hd , if the ratio ρ lies in the interval: 3/2 < ρ < 2. With the pentaquark mass difference in (5.6) and (K cq )3 = 67 MeV, one obtains the value ρ 1.92, which is in the unphysical interval (the corresponding value of the spin-spin coupling ratio r (2) hd −10.4). The situation does not improve if we assign J P = 3/2 + to the third state. In this case, it has a higher mass than the J P = 1/2 + pentaquark by 3A t 10 MeV, thereby the r.h.s. of (5.6) gets an additional factor −3A t . There are again two solutions as in (5.7) but in r hd becomes positive but too large, r (2) hd 16, which is again unacceptable. So, the assignment of the third state, P c (4312) + , with the orbitally-excited pentaquark is not tenable in the Model II.
It should be noted that the assignment of all the three observed states with the P -wave pentaquarks was also discussed in ref. [65]. In contrast to the one-gluon-exchange model, from which the effective Hamiltonian used in our paper follows, the effective Hamiltonian in [65] is based on a confining potential and an attractive flavor-spin interaction, resulting from the underlying chiral dynamics. Under the assumption that the flavor-spin interaction is dominant, the three LHCb states are then identified with the lowest-mass orbitallyexcited states. Moreover, in this model, the lowest S-wave state lies above the P -wave states. This obviously contradicts the predictions of our model discussed above.
We remark that the assignments we have discussed satisfy the heavy-quark symmetry, which means that the spin S ld = 0 of the light diquark in Λ b -baryon is conserved in the decay to the charged pentaquark and K − -meson. That is why, only pentaquarks with the "good" light diquark are considered. As the heavy-quark symmetry is not exact, pentaquarks with  Table 8. Masses of the hidden-charm unflavored pentaquarks (in MeV) and their comparison with the results presented in [51,52]. For the P -wave pentaquarks, the values of the orbital, spin-orbit and tensor couplings are taken from table 9.
S ld = 1 can also be produced in the Λ b -baryon decays, but being suppressed by 1/m b , they are rather unlikely.
Masses of the hidden-charm unflavored pentaquarks are presented in table 8, and compared with the results obtained in [51,52]. For the P -wave pentaquarks, the values of the orbital, spin-orbit and tensor couplings are taken from table 9. The threshold for the observed pentaquarks in the P + c → J/ψ + p decay mode is M thr J/ψ p = m J/ψ + m p = 4035.17 MeV [3]. With the masses given in table 8, there are two states with the "good" light diquark (S ld = 0), the J P = 1/2 − state with the mass M = 3830 MeV and the J P = 1/2 + state with the mass M = 4030 MeV, which lie below the M thr J/ψ p threshold. Also, the third state having J P = 3/2 + , with the mass M = 4040 MeV, may also lie below M thr J/ψ p . There are also two states, J P = 1/2 − and J P = 3/2 − , with the "bad" light diquark (S ld = 1), having the degenerate masses M = 4026 MeV, which also lie below the JHEP10(2019)256 B A t A ld b 207 ± 20 3.4 ± 1.1 13.5 ± 0.4 13.3 ± 0.5 Table 9. Input values (in MeV) for the couplings required for the mass estimations of the P -wave pentaquarks. For the coupling B, the error is taken from (C.13).
M thr J/ψ p threshold. One of these states (with the mass M = 3830 MeV) is even below the threshold for the decay P + c → η c + p with M thr ηc p = m ηc + m p = 3922 MeV [3]. They will decay through the annihilation of the cc-pair into light hadrons or a lepton pair, and hence will be narrower the P + c -resonances observed. It is conceivable that a dedicated search of narrow states in the LHCb data with improved statistics may reveal the existence of these pentaquark states. It is interesting to note that the virtual J/ψ-meson decay will also lead to an off-shell dilepton pair. Hence, a search for a resonance in the p µ + µ − (as well as p e + e − ) may result into observing narrow structures. In addition, the other states shown in table 8 are also reachable in the Λ b → J/ψ p K − decay, in particular, our model predicts a peak with the mass M = (4150 ± 29) MeV below the lowest-mass observed pentaquark, P c (4312) + , and the other peak with a mass near 4350 MeV, being a combination of two J P = 1/2 + and J P = 3/2 + pentaquarks with close masses. The measured spectrum by the LHCb collaboration [1] has indications for such states, though more data are needed for a definite conclusion.

Comparisons with the dynamical diquark model
Before presenting the strange pentaquark spectrum, we compare some of our estimates in the non-strange sector with those obtained in the dynamical diquark model [54].
The common feature of our model and the dynamical diquark model is the assumption of point-like diquarks in the color-antitriplet representation, with the pentaquarks built from a diquark and a color-triplet triquark. However, the two approaches differ in their dynamical details and in assigning the flavors of the quarks in the diquarks and triquarks. In our model, we assume that the light diquark [q q ]3 is bound to the doubly-heavy triquark [c3 [cq]3] 3 and the orbital excitations in the pentaquark are due to the light-diquark orbital excitations (the doubly-heavy triquark is assumed to be in the S-wave). In the dynamical diquark model [54], both the diquark [cq]3 and the triquark [c3 [q q ]3] 3 are heavy objects, and orbital excitations are between them [54]. More importantly, in our analysis we take into account all sizable spin-spin, spin-orbit, orbital and tensor interactions between the charm antiquark and the diquarks, except for the spin-spin interaction between the constituents of the doubly-heavy triquark and the light diquark, which are assumed to be strongly suppressed. All these interactions are omitted in the current formulation of the dynamical diquark model and left for further study [54]. So, a lot of pentaquark states which differ in their masses in our model (see, for example, tables 2 and 4) are mass degenerate in the dynamical diquark model, which makes a numerical comparison rather difficult.
We also differ in the J P assignments of the states. In the dynamical diquark model [54], two possible assignments are suggested: the first one identifies the two pentaquark states

JHEP10(2019)256
with larger masses, P c (4440) + and P c (4457) + , with the states from the first orbitallyexcited multiplet called Σ + g (1P ) (see table 7 in [54]). With this assignment, there is no place in the mass spectrum for the third narrow state, P c (4312) + , observed recently by the LHCb collaboration. This is not surprising, as also remarked in [54], since the strength of the orbital excitation in this assignment is about ∼ 400 MeV for the states with the principle quantum number n = 1, in accord with the corresponding excitation energy in charmonia. The observed mass splitting ∼ 140 MeV is, therefore, impossible to accommodate with these assignment. So, the second possibility discussed in [54] is to put the lowest-mass state, P c (4312) + , in the first orbitally-excited multiplet Σ + g (1P ), and in this case the two other states P c (4440) + and P c (4457) + are the members of the radially excited multiplet Σ + g (2S). The suggested spin-parity of the P c (4380) + -state is J P = 5/2 + (this means that the bottom of the Σ + g (1P )-multiplet is fixed to the mass of the P c (4312) +pentaquark) while one of the P c (4440) + -and P c (4457) + -states has J P = 3/2 − . In our model, all three observed pentaquark states have the principle quantum number n = 1 and have the indicated parities: J P = 3/2 − for P c (4312) + and J P = 3/2 + and J P = 5/2 + for P c (4440) + and P c (4457) + , respectively. We are able to explain quantitatively the mass splitting in the higher mass pentaquarks by the spin-orbit interaction, and the difference between the lower-and higher-mass states, to a large extent, by the orbital interaction. The strength of the latter interaction was determined as ∼ 200 MeV, which is smaller by ∼ 100 MeV compared to charmonia, but still acceptable as the size of the pentaquark and, hence, the distance from the doubly-heavy triquark to the light diquark, is anticipated to be significantly larger than the hadronic size of the charmonium system.
Thus, despite the similarity in the general approach about the constituents (colored diquark and triquark) of the pentaquarks, the two model are different in details and lead to different spectra. Both models require estimates of the diquark-size effects on the spectroscopy, a point also discussed later in this paper.

Mass predictions for the strange hidden-charm pentaquarks
Including the s-quark, the hidden-charm pentaquark spectrum becomes much richer. Depending on the strange content of the pentaquarks, they can be classified into singly-, doubly-, and triple-strange pentaquarks and their SU(3) F -multiplets are in one to one correspondence with the ordinary and strange baryons [51]. The mass estimates given below assume mass-degenerate u-and d-quarks. The constituent masses of the quarks are accepted as specified in subsection 5.1. The diquark masses and spin-spin couplings are not SU(3) F -invariant, and their breaking is taken into account (see tables 6 and 7). For the P -wave pentaquarks, the orbital, spin-orbit and tensor couplings are assumed equal for all three light quarks and are taken from table 9.
The masses of the singly-strange pentaquarks are presented in tables 10 and 11 and compared with the results of [51,52]. The difference between the entries in the two tables lies in the s-quark attachment either to the heavy diquark or to the light diquark, i.e., the pentaquark structures are (c3 [

cs]3[qq ]3) and (c3[cq]3[sq ]
3), respectively, with q ( ) being u-or d-quark. One should also remember that the states presented in tables 10 and 11 with the same set of quantum numbers: spin-parity J P and isospin I, having strangeness  Table 10.
Masses of the hidden-charm strange pentaquarks (in MeV) with a strange-heavy diquark, i.e. having the structure (c3[cs]3[qq ]3), where q ( ) is u-or d-quark, and their comparison with the results presented in [51,52]. For the P -wave pentaquarks, the values of the orbital, spin-orbit and tensor couplings are taken from table 9. S = −1, are experimentally indistinguishable and, hence, mix with each other. Thus, the mass eigenstates, resulting from this mixing, and the weak interaction eigenstates, in which the decay of the Ξ b -baryon would lead dominantly to the states with the internal quantum number (c3[cq]3[sq ]3), are different. In the present analysis, we have not taken into account this mixing.
As for the discovery modes of such pentaquarks in weak decays of bottom baryons, they were studied in [51] in the heavy-quark symmetry limit, assuming SU(3) F -symmetry. Both Cabibbo-allowed b → ccs (with ∆I = 0, ∆S = −1) and Cabibbo-suppressed b → ccd (with ∆I = 1/2, ∆S = 0) transitions were considered. Among the possible modes of the weakly-decaying b-baryons, the most promising modes in searching for the hidden-charm singly-strange pentaquarks at the LHC are: Assuming a similarity in the unflavored and strange pentaquark production, we anticipate the observation of the strange partners of the known pentaquarks. We have marked the masses of these states by enclosing them in solid boxes in table 11. The threshold for the strange pentaquarks in the P 0 Λ → J/ψ Λ decay mode is M thr J/ψ Λ = m J/ψ + m Λ = 4212.58 MeV [3]. With the masses given in table 11, there is only one state with the "good" light diquark (S ld = 0), the J P = 1/2 − with the mass M = (4112 ± 32) MeV, which lies 100 MeV below the M thr J/ψ Λ threshold. This will result into a state narrower than one should expect due to a strong decay, with the cc-pair annihilating into light hadrons or a lepton pair. This makes a good case for a search of narrow Λ µ + µ − (and Λ e + e − ) resonant structure. The thresholds in other two modes P +,0 Σ → J/ψ Σ +,0 are higher:  these states contain the "bad" light diquark and their production in the Ξ b -baryon decays is suppressed by heavy quark symmetry.
The masses of doubly-strange pentaquarks are presented in tables 12 and 13 and compared with the results of [51,52]. The difference between the entries in these tables is as follows: the states from table 12 contain both the strange heavy diquark and the strange light diquark, while those from table 13 have hidden-charm unflavored triquark and doublystrange light diquark, i.e., the pentaquark structures are (c3[cs]3[sq]3) and (c3[cq]3{ss}3), respectively, with q being u-or d-quark. The number of states with the doubly-strange light diquark is smaller 5 because such a diquark, being flavor-symmetric, can be only the 5 In refs. [51,52] similar analysis was done for the hidden-charm doubly-strange pentaquarks and mass predictions for the states with the doubly-strange "good" light diquark were presented erroneously. The corresponding entries should be excluded from tables V and VI (the last values in the columns entitled PX 2 and PY 2 , respectively) in [51] and from tables V (the last values in the columns entitled PX 1 and PX 3 ), VI (similar values in the columns PY 1 and PY 3 ), and VII (the last value in the columns entitled PY 7 ) in [52].  "bad" diquark with the spin S ld = 1. So, this means that the pentaquarks with the "good" light diquark from table 12 remain unmixed while the states with the "bad" light diquark mix, in general, with the corresponding states (with the same set of quantum numbers: spin-parity J P and isospin I) from table 13. In this paper, we neglect this mixing. Based on the analysis in [51], we specify the most promising decay mode of the Ω bbaryon for searching of hidden-charm doubly-strange pentaquarks at the LHC as Ω − b → P 0 Ξ 10 K − → J/ψ Ξ 0 K − . The threshold for the doubly-strange pentaquarks in the P 0 Ξ 10 → J/ψ Ξ 0 decay mode is M thr J/ψ Ξ 0 = m J/ψ +m Ξ 0 = (4628.70±0.32) MeV [3]. With the masses given in table 13, there are two mass degenerate states, having J P = 1/2 − and J P = 3/2 − , with the mass M = 4511 MeV, which lie below the M thr J/ψ Ξ 0 threshold. Also, the third orbitally-excited state with the spin-parity J P = 1/2 + and the mass M = 4629 MeV could also lie below this threshold. All the other states presented in table 13 can be produced in the Ω b -baryon decays and their observation will herald the era of the as yet unknown field of doubly-strange pentaquarks.

JHEP10(2019)256
The masses of the triple-strange pentaquarks are presented in table 14. As these pentaquarks do not contain u-or d-quarks, they are isospin singlets. All these states have the doubly-strange "bad" light diquark only, i.e., the pentaquark structure is (c3[cs]3{ss}3). They can be produced at the LHC in the weak decays of the Ω b -baryons, in particular,  due to the b → ccs transition which gives rise to the quasi-two-body decay mode Ω − b → P − Ω 10 φ → J/ψ Ω − φ. The threshold for the triple-strange pentaquarks in the P − Ω 10 → J/ψ Ω − decay is M thr J/ψ Ω − = m J/ψ + m Ω − = (4769.35 ± 0.29) MeV [3]. With the masses given in table 14, there are two mass degenerate states, having J P = 1/2 − and J P = 3/2 − , with the mass M = (4642 ± 31) MeV, which lie ∼ 100 MeV below the M thr J/ψ Ω − threshold. Also, the third orbitally-excited state with the spin-parity J P = 1/2 + and the mass M = (4761 ± 37) MeV could also lie below this threshold. All the other states presented in table 14 can be produced in the Ω b -baryon decays. The other possibility is to produce such pentaquarks promptly, but their production is highly suppressed due to their triple-strange content.

Pentaquark decay widths
In the compact diquark picture, the quarks in a diquark are bound and not free. In the present context it means that there is a barrier (or bound-state effect) which reduces the probability of thec-quark and the charm quark in the [uc]-diquark to form a charmonium state. This is seen also in the decays of the X, Y, Z states, which are tetraquark candidates in the compact diquark picture [48,54,[66][67][68][69]. In the case of the hidden-charm tetraquarks, it is the radius of the tetraquark R 4q ∼ 1 − 2 fm (or, perhaps somewhat larger) compared to the size of the compact heavy-light diquark, typically R Qq ∼ 0.5 fm, which acts as a JHEP10(2019)256 barrier to form a J/ψ-or η c -meson. For example, it is argued in [66], that the relative ratio λ is expected to be λ ≡ R 4q /R Qq ≥ 3. At long distance (but within the confined tetraquark radius), the potential is attractive, which changes as the internal structure of the diquark (and anti-diquark) gets resolved. The tunneling probability depends on the mass of the quark, with the probability exponentially suppressed the heavier the quark is. This can be expressed through the semi-classical approximation of the tunneling amplitude where M is the quark mass and the energy E is typically 100 MeV [66].
This suppresses the formation of the charmonium states, leading to a much reduced P c → J/ψ p decay width. In all these cases, the size of the diquark R Qq plays an important role in reducing the decay widths, though the extent of this suppression may be different in the tetraquark and pentaquark cases. To work it out requires a dynamical theory, which is currently not at hand. As some of the states are orbitally excited states, anticipated in the compact diquark picture, the decay widths are further reduced due to the angular momentum barrier. Diquark-size effects are also present in the tetraquark and pentaquark mass estimates in the diquark model, and they have been calculated in a particular case, namely double-bottom tetraquark [70].
For four quarks (q, q , q and c) and the charm antiquark, present in a pentaquark P c , there are different alternatives for clustering in color neutral states. The first possibility is the triquark-diquark alternative: From the color algebra point of view, these states are related, Ψ D 2 = −Ψ D 1 , but other internal dynamical properties can be different. The square parenthesis on the r.h.s. of the second equality in each equation denote color antisymmetrization of the constituents. The color connection of the quarks in Ψ D 1 is assumed in the present paper, while the Ψ D 2 color structure is employed in the dynamical diquark model of multiquark exotic hadrons [50,54,68]. The other color-singlet combinations are the meson-baryon alternatives: Among these four combinations, only Ψ H 1 and Ψ H 2 satisfy the heavy-quark-symmetry condition, which implies that the light [q q ]-diquark is transmitted intact, retaining its spin quantum number, from the b-baryon decays to the pentaquark.

JHEP10(2019)256
It is possible to reduce the states in the color configurations (6.1) or (6.2) to a linear combination of the states in the set (6.3)-(6.6). Keeping the color of the light diquark unchanged, the convolution of two Levi-Civita tensors entering the triquark in (6.1) gives the following decomposition: as mentioned earlier. However, a color reconnection is not enough to reexpress a pentaquark operator as a direct product of the meson and baryon operators. As hadrons are determined by their spin-parities, one should also project spins of the quarks and diquarks onto definite hadronic spin states. This requires to know the Dirac structure of the pentaquark operators to undertake the Fierz transformations in the Dirac space. This is a rather involved problem, already for the S-wave pentaquarks, and requires a dedicated study, to which we hope to return in a future publication. We exemplify it here by considering the three S-wave pentaquarks from table 2 with a light "good" diquark, following similar considerations for the hidden-charm tetraquarks [66]. The heavy diquark in them, being also in S-wave, can have both spins, S hd = 0 and S hd = 1. We further simplify the problem by assuming that the charm quark and antiquark are non-relativistic particles, keeping them as point-like objects at different positions in the Minkowski space-time. The diquark-diquark-antiquark operators of the first state from table 2, containing the heavy-diquark with S hd = 0, after color reconnection (the spin structure is determined by eq. (6.1)) and making the color indices explicit, are as follows: wherec(x), c(y) and q(y) are the charm-antiquark, charm-quark and light-quark spinors in the non-relativistic limit, respectively, i and k are the color indices, the light diquark operator is denoted as d 0 (z), and the usual notation for the product of matrices is employed. With the help of the Fierz transform for the Pauli matrices [15]: we reduce the states (6.8) and (6.9) to the forms: 12) where the quark-antiquark and quark-diquark combinations are properly normalized in the color space.
The transition to hadrons needs specifying the flavors of all the quarks. If, as an example, we consider pentaquarks with the same quark content as the observed ones, this

JHEP10(2019)256
immediately specifies q = u and d 0 = [u C γ 5 d], being the color-antitriplet scalar state. 6 Neglecting the position dependence of the operators, the quark operators (6.11) and (6.12) can be rewritten in terms of the hadrons: where A and B (a and b) are non-perturbative coefficients associated with different barrier penetration amplitudes for different light-quark (heavy-quark) spin configurations. They are equal in the limit of naive Fierz coupling. As mentioned in subsection 5.2, the predicted mass of this state, M = (3830 ± 34) MeV, is below the J/ψ p and η c p thresholds, so its decay to these states is forbidden by phase space, and it will decay through the charmanticharm pair annihilation, having a width much smaller than the widths of the newly observed pentaquarks in table 1. The decays of this pentaquark to the Λ + cD 0 and Λ + cD * 0 final states are also forbidden by phase space, as these channels have higher thresholds than the proton-charmonium modes. Along similar lines, assigning q = s, but keeping the light diquark d 0 as above, we get the first state with the internal spins S hd = S ld = 0 from table 11. In terms of the physical hadrons, the quark operators (6.11) and (6.12) can be written as follows: where A s and B s (a s and b s ) are non-perturbative coefficients associated with barrier penetration amplitudes. The predicted mass of this state, M = (4112 ± 32) MeV, (see table 11 in subsection 5.4) is below the J/ψ Λ and η c Λ thresholds, so its decay to these states is forbidden by phase space and it will decay through the charm-anticharm pair annihilation. As the thresholds of Λ + cD 0 s and Λ + cD * 0 s are higher, the decays of this pentaquark to these final states are also forbidden by phase space.
In the same approximation, the diquark-diquark-antiquark operators of the other states from table 2, containing the heavy-diquark with S hd = 1, after color reconnection (the spin structure is determined by eq. (6.1)) and making the color indices explicit, are as follows: These operators, being the direct product of a spinor and a vector, are required to be decomposed into two irreducible representations, corresponding to the spin J = 1/2 and J = 3/2 pentaquarks. The projection onto these states is done after Fierz transform is

JHEP10(2019)256
performed, for which one uses the relation [15]: With this, the states (6.17) and (6.18) are reduced to: where the quark-antiquark and quark-diquark combinations are properly normalized in the color space. Neglecting as above the position dependence, the quark operators (6.20) and (6.21) of the pentaquarks with the light quark q = u and light diquark d 0 = [u C γ 5 d] can be rewritten in terms of the hadronic states: where A , B and C (a , b and c ) are non-perturbative coefficients associated with different barrier penetration amplitudes. They are equal in the limit of naive Fierz coupling. These operators are the linear combination of two operators corresponding to baryons with spins J = 1/2 and J = 3/2. They can be separated by using the equation: with the following result: This analysis shows that the unflavored pentaquark with the spin-parity J P = 3/2 − is mainly decaying either to J/ψ p final state, in which it was observed due to our assignment, or to Λ + cD * 0 state.

JHEP10(2019)256
If instead we assign q = s, but keeping the light diquark d 0 as above, we get two states with spins S hd = 1 from table 11 and the corresponding operators are as follows: It is easy to recognize that the main decays in which a singly-strange J P = 3/2 − pentaquark should be searched, are the J/ψ Λ and Λ + cD * 0 s final states. Analysis of the other pentaquark states can be done along similar lines.
From the foregoings in this section, it is evident that the color representations and color reconnections, on one hand, and localization of the diquarks and quarks, on the other, are essential to disentangle the dynamics (production and decays) of multiquark hadrons. We have worked out some illustrative cases here, in the color-antitriplet diquark representation, imposing heavy-quark symmetry constraints in the decays of b-baryons.
We end this section with a remark about the composition, localisation, and bindings of the various building blocks of a pentaquark in our model. In the compact diquark picture which we follow, the quarks in a diquark are bound and not free, and they have a hadronic size, which could be as large as several tenths of a Fermi. This applies to both the light [q q ]-diquark as well as the heavy-light [qc]-diquark. In this respect, the compact pentaquark model differs from the hadrocharmonium model [55], in which the light degrees of freedom are considered as a cloud. We implicitly assume that the diquark sizes are smaller than that of the triquark, which in turn is smaller than the pentaquark size. So, the compact pentaquark model is essentially a three-scale problem. The light diquark is localized at some distance from the doubly-heavy triquark and an overlap of their wave functions can be described within the one-gluon-exchange approach. The pentaquark stabilization is essentially determined by the triquark stabilization, which is accomplished by assuming that there is a barrier/bound-state-effect which reduces the probability of thē c and the charmed quark in the heavy [qc]-diquark to annihilate each other. This barrier is essential, as otherwise the c-andc-quarks will immediately annihilate. The triquark topology is not too different than that of a tetraquark consisting of [qc]-and [q c]-diquarks, where the tetraquark is stabilized due to a barrier between the diquark and the anti-diquark shells [48,49].
As a final remark, we note that in the diquark models light diquarks [qq ] are ubiquitous. They are present in baryons, such as Λ b , and the hidden-charm pentaquarks, discussed in detail in our paper, but also in the excited baryons, such as Ω * c = c{ss} [59], and in doubly-heavy tetraquarks, such asbb[qq ], discussed in a number of papers on stable tetraquarks [41][42][43][44][45]. How their hadronic sizes impact on properties of a given multiquark state depends on an underlying dynamical model. In the constituent diquark model, which

JHEP10(2019)256
we use, the size of the light diquark influences various parameters of the effective Hamiltonian. We have ignored the differences which emerge from intrinsically different embeddings of the diquarks in the baryons, tetraquarks and pentaquarks. This is an approximation, as the hadrons in question are not expected to have a universal hadronic size, and it limits the precision of the pentaquark masses presented here. A detailed dynamical scheme for pentaquarks remains to be developed -an issue which we hope to address in the future.

Conclusions
We have presented the mass spectrum of the hidden-charm pentaquark states (ccqq q ), where q, q , and q are light u-, d,-and s-quarks, using isospin symmetry for the u-and d-quarks. In doing this, we have used an effective Hamiltonian, based on a doubly-heavy triquark -light diquark picture, shown in figure 1. Apart from the constituent diquark and quark masses, the Hamiltonian incorporates dominant spin-spin, spin-orbit, orbital and tensor interactions. As the first step, we present analytical expressions for the masses of the states having well-defined quantum numbers, making explicit the transformations needed in the relations among various angular momentum basis vectors. This formalism is the main theoretical contribution of our paper.
For the numerical estimates of the pentaquark masses, we use Model II by Maiani et al. [56], developed to calculate the X, Y , and Z tetraquark mass spectrum in the compact diquark model, in which only the dominant spin-spin couplings are taken into account. Since not all parameters in the effective Hamiltonian can be uniquely determined, as there are currently only three pentaquark states observed so far, and also their spinparity assignments are not known, we use the results of the known heavy tetraquarks and heavy baryons to fix them. Due to the paucity of data on pentaquarks, and the inference on the model parameters from other systems, it is difficult to be quantitatively precise. A remark about the accuracy on our mass estimates due to the neglect of the diquark size effects is in order here. While we have neglected this, it has been estimated recently for the double-bottom hadrons (baryons and tetraquarks) by taking into account gluon exchanges between the spectator quarks and the two heavy quarks in the diquark. Using the non-relativistic potential model to estimate the matrix element of the resulting dominant operator leads to a contribution ∆M Ξ bbq 30 MeV, with a similar contribution to the mass shift of the corresponding tetraquark [70]. The corresponding effect for the hidden-charm tetraquarks and pentaquarks is not yet at hand, but is presumably of a similar order of magnitude. Thus, there is an intrinsic diquark-size dependent uncertainty in our mass estimates, which we hope can be estimated in an improved dynamical diquark framework. In view of this, the masses of the pentaquarks in the compact diquark approach presented here, while being in the right ball-park of the current experimental measurements, require an improved theoretical framework and more data to be precise.
The resulting pentaquark spectrum is very rich, as the states have both spin-0 and spin-1 light diquarks. Since the current experimental effort is dominated by searches in Λ bdecays, we impose heavy quark symmetry, in which the light-diquark spin in Λ b -baryon is a good quantum number. This, in turn, implies that only pentaquarks with the spin-0 light JHEP10(2019)256 diquarks are anticipated in Λ b -decays. This reduces the spectrum of the observable states greatly, yet there are more states present in the spectrum than have been observed, and hence we are not able to assign the J P quantum numbers to the observed states uniquely. We have argued, based on theoretical consistency, why some spin-parity assignments for the new pentaquarks can be ruled out in the compact diquark model. Yet, by virtue of the fact that both S-and P -state pentaquarks are anticipated in such models, predicting both positive-and negative-parity states, the spectrum provides a generic discrimination between the compact pentaquarks and the ones foreseen in hadron molecules, in which the observed pentaquarks have all the negative parity. Thus, it is crucial to determine the J P quantum numbers of the observed pentaquarks. Likewise, the remaining states in the pentaquark spectrum, whose masses we have presented, should be searched for. We find that some of them may lie in mass below their strong thresholds, defined by J/ψ p (or η c p). Hence, they will decay via the cc-pair annihilation into light hadrons and into a lepton pair. Thus, the mode p µ + µ − may reveal the presence of these additional pentaquark states, and we urge the LHCb collaboration to search for them.
The hidden-charm pentaquarks having a strange quark in their valence composition are at present an experimentally virgin territory. While they can be produced also promptly, like their non-strange counterparts, only search strategies based on weak decays of the b-baryons with a strange quark are expected to bear fruits. The strange b-baryons which decay weakly are Ξ 0 b (bus), Ξ − b (bds), and Ω − b (bss). Of these, the decays of Ω − b (bss) will produce pentaquarks with a spin-1 light diquark, in accordance with the heavy quark symmetry, thus they allow to investigate the SU(3) F -decuplet pentaquark sector. However, the production of such b-baryons is suppressed due to the excitation of the ss-pair from the vacuum (for producing Ξ 0 b and Ξ − b in a b-quark initiated jet) and twice suppressed (for the Ω − b -baryon). As an example, the product branching ratio B( [3]. This, compared with the product branching ratio B( , is about 20 times less probable. Thus, the discovery of hidden-charm strange pentaquarks requires a lot more data, but they remain a firm prediction, as in the compact diquark model complete SU(3) F -multiplets are anticipated. Hence, their non-observation in experiments would seriously undermine such models. In the foreseeable future, experiments at the LHC, in particular LHCb, remain our best hope for their discovery. In the long-term, the entire hidden-charm pentaquark sector can be studied at a Tera-Z factory. A Corrections due to spin-spin interactions between the constituents of doubly-heavy triquark and light diquark The matrix elements of the operators in eq. (3.4) have in general six different coupling strengths. For the pentaquarks, which contain the light u-and d-quarks only, isospin symmetry among the light quarks can be used to reduce the number of couplings: (K cq )3 = (K cq )3 = (K cq )3, (K qq )3 = (K qq )3 = (K qq )3, andKc q =Kc q =Kc q . So, the expression for H t−ld SS in the effective Hamiltonian responsible for the spin-spin interactions between the (anti)quarks in the diquark and triquark, given in eq. (3.4), is greatly simplified: where the definition of the light-diquark spin operator, S ld = S q + S q , is used. We start by calculating the matrix elements of the spin-spin operators from the first term, 2Kc q (Sc · S ld ), which contains the spins of the charm antiquark and light diquark. Again, it is convenient to transform the vector basis. To do so, we denote the pentaquark state The required matrix elements can be written in the form: where the charm-antiquark spin, Sc = 1/2, is substituted. We calculate now the matrix elements of the last two operators in eq. (A.1). Let us start from the operator 2 (K cq )3 (S c · S ld ) which works on the wave functions of the charm quark in the heavy diquark and light diquark. To find it, we need to couple the heavy and light diquark first and then couple the light diquark with the charm quark -the constituent of heavy-light diquark. This requires a double transformation of the original pentaquark vector state |S hd , S t , L t ; S ld , L ld ; S, L J , which we denote as |S ld , (S hd , Sc) St ; S , to the vector JHEP10(2019)256 state |(S ld , S hd ) S 1 , Sc; S ; in the second transformation these basis states, for simplicity called as |S ld , (S c , S q ) S hd ; S 1 , are changed to the required ones, |(S ld , S c ) S 2 , S q ; S 1 . Specifically: The matrix elements of interest are as follows: where the quark and antiquark spins, S c = S q = Sc = 1/2, are substituted. Note that the matrix elements of the last operator in (A.1) 2 (K qq )3 (S q · S ld ) differ from the above matrix elements by the replacement (K cq )3 → (K qq )3 and the interchange S c = 1/2 ↔ S q = 1/2. It means that, up to the coupling strengths, they coincide. So, the contribution of these two operators enters the matrix elements with the factor (K cq )3 + (K qq )3. We continue by considering the spin-spin interaction between the strange hidden-charm triquark and unflavored (i.e., having u-and d-quarks) light diquark. In this case, the expression given in eq. (3.4) of the effective Hamiltonian is again simplified in the isospinsymmetry limit: where S s and S ld = S q + S q are the spin operators of the s-quark in the triquark and the light-diquark, respectively. Comparing (A.6) with (A.1), it is easy to recognize that they differ by the factor entering the last terms, (K sq )3 and (K qq )3. So, the matrix elements (A.3) and (A.5) remain unchanged, but to get the contribution of the last operator in (A.6) one should simply change the coupling: (K qq )3 → (K sq )3. Taking into account the statement in the paragraph above, the second and third terms in (A.6) are equal (modulo couplings) and contribute with the coefficient (K cq )3 + (K sq )3 in the mass formulae.

JHEP10(2019)256
The same simple form of the spin-spin interaction is also hold for the hidden-charm pentaquark with the strangeness S = −2 which consists of the unflavored doubly-heavy triquark and light diquark containing two s-quarks: So, the matrix elements (A.3) and (A.5), up to spin-spin couplings, remain unchanged and their contributions to the mass formulae are accompanied by the couplingsKc s and (K cs )3 + (K sq )3, respectively. Similar simple effective Hamiltonian for the spin-spin interactions is present in the hidden-charm pentaquarks with the strangeness S = −3, which contains three s-quarks: The matrix elements (A.3) and (A.5) can be safely used in this case also and their contributions are entering the mass formulae with the couplingsKc s and (K cs )3 + (K ss )3, respectively. The next flavor state which corresponds to the hidden-charm strange pentaquark consists of the the unflavored doubly-heavy triquark and strange light diquark. The corresponding spin-spin effective Hamiltonian (3.4) does not simplified and all the six terms, each of which is accompanied by its own coupling, should be considered separately: where we identify q = s and apply the isospin symmetry to the spin-spin couplings. A calculation of these operators requires a recoupling of four angular momenta. In general, this transformation scheme is expressed through the 9j-symbol which is defined as follows [60]: (j 1 , j 2 ) j 12 , (j 3 , j 4 ) j 34 , j|(j 1 , j 3 ) j 13 , (j 2 , j 4 ) j 24 , j = (2j 12 + 1) (2j 34 + 1) (2j 13 + 1) (2j 24 + 1) (A.10) To demonstrate its necessity, let us consider the first term in (A.9). In our basis of states the charm antiquark belongs to the doubly-heavy triquark while the light q -quark is inside the light diquark. The original pentaquark vector state |S hd , S t , L t ; S ld , L ld ; S, L J , which we denote for simplicity as |(Sc, S hd ) St , (S q , S s ) S ld ; S , should be transformed to the new state with the recoupling momenta |(Sc, S q ) Sc q , (S hd , S s ) S hd,s ; S , and the 9j-symbol is responsible for this transform:

JHEP10(2019)256
So, the matrix element of the operator 2Kc q (Sc · S q ) sandwiched between the states with different triquark S ( ) t and light-diquark S ( ) ld vector states can be written as a double sum (here, we come back to the original notation for the vector state): J S hd , S t , L t ; S ld , L ld ; S, L|2Kc q (Sc · S q )|S hd , S t , L t ; S ld , L ld ; S, L J =Kc q (2S t + 1) (2S t + 1) (2S hd + 1) (2S hd + 1) Sc q ,S hd,s (2Sc q + 1) (2S hd,s + 1) where the spins Sc = S q = S s = 1/2 of quarks were subsituted. In a similar manner one can get the contribution of the operator 2Kc s (Sc · S s ) with the result: Up to the coupling,Kc s , this equation coincides with (A.12), so their common contribution contains the factorKc q +Kc s .
With the other four operators the situation is a little bit more complicated as spin-spin interactions relate quarks from different diquarks, moreover, the spin of the heavy diquark is coupled with the spin of the charm antiquark inside the triquark system. So, at the first step one needs to decouple the heavy diquark from the charm antiquark and couple it with the light diquark which generates the tetraquark system. If we denote the vector state as |S ld , (S hd , Sc) St ; S , we need to recouple in a way to get the new basis |(S ld , S hd , ) S T Sc; S which is possible to do with the help of the transformation (4.4): At the second step, it is necessary to recouple quark spins inside the tetraquark. So, for this transformation the 9j-symbol (A.10) is required. To do this, the introduced pentaquark vector state |(S ld , S hd ) S T , Sc; S is more convenient to rewrite in the form

JHEP10(2019)256
|(S c , S q ) S hd , (S q , S s ) S ld ; S T and then perform the transformations: Now, we can sandwich the operators 2 (K cq )3 (S c · S q ) and 2 (K sq )3 (S s · S q ) by the vector states with varying spins of the heavy, S hd , and light, S ld , diquarks as well as of the doublyheavy triquark, S t . The resulting matrix element of these two operators is the triple sum: The matrix elements of the rest two operators 2 (K cs )3 (S c · S s ) and 2 (K qq )3 (S q · S q ) are calculated similarly with the following result: J S hd , S t , L t ; S ld , L ld ; S, L|2 (K cs )3 (S c · S s ) + 2 (K qq )3 (S q · S q )|S hd , S t , L t ; S ld , L ld ; S, L J = (K cs )3 + (K qq )3 (2S t + 1) (2S t + 1) (2S hd + 1) (2S hd + 1) (2S ld + 1) (2S ld + 1) S T ,Scs,S qq (2S T + 1) (2S cs + 1) (2S qq + 1) S cs (S cs + 1) − 3 2 In general, eqs. (A.16) and (A.17) are the matrices which mix all the states having the same total spin, S, and the same parity. The last flavor state which corresponds to the hidden-charm strange pentaquark with the strangeness S = −2 consists of the the strange doubly-heavy triquark and strange light diquark. The corresponding spin-spin effective Hamiltonian (3.4) does not get simplified in this case also, and all the six terms, each of which is accompanied by its own coupling, should be considered separately: +2Kc s (Sc · S s ) + 2 (K cs )3 (S c · S s ) + 2 (K ss )3 (S s · S s ),

JHEP10(2019)256
where we make the substitution q = s, q = q, and q = s . The above consideration holds also here, and one can use the expressions (A.13), (A. 16), and (A.17) with the obvious replacements (K sq )3 → (K ss )3 in (A.16) and (K qq )3 → (K sq )3 in (A.17). The contributions discussed in this appendix are relevant for the Model I by Maiani et al. [35], where such (subdominant) spin-spin interactions are included in the effective Hamiltonian. We presented the formalism here. They are neglected in the numerical analysis given in this paper, as we use the Model II [56] for our mass estimates. B Mass derivations for the P -wave pentaquark states B.1 The states with J P = 3/2 + and triquark spin S t = 3/2 The mass matrix of the states with the spin-parity J P = 3/2 + and triquark spin S t = 1/2 is as follows: Before calculating the mass spectrum, we make the three matrices above to be traceless: There is a standard procedure to find the eigenvalues, λ i , i = 1, 2, 3, of the traceless symmetric (3 × 3) matrix. Denoting such a matrix as its characteristic equation

JHEP10(2019)256
is expressed in terms of three rotational invariants: (B.6) Since real eigenvalues are of our interest, two non-trivial invariants, I 2 and I 3 , must satisfy the following conditions: If these conditions are fulfilled, we can define the angle: through which the sorted eigenvalues (λ 1 > λ 2 > λ 3 ) can be written as follows: After these eigenvalues are obtained, the masses are calculated easily: To derive the expressions for the masses of the four states with the spin-parity J P = 3/2 + and triquark spin S t = 1/2 from table 5, one needs to find the eigenvalues of the following non-diagonal symmetric (4 × 4) matrix: Noting that the last three matrices are block-diagonal and their non-trivial (2 × 2) blocks are identical. A diagonalization of these matrices is reduced to the diagonalization of the blocks: where the following parameter is introduced: (B.14) Coming back to the (4 × 4) mass matrix (B.12), it is necessary to point out that the transformation to the diagonal form of the sum of the last three matrices in (B.12) can be done by using the orthogonal matrix: For the first four states with the spin-parity J P = 1/2 + and triquark spin S t = 1/2 from table 5, the mass matrix is again the non-diagonal symmetric (4 × 4) matrix: The rotational matrix (B.15) with some fixed value of the angle α transforms the sum of the last three matrices to the diagonal form. As in the case of J P = 3/2 + pentaquarks, the diagonalization procedure can be applied to the (2 × 2) blocks: where the following parameter is introduced: µ (1/2) = 1 30 25 (7A t + 2A ld + 3b) 2 + 2 (20A t + 20A ld + 21b) 2 .

(B.23)
After this transformation is done, the (4 × 4) mass matrix (B.21) takes the form:  Table 15. Experimental data from the LHCb [63] and Belle [64] collaborations. (Masses are given in MeV.) All the measurements have the uncertainty +0.3 −0.5 MeV coming from the mass of the groundstate Ξ + c -baryon, which is common in both analysis. The spin-parity, J P , assignment assumes that all the states are orbitally-excited baryons as suggested in [59]. orthogonal transformations and we get the following set of masses: C χ 2 -analysis of orbitally-excited Ω * c -baryons Existing experimental data on the masses of orbitally-excited Ω * c -baryons were obtained by the LHCb [63] and Belle [64] collaborations and collected in table 15.
The theoretical expressions for the orbitally-exciting Ω * c -baryon masses follow from the effective Hamiltonian [59]: 7 H eff = m c + m {ss} + 2κ ss (S s 1 · S s 2 ) + B Q 2 L 2 + V SD , (C.1) In eq. (C.1), m c and m {ss} are the masses of the c-quark and the spin-1 {ss}-diquark, respectively, κ ss is the spin-spin coupling of the quarks in the diquark, L is the orbital JHEP10(2019)256 angular momentum of the Ω * c -baryon, and B Q is the orbital coupling. The coefficients a 1 and a 2 are the strengths of the spin-orbit terms involving the spin of the diquark S {ss} and the charm-quark spin S c , respectively, c is the strength of the spin-spin interaction between the diquark and the charm quark, and S 12 /4 represents the tensor interaction, defined by where n = r/r is the unit vector in the direction from charm quark to the doubly-strange diquark.
One should remember that several required internal parameters are fixed: charmquark spin S c = 1/2, diquark spin S {ss} = 1, and Ω * c -baryon orbital angular momentum L = 1. So, in the L − S coupling scheme, two possible values of the total spin S = 1/2 and S = 3/2, after their coupling to L, allow to get five P -wave states: two with the total angular momentum J = 1/2, two with J = 3/2 and the last one with J = 5/2 (see table 15). In this scheme the spin-dependent part, V SD , is represented by the block-diagonal matrix, each block of which corresponds to the states with a fixed value of J [59]: After diagonalizing the matrices and adding the common mass term   Table 16. The χ 2 -analysis of the orbitally-excited Ω * c -baryon masses based on the measurements by the LHCb and Belle collaborations. This analysis is performed based on the LHCb data alone and on the combined data set from both the LHCb and Belle collaborations presented in table 15. , a 1 , a 2 , b, c} in these five equations. If there are five experimentally measured states which can be assigned with the P -wave Ω c -baryons, all five variables, in general, can be determined but their values do not always satisfy physical requirements.
After the assignment is done as specified in table 15, it is possible to perform a χ 2analysis of the data. We start from the LHCb data alone as all five states are present there and then add the Belle measurements which allow to find χ 2 min for four degrees of freedom. The best-fit values and 1σ uncertainties of free parameters are presented in table 16. There are no free degrees of freedom for the fit based on the LHCb data alone (the first row in table 16) but combining them with the Belle data results χ 2 min /ndf = 0.87/4. The correlation matrices in both cases are as follows: