The ¯ K N Interaction in Higher Partial Waves

We present a chiral ¯ KN interaction model that has been developed and optimized in order to account for the experimental data of inelastic ¯ KN reaction channels that open at higher energies. In particular, we study the effect of the higher partial waves, which originate directly from the chiral Lagrangian, as they could supersede the role of high-spin resonances employed in earlier phenomenological models to describe meson-baryon cross sections in the 2 GeV region. We present a detailed derivation of the partial wave amplitudes that emerge from the chiral SU(3) meson-baryon Lagrangian up to the d-waves and next-to-leading order in the chiral expansion. We implement a nonperturbative unitarization in coupled channels and optimize the model parameters to a large pool of experimental data in the relevant energy range where these new contributions are expected to be important. The obtained results are encouraging. They indicate the ability of the chiral higher partial waves to extend the description of the scattering data to higher energies and to account for structures in the reaction cross-sections that cannot be accommodated by theoretical models limited to the s-waves.


Introduction
Over the last forty years, strong interaction processes in the low-energy domain have been described with great success by chiral perturbation theory (χPT), an effective field theory with hadron degrees of freedom, which respects all the symmetries of the underlying theory of the strong interaction, quantum chromodynamics (QCD), particularly its spontaneously broken chiral symmetry [1].
In the early years of χPT many studies were focused on the low-energy pion interactions with pions and nuclei, although extensions to the strangeness sector, invoking flavor SU(3) symmetry, were also explored. In particular, theKN interaction was especially challenging because the description of the scattering amplitude cannot be properly achieved within a perturbative treatment, owing to the presence of a resonance, the Λ(1405), located 27 MeV below theKN threshold.
Another challenge in hadron physics over the past few decades has been to disentangle the nature of some hadron resonances that do not fit well in a conventional quark model description, according to which baryons are composed by three quarks and mesons by a quark-antiquark pair [2]. For many years, a large amount of theoretical and experimental activity has been devoted to the discovery and characterization of exotic mesons and baryons. Early evidence was provided, precisely, by the Λ(1405) resonance, whose mass was systematically predicted to be too high by quark models, but it found a better explanation if it was described as aKN quasi-bound state. This was already pointed out in the late fifties [3,4] and was later corroborated by models that built the meson-baryon interaction from a chiral effective Lagrangian, derived respecting the symmetries of QCD, In the present work we consider the s-and p-wave contributions of the scattering amplitudes but, differently to Reference [40], we also employ the data of inelastic reactions opening at higher energies to constrain our models. This implies that the NLO Lagrangian, which we take from Reference [37], must also consider terms that have usually been disregarded owing to their negligible effect at lower energies. We will show that the new terms of the NLO Lagrangian are relevant, being in fact strongly intertwined with the p-wave components of theKN interaction.
Our paper is organized as follows. In Section 2 we develop the chiral unitary formalism, show the LO and NLO contributions of the chiral Lagrangian, give the explicit expressions of the corresponding interaction kernels-with their complete momentum dependence structure-and, finally, describe the procedure that allows one to obtain the partial-wave components of the scattering amplitude. The models studied in this work are presented in Section 3, together with the data employed in the fits. Section 4 is devoted to the discussion of our results, with a special focus on the role played by the p-wave amplitudes and the new terms of the NLO Lagrangian. A few concluding remarks are given in Section 5.

