The $P_{cs}(4459)$ pentaquark from a combined effective field theory and phenomenological perspective

The observation of the $P_{cs}(4459)$ by the LHCb collaboration adds a new member to the set of known hidden-charm pentaquarks, which includes the $P_c(4312)$, $P_c(4440)$ and $P_c(4457)$. The $P_{cs}(4459)$ is expected to have the light-quark content of a $\Lambda$ baryon ($I=0$, $S=-1$), but its spin is unknown. Its closeness to the $\bar{D}^* \Xi_c$ threshold -- $4478\,{\rm MeV}$ in the isospin-symmetric limit -- suggests the molecular hypothesis as a plausible explanation for the $P_{cs}(4459)$. While in the absence of coupled-channel dynamics heavy-quark spin symmetry predicts the two spin-states of the $\bar{D}^* \Xi_c$ to be degenerate, power counting arguments indicate that the coupling with the nearby $\bar{D} \Xi_c'$ and $\bar{D} \Xi_c^*$ channels might be a leading order effect. This generates a hyperfine splitting in which the $J=\tfrac{3}{2}$ $\bar{D}^* \Xi_c$ pentaquark will be lighter than the $J=\tfrac{1}{2}$ configuration, which we estimate to be of the order of $5-15\,{\rm MeV}$. We also point out an accidental symmetry between the $P_{cs}(4459)$ and $P_c(4440/4457)$ potentials. Finally, we argue that the spectroscopy and the $J/\psi \Lambda$ decays of the $P_{cs}(4459)$ might suggest a marginal preference for $J = \tfrac{3}{2}$ over $J = \tfrac{1}{2}$.


Introduction
The discovery by the LHCb collaboration of three hiddencharm pentaquarks [1] -the P c (4312), P c (4440) and P c (4457) -has triggered intense theoretical efforts to decode their nature, in particular whether they are molecular [2][3][4][5][6][7][8][9][10][11] or not [12][13][14][15][16]. Recently a new hidden-charm pentaquark has been found [17] -the P cs (4459) 0 -which we will simply denote as P cs in this work. This pentaquark has been observed in the J/ψΛ channel, from which it can be deduced that its quark content is ccsqq with q = u, d. Its mass and width are M P cs = 4458.8 ± 2.9 +4.7 −1.1 MeV , Γ P cs = 17.3 ± 6.5 +8.0 but the statistical significance of the signal is merely 3.1 σ. Besides, its spin and parity have not been determined yet. It is also worth noticing that predictions of P c and P cs pentaquarks [18][19][20] have been there long before their eventual observation. The P cs pentaquark lies a few MeV below theD * Ξ c threshold -4478.0 MeV in the isospin symmetric limit -suggesting a strong molecular component [21][22][23]. However there are at least other two nearby thresholds: theDΞ ′ c andDΞ * c ones at 4446.0 and 4513.2 MeV, respectively (i.e. 32.0 and 35.2 MeV away from theD * Ξ c threshold). If the spin of the P cs pentaquark is J = 1 2 ( 3 2 ), it will mix with theDΞ ′ c (DΞ * c ) channel, which will result into a molecular picture more complex Correspondence to: mpavon@buaa.edu.cn than that of the P c pentaquarks (i.e. that of a single channel DΣ c orD * Σ c molecule). Here we will consider how the aforementioned coupled channel dynamics affects the spectrum of a molecular P cs . If we consider the possible isoscalarDΞ ′ c and DΞ * c molecular states, one quickly realizes that owing to SU(3)flavor and heavy quark spin symmetry (HQSS) it is possible to make predictions [24]. For theD * Ξ system this is not the case though and we will have to resort to phenomenology to relate its interaction with the already known non-strange molecular pentaquarks. If this is done, the molecular description of the P c (4312/4440/4457) and P cs (4459) pentaquarks turns out to be coherent, as we will explain in the following lines.
The manuscript is organized as follows: in Sect. 2 we briefly explain the non-relativistic effective field theory we will use to describe the molecular pentaquarks. In Sect. 3 we discuss the symmetry constraints of the pentaquarks. Sect. 4 is devoted to the power counting of the coupled channels affecting the P cs . In Sect. 5 we explain how to estimate the low energy constants of the effective field theory from meson-exchange saturation. In Sect. 6 we will show an accidental symmetry between the potentials of the P c (4440/4457) and P cs pentaquarks. In Sect. 7 we discuss the size of the hyperfine splitting between the J = 1 2 , 3 2 P cs pentaquarks. In Sect. 8 we consider the decay of the P cs pentaquark into J/ψΛ depending on its spin. Finally, we summarize our conclusions in Sect. 9 and explain a few technicalities in Appendices A and B.

Effective field theory description
Before explaining how symmetries inform the pentaquark spectrum, first we will briefly explain the effective field theory (EFT) formalism we follow. We will describe interactions among heavy hadrons with a non-relativistic contact-range potential of the type with C an unknown coupling constant, where this coupling can be further decomposed into a sum of irreducible components C = R λ R C R , with R denoting some quantum-number / representation, λ R some coefficient / operator and C R the particular coupling that applies in each case. This type of contactrange potential often appears in lowest-(or leading-) order EFT descriptions of hadron-hadron interactions (concrete examples with full derivations can be found in Refs. [25][26][27][28] for antimesonmeson molecules and in Ref. [29] for pentaquarks). Of course this is true provided that the one-pion-exchange potential, which is the longest range piece of the hadron-hadron interaction, is weak and thus subleading [30,31] (otherwise it should be included at lowest-order). The previous contact-range potential is singular though and has to be regularized, which we do by introducing a regulator function f (x) and a cutoff Λ, i.e.
where the coupling now depends on the cutoff C = C(Λ). For the regulator we will choose a Gaussian, f (x) = e −x 2 , and for the cutoff we will use the range Λ = 0.5 − 1.0 GeV. Finally this potential is included in a dynamical equation, such as Schrödinger or Lippmann-Schwinger, for obtaining predictions. If we choose Lippmann-Schwinger and are interested in poles of the scattering amplitude, i.e. bound/virtual states or resonances, we can simply solve where φ is the vertex function, which is defined as the the wave function Ψ times the propagator (φ(p) = [M th + p 2 /(2µ) − M mol ] Ψ (p)), V the potential, M th the mass of the threshold (i.e. the sum of the masses of the two hadrons comprising a molecular candidate), µ their reduced mass and M mol the mass of the hadronic molecule we want to predict.

