A hidden-charm $S=-1$ pentaquark from the decay of $\Lambda_b$ into $J/\psi \eta \Lambda$

The hidden charm pentaquark state $P_c(4450)$ observed recently by the LHCb collaboration may be of molecular nature, as advocated by some unitary approaches that also predict pentaquark partners in the strangeness $S=-1$ sector. In this work we argue that a hidden-charm strange pentaquark could be seen from the decay of the $\Lambda_b$, just as in the case of the non-strange $P_c(4450)$, but looking into the $J/\psi \eta \Lambda$ decay mode and forming the invariant mass spectrum of $J/\psi\Lambda$ pairs. In the model presented here, which assumes a standard weak decay topology and incorporates the hadronization process and final state interaction effects, we find the $J/\psi \eta \Lambda$ final states to be populated with similar strength as the $J/\psi K^- p$ states employed for the observation of the non-strange pentaquark. This makes the $\Lambda_ b \to J/\psi \eta \Lambda$ decay to be an interesting process to observe a possible strange partner of the $P_c(4450)$ . We study the dependence of the $J/\psi \Lambda$ mass spectra on various model ingredients and on the unknown properties of the strange pentaquark.


Introduction
The LHCb collaboration reported recently two exotic structures in the invariant J/ψp mass spectrum of the Λ 0 b → J/ψK − p process. These pentaquark states were named P c (4380), with a mass of 4380 ± 8 ± 29 MeV and a width of 205±18±86 MeV, and P c (4450), with a mass a e-mail: feijoo@icc.ub.edu b e-mail: vladimir@icc.ub.edu c e-mail: angels.ramos@icc.ub.edu d e-mail: Eulogio.Oset@ific.uv.es of 4449.8 ± 1.7 ± 2.5 MeV and a width of 39 ± 5 ± 19 MeV [1,2]. Hidden charm baryon states with similar characteristics of the states reported had already been predicted, employing a molecular picture [3,4,5,6,7] or a quark model approach [8,9]. A list of early references on pentaquark states can be seen in Ref. [10]. The CERN discovery triggered a large number of theoretical works trying to give an explanation for the two reported states. The molecular picture was invoked in [11,12,13,14], the diquark picture in [15,16,17,18,19], QCD sum rules were used in [20,21], and the soliton model was employed in [22]. It has also been argued that the observed enhancement could be due to kinematical effects or triangular singularities [23,24,25]. Suggestions of different reactions to observe the pentaquarks have been reported [26,27,28,29,30], while explicit decay modes to elucidate their structure have also been studied in [31,32]. Further discussions on the issue and the nature of the two P c states can be seen in Refs. [33,34] and particularly in the recent detailed review of Ref. [35].
In what respects the present work, we recall that a theoretical study of the Λ 0 b → J/ψK − p reaction was done in [36], prior to the experimental study of Ref. [1], predicting the contribution of the tail of the Λ(1405) in the K − p invariant mass distribution. The analysis of [1] contains such a contribution in agreement in shape with the predictions, where the absolute normalization is unknown. Moreover, it was shown in [12] that the distributions in the pentaquark channel, i.e. in the invariant J/ψp mass spectrum of [1], could be explained via the incorporation of the hidden charm N * states predicted in [3,4,5,6], which are molecular states mostly made fromD * Σ c orD * Σ * c components and having a small coupling to J/ψp, one of their open decay channels. It is unlikely that there are no partners of the states found in [1], and indeed, in [3,4] states of spin-parity 3/2 − with hidden charm but strangeness S = −1 were predicted, mostly made ofD * Ξ c orD * Ξ c , decaying into J/ψΛ. In view of this, the decay of Ξ − b into J/ψK − Λ was suggested in [37] as a suitable reaction to find a hidden charm strange state. Predictions for theKΛ and J/ψΛ mass distributions were done, and, playing with uncertainties, it was shown that a clear peak in the J/ψΛ mass distribution should show up. This reaction is presently being considered by the LHCb collaboration. However, since there is a much smaller statistics in the production of Ξ − b than that of Λ b [38] , it is interesting to explore alternative reactions to observe this strangeness S = −1 hidden-charm pentaquark. In the present paper we suggest to employ the J/ψ ηΛ decay mode of the Λ b . Since the ηΛ pair is populated with a weight √ 2/3 relative to the K − p pair in the primary [36], the Λ b → J/ψ ηΛ decay rate should be similar as that found for J/ψK − p final states in the study of the non-strange pentaquark, and the new strange state should be looked for in the J/ψΛ mass distribution instead of the J/ψp one. We note that the possible existence of an strange S = −1 pentaquark partner was also studied in [39] from the non-strange decay mode Λ b → J/ψK 0 Λ, which is one of the coupled channels of the decay Λ b → J/ψπ − p from which, even if it is more Cabbibo suppressed than the Λ b → J/ψK − p process, a possible signal of the P c (4450) may also have been seen [2,40]. The study of [39] explored the effect of different weak decay amplitudes to produce either a J P = 1/2 − or a J P = 3/2 − strange pentaquark. In this work, we will also take these possibilities into account. This paper is organized as follows. In Sect. 2 we present our formalism for the Λ b → J/ψ ηΛ decay, describing the weak transition process and the implementation of final state interactions. Our results are shown in Sect. 3, where the spectra of both ηΛ and J/ψΛ states can be seen and a discussion on their dependence on various uncertain parameters of our model can be found. We shall show that, even within uncertainties, the signal for a strangeness S = −1 hidden charm pentaquark remains as a clear peak in the J/ψΛ mass distribution. Our conclusions are summarized in Sect. 4.