Formalism
Chiral unitary approaches (UChPT) have shown to be a powerful tool to treat the meson-baryon scattering at energies around resonances. These nonperturbative schemes prevent plain chiral perturbation theory from non converging and guarantee the unitarity and analyticity of the scattering amplitude. In the present work, unitarity is implemented by solving the Bethe-Salpeter (BS) equation with coupled channels. Following [6,41], the interaction kernel is conveniently split into its on-shell contribution and the corresponding off-shell one. The off-shell part gives rise to a tadpole-type diagram which can be reabsorbed into renormalization of couplings and masses and, hence, be omitted from the calculation. This method permits factorizing the interaction kernel and the scattering amplitude out of the integral equation, transforming a complex system of coupled integral equations into a simple system of algebraic equations that, in matrix form, reads where V ij is the driving kernel derived from the chiral Lagrangian, T ij is the corresponding scattering amplitude for the transition from an i channel to a j one, and G l is the loop function of the intermediate channel l, which reads with M l and m l being the baryon and meson masses of the channel. As this function diverges logarithmically, a dimensional regularization scheme is applied that leads to the expression: The subtraction constants a l replace the divergence for a given dimensional regularization scale µ, which is taken to be 1 GeV in the present work. These constants are unknown parameters to be fitted to the experimental data. For the S = −1 meson-baryon interaction they amount to ten, one per channel, although it is quite common to reduce them to six taking into account isospin symmetry arguments, as we shall also do.
The kernel employed in the chiral unitary approaches is derived from the SU(3) effective chiral Lagrangian, which provides the fundamental blocks of the interaction that preserve the symmetries of QCD employing hadron fields as the relevant degrees of freedom. Following a power counting scheme, these constituent pieces are arranged in the expansion by order of relevance (see for instance [42] for a more detailed explanation). At leading order (LO), the most general Lagrangian can be expressed as: where M 0 is the common baryon octet mass in the chiral limit, the constants D, F denote the axial vector couplings of the baryons to the mesons, and the symbol · stands for the trace in flavor space. The baryon octet field (N, Λ, Σ, Ξ) is denoted by B, while the pseudoscalar meson octet field φ (π, K, η) enters in [28]) with f being the pseudoscalar decay constant that acts as typifying scale factor in the expansion in powers of momentum. Finally, [D µ , B] contains the covariant derivative that accounts for the local character of the chiral transformation of u, and it is defined as: with Γ µ = [u † , ∂ µ u]/2 being the chiral connection. The next-to-leading order (NLO) contributions are given by: where M N stands for the nucleon mass, while the coefficients b D , b F , b 0 , d i (i = 1, . . . , 4), g i (i = 1, . . . , 4) and h i (i = 1, . . . , 4) are the low-energy constants (LECs) at this order. Despite the symmetries of the underlying theory cannot fix these constants, some of them can be constrained by the mass splitting of baryons, the pion-nucleon sigma term and from the strangeness content of the proton [43], and some others using data from low-energy meson-baryon interactions, such as the isospin even πN s-wave scattering length [44] and the isospin zero kaon-nucleon s-wave scattering length [45]. However, we treat them as free parameters in the fitting procedure, as it is usually done in the literature. The reason lies in the coupled-channel unitarization that produces amplitudes that go beyond tree level. The quantity χ + = 2B 0 (u † Mu † + uMu) explicitly breaks chiral symmetry via the quark mass matrix M = diag(m u , m d , m s ), while B 0 = − 0|qq|0 / f 2 relates to the order parameter of spontaneously broken chiral symmetry. Before proceeding further, we would like to make a few remarks on the NLO Lagrangian of Equation (6), taken from [37], where the authors study the KN scattering amplitude in a chiral unitary approach looking for a possible broad resonance in the S = +1 sector. Formally speaking, the expression differs from the one we employed in our previous works [28][29][30] in the terms with coefficients g i (i = 1, . . . , 4) and h i (i = 1, . . . , 4) (from now on g-terms and h-terms), which naturally appear from an extension of the SU(2) chiral Lagrangian introduced in [8,46,47], but have been commonly discarded in the literature studying S = −1 meson-baryon scattering within the UChPT framework owing to their nonsignificant role at low energies. Nevertheless, as our models explore higher energies, it is reasonable to consider these additional contributions. We will show that they play a relevant role mainly in theKN transitions to the ηΛ, ηΣ 0 , and KΞ channels. We also note that the b D , b F , b 0 , d i and h-terms are invariant under the hermitian conjugation (h.c.) transformation in Equation (6); hence, their h.c. contributions can be reabsorbed in the LECs of the corresponding original terms. In contrast, the h.c. of the g-terms produces new structures in the Lagrangian, which have to be considered explicitly. Specially remarkable is the cancellation of the g 3 monomial, which implies the reduction in the parameter space in the derived models. Finally, it should be pointed out that other versions for SU(3) Lagrangians at O(q 2 ) are presented in [48,49].
Not without a tedious calculation, one can derive the amplitudes for the φ i B s i → φ j B s j processes (with the incoming and outgoing spins s, s , respectively) which are schematically represented in Figure 1. Diagram (i) corresponds to the Weinberg-Tomozawa (WT) term that can be written in momentum space as: The indices (i, j) cover all possible initial and final channels in the S = −1, Q = 0 sector. The matrix of coefficients C ij can be found in Table VII of Reference [28]. The normalization factor N is defined as N = (M + E)/(2M), with M and E being, respectively, the mass and energy of the baryon in the corresponding channel. The two-component Pauli spinor of a baryon with spin projection s is denoted by χ s . The symbol for the spin, only present in the spinors, should not be confused with √ s that represents the total energy of the meson-baryon system in the center-of-mass (CM) frame. The CM three-momentum of the meson (baryon) is given by q (− q). In Figure 1, the Born contributions are represented by diagrams (ii) (direct Born term) and (iii) (crossed Born term), the vertices of which are obtained from the D and F terms of Equation (4). The analytical form of the direct Born term is given by: and, similarly, the crossed Born term reads: where the label k refers to the intermediate baryon involved in the process. Here, u stands for the Mandelstam variable defined by u = (p i µ − q j µ ) 2 , with q j µ being the fourmomentum of the final meson and p i µ of the initial baryon. The coefficients C (Born) xy,z , which depend on the axial vector constants D and F, can be found in Apendix A of Reference [10].
Taking into account only the tree-level contributions at NLO (diagram (iv) in Figure 1), one obtains: where the D ij and L ij coefficients are constructed from combinations of the NLO parameters b 0 , b D , b F , d 1 , d 2 , d 3 and d 4 and are given in Table VIII of Reference [28], while the g ij and h ij matrices are listed in Appendix A.
One should bear in mind that the direct incorporation of the amplitudes given by Equations (7)-(10) into Equation (1) is simply not possible since such interaction kernels are composed of a mixture of contributions with different angular momenta. Thus, it is convenient to express the T-matrix in terms of the spin-nonflip and spin-flip parts.
where θ is the CM angle between the inital and final meson momenta andn = q j × q i /| q j × q i | is the normal vector to the scattering plane. The functions f ( √ s, θ) and g( √ s, θ) can be expanded in Legendre polynomials as: with f l and g l being the projections of the sum of all the above defined kernels onto P l (cos θ) and sin θ dP l (cos θ) d cos θ , respectively. The amplitudes can be redefined as: to ensure that the quantum numbers of spin ( 1 2 ), orbital angular momentum (l), and total angular momentum (J) are preserved in the unitarization procedure implemented by the BS equations. Following the notation of Reference [38], each unitarized J-scattering amplitude should be calculated by a new version of Equation (1), which in matrix form reads: where the amplitudes f tree l± are obtained from Equations (11)-(13) but employing, for the amplitude T ij (s, s ), the sum of the interaction kernels shown in Equations (7)- (10).
The aim of this work is the study of the effects caused by the inclusion of higher partial waves on the physical observables in this sector, particularly on the K − p cross-sections. The most general expression of the differential cross-sections for a given with where the first factor averages over the initial baryon spin projections, giving 1/2 for this particular case, and where we have also summed over all possible final baryon spin projections.
Since the scope of the present study does not go beyond d-wave contributions, focusing mostly on p-wave effects, the expression for the total cross-section incorporating such partial waves can be written as:

Models and Data Treatment
The availability of a large collection of experimental data makes the S = −1 sector an interesting benchmark to test meson-baryon EFTs, not only for checking their predictive power but also to extract information about the LECs, especially those beyond LO. Constructing a model with higher partial waves capable to describe the phenomenology of the S = −1 meson-baryon interaction provides, indeed, a tool to explore other sectors lacking experimental data, like S = −2, S = −3, and even to check the validity of SU(3) symmetry upon comparing the resulting parameters to those of existing chiral models that describe the meson-baryon interaction in the S = 0 sector.

Models
In this section, we introduce the models derived for this study and discuss the fitting procedure. Since the present study is the natural extension of the one presented in [30], we include a description of the old model as a reminder and to provide a better comprehension of the p-wave effects.
• s-wave (old): This first model corresponds to the fit called WT+Born+NLO carried out in [30]. It was constructed by adding the interaction kernels derived from the Lagrangian up to NLO and neglecting the hand g-terms. We limited this model to the s-wave contributions. • s-wave: The second model improves upon the first one by incorporating the novel hand g-terms that come from the NLO interaction kernel and, as in the first model, only the s-wave contribution is taken into account. • s+p-waves: This model employs the same Lagrangian as the s-wave model, but it also incorporates the p-wave contributions.

