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

TL;DR

This paper models the spectrum of the hidden charm pentaquark $uudcar c$ using an SU(4) flavor-spin hyperfine interaction, predicting the ordering of positive and negative parity states and comparing with LHCb experimental data.

Contribution

It introduces a detailed constituent quark model with specific orbital configurations to predict pentaquark state energies and parities, aligning theoretical results with experimental observations.

Findings

Lowest positive parity states are below certain negative parity states.

Positive parity states with specific quantum numbers are identified.

Model results are consistent with LHCb pentaquark structures.

Abstract

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.

Tables

Table 1. Table 1: Masses and decay widths of the 2019 LHCb resonances Skwarnicki ; Aaij:2019vzc .

Resonance	Mass (MeV)	Width (MeV)
$P_c^{+}\left(4312\right)$	4311.9$\pm$ 0.7 ${}_{-0.6}{}^{}{}_{}{}^{+6.8}$	9.8$\pm$ 2.7 ${}_{-4.5}{}^{}{}_{}{}^{+3.7}$
$P_c^{+}\left(4440\right)$	4440.3$\pm$ 1.3 ${}_{-4.7}{}^{}{}_{}{}^{+4.1}$	20.6$\pm$ 4.9 ${}_{-10.1}{}^{}{}_{}{}^{+8.7}$
$P_c^{+}\left(4457\right)$	4457.3$\pm$ 0.6 ${}_{-1.7}{}^{}{}_{}{}^{+4.1}$	6.4$\pm$ 2.0 ${}_{-1.9}{}^{}{}_{}{}^{+5.7}$

Table 2. Table 2: The simplified interaction - V χ subscript 𝑉 𝜒 V_{\chi} / C χ subscript 𝐶 𝜒 C_{\chi} , Eq. ( 1 ), integrated in the flavor-spin ( F ​ S 𝐹 𝑆 FS ) space for four quark subsystems in three different states defined in column 2 in the F ​ S 𝐹 𝑆 FS coupling. The notation is: O 𝑂 O for orbital, C 𝐶 C for color, F 𝐹 F for flavor, S 𝑆 S for spin.

State	Definition	L	S
$|1⟩$	$({\left[31\right]}_O{\left[211\right]}_C{\left[1^4\right]}_{OC};{\left[22\right]}_F{\left[22\right]}_S{\left[4\right]}_{FS})$	1	0	27
$|2⟩$	$({\left[31\right]}_O{\left[211\right]}_C{\left[1^4\right]}_{OC};{\left[31\right]}_F{\left[31\right]}_S{\left[4\right]}_{FS})$	1	1	21
$|3⟩$	${\left({\left[4\right]}_O{\left[211\right]}_C\left[211\right]\right]}_{OC};{\left[211\right]}_F{\left[22\right]}_S{\left[31\right]}_{FS})$	0	0	15

Table 3. Table 3: The hyperfine interaction V χ subscript 𝑉 𝜒 V_{\chi} , Eq. ( 7 ), integrated in the flavor-spin space, for the quark subsystem u ​ u ​ d ​ c 𝑢 𝑢 𝑑 𝑐 uudc I 𝐼 I =1/2, S 𝑆 S = 0. V γ q a ​ q b subscript superscript 𝑉 subscript 𝑞 𝑎 subscript 𝑞 𝑏 𝛾 {{V}^{q_{a}q_{b}}_{\gamma}} are defined in Eq. ( 10 ) where the upper index q a ​ q b subscript 𝑞 𝑎 subscript 𝑞 𝑏 {q_{a}q_{b}} indicates the flavor of the interacting quark pair .

State	$V_{\chi }$
$|{\left[31\right]}_O{\left[22\right]}_F{\left[22\right]}_S{\left[4\right]}_{FS}⟩$	15 $V_{\pi }-V_{\eta }^{uu}-2V_{{\eta }^{\prime }}^{uu}+12V_D^{uc}-\frac{1}{2}{V}_{{\eta }_c}^{uu}+\frac{3}{2}{V}_{{\eta }_c}^{uc}-2V_{{\eta }^{\prime }}^{uc}$
$|{\left[31\right]}_O{\left[31\right]}_F{\left[31\right]}_S{\left[4\right]}_{FS}⟩$	3 $V_{\pi }+V_{\eta }^{uu}+\frac{1}{3}{V}_{{\eta }^{\prime }}^{uu}+14V_D^{uc}+\frac{1}{2}{V}_{{\eta }_c}^{uu}+\frac{5}{2}{V}_{{\eta }_c}^{uc}-\frac{10}{3}{V}_{{\eta }^{\prime }}^{uc}$
$|{\left[4\right]}_O{\left[211\right]}_F{\left[22\right]}_S{\left[31\right]}_{FS}⟩$	7 $V_{\pi }-\frac{7}{9}{V}_{\eta }^{uu}-\frac{14}{9}{V}_{{\eta }^{\prime }}^{uu}++\frac{22}{3}{V}_D^{uc}-\frac{7}{18}{V}_{{\eta }_c}^{uu}+\frac{11}{6}{V}_{{\eta }_c}^{uc}-\frac{22}{9}{V}_{{\eta }^{\prime }}^{uc}$

