Spectrum of the $uudc \bar c$ hidden charm pentaquark with an SU(4) flavor-spin hyperfine interaction

We study a few of the lowest states of the pentaquark $uudc\overline{c}$, of positive and negative parity, in a constituent quark model with an SU(4) flavor-spin hyperfine interaction. For positive parity we introduce space wave functions of appropriate permutation symmetry with one unit of orbital angular momentum in the internal motion of the four-quark subsystem or an orbital excitation between the antiquark and the four quark subsystem which remains in the ground state. We show that the lowest positive parity states $1/2^+, 3/2^+$ are provided by the first alternative and are located below the $1/2^-$ and the $1/2^+$ states with all quarks in the ground state. We compare our results with the LHCb three narrow pentaquark structures reported in 2019.

Recently the LHCb Collaboration have updated their analysis of the Λ 0 b → J/ψK − p decay [5,6] and observed three narrow peaks with masses and widths as given in Table 1, where the entirely new P + c (4312) has a 7.3σ statistical significance and the other two resonances replace the previous P + c (4450) now resolved at 5.4σ significance. The broad P + c (4380) resonance observed in 2015 awaits a e-mail: fstancu@ulg.ac.be confirmation and a proper identification of the spin-parity is necessary in all cases. The mass of P + c (4312) lies a few MeV below the Σ + c D 0 threshold (4318 MeV) and the masses of P + c (4440) and P + c (4457) lie a few MeV below the Σ + c D * 0 threshold (4460 MeV). Proximity of the Σ + c D 0 and Σ + c D * 0 thresholds to these narrow peaks made their interpretation as molecular S-wave of Σ + c D 0 or Σ + c D * 0 quite natural [7,8,9,10,11,12]. In such an interpretation, they all acquire a negative parity. In Ref. [13] it is shown that the usual molecular scenarios cannot explain the production rate of P c states in Λ b decays. Moreover, the paper points out the need to couple the Σ + c D * 0 and the Λ c (2595)D channels due to the very close proximity of their thresholds. As a result, the second and the third resonances of Table 1 acquire opposite parities, namely J P (4440) = 3 2 − and J P (4457) = 1 2 + and respectively.
The 2019 LHCb pentaquarks have also been analysed in the compact diquark model [14]. This analysis has been extended in Ref. [15], and presently includes a spin-spin, a spin-orbit and a tensor interaction, obtained from the integration in the color space of the one gluon exchange (OGE) interaction between quarks [16]. The basis contains only the color3 configuration, therefore is smaller than the basis used in the present paper. Also there is a belief that the spin 0 diquarks are more tightly bound than the spin 1 diquarks. The lowest positive and negative parity states were calculated and it was found that the lowest state has negative parity.
Here we explore the spectrum of pentaquarks within a quark model [17,18,19], which has a flavor dependent hyperfine interaction. Contrary to the OGE model result, the lightest pentaquark has positive parity, as shown below. As a general feature, the chromomagnetic (color-spin) and the flavor-spin interactions predict contradictory results regarding the parity of open charm exotic systems, as the tetraquark uucc, the pentaquarks uuddc, uudsc and the hexaquark uuddsc [20].
In Refs. [17,19] the hyperfine splitting in hadrons is obtained from the short-range part of the Goldstone boson exchange interaction between quarks. The main merit of the model is that it reproduces the correct ordering of positive and negative parity states in the baryons spectrum in contrast to the OGE model. On the other hand it does not apply to the hyperfine splitting in mesons because it does not explicitly contain a quark-antiquark interaction.
Here we extend the model of Refs. [17,19] from SU(3) to SU(4) in order to incorporate the charm quark. The extension is made in the spirit of the phenomenological approach of Ref. [21] where, in addition to Goldstone bosons of the hidden approximate chiral symmetry of QCD, the flavor exchange interaction is augmented by additional exchange of D mesons between u, d and c quarks and of D s mesons between s and c quarks. The model provided a satisfactory description of the heavy flavor baryons known at that time. With the SU(4) extension we give a more fundamental ground to the work of Ref. [17] and test the corresponding approximations. In the following we shall use the name of flavor-spin (FS) model for the SU(4) extension. The extended flavor exchange allows a complete antisymmetrization of quarks even for quarks of widely different mass, since the Young diagrams (or partitions) serve to label the irreducible representations both of the permutation and of the SU(4) group [22].
Presently we study the pentaquarks of structure uudcc in order to see whether or not the FS model can accommodate the new resonances observed at LHCb [6]. The approach is similar to that previously used in Ref. [23] for the positive parity uuddc and uudsc pentaquarks. That work preceded the H1 experiment [24] which has found evidence, so far not confirmed by other experiments, for a narrow baryon resonance in the D * ± p invariant mass, interpreted as a uuddc pentaquark. The mass prediction made in Ref. [23] is not far from the observation of the H1 experiment.
The parity of the pentaquark is given by P = (−) ℓ + 1 . For the lowest positive parity states (see below) the subsystem of four quarks is defined to carry an angular momentum ℓ = 1 in their internal motion. This implies that this subsystem must be in a state of orbital symmetry [31] O . Although the Table 2 The simplified interaction -V χ /C χ , Eq. (1), integrated in the flavor-spin (FS ) space for four quark subsystems in three different states defined in column 2 in the FS coupling. The notation is: O for orbital, C for color, F for flavor, S for spin.

