On the hidden charm pentaquark-s in $\Lambda_b \to J/\psi K^- p$ decay

In a previous work we presented a theoretical analysis of the $\Lambda_b \to J/\psi K^- p$ reaction from where a recent experiment by the LHCb collaboration at CERN claimed the existence of two hidden charm pentaquarks, $P_c(4380)^+$ and $P_c(4450)^+$. In that work we focused only on the $\Lambda(1405)$ and $P_c(4450)^+$ signals and discussed the possible explanation of this pentaquark state within the picture of a dynamical meson-baryon molecule made up mostly from $\bar D^* \Sigma_c$ and $\bar D^* \Sigma^*_c$ components. In the present work we improve upon the previous one by considering the total $K^- p$ and $J/\psi p$ data including all the relevant resonances contributing to the spectra, and discuss the possible nature of both $P_c(4380)^+$ and $P_c(4450)^+$. We also discuss several important topics, like the effect of the contact term in the reaction, the viability of reproducing the data without the $P_c(4380)^+$ and the possible quantum numbers assignment to these pentaquarks.


INTRODUCTION
In the early beginning of the quark model Gell-Mann and Zweig already mentioned that apart from the standard qq and qqq configurations for mesons and baryons respectively, there could exist some multiquark configurations [1,2]. Concrete calculations for the case of pentaquarks were done by Hogaasen and Sorba [3] and Strotmann [4]. Much excitement was generated by the claims of an observation of a pentaquark state, Θ + , at Spring8/Osaka in a photonuclear reaction [5]. A similar method was used to analyze the γd → pnK + K − reaction where again a claim was made for the Θ + pentaquark [6]. After a period of excitement where the peak was observed in most laboratories, searches with better statistics and analysis methods started to report negative results and the issue was closed. A detailed report on this issue can be seen in [7]. An important work to clarify the issue was done in [8,9], were it was shown that the experimental peak observed in [6] was a consequence of the analysis method in [6], where only the K + and K − were observed and the K + n invariant mass was constructed with a prescription for the unmeasured n momentum, that was incorrect and artificially produced a peak in the "Θ + " region.
With this precedent, the claim of two pentaquark states in the LHCb experiment [10,11] should have been taken with caution, but given the thoroughness of the experimental analysis, the result was difficult to challenge. Yet, issues concerning the implementation of unitarity in the experimental analysis, and the lack of a tree level contribution, which is unavoidable from a theoretical point of view, were raised [12]. One of the aims of the present paper is to discuss in detail these issues and show how some accidental circumstances make the present exper-imental analysis overcome this problem, thus providing extra support to the experimental claims.
The precedent of the former "Θ + " pentaquark unjustified claims has not prevented a wave of excitement among theoreticians, who have proposed a variety of possible explanations for these two states. One of the reasons for it is that predictions of hidden charm baryon states had been made before. Indeed, in [13,14] baryon states of hidden charm were found in the study of the interaction of theDΣ c -DΛ c ,D * Σ c -D * Λ c coupled channels as the main building blocks, together with the η c N and J/ψN states, plus decay channels in the light sector. Related studies were done in [15], where bound states ofDΣ c andD * Σ c were also found. In [16], using an admixture of SU (6) and Heavy Quark Spin symmetry, HQSS, states of hidden charm similar to those predicted in [13,14] were also found. Further studies were done and in [17] similar results to those of [13,14] were found, using HQSS and the local hidden gauge approach as tools to evaluate the matrix elements of the interaction. A quark model was used in [18], where some hidden charm baryons states were also obtained. All these works share similar qualitative results, but differ in the predictions of the masses of the particles by as much as 200 MeV up in [18] to 200 MeV down in [16] with respect to those found in [13,14,17]. Further work on this line is done in [19], using an admixture of Vector-Baryon and Pseudoscalar-Baryon states in coupled channels which allow one to have a better hold on the decay width of the states.