Light-flavor and heavy-quark symmetries
Symmetry constrains the potential binding the molecular pentaquarks. If we begin by considering the three known P c pentaquarks, in the molecular picture they are thought to bē DΣ c andD * Σ c bound states. From the SU(3)-flavor perspective the P c 's are composed of a triplet charmed antimeson and a sextet charmed baryon, which together can couple into the octet and decuplet representations of SU(3), i.e. 3 ⊗ 6 = 8 ⊕ 10. The flavor structure of the potential is thus with H c = D, D * or D s , D * s and S c = Σ c , Σ * c , Ξ ′ c , Ξ * c or Ω c , Ω * c representing an arbitrary charmed meson or baryon, C O and C D the octet and decuplet couplings and λ O and λ D a coefficient that depends on the particular antimeson-baryon configuration considered (they are explained in detail in Ref [24]). From the HQSS perspective the potential between two heavy hadrons can only depend on the spin of the light quarks inside them. For the triplet charmed meson and sextet charmed baryon the light-spins are S L = 1 2 and S L = 1, respectively, which couple to 1 2 ⊗ 1 = 1 2 ⊕ 3 2 . However it is more compact to express the light-quark spin structure of the potential in terms of light-spin operators: with C a and C b couplings that represent the spin-independent and spin-dependent pieces of the potential, respectively, and σ L and S L the spin-operators for the light-spin degrees of freedom within the charmed meson and baryon (for the notation in terms of light-spin check for instance Ref. [32], while the channel-bychannel potential can be found in Ref. [29]). From the SU(3)-flavor and HQSS structure we have just explained it is already possible to derive the existence ofDΞ ′ c andDΞ * c molecular states [24]. First we notice that the standard molecular interpretation of the P c (4312) pentaquark is that it is a I = 1 2D Σ c bound state. Thus the decomposition of the P c (4312) potential is i.e. the octet SU(3)-representation and the HQSS part of the potential that is independent of the spin of the light-quarks. Any other molecular pentaquark with the same decomposition will have the same potential as the P c (4312) and consequently, will be likely to have a similar binding energy. Among these pentaquarks we have the I = 0DΞ ′ c andDΞ * c systems, for which the potential reads where again only the octet, spin-independent piece of the contactrange potential (C O a ) is involved. From now on we will simply write C a = C O a , as the octet configuration is the only one we are considering in this work.
If we now consider the P cs , its most natural molecular interpretation isD * Ξ c . This two-body system is not connected tō DΞ ′ c ,DΞ * c andDΣ c neither by SU(3)-flavor nor HQSS symmetries. From SU(3)-flavor symmetry, the P cs pentaquark contains a triplet charmed antimeson and antitriplet charmed baryon and is a combination of a singlet and an octet, i.e. 3 ⊗3 = 1 ⊕ 8. The concrete flavor structure of the potential is unessential though, as we are only considering the I = 0, S = −1 sector (i.e.D * Ξ c ). Regarding HQSS, the antitriplet charmed baryon contains a diquark with S L = 0, from which we expect a trivial light-spin structure owing to 1 2 ⊗ 0 = 1 2 . The potential reads with no spin dependence whatsoever and T c = Λ c , Ξ c representing a generic antitriplet charmed baryon. In addition to this, thē H c T c andH c S c systems can couple by means of a transition potential of the type with E b a coupling, σ L the spin-operator for the light-quark within the charmed meson and ǫ L the polarization vector of the light-diquark in the sextet charmed baryon. The couplings D a and E b can be further decomposed in isospin and flavor representations, but this is not necessary for the set of molecules we are considering. Putting the pieces together for theD * Ξ c , if we consider the coupled channel bases B(J = 1 2 ) = {DΞ ′ c ,D * Ξ c } and B(J = 3 2 ) = {D * Ξ c ,DΞ * c } we will have the following potentials: By including these potentials in a bound state equation such as the coupled-channel extension of Eq. (4) we can calculate the mass of the P cs .