State
Definition  (1), suggests that the [31] O symmetry leads to a stable pentaquark of positive parity and the lowest state of symmetry [4] O is unbound.
In the exact SU(4) limit, like in the exact SU(3) limit, the flavor-spin interaction introduced in the next section takes the following form with λ F i the Gell-Mann matrices, σ i the Pauli matrices and C χ an equal strength constant for all pairs.
Let us consider two totally antisymmetric states with the orbital symmetry [31] O , allowed by Pauli principle. They are written in the flavor-spin (FS) coupling scheme, as given in Table 2. They are consistent with the isospin and J P combinations of Table 2 of Ref. [25].
The expectation value of (1), as shown in Table 2 (see the next section for the derivation of the last column), is −27 C χ for | 1 and −21 C χ for | 2 .
First we consider the lowest state, i.e. | 1 . The quarkantiquark interaction is neglected in the description of mesons as well, so that the meson Hamiltonian contains a kinetic and a confinement term only. If unstable against strong decays, the pentaquark would spilt into a baryon q 3 and a meson QQ, where Q stands for a heavy flavor. Then, in discussing the stability we introduce the quantity where E defines the energy (mass) of a system with the content defined in the bracket. In our schematic estimate, we suppose that the confinement energy roughly cancels out in ∆E. Then, the kinetic energy contribution to ∆E is ∆KE = 5/4 ω in a harmonic oscillator model [26]. Table 2 shows that the FS contribution to the state | 1 is −27 C χ . For the nucleon or Λ c it is -27/2 C χ (see next section). Then we have ∆V χ = −27/2 C χ . With ω ≈ 157 MeV and C χ ≈ 20 MeV [21] this gives i.e. a considerable binding. This is to be contrasted with the negative parity pentaquarks similar to those containing light quarks, studied in Ref. [27]. For the lowest negative parity state | 3 of uudcc one gets i.e. instability. Thus the FS interaction overcomes the excess of kinetic contribution in [31] O and generates a lower expectation value for [31] O than for [4] O . The paper is organized as follows. In Sec. 2 we introduce the model Hamiltonian and generalize the twobody matrix elements of the FS interaction from SU(3) to SU(4). Sec. 3 describes the orbital part of the four quark subsystem constructed to be translationally invariant both for positive and negative parity states. Sec. 4 summarizes a few analytic details. Sec. 5 exhibits the numerical results for masses and distances between quarks/antiquarks. A qualitative discussion of the decay widths is given in Section 6. A comparison with previous studies of hidden charm pentaquark based on various versions of the FS model is made in Sec. 7. In the last Section we draw some conclusions. Appendix A gives the matrices λ F of SU(4). Appendix B is a reminder of useful group theory formulae for SU(n). Appendix C exhibits the variational solution for the baryon masses relevant for the present study. In Appendix D we construct the explicit form of flavor states of content uudc in the Young Yamanouchi basis for specific irreducible representations [ f ] F .