FORMALISM
The core of the present analysis is the identification of the most relevant mechanisms contributing to the amplitude to describe the Λ b → J/ψK − p experimental data. We do not intend to obtain a better fit than the one carried out in the experimental analysis [10], which included many different Λ resonances in addition to the pentaquarks, considering possible different quantum number assignments for them and all angular dependences relative to the decay products of the J/ψ. Our intention is to implement a good enough amplitude, but as simple as possible in order to discuss the relevant theoretical issues regarding the role played by the different contributions to the Λ b → J/ψK − p, specially the pentaquarks and the dominant Λ resonances. Nevertheless, regarding the K − p and J/ψ p invariant mass distributions our approach is quite accurate since the angular information mentioned above, considered by the experimental analysis, is integrated in these observables. The procedure followed in the present work is also rather different from the one followed in [10]. Instead of using the helicity formalism, we construct explicit amplitudes suited to the excitation of the different resonances and their quantum numbers. In addition, we explicitly consider the K − p formation at the tree level and its interaction with coupled channels in s-wave.
In ref. [21] only the dynamically generated resonances, Λ(1405) and P c (4450), were included, and it was shown in [48] that the relevant mechanisms in the Λ b → J/ψ K − p decay are those depicted in Figs. 1(a) and 1(b). The K − p is produced at the quark level as depicted in Fig. 1(a). It is also interesting to recall that in [48] it was also found that the J/ψπΣ states were not produced at this tree level. First the weak interaction produces the cc state (J/ψ) and an s quark and then the u, d and s remaining quarks hadronize into a meson-baryon pair. Note that the u and d quarks are spectators during the entire process, as justified in ref. [48]. Within the framework of the chiral unitary approach, the Λ(1405) resonance, of crucial role in the present decay, is produced dynamically from the final K − p state interaction implementing unitarity in the different meson baryon channels with the Λ(1405) quantum numbers. In this way, and without the need to include the Λ(1405) resonances as explicit degrees of freedom, two poles were obtained for the Λ(1405) resonance [49]. Recently the positions of the poles were more finely obtained at √ s = 1352 − 48i MeV and 1419 − 29i MeV [50], with the lowest mass pole coupling mostly to πΣ and the highest mass one toKN . Therefore its implementation into the Λ b → J/ψ K − p decay is depicted in Fig. 1(b), and it is given by where the s-wave, isospin 0,K-nucleon unitarized scattering amplitude from ref. [50], (thick circle in Fig. 1(b)). The parameter α 1 is a free parameter, to be fitted later on, accounting for the elementary production process, Fig. 1(a). Note that in addition to the previous amplitude we have to add to the total amplitude the J/ψ K − p tree level contact term contribution, Fig. 1(a), which is shown in ref. [48] to have the same weight, α 1 , as the term in Eq. (1). This is a non-trivial result. Thus the contribution of the mechanisms in Fig. 1(a) and 1(b) is Note that this tree level contribution, accounted for by the 1 addend in Eq. (2), interferes with the amplitude of Eq. (1). Although constant nonresonant terms were used in the fit of [10], the output of the fit did not return a significant contribution of these terms. We shall try to understand this feature, and this issue will be widely discussed in the results section. As mentioned in the Introduction, in refs. [13,17,19], several poles were obtained (see table II of ref. [17] and table 8 of ref. [19]) when implementing unitarity in coupled channels in s-wave considering the channels J/ψN , η c N , DB andD * B, with B baryon charmed states belonging to the 20 representations of SU (4) with J P = 1/2 − and 3/2 − . Exploring the possibility that some of these poles (or a mixture of several ones) could correspond to the experimentally found P c pentaquarks is one of the aims of the present work. In such a case, the production of the pentaquark in our model would proceed through the mechanism depicted in Fig. 1(c): The J/ψ p pair initially produced undergoes final state interaction, which is accounted for by the J/ψ p → J/ψ p unitarized scattering amplitude represented by the thick circular dot in Fig. 1(c). In this particular case, the J/ψ p → J/ψ p unitarized scattering amplitude resembles very much a Breit-Wigner [17], therefore, and for the numerical evaluation carried out in the present work, the mechanism of Fig. 1(c) can be effectively accounted for by a term proportional to with G J/ψp the J/ψp loop function and m Pc (Γ PC ) the mass (width) of either of the pentaquarks. The term in Eq.
(3) has to be multiplied by a momentum structure which depends on the different possible quantum number assignment of the pentaquarks. (See Appendix for explicit details). The pole positions of the amplitudes obtained in refs. [13,17,19] provide directly m Pc and Γ Pc , but with uncertainties in the mass of the order of 200 MeV [19]. Therefore we fine tune these values to the experimental results of ref. [10], m Pc = 4380 MeV and Γ Pc = 205 MeV for the lowest pentaquark and m Pc = 4449.8 MeV and Γ Pc = 40 MeV for the highest one, which lie indeed in between the different values obtained in refs. [13,17,19]. The coupling of the dynamically generated resonance to J/ψ p, g J/ψ p in Eq. (3), was determined to be of the order of 0.5 in refs. [13,17,21]. However, we have included in the amplitude of Eq. (3) a free parameter α i , (i = 2 for P c (4450) and i = 3 for P c (4380)), to better fit the experimental data, which should have a natural value of the order of 1. The presence of the J/ψp loop function, G J/ψp , in Eq. (3) has to be included when assuming the resonances to be dynamically generated, since in such a case it is always produced by the scattering of an initial J/ψp pair. Nevertheless, the G J/ψp factor has little impact in the global fit. We can also consider the scenario where one or both pentaquarks have J P different to 1/2 − or 3/2 − which are the quantum numbers of the states generated by the chiral unitary approach [13,17,19]. In this scenario, and if the J/ψp pair was in p-wave, the corresponding pentaquark would carry J P = 1/2 + , 3/2 + or 5/2 + and its contribution should be added to the previous terms leading to the Λ(1405) (see Appendix). Nonetheless, and as we will explain below, since the qualitative output of our fits are similar irrespective of the quantum number assignments for the pentaquarks, we will only work out the 5/2 + case (in addition to the previous 1/2 − and 3/2 − ) for simplicity since only one possible partial wave for the kaon is possible, (L ′ = 2), in that case, (see Appendix , table I, for details). In this case the mechanism for the pentaquark production would be simply that depicted in Fig. 2a and, in order to establish its amplitude, we must remove G J/ψp in Eq. (3) and multiply it by a momentum dependence according to its quantum numbers as shown in Eqs. (12), (17) and (26).
On the other hand, in the analysis of ref. [21] no further Λ resonances, (in addition to the Λ(1405)), were considered since that study was based only on the Λ(1405) and the pentaquark filtered signal. However, since in the present work we aim at reproducing the total K − p and J/ψ p invariant mass, we must add the relevant Λ contributions to the process, Fig. 2b. In ref. [10], 13 different Λ resonances were included in the fit. In our analysis it is enough to consider only those which gave a sizable contribution to the final cross section in [10] which, in addition to the Λ(1405) discussed above, turn out to be Λ(1520) The Λ resonances (except the Λ(1405) which is explained above, see Eq. (1)) are parametrized by a Breit-Wigner shape with Flatté parametrization of the width by up to a factor containing the spin and momentum structure as explained in the Appendix. The widths of the Λ resonances have been taken energy dependent in the following way: with Γ oi the on shell width, and f i the branching ratio of the Λ i into the dominant decay channels, j, obtained from the PDG [51]. In Eq. (5) B(L, q, q o ) is the Blatt-Weisskopf penetration factor for L-wave [52], where q (q o ) is the K or p momentum at the K − p center of mass frame at M K − p (m i ). Furthermore contact terms with the different spin and angular momentum of the J/ψ and kaon are also considered as explained in the Appendix. These contact terms were found to be negligible in the experimental fit [10] but we will explain in the results section that some of them could play indeed an important role.
Altogether, the Λ b → J/ψK − p differential decay rate is given by The interferences between the different Λ mechanisms and the pentaquarks and contact terms depend on the different quantum numbers of the particular Λ resonances and the pentaquarks. The total amplitude squared |T (M K − p M J/ψp )| 2 (of course averaged over initial spins and added over final ones) is given in the Appendix for the different cases. The α i and C i in Eqs. (22), (23), (24), (25) and (26) and are free parameters to be fitted in our analysis of the following section.