Power counting and coupled channel dynamics
EFTs are expected to be power series in terms of the expansion parameter (Q/M), where Q and M represent the characteristic low and high energy scales of the system, respectively. For molecular pentaquarks Q is of the order of the pion mass (m π ≃ 140 MeV) or the wave number of the bound state (i.e. γ = 2µB 2 ∼ 206 MeV for the P cs as aD * Ξ c molecule), while M will be of the order of the vector meson mass (m ρ ≃ 770 MeV). This suggests the expansion parameter where we have identified the wave number γ with the light scale Q. We will now compare this number with the expected size of coupled channel effects.
If we are interested in the mass difference between the J = 1 2 and 3 2D * Ξ c bound states, i.e. the hyperfine splitting, the relevant coupled channels are of theH c T c -H c S c type, i.e. Eq. (12), which can break the spin degeneracy. There areH c T c -H c T c coupled channel effects too (e.g.D * s Λ c -D * Ξ c ), but they do not generate a dependence on the light spin. For the coupled channel dynamics relevant to theD * Ξ c system (independently of whether they generate spin dependence), their expected size with respect to the diagonal interaction is [30,31] respectively, where B 2 = 19.2 MeV is the binding energy of a molecular P cs and ∆ CC the mass gap of the listed coupled channels. This indicates that only theDΞ ′ c andDΞ * c channels are expected to be larger than the size of subleading corrections. The next channel in importance,D * s Λ c , does not break spin degeneracy, as previously mentioned, and in addition its size is subleading. Finally, though theD * Ξ * channel will indeed contribute to the hyperfine splitting, its size is strongly suppressed with respect toDΞ ′ c andDΞ * c and thus we will not take it into account.
For analyzing the possible impact of the coupled channel dynamics, we will do the following calculation (i) Consider the P cs pentaquark to be a J = 1 2 or 3 2 molecule, which in analogy with Ref. [5] we will call scenarios A and B, respectively. (ii) Consider different E b /D a coupling ratios: with this ratio fixed, the D a coupling can be determined from the P cs pentaquark (and the C a from the P c (4312) one). Then we check how the hyperfine splitting changes with this ratio.
The result of these calculations is shown in Fig. 1 for scenarios A and B and a cutoff Λ = 0.5 −1.0 GeV. The hyperfine splitting grows quickly with the E b /D a ratio and it is sizable even for small ratios. This can be understood from the power counting of contact-range theories [34], in which a coupling generating a bound state near threshold is fine-tuned, thus explaining how the effect of a comparatively small E b is amplified by the fact that D a can generate a molecular P cs . We also notice that the for the same E b /D a ratio the hyperfine splitting will be considerably larger in scenario A, which has to do with the fact that D a is also larger in this scenario: coupled-channel dynamics require M(J = 1 2 ) > M(J = 3 2 ), which in turn forces D a to be larger in scenario A if E b 0.
However, without being able to estimate the E b /D a coupling it will be not possible to know the hyperfine splitting. From power counting arguments the size of each of these couplings will be [34] where the superscript (R) refers to a renormalized coupling (we elaborate below): D (R) a is said to be enhanced (i.e. its size is larger than expected owing to the existence of a bound state close to threshold), while E (R) b is natural (i.e. its size can be determined from standard or naive dimensional analysis arguments). The renormalized couplings D (R) a and E (R) b (which is the type of couplings for which the arguments of Ref. [34] were originally developed) refer loosely speaking to the parts of the couplings that do not depend on the cutoff. But here we are working instead with the bare (or running) couplings D a = D a (Λ) and E b = E b (Λ), which explicitly depend on the cutoff. Nonetheless the previous power counting estimates apply to the bare couplings for specific cutoff ranges: (i) for couplings of natural size (e.g. E b ), this will be the case irrespectively of whether the cutoff is soft (Λ ∼ Q) or hard (Λ ∼ M), while (ii) for couplings of unnatural size (e.g. D a ) the enhancement ideally requires a soft cutoff Λ ∼ Q, with the size of the coupling reverting to its natural size as the cutoff becomes harder [35,36]. However, in the molecular pentaquarks the separation of scales is far from perfect: as explained in Appendix A for the P cs pentaquark and the Gaussian regulator we use here, D a (Λ) coincides with its power counting estimation for Λ ∼ 0.9 GeV (i.e. within the cutoff range we use). We thus expect that in a first approximation the previous relations will hold in the Λ = (0.5 − 1.0) GeV range, i.e.
from which we get E b /D a ∼ Q/M ∼ 0.27, yielding an estimated hyperfine splitting of ∆M A ∼ 12 − 35 and ∆M B ∼ 7 − 15 MeV in scenarios A and B, respectively, where there is still a noticeable cutoff dependence. It is nonetheless possible to improve over the previous picture by including a renormalization factor to better connect the bare and running couplings: where we discuss the derivation of the factor F in Appendix A. Concrete calculations show that F = (1.47 − 0.91) in the Λ = 0.5 − 1.0 GeV cutoff window used in this work, leading to the hyperfine splittings ∆M A ∼ 27−40 and ∆M B ∼ 12−13 MeV in scenarios A and B, respectively, which happen to display less cutoff dependence (though they are still of the same order as our original estimation). It is worth noticing that the reason behind these elaborations is that we are applying power counting arguments to non-observable quantities (which are allowed to have a strong dependence on the cutoff). Actually, besides the standard P cs (4459) single peak found in the J/ψΛ invariant mass distribution, the LHCb collaboration also reports a second possible two-peak solution [17] involving two P cs pentaquarks with masses which, if they are both to be interpreted primarily asD * Ξ c bound states, will result in the hyperfine splitting ∆M = 12.9 ± 4.6 MeV. In principle this is compatible with scenario B. But both scenarios A and B use the standard single peak solution as the reference input. Had we determined the D a and E b couplings from the two-peak solution instead, then the E b /D a ratio would have been for the cutoff range Λ = 0.5−1.0 GeV. This ratio is in fact compatible with the power counting estimation of Q/M ∼ 0.27. In , which are expected to be LO effects in the EFT description. This description contains three independent couplings C a , D a , E b , where the first one (C a ) is determined from the P c (4312) pentaquark and SU(3)-flavor symmetry. For the other two couplings we do as follows: scenario A (B) assumes that the observed P cs pentaquark is the J = 1 2 ( 3 2 ) configuration; then for a given E b /D a ratio we determine D a from the P cs mass. Finally, we show the hyperfine splitting ∆M(P cs ) = M(P cs , J = 1 2 ) − M(P cs , J = 3 2 ) in both scenarios. the following lines we will resort to phenomenological information for further elucidating the E b /D a ratio.