The Hamiltonian
The nonrelativistic FS Hamiltonian has the general form [17] with m i and p i denoting the quark masses and momenta respectively and r i j the distance between the interacting quarks i and j. The Hamiltonian contains the internal kinetic energy and the linear confining interaction The hyperfine part V χ (r i j ) has a flavor-spin structure which presently is extended to SU(4), similarly to Refs. [28] and [21], except that in the latter, the term V η ′ λ 0 i · λ 0 j of Eq. (7) has been ignored. One has with λ F i (F = 1,2,...,15) defined in Appendix A and λ 0 i = √ 2/3 1, where 1 is the 4 × 4 unit matrix. Note that the norm of λ 0 i is not that satisfying the general property tr(λ a λ b ) = 2 δ ab if 1 is the 4 × 4 unit matrix [23], but we maintain the norm from the SU(3) version of the flavor-spin interaction in order to use the same coupling constant for η and η ′ , which gave a good spectrum for strange baryons. At this stage it is useful to mention the study of Ref. [29] where it has been shown that η ′ has a significant charm component.
In the SU(4) version the interaction (7) contains γ = π, K, η, D, D s , η c and η ′ meson-exchange terms. Every V γ (r i j ) is a sum of two distinct contributions: a Yukawa-type potential containing the mass of the exchanged meson and a short-range contribution of opposite sign, the role of which is crucial in baryon spectroscopy [30]. For a given meson γ the meson exchange potential is In the present calculations we use the parameters of Ref. [30] to which we add the µ D mass and the coupling constant . These are The meson masses correspond to the experimental values from the Particle Data Group [31]. There is no K-and D sexchange in the uudcc pentaquark and, as discussed in the following, we ignore η c -exchange. Therefore the K, D s and η c masses are not needed in these calculations.
For the quark masses we take the values determined variationally in Ref. [28]  After integrating in the flavor space, the two-body matrix elements which generalize Eq. (3.3) of Ref. [17] to SU(4), Table 3 The hyperfine interaction V χ , Eq. (7), integrated in the flavor-spin space, for the quark subsystem uudc I =1/2, S = 0. V qaqb γ are defined in Eq. (10) where the upper index q a q b indicates the flavor of the interacting quark pair .
one obtains a result containing new contributions due to charm. This is In Eqs. (10) the pair of quarks i j is either in a symmetric [2] F or in an antisymmetric [11] F flavor state and the isospin I is defined by the quark content. The upper index of V exhibits the flavor of the two quarks interchanging a meson specified by the lower index. In order to keep close to the notations of Ref. [17] the upper index of π and K is not indicated. Obviously, in every sum/difference of Eq. (10) the upper index is the same for all terms. Note that the term 1 2 V uc η c was missing in Ref. [28]. In practice it can be neglected but theoretically it is important because it recovers the exact SU(4) limit exhibited in Table  2.
Using the flavor wave functions given in Appendix D and the expressions (10) one can calculate the matrix elements of the flavor-spin interaction (7) for four quark states. They are presented in Table 3. To recover the SU(4) limit of Table 2 for the uudc configuration, one has to take In the exact SU(4) limit the states [2] F and [11] F of Eq. (10) containing u, d, c quarks become -3/2 C χ and 5/2 C χ respectively. Then, making use of the spin matrix elements σ i · σ j = 1 and σ i · σ j = -3 for S = 1 and S = 0 respectively, one finds that the contribution of the FS interaction to the baryon masses are V χ (N) = V χ (Λ c ) = -27/2 C χ and V χ (∆) = -9/2 C χ , which are close to the SU(3) limits -14 C χ and -4 C χ respectively. As a matter of fact one can also see that the SU(4) limit given in -C χ units, as listed in Table 2, namely 27, 21 and 15 are quite close to the SU(3) limit (three flavors) of 28, 64/3 and 16 respectively, which can easily be obtained using Appendix B.
The exact SU(3) or SU(4) limits are useful in indicating the sequence of positive and negative parity states in the spectrum, which is maintained at a broken unitary symmetry.