Formalism
The study of the Λ b → J/ψ ηΛ decay follows the same approach as that presented in [36] for Λ b → J/ψ K − p. At quark level, both processes are identical and proceed through the transition diagram depicted in Fig. 1, where we can see the W -exchange weak process transforming the b quark into ccs, followed by the hadronization of a pair of quarks which eventually produces a meson and a baryon, in addition to the J/ψ. The process depicted in Fig. 1 assumes that the elementary weak decay involves only the b quark of the Λ b , while the u and d quarks remain as spectators, the reason being that one expects one-body operators in a microscopical evaluation to have larger strength than twoor multi-body operators. According to this assumption, since the Λ b has isospin I = 0, so does the spectator ud pair, which, combined with the s quark after the weak decay, can only form I = 0 Λ states. The findings of the experimental analysis of Ref. [1] clearly support this hypothesis.
For the hadronization process we introduce aqq pair between two quarks with the quantum numbers of the vacuum,ūu +dd +ss. The dominant contribution of the hadronization preserves the spectator role of the ud pair, which ends up into the final baryon, and requires the involvement of the s quark, which ends up into the final meson. Any other topology that would bring the u or d quark into the final meson requires a large momentum transfer that supresses the mechanism. If, in addition, we wish to have the meson-baryon pair in s-wave, it will have negative parity, forcing the s quark prior to hadronization to have also this parity and thus be in an excited state. Since in the final K − or η mesons the s quark is in its ground state, this also implies that the s quark produced immediately after the weak process must participate actively in the process of hadronization, which proceeds as shown in Fig. 1. A further discussion on the reduced size of other alternative mechanisms can be found in [41].  The technical way to implement the hadronization and produce meson-baryon pairs in the final state follows the same steps as in [42,43,44] for meson decays and in [36,45] for the Λ b decay. The flavor decomposition of the Λ b state is: which becomes, after the weak process or, upon hadronization, which can be written in terms of the qq matrix P , as where P =   uū ud us dū dd ds sū sd ss Writing the matrix P in terms of the meson states, P → φ, where the η, η mixing [46] has been assumed, the hadronized state becomes: By the former equation one obtains the mixed antisymmetric representation of the octet of baryons and taking the results of [47] (see also [36]) one finds the final representation where we have omitted the |η Λ contribution because of the large mass of the η meson [36]. The final step consists in taking into account the final state interaction of the meson-baryon pairs. The amplitude for the Λ b → J/ψ ηΛ decay will then be built from the diagrams of Fig. 2, where we can see the direct tree-level process, depicted by diagram (a), the finalstate interaction contribution of the meson-baryon pair into ηΛ production (b), and the final-state J/ψΛ → J/ψΛ interaction (c). The corresponding amplitude can be written as: Fig. 2 Diagrammatic representation of the decay amplitude for Λ b → J/ψ ηΛ: a) tree level, b) the ηΛ production through the coupled channel interaction of the initially produced ηΛ andKN meson-baryon pairs, c) J/ψΛ → J/ψΛ interaction.
where the weights h i , obtained from Eq. (8), are: and where G i , with i = K − p,K 0 n, ηΛ, denotes the meson-baryon loop function, chosen in accordance with the model for the scattering matrix t i,ηΛ [48]. Similarly, we take the loop function G J/ψΛ employed in the model of [3,4] on which, as we will show below, we base our prescription for t J/ψΛ,J/ψΛ . The factor V p , which includes the common dynamics of the production of the different pairs, is unknown and we take it as constant, see Ref. [45] for a more detailed argumentation. At this point it is worth mentioning that the model for the final state interaction in the ηΛ channel, which is briefly described in the next section, generates some resonances dynamically, like the Λ(1405) or the Λ(1670), that are either below or at the edge of the threshold of ηΛ invariant masses M ηΛ accessible from the decay of the Λ b . We therefore would like to explore the possibility of adding to the amplitude the explicit contribution of some Λ * which lies in the relevant M ηΛ region, essentially between 1700 MeV and 2500 MeV, and might couple sensibly to ηΛ states, as represented diagrammatically in Fig. 3. One state with these characteristics is listed in the PDG compilation [49], the one star Λ(2000), having a width Γ ∼ 100 − 300 MeV and a branching ratio to the decay into ηΛ of (16 ± 7)%. The recent unitary multichannel model forKN scattering, with parameters fitted to partial waves up to J = 7/2 and up to 2.15 GeV of energy, also gives an s-wave J P = 1/2 − Λ state with similar mass and width properties [50].
We note that, since our model will rely on the strange pentaquark predicted in Refs. [3,4] at an energy around 4450 MeV, which couples strongly toD * 0 Ξ c states, one Fig. 3 Diagrammatic representation of an s-wave resonance contribution to the Λ b → J/ψ ηΛ decay amplitude should consider the possibility that the influence of this resonance in the final J/ψΛ mass distribution could also be due to the creation of a virtualD * ηΞ c state in a first step of the Λ b decay, through the mechanism of Fig. 4, followed by multiple interactions to generate the resonance, which would eventually decay into a J/ψΛ pair in the final state, represented by the diagrams of Fig. 5. However, this configuration requires a different topology, as seen in Fig. 4, in which the ud quarks of the Λ b do not act as a coupled spectator pair. Although it is hard to quantify the size of the amplitude of Fig. 1 with respect to that of Fig. 4, the fact that in this later case one of the spectator quarks ends up in the charmed meson and the other one goes to the baryon makes us believe that the corresponding amplitude will be reduced. We will therefore assume the dominance of the mechanism of Fig. 1 over that of Fig. 4 by a factor of two or more and will give predictions for different relative signs of the two processes. We note that the lowest order contribution to the Λ b → J/ψ ηΛ decay induced by virtualD * ηΞ c states is the amplitude of Fig. 5 (a). We have checked, by explicit numerical evaluation, that the next-order contribution of Fig. 5 (b), involving the additional final state interaction of the ηΛ pair, gives a negligible correction, hence it will be ignored in the results presented in Sect. 3. Adding the s-wave resonant contribution of Fig. 3 and the process initiated by an intermediateD * ηΞ c state followed by final state interactions leading to a J/ψΛ pair and an η meson represented by the diagram of Fig. 5 (a), the final amplitude for Λ b → J/ψ ηΛ de- cay, producing a strange pentaquark with J P = 1/2 − becomes: where α is a dimensionless parameter that determines the strength of the s-wave resonant mechanism, while β controls the strength of Λ b decaying virtually intō D * ηΞ c , relative to its decay into J/ψ ηΛ. The amplitudes of Eqs. (9) and (12) come with the matrix element m Λ | σ | m Λ b , tied to the s-wave character assumed for the weak decay vertex and accounting for the spin 1/2 of the decaying Λ b , and the spins 1/2 and 1 of the emitted Λ and J/ψ meson, respectively. Moreover, as will be recalled in the next subsections, the J/ψΛ (and ηΛ) interaction models are also taken in s-wave, hence the spin values of the J/ψΛ pair could in principle be J = 1/2 or J = 3/2. However, only the J = 1/2 case is allowed to match the spin 1/2 of the decaying Λ b . In fact the σ operator projects the J/ψΛ system into J = 1/2 as shown in [51].
Our general strategy is to assume that the decay process proceeds involving the smallest possible angular momentum at the vertices. Therefore, in order to produce the strange pentaquark with J = 3/2 it is necessary to implement at least a p-wave contribution in the weak decay mechanism. A p-wave operator of the form was considered in [39], where it was also shown that, if the peak seen in [1] corresponds to the molecular states with J P = 3/2 − generated in [3,4] from the scalarvector meson interaction in s-wave, the momentum q in the former equation must be that of the η-meson, which is taken in the rest frame of the ηΛ system. A decomposition of the p-wave vertex in terms of two operators, that project, respectively, over the spin J = 3/2 and J = 1/2 of the two-body J/ψΛ system, was also given in [39]: The J = 3/2 pentaquark will then be generated by the final state interaction of the J/ψΛ pair initiated by the p-wave decay vertex of Eq. (13), in a process of the type represented by the diagram of Fig. 2(c), or also from the virtual excitation of intermediateD * ηΞ c states followed by multiple interactions leading to J/ψΛ pairs in the final state, as seen in Fig. 5(a). Note that in this later case the p-wave decay vertex will have the same structure as that of Eq. (13), but with a strength constant B , and will not act at tree-level.
Considering the s-wave and p-wave contributions, the amplitude that allows for the appearance of a pentaquark with J = 3/2 is then given by: where the term proportional to V p stands for the contribution of the s-wave weak decay amplitude 1 , the next row contains the contribution of the p-wave tree level term and the last two rows correspond to the final state interaction contributions that generate the J = 3/2 pentaquark initiated by J/ψΛ states (proportional to B) or byD * Ξ c states (proportional to B ). The former equation can be cast schematically as: Finally, the double differential cross-section for the Λ b → J/ψ ηΛ decay process reads [36]: where, after performing the sum over final spins and polarizations and the average over initial spins (see appendix in [39]), one has: with M being that of Eqs. (9) or (12), corresponding to an s-wave weak vertex and, hence, producing a pentaquark with J = 1/2, or with M being that of Eq. (16), corresponding to a weak vertex that also has a p-wave term and, hence, making the production of a pentaquark with J = 3/2 possible. Fixing the invariant mass M J/ψΛ , one can integrate over M ηΛ in order to obtain dΓ/dM J/ψΛ . In this case, the limits are given by: and where Similar formulas are obtained when fixing the invariant mass M ηΛ and integrating over M J/ψΛ to obtain dΓ/dM ηΛ : where