Table 4. Table 4: Lowest positive and negative parity u ​ u ​ d ​ c ​ c ¯ 𝑢 𝑢 𝑑 𝑐 ¯ 𝑐 uudc\bar{c} pentaquarks of quantum numbers I , J P 𝐼 superscript 𝐽 𝑃 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.

State	$I$	$J^P$	Variational parameters		Mass	$⟨r_{12}⟩$	$⟨r_{45}⟩$
			a (fm-2)	b (fm-2)	(GeV)	fm	fm
$|1⟩|\overline{c}⟩$	$\frac{1}{2}$	${\frac{1}{2}}^{+}$, ${\frac{3}{2}}^{+}$	2.055	1.079	4.273	0.875	1.018
$|2⟩|\overline{c}⟩$	$\frac{1}{2}$, $\frac{3}{2}$	${\frac{1}{2}}^{+}$,  ${\frac{3}{2}}^{+}$,  ${\frac{5}{2}}^{+}$	1.541	1.027	4.453	1.010	1.085
$|3⟩|\overline{c}⟩$	$\frac{1}{2}$	${\frac{1}{2}}^{-}$	1.284	1.079	4.487	0.996	1.062

Table 5. Table 5: Variational solution for the baryon ground state masses with the flavor-spin interaction of Sec. 2 . Column 1 gives the baryon, column 2 the isospin, column 3 the spin and parity column 4 the mass with V η ′ u ​ c = 0 subscript superscript 𝑉 𝑢 𝑐 superscript 𝜂 ′ 0 V^{uc}_{\eta^{\prime}}=0 and V η ′ u ​ c ≠ 0 subscript superscript 𝑉 𝑢 𝑐 superscript 𝜂 ′ 0 V^{uc}_{\eta^{\prime}}\neq 0 , column 5 the variational parameter and the last column the experimental mass.

Baryon	$I$	$J^P$	Calculated Mass (GeV)		a(fm-2)	Exp.mass (GeV)
			$V_{{\eta }^{\prime }}^{uc}=0$	$V_{{\eta }^{\prime }}^{uc}\ne 0$
$N$	$\frac{1}{2}$	${\frac{1}{2}}^{+}$	0.960		2.594	0.940
$\Delta$	$\frac{3}{2}$	${\frac{3}{2}}^{+}$	1.304		2.594	1.232
${\Lambda }_c$	$0$	${\frac{1}{2}}^{+}$	2.180		2.055	2.283
${\Sigma }_c$	$1$	${\frac{1}{2}}^{+}$	2.346	2.434	2.055	2.455
${\Sigma }_c^{\ast }$	$1$	${\frac{3}{2}}^{+}$	2.482	2.438	2.055	2.530

Equations

$$V_{\chi }=-C_{\chi }\sum _{i<j}{\lambda }_i^F\cdot {\lambda }_j^F{\mathbf{\sigma }}_i\cdot {\mathbf{\sigma }}_j$$

$$\Delta E=E(qqqQ\overline{Q})-E(q^3)-E(Q\overline{Q})$$

$$\Delta E=\frac{5}{4}\hslash \omega -27/2C_{\chi }\approx -74\text{mbox}MeV,$$

$$\Delta E=3/4\hslash \omega -3/2C_{\chi }\approx 88\text{mbox}MeV,$$

$$H$$

$$V_{\text{conf}}(r_{ij})=-\frac{3}{8}{\lambda }_i^c\cdot {\lambda }_j^{c}C{r}_{ij}.$$

$$V_{\chi }(r_{ij})$$

$$V_{\gamma }(r)$$

$$\frac{g_{\pi q}^2}{4\pi }=\frac{g_{\eta q}^2}{4\pi }=\frac{g_{Dq}^2}{4\pi }=0.67,\frac{g_{{\eta }^{\prime }q}^2}{4\pi }=1.206,$$

$$r_0=0.43\text{mbox}fm,\alpha =2.91\text{mbox}{fm}^{-1},C=0.474\text{mbox}{fm}^{-2},$$

$${\mu }_{\pi }=139\text{mbox}MeV,{\mu }_{\eta }=547\text{mbox}MeV,{\mu }_{{\eta }^{\prime }}=958\text{mbox}MeV,$$

$${\mu }_D=1867\text{mbox}MeV.$$

$$m_{u,d}=340\text{mbox}MeV,m_c=1350\text{mbox}MeV.$$