Orbital space
This section follows closely the derivation of the orbital part of the four-quark wave function of symmetry [31] O as presented in Ref. [23].
One needs four internal Jacobi coordinates for pentaquarks. The convenient choice is where 1,2,3 and 4 are the quarks and 5 the antiquark so that t gives the distance between the antiquark and the center of mass coordinate of the four-quark subsystem. Then one has to express the orbital wave functions of the four-quark subsystem in terms of the internal coordinates x, y, z for the specific permutation symmetry [31] O . For the lowest positive parity states we have considered an s 3 p structure for [31] O and extended Moshinski's method [32], from three to four particles, to construct translationally invariant states of definite permutation symmetry. In this way we have obtained three independent states denoted below by ψ i , defined as Young-Yamanouchi basis vectors [22] of the irreducible representation [31] O in terms of shell model states r |nℓm where n = 0, ℓ = 1, and, for simplicity, we took m = 0 everywhere. They read One can see that the ℓ = 1 orbital excitation is located in the relative motion between quarks defined by the internal coordinates x, y and z. The states ψ 1 and ψ 2 are the generalization, from three to four quarks, of ψ λ ℓ0 and ψ 3 the generalization, from three to four quarks, of ψ ρ ℓ0 radial wave functions of mixed symmetry [21] used in baryon spectroscopy. Thus, a more appealing notation could be ψ λ 1 ℓ0 , ψ λ 2 ℓ0 and ψ ρ ℓ0 instead of ψ i with i = 1,2,3. In this picture there is no excitation in the relative motion between the cluster of four quarks and the antiquark defined by the coordinate t. Then the pentaquark orbital wave functions ψ 5 i are obtained by multiplying each ψ i from above by the wave function t |000 which describes the relative motion between the four-quark subsystem and the antiquark c. Assuming an exponential behavior we introduce two variational parameters, a for the internal motion of the four-quark subsystem and b for the relative motion between the subsystem qqqc and c. We explicitly have where N = 2 3/2 a 11/4 b 3/4 3 1/2 π 5/2 We also need the orbital wave function of the lowest negative parity state described by the s 4 configuration of symmetry [4] O which is with

Analytic details
The kinetic energy part T of the Hamiltonian (5) can be calculated analytically. For the positive parity states of Table  2 its expectation value becomes and where q = u, d and Q = c. Here, we have m q = 340 MeV and m c = 1350 MeV, as defined by Eq. (9). Taking m u = m d = m Q = m and setting a = b, one can recover the identical particle limit T = 7 2 ω with ω = 2 a 2 /m, see [26]. Note that one can obtain the desired units of all quantities involved, by muliplying/dividing by c ≃ 197.33 MeV fm at the correct power.
By integrating in the color space, the expectation value of the confinement interaction (6) becomes where the coefficients 6 and 4 account for the number of quark-quark and quark-antiquark pairs, respectively, and where i, j = 1,2,3,4,5 (i j). An analytic evaluation gives as well as for r i j with i, j = 1,2,3,4 and and likewise for r j5 where j = 1,2,3.
One also needs to derive the expressions of r 12 and r 45 in the case where the four quarks are in the s 4 configuration. These are and The matrix elements of the spin-flavor operators of Eq. (7) have been calculated using the fractional parentage technique described in Ref. [22] based on Clebsch-Gordan coefficients of the group S 4 [33]. In this way, each matrix element reduces to a linear combination of two-body matrix elements of either symmetric or antisymmetric states for which Eq. (10) can be used to integrate in the spin-flavor space. For positive parity states with one unit of orbital excitation the resulting linear combinations contain orbital two-body matrix elements of type ss V q a q b γ ss , sp V q a q b γ sp and sp V q a q b γ ps where γ = π, D, η, η ′ and q a q b is a pair of quarks from Eq. (10).
It is useful to note that the ratio of the orbital matrix elements of the D meson exchange interaction between the pair uc of quarks, denoted here by V uc D and the matrix elements of the π meson exchange interaction V π , see Eq. (10), turned out to be about 0.25, in agreement with the mass ratio m u,d /m c of Eq. (9).