Meson exchange saturation
The problem we have is that there are three couplings (C a , D a , E b ) of which we can only determine two (C a from the P c (4312) and D a or E b from the P cs (4459)). Yet, if we use phenomenology it might be possible to find relations among these couplings and thus determine the three of them. In particular we will focus on light-meson saturation, i.e. the idea that the contact-range couplings of a given EFT are saturated by the exchange of light mesons [37,38]. Here we choose the novel saturation procedure of Ref. [39], which we explain below.
Standard saturation maps the finite-range S-wave potential generated by the exchange of a light-meson, V M ( q 2 ) (with q = p ′ − p the exchanged momentum, where for S-wave we can express the potential as a function of q 2 ), into a contact-range coupling C sat by taking the limit which is expected to work for Λ close to the mass m of the exchanged meson. However, if the potential vanishes in this limit, we will obtain C sat = 0. For instance, the potentials generate exactly the same finite-range potential in r-space, namely the Yukawa potential where the difference between the two is a distribution Of course, at this point we have to discuss the impact of formfactors, which modify the light-meson exchange potentials as follows where F M is the aforementioned form-factor (with the most used parametrizations being of the multipolar type) while Λ M is the form-factor cutoff, which should not be confused with the EFT cutoff Λ. For a local form-factor the resulting r-space potentials will be also local, and the Dirac-delta will acquire a finite-size: where δ (3) F represents a Dirac-delta that has been already smeared out by the form-factor. Yet, the characteristic scale Λ M of these finite-range effects is expected to be larger than the mass of the exchanged meson, i.e. Λ M > m (otherwise the effect of said exchanged light-meson will be washed out by the form factors). For instance, for the Bonn-B potential [40] Λ σ = 1.9 − 2.0 GeV and Λ ρ/ω = 1.85 GeV for the M = σ and ρ/ω mesons, respectively, while for the CD-Bonn potential [41] we have Λ σ = 2.5 GeV, Λ ρ = 1.31 GeV and Λ ω = 1.5 GeV.
Thus for the range of cutoffs in which saturation is expected to work we have (Λ ∼ m) < Λ F , which implies that the previous Dirac-delta is unimportant: independently of whether the potential is derived from derivative interactions (V ′ M ) or not (V M ), the potentials at mr ∼ 1 will be similar and thus the saturation of the couplings should follow suit. That is, if the renormalization scale is similar to the exchanged meson mass, the potentials V M and V ′ M are expected to lead to approximately the same saturated coupling. This is achieved with the convention i.e. by extracting the residue of the potential at q 2 = −m 2 , which effectively recovers the expectations from the mr ∼ 1 behavior of the potential. Ref. [39] explicitly checked this method with the one-pion-exchange potential (with an arbitrary coupling strength) as a specific example. In Appendix B we include a detailed comparison between the standard saturation procedure of Ref. [38] and the one presented here for the particular case of the nucleon-nucleon system, which indicates that both saturation methods yield comparable results. Now, if we consider the scalar meson σ, in the non-relativistic limit it generates a spin-independent potential which can contribute to the saturation of C a and D a (but not C b or E b ): where m S the mass of the sigma meson and g σi its coupling, and the indices i, j = 1, 2, 3 referring to theD ( * ) , Ξ ( ′ / * ) c and Ξ c , respectively. Independently of the saturation method we will obtain with F a = C a , D a the generic name for the spin-independent couplings. The proportionality constant is in principle unknown, but we will assume it to be similar for all the couplings. We remind that saturation is expected to work for cutoffs close to the mass of the meson being exchanged, Λ ∼ m S in this case. The vector mesons (ρ and ω) generate a more complex potential which can be expanded in a multipole expansion similar to the one we have for electromagnetic interactions. There are electric-and magnetic-type components (indicated by the subscripts a and b) with and O L12 = 0, σ L · ǫ L or σ L · S L forH c T ,H c T c -H c S c and H c S c , respectively, where the dots represent S-to-D-wave components (we assume they will not appreciably contribute to the saturation of the EFT couplings). In the vector-exchange potential the i, j indices refer to the different hadrons involved, g Vi to the electric-type couplings, f Vi to the magnetic-like ones, m V to the vector meson mass and M a typical hadronic mass scale. The V Va component contribute to the C a and D a couplings as where F a = C a , D a , with no difference between the standard and modified saturation procedures, i.e. Eqs. (23) and (29). For the V Vb component to contribute to the C b and E b couplings in a non-trivial way we have to use the modified saturation procedure of Ref. [39] (i.e. Eq. (29)) in which case we obtain where Here we notice that with the standard saturation method we would have arrived to This is what happens for instance in Refs. [8,33], which in principle consider the complete set of LO interactions we have here, i.e. Eqs. (6), (11) and (12), but set the F b -type couplings to zero leading to pentaquark predictions without sizable hyperfine splittings. Yet, as we will see, from the point of view of phenomenology the choice of Refs. [8,33] is completely justified: the size of the F b -type of coupling is expected to be considerably smaller than the F a -type, which prompts the approximation F b = 0. In contrast from the EFT point of view the inclusion of the F b couplings is justified in terms of power counting.
The only thing left is to determine the couplings: for the σ we will use the linear-sigma model [42] and the quark model (as used in Ref. [43]), from which we get that the coupling of the sigma to the nucleon is g σNN = √ 2M N / f π ∼ 10.2 (for the f π ≃ 132 MeV normalization). For the charmed mesons with only one light-quark we end up with g σ1 = g σNN /3 ≃ 3.4, while for the charmed baryons with two light-quarks we have g σ2 = g σ3 = 2 g σNN /3 ≃ 6.8. We note that here we have assumed the same coupling of the sigma with the q = u, d, s quarks, as happens for instance in the quark model of Ref. [44], where we notice that this pattern can either arise from a singlet sigma (not necessarily a realistic assumption) or alternatively from a negligible coupling of the octet component of the sigma to the light quarks, see for instance Ref. [45] for a more detailed discussion. In either case, this runs counter to the standard expectation that the strange and non-strange components of light mesons would decouple (owing to the OZI rule), leading in the case of the sigma to a vanishing coupling to the strange quark. However it has been argued that the OZI rule does not work so well in the 0 ++ sector [46][47][48][49], i.e. a nonsinglet sigma might potentially have a non-negligible coupling with the strange quark. Besides, the singlet and octet mixing angle of the sigma might very well be far away from the angle decoupling strange and non-strange components [50]. Yet, even though the sigma coupling to the strange quark might not be necessarily suppressed, it might not be as strong as with the light-quarks. This can generate SU(3)-flavor breaking effects which we will discuss later, where the naive expectation will be that the contribution of the sigma to the contact-range couplings will be weaker in theDΞ ′ c andD * Ξ c molecules than in theDΣ c one. For further details we refer to the discussion around Eq. 73, where we notice that the previous expectation seems to be challenged by the experimental information currently available 1 .
For the vector mesons, we have electric-and magnetic-like couplings, for which we will resort to Sakurai's universality and vector meson dominance [51][52][53], i.e. the mixing of the neutral vector mesons with the electromagnetic current, which can be encapsulated in the substitution rules [54,55] where e is the proton charge, g = m V /2 f π ≃ 2.9 the universal rho coupling constant (with m V the vector meson mass and f π ≃ 132 MeV the pion weak decay constant), A µ the photon field, ρ 3 µ the neutral rho field (the superscript refers to the isospin index, a = 3 for the neutral component), ω µ the omega field and µ a Lorentz index. The parameter β indicates the degree of vector meson dominance: if the electromagnetic couplings of the heavy hadrons were to be completely dominated by the vector mesons, we will have β = 1. Otherwise it will be β < 1. Usually β ≃ 0.9 is estimated [56][57][58] (which can be traced back to the ratio of the couplings required for the ρ → γγ and ρ → ππ processes, i.e. √ 2 f π / f ρ ∼ 0.9). By applying these substitution rules to the interaction of the heavy hadrons with the vector mesons, we will get the electromagnetic interaction for the light quarks within the heavy hadrons and from this, we can determine the g Vi and f Vi couplings. In practical terms this means that g Vi and f Vi are proportional to the light-quark contribution of the electric charge and magnetic moments of the heavy hadrons, respectively. For the E0 couplings we get g V1 = g V2 = g V3 = βg. For the M1 couplings, we will first make the decomposition f V = κ V g V and use the choice M = m N for the mass scale with m N ≃ 940 MeV the nucleon mass. From this κ V = 3 2 (µ u /µ n.m. ), with µ u /µ n.m. the magnetic moment of the u-quark within light-diquark pair inside the Ξ + c /Ξ ( ′ / * )+ c charmed baryon expressed in units of the nuclear magneton (µ n.m. ). If we make use of the quark model a second time, we obtain κ V1 = κ V2 = 3 2 µ u /µ n.m. ≃ 2.9 for the charmed antimeson and sextet strange charmed baryon. For the antitriplet charmed baryon we have instead κ V3 = 0, which is a consequence of the two light-quarks within the Ξ c baryon being in a spin-0 configuration. The E b coupling involves a M1 antitriplet to sextet transition for which κ V (3 → 6) = 3 2 µ(3 → 6)/µ n.m. = 3 2 (µ u /µ n.m. ) ≃ 2.9, where µ(3 → 6) refers to the uquark magnetic moment in the3 → 6 light-diquark transition 2 within Ξ + c → Ξ ( ′ , * )+ c . 1 In particular D a (as obtained from the P cs (4459)) seems to be more attractive than C a (as obtained from the P c (4312)), see Eq. (77), which runs counter with the expectation that sigma exchange should be weaker in the former case. Yet more accurate experimental measurements are needed to confirm whether this is really the case. 2 For translating the magnetic moment of the light-diquark into the one for the heavy baryons we can use the relations µ q (Ξ ′ c ) = 2 3 µ q (6) and µ q (Ξ c → Ξ ′ c ) = ±µ q (3 → 6)/ √ 3, with µ q (X) the magnetic moment of a particular light-quark q within "X" (X being a baryon or a light-diquark in the antitriplet or sextet configuration) and with the sign of the3 → 6 transition depending on the ordering of the light quarks in the3 flavor wavefunction (which is flavor antisymmetric).
With the previous couplings and setting m S ≃ 550 MeV and m V = (m ρ + m m ω )/2 ≃ 775 MeV, we arrive to the ratios: The first one of these ratios was checked in Ref. [54], where it was respected at the 30% level. If we determine C a and C b as in Ref. [5], we obtain |C b /C a | = 0.16 for Λ = 0.75 GeV, which reproduces the absolute magnitude of Eq. (38) at the 25% level. However we note that in Ref. [5] the sign of C b depends on which are the spin of the P c (4440) and P c (4457) pentaquarks: the sign is correctly reproduced if the P c (4457) is J = 1 2 . The second one (C a = D a ) appears for instance in Ref. [18], which predicted the P cs , and Refs. [19,33,59] as a consequence of the universality of the vector meson coupling. The recent work of Ref. [60] uses the chiral quark model to guess the contactrange couplings, leading to C a ≃ D a too. The third relation has not been previously used, as E b is usually set to zero (e.g. Ref. [33], where their E b equivalent is calledμ 24 ). Finally it is worth noticing the following where this relation actually does not depend so much on saturation being correct or accurate, but rather on the fact that the light-meson exchange potentials are identical under the approximations we have made. In the following lines we will explore the consequences of this relation.