$$V_{ij}={\bm{\sigma}}_{i}\ \cdot{\bm{\sigma}}_{j}\,\left\{\begin{array}[]{cl}V_{\pi}+\frac{1}{3}V^{uu}_{\eta}+\frac{1}{6}V^{uu}_{\eta_{c}},&\hskip 8.5359pt\mbox{ $[2]_{F},I=1$}\\ 2V_{K}-\frac{2}{3}V^{us}_{\eta},~{}~{}~{}2V^{uc}_{D}-\frac{1}{2}V^{uc}_{\eta_{c}}&\hskip 8.5359pt\mbox{ $[2]_{F},I=\frac{1}{2}$}\\ 2V^{sc}_{D_{s}}-\frac{1}{2}V^{sc}_{\eta_{c}}&\hskip 8.5359pt\mbox{ $[2]_{F},I=0$}\\ \frac{4}{3}V^{ss}_{\eta}+\frac{3}{2}V^{cc}_{\eta_{c}}&\hskip 8.5359pt\mbox{ $[2]_{F},I=0$}\\ -2V^{sc}_{D_{s}}-\frac{1}{2}V^{sc}_{\eta_{c}}&\hskip 8.5359pt\mbox{ $[11]_{F},I=0$}\\ -2V_{K}-\frac{2}{3}V^{us}_{\eta},~{}~{}~{}-2V^{uc}_{D}-\frac{1}{2}V^{uc}_{\eta_{c}}&\hskip 8.5359pt\mbox{ $[11]_{F},I=\frac{1}{2}$}\\ -3V_{\pi}+\frac{1}{3}V^{uu}_{\eta}+\frac{1}{6}V^{uu}_{\eta_{c}},&\hskip 8.5359pt\mbox{ $[11]_{F},I=0$}\end{array}\right.$$

$$\displaystyle\begin{array}[]{c}\bm{x}\ =\ {\bm{r}}_{1}\ -\ {\bm{r}}_{2}\ ,\\ \bm{y}\ =\ {\left({{\bm{r}}_{1}\ +\ {\bm{r}}_{2}\ -\ 2{\bm{r}}_{3}}\right)/\sqrt{3}},\\ \bm{z}\ =\ {\left({{\bm{r}}_{1}\ +\ {\bm{r}}_{2}\ +\ {\bm{r}}_{3}\ -\ 3{\bm{r}}_{4}}\right)/\sqrt{6}},\\ \bm{t}\ =\ {\left({{\bm{r}}_{1}\ +\ {\bm{r}}_{2}\ +\ {\bm{r}}_{3}+\ {\bm{r}}_{4}-\ 4{\bm{r}}_{5}}\right)/\sqrt{10}},\end{array}$$

$$\displaystyle{\psi}_{1}=\begin{array}[]{c}$\framebox{1}\framebox{2}\framebox{3}$\\ $\framebox{4}$\end{array}=\left\langle{\vec{x}\left|{000}\right.}\right\rangle\left\langle{\vec{y}\left|{000}\right.}\right\rangle\left\langle{\vec{z}\left|{010}\right.}\right\rangle$$

$$\displaystyle{\psi}_{2}=\begin{array}[]{c}$\framebox{1}\framebox{2}\framebox{4}$\\ $\framebox{3}$\end{array}=\left\langle{\vec{x}\left|{000}\right.}\right\rangle\left\langle{\vec{y}\left|{010}\right.}\right\rangle\left\langle{\vec{z}\left|{000}\right.}\right\rangle$$

$$\displaystyle{\psi}_{3}=\begin{array}[]{c}$\framebox{1}\framebox{3}\framebox{4}$\\ $\framebox{2}$\end{array}=\left\langle{\vec{x}\left|{010}\right.}\right\rangle\left\langle{\vec{y}\left|{000}\right.}\right\rangle\left\langle{\vec{z}\left|{000}\right.}\right\rangle$$

$${\psi }_1^5=N\exp \left[-\frac{a}{2}\left(x^2+y^2+z^2\right)-\frac{b}{2}{t}^2\right]z{Y}_{10}\left(\hat{z}\right)$$

$${\psi }_2^5=N\exp \left[-\frac{a}{2}\left(x^2+y^2+z^2\right)-\frac{b}{2}{t}^2\right]y{Y}_{10}\left(\hat{y}\right)$$

$${\psi }_3^5=N\exp \left[-\frac{a}{2}\left(x^2+y^2+z^2\right)-\frac{b}{2}{t}^2\right]x{Y}_{10}\left(\hat{x}\right)$$

$$N=\frac{2^{3/2}{a}^{11/4}{b}^{3/4}}{3^{1/2}{\pi }^{5/2}}$$

$${\varphi }_0=N_0\exp \left[-\frac{a}{2}\left(x^2+y^2+z^2\right)-\frac{b}{2}{t}^2\right],$$

$$N_0=(\frac{a}{\pi }{)}^{9/4}(\frac{b}{\pi }{)}^{3/4}.$$

$$\displaystyle\begin{array}[]{lcl}\left\langle{T}\right\rangle&=&\ {\frac{1}{3}}\ \left[{\left\langle{{\psi}_{1}^{5}\left|{T}\right|{\psi}_{1}^{5}}\right\rangle\ +\ \left\langle{{\psi}_{2}^{5}\left|{T}\right|{\psi}_{2}^{5}}\right\rangle\ +\ \left\langle{{\psi}_{3}^{5}\left|{T}\right|{\psi}_{3}^{5}}\right\rangle}\right]\\ \\ &=&\ {\hbar}^{2}\ \left({{\frac{11}{2{\mu}_{1}}}\ a\ +\ {\frac{3}{{2\mu}_{2}}}\ b}\right),\end{array}$$

$$\displaystyle{\frac{4}{{\mu}_{1}}}\ =\ \frac{3}{{m}_{q}}\ +\ \frac{1}{{m}_{Q}},\$$

$$\frac{5}{{\mu }_2}=\frac{1}{{\mu }_1}+\frac{4}{m_Q},$$

$$⟨V_{conf}⟩=\frac{C}{2}\left(6⟨r_{12}⟩+4⟨r_{45}⟩\right)$$

$$⟨r_{ij}⟩=\frac{1}{3}\left[⟨{\psi }_1^5\left|r_{ij}\right|{\psi }_1^5⟩+⟨{\psi }_2^5\left|r_{ij}\right|{\psi }_2^5⟩+⟨{\psi }_3^5\left|r_{ij}\right|{\psi }_3^5⟩\right],$$

$$⟨r_{12}⟩=\frac{20}{9}\sqrt{\frac{1}{\pi a}},$$

Full text

††journal: Eur. Phys. J. C

1 Introduction

A considerable wave of interest has been raised by the 2015 observation of the LHCb collaboration of two resonances $P^{+}_{c}(4380)$ ($\Gamma=205\pm 18\pm 86$ MeV) and $P^{+}_{c}(4450)$ ($\Gamma=39\pm 5\pm 19$ MeV) in the $\Lambda^{0}_{b}\rightarrow J/\psi K^{-}p$ decay, interpreted as hidden charm pentaquarks of structure $uudc\bar{c}$ Aaij:2015tga . The preferred quantum numbers, which gave the best solution, were $3/2^{-},5/2^{+}$. However, the quantum numbers $3/2^{+}$, $5/2^{-}$ and $5/2^{-}$, $3/2^{+}$ were also acceptable. Various possible interpretations of the 2015 LHCb resonances as kinematical effects, molecular states or compact pentaquarks were reviewed, see for example Refs. Ali:2017jda , Karliner:2017qhf and Liu:2019zoy .

Recently the LHCb Collaboration have updated their analysis of the $\Lambda^{0}_{b}\rightarrow J/\psi K^{-}p$ decay Skwarnicki ; Aaij:2019vzc 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\sigma$ statistical significance and the other two resonances replace the previous $P^{+}_{c}(4450)$ now resolved at $5.4\sigma$ significance. The broad $P^{+}_{c}(4380)$ resonance observed in 2015 awaits 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 $\Sigma^{+}_{c}\overline{D}^{0}$ threshold (4318 MeV) and the masses of $P^{+}_{c}(4440)$ and $P^{+}_{c}(4457)$ lie a few MeV below the $\Sigma^{+}_{c}\overline{D}^{*0}$ threshold (4460 MeV). Proximity of the $\Sigma^{+}_{c}\overline{D}^{0}$ and $\Sigma^{+}_{c}\overline{D}^{*0}$ thresholds to these narrow peaks made their interpretation as molecular S-wave of $\Sigma^{+}_{c}\overline{D}^{0}$ or $\Sigma^{+}_{c}\overline{D}^{*0}$ quite natural Guo:2019kdc ; Guo:2019fdo ; Xiao:2019mst ; Xiao:2019aya ; Shimizu:2019ptd ; Cheng:2019obk . In such an interpretation, they all acquire a negative parity. In Ref. Burns:2019iih it is shown that the usual molecular scenarios cannot explain the production rate of $P_{c}$ states in $\Lambda_{b}$ decays. Moreover, the paper points out the need to couple the $\Sigma^{+}_{c}\overline{D}^{*0}$ and the $\Lambda_{c}(2595)\overline{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)={\frac{3}{2}}^{-}$ and $J^{P}(4457)={\frac{1}{2}}^{+}$ and respectively.

The 2019 LHCb pentaquarks have also been analysed in the compact diquark model Ali:2019npk . This analysis has been extended in Ref. Ali:2019clg , 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 DeRujula:1975qlm . The basis contains only the color $\bar{3}$ 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 Glozman:1995fu ; Robson:1988mp ; Glozman:1997ag , 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 $uu\bar{c}\bar{c}$, the pentaquarks $uudd\bar{c}$, $uuds\bar{c}$ and the hexaquark $uuddsc$ Stancu:1998pr .

In Refs. Glozman:1995fu ; Glozman:1997ag 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. Glozman:1995fu ; Glozman:1997ag 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. Glozman:1995xy 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. Glozman:1995fu 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 Stancu:1991rc .

Presently we study the pentaquarks of structure $uudc\bar{c}$ in order to see whether or not the FS model can accommodate the new resonances observed at LHCb Aaij:2019vzc . The approach is similar to that previously used in Ref. Stancu:1998sm for the positive parity $uudd\bar{c}$ and $uuds\bar{c}$ pentaquarks. That work preceded the H1 experiment Aktas:2004qf which has found evidence, so far not confirmed by other experiments, for a narrow baryon resonance in the $D^{*\pm}p$ invariant mass, interpreted as a $uudd\bar{c}$ pentaquark. The mass prediction made in Ref. Stancu:1998sm is not far from the observation of the H1 experiment.

The parity of the pentaquark is given by P$\ ={\left({-}\right)}^{{\ell\ +\ 1}}$. For the lowest positive parity states (see below) the subsystem of four quarks is defined to carry an angular momentum $\ell$ = 1 in their internal motion. This implies that this subsystem must be in a state of orbital symmetry ${\left[{31}\right]}_{O}$. Although the kinetic energy of such a state is higher than that of the totally symmetric ${\left[{4}\right]}_{O}$ state, a crude estimate based on the simplified interaction (1), suggests that the ${\left[{31}\right]}_{O}$ symmetry leads to a stable pentaquark of positive parity and the lowest state of symmetry ${\left[{4}\right]}_{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

[TABLE]

with $\lambda_{i}^{F}$ the Gell-Mann matrices, $\bm{\sigma}_{i}$ the Pauli matrices and $C_{\chi}$ an equal strength constant for all pairs.

Let us consider two totally antisymmetric states with the orbital symmetry ${\left[{31}\right]}_{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. Burns:2015dwa .

The expectation value of (1), as shown in Table 2 (see the next section for the derivation of the last column), is $-27\ C_{\chi}$ for $\left|{\left.{1}\right\rangle}\right.$ and $-21\ C_{\chi}$ for $\left|{\left.{2}\right\rangle}\right.$.

First we consider the lowest state, i.e. $\left|{\left.{1}\right\rangle}\right.$. The quark-antiquark 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 $Q\overline{Q}$, where $Q$ stands for a heavy flavor. Then, in discussing the stability we introduce the quantity

[TABLE]

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 $\Delta E$. Then, the kinetic energy contribution   to   $\Delta E$   is   $\Delta KE\ =\ 5/4\ \hbar\omega$ in a harmonic oscillator model HO .

Table 2 shows that the FS contribution to the state $\left|{\left.{1}\right\rangle}\right.$ is $-27\ C_{\chi}$. For the nucleon or $\Lambda_{c}$ it is - 27/2 $C_{\chi}$ (see next section). Then we have $\Delta V_{\chi}=\ -27/2\ C_{\chi}$. With $\hbar\omega\approx$ 157 MeV and $C_{\chi}\approx$ 20 MeV Glozman:1995xy this gives

[TABLE]

i.e. a considerable binding. This is to be contrasted with the negative parity pentaquarks similar to those containing light quarks, studied in Ref. Genovese:1997tm . For the lowest negative parity state $\left|{\left.{3}\right\rangle}\right.$ of $uudc\bar{c}$ one gets

[TABLE]

i.e. instability. Thus the FS interaction overcomes the excess of kinetic contribution in ${\left[{31}\right]}_{O}$ and generates a lower expectation value for ${\left[{31}\right]}_{O}$ than for ${\left[{4}\right]}_{O}$.

The paper is organized as follows. In Sec. 2 we introduce the model Hamiltonian and generalize the two-body 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. A gives the matrices $\lambda^{F}$ of SU(4). B is a reminder of useful group theory formulae for SU(n). C exhibits the variational solution for the baryon masses relevant for the present study. In D we construct the explicit form of flavor states of content $uudc$ in the Young Yamanouchi basis for specific irreducible representations $[{f}]_{F}$.

2 The Hamiltonian

The nonrelativistic FS Hamiltonian has the general form Glozman:1995fu

[TABLE]

with $m_{i}$ and ${\bm{p}}_{i}$ denoting the quark masses and momenta respectively and $r_{ij}$ the distance between the interacting quarks $i$ and $j$. The Hamiltonian contains the internal kinetic energy and the linear confining interaction

[TABLE]

The hyperfine part $V_{\chi}(r_{ij})$ has a flavor-spin structure which presently is extended to SU(4), similarly to Refs. Pepin:1996id and Glozman:1995xy , except that in the latter, the term $V_{\eta^{\prime}}\lambda^{0}_{i}\cdot\lambda^{0}_{j}$ of Eq. (7) has been ignored. One has

[TABLE]

with $\lambda^{F}_{i}$ ($F$ = 1,2,…,15) defined in A and $\lambda^{0}_{i}=\sqrt{2/3}~{}{\bf 1}$, where $\bf 1$ is the $4\times 4$ unit matrix. Note that the norm of $\lambda^{0}_{i}$ is not that satisfying the general property tr$(\lambda^{a}\lambda^{b})$ = 2 $\delta^{ab}$ if $\bf 1$ is the $4\times 4$ unit matrix Stancu:1998sm , but we maintain the norm from the SU(3) version of the flavor-spin interaction in order to use the same coupling constant for $\eta$ and $\eta^{\prime}$, which gave a good spectrum for strange baryons. At this stage it is useful to mention the study of Ref. Shuryak:1997xd where it has been shown that $\eta^{\prime}$ has a significant charm component.

In the SU(4) version the interaction (7) contains $\gamma=\pi,K,\eta,D,D_{s},\eta_{c}$ and $\eta^{\prime}$ meson-exchange terms. Every $V_{\gamma}(r_{ij})$ 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 Glozman:1996wq . For a given meson $\gamma$ the meson exchange potential is

[TABLE]

In the present calculations we use the parameters of Ref. Glozman:1996wq to which we add the $\mu_{D}$ mass and the coupling constant $\frac{g_{Dq}^{2}}{4\pi}$. These are

[TABLE]

[TABLE]

[TABLE]

[TABLE]

The meson masses correspond to the experimental values from the Particle Data Group Tanabashi:2018oca . There is no $K$- and $D_{s}$-exchange in the $uudc\bar{c}$ pentaquark and, as discussed in the following, we ignore $\eta_{c}$-exchange. Therefore the $K$, $D_{s}$ and $\eta_{c}$ masses are not needed in these calculations.

For the quark masses we take the values determined variationally in Ref. Pepin:1996id

[TABLE]

They were adjusted to reproduce the average mass ${\overline{M}}=(M+3M^{*})/4$ = 2008 MeV of the $D$ mesons. With these values one obtains an average mass of 2985 MeV for $J/\psi$ and $\eta_{c}$.

After integrating in the flavor space, the two-body matrix elements which generalize Eq. (3.3) of Ref. Glozman:1995fu to SU(4), one obtains a result containing new contributions due to charm. This is

[TABLE]

In Eqs. (10) the pair of quarks $ij$ 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. Glozman:1995fu the upper index of $\pi$ 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 $\frac{1}{2}V^{uc}_{\eta_{c}}$ was missing in Ref. Pepin:1996id . 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 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 $V_{\pi}$ = $V^{uu}_{\eta}$ = $V^{uu}_{\eta_{c}}$ = $V^{uc}_{D}$ = $V^{uc}_{\eta_{c}}$ = - ${C}_{\chi}$ and $V^{uu}_{\eta^{\prime}}$ = $V^{uc}_{\eta^{\prime}}$ = 0.

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}_{\chi}$ and 5/2 ${C}_{\chi}$ respectively. Then, making use of the spin matrix elements $\langle\sigma_{i}\cdot\sigma_{j}\rangle$ = 1 and $\langle\sigma_{i}\cdot\sigma_{j}\rangle$ = - 3 for $S=1$ and $S=0$ respectively, one finds that the contribution of the FS interaction to the baryon masses are $V_{\chi}(N)$ = $V_{\chi}(\Lambda_{c}$) = - 27/2 $C_{\chi}$ and $V_{\chi}(\Delta)$ = - 9/2 $C_{\chi}$, which are close to the SU(3) limits

• 14 $C_{\chi}$ and - 4 $C_{\chi}$ respectively. As a matter of fact one can also see that the SU(4) limit given in - $C_{\chi}$ 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 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.

3 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. Stancu:1998sm .

One needs four internal Jacobi coordinates for pentaquarks. The convenient choice is

[TABLE]

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 $\vec{x},\vec{y},\vec{z}$ for the specific permutation symmetry ${\left[{31}\right]}_{O}$. For the lowest positive parity states we have considered an $s^{3}p$ structure for ${\left[{31}\right]}_{O}$ and extended Moshinski’s method MO69 , 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 $\psi_{i}$, defined as Young-Yamanouchi basis vectors Stancu:1991rc of the irreducible representation ${\left[{31}\right]}_{O}$ in terms of shell model states $\left\langle{\vec{r}\left|{n\ell m}\right.}\right\rangle$ where $n$ = 0, $\ell$ = 1, and, for simplicity, we took m = 0 everywhere. They read

[TABLE]

[TABLE]

[TABLE]

One can see that the $\ell$ = 1 orbital excitation is located in the relative motion between quarks defined by the internal coordinates $\vec{x}$, $\vec{y}$ and $\vec{z}$. The states ${\psi}_{1}$ and ${\psi}_{2}$ are the generalization, from three to four quarks, of $\psi^{\lambda}_{\ell 0}$ and ${\psi}_{3}$ the generalization, from three to four quarks, of $\psi^{\rho}_{\ell 0}$ radial wave functions of mixed symmetry $[21]$ used in baryon spectroscopy. Thus, a more appealing notation could be $\psi^{\lambda_{1}}_{\ell 0}$, $\psi^{\lambda_{2}}_{\ell 0}$ and $\psi^{\rho}_{\ell 0}$ instead of ${\psi}_{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 $\vec{t}$. Then the pentaquark orbital wave functions $\psi_{i}^{5}$ are obtained by multiplying each $\psi_{i}$ from above by the wave function $\left\langle{\vec{t}\left|{000}\right.}\right\rangle$ which describes the relative motion between the four-quark subsystem and the antiquark $\overline{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 $\overline{c}$. We explicitly have

[TABLE]

[TABLE]

[TABLE]

where

[TABLE]

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

[TABLE]

with

[TABLE]

4 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

[TABLE]

with

[TABLE]

and

[TABLE]

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 $\left\langle{T}\right\rangle\ =\ {\frac{7}{2}}\ \hbar\omega$ with $\hbar\omega\ =\ 2\ a{\hbar}^{2}/m$, see HO . Note that one can obtain the desired units of all quantities involved, by muliplying/dividing by $\hbar c\simeq$ 197.33 MeV fm at the correct power.

By integrating in the color space, the expectation value of the confinement interaction (6) becomes

[TABLE]

where the coefficients 6 and 4 account for the number of quark-quark and quark-antiquark pairs, respectively, and

[TABLE]

where $i,j$ = 1,2,3,4,5 ($i\neq j$). An analytic evaluation gives

[TABLE]

as well as for ${r}_{ij}$ with $i,j$ = 1,2,3,4 and

[TABLE]

and likewise for ${r}_{j5}$ where $j$ = 1,2,3.

One also needs to derive the expressions of $\left\langle{{r}_{12}}\right\rangle$ and $\left\langle{{r}_{45}}\right\rangle$ in the case where the four quarks are in the $s^{4}$ configuration. These are

[TABLE]

and

[TABLE]

The matrix elements of the spin-flavor operators of Eq. (7) have been calculated using the fractional parentage technique described in Ref. Stancu:1991rc based on Clebsch-Gordan coefficients of the group $S_{4}$ Stancu:1999qr . 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 $\left\langle{ss\left|{{V}^{q_{a}q_{b}}_{\gamma}}\right|ss}\right\rangle\ ,\ \left\langle{sp\left|{{V}^{q_{a}q_{b}}_{\gamma}}\right|sp}\right\rangle$ and $\left\langle{sp\left|{{V}^{q_{a}q_{b}}_{\gamma}}\right|ps}\right\rangle$ where $\gamma=\pi,D,\eta,\eta^{\prime}$ 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 $\pi$ meson exchange interaction $V_{\pi}$, see Eq. (10), turned out to be about 0.25, in agreement with the mass ratio $m_{u,d}/m_{c}$ of Eq. (9).

5 Results

The four quark states described in the previous sections couple to the antiquark to a total angular momentum $\vec{J}\ =\ \vec{L}\ +\ \vec{S}\ +\ {\vec{s}}_{Q}$, with $\vec{L}$, $\vec{S}$ the angular momentum and spin of the four-quark cluster and $\vec{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 $|\bar{c}\rangle$. The orbital wave functions of the pentaquark are given by Eqs. (25)-(30) containing the variational parametres $a$ and $b$. We have neglected the contribution of $V^{uu}_{\eta_{c}}$ and of $V^{uc}_{\eta_{c}}$, because little $u\bar{u}$ and $d\bar{d}$ is expected in ${\eta_{c}}$. We have also neglected $V^{uc}_{\eta^{\prime}}$ assuming a little $c\bar{c}$ component in $\eta^{\prime}$, which means that we took

[TABLE]

The optimal values found for the parameters $a$ and $b$ are indicated in each case. The eigenvalues of $|1\rangle|\bar{c}\rangle$ and $|2\rangle|\bar{c}\rangle$ states are degenerate for the allowed values of $J$ in each case. Moreover, for the wave function $|2\rangle|\bar{c}\rangle$ the states with $J^{P}=\frac{1}{2}^{+}$ and $\frac{3}{2}^{+}$ have multiplicity 2. Presently there are more states than the observed ones. In the case of pentaquarks with heavy flavor, Ref. Ali:2019npk argues that only part of the spectrum is reachable in $\Lambda_{b}$ decays.

The SU(4) FS model described in Sec. 2 gives a value for the mass of the state $|1\rangle|\bar{c}\rangle$ of 4273 MeV, close to the mass of the $P^{+}_{c}(4312)$ resonance and supports the assignment $J^{P}=\frac{1}{2}^{+}$ or $\frac{3}{2}^{+}$ for it. The state $|2\rangle|\bar{c}\rangle$ has a mass of 4453 MeV close to that of the $P^{+}_{c}(4440)$ and $P^{+}_{c}(4457)$ resonances and supports the quantum numbers $J^{P}=\frac{1}{2}^{+},\frac{3}{2}^{+}$ or $\frac{5}{2}^{+}$.

As discussed above, in the exact limit one obtains binding for $|1\rangle|\bar{c}\rangle$, 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\eta_{c}(3920)$ or $NJ/\psi(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 $\eta$, $\eta^{\prime}$ 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\rangle|\bar{c}\rangle$ 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 $\left\langle{{r}_{12}}\right\rangle$ and $\left\langle{{r}_{45}}\right\rangle$ from Eqs. (38)-(42). They are exhibited in Table 4. In the $|1\rangle|\overline{c}\rangle$ 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 $\left\langle{{r}_{12}}\right\rangle$ and $\left\langle{{r}_{45}}\right\rangle$ 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 deuteron-like. 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.

6 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. Burns:2015dwa 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/\psi$ final state. For the lowest state, $|{{\left[{31}\right]}_{O}\ {\left[{22}\right]}_{F}{\left[{22}\right]}_{S}{\left[{4}\right]}_{FS}}\rangle$, the overlap with $q^{3}$ + $Q\bar{Q}$ may be small. This problem deserves a separate investigation.

7 Comparison with other studies

The basis (25) - (27) is entirely consistent with the parity considerations made in the context of the tetrahedral group $T_{d}$ Bijker:2003pm , isomorphic with the permutation group $S_{4}$. The states $|1\rangle$ and $|2\rangle$ of Table 2 are states of symmetry $1^{-}_{A_{1}}$ in the notation of Ref. Bijker:2003pm . Likewise the state $|3\rangle$ of the same table has symmetry $0^{+}_{A_{1}}$. There is agreement between our findings and the conclusion of Ref. Bijker:2003pm that the parity of a pentaquark depends on the interplay between the flavor-spin 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 four-quark subsystem containing a charmed quark, as it was the case for pentaquarks containing only light flavors like $uudd$ Stancu:1998sm , so that the lowest state has positive parity in both cases, either for light or light+heavy quarks.

The authors of Ref. Yuan:2012wz have studied the $qqqc\bar{c}$ pentaquark prior to the 2015 LHCb observation of pentaquarks, in three different models, including the flavor-spin (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. Yuan:2012wz 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

[TABLE]

where

[TABLE]

Then the expectation value of the kinetic energy becomes

[TABLE]

with

[TABLE]

and

[TABLE]

The expectation value of the confinement interaction is given by Eq. (37) with

[TABLE]

and

[TABLE]

The hyperfine interaction integrated in the flavor-spin space is given by the third row of Table 3 because the symmetry state corresponds to $|{{\left[{4}\right]}_{O}\ {\left[{211}\right]}_{F}{\left[{22}\right]}_{S}{\left[{31}\right]}_{FS}}\rangle$.

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 GeV2, $b$ = 0.037 GeV2) 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. Yuan:2012wz is lower than the lowest $1/2^{+}$ state.

In Ref. Yang:2015bmv the spectrum of the $uudc\bar{c}$ 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. Giron:2019cfc . 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.

8 Conclusions

We have calculated the lowest part of the mass spectrum of the hidden charm pentaquark $uudc\bar{c}$, 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. Yuan:2012wz , 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 flavor-spin model is that it introduces an isospin dependence of the pentaquark states, necessary to discriminate between decay channels, like for tetraquarks Giron:2019cfc .

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. Glozman:1997jy 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. Ali:2019npk 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 $\ell=1$ states to be higher than the $\ell=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 Eides:2019tgv . In that model the new $P^{+}_{c}(4312)$ resonance is interpreted as a bound state of $\chi_{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 $uudc\bar{c}$ 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 triquark — light diquark model, the molecular picture or the hadronic molecular picture.

Appendix A The SU(4) generators

Here we reproduce the $\lambda^{F}$ matrices which are the SU(4) generators in the fundamental representation of SU(4) Stancu:1991rc . 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.

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

Appendix B Exact SU(4) limit

The expectation value of the operator (1) displayed in Table 2, can be checked with the following formula Ortiz-Pacheco:2018ccl

[TABLE]

where $n$ is the number flavors and $k$ the number of quarks, here $n$ = 4 and $k$ = 4. $C_{2}^{SU(n)}$ is the Casimir operator eigenvalues of $SU(n)$ which can be derived from the expression Stancu:1997dq :

[TABLE]

where $f_{i}^{\prime}=f_{i}-f_{n}$, for an irreducible representation given by the partition $[f_{1},f_{2},...,f_{n}]$. Eq. (120) has been previously used for $n$ = 3 and $k$ = 6 in Ref. Stancu:1997dq .

Appendix C The baryons

The masses of baryons related to the present study were calculated variationally using a radial wave function of the form $\phi\propto exp[-\frac{a}{2}(x^{2}+y^{2})]$ containing the variational parameter $a$ and the coordinates $x$ and $y$ defined by Eq. (15).

The results are indicated in Table 5 together with the experimental masses. One can see that condition $V^{uc}_{\eta^{\prime}}\neq 0$ is disfavored because in this case the mass difference between $\Sigma_{c}^{*}$ and $\Sigma_{c}$ is negligible as compared to the experiment. This quantitatively justify the condition $V^{uc}_{\eta^{\prime}}=0$, imposed in the case of pentaquarks. (Note that the quantity $V^{uc}_{\eta^{\prime}}$ is not present in the mass of $\Lambda_{c}$.)

The calculated charmed baryon masses with $V^{uc}_{\eta^{\prime}}=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 $\Xi^{++}_{cc}({\frac{1}{2}})$ and $\Xi^{++}_{cc}({\frac{3}{2}})$. They are 3386 MeV and 3520 MeV respectively, when $V^{uc}_{\eta^{\prime}}$ = 0. From experimental point of view the situation remains controversial. Recently the LHCb collaboration announced the obsevation of the double charm baryon $\Xi^{++}_{cc}$ in the $\Lambda^{+}_{c}K^{-}\pi^{+}\pi^{+}$ mass spectrum Aaij:2017ueg . Its mass of 3621 MeV is at variance with the 3519 MeV value found earlier by the SELEX Collaboration Mattson:2002vu . 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. Stancu:1991rc 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

[TABLE]

[TABLE]

For the irrep $[31]$ there are three basis vectors

[TABLE]

[TABLE]

[TABLE]

Note that there is a conflict with the basis vectors (A.9)-(A.11) of Ref. Yuan:2012wz which do not form a Young Yamanouchi basis, thus, although orthogonal, they do not satisfy the correct permutation symmetry.

In the same way the three basis vectors of the irrep $[211]$ were obtained as

[TABLE]

[TABLE]

[TABLE]

In the present model it is useful to rewrite the above wave functions in terms of products of symmetric $\phi_{[2]}(q_{a}q_{b})=(q_{a}q_{b}+q_{b}q_{a})/\sqrt{2}$ or antisymmetric $\phi_{[11]}(q_{a}q_{b})=(q_{a}q_{b}-q_{b}q_{a})/\sqrt{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.

In particular the states (122) and (123) become

[TABLE]

and

[TABLE]

where (131) obviously has isospin $I=1/2$ which means that the pairs 12 and 34 in (130) have to couple to the same isospin value as well.

For the irrep $[31]$ one has to use the Rutherford-Young-Yamanouchi basis, which symmetrizes or antisymmetrizes the last two particles Stancu:1991rc . Accordingly the basis vectors (124) and (125) transform into

[TABLE]

where the pair 34 is in a symmetric state, and

[TABLE]

where the pair 34 is in an antisymmetric state.

The state (126) can simply be rewritten as

[TABLE]

where the pair 12 is in an antisymmetric state and 34 in a symmetric state. The states (132)-(134) can have $I$ = 1/2 or 3/2.

For the irrep $[211]$ the vector (127) can be directly rewritten as

[TABLE]

and the mixture of the other two in the Rutherford-Young-Yamanouchi basis reads

[TABLE]

[TABLE]

The states (135)-(137) have $I$ = 1/2.

Acknowledgements.

I thank Ileana Guiasu for stimulating discussion. This work has been supported by the Fonds de la Recherche Scientifique - FNRS, Belgium, under the grant number 4.4503.19.