$$\Xi _c$$ and $$\Xi _b$$ excited states within a $$\mathrm{SU(6)}_{\mathrm{lsf}}\times $$HQSS model

We study odd parity J = 1 / 2 and J = 3 / 2 Ξ c resonances using a unitarized coupled-channel framework based on a SU (6) lsf × HQSS-extended Weinberg–Tomozawa baryon–meson interaction, while paying a special attention to the renormalization procedure. We predict a large molecular Λ c K ¯ component for the Ξ c (2790) with a dominant 0 - light-degree-of-freedom spin configuration. We discuss the differences between the 3 / 2 - Λ c (2625) and Ξ c (2815) states, and conclude that they cannot be SU(3) siblings, whereas we predict the existence of other Ξ c -states, one of them related to the two-pole structure of the Λ c (2595) . It is of particular interest a pair of J = 1 / 2 and J = 3 / 2 poles, which form a HQSS doublet and that we tentatively assign to the Ξ c (2930) and Ξ c (2970) , respectively. Within this picture, the Ξ c (2930) would be part of a SU(3) sextet, containing either the Ω c (3090) or the Ω c (3119) , and that would be completed by the Σ c (2800) . Moreover, we identify a J = 1 / 2 sextet with the Ξ b (6227) state and the recently discovered Σ b (6097) . Assuming the equal spacing rule and to complete this multiplet, we predict the existence of a J = 1 / 2 Ω b odd parity state, with a mass of 6360 MeV and that should be seen in the Ξ b K ¯ channel.

The study of heavy baryons with charm or bottom content has been the subject of much interest over the past years in view of newly discovered states [[1]]. In particular, there has been a tremendous effort to understand the nature of the experimental states within conventional quarks models, QCD sum-rules frameworks, QCD lattice analysis or molecular baryon–meson models (see Refs. [[2]–[7]] for recent reviews).

The attention has been recently revived by the experimental observation of several excited states. Recent detections have been reported by the LHCb Collaboration regarding five ${\mathrm{\Omega }}_c$ excited states in the ${\mathrm{\Xi }}_c^{+}{K}^{-}$ spectrum in pp collisions [[8]], and the excited ${\mathrm{\Xi }}_b(6227)$ state in ${\mathrm{\Lambda }}_b^0{K}^{-}$ and ${\mathrm{\Xi }}_b^0{\pi }^{-}$ invariant mass spectra also in pp collisions [[9]]. Moreover, the Belle Collaboration has confirmed the observation of four of the excited ${\mathrm{\Omega }}_c$ states [[10]], and detected the ${\mathrm{\Xi }}_c(2930)$ state in its decay to ${\mathrm{\Lambda }}_c^{+}{K}^{-}$ in $B^{-}\to K^{-}{\mathrm{\Lambda }}_c^{+}{\overline{\mathrm{\Lambda }}}_c^{-}$ decays [[11]].

In view of these new observations, a large theoretical activity has been indeed triggered, in particular within dynamical approaches based on a molecular description of these states. Starting from the newly observed ${\mathrm{\Omega }}_c$ states, the molecular models of Refs. [[12]] have been reanalyzed in view of the new discoveries. While in Ref. [[14]] two ${\mathrm{\Omega }}_c$ resonant states at 3050 MeV and 3090 MeV with $J^P=1/2^{-}$ were obtained, being identified with two of the experimental states, in Ref. [[15]] two $J^P=1/2^{-}$ ${\mathrm{\Omega }}_c$ states and one $J^P=3/2^{-}$ ${\mathrm{\Omega }}_c$ were determined within an extended local hidden gauge approach, the first two in good agreement with [[14]]. Other theoretical works also examined the ${\mathrm{\Omega }}_c$ sector, trying to explain the extra broad structure observed by the LHCb around 3188 MeV [[8]]. In Ref. [[16]] it was shown that this bump could be interpreted as the superposition of two $D\mathrm{\Xi }$ bound states, whereas in Ref. [[17]] a loosely bound molecule of mass 3140 MeV was determined.

Excited ${\mathrm{\Xi }}_c$ states below 3 GeV and the excited ${\mathrm{\Xi }}_b$ state found experimentally [[1]]. We show the assigned $J^P$ (when possible), the mass $\mathrm{M}$ and width $\mathrm{\Gamma }$ , as well as the decay channels

With regards to ${\mathrm{\Xi }}_c$ , the theoretical analysis based on the local hidden gauge formalism has shown that not only the ${\mathrm{\Xi }}_c(2930)$ can have a molecular interpretation, but also other ${\mathrm{\Xi }}_c$ states around 3 GeV reported in the PDG [[1]]. In particular, the ${\mathrm{\Xi }}_c(2790)$ would be a $J^P=1/2^{-}$ molecular state, whereas ${\mathrm{\Xi }}_c(2930)$ , ${\mathrm{\Xi }}_c(2970)$ , ${\mathrm{\Xi }}_c(3055)$ and ${\mathrm{\Xi }}_c(3080)$ could be described as molecules with either $1/2^{-}$ or $3/2^{-}$ [[18]]. On the other hand, the same model has produced two states for ${\mathrm{\Xi }}_b(6227)$ with masses close to the experimental one with similar widths, being the spin–parity assignment either $1/2^{-}$ or $3/2^{-}$ [[18]]. The ${\mathrm{\Xi }}_b$ state has been also studied within a unitarized model that uses the leading-order chiral Lagrangian in Refs. [[19]], identifying the ${\mathrm{\Xi }}_b(6227)$ state as a S-wave ${\mathrm{\Sigma }}_b\overline{K}$ molecule, with a preferred $1/2^{-}$ spin–parity assignment [[20]].