Final interaction models
In this section we briefly describe the theoretical models employed to obtain the amplitudes t i,ηΛ , t J/ψΛ,J/ψΛ and tD * Ξ c ,J/ψΛ , which account for the final state interaction effects. The S = −1 meson-baryon amplitude with ηΛ in the final state appearing in diagram (b) of Fig. 2 is determined from the coupled-channel unitary model of Ref. [48], developed with the aim of improving upon the knowledge of the chiral interation at next-to-leading order (NLO). The parameters of the model were fitted to a large set of experimental scattering data [52], as well as to branching ratios at threshold [53], and to the precise SIDDHARTA value of the energy shift and width of kaonic hidrogen [54]. Differently to other works, as e.g. [55,56,57], the model was also constrained to reproduce the K − p → K + Ξ − , K 0 Ξ 0 reactions, since they are especially sensitive to the NLO terms. The work of Ref. [48] also investigated the influence of high spin hyperon resonances on the K − p → KΞ amplitudes, finding that the resonant terms helped in improving the description of the scattering data and produced more precise values of the low energy constants of the chiral unitary model.
More especifically, the meson-baryon amplitudes of Ref. [48] are built from a kernel obtained from the SU(3) chiral Lagrangian up to NLO: where and with The indices i, j stand for any of the ten meson-baryon channels in the neutral S = −1 sector: K − p,K 0 n, π 0 Λ, π 0 Σ 0 , π − Σ + , π + Σ − , ηΛ, ηΣ 0 , K + Ξ − and K 0 Ξ 0 , while M i , M j and E i , E j are the masses and energies, respectively, of the baryons involved in the transition, and k i,µ , k µ j are the corresponding meson four-momenta. The lagrangian is written in terms of SU(3) coefficients C ij , the pion decay constant f and other low energy constants embedded in the matrices D ij and L ij of the NLO term, which can be found, for example, in the appendices of Ref. [48]. Chiral unitary amplitudes are obtained by solving the Bethe-Salpeter equation in its on-shell factorized form: where the meson-baryon loop function G l is obtained employing dimensional regularization where M l and m l are the baryon and meson masses of the "l" channel, the regularization scale µ is taken to be 1 GeV , and a l are subtraction constants, which, together with the low energy parameters of the lagrangian, were determined from fits to data performed in Ref. [48]. We will employ the set of parameters corresponding to the model named "NLO*" there, which will be referred to as Model 1 here. These results will be compared to those obtained with another set of parameters -named "WT (no KΞ)" in Ref. [48] and referred to as WT in the present paper-obtained from a fit that employs the lowest order Weinberg-Tomozawa (WT) term without taking into account the experimental data corresponding to the KΞ channels, as most of the works in this field. The chiral approach was complemented with the explicit inclusion, in the K − p → K 0 Ξ 0 , K + Ξ − amplitudes, of two high spin resonances, Σ(2030) and Σ(2250), selected from the possible candidates listed in the PDG [49] and in accordance to other resonance-based models [58,59]. The spin and parity J π = 7/2 + of the Σ(2030) are well established. Those of the Σ(2250) are not known, but the choice J π = 5/2 − was adopted out of the two most probable assignments, 5/2 − or 9/2 − . The fit of the model that includes the resonances, named "NLO+RES" in [48] and Model 2 here, determines not only the low energy parameters and subtraction constants but also the couplings, masses, widths and form-factor cut-offs of the resonances. More details on the implementation of the resonant terms can be found in Ref. [48].
With respect to the final state interaction in the J/ψΛ sector, represented by the diagrams of Fig. 2(c) and Fig. 5(a), we recall that two states with strangeness and hidden charm with J P = 3/2 − and I = 0 were found in Refs. [3,4] as meson-baryon molecules, having pole positions √ s = 4368 − 2.8i and √ s = 4547 − 6.4i and coupling to J/ψΛ states with strength | g J/ψΛ |= 0.47 and 0.61, respectively. The magnitude of each of these couplings is relatively small compared to the coupling of the pole to the main meson-baryon component, which for the lower energy pole isD * Ξ c , with | gD * Ξc |= 3.6, while for the higher energy one isD * Ξ c , with | gD * Ξ c |= 2.6. In any case, | g J/ψΛ | is large enough to create a peak in the mass distribution, as we shall see. As candidate for the strangeness −1 pentaquark, we will consider the state at higher energy since its mass is close to the non-strange pentaquark found in [1]. One must however accept that the mass obtained for these states has uncertainties since, unlike in other sectors, one does not have any experimental data to constrain the parameters of the theory. We therefore take the nominal value of about M R = 4550 MeV for the mass of the strange pentaquark and will explore the stability of our results to variations of this mass. We shall take Γ R = 10 MeV in agreement with the findings of [3,4]. Our explorations are implemented employing the following Breit-Wigner representation for the t J/ψΛ,J/ψΛ and tD * Ξ c ,J/ψΛ amplitudes Then the production of the resonance is done through the J/ψΛ → J/ψΛ andD * Ξ c → J/ψΛ amplitudes, parametrized through the expressions given in Eqs. (35) and (36), as seen in diagrams of Fig. 2(c) and Fig. 5(a), respectively, as well as in Eqs. (9), (12) or (16). The values of the couplings are g J/ψΛ = −0.61 − 0.06i and gD * Ξ c = 2.61 − 0.13i. The loop functions G J/ψΛ and GD * Ξ c appearing in these equations are taken from [3,4], where a dimensional regularization method with a scale µ = 1000 MeV was employed, using subtraction constants a J/ψΛ = aD * Ξ c = −2.3.