Accidental symmetry in the pentaquark potential
Light-meson exchange actually suggests a very interesting relation between the P c (4440), P c (4457) and the two spin states ofD * Ξ c . The S-wave light-meson exchange potential in thē H c S c system is which for the J = 1 2 , 3 2D * Σ c systems reads from which we expect the hyperfine splitting to be proportional to In comparison for the coupledH c T c -H c S c system, the corresponding potential reads where if we consider vector meson exchange, vector meson dominance, the quark model relations for the charmed baryon magnetic moments and SU(3) symmetric sigma exchange, we will have This will receive small corrections from η-exchange (that work in the direction of making |W b | > |V b |), which we will ignore as they are small. Now, in the limit where theDΞ ′ c -DΞ c andDΞ * c -D * Ξ c thresholds are degenerate and have the same mass, there will be two eigenvalues for this potential, which correspond to the linear combinations with potential eigenvalues In particular, in this limit the hyperfine splitting between the "± states" will be that is, expected to be similar to the hyperfine splitting between the J = 1 2 and 3 2D * Σ c states. Of course, this accidental symmetry in the potential will be broken by the fact that there is a mass gap between the involved coupled channels. The effect of this mass gap will be to decrease the hyperfine splitting and to force the "±" sign in Eq. (50) as to make the state corresponding to the lower threshold the most attractive configuration. If we compare the characteristic momentum scale of the coupled channel exchange potential (i.e. the vector meson mass m V ) and the momentum scale of the coupled channel dynamics (Λ CC = 2µ∆ CC ≃ 266 and 279 MeV), the ratio is 0.34 and 0.36 for theDΞ ′ c andDΞ * c channels, respectively. Together they add to 0.70, i.e. we expect the hyperfine splittings to be about 30% of the expected value were not to be a mass gap. For the J = 1 2 ( 3 2 ) configuration, theD * Ξ c threshold is heavier (lighter) than theDΞ ′ c (DΞ * c ) one, which forces the most repulsive (attractive) sign configuration in Eq. (51).
As a consequence, if the previous approximations hold, the hyperfine splitting of the twoD * Ξ c pentaquarks will be similar to the one of the twoD * Σ c pentaquarks, i.e.