Results
The four quark states described in the previous sections couple to the antiquark to a total angular momentum J = L + S + s Q , with L, S the angular momentum and spin of the four-quark cluster and s Q the spin of the heavy antiquark.
We have searched for variational solutions for the pentaquark described by the Hamiltonian defined in Sec. 2. The numerical results are presented in Table 4. The wave functions are the product of the four quarks subsystem states of flavor-spin structure defined in Table 3 and the charm antiquark wave function denoted by |c . The orbital wave functions of the pentaquark are given by Eqs. (15)- (20) containing the variational parametres a and b. We have neglected the contribution of V uu η c and of V uc η c , because little uū and dd is expected in η c . We have also neglected V uc η ′ assuming a little cc component in η ′ , which means that we took The optimal values found for the parameters a and b are indicated in each case. The eigenvalues of |1 |c and |2 |c states are degenerate for the allowed values of J in each case. Moreover, for the wave function |2 |c the states with J P = 1 2 + and 3 2 + have multiplicity 2. Presently there are more states than the observed ones. In the case of pentaquarks with heavy flavor, Ref. [14] argues that only part of the spectrum is reachable in Λ b decays.
The SU(4) FS model described in Sec. 2 gives a value for the mass of the state |1 |c of 4273 MeV, close to the mass of the P + c (4312) resonance and supports the assignment J P = As discussed above, in the exact limit one obtains binding for |1 |c , see Eq. (3), but this is not the case for the variational solution in the broken SU(4), which is larger than the lowest thresholds, Nη c (3920) or N J/ψ(4040). It is not surprising, because in the exact limit the matrix elements are equal for all exchanged mesons, which means that the matrix elements of η, η ′ and D exchanges are overestimated, thus introducing more attraction. However the SU(4) exact limit brings useful information about the relative position of various states.
The lowest negative parity state |3 |c is located above the positive parity states, as expected. It has a mass of 4487 MeV close that of the P + c (4440) or the P + c (4457) resonances. It supports the quantum number J P = 1/2 − . It may be useful to look at the compactness of the present pentaquarks by estimating the distances r 12 and r 45 from Eqs. (25)- (29). They are exhibited in Table 4. In the |1 |c state, which is the lowest, one can see that at equilibrium the four quarks are clustered together and the antiquark is somewhat far apart. For the other states the distances r 12 and r 45 are comparable.
The values of the above convenient relative distances were exhibited in order to see whether or not some clustering of four-quark formation is possible. The effect is quite small. One can use alternative relative coordinates convenient either for a molecular-type structure, or for the hadrocharmonium picture. They are both deuteronlike. Unitary transformations can relate various coordinate systems. This can be considered in the future.
Also, one could introduce a third variational parameter, different from a or b, associated to the coordinate z where the charm quark is involved. In that case the analytic work would be more involved. It was found that the energy levels vary slowly with the present variational parameters, so that the minima are very shallow. In view of this observation the spectrum is expected to be a little lowered by the introduction of a new parameter, but the parity sequence will remain the same. However a more sophisticated variational solution would be important for the wave function, thus for studying decay properties numerically. More theoretical work should be done, in wait for more precise data.

Strong decay widths
Assuming a coupling to open channels, which could provide estimates of widths (and level shifts), one can qualitatively discuss the strong decay widths in the context of compact pentaquarks. First, for P-waves the centrifugal barrier will reduce the widths relative to S-waves. Thus, in the present model, the lowest state, having positive parity, will have a smaller width than the negative parity ones, contrary to the arguments of Ref. [25] inspired by OGE or molecular models. Next, the narrowness of the lowest state is related to the overlap of the color-flavor-spin-orbital wave function with an p + J/ψ final state. For the lowest state, [4] FS , the overlap with q 3 + QQ may be small. This problem deserves a separate investigation. Table 4 Lowest positive and negative parity uudcc pentaquarks of quantum numbers I, J P and symmetry structure defined in Table 2. Column 1 gives the state, column 2 the isospin, column 3 the parity and total angular momentum, column 4 the optimal variational parameters associated to the wave functions defined in Sec. 3, column 5 the calculated mass and the last two columns the distances between quarks/antiquarks.