Over the past years, a unitarized coupled-channel scheme has been developed in Refs. [[21]–[27]] that implements heavy-quark spin symmetry (HQSS), which is a proper QCD symmetry that appears when the quark masses, such as that of the charm or bottom quark, become larger than the typical confinement scale. This approach is based on a consistent ${\mathrm{SU}(6)}_{\mathrm{lsf}}\times \mathrm{HQSS}$ extension of the Weinberg–Tomozawa (WT) $\pi N$ interaction, where "lsf" stands for light quark–spin–flavor symmetry, respectively. Within this framework, it has been identified a two-pole pattern for the ${\mathrm{\Lambda }}_c(2595)$ resonance[1] [[21], [23]], similar to the $\mathrm{\Lambda }(1405)$ [[29]–[31]]. The same scheme has also generated dynamically the ${\mathrm{\Lambda }}_b(5912)$ and ${\mathrm{\Lambda }}_b(5920)$ narrow resonances, discovered by LHCb [[32]], which turn out to be HQSS partners, naturally explaining their approximate mass degeneracy [[24]].

More recently, the work of Ref. [[23]] has been revisited in view of the newly discovered ${\mathrm{\Omega }}_c$ states, paying a special attention to the renormalization procedure used in the unitarized coupled-channel model and its impact on the ${\mathrm{\Omega }}_c$ sector. In Ref. [[33]] it was shown that some (probably at least three) of the ${\mathrm{\Omega }}_c$ states experimentally observed by LHCb would have $1/2^{-}$ or $3/2^{-}$ .

The discovery of the ${\mathrm{\Xi }}_c(2930)$ and ${\mathrm{\Xi }}_b(6227)$ has stimulated and motivated further research along this line. In the present work we follow a similar procedure as described in [[33]] and study the possible molecular interpretation of those states, revisiting the previous works on the ${\mathrm{\Xi }}_c$ [[23]] and ${\mathrm{\Xi }}_b$ [[24]] sectors. However, in the ${\mathrm{\Xi }}_c$ sector we do not restrict ourselves to the recent ${\mathrm{\Xi }}_c(2930)$ observation, but analyze all excited ${\mathrm{\Xi }}_c$ states found experimentally with masses up to 3 GeV [[1]]. The four excited ${\mathrm{\Xi }}_c$ states with masses below 3 GeV and the ${\mathrm{\Xi }}_b(6227)$ are collected in Table 1, showing the assigned spin–parity $J^P$ (when possible) as well as masses, widths, and decay channels. In this work we pay a special attention to the dependence on the renormalization scheme as well as to the flavor-symmetry content of the ${\mathrm{SU}(6)}_{\mathrm{lsf}}\times \mathrm{HQSS}$ model, as we determine the possible HQSS partners and siblings among the experimental states while predicting new ones. Thus, we follow the discussion of Ref. [[23]] on its spin–flavor symmetry breaking pattern. Flavor SU(4) is not a good symmetry in the limit of a heavy charm quark, for this reason, instead of the breaking pattern SU(8) $\supset$ SU(4), we consider the pattern SU(8) $\supset$ SU(6), since the light spin–flavor group [SU(6)] is decoupled from heavy-quark spin transformations. This allows us to implement HQSS in the analysis and to unambiguously identify the corresponding multiplets among the resonances generated dynamically. At the same time, we are also able to assign approximate heavy [SU(8)] and light [SU(6)] spin–flavor multiplet labels to the states.

This work is organized as follows. In Sect. 2 we present the ${\mathrm{SU}(6)}_{\mathrm{lsf}}\times \mathrm{HQSS}$ extension of the WT interaction, while in Sect. 3 we show our results for the ${\mathrm{\Xi }}_c$ and ${\mathrm{\Xi }}_b$ states,[2] respectively, and the possible experimental identification. Finally, in Sect. 4 we present our conclusions, emphasizing the possible classification of these experimental states according to the flavor-symmetry content of the scheme, while predicting new observations.

We consider the sector with charm $C=1$ , strangeness $S=-1$ and isospin $I=1/2$ quantum numbers, where the ${\mathrm{\Xi }}_c(2930)$ state has been observed by the Belle Collaboration [[11]]. Also, we examine the bottom $B=-1$ , strangeness $S=-1$ and isospin $I=1/2$ , where the ${\mathrm{\Xi }}_b(6227)$ has been found [[9]]. In order to do so, we revise the results in Ref. [[23]] for the ${\mathrm{\Xi }}_c$ states and in Ref. [[24]] for ${\mathrm{\Xi }}_b$ ones.

In the case of the $C=1$ , $S=-1$ and $I=1/2$ sector, the building-blocks are the pseudoscalar ( $D_s,D,K,\pi ,\eta ,\overline{K}$ ) and vector ( $D_s^{\ast },D^{\ast },K^{\ast },\rho ,\omega ,{\overline{K}}^{\ast },\varphi$ ) mesons, and the spin-1/2 ( $\mathrm{\Lambda }$ , $\mathrm{\Sigma }$ , $\mathrm{\Xi }$ , ${\mathrm{\Lambda }}_c$ , ${\mathrm{\Sigma }}_c$ , ${\mathrm{\Xi }}_c$ , ${\mathrm{\Xi }}_c^{\prime }$ , ${\mathrm{\Omega }}_c$ ), and spin-3/2 ( ${\mathrm{\Sigma }}_c^{\ast }$ , ${\mathrm{\Xi }}_c^{\ast }$ , ${\mathrm{\Omega }}_c^{\ast }$ ) charmed baryons [[21], [23]]. For the bottom sector $B=-1$ , $S=-1$ and $I=1/2$ , one can substitute the c quark by a b quark, and we have the pseudoscalar ( ${\overline{B}}_s,\overline{B},K,\pi ,\eta ,\overline{K}$ ) and vector ( ${\overline{B}}_s^{\ast },{\overline{B}}^{\ast },K^{\ast },\rho ,\omega ,{\overline{K}}^{\ast },\varphi$ ) mesons, and the spin-1/2 ( $\mathrm{\Lambda }$ , $\mathrm{\Sigma }$ , $\mathrm{\Xi }$ , ${\mathrm{\Lambda }}_b$ , ${\mathrm{\Sigma }}_b$ , ${\mathrm{\Xi }}_b$ , ${\mathrm{\Xi }}_b^{\prime }$ , ${\mathrm{\Omega }}_b$ ), and spin-3/2 ( ${\mathrm{\Sigma }}_b^{\ast }$ , ${\mathrm{\Xi }}_b^{\ast }$ , ${\mathrm{\Omega }}_b^{\ast }$ ) baryons [[24]]. All baryon–meson pairs with $(C=1/B=-1,S=-1,I=1/2)$ quantum numbers span the coupled-channel space for a given total angular momentum (J).