RESULTS AND DISCUSSION
We have carried out different fits to the experimental data [10] considering the following possibilities for the spin-parity of the pentaquarks (J P A , J P B ) where J P A and J P B stand for the spin-parity of the P c (4380) and P c (4450) respectively: and (5/2 + , 5/2 + ). Nevertheless, and advancing some results, neither of them provides a remarkably better fit than the rest, and all of them produce qualitatively similar results. However, experimentally [10] a best fit was obtained for J P assignments of (3/2 − , 5/2 + ), and also acceptable fits were obtained for (3/2 + , 5/2 − ) and (5/2 + , 3/2 − ) and "other combinations are less likely", but there was a clear preference for two states with opposite parity [10]. The experimental analysis [10] is more complete than the one carried out in the present work since it takes into consideration more Λ resonances and, specially, all angular dependences relative to the decay products of the J/ψ, which we do not consider. However the angular dependence is not relevant for the invariant mass distributions, which are the observables we fit, therefore the extra conclusions drawn by the experimental analysis [10] about the J P of the pentaquarks must certainly come from that extra information used in the experimental analysis, (note that there are more than 150 parameters fitted in [10] while there are just 19 in our analysis). Yet, it is worth stressing again that it is not the intent of the present work to improve or even be on a par with the already good experimental analysis. Our focus is on specific theoretical issues for which the approach followed here is good enough. Therefore, since the qualitative discussion is similar irrespectively of the J P pentaquark as-signments in our model, we will discuss upon results for (3/2 − , 3/2 − ) assignment, one of the chiral unitary options, unless other case is explicitly stated.
First we show in the (a) and (b) panels of Fig. 3 the result of the fit to the experimental data [10] and the individual contributions of the different resonances. (The 1 + GT label stands for the term 1 + (2) and (22) which essentially produces the Λ(1405), as explained in the previous section). We can see that the global fit is quite fair, given the simplified version of the model compared to the analysis done in the experimental work [10]. Note specially the important strength of both pentaquarks, P c (4380) and P c (4450) in the J/ψp mass distribution. We have only fitted the K − p mass distributions up to M K − p < 2 GeV and the J/ψp up to 4.8 GeV since our chiral unitary model for the Λ(1405) cannot be extrapolated to further higher energies and the reduced range is preferable and sufficient for our discussion. Panels (c) and (d) in Fig. 3 reflect data and results for J/ψp mass distribution implementing the kinematic cuts 1.70 GeV < M K − p < 2 GeV and M K − p > 2.0 GeV respectively. (See figure 8 of ref. [10]). The data in panels (c) and (d) are not fitted, thus the curves therein are output of the calculation. It is worth noting that the experimental data for the K − p and J/ψ p invariant mass distributions are not corrected for experimental setup acceptance in ref. [10]. However, in the experimental paper, phase-space curves for Fig. 3a and Fig. 3b are provided, which of course are affected by the acceptance. Therefore, comparing those curves to the corresponding theoretical phase-space three body distribution, we have renormalized each experimental datum such that the phase-space agrees with the actual one and such that the areas below both invariant mass distributions are the same (i.e. same total Λ b width). This acceptance correction cannot be performed in the theoretical results in panels (c) and (d) since the acceptance modified phasespace is not provided for such kinematic cuts in [10]. This is one of the reasons of the rough agreement between our theoretical calculation and experiment in panels (c) and (d). However, the main reason is that those cuts filter events in the higher part of the K − p spectrum which are little relevant in the global fit, or actually are not fitted at all. Nevertheless, it is worth noting that our results of Fig. 3c and d will be used just qualitatively in the coming discussions below.
In ref. [21] it was pointed out that the experimental support for the existence of the P c (4380) state was not as clear as for the P c (4450) one. The reason was that the Argand plot for the P c (4380) is not as clean as that of the narrow P c (4450) [10]. On the other hand there is not a clear bump for the P c (4380) in the total J/ψp experimental invariant mass distribution. To shed some light into this issue we next carry out a fit removing the P c (4380) term. The result is shown in Fig. 4a   This fit without P c (4380) is just slightly worse compared to that in Figs. 3a and b, (about a 20% bigger χ 2 /dof), but it is specially very similar or even better in the lowest region of the J/ψp mass distribution. We have traced the reason for the good agreement in the low J/ψp mass region when removing the P c (4380) to the J P = 5/2 + contact term (the 1 addend in Eq. (26)). (The shape of this contact term by itself corresponds to the label "only C 5/2 " in the figures.) We see that with a slight increase of this contact term when the P c (4380) contact term is removed (compare curves "only C 5/2 " between figures Figs. 3b and Fig. 4b), a similar effect than the one of the P c (4380) can be mostly accommodated. Very similar results and conclusions are obtained for the other (J P A , J P B ) possibilities studied in the present work. Therefore we must conclude that the fit to only the invariant mass distributions is not enough to draw firm conclusions about the existence of the P c (4380) state or the spin-parity of the pentaquark states, since many different (J P A , J P B ) possibilities yield similar results. The reason why the experimental analysis get (3/2 − , 5/2 + ) as the best option and (3/2 + , 5/2 − ) and (5/2 + , 3/2 − ) also acceptable and "other combinations are less likely" [10] must then be traced to some other observable be-  yond just the K − p and J/ψp mass distribution, like the angular dependence of the J/ψ decay products, etc.
On the other hand, if we next direct the attention to the phase space distribution in Fig. 4c, we see that the kinematic cut implemented in that plot causes by itself a sharp kink around M J/ψp ∼ 4.2 GeV which is responsible for some of the accumulation of strength seen in the experimental data around that region. Therefore this apparent accumulation of strength at low M J/ψp is not an indication by itself of the existence of a nearby resonance, like the P c (4380).
In order to illustrate the similarity of the fits among different J P pentaquark assignments mentioned at the beginning of this section, we show in Fig. 5 the result of the fit for the (3/2 − , 5/2 + ) case which is the case for the best fit in the experimental work [10]. In our case, however, we get a fit with a very similar χ 2 /dof to that in Fig. 3a and b. Let us next discuss one of the most important issues to be addressed in the present work, that is the effect of the contact term contribution of Fig. 1(a), represented by the 1 addend in Eqs. (2) and (22). Although a nonresonant term of this type was considered in the fit of [10], the output returned a negligible contribution from this source, and, in principle, it could play an important role in the results. In the experimental fit [10], a Breit-Wigner shape is considered for the different resonances but for the Λ(1405) the width has two components to account for the Flatté effect. However, as explained above in the formalism section, in the chiral uni- tary approach this resonance is generated (actually two of them) from the K − p interaction and accounted for by the 1 (2) and (22), (we will just call it 1 + GT term in the following). Note that we include the 1 addend for the contact term which seems to be negligible in the final fit of [10]. Indeed, we see in Fig. 3a and b that the contribution of the 1 addend by itself (curves labeled "only 1") is very large and it is crucial to produce the final Λ(1405) strength and shape in the invariant mass distributions from its interference with the GT term (curves labeled "only GT"), which is also very large by itself and with a shape very different from a Λ(1405) resonance. However, we are going to show that a fortuitous combination of facts renders the approach of [10] very similar to the fully unitary approach that we follow here.
Elaborating on this latter issue, let us recall that the two Λ(1405) are basically obtained from the interaction of the coupled channels πΣ andKN . We shall call T ij the transition matrices from channels 1 (πΣ) and 2 (KN ). Actually what one would expect to be approximately a Breit-Wigner is T 11 itself, not (1 + GT ) 11 . Therefore, in order to mimic the experimental approach to the Λ(1405) we have performed a different fit implementing the substitution where m R , Γ o and α are adjusted to approximately reproduce the t I=0 KN,KN amplitude in the Λ(1405) resonance region and Θ is the step function. The term proportional to α in Eq. (8) is included in order to account for the Flatté effect, which is also incorporated in the analysis of ref. [10]. The result of the fit is shown in Fig. 6. We can see that the result is very similar to Fig. 3a,b, in spite of the fact that one is neglecting the contact term in the latter analysis.
This certainly requires a detailed explanation: If we plot (1 + GT ) 11 and −T 11 we see in Fig. 7 that they are very approximately proportional, (in that figure −T 11 has been multiplied by a constant factor such that the absolute values are the same at the Λ(1405) mass). Actually the curves for the absolute value of (1 + GT ) 11 and −T 11 are almost the same (up to a global constant). This is not true in a general case for a unitarized scattering amplitude since, in the chiral unitary approach and for two coupled channels, we have However, it turns out that V 12 is about a factor 3 smaller than V 11 and V 22 and therefore, taking into account that V 11 is negative and smooth in the energy region considered, we have An example where the approximation in Eq. (10) does not hold is the J/ψp scattering that produces the pentaquarks in the chiral unitary approach. In this case the  potential V J/ψp,J/ψp is much smaller than those for the other channels [17] (mostlyD * Σ c andD * Σ * c ). Let us call channel 1 the J/ψp and assume for simplicity that there was only another channel,D * Σ c , (channel 2). In this case the potential matrix element V 12 is of the same order of magnitude than V 22 and T 12 is larger than T 11 [17]. Therefore (1 + GT ) 11 is not proportional to T 11 in the J/ψp scattering case.
Coming back to the meson-baryon interaction in I = 0, L = 0, that produces the Λ(1405), one could naively think, in a first impression, that even if V 12 was not small, Eq. (10) should also hold next to a pole since T 11 and T 12 contain both the singularity of the pole which would factorize out. However this is true if only one pole is present, but as already mentioned above, there are two poles associated to the Λ(1405) resonance. If we call these two poles A and B, respectively, we would have that which are indeed not proportional since the couplings of the different poles, g i , to the different channels are different, (g A 1 = 2.71, g B 1 = 2.96, g A 2 = 3.06 and g B 2 = 1.96 [50]). In Fig. 8 we can see that T 11 and T 12 are actually very different.
In Fig. 9 we show the comparison between (1 + GT ) 11 and the Breit-Wigner parameterization of the Λ(1405) of Eq. (7). In this figure we can see the similarity between both amplitudes, which makes the fit almost equivalent using any of them.
In conclusion, when fitting the Λ(1405) with a Breit-Wigner with a Flatté width, as the experimental analysis does [10], and neglecting the contact term, the result is equivalent to having considered 1 + GT as in the chiral unitary approach. But this is true in this particular case by chance. Note that it is a combination of the K − p tree level plus rescattering (1 + GT ) that makes the two approaches equivalent. As mentioned, this is not trivial. Things could have also been different if, instead of having the K − p at tree level and not πΣ, one had the opposite  situation, with πΣ at tree level and not K − p. We would then have something similar to the "only GT " contribution of Fig. 3(a), instead of the (1 + GT ), and the spectra of K − p would be different.