|M(D *
However if the effect is to be reduced by a 70%, as suggested by the scale comparison, we will end up with a 5.1 MeV splitting. Though the sign of the hyperfine splitting might be protected owing to power counting and the nature of the coupled channel dynamics, its size will be diminished owing to the finite mass gaps between the channels. Besides, the uncertainties in the couplings of the light-meson exchange potential are also large. In the following we will see how this expected effect holds when compared with different error sources. For doing the explicit calculations we first obtain the C a and C b couplings from the masses of the P c (4440) and P c (4457) pentaquarks, i.e. the calculation of Ref. [5], where we notice that for the hyperfine splitting it does not matter which spin is each pentaquark 3 for Λ = 0.5 (1.0) GeV, with the hyperfine splitting ∆M(P cs ) = 4.6 (4.7) MeV , which indeed indicates a reduction of the coupled channel effects owing to the finite mass gap between the channels, and where from now on we define Yet, even if this accidental symmetry is greatly reduced in the hyperfine splittings, it is worth noticing that the predictions we obtain from using the effective C a and C b couplings describing the P c (4440) and P c (4457) pentaquarks are basically compatible with the experimental mass of the observed P cs pentaquark.

The hyperfine splitting
Now we will analyze the possible size of the hyperfine splittings with the (admittedly approximate) information we have derived from light-meson exchange saturation. We begin with the simplest of the relations, that is: and determine C a from reproducing the P c (4312) pentaquark, which yields C a = −1.19 (−(2.17 − 0.80)) fm 2 for Λ = 0.75 (0.5 − 1.0) GeV, where we will use a central value of the cutoff close to the expected scale at which saturation works (i.e. Λ ∼ m V ), while we will still consider the 0.5 − 1.0 GeV cutoff range for estimating the regulator uncertainties. If we set E b = 0, it will lead to M(P cs , J = 1 2 , 3 2 ) = 4467.5 (4466.9 − 4468.1) MeV , (59) that is, we obtain two degenerate P cs pentaquarks. If we allow for E b 0 we can effectively fit one of the pentaquarks to the experimental mass. We obtain 4 .: where it should be noticed that scenario B is automatically chosen, as for C a /D a ∼ 1 the mass of the P cs pentaquark can only be reproduced if J = 3 2 . The hyperfine splitting will be ∆M(P cs ) = 14.3 (12.7 − 16.5) MeV , which is definitely larger than the estimation from the phenomenological symmetry in the pentaquark potential. The ratio E b /D b = 0.34 (0.24 − 0.53) will also be larger than the expectation from saturation, Eq. (40).
Alternatively we can assume that the saturation relations in Eqs. (39) and (40) and the hyperfine splitting is which is closer to the one we derived from the phenomenological symmetry, i.e. Eq. (56). Now there are a series of (potentially large) uncertainties related to the previous relations. The most obvious one is the existence of SU(3)-breaking corrections between the coupling in theDΣ c (strangeness S = 0) and theDΞ ′ c ,DΞ * c (strangeness S = −1) systems, i.e. between the P c (4312) and P ′ cs , P * cs pentaquarks. Their couplings, which are identical if SU(3)-flavor symmetry is exactly preserved, will differ by a correction where δC / F a indicates the correction to the C a coupling. From a comparison with the pion and kaon weak decay constants, f π ≃ 130 MeV and f K ≃ 160 MeV, we expect the size of the SU(3)-breaking effects to be of the order of At the end ("Two-peak solution") we compare with the two-peak fit of Ref. [17], from which an hyperfine splitting of 12.9 ± 4.6 MeV is expected. If this figure turns out to be confirmed by future experimental studies, this will discard the "EFT A", "EFT+RGA A", "Accidental" and "Saturation II" estimations.
Alternatively, we can consider this problem from the point of view of chiral symmetry 5 , in which the C a coupling can be decomposed into quark-mass independent and quark-mass dependent pieces [61]: the quark-mass dependent piece comes from terms in the Lagrangian where the quark-mass matrix is inserted between the hadron fields, with these terms expected to be subleading (as in this case we are expanding around the massless quark limit). This quark-mass dependence can be translated into a quadratic dependence on the mass of the pseudoscalar Nambu-Goldstone bosons, which can be schematically written as where C [n,m] indicates n (m) insertions of the quark mass matrix between the charmed antimeson (baryon) fields, m π and m K are the pion and kaon masses and Λ χ ∼ 1 GeV is the chiral symmetry breaking scale. Rearranging the terms, the difference can be rewritten in terms of a chiral symmetry breaking (χS B) correction, which takes the form 5 Notice that the EFT we are using here is not a pionless EFT, but rather a pionful EFT for which pion (and pseudoscalar meson) exchanges are considered to be subleading and thus not explicitly included in the leading order description. Thus chiral symmetry considerations can play an explicit role.
From the arguments about quark-mass dependence shown above, the size of this correction is expected to be of the order where it is interesting to notice that its size is identical to the previous estimation based on f π and f K . The experience in the light-baryon sector suggests breakings larger than the previous estimation [62], where repulsion increases with the number of strange quarks. The recent discovery of the Z cs [63] allows for a comparison of the couplings required to reproduce the Z c /Z * c [64] and Z cs [65] as virtual states, with attraction apparently increasing with the strange quark content (though the uncertainties are large).
Independently of the derivation, if SU(3)-breaking corrections reach a certain level they would lead to unboundDΞ ′ c and DΞ * c systems. In particular this will happen for for Λ = 0.75(0.5 − 1.0) GeV, where we indicate which pentaquark we are referring to in the parentheses. In this regard the previous 20% estimations indicate that the SU(3)-symmetry partner of the P c (4312) pentaquark is still likely to bind. The eventual observation of aDΞ ′ c molecular pentaquark could bring light to this issue.
A significant effect which might influence the size of SU(3)flavor breaking corrections is the nature of sigma meson exchange: if the sigma were not to couple with the strange meson, saturation will suggest which is in the limit between binding and not binding for P ′ cs (though in the no-binding case, the P ′ cs pentaquark might still survive as a virtual state). Additionally, in this scenario from saturation we will expect which would imply that the P cs could also very well be close to not binding (except for the increase in the relative strength of E b , which would help in the J = 3 2 configuration). With these coupling ratios, the masses of the molecular P cs pentaquarks would be with the J = 1 2 state in the second Riemann sheet with respect to the D * Ξ c threshold (i.e. it is a shallow virtual state / resonance). That is, only the J = 3 2 state is an actual bound state. In this later case the hyperfine splitting is compatible with zero and can even change signs as the J = 1 2 is allowed to be a virtual state. However this large SU(3)-breaking in the direction of making theD * Ξ c less bound does not reproduce the experimental mass of the P cs pentaquark (unless we allow for E b /D a ∼ 1.23, which is a considerable deviation from the saturation relations and would lead to a hyperfine splitting of 50.2 MeV). Thus it is unlikely that SU(3)-breaking would be as large as in Eq. (73), at least if its effect is to reduce attraction in theDΞ ′ system.
Finally, it is possible to make a comparison with the twopeak fit included in Ref. [17], which leads to two P cs pentaquarks with masses 4454.9 ± 2.7 MeV and 4467.8 ± 3.7 MeV, respectively. Concrete calculations assuming that the lighter (heavier) P cs is a J = where the central value and the spreads correspond to Λ = 0.75 (0.5 − 1.0) GeV, as usual. If the two-peak fit ends up being confirmed in future studies, the previous indicates a bit more attraction than expected for the D a coupling (but compatible with the errors we would expect for a phenomenological determination of the D a /C a ratio) and that the size of E b seems to be underestimated by the phenomenological arguments we provide (yet compatible with the power counting estimates we proposed). We summarize the different estimations we have considered along this work in Table 1. These indicate that, though it is not possible to determine the hyperfine splitting accurately from theory alone, power counting arguments and phenomenological approximations suggest it might be in the 5 − 15 MeV range (i.e. compatible with the two-peak fit in [17]), with the J = 3/2 pentaquark being the lighter state. For comparison, Refs. [18,19] predict degenerate P cs pentaquarks. Meanwhile, a recent calculation in the one-boson-exchange model generates a ∆M(P cs ) = (2.4 − 20.0) MeV splitting, which also comes from theDΞ ′ c -D * Ξ c andD * Ξ c -DΞ * c coupled channel dynamics [66]. In contrast Ref. [60] obtains ∆M(P cs ) = −6.0 MeV from two-pion-exchange (TPE). This is interesting as the naive expectation would be that TPE is of order (Q/M) 2 in the EFT expansion and we might expect it to play a minor role. Thus, the role of TPE might indeed deserve further attention in the future. However there is the practical limitation that this calculation will also require (Q/M) 2 corrections to the contact-range potential, i.e. more unknown parameters.
Regarding which is the spin of the P cs pentaquark, from the different predictions in Table 1 it seems that the J = 3 2 DΞ * c configuration might provide a better match to its experimental mass, though uncertainties are too large to draw definite conclusions. It is also important to stress that the experimental determination of resonance masses usually relies in using the Breit-Wigner parametrization. Other parametrizations might yield different masses, as happened with the P c (4312) [9], the Z c (3900) / Z c (4020) [64], and the Z cs (3985) [65]. Thus the mass of a molecular P cs might not coincide with the experimental determination of Ref. [17].

Decay into J/ψΛ
It is interesting to notice that the observation of the P cs (4459) in the J/ψΛ channel might provide further circumstantial evidence of its spin. If we decompose the J = 1 2 , 3 2D * Ξ c system into its heavy-and light-quark spin components, we find which is to be compared with where S H and S L refer to the heavy-and light-quarks spin. If the decay preserves HQSS [67] the expect the following relation between matrix elements: which for degenerateD * Ξ c states implies that the partial decay widths of the J = 1 2 and 3 2 configurations will show a 1 : 4 ratio. In fact phenomenological calculations seems to support these ratios, with Ref. [33] giving a 1 : 1.78 ratio for the amplitudes / couplings and Ref. [68] yielding 1 : 4.35 for the partial decay widths. Of course, this does not determine the spin of the P cs (4459), but nonetheless indicates that, ceteris paribus, the probability of discovering the J = 3 2 molecule in the J/ψΛ invariant mass distribution might be larger than for its J = 1 2 partner. But this conclusion is dependent on the production rates from the Ξ b decays, which have been recently investigated in Ref. [69], suggesting that the production rate of a J = 3 2D * Ξ c pentaquark would be 4.9 times the one for its J = 1 2 partner. Ideally, it would be possible to adapt the methods of Refs. [70,71] (originally formulated for the three P c pentaquarks) to analyze the invariant mass distribution data of the new P cs .

Conclusions
The P cs (4459) is the latest piece of the pentaquark puzzle. Its closeness to theD * Ξ c threshold suggests that it might be a bound state of these hadrons. Then the question is what is the connection of the P cs (4459) with the well-known P c (4312), P c (4440) and P c (4457) in the molecular hypothesis. At first sight the answer is unclear: theD * Ξ c system is not directly connected by SU(3)-flavor and HQSS with theDΣ c andDΣ * c systems, which are the usual molecular interpretations of the three P c pentaquarks. However if we resort to phenomenological arguments then we can bridge the gap between the new P cs and the previous P c 's, resulting in a coherent description of these four pentaquarks. From vector meson dominance and the quark model, we point out a possible accidental symmetry between the potentials in the P c and P cs sectors, though owing to its phenomenological nature large deviations are to be expected.
There are two possible spin configurations for a molecular P cs , which in principle should be degenerate. In this regard the explicit inclusion of theDΞ ′ c -D * Ξ c andD * Ξ c -DΞ * c coupled channel dynamics, which is required by power counting arguments, breaks this degeneracy and thus might have important implications for spectroscopy. This mechanism generates a sizable hyperfine splitting which we estimate to be in the 5 − 15 MeV range. Incidentally, this estimation is in line with the proposed two-peak solution in [17].
In general the predicted mass of the J = 3 2 P cs pentaquark are closer to its experimental value than its J = 1 2 partner, which might be interpreted as favoring the former spin assignment. But theoretical errors in the masses make it unpractical to determine the spin of the new pentaquark from spectroscopy alone. In this regard, the partial decay widths of the two spin configurations to J/ψΛ, where the P cs have been discovered, approximately differ by a factor of four, making the J = 3 2 configuration considerably easier to detect in this channel. Of course this is only true provided all other effects are similar. Independently of their spin, the existence of two possibleD * Ξ c molecules tends to hold up well within the expected uncertainties of a phenomenological approach. As happened with the original P c (4450), future experiments could further determine whether there are really twoD * Ξ c states and which is their mass difference.

Acknowledgments
We would like to thank Li-Sheng Geng

A Power counting arguments with running coupling constants
In this work we have made use of power counting estimations of the size of the contact-range couplings appearing in the EFT description of the molecular pentaquarks. However the contact-range couplings are cutoff dependent and the aforementioned estimations were originally formulated for the renormalized couplings [34], the definition of which is scheme dependent. Yet, it has been argued that these estimations do indeed apply to running couplings if the cutoff is sufficiently soft [35,36]. But the previous arguments are qualitative in nature: here we will see how to modify dimensional estimations in a quantitative manner for their use with running couplings. This will be particularly useful for the dimensional estimation of the ratio of the D a and E b coupling constant employed in Sect. 4 (check Eq. (17)).
The easiest example will be a contact-range theory with only one channel, in which we determine the coupling from the condition of reproducing the two-body scattering amplitude. For instance, if we consider the scattering of aD * antimeson and a Ξ c baryon, the T-matrix can be written as where f (x) is the regulator we use for the contact-range potential and with I 0 (k, Λ) the loop function: the exact evaluation of which depends on the details of the regulator. If we want the T-matrix to be (exactly) cutoff independent at a given reference momentum k R , this will lead to the condition where Λ 1 and Λ 2 are two different cutoffs. If we choose to renormalize the scattering amplitude at the bound state pole, i.e. at k R = iγ R = i 2µB 2 , the D a coupling will be given by This coupling will exactly reproduce its power counting estimation for a privileged cutoff Λ * where for a Gaussian regulator Λ * ≃ 4.3 2µB 2 , which for the particular case of the P cs (4459) yields Λ * ≃ 886 MeV. However, if the cutoff is different from this privileged value, the dimensional estimations will have to be corrected as follows It happens that Λ * ∼ 0.9 GeV, which implies that corrections will be small for the range of cutoffs considered in this work (i.e. Λ = 0.5 − 1.0 GeV). We notice that, in contrast with what is expected in Refs. [35,36], the estimation of Λ * in the P cs pentaquark is rather large and can be hardly considered to be a soft scale. However Refs. [35,36] deal with the two-nucleon system and the deuteron happens to be considerably more shallow than the P cs . Indeed, repeating the previous arguments for the deuteron yields Λ * ≃ 196 MeV, which is of the order of the pion mass and more in line with the expectations of Refs. [35,36]. The inclusion of coupled channel effects can be taken into account by considering the matrix version of the previous renormalization group equation which ensures that the T-matrix is cutoff-independent at the renormalization point k = k R and where k ′ R refers to the reference momentum in the second channel. If we write the previous equations in coefficients, we end up with where det(V C ) = C a D a − E 2 b is the determinant of the contactrange potential. Combining the last two equations, we end up with If for simplicity we assume the same Λ * as in the uncoupled channel case and that for Λ * all the power counting estimations are followed as expected (i.e. |D a | = 2π µγ R , |C a | = 2π µγ ′ R and |E b | = 2π µM , with γ R , γ ′ R the wave number in each of the coupled channels and M = m ρ ), we will end up with a correction factor of for Λ = 0.5−1.0 GeV, which happens to be close to the original estimation. However this small change reduces the cutoff dependence of the hyperfine splitting (which is now a renormalizable quantity), yielding 27 − 40 and 12 − 13 MeV in scenarios A and B respectively (check the discussion around Eq. (19) in the main text).

B Standard and novel saturation procedures in the two-nucleon system
In this appendix we compare the standard and modified saturation procedures for the particular case of the S-wave twonucleon contact-range interactions. Standard saturation is known to work well when comparing the EFT contact-range couplings with a series of OBE potentials [38]. Thus the natural question at this point is whether this is still the case with the novel saturation method of Ref. [39].
However nucleons are not as heavy as the charmed mesons and baryons, which means that relativistic corrections are often included in the light-meson exchange potentials. This implies that the saturation formulas we have to use are somewhat more involved than the ones we obtained in Sect. 4. Actually, the corresponding formulas for the standard saturation method (which also include form factors) can be found in Ref. [38]. It happens that the novel saturation method merely generates and additive factor in the standard saturation relations of Ref. [38]: Here we notice that, for potentials of the q 2 /(m 2 + q 2 ) type, the novel saturation procedure can be encapsulated in the following substitution rule which merely amounts to the inclusion of an additive term in the potential to manually remove the Dirac-delta term, thus justifying the rule in Eq. (95). Now in the two-nucleon system we also find relativistic corrections that follow the general form k 2 q 2 + m 2 = − 1 4 q 2 q 2 + m 2 + 1 4 which we have rewritten as the sum of a purely local and nonlocal terms (i.e. the terms proportional to q 2 and ( p 2 + p ′2 ), respectively). For the saturation of the local term we will use again the substitution rule of Eq. (96). From the previous, the modifications for a scalar, pseudoscalar and vector meson are where M N is the nucleon mass, g S (g P ) the coupling constant for the scalar (pseudoscalar) meson, g V and f V the electric-and magnetic-like couplings for the vector meson, σ 1(2) the spin operators for nucleons 1(2) and F M ( q 2 ) the form-factors, which accept the expansion with Λ M the form-factor cutoff for the particular meson M under consideration. If the exchanged meson is an isovector (e.g. the ρ), the isospin operator τ 1 · τ 2 will have to be included in the previous expressions. If we consider the momentum expansion of the S-wave contact-range potential p ′ |V C |p = C 0 + C 2 (p ′2 + p 2 ) + . . . , the novel saturation procedure generates the following modifications for the 1 S 0 and 3 S 1 partial waves in the two-nucleon system (where we have used the spectroscopic notation 2S +1 L J with S , L, J the intrinsic, orbital and total spin, respectively). For a scalar meson, we have to add the terms while for a pseudoscalar meson we will have For vector meson exchange we add the terms Finally, Ref. [38] uses the following definition for the coupling constantsĈ = 4π C 0 and C = 4π C 2 , which we will also use in what follows for a more convenient comparison. Putting all the pieces together, for the particular case of the Bonn B potential [40] we obtain C novel 1 S 0 =Ĉ std 1 S 0 + δĈ1 S 0 = (−0.117 − 0.134) · 10 4 GeV −2 = −0.251 · 10 4 GeV −2 , C novel 1 S 0 = C std 1 S 0 + δC1 S 0 = (1.276 + 0.178) · 10 4 GeV −4 = 1.454 · 10 4 GeV −4 , for the singlet, to be compared withĈ N 2 LO 1 S 0 = {−0.160, −0.158 · 10 4 GeV −2 and C N 2 LO 1 S 0 = {1.134, 1.135} · 10 4 GeV −4 (i.e. the results for the EFT couplings in the two-nucleon system at nextto-next-to-leading (N 2 LO) as determined in Ref. [38], where the brackets indicate their expected variation within the formalism of the aforementioned reference), while for the triplet we getĈ novel 3 S 1 =Ĉ std 1 S 0 + δĈ1 S 0 = (−0.101 − 0.091) · 10 4 GeV −2 = −0.192 · 10 4 GeV −2 , C novel 3 S 1 = C std 1 S 0 + δC1 S 0 = (0.660 + 0.178) · 10 4 GeV −4 = 0.838 · 10 4 GeV −4 , and C N 2 LO 3 S 1 = {0.637, 0.587} · 10 4 GeV −4 for EFT. From this, in the case of the singlet channel with the Bonn-B potential, the novel saturation method underperforms the standard one, while for the triplet channel their deviations with respect the EFT couplings are similar. However this comparison is potentialdependent: Ref. [38] considers a total of six phenomenological potentials, while here we limit ourselves to Bonn-B, i.e. the easiest one on which to apply saturation.