Results
We start this section by presenting, in Fig. 6, the invariant mass distributions of J/ψΛ states produced in the decay Λ b → J/ψ ηΛ, obtained from the simpler s-wave weak decay approach of Eq. (9) and for three different models of the S = −1 ηΛ interaction [48]: one that only considers the lowest-order WT term of the Lagrangian (dotted line) and two other models, Model 1 (dashed line) and Model 2 (solid line), that incorporate the next-to-leading order terms and, in the case of the later one, the effect of higher spin resonances. Please note that although the Model 2 includes the additional contribution of the two high-spin resonances, they do not contribute directly to the studied I = 0 decay due to their I = 1 nature, but their inclusion does modify the parameters of Model 2 with respect to those of Model 1. For the three models of the ηΛ interaction, the peak of the pentaquark is clearly observed at 4550 MeV, the value of the mass M R employed in the parametrization of Eqs. (35) and (36). However the overall strength is enhanced for the NLO models, which also show a different interference pattern with the non-resonant background to that of the lowest-order WT model.
The invariant mass distribution of ηΛ pairs is shown in Fig. 7, where the J/ψΛ resonant structure has disappeared since the invariant M J/ψΛ masses have been integrated out. The ηΛ invariant mass distributions have essentially the same shape as that of the distributions shown already in Fig. 6 of Ref. [45], which did not consider the additional contribution associated to the hidden charm strange pentaquark. We see a broad peaked shape, associated to the NLO terms of the Lagrangian, which is more pronounced in the case of Model 1. As already discussed in Ref. [45], this structure is not associated to any resonant state since it appears at different energies in different channels and no pole in the complex plane was found either.
In the following, results will be presented for only one model of the strong ηΛ interaction, chosen to be Model 2 as it provides a better account of the scattering observables [48]. The J/ψΛ invariant mass distributions displayed in Figs. 8 and 9, for different values of the pentaquark coupling to J/ψΛ and for different values of the pentaquark mass, respectively, show obvious trends. From Fig. 8 we can conclude that the pentaquark could be seen over the background even if its coupling to J/ψΛ states were as low as | g J/ψΛ |= 0.48. The unitary approaches of Refs. [3,4,5,6] predict values for this coupling in between 0.5 − 1.0, which make us believe that the strange pentaquark could leave a clear signature in the J/ψΛ mass spectrum. mass M R do not practically change the aspect of the ηΛ invariant mass spectrum. In Fig. 10 we explore the effect of including the additional effect of a Λ(2000) s-wave resonance coupling to ηΛ states. The unknown coupling strength α of Eq. (12) is varied such that it produces a clearly visible change in the spectrum of ηΛ invariant masses over what we obtain in the absence of this contribution, as seen in the bottom panel of Fig. 10. In the top panel we observe that the inclusion of the Λ(2000) on the J/ψΛ pair distribution, where the ηΛ invariant masses have been integrated out, essentially enhances the strength while keeping the same shape for the different values of α.
We next explore the influence of the strange pentaquark being initiated by the excitation of a virtual D * 0 ηΞ c state, followed by the multiple scattering of D * 0 Ξ c leading to a final J/ψΛ pair and an η meson. As discussed in Sect. 2.1, the topology for this decay should lead to a reduced amplitude with respect to that of the J/ψ ηΛ case. We implement this phenomenologically through the parameter β, as seen in Eq. (12), which is given the values −0.5, −0.25, 0.0, 0.25, and 0.5 accounting also for different relative sign cases. The results obtained with the negative values are displayed in Fig. 11 and those with the positive values in Fig. 12. As seen in the top panel of Fig. 11, the pentaquark signal for the negative values of β gets somewhat reduced with respect to the case in which the virtual excitation ofD * 0 ηΞ c states is omitted, indicating a destructive interference with the direct excitation of J/ψ ηΛ states. However, the signal is still clearly visible over the background. The situation is completely different for the positive values of β. As seen on the top panel of Fig. 12, the signal of the pentaquark is tremendously enhanced due to a constructive interference between both mechanisms and to the fact that the coupling strength of the pentaquark toD * 0 Ξ c states is a factor 4 larger than that to J/ψΛ. In both positive and negative β cases, the changes seen in the ηΛ invariant mass spectra are relatively minor, as can be seen in the bottom panels of Figs. 11 and 12.
Up to here, we have been discussing the results assuming the pentaquark to have J P = 1/2 − , which can then be produced by an s-wave mechanism for the Λ b decay. In the case of a J P = 3/2 − pentaquark, which is another of the possibilities for the states predicted in [3,4], it is necessary to implement at least an additional p-wave contribution, as that of Eq. (13). Our results for this case are displayed in Fig. 13. The dotted line represents the case in which only the s-wave contribution is kept, producing a J/ψΛ pair in 1/2 − . Since the pentaquark is now assumed to have J P = 3/2 − , it does not show in that J/ψΛ invariant mass spectrum, which reduces to a structureless background. We could have included, as in the study of the J = 1/2 pentaquark case, the Λ(2000) s-wave resonance contribution, but we have omitted this effect in the present J = 3/2 discussion because, although it would be seen as an additional structure in the ηΛ invariant mass distribution, it would simply contribute with a practically constant strength to the spectrum of J/ψΛ pairs, similarly to what we have found in Fig. 10.
When we add the p-wave vertex of Eq. (13), we obtain the distributions displayed by the dashed curves in Fig. 13. The size of the coupling constant, B = 0.001 MeV −1 , has been chosen so that the p-wave contribution has a visible effect over the s-wave J/ψΛ and ηΛ invariant mass distributions. The J/ψΛ spectrum, shown in the top panel of Fig. 13, presents a dip at the pentaquak mass, which comes from the interference between the tree level and the J/ψΛ final state interaction terms, displayed by Figs. 2(a) and (c), respectively, as can also be seen in Eq. (16). This is the same behavior as that observed in the study of the strange pentaquark from the Λ b → J/ψK 0 Λ decay in [39]. In the present work, we also incorporate the excitation of the pentaquark from the multiple scattering ofD * 0 Ξ c pairs produced in the virtual Λ b →D * 0 ηΞ c decay, which proceeds also in p-wave with a strength B . This is a necessary consideration if one wants to interpret the pentaquark as the state emerging from the interaction of D * 0 Ξ c and its related coupled states. If we now assume a ratio between the Λ b →D * 0 ηΞ c and Λ b → J/ψ ηΛ amplitudes of B /B = 0.5, we obtain the solid curve, where the dip has turned into a wiggled shape. When the sign of B is opposed to that of B, we find a similar behavior, although in a reflected way, as depicted by the dot-dashed curve. In either case, a visible pentaquark signal is obtained, as can be seen more clearly in the inset of Fig. 13.