The S-wave tree level amplitudes between two baryon–meson channels are given by the ${\mathrm{SU}(6)}_{\mathrm{lsf}}\times \mathrm{HQSS}$ WT kernel,

1 $$\begin{array}{c}\hfill V_{\mathit{ij}}^J(s)=D_{\mathit{ij}}^J\frac{2\sqrt{s}-M_i-M_j}{4f_i{f}_j}\sqrt{\frac{E_i+M_i}{2M_i}}\sqrt{\frac{E_j+M_j}{2M_j}}.\end{array}$$

Graph

The $M_i$ and $m_i$ are the masses of the baryon and meson in the i channel, respectively, and $E_i$ is the center-of-mass energy of the baryon in the same channel,

2 $$\begin{array}{c}\hfill E_i=\frac{s-m_i^2+M_i^2}{2\sqrt{s}}.\end{array}$$

Graph

The projection onto definite J is done as explained in Refs. [[23], [34]]. In [[34]], only $J=1/2$ baryons and pseudoscalar mesons are considered. The inclusion of spin $J=3/2$ baryons and vector mesons was discussed in Appendix A of Ref. [[23]], using results derived in [[35]]. The hadron masses, meson decay constants, $f_i$ , and $D_{\mathit{ij}}^J$ matrices are taken from Refs. [[23]], where the underlying ${\mathrm{SU}(6)}_{\mathrm{lsf}}\times$ HQSS group structure of the interaction has been considered.

Starting from $V_{\mathit{ij}}^J$ , we solve the Bethe–Salpeter equation (BSE) in coupled channels,

3 $$\begin{array}{c}\hfill T^J(s)=\frac{1}{1-V^J(s)G^J(s)}{V}^J(s),\end{array}$$

Graph

where the $G^J(s)$ is a diagonal matrix that contains the different baryon–meson loop functions $G_i$ ,

4 $$\begin{array}{c}\hfill G_i(s)=i2M_i\int \frac{d^{4}q}{{(2\pi )}^4}\frac{1}{q^2-m_i^2+iϵ}\frac{1}{{(P-q)}^2-M_i^2+iϵ},\\ \hfill \end{array}$$

Graph

with P the total momentum of the system such that $P^2=s$ . We omit the index J from here on for simplicity. The bare loop function is logarithmically ultraviolet (UV) divergent and needs to be renormalized. This can be done by separating the divergent and finite parts of the loop function,

5 $$\begin{array}{c}\hfill G_i(s)={\overline{G}}_i(s)+G_i(s_{i+}),\end{array}$$

Graph

with the finite part of the loop function, ${\overline{G}}_i(s)$ , given in Refs. [[33]]. The divergent contribution of the loop function, $G_i(s_{i+})$ in Eq. (5) needs to be renormalized.

On the one hand, this can be done by one subtraction at certain scale ( $\sqrt{s}=\mu$ )

6 $$\begin{array}{c}\hfill G_i^{\mu }(s)={\overline{G}}_i(s)-{\overline{G}}_i({\mu }^2),\end{array}$$

Graph

where $\mu =\sqrt{m_{\mathit{th}}^2+M_{\mathit{th}}^2}$ , with $m_{\mathit{th}}$ and $M_{\mathit{th}}$ the masses of the meson and the baryon, respectively, that belong to the channel with the smallest threshold for a given (C, S, I) or (B, S, I) sectors. This common scale $\mu$ is chosen to be independent of total angular momentum J [[13], [36]], and it is the approach used in the previous works of Refs. [[23]].

On the other hand, as discussed in our recent paper [[33]], we could also use a sharp-cutoff regulator $\mathrm{\Lambda }$ in momentum space, so that

7 $$\begin{array}{c}\hfill G_i^{\mathrm{\Lambda }}(s)={\overline{G}}_i(s)+G_i^{\mathrm{\Lambda }}(s_{i+}),\end{array}$$

Graph

where $G_i^{\mathrm{\Lambda }}(s_{i+})$ is given in Refs. [[33], [37]].

Graph: Fig. 1 Evolution of the masses and widths of the dynamically generated Ξc states as we vary the renormalization scheme from using a subtraction constant to a common cutoff of Λ=1090 MeV. The cross symbolizes the position of the states in the subtraction constant scheme [or dimensional regularization (DR)] [[23]], while the triangle indicates the mass and width of the same states for the cutoff scheme

Note that if one uses channel-dependent cutoffs, the one-subtraction renormalization scheme is recovered by choosing ${\mathrm{\Lambda }}_i$ in each channel in such a way that

8 $$\begin{array}{c}\hfill G_i^{{\mathrm{\Lambda }}_i}(s_{i+})=-{\overline{G}}_i({\mu }^2).\end{array}$$

Graph

However, we employ a common UV cutoff for all baryon–meson loops within reasonable limits. In this manner, we avoid a fictitious reduction or enhancement of any baryon–meson loop by using an unappropriated value of the cutoff [[28], [38]], as well as we prevent an arbitrary variation of the subtraction constants, as we correlate all of them with a UV cutoff [[33]].