Comparison with other studies
The basis (15) - (17) is entirely consistent with the parity considerations made in the context of the tetrahedral group T d [34], isomorphic with the permutation group S 4 . The states |1 and |2 of Table 2 are states of symmetry 1 − A 1 in the notation of Ref. [34]. Likewise the state |3 of the same table has symmetry 0 + A 1 . There is agreement between our findings and the conclusion of Ref. [34] that the parity of a pentaquark depends on the interplay between the flavorspin interaction and the contribution of the orbital excitation to the mass. Here we found that the flavor-spin attraction compensates the excess of kinetic energy brought in by one unit of orbital excitation in the internal motion of the fourquark subsystem containing a charmed quark, as it was the case for pentaquarks containing only light flavors like uudd [23], so that the lowest state has positive parity in both cases, either for light or light+heavy quarks.
The authors of Ref. [35] have studied the qqqcc pentaquark prior to the 2015 LHCb observation of pentaquarks, in three different models, including the flavorspin (FS) model. The orbital wave function of the four-quark subsystem has symmetry [4] O for both parities. A symmetric state with an s 3 p structure describes the spurious centre of mass motion, while the state of orbital symmetry [31] O derived in Sec. 3 is translationally invariant. Although the radial wave function was not specified, one can infer that the positive parity states of Ref. [35] were obtained by including a unit of orbital angular momentum in the relative motion between the four-quark subsystem and the antiquark. In this case the orbital wave function takes the form where Then the expectation value of the kinetic energy becomes and The expectation value of the confinement interaction is given by Eq. (24) with and The hyperfine interaction integrated in the flavor-spin space is given by the third row of Table 3  The attraction brought by the flavor-spin interaction in the four-quark state [31] FS is not strong enough to compensate the excess of kinetic energy due to the orbital excitation of the antiquark relative to the four quark subsystem. In our case we obtain a mass of 4573 MeV (a = 0.040 GeV 2 , b = 0.037 GeV 2 ) which corresponds to J P = 1/2 + or 3/2 + , higher than the mass 4487 GeV of the 1/2 − state of Table 4. This can explain why the lowest negative parity state 1/2 − in the calculations of Ref. [35] is lower than the lowest 1/2 + state.
In Ref. [36] the spectrum of the uudcc is studied in a chiral quark model where the Hamiltonian contains both the chromomagnetic and the flavor-spin interactions. For the flavor-spin part the symmetry is restricted to SU(3) with the replacement of the s quark by the c quark which, according to Eq. (10), seems to be reasonable. However, the lowest states have negative parity and the results for positive parity states are not shown, claiming that they are unbound. A possible interpretation is that the chromomagnetic dominates over the flavor-spin interaction in those calculations.
However, some support in favor of a flavor-spin interaction in exotics is given by Ref. [37]. Although restricted to tetraquarks described in a diquark basis, it has been shown that the mass difference observed in some charmed tetraquarks can be explained by an isospin dependent interaction, which is an SU(2) flavor-spin interaction. An important remark is that the contribution of the flavor-spin part relative to the one gluon exchange part must be adequate, irrespective of the basis used in the calculations.

Conclusions
We have calculated the lowest part of the mass spectrum of the hidden charm pentaquark uudcc, within the SU(4) version of the flavor-spin model introduced in this work. This is a model which provides an internal structure for pentaquarks, contrary to the models mentioned in the introduction. An important achievement is that for positive parity states we have made a clear distinction between the case where a unit of angular momentum is located in the internal motion of the four-quark subsystem and the previously schematic case used in Ref. [35], where the angular momentum has to be located in the relative motion between the four-quark subsystem and the antiquark in order to have translationally invariant states (no center of mass motion) of specific flavor-spin symmetry. We have presented analytic expressions for the kinetic, confinement and the flavor-spin parts of the Hamiltonian for a simple variational solution. In the numerical estimates we have used model parameters from earlier studies of pentaquarks.
The calculated masses are in the range of the presently observed P + c (4312), P + c (4440) and P + c (4440) resonances of the 2019 LHCb data both for positive and negative parity states.
One should stress that an important feature of the flavorspin model is that it introduces an isospin dependence of the pentaquark states, necessary to discriminate between decay channels, like for tetraquarks [37].
The radial form (8) favors the p states over the s states due to the second term, which makes the FS interaction closer to the SU(4) limit. Therefore it would be interesting to find out how the pentaquark spectrum depends on the particular radial form of the FS hyperfine interaction by choosing the more realistic version of Ref. [38] with a cut-off parameter for each individual meson exchange. Based on symmetry arguments we however expect a similar conclusion for other forms of the hyperfine interaction.
In addition, one should look for more suitable coupling constants related to the presence of the charmed quark.
Presently, the most important conclusion is that the FS model predicts positive parity for the lowest state in contrast to the widely used OGE model, with or without correlated quarks/antiquarks. Among the recent studies interpreting the 2019 LHCb data, positive parity states have been considered in Ref. [14] within a doubly-heavy triquark -light diquark model. In that model, besides the spin-spin chromomagnetic interaction, the Hamiltonian incorporates spin-orbit and orbital interactions. The fitted parametres lead to a positive contribution of the angular momentum dependent terms which make the ℓ = 1 states to be higher than the ℓ = 0 states. This means that the lowest state has a negative parity, contrary to our results.
The 2019 LHCb pentaquarks have also been studied as hadrocharmonium states [39]. In that model the new P + c (4312) resonance is interpreted as a bound state of χ 0 (1P) and a nucleon with I = 1/2 and J P = 1/2 + . The other two resonances P + c (4440) and P + c (4457), have a negative parity. The parity sequence in the spectrum is closer to the one produced by the flavor-spin model.
The main issue of the present work is the parity sequence in the pentaquark uudcc spectrum. Thus, a crucial test is to experimentally determine the spin and parity of the 2019 LHCb resonances in order to discriminate between different interpretations, in particular between the OGE and the flavor-spin quark models, the doubly-heavy triquarklight diquark model, the molecular picture or the hadronic molecular picture.