SUMMARY AND CONCLUSIONS
We have addressed some interesting issues from a theoretical point of view regarding the possible pentaquark states that show up in the Λ b → J/ψK − p decay from the LHCb experimental data [10]. In order to model theoretically the process, we have improved over a previous work [21], implementing the K − p and J/ψp final state interaction and including explicitly the relevant Λ resonances not generated in the former rescattering. The K − p interaction generates dynamically the double Λ(1405) pole structure implementing unitarity withKN and πΣ coupled channels in s-wave. In the J/ψp final state the possibility of the pentaquark to be dynamically generated is also allowed. With this simplified model, as the number of Λ resonances included with respect to the experimental analysis is concerned, but accurate enough to reproduce the experimental data, we perform several fits to the LHCb invariant mass distributions [10] for the different spin-parity possibilities for the pentaquarks.
One of the conclusions obtained is that, with the only fit to the K − p and J/ψp mass distributions, the existence of the P c (4380) state cannot be undoubtedly claimed (unlike the P c (4450) state), since we get not much worse results removing this pentaquark from the fit. We have traced the origin of this similarity in the results with or without P c (4380) to the effect of a nonresonant term which provides some of the strength in the absence of the P c (4380) resonance. Furthermore, we also obtain similar results for the different possibilities of the spin-parity of the pentaquarks. Therefore the claims regarding the existence of the P c (4380) pentaquark and the spin-parity assignments of both pentaquarks made in the experimental analysis [10] cannot be inferred just from a fit to the K − p and J/ψp mass distributions. Thus, it would be most welcomed if the experimental group could singled out the observables that show unambiguously the existence of the P c (4380) pentaquark and the spin-parity of both of them.
On the other hand, we have widely discussed the important role played by the tree level contact elementary Λ b → J/ψK − p production and the K − p rescattering. The interference of these two terms gives rise to the Λ(1405) resonance in our case, and we showed that this procedure, implementing unitarity in coupled channels, produced results in remarkable agreement with those of the analysis of [10], where a Λ(1405) resonance Breit-Wigner term (accounting for Flatté effect) was introduced. We have shown that this agreement is not trivial or general, but occurs in the present case. In this Appendix we explicitly evaluate the contribution of the different mechanisms to the total scattering amplitude in Eq. 6 which depend on the different partial waves and spin possibilities of the different resonances considered.
For the different spin and angular momentum we will follow the nomenclature of fig. 10, i.e., s, L and J P are the spin, orbital angular momentum and total spin-parity respectively of the J/ψp pair and L ′ and J ′P ′ stand for the K − orbital angular momentum and total spin-parity of the pK − system.
If we assume that the dominant contribution to a given partial wave is given by the smallest allowed orbital angular momentum, we can write the relevant terms of the Λ b → J/ψK − p amplitude of interest in the present work as The term S 1 2 − term is given by S 1 2 − = σ · ǫ where ǫ stands for the J/ψ polarization vector. This operator projects over J/ψp and K − s-waves [53]. The structures P 1 2 − and P 3 2 − are obtained in ref. [54] by imposing orthonormality to S 1 2 − and are given by [55]: where k is the K − momentum. Note that the Λ(1600) and Λ(1810) contribute both to P 3 2 − and P 1 2 − . Since the first Λ b J/ψΛ vertex is s-wave and the ΛpK − vertex is pwave, the actual structure for the Λ(1600) and Λ(1810) I: Quantum numbers of the different contributions to the Λ b → J/ψK − p amplitude with the lowest possible value of orbital angular momentum for a given value of J P and the J ′P ′ used in the present work. cases is m p | σ · k σ · ǫ|m Λ b which in terms of P 3 2 − and P 1 2 − reads Although not needed in the present work, the J ′P ′ = 3 2 + combination is given by with S † the spin 1/2 to 3/2 transition operator. On the other hand, D 3 2 − is given by since it accounts for the K − in d-wave and J/ψp in swave. Finally, the expression for the D 5 2 + is where p is the J/ψ momentum. Note that, while the (k i k j − k 2 δ ij /3) term in Eq. (17) is purely d-wave, the i( σ × ǫ) i p j term gives contribution to both This is the reason why we consider Eq.(17) to account for the J P = 5 2 + case.
In the previous equations J P A and J P B stand for the spin-parity of the P c (4380) and P c (4450) respectively. The Kronecker deltas are introduced to account for the different possibilities of quantum numbers of the two pentaquarks which, according to their J P , contribute to the corresponding partial wave amplitude. For instance, if one wanted to consider the (J P A , J P B ) = (1/2 − , 5/2 + ) case then the P c (4380) propagator would contribute to the a and c coefficient, while the P c (4380) propagator would contribute only to the f coefficient.
Note the extra inclusion of contact terms (the 1 addends) in b, c and f coefficients to account for possible contact terms with those quantum numbers.
The α i and C i coefficients are complex in general. Therefore, taking into account that there are some unobservable arbitrary global phases, there are 19 free parameters to be fitted in the general case.