Conclusions
The recent finding of two structures in the J/ψp invariant mass distribution of the Λ b → J/ψK − p decay, associated to two pentaquark states, together with its plausible explanation in terms of a previously predicted hidden charm baryon molecular state, prompted us to study the decay of the Λ b into J/ψ ηΛ final states. The Λ b → J/ψ ηΛ decay, being a coupled channel of the Λ b → J/ψK − p one, will occur with similar strength and one could observe, in the J/ψΛ invariant mass spectrum, possible strange partners of the two non-strange pentaquark states reported by the LHCb collaboration.
We recall that when the hidden charm N * resonances were theoretically predicted as molecular states in several unitary approaches, some partner hidden charm strange Λ * states were also found. We have taken advantage of this finding and have predicted what signal should one of these states leave in the ηΛ and J/ψΛ invariant mass distributions of the Λ b → J/ψ ηΛ reaction. We have found that, taking the values of the couplings of the hidden charm Λ * state to theD * 0 Ξ c and J/ψΛ channels obtained in the unitary approaches, one should observe clear and sizable peaks in the J/ψΛ mass distribution of the Λ b → J/ψ ηΛ decay. We have also used for this purpose the interaction of ηΛ with its coupled channels by means of a recent chiral unitary approach at next to leading order.
We have studied the dependence of our results on reasonable changes in the parameters of the models involved in our description of the process, as well as on the unknown properties of the speculated hidden charm strange pentaquark. We have observed that, while there appear changes in the position of the peak and in the shapes of the distributions, a resonance signal in the J/ψΛ invariant mass spectrum is clearly seen in all the cases. This gives us confidence that such an experiment should result into a successful proof of the existence of this new state and encourage the experimental analysis of this decay channel, for which our theoretical study predicts a similar strength than for the Λ b → J/ψK − p reaction already analyzed by LHCb.