Appendix A: The SU(4) generators
Here we reproduce the λ F matrices which are the SU(4) generators in the fundamental representation of SU(4) [22]. Implementing them in Eq. (7) one can obtain the two-body matrix elements of Eq. (10) for each pair of quarks of a given flavor.
Appendix B: Exact SU(4) limit The expectation value of the operator (1) displayed in Table  2, can be checked with the following formula [40] i< j where n is the number flavors and k the number of quarks, here n = 4 and k = 4. C S U(n) 2 is the Casimir operator eigenvalues of S U(n) which can be derived from the expression [41] : where f ′ i = f i − f n , for an irreducible representation given by the partition [ f 1 , f 2 , ..., f n ]. Eq. (B.9) has been previously used for n = 3 and k = 6 in Ref. [41].

Appendix C: The baryons
The masses of baryons related to the present study were calculated variationally using a radial wave function of the form φ ∝ exp[− a 2 (x 2 + y 2 )] containing the variational parameter a and the coordinates x and y defined by Eq. (11).
The results are indicated in Table 5 together with the experimental masses. One can see that condition V uc η ′ 0 is disfavored because in this case the mass difference between Σ * c and Σ c is negligible as compared to the experiment. This quantitatively justify the condition V uc η ′ = 0, imposed in the case of pentaquarks. (Note that the quantity V uc η ′ is not present in the mass of Λ c .) The calculated charmed baryon masses with V uc η ′ = 0 are somewhat smaller than the experimental values. By increasing the charmed quark mass from m c = 1.35 GeV to m c = 1.45 GeV the agreement with the experiment is much better. If the same increase is taken into account in calculating the pentaquark masses, the spectrum is shifted to larger masses but still remains at the edge of the desired experimental range. However the conclusion remains unchanged: the lowest state has positive parity.
In addition, we have calculated the masses of two doubly charmed baryons, namely Ξ ++ cc ( 1 2 ) and Ξ ++ cc ( 3 2 ). They are 3386 MeV and 3520 MeV respectively, when V uc η ′ = 0. From experimental point of view the situation remains controversial. Recently the LHCb collaboration announced the obsevation of the double charm baryon Ξ ++ cc in the Λ + c K − π + π + mass spectrum [42]. Its mass of 3621 MeV is at variance with the 3519 MeV value found earlier by the SELEX Collaboration [43]. The spin and parity are not known experimentally.

Appendix D: The flavor wave functions
Here we give the explicit form of flavor states of content uudc in the Young Yamanouchi basis of irreducible representations (irrep) [ f ] F appearing in the wave functions. We assume identical particle and use the method of Ref. [22] to derive these basis vectors. The order of particles is always 1234 in every term below and the permutation symmetry of each state is defined by its Yamanouchi symbol, which is a compact notation for a Young tableau. For a tableau with n particles it is defined by Y = (r n , r n−1 , ..., r 1 ) where r i represents the row of the particle i.
For the irrep [22] there are two basis vectors For the irrep [31] there are three basis vectors |[31] F 2111 = 1 6 (3uudc + 3duuc + 3uduc − cudu − cuud − cduu − ucdu − dcuu In the present model it is useful to rewrite the above wave functions in terms of products of symmetric φ [2] (q a q b ) = (q a q b + q b q a )/ √ 2 or antisymmetric φ [11] (q a q b ) = (q a q b − q b q a )/ √ 2 quark pair states for the pairs 12 and 34. This allows a forward calculation of the flavor integrated matrix elements (10) and in addition one can easily read off the isospin of the corresponding wave function.