The poles of the T matrix describe the odd-parity dynamically-generated ${\mathrm{\Xi }}_c$ and ${\mathrm{\Xi }}_b$ states, which appear in the first and second Riemann sheets (FRS and SRS). Poles of the scattering amplitude on the FRS below threshold are bound states, whereas poles on the SRS below the real axis and above threshold are resonances. The mass and the width of the bound state/resonance can be found from the position of the pole on the complex energy plane. At the complex pole, the T-matrix is given by

9 $$\begin{array}{c}\hfill T_{\mathit{ij}}(s)=\frac{g_i{g}_j}{\sqrt{s}-\sqrt{s_R}}+\cdots ,\end{array}$$

Graph

where $\sqrt{s_R}=M_R-\mathrm{i}{\mathrm{\Gamma }}_R/2$ , with $M_R$ the mass and ${\mathrm{\Gamma }}_R$ the width of the state, and $g_i$ is the complex coupling of the state to the channel i, which is determined from Cauchy's theorem of residues. Thus, the dimensionless couplings $g_i$ are obtained by first assigning an arbitrary sign to one of them ( $g_1$ ), so

10 $$\begin{array}{c}\hfill g_1^2=\underset{\sqrt{s}\to \sqrt{s_R}}{lim}(\sqrt{s}-\sqrt{s_R})T_{11}(s).\end{array}$$

Graph

The other couplings are calculated as,

11 $$\begin{array}{c}\hfill g_j=g_1\underset{\sqrt{s}\to \sqrt{s_R}}{lim}\frac{T_{1j}(s)}{T_{11}(s)},\end{array}$$

Graph

and can be used to analyze the contribution of each baryon–meson channel to the generation of the state.

As mentioned in the Introduction, the first observation of the ${\mathrm{\Xi }}_c(2930)$ state was reported by the Belle Collaboration in Ref. [[11]]. This state was observed through its decay to ${\mathrm{\Lambda }}_c^{+}{K}^{-}$ with no assigned quantum numbers. Besides this recently discovered state, there are other three ${\mathrm{\Xi }}_c$ excited states with energies below 3 GeV [[1]]. As seen in Table 1, the $1/2^{-}$ ${\mathrm{\Xi }}_c(2790)$ state decays into ${\mathrm{\Xi }}_c^{\prime }\pi$ , whereas the $3/2^{-}$ ${\mathrm{\Xi }}_c(2815)$ decays into ${\mathrm{\Xi }}_c^{\prime }\pi$ and has also the decay chain ${\mathrm{\Xi }}_c^{\ast }\pi$ , followed by ${\mathrm{\Xi }}_c^{\ast }\to {\mathrm{\Xi }}_c\pi$ [[40]]. Also, a ${\mathrm{\Xi }}_c(2970)$ with unknown quantum numbers has been observed decaying into ${\mathrm{\Lambda }}_c^{+}\overline{K}\pi$ , ${\mathrm{\Sigma }}_c\overline{K}$ , ${\mathrm{\Xi }}_{c}2\pi$ , ${\mathrm{\Xi }}_c^{\prime }\pi$ and ${\mathrm{\Xi }}_c^{\ast }\pi$ .

Graph: Fig. 2 Evolution of the masses and widths of the dynamically generated Ξc states, as we vary the cutoff from Λ=1 GeV (triangles) to Λ=1.2 GeV (crosses). In a, b, the c1 (c8) state becomes virtual above (below) the ΛcK¯ (Ξc∗π) threshold. The squares and their associated errorbars show the masses and widths of the experimental Ξc(2790) and Ξc(2930) (a) and Ξc(2815) and Ξc(2970) (b) together with their experimental errors. The spin–parity of both Ξc(2930) and Ξc(2970) resonances are not experimentally determined [[1]], and we have displayed them here just for illustrative purposes

Masses and widths of the $c_1$ to $c_9$ states with $J=1/2$ or $J=3/2$ and odd parity in the $C=1$ , $S=-1$ and $I=1/2$ sector, together with the couplings (in modulus) to the dominant baryon–meson channels ( $g>1$ ) and the couplings to the decay channels reported experimentally for the ${\mathrm{\Xi }}_c$ states. All results have been obtained for $\mathrm{\Lambda }=1150$ MeV. We also indicate the possible experimental identification as well as the ${\mathrm{SU}(6)}_{\mathrm{lsf}}\times$ HQSS, SU(6) and SU(3) irreducible representations of these states. We use the notation ${\mathbf{R}}_{\mathbf{2}{\mathbf{J}}_{\mathbf{C}}+\mathbf{1},\mathbf{C}}$ , where $\mathbf{R}$ is the SU(6) irreducible representation (irrep) label (for which we use the dimension), $J_C$ is the spin carried by the quarks with charm (1/2 in all cases) and C the charm content (1 in all cases). In addition, we also use ${\mathbf{r}}_{\mathbf{2}\mathbf{J}+\mathbf{1}}$ , where $\mathbf{r}$ is the SU(3) irrep, with J the total angular momentum of the state (see Ref. [[23]] for details)

