Ξc and Ξb excited states within a SU(6)lsf×HQSS model

SU(6)lsf×HQSS model J. Nieves, R. Pavao, and L. Tolos 3, 4 Instituto de F́ısica Corpuscular (centro mixto CSIC-UV), Institutos de Investigación de Paterna, Aptdo. 22085, 46071, Valencia, Spain Institut für Theoretische Physik, University of Frankfurt, Max-von-Laue-Str. 1, 60438 Frankfurt am Main, Germany Frankfurt Institute for Advanced Studies, University of Frankfurt, Ruth-Moufang-Str. 1, 60438 Frankfurt am Main, Germany Institute of Space Sciences (CSIC-IEEC), Campus Universitat Autònoma de Barcelona, Carrer de Can Magrans, s/n, 08193 Cerdanyola del Vallès, Spain (Dated: November 15, 2019)


I. INTRODUCTION
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][3][4][5][6][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 Ω c excited states in the Ξ + c K − spectrum in pp collisions [8], and the excited Ξ b (6227) state in Λ 0 b K − and Ξ 0 b π − invariant mass spectra also in pp collisions [9]. Moreover, the Belle Collaboration has confirmed the observation of four of the excited Ω c states [10], and detected the Ξ c (2930) state in its decay to Λ + c K − in [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 Ω c states, the molecular models of Refs. [12,13] have been reanalyzed in view of the new discoveries. While in Ref. [14] two Ω 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 − Ω c states and one J P = 3/2 − Ω c were determined within an extended local hidden gauge approach, the first two in good agreement with [14]. Other theo-retical works also examined the Ω 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Ξ bound states, whereas in Ref. [17] a loosely bound molecule of mass 3140 MeV was determined.
With regards to Ξ c , the theoretical analysis based on the local hidden gauge formalism has shown that not only the Ξ c (2930) can have a molecular interpretation, but also other Ξ c states around 3 GeV reported in the PDG [1]. In particular, the Ξ c (2790) would be a J P = 1/2 − molecular state, whereas Ξ c (2930), Ξ c (2970), Ξ c (3055) and Ξ 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 Ξ 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 Ξ b state has been also studied within a unitarized model that uses the leading-order chiral Lagrangian in Refs. [19,20], identifying the Ξ b (6227) state as a S-wave Σ bK molecule, with a preferred 1/2 − spin-parity assignment [20].
More recently, the work of Ref. [23] has been revisited in view of the newly discovered Ω c states, paying a special attention to the renormalization procedure used in the unitarized coupled-channel model and its impact on the Ω c sector. In Ref. [33] it was shown that some (probably at least three) of the Ω c states experimentally observed by LHCb would have 1/2 − or 3/2 − .
The discovery of the Ξ c (2930) and Ξ 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 Ξ c [23] and Ξ b [24] sectors. However, in the Ξ c sector we do not restrict ourselves to the recent Ξ c (2930) observation, but analyze all excited Ξ c states found experimentally with masses up to 3 GeV [1]. The four excited Ξ c states with masses below 3 GeV and the Ξ b (6227) are collected in Table I, 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 SU(6) lsf × 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) ⊃ SU(4), we consider the pattern SU(8) ⊃ 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 Section II we present the SU(6) lsf × HQSS extension of the WT interaction, while in Section III we show our results for the Ξ c and Ξ b states 2 , respectively, and the possible experimental identification. Finally, in Section IV 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.  I: Excited Ξ c states below 3 GeV and the excited Ξ b state found experimentally [1]. We show the assigned J P (when possible), the mass M and width Γ, as well as the decay channels.
We consider the sector with charm C = 1, strangeness S = −1 and isospin I = 1/2 quantum numbers, where the Ξ 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 Ξ b (6227) has been found [9]. In order to do so, we revise the results in Ref. [23] for the Ξ c states and in Ref. [24] 1 The details of this double pole structure, generated by the Σcπ, N D and N D * coupled-channels dynamics, depend strongly on the adopted renormalization scheme, which could considerably enhance the role played by the two latter channels around the resonance energy. This is discuss in detail in Ref. [28]. 2 From now on we refer to excited Ξc and Ξ b independently of 1/2 − or 3/2 − spin-parity assignment for Ξ 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, π, η,K) and vector (D * s , D * , K * , ρ, ω,K * , φ) mesons, and the spin-1/2 (Λ, Σ, Ξ, Λ c , Σ c , Ξ c , Ξ c , Ω c ), and spin-3/2 (Σ * c , Ξ * c , Ω * c ) 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 (B s ,B, K, π, η,K) and vector (B * s ,B * , K * , ρ, ω,K * , φ) mesons, and the spin-1/2 (Λ, Σ, Ξ, [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 SU(6) lsf × HQSS WT kernel, (1) 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, The hadron masses, meson decay constants, f i , and D J ij matrices are taken from Ref. [23,24], where the underlying SU(6) lsf × HQSS group structure of the interaction has been considered. Starting from V J ij , we solve the Bethe-Salpeter equation (BSE) in coupled channels, where the G J (s) is a diagonal matrix that contains the different baryon-meson loop functions G i , 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, with the finite part of the loop function, G i (s), given in Refs. [34] and [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 ( where µ = m 2 th + M 2 th , with m th and M 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 µ is chosen to be independent of total angular momentum J [13,35], and it is the approach used in the previous works of Refs. [23,24].
On the other hand, as discussed in our recent paper [33], we could also use a sharp-cutoff regulator Λ in momentum space, so that where G Λ i (s i+ ) is given in Refs. [33,36]. Note that if one uses channel-dependent cutoffs, the one-subtraction renormalization scheme is recovered by choosing Λ i in each channel in such a way that However, we employ a common UV cutoff for all baryonmeson 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,37,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 Ξ c and Ξ 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. Close to the pole, the T -matrix behaves as where √ s R = M R − i Γ R /2, with M R the mass and Γ R the width of the state, and g i is the complex coupling of the state to the channel i. The dimensionless couplings g i are obtained by first assigning an arbitrary sign to one of them (g 1 ), so The other couplings are calculated as, and can be used to analyze the contribution of each baryon-meson channel to the generation of the state.

A.
Ξc excited states 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.
As mentioned in the Introduction, the first observation of the Ξ c (2930) state was reported by the Belle Collaboration in Ref. [11]. This state was observed through its decay to Λ + c K − with no assigned quantum numbers. Besides this recently discovered state, there are other three Ξ c excited states with energies below 3 GeV [1]. As seen in Table I, the 1/2 − Ξ c (2790) state decays into Ξ c π, whereas the 3/2 − Ξ c (2815) decays into Ξ c π and has also the decay chain Ξ * c π, followed by Ξ * c → Ξ c π [39]. Also, a Ξ c (2970) with unknown quantum numbers has been ob-served decaying into Λ + cK π, Σ cK , Ξ c 2π, Ξ c π and Ξ * c π. We start by revising the results Ref. [23] in the Ξ c sector in order to understand whether the experimental states can be accommodated in our model. The widths of our Ξ 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 . . . 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 Ξ c states using one subtraction constant (DR) are below or close to the experimental Ξ c (2790) or Ξ c (2815) states, while being far below in mass with respect to Ξ c (2930) or Ξ c (2970).
We then study the effect on masses and widths of the renormalization procedure so as to determine whether any our Ξ 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 Ω 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 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 Λ = 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 Λ cK threshold. Once we have identified our Ξ 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 II 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 Ξ c states. All these results are obtained for Λ = 1150 MeV. In this table we also indicate the SU(6) lsf ×HQSS, SU(6) and SU(3) ir-reducible 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 Λ = 1000 MeV to Λ = 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 Ξ c (2790) could be one of the c 1 , c 3 , c 6 or even the c 5 states. The identification with the Ξ c (2790) is possible because these states couple to Ξ c π, although this baryon-meson channel is not the dominant one for their dynamically generation, as seen in Table II for a Λ = 1150 MeV, except for c 5 . Indeed, this latter feature of c 5 disfavors its identification with the Ξ c (2790). This is because it would become too broad (Γ ≥ 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 Γ Ξc(2790) ∼ 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 Ξ 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 Λ cK threshold. At the same time, the state has large Λ cK and small Ξ c π couplings (see Table II), respectively, which explains its small experimental width despite being placed well above the latter threshold, and it is natural to think that the Λ cK channel should play an important role in the dynamics of the Ξ 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 π ldof = 0 − spin-parity quantum numbers [28]. Thus with this identification, this first odd parity excited Ξ 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 [39]. This is to say, the Ξ c (2790) will not be a constituent quark model λ−mode excited state [28,40] with j π ldof = 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 ldof in the Ξ c (2790) were predominantly 1 − , one would expect a larger width for this resonance, since its decay to the open channel Ξ c π 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 Λ cK component for the Ξ c (2790) that will have then a dominant j π ldof = 0 − configuration. Our present Ξ c (2790) identification with the c 1 pole differs from the previous assignments in Ref. [23], where the Ξ c states were obtained using the one subtraction renormalization scheme. In this previous work, the c 7 state was assigned to Ξ c (2790) due to its closeness in energy and the sizable Ξ c π coupling within the DR scheme. GeV (triangles) to Λ = 1.2 GeV (crosses). In Fig. 2a (Fig. 2b), the c 1 (c 8 ) state becomes virtual above (below) the Λ c K (Ξ * c π) threshold. The squares and their associated errorbars show the masses and widths of the experimental Ξ c (2790) and Ξ c (2930) (Fig. 2a) and Ξ c (2815) and Ξ c (2970) (Fig. 2b) 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.  II: 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 Ξ c states. All results have been obtained for Λ = 1150 MeV. We also indicate the SU(6) lsf × HQSS, SU(6) and SU(3) irreducible representations of these states. We use the notation R 2J C +1,C , where 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 r 2J+1 , where r is the SU(3) irrep, with J the total angular momentum of the state (see Ref. [23] for details).
With regards to the recently discovered Ξ 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 Λ cK channel, although not dominantly as seen in Table II for a Λ = 1150 MeV. The assignment to the c 6 pole is, however, disfavored because of the mass difference between this state and the experimental Ξ c (2930). As for c 7 , the small Λ cK coupling of this state makes also somehow doubtful its identification with the Ξ 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 Ξ c (2815) with c 4 or c 8 . These states couple to Ξ * c π in S−wave, although for c 4 , couplings to other baryon-meson states (Σ * cK , Λ cK * or Σ * cK * ) are larger as seen in Table II. The experimental Ξ c (2815) is quite narrow, Γ Ξc(2815) ∼ 2 − 3 MeV, despite the Ξ * c π threshold being around 30 MeV below its mass. This hints to a subdominant Ξ * c π 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 Ξ c (2815) resonance.
As for Ξ c (2970), assuming that it has J = 3/2 − , we could identify it with the c 9 state for values of the cutoff around Λ 1.1 GeV. In this case, we have to take into account that this state couples to Λ cK * and Σ * cK , and Ξ * c π (though not dominantly), and those baryon-meson channels can decay into Λ cK π and Ξ c ππ, 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 I. Compared to the results of Ref. [23], the Ξ c (2815) was identified there with c 9 , assuming that Ξ c (2790) and Ξ 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 Ξ c (2790) resonance belongs to an J = 1/2 SU(3) antitriplet irrep, and it would be the SU(6) lsf ×HQSS (see Table II) partner of a narrow Λ * c state discussed in [21,23,28]. This latter state has large (small) N D and N D * (Σ c π) couplings, and depending on the renormalization scheme (one-subtraction or UV cutoff), it is part of a double pole pattern for the Λ c (2595), similar to that found for the Λ(1405) within unitarized chiral models [29,30,[41][42][43][44][45][46] (see related review in [1]), or it is located in the region of 2.8 GeV close to the N D threshold [28].
On the other hand, the c 3 pole belonging to (168, 21, 6 2 ) representation and the c 4 of the (168, 21, 6 4 ) form a (j π ldof = 1 − )−HQSS doublet. As mentioned earlier, the c 4 can be identified with the Ξ c (2815), but we note that the Ξ c (2815) is not the sibling of the Λ c (2625) because of the different coupling strengths to Ξ * c π and Σ * c π, respectively. Whereas Ξ c (2815) weakly couples to Ξ * c π, the Λ c (2625) strongly does to Σ * c π. However, this latter state is narrow because the Σ * c π channel is closed (located around 30 MeV above the mass of the resonance). Indeed, recently it has been argued that the Λ(2625) is probably a constituent three quark state [28,40].
The features of the J = 1/2 counterpart of c 5 in the Λ 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 Σ c π pair. Within the UV cutoff RS, this state can be easily moved below the Σ c π threshold and be identified with the narrow Λ c (2595) [19]. In the DR scheme advocated in Ref. [23], this broad state, together with the j π ldof = 0 − narrow state mentioned above in the discussion of the Ξ c (2790), gives rise to a double pole structure for the Λ c (2595).
Within the UV cutoff renormalization scheme examined here, the (c 7 , c 9 ) HQSS-doublet might correspond to the experimental Ξ c (2930) and Ξ c (2970) states. The c 7 state, that we have tentatively assigned to the Ξ c (2930), exhibits (Table II) moderate couplings to Ξ c π and Λ cK , small ones to Ξ c π and Σ cK , and finally large couplings to ΛD ( * ) , ΣD ( * ) and Σ * D * . It belongs to a SU(3) sextet, where there is also a Ω c state. The latter corresponds to the one labeled as d in our previous study of the Ω c odd-parity resonances [33], where it was tentatively assigned either to the Ω c (3090) or the Ω c (3119) observed by the LHCb Collaboration in the Ξ cK 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 − Σ c state around 2800 MeV that will complete the sex-tet. The Σ c (2800) clearly fits into this picture since it is observed in the Λ c π channel [1].
Recently there has been an analysis of the Ξ 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 Ξ c states can be identified with the experimental Ξ c (2790), Ξ c (2930), Ξ c (2970), Ξ c (3055) and Ξ c (3080). Whereas the Ξ c (2790) would be a 1/2 − state, the Ξ c (2930), Ξ c (2970), Ξ c (3055) and Ξ c (3080) could be either 1/2 − or 3/2 − ones. Compared to this approach, our model identifies the experimental Ξ c (2790) and Ξ c (2930) as 1/2 − states, and the Ξ c (2815) and Ξ 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 * and light vector mesons with baryons are not completely fixed by HQSS or chiral symmetries, thus allowing for different assumptions.

B. Ξ b excited states
With regards to the bottom sector, the Ξ b (6227) resonance has been recently measured by the LHCb experiment [9], with Γ Ξ b (6227) ∼ 18 MeV. Its quantum numbers, though, remain unknown, whereas the observed decay channels are Λ 0 b K − and Ξ 0 b π − (see Table I). We start again by revising the previous results of Ref. [24] with B = −1, S = −1, I = 1/2 (Ξ 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 baryonmeson channels as well those for the experimental decay channels of Ξ b (6227) are shown in Table III. We obtain nine states, which are the bottom counterparts of the Ξ 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 Λ 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 III. 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 SU(6) lsf ×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.
In Fig. 3 we display the evolution of the masses and widths of the dynamically generated Ξ b states as we vary the cutoff from Λ = 1000 MeV (triangles) to Λ = 1400 MeV, for J = 1/2 (upper plot) and J = 3/2 (lower plot). The square and its bar represent the position of the Ξ b (6227) resonance, and the error for its mass and width, respectively. We show the experimental result (Ξ b (6227)) for both J = 1/2 and J = 3/2 because its quantum numbers have not been determined yet. Additionally, in Table IV, we collect the masses and the widths of the b 1 to b 9 states with J = 1/2 or J = 3/2, together with the couplings to the dominant baryon-meson channels and the couplings to the decay channels of the Ξ b (6227), for Λ = 1150 MeV as in the charm sector. We also indicate the SU(6) lsf × HQSS, SU(6) and SU(3) irreducible representations of these states.
We might try now to assign the experimental Ξ b (6227) to any of our states, while determining the negative parity baryons with B = −1 belonging to the same 3 * and 6 SU(3) representations. The observed decay modes, , of the resonance support that this state should have 1/2 − spin-parity, assuming S−wave. Moreover, the j π ldof = 0 − component should be also quite relevant, which according to the couplings collected in Table IV makes plausible its identification either with the b 1 or b 2 states. The evolution displayed in the upper plot of (LHCb [32]) would be part of the 3 * 2 and 3 * 4 multiplets forming a HQSS-doublet [24]. These antritriplets should be completed by another HQSS-doublet of Ξ b and Ξ * b states, b 5 and b 8 , that according to Fig. 3 and Table IV should have masses around 6250 MeV and could be seen in the Σ ( * ) bK and Ξ ( * ) b π modes. Coming back to the Ξ b (6227), it belongs to a j π ldof = 0 − −sextet that should be completed by J = 1/2 Σ b and Ω b states. The recent Σ b (6097) resonance seen by the LHCb Collaboration [47] in the Λ b π channel nicely fits in this multiplet. Relying again in the equal spacing rule, we could foresee the existence of a J = 1/2 Ω b odd parity state with a mass of around 6360 MeV that should be observed in the Ξ bK channel. Some molecular Ω b states were predicted previously in Ref. [48], but all of them above 6.4 GeV.

IV. CONCLUSIONS
In this work we have explored the possible molecular interpretation of several experimental excited Ξ c and Ξ b states. We have used a coupled-channel unitarized model, that is based on a SU(6) lsf ×HQSS-extended WT baryon-meson interaction, within the on-shell approximation. We have paid a special attention to the depen-dence of our predictions on the renormalization scheme, so as to assess the robustness of our results.
We have presented a molecular interpretation for the experimental Ξ c (2790), Ξ c (2815), Ξ c (2930), Ξ c (2970) and Ξ b (6227) states, and have predicted the spin-parity quantum numbers of the latter three resonances. We have found that the Ξ c (2790) state has a large molecular Λ cK component, with a dominant j π ldof = 0 − configuration, and discussed the differences between the 3/2 − Λ c (2625) and Ξ c (2815) states, finding that they cannot be SU(3) siblings. We have also predicted the existence of other Ξ c −states, not experimentally detected yet, being two of them siblings of the two poles that form the Λ c (2595). Interestingly, the recently discovered Ξ c (2930) and Ξ c (2970) are found to be HQSS partners.
The flavor-symmetry content of the framework has also allowed us to understand the nature of the Σ c (2800) and Σ b (6097) states, for which we have determined their spinparity. Moreover, we have predicted several states, some of them displayed in Figs. 4 and 5 (marked with a ? symbol). Among them, we stress the Ω b (6360) state, with a dominant Ξ bK contribution, in the sextet where the Σ b (6097) and Ξ b (6227) are located, together with the Ξ b (6240) and Ξ * b (6240) states, partners of the HQSS doublet Λ b (5912) and Λ b (5920) discussed in [24].