Data Treatment and Fitting Procedure
The parameters in our models originate from the chiral Lagrangian as well as from the employed regularization scheme. One finds 17 LECs in accordance with Equations (7)-(10), namely the meson decay constant f , the axial vector couplings D and F, and the NLO coefficients b 0 , b D , b F , d 1 , d 2 , d 3 , d 4 , g 1 , g 2 , g 4 , h 1 , h 2 , h 3 , and h 4 . It is a well-known fact that, in the UChPT models, the f parameter effectively takes larger values than the experimental one. Literature offers plenty of miscellaneous choices, but, in our particular case, in order to allow for a sort of average over all mesons involved in the various coupled channels, this parameter is constrained to vary from f = f π to f = 1.3 f π . The axial vector couplings D and F are allowed to vary within 15% of their canonical values in order to accommodate the spread of values found in the literature [29]. The NLO LECs are treated as completely free parameters in the fit. Additionally, the dimensional regularization introduces six subtraction constants, a πΣ , aK N , a πΛ , a ηΣ , a ηΛ , and a KΞ that are also parameters of the fit. However, these parameters are restricted to their natural-sized values, which should, as discussed in Reference [30], lie in the range [−10, 10] × 10 −3 for µ = 1 GeV. Table 1 displays the distribution per observable of the 224 experimental points employed in the current fitting procedure (the same as in [30]), most of them coming from the total cross-section for K − p scattering into different final channels [50][51][52][53][54][55][56][57][58][59][60][61][62][63][64][65]. Apart from this, the measured branching ratios of cross-section yields [66,67] are also fitted. Such yields can be defined from the elastic and inelastic K − p cross-sections evaluated at threshold: The precise measurement of the energy level shift and width of the atomic 1s state in kaonic hydrogen by the SIDDHARTA Collaboration [20] implied a substantial change in constraining the theoretical models. We incorporate this input in our fits via the K − p scattering length, which is obtained from the K − p scattering amplitude at threshold as: The scattering length is connected to the shift and width of kaonic hydrogen by means of the second-order corrected Deser-type formula [68]: where α is the fine-structure constant and µ r the reduced mass of the K − p system.
The minimizing criteria taken for the fitting procedure are based on the χ 2 per degree of freedom (χ 2 d.o.f. ). We avoid the use of the the standard definition of χ 2 d.o.f. since it favors observables with a larger number of experimental points over those with a smaller number. To get around this misleading effect, we adopt the method used in [9,23,24,26,[28][29][30], consisting of the assignment of equal weights to the different measurements by defining a renormalized χ 2 d.o.f. as: with Here, y exp i,k , y th i,k , and σ i,k are the experimental value, the theoretical prediction, and the experimental uncertainty of the ith point of the kth observable, respectively. Each observable has a total of n k points; K is the total number of observables, and p denotes the number of free fit parameters. In Equation (20), the renormalization is incorporated by averaging the χ 2 per degree of freedom over the different experiments.

Results and Discussion
As mentioned in the introduction, our previous works [28][29][30] studied the influence of the Born and NLO contributions, especially on the inelastic processes that open up at higher energies, namely the K − p → ηΛ, ηΣ 0 , K + Ξ − , K 0 Ξ 0 reactions. In the present work we include new ingredients, such as the hand g-terms of the NLO Lagrangian and the p-wave contributions, which are expected to play a relevant role at these energies. The effects of these new contributions on the model parameters are displayed in Table 2, where the parameters of the s-wave (old) model, obtained in Reference [30], are also included for comparison.
The parameters of the new s-wave model are displayed in the second column of Table 2. Although this fit has the lowest χ 2 d.o.f. of 0.77, we observe an overall discrepancy in the values of the parameters and in their accuracy compared to those obtained in our previous s-wave (old) study. It is important to stress that the error bars can not be directly compared because of the different methods employed to calculate them, as mentioned in the caption of Table 2. Despite this fact, some qualitative comparison can be made, as we discuss below.
Although we get natural-sized values for all the subtraction constants obtained with the two s-wave models, the parameter sets are very different. Moreover, the s-wave model provides larger values of the d i parameters compared to those of the s-wave (old) model, quite high values of the new g i 's and h i 's, as well as a small value of f / f π , almost at the lower edge, together with values of the axial vector constants at their upper edges. A plausible explanation is that the gand h-terms of the new model can play a role in accommodating the data at higher energies, where these terms are more important, and the fitting procedure forces them to be sizable at the expense of producing extreme values for the other parameters.
It is therefore expected that the incorporation of the p-waves in the fit, acting also more importantly at higher energies [36], might bring the parameters of the gand h-terms within a more natural range. This is indeed what happens for the s+p-waves model, whose parameters are compiled in the first column of Table 2. The parameters b i 's and d i 's become also more moderate in size, except for the value of d 4 , which shows big differences among the different models. We note that the volatility of this NLO LEC was already discussed in [30], and it was tied to the fact that this parameter only plays a role in the elastic transitions of the sector. Finally, we observe that the fit of the new s+p-waves model produces a value of f / f π very close to that of the old model, together with exactly the same values for D and F.  (20), for the different models described in the text. The subtraction constants are taken at a regularization scale µ = 1 GeV. The error bars in the parameters of the s-wave (old) are determined as explained in [30], while those of s+p-waves and s-wave are directly provided by the MINUIT [  The results of the threshold observables for the previous models are collected in Table 3, where they can be compared to the corresponding experimental values also included there. There is no substantial change in the reproduction of the experimental values when comparing the new models to the s-wave (old) one. In Figure 2, the calculated total cross-sections of K − p scattering to all channels of the S = −1 sector are compared to experimental data. As can be seen from the first three panel rows, the agreement of all models in reproducing the threshold observables is also reflected on the total cross-sections of the classical processes, namely K − p → K − p,K 0 n, π − Σ + , π + Σ − , π 0 Σ 0 , and π 0 Λ, which are studied at energies close to the K − p threshold. It is interesting to comment on the little enhancement of strength seen for the s+p-waves model at CM energies around 1500 MeV in the K − p → π 0 Λ cross-section. In this case, since the corresponding threshold is located several tens of MeV below the other channels, the momentum of the outgoing particles can reach higher values, thereby providing a more relevant role to the p-wave contribution. In contrast to what we see in the conventional channels, larger differences among the models are observed in the cross-sections of the K − p → ηΛ, ηΣ 0 , K + Ξ − , K 0 Ξ 0 processes, shown in the two bottom panel rows of Figure 2. In the case of the η channels, one can clearly see that all models describe qualitatively well the data, exhibiting minor discrepancies at energies below 1800 MeV. The differences among the models become appreciable at higher energies. From the study in [23,30], we know that the ηΛ and ηΣ 0 reactions are sensitive to the NLO (and Born) terms, to the extent that one cannot reproduce properly the experimental data in the case of ηΛ channel at this energy range, unless such terms are taken into account. Consequently, the disagreement observed between the s-wave (old) and s-wave models can be directly attributed to effects induced by the incorporation of the new gand h-terms. The additional p-wave contribution in the s+p-waves model, which gains significance as the momenta of particles involved increase, seems to provide more strength to the K − p → ηΛ cross-section around 1900 MeV and to the K − p → ηΣ 0 one from 2000 MeV on.
In the case of the KΞ channels, there are appreciable differences between the models in the whole energy range. In particular, the s+p-waves model for the K 0 Ξ 0 cross-section slightly overestimates the experimental data below 2000 MeV while remaining somewhat below the experimental structure that appears at 2100 MeV. On the contrary, this model offers the best description for the experimental points located at higher energies for both the K 0 Ξ 0 and K + Ξ − total cross-sections. This is in opposition to the s-wave and s-wave (old) models that, in such a region, accommodate the experimental points of the K 0 Ξ 0 and K + Ξ − cross-sections, respectively, at the edges of their errors bars. In general, the s+p-waves model is smoothly averaging over the whole set of points, and the s-wave models are able to provide a more marked structure around 2100 MeV.
For a deeper understanding of the role of the p-waves, in Figure 3 we present the individual partial-wave contributions to the K − p processes included in the s+p-waves model, namely the J P = 1 2 − (s-wave), 1 2 + and 3 2 + (p-wave) channels. Turning first to the classical processes (top six panels), we notice that p-wave effects are only noticeable in the π 0 Λ production process at the upper edge of the energy range, as discussed in the previous paragraph. The other cross-sections are essentially built from s-wave contributions at the energies explored here.
The situation becomes much more interesting for those reaction channels that open up at higher energies. Examining the K − p → ηΛ cross-section, one can observe that the s-wave J P = 1 2 − contribution (dotted red line) is dominant just above threshold due to the Λ(1670) resonance, generated dynamically in this partial wave. Then, as the energy increases, its role is moderately losing relevance against the J P = 3 2 + component (dash-dotted violet line), which is the main contribution in the region ranging from 1750 to 2100 MeV. As a matter of fact, the predominant character of the J P = 3 2 + component is another proof of the need to incorporate NLO terms in the kernel, as it can be seen that the WT and the Born s-channel terms do not contribute to this partial wave. Moreover, it is interesting to note that the present J P = 3 2 + contribution may replace the role of the explicit resonance Λ(1890) with these quantum numbers that was incorporated, in addition to the s-wave chiral amplitudes, in one of the models studied in [30] producing a similar description of the ηΛ cross-section above 1900 MeV. Concerning the J P = 1 2 + partial-wave (dashed green line), it is hardly contributing to the cross-section, which is in contrast to the coupled channel partial wave analysis carried out in Reference [22]. There, the corresponding crosssection contains a contribution coming from the J P = 1 2 + channel comparable in size to that of the J P = 3 2 + one, the latter having a similar strength as ours. As a last comment, we emphasize that the information encoded in this process can be very valuable, as it acts as an I = 0 filter, thereby helping to reduce the potential ambiguities in the amplitudes [30,31].  Next, we focus on another isospin filtering process, namely the K − p → ηΣ 0 reaction, which proceeds in the I = 1 channel. In this case, the s-wave J P = 1 2 − contribution is clearly the dominant one over the whole range of energies explored. The p-wave contributions start being noticeable from 1850 MeV, representing 20% of the total cross-section at most. However, one should be very careful about drawing conclusions because of the lack of experimental data in this channel beyond 1800 MeV.
The partial-wave decomposition of the KΞ total cross-section provides very valuable information about the relevance of higher partial waves in our model. On the one hand, the K 0 Ξ 0 production discloses the fundamental role played by the J P = 1 2 + contribution above 2200 MeV and how the low-energy regime is dominated by the J P = 1 2 − wave.
The J P = 3 2 + component provides a non-negligible strength around 2100 MeV, which is however not enough to reproduce the sizable experimental structure. On the other hand, in the K − p → K + Ξ − process, the J P = 1 2 − and J P = 1 2 + contributions are moderate and roughly constant, being the J P = 3 2 + component the one that governs the description of the experimental data. A subdominance of the s-wave component was also obtained for the π + p → K + Σ + process in Reference [36], where a similar chiral scheme up to NLO for the S = 0 sector was employed. The resonant model of Reference [22] produces K 0 Ξ 0 and K + Ξ − cross-sections that are dominated by the 3 2 + contribution with a sizable influence of The incorporation of the d-wave components of the interaction is the natural next step. It can be easily verified that these partial waves emerge from the crossed Born and the NLO kernels. It seems that a model containing the s-, p-and d-wave components simultaneously would provide the most realistic information about the NLO LECs of the Lagrangian (6). However, in practice, reaching a proper parametrization of the LECs when the d-waves are included is highly non-trivial, given the interdependence of the parameters contributing to the various partial waves at tree level, together with the complications of the coupled-channel scheme. Preliminary attempts have led us to the conclusion that a proper and controlled inclusion of higher partial waves in the fit requires the extension of the data set to include the available total and differential cross-sections at higher energies, an enormous task that is beyond the scope of the present work.
However, in order to illustrate the relevance of the d-waves, we take the parameters of the s+p-waves model and present, in Figure 4, the J P = 3 2 − and J P = 5 2 − contributions to the cross-sections of the K − p → ηΛ, ηΣ 0 , K + Ξ − , K 0 Ξ 0 reactions, which are the ones expected to be more substantially affected by the higher partial waves. As can be noticed from all panels in Figure 4, the d-wave contributions start to have a non-negligible effect about 100 MeV above the threshold of the corresponding reaction. In the cases of the ηΛ and K + Ξ − cross-sections, these partial waves contribute in a very remarkable way. This demonstrates that the d-wave terms seem to be an essential ingredient to describe the K − p inelastic cross-sections at higher energies, and, consequently, a fitting procedure including the effect of the d-wave terms needs to be performed. From a purely phenomenological point of view, it seems that the J P = 3   Table 2). See text for a detailed explanation of the content.
Finally, we analyze the pole content of the scattering amplitudes derived from the s+p-waves and s-wave models and compare with the corresponding pole content from the s-wave (old) model. The first thing to be mentioned is that we cannot dynamically generate any state, neither with J P = 1 2 + nor with J P = 3 2 + , from the s+p-waves model, in spite of the existence of some Λ * and Σ * resonances with these quantum numbers in the range covered by the higher-energy inelastic channels. According to the PDG compilation [70], these resonances decay mostly intoKN, πΣ, πΛ or to more exotic two-and three-body products. Thus, the inclusion of higher-energy experimental data from the K − p →KN, πΣ and πΛ processes in future fits could favor the appearance of poles in the scattering amplitude, since there are clear structures reflecting resonant signals in the corresponding cross sections.
The results of the pole position of the resonances with J P = 1 2 − , together with their couplings to the different channels, are compiled in Table 4 for the s+p-waves (top), the swave (middle) and the s-wave (old) (bottom) models. Table 4. Comparison of the pole positions between the models: s+p-waves, s-wave (old) and s-wave (in MeV) and the corresponding modulus of the couplings |g i | found in I = 0 and I = 1 channels.   In the I = 0 channel, all models generate the double-pole structure of Λ(1405) and a pole corresponding to Λ(1670). As can be appreciated from the panel on the top, the broader pole of Λ(1405) has been shifted 20 MeV towards lower energies compared to the corresponding resonance found in the s-wave (old) amplitudes. Despite the location of this first pole is not well established, given the scattered positions found in the literature, the new position of this pole is in better agreement with most of the studies devoted to this topic (see e.g., [27,30]). This is also the case for the corresponding pole in the s-wave model of the middle panel, but in this case it is shifted downwards by 55 MeV. The widths of these broad poles behave in a contrary way, while that of s+p-waves model gets narrower, compared to that of the s-wave (old) model, the one obtained from the s-wave model is even larger. All these broad poles couple most notably to πΣ andKN states, the coupling to πΣ being slightly higher. All models produce a higher energy and narrower Λ(1405) state. While the position in energy is similar for the three models, the pole found by the s-wave model is narrower by more than 10 MeV compared to those of found by the two other models. The pattern of the couplings for the narrow Λ(1405) in the three models is in accordance to most of the works in literature, meaning that this pole couples mostly to theKN channel but has a sizable coupling to the πΣ one. As for the Λ(1670) resonance, the most remarkable fact is that the s+p-waves model produces a width that is half of that obtained by the other two models and in much better agreement with the value quoted by the PDG [70].

s+p-waves
Concerning the I = 1 sector, the s-wave model does not generate any state, whereas the s+p-waves model produces a wide resonance about 100 MeV below in energy and substantially wider than that found by the s-wave (old) model. These states can hardly be identified with any of the Σ resonances listed in the PDG with J P = 1/2 − [70]. For instance, the observed 1-star Σ(1620) resonance is substantially narrower and has been seen to decay into πΣ states with twice as much probability than into πΛ states, a decay pattern that cannot be reproduced by neither of the I = 1 states listed in Table 4 according to the size of their couplings to the different meson-baryon channels. The identification with the Σ(1750) resonance is also discarded, as it has been seen to decay more strongly into ηΣ states and, with a somewhat smaller branching ratio, into πΛ and πΣ. We note, however, that the large width of the I = 1 states found here hinders their experimental identification in reactions where they might be produced.

Conclusions
We have performed a new study of the meson-baryon interaction in the S = −1 sector including both s-and p-waves, aimed at improving our knowledge about the NLO terms of the chiral SU(3) Lagrangian. We pay special attention to the processes that are very sensitive to these subleading terms of the Lagrangian, such as the K − p → K + Ξ − , K 0 Ξ 0 , ηΛ, ηΣ 0 reactions. This work presents a detailed derivation of the formalism needed to obtain the higher partial waves of the scattering amplitudes, beyond the s-wave component usually considered in the literature.
Our previous models evolved in the direction of implementing systematic improvement to obtain a more reliable determination of the NLO terms of the Lagrangian, namely the b i and d i LECs, culminating into the s-wave (old) set obtained in [30]. In the present work we have first improved upon the later model by incorporating the novel gand h-terms from the NLO interaction kernel. This has led to the new s-wave set of parameters, which reproduces experimental data very well and achieves very low χ 2 d.o. f . of 0.77. However, we have observed that the LECs of this s-wave model, both the LO and NLO ones, are qualitatively very different from our earlier estimations.
We have then checked how stable is this new s-wave set by incorporating into the model the p-wave contributions, which strongly depend on the novel gand h-terms but also on the rest of NLO parameters. With respect to the s-wave results, the s+p-waves model also produces a very good fit, although with a slightly worse χ 2 d.o. f . value of 0.86 and a rather different set of LECs. The latter fact is an indication of the clear interconnection between the higher partial waves and the higher orders of the chiral expansion, as both contributions acquire importance at higher energies. It is interesting to note, however, that the size of the LECs obtained with the s+p-waves model is smaller and closer to that of our earlier s-wave (old) model. This makes us believe that the s+p-waves model developed in the present work constitutes a very good starting point for subsequent implementations of higher partial waves. In this work, we have shown that a notable fraction of the crosssections of the K − p → K + Ξ − , K 0 Ξ 0 , ηΛ, ηΣ 0 reactions comes indeed from the novel p-wave contributions, thereby making it clear that higher partial waves cannot be considered as minor corrections in these processes.
The main conclusion of this study is that one does need to incorporate higher partialwave contributions to the scattering amplitudes in order to properly describe the K − p inelastic reactions opening up at higher energies. This is in contrast to our previous works, where the experimental data of such processes were reasonably reproduced, employing pure s-wave scattering amplitudes. In other words, our previous models have effectively overestimated the lowest partial-wave contributions, thereby masking the physics of higher partial waves behind the values of the NLO parameters. As for the relevance of the new gand h-terms considered in this work, compared to the rest of NLO contributions, we cannot say anything conclusive, as we have obtained very different values of the parameters in the two new fits presented here. What seems to be clear is that they provide a non-negligible contribution to the scattering amplitudes, especially far enough from thresholds, because of their strong dependence on momenta. Therefore, it is natural to consider the gand hterms when implementing higher partial waves in the scattering amplitudes, since both contributions acquire a relevant role at higher energies.
As a final remark we would like to comment that the s+p-waves model is the natural continuation of our previous ones [28][29][30]. Since it employs the same set of data points in the fit, it provides a clear quantification of the effect of the new ingredients, namely the gand h-terms and the p-waves. This model should then be considered as a promising reference point upon which additional partial-wave contributions, in particular the d-waves, need to be added. However, this challenging study, which is already in progress, requires including more experimental data points of total and differential cross-sections within a wider range of energies in the fits. As can be clearly seen in Table 2, the determination of LECs becomes an extensive nonlinear optimization problem where the relatively large number of parameters makes it challenging to determine their optimal values. Thus, it is also crucial to establish a credible program to quantify the magnitude of statistical and systematic uncertainties of the LECs and to study their propagation into experimental observables. This can be achieved by employing advanced mathematical optimization tools, as introduced, e.g., in [71,72]. We believe that this is the proper procedure to have a realistic determination of the NLO parameters of the chiral SU(3) Lagrangian. Funding: V.M. and A.R. acknowledge support from the Ministerio de Ciencia, Innovación y Universidades under the project CEX2019-000918-M of ICCUB (Unidad de Excelencia "María de Maeztu"), and from the Ministerio de Economía y Competitividad, with additional European Regional Development Funds, under the contract FIS2017-87534-P. A.F. acknowledges support from the from the Spanish Ministerio de Ciencia, Innovación y Universidades under the project PID2020-112777GB-I00, from the Ministerio de Economía y Competitividad and the European Regional Development Fund, under contract FIS2017-84038-C2-1-P, and from the Generalitat Valenciana under contract PROME-TEO/2020/023. V.M., A.R. and A.F. also acknowledge support from the EU STRONG-2020 project under the program H2020-INFRAIA-2018-1, grant agreement no. 824093. The work of D.G. was supported by the Czech Science Foundation, GAČR grant no. 19-19640S.

Conflicts of Interest:
The authors declare no conflict of interest.