We start by revising the results Ref. [[23]] in the ${\mathrm{\Xi }}_c$ sector in order to understand whether the experimental states can be accommodated in our model. The widths of our ${\mathrm{\Xi }}_c$ states as a function of their masses in the $J=1/2$ and $J=3/2$ sectors are shown in the upper and lower plots of Fig. 1, respectively, together with different baryon–meson thresholds, to which they can couple. The dynamically generated states of Ref. [[23]] are displayed with a cross and the "DR" legend, as those have been obtained using one subtraction at certain scale or dimensional regularization. In what follows, we label the states as $c_1\cdots c_9$ , and they correspond to those given in Table V of Ref. [[23]]. They have either $J^P=1/2^{-}$ or $J^P=3/2^{-}$ and are ordered by their mass position. Hence, $c_1$ ( $c_9$ ) corresponds to the lightest (heaviest) state of mass 2699.4 MeV (2845.2 MeV), among those quoted in the mentioned table, where their SU(6) and SU(3) quantum numbers are also given. We observe that the masses of our ${\mathrm{\Xi }}_c$ states using one subtraction constant (DR) are below or close to the experimental ${\mathrm{\Xi }}_c(2790)$ or ${\mathrm{\Xi }}_c(2815)$ states, while being far below in mass with respect to ${\mathrm{\Xi }}_c(2930)$ or ${\mathrm{\Xi }}_c(2970)$ .

We then study the effect on masses and widths of the renormalization procedure so as to determine whether any our ${\mathrm{\Xi }}_c$ can be identified with a experimental state while assessing the dependence on the renormalization scheme, which might be significant (as shown in Ref. [[28]]). We proceed as described in Ref. [[33]] for the ${\mathrm{\Omega }}_c$ states, where we explore a different renormalization procedure, the cutoff scheme. In order to do so, we first need to determine how the masses and widths of our dynamically generated states change as we adiabatically move from the one subtraction renormalization scheme to the cutoff one. Thus, we change the loop functions by

12 $$\begin{array}{c}\hfill G_i^{\mathrm{\Lambda }}(s)={\overline{G}}_i(s)-(1-x){\overline{G}}_i({\mu }^2)+x{\overline{G}}_i^{\mathrm{\Lambda }}(s),\end{array}$$

Graph

where we slowly evolve x from 0 to 1 while following the evolution of the states, as seen in Fig. 1. The $c_1$ to $c_9$ states for a $\mathrm{\Lambda }=1090$ MeV are shown with a triangle. We find that most of these states move to higher energies, except for $c_2$ , $c_5$ and $c_8$ , whereas getting closer to the experimental values. Note that fot this cutoff, the $J^P=1/2^{-}$ $c_1$ state become virtual above the ${\mathrm{\Lambda }}_c\overline{K}$ threshold.

Once we have identified our ${\mathrm{\Xi }}_c$ states in the cutoff scheme, we can assess the dependence of our results on this regulator, as well as their possible experimental identification. In Fig. 2 we show the evolution of the $c_1$ to $c_9$ states as we vary the cutoff from 1 GeV (triangles) to 1.2 GeV (crosses), and we also display different two-body thresholds. Moreover, in Table 2 we show masses and widths of the $c_1$ to $c_9$ states with $J=1/2$ or $J=3/2$ , together with the couplings to the dominant baryon–meson channels ( $g>1$ ) and the couplings to the decay channels reported experimentally for the ${\mathrm{\Xi }}_c$ states. All these results are obtained for $\mathrm{\Lambda }=1150$ MeV. In this table we also indicate the possible experimental identification as well as the ${\mathrm{SU}(6)}_{\mathrm{lsf}}\times$ HQSS, SU(6) and SU(3) irreducible representations (irreps), to which the $c_1$ to $c_9$ states belong (see Ref. [[23]] for group-structure details).

As we evolve the cutoff value from $\mathrm{\Lambda }=1000$ MeV to $\mathrm{\Lambda }=1200$ MeV, that is, from the right to left in Fig. 2, we observe that some of our dynamically generated resonances can be identified with the experimental states attending to the complex energy position.

In the $J^P=1/2^{-}$ sector, attending to the position in the complex plane, we observe that the ${\mathrm{\Xi }}_c(2790)$ could be one of the $c_1$ , $c_3$ , $c_6$ or even the $c_5$ states. The identification with the ${\mathrm{\Xi }}_c(2790)$ is possible because these states couple to ${\mathrm{\Xi }}_c^{\prime }\pi$ , although this baryon–meson channel is not the dominant one for their dynamically generation, as seen in Table 2 for a $\mathrm{\Lambda }=1150$ MeV, except for $c_5$ . Indeed, this latter feature of $c_5$ disfavors its identification with the ${\mathrm{\Xi }}_c(2790)$ . This is because it would become too broad ( $\mathrm{\Gamma }\ge 70$ MeV) for UV cutoffs of around 1 GeV, that would lead the $c_5$ resonance to have masses closer to the experimental one, as seen in Fig. 2. In addition in the DR scheme, the mass of the $c_5$ state is close to 2790 MeV, but its width is approximately of 84 MeV [[23]] (see also Fig. 1), while experimentally ${\mathrm{\Gamma }}_{{\mathrm{\Xi }}_c(2790)}\sim 10$ MeV.

Looking at the behavior of the $c_1$ , $c_3$ , $c_6$ poles with the UV cutoff in Fig. 2, it seems natural to assign the ${\mathrm{\Xi }}_c(2790)$ to the $c_1$ pole. This state has a width of the order of 10 MeV for UV cutoffs in the region of 1.2 GeV, where it is located below the ${\mathrm{\Lambda }}_c\overline{K}$ threshold. At the same time, the state has large ${\mathrm{\Lambda }}_c\overline{K}$ and small ${\mathrm{\Xi }}_c^{\prime }\pi$ couplings (see Table 2), respectively, which explains its small experimental width despite being placed well above the latter threshold, and it is natural to think that the ${\mathrm{\Lambda }}_c\overline{K}$ channel should play an important role in the dynamics of the ${\mathrm{\Xi }}_c(2790)$ given its proximity to that threshold. Note that the light degrees of freedom (ldof) in the inner structure of the $c_1$ are predominantly coupled to $j_{\mathit{ldof}}^{\pi }=0^{-}$ spin–parity quantum numbers [[28]]. Thus with this identification, this first odd parity excited ${\mathrm{\Xi }}_c$ state would not have a dominant configuration consisting of a spinless light diquark and a unit of angular momentum between it and the heavy quark, as argued for instance in the Belle Collaboration paper [[40]]. This is to say, the ${\mathrm{\Xi }}_c(2790)$ will not be a constituent quark model $\lambda$ -mode excited state [[28], [41]] with $j_{\mathit{ldof}}^{\pi }=1^{-}$ and hence it will not form part of any HQSS doublet, thus making the assignment to $c_3$ unlikely. Actually, if the spin–parity quantum numbers for the $\mathit{ldof}$ in the ${\mathrm{\Xi }}_c(2790)$ were predominantly $1^{-}$ , one would expect a larger width for this resonance, since its decay to the open channel ${\mathrm{\Xi }}_c^{\prime }\pi$ is HQSS allowed. This is precisely the situation for the $c_3$ that is broader than the experimental state. In summary, we conclude a large molecular ${\mathrm{\Lambda }}_c\overline{K}$ component for the ${\mathrm{\Xi }}_c(2790)$ that will have then a dominant $j_{\mathit{ldof}}^{\pi }=0^{-}$ configuration.

Our present ${\mathrm{\Xi }}_c(2790)$ identification with the $c_1$ pole differs from the previous assignments in Ref. [[23]], where the ${\mathrm{\Xi }}_c$ states were obtained using the one subtraction renormalization scheme. In this previous work, the $c_7$ state was assigned to ${\mathrm{\Xi }}_c(2790)$ due to its closeness in energy and the sizable ${\mathrm{\Xi }}_c^{\prime }\pi$ coupling within the DR scheme.

With regards to the recently discovered ${\mathrm{\Xi }}_c(2930)$ , if we assume that this state has $J^P=1/2^{-}$ , we could identify it either with our $c_6$ or $c_7$ states, as they both couple to the ${\mathrm{\Lambda }}_c\overline{K}$ channel, although not dominantly as seen in Table 2 for a $\mathrm{\Lambda }=1150$ MeV. The assignment to the $c_6$ pole is, however, disfavored because of the mass difference between this state and the experimental ${\mathrm{\Xi }}_c(2930)$ . As for $c_7$ , the small ${\mathrm{\Lambda }}_c\overline{K}$ coupling of this state makes also somehow doubtful its identification with the ${\mathrm{\Xi }}_c(2930)$ . In the case of our $c_2$ and $c_5$ states, we should mention that we do not have any clear experimental candidate at this point for the $c_5$ dynamically generated $J=1/2$ state, whereas the $c_2$ state becomes broad and appears below 2650 MeV, thus not allowing for any reasonable experimental assignment.

For $J^P=3/2^{-}$ , the analysis of the evolution of the different states in Fig. 2 allows for the identification of the experimental ${\mathrm{\Xi }}_c(2815)$ with $c_4$ or $c_8$ . These states couple to ${\mathrm{\Xi }}_c^{\ast }\pi$ in S-wave, although for $c_4$ , couplings to other baryon–meson states ( ${\mathrm{\Sigma }}_c^{\ast }\overline{K}$ , ${\mathrm{\Lambda }}_c{\overline{K}}^{\ast }$ or ${\mathrm{\Sigma }}_c^{\ast }{\overline{K}}^{\ast }$ ) are larger as seen in Table 2. The experimental ${\mathrm{\Xi }}_c(2815)$ is quite narrow, ${\mathrm{\Gamma }}_{{\mathrm{\Xi }}_c(2815)}\sim 2-3$ MeV, despite the ${\mathrm{\Xi }}_c^{\ast }\pi$ threshold being around 30 MeV below its mass. This hints to a subdominant ${\mathrm{\Xi }}_c^{\ast }\pi$ molecular component in the inner structure of this resonance. Moreover, looking at the dependence of the $J^P=3/2^{-}$ pole masses and widths with the UV cutoff displayed in Fig. 2, it seems reasonable to assign the $c_4$ state to the ${\mathrm{\Xi }}_c(2815)$ resonance.

As for ${\mathrm{\Xi }}_c(2970)$ , assuming that it has $J=3/2^{-}$ , we could identify it with the $c_9$ state for values of the cutoff around $\mathrm{\Lambda }\simeq 1.1$ GeV. In this case, we have to take into account that this state couples to ${\mathrm{\Lambda }}_c{\overline{K}}^{\ast }$ and ${\mathrm{\Sigma }}_c^{\ast }\overline{K}$ , and ${\mathrm{\Xi }}_c^{\ast }\pi$ (though not dominantly), and those baryon–meson channels can decay into ${\mathrm{\Lambda }}_c\overline{K}\pi$ and ${\mathrm{\Xi }}_c\pi \pi$ , respectively. Nevertheless, the predicted width would be significantly smaller than the range of 20-30 MeV quoted in the PDG [[1]] and shown in Table 1. Compared to the results of Ref. [[23]], the ${\mathrm{\Xi }}_c(2815)$ was identified there with $c_9$ , assuming that ${\mathrm{\Xi }}_c(2790)$ and ${\mathrm{\Xi }}_c(2815)$ were the $c_7$ and $c_9$ HQSS partners.

In fact, we observe several HQSS partners among our states as well as possible siblings within the same SU(3) representation. The ${\mathrm{\Xi }}_c(2790)$ resonance belongs to an $J=1/2$ SU(3) antitriplet irrep, and it would be the ${\mathrm{SU}(6)}_{\mathrm{lsf}}\times$ HQSS (see Table 2) partner of a narrow ${\mathrm{\Lambda }}_c^{\ast }$ state discussed in [[21], [23], [28]]. This latter state has large (small) ND and $N{D}^{\ast }$ ( ${\mathrm{\Sigma }}_c\pi$ ) couplings, and depending on the renormalization scheme (one-subtraction or UV cutoff), it is part of a double pole pattern for the ${\mathrm{\Lambda }}_c(2595)$ , similar to that found for the $\mathrm{\Lambda }(1405)$ within unitarized chiral models [[29], [42]–[47]] (see related review in [[1]]), or it is located in the region of 2.8 GeV close to the ND threshold [[28]].

Masses and widths of the $b_1$ to $b_9$ states with $J=1/2$ or $J=3/2$ in the $B=-1$ , $S=-1$ and $I=1/2$ sector, together with the couplings to the dominant baryon–meson channels and the couplings to the experimental decay channels of the ${\mathrm{\Xi }}_b(6227)$ , using one-subtraction renormalization, as in Table IV of Ref. [[24]]. We also indicate the ${\mathrm{SU}(6)}_{\mathrm{lsf}}\times$ HQSS, SU(6) and SU(3) irreducible representations of these states, as in Table 2. States with $\ast$ are virtual states. Note that the $b_8^{\ast }$ lies in the real axis, but in a sheet that is not connected to the physical sheet, thus we are not showing the couplings indicating "R.S (real sheet) not connected"

On the other hand, the $c_3$ pole belonging to ( $\mathbf{168},\mathbf{21},{\mathbf{6}}_{\mathbf{2}}$ ) representation and the $c_4$ of the ( $\mathbf{168},\mathbf{21},{\mathbf{6}}_{\mathbf{4}}$ ) form a $(j_{\mathit{ldof}}^{\pi }=1^{-})$ -HQSS doublet. As mentioned earlier, the $c_4$ can be identified with the ${\mathrm{\Xi }}_c(2815)$ , but we note that the ${\mathrm{\Xi }}_c(2815)$ is not the sibling of the ${\mathrm{\Lambda }}_c(2625)$ because of the different coupling strengths to ${\mathrm{\Xi }}_c^{\ast }\pi$ and ${\mathrm{\Sigma }}_c^{\ast }\pi$ , respectively. Whereas ${\mathrm{\Xi }}_c(2815)$ weakly couples to ${\mathrm{\Xi }}_c^{\ast }\pi$ , the ${\mathrm{\Lambda }}_c(2625)$ strongly does to ${\mathrm{\Sigma }}_c^{\ast }\pi$ . However, this latter state is narrow because the ${\mathrm{\Sigma }}_c^{\ast }\pi$ channel is closed (located around 30 MeV above the mass of the resonance). Indeed, recently it has been argued that the $\mathrm{\Lambda }(2625)$ is probably a constituent three quark state [[28], [41]].

As for the $J=1/2$ $c_5$ and the $J=3/2$ $c_8$ states, those form part of a SU(6) $\mathbf{15}$ -plet, belonging to the ${\mathrm{SU}(6)}_{\mathrm{lsf}}\times$ HQSS ( $\mathbf{168},\mathbf{15},{\mathbf{3}}_{\mathbf{2}}$ ) and ( $\mathbf{168},\mathbf{15},{\mathbf{3}}_{\mathbf{4}}$ ) irreps [[23]]. They form a HQSS doublet with $j_{\mathit{dof}}^{\pi }=1^{-}$ and hence have large couplings to ${\mathrm{\Xi }}_c^{\prime }\pi$ and ${\mathrm{\Xi }}_c^{\ast }\pi$ , respectively. Indeed, as a good approximation, they are dynamically generated by the charmed baryon–Goldstone boson interactions. These moderately broad states are in the SU(3) ${}_{2J+1}$ ${\mathbf{3}}_{\mathbf{2}}^{\ast }$ and ${\mathbf{3}}_{\mathbf{4}}^{\ast }$ irreps, which should be completed by one $J=1/2$ and one $J=3/2$ ${\mathrm{\Lambda }}_c$ resonances stemming from the ${\mathrm{\Sigma }}_c\pi$ and ${\mathrm{\Sigma }}_c^{\ast }\pi$ chiral interactions [[19], [28]], neglecting higher energy channels. The $J=3/2$ sibling is, however, not the ${\mathrm{\Lambda }}_c(2625)$ . As mentioned before, the ${\mathrm{\Lambda }}_c(2625)$ is probably a quark model ( $\lambda$ -mode excitation) state [[28], [41]]. Another resonance with mass and width of around 2.7 GeV and 60 MeV [[19], [28]], that has not been discovered yet, would then be the SU(3) sibling of the $c_8$ state.

The features of the $J=1/2$ counterpart of $c_5$ in the ${\mathrm{\Lambda }}_c$ sector are much more uncertain and depend on both the employed renormalization scheme and on the interplay between quark-model and baryon–meson degrees of freedom [[28]]. Thus, for instance neglecting the latter, it would appear around 2.6 GeV with a large width of 60-80 MeV because its sizable coupling to the ${\mathrm{\Sigma }}_c\pi$ pair. Within the UV cutoff RS, this state can be easily moved below the ${\mathrm{\Sigma }}_c\pi$ threshold and be identified with the narrow ${\mathrm{\Lambda }}_c(2595)$ [[19]]. In the DR scheme advocated in Ref. [[23]], this broad state, together with the $j_{\mathit{ldof}}^{\pi }=0^{-}$ narrow state mentioned above in the discussion of the ${\mathrm{\Xi }}_c(2790)$ , gives rise to a double pole structure for the ${\mathrm{\Lambda }}_c(2595)$ .

Within the UV cutoff renormalization scheme examined here, the $(c_7,c_9)$ HQSS-doublet might correspond to the experimental ${\mathrm{\Xi }}_c(2930)$ and ${\mathrm{\Xi }}_c(2970)$ states. The $c_7$ state, that we have tentatively assigned to the ${\mathrm{\Xi }}_c(2930)$ , exhibits (Table 2) moderate couplings to ${\mathrm{\Xi }}_c\pi$ and ${\mathrm{\Lambda }}_c\overline{K}$ , small ones to ${\mathrm{\Xi }}_c^{\prime }\pi$ and ${\mathrm{\Sigma }}_c\overline{K}$ , and finally large couplings to $\mathrm{\Lambda }{D}^{(\ast )}$ , $\mathrm{\Sigma }{D}^{(\ast )}$ and ${\mathrm{\Sigma }}^{\ast }{D}^{\ast }$ . It belongs to a SU(3) sextet, where there is also a ${\mathrm{\Omega }}_c$ state. The latter corresponds to the one labeled as d in our previous study of the ${\mathrm{\Omega }}_c$ odd-parity resonances [[33]], where it was tentatively assigned either to the ${\mathrm{\Omega }}_c(3090)$ or the ${\mathrm{\Omega }}_c(3119)$ observed by the LHCb Collaboration in the ${\mathrm{\Xi }}_c\overline{K}$ mode [[8]]. This is in fact consistent with what one might expect from its $c_7$ -sibling couplings. Assuming the equal spacing rule we could predict the possible existence of a $J=1/2^{-}$ ${\mathrm{\Sigma }}_c$ state around 2800 MeV that will complete the sextet. The ${\mathrm{\Sigma }}_c(2800)$ clearly fits into this picture since it is observed in the ${\mathrm{\Lambda }}_c\pi$ channel [[1]].

Graph: Fig. 3 Evolution of the masses and widths of the dynamically generated Ξb states, as we vary the cutoff from Λ=1000 MeV (triangles) to Λ=1400 MeV (crosses), with J=1/2 (upper panel) and J=3/2 (lower panel). The square and its bars represent the position of the Ξb(6227) resonance, and its errors in mass and width, respectively. We show the experimental result for both values of J due to its unknown quantum numbers. In b, the last five thresholds (not labelled in the figure because they are too close to each other) are: Ξb∗η(6492.45MeV), ΣB∗(6518.35MeV), ΩbK(6564.68MeV), Ξbρ(6565.04MeV) and Ξbω (6572.12MeV)

Recently there has been an analysis of the ${\mathrm{\Xi }}_c$ sector within a baryon–meson molecular model based on local hidden gauge that implements the interaction between the $1/2^{+}$ and $3/2^{+}$ ground-state baryons with $0^{-}$ and $1^{-}$ mesons [[18]]. The authors have found that five of their dynamically generated ${\mathrm{\Xi }}_c$ states can be identified with the experimental ${\mathrm{\Xi }}_c(2790)$ , ${\mathrm{\Xi }}_c(2930)$ , ${\mathrm{\Xi }}_c(2970)$ , ${\mathrm{\Xi }}_c(3055)$ and ${\mathrm{\Xi }}_c(3080)$ . Whereas the ${\mathrm{\Xi }}_c(2790)$ would be a $1/2^{-}$ state, the ${\mathrm{\Xi }}_c(2930)$ , ${\mathrm{\Xi }}_c(2970)$ , ${\mathrm{\Xi }}_c(3055)$ and ${\mathrm{\Xi }}_c(3080)$ could be either $1/2^{-}$ or $3/2^{-}$ ones. Compared to this approach, our model identifies the experimental ${\mathrm{\Xi }}_c(2790)$ and ${\mathrm{\Xi }}_c(2930)$ as $1/2^{-}$ states, and the ${\mathrm{\Xi }}_c(2815)$ and ${\mathrm{\Xi }}_c(2970)$ as $3/2^{-}$ . The different assignment is mainly due the distinct renormalization scheme used in the two approaches as well as the fact the interactions involving D and $D^{\ast }$ and light vector mesons with baryons are not completely fixed by HQSS or chiral symmetries, thus allowing for different assumptions.

With regards to the bottom sector, the ${\mathrm{\Xi }}_b(6227)$ resonance has been recently measured by the LHCb experiment [[9]], with ${\mathrm{\Gamma }}_{{\mathrm{\Xi }}_b(6227)}\sim 18$ MeV. Its quantum numbers, though, remain unknown, whereas the observed decay channels are ${\mathrm{\Lambda }}_b^0{K}^{-}$ and ${\mathrm{\Xi }}_b^0{\pi }^{-}$ (see Table 1).

We start again by revising the previous results of Ref. [[24]] with $B=-1$ , $S=-1$ , $I=1/2$ ( ${\mathrm{\Xi }}_b$ sector). Masses and widths of the dynamically generated states within our model using the DR scheme, together with their irreps, spins and couplings to the dominant baryon–meson channels as well those for the experimental decay channels of ${\mathrm{\Xi }}_b(6227)$ are shown in Table 3. We obtain nine states, which are the bottom counterparts of the ${\mathrm{\Xi }}_c$ ones discussed in the previous subsection. Compared to Ref. [[24]], we report here five more poles, since in that reference only SU(3) flavor partners of ${\mathrm{\Lambda }}_b$ states were searched (members of antitriplet irreps). Also, two of them, the state at 6035 MeV with $J=1/2$ and the one at 6043 MeV with $J=3/2$ were wrongly assigned in Ref. [[24]] to the SU(6) 15 representation. Instead, their should belong to the SU(6) 21 representation, as seen in Table 3. Moreover, there is a state at 6073 MeV in Table IV in Ref. [[24]] that does not appear in our present calculation. The differences between of them are due to the difficulty in determining the number of states and their representations as we break the ${\mathrm{SU}(6)}_{\mathrm{lsf}}\times$ HQSS symmetry to SU(3) in the bottom sector, as almost all states have zero width and states with widths closer to zero are more difficult to follow in the complex energy plane.

As Table 2, but for the ${\mathrm{\Xi }}_b$ sector ( $\mathrm{\Lambda }=1150$ MeV)