Triangular singularity and a possible resonance in the decay
Abstract
We study the decay by considering a triangle singularity mechanism. In this mechanism, the decays into the , the decays into the (or ), and then the (or ) interact to produce the in the final state. This mechanism produces a peak structure around MeV. In addition, the possibility that there is a hiddenstrange pentaquarklike state is also considered by taking into account the final state interactions of , , and . We conclude that it is difficult to search for the hiddenstrange analogue of the states in this decay. However, we do expect nontrivial behavior in the invariant mass distribution. The predictions can be tested by experiments such as BESIII, LHCb and BelleII.
I Introduction
In 2015, two hiddencharm pentaquarklike structures, and , were observed in the invariant mass spectrum via the decay by the LHCb Collaboration Aaij:2015tga (). After they were observed, the two resonances were investigated within multiple theoretical schemes with the aim to explain their nature (for more details and references, see the recent reviews Chen:2016qju (); Guo:2017jvc ()). The existence of pentaquarks with hidden charm in that mass region was already predicted in Refs. Wu:2010jy (); Wu:2010vk (); Wang:2011rga (); Yang:2011wz (); Xiao:2013yca () by studying the interactions of anticharm mesons and charm baryons using different models. Furthermore, it was pointed out in Ref. Guo:2015umn () (see also Refs. Liu:2015fea (); Guo:2016bkl ()) that a triangle singularity is located very close to the threshold, GeV, and thus at the mass. Such a singularity could produce a narrow peak mimicking the behavior of a narrow resonance, which requires the and proton to rescatter in an wave into the final state Bayar:2016ftu (). This would require quantum numbers or for the peak. Notice that although the and were reported as the most preferred quantum numbers in the original LHCb publication Aaij:2015tga (), remains one of the favored possibilities in a later experimental analysis using an extended model Jurik:2016bdm (). Clearly, further investigations on the and structures, in particular from more processes and more experiments, are needed. Since the structures were observed in the decay mode , it is natural to expect that these states, were they hadronic resonances, can be produced in photoproduction process where they will appear as channel resonances Wang:2015jsa (); Kubarovsky:2016whd (); Meziani:2016lhg (); Blin:2016dlf ().
Analogous to the hiddencharm pentaquark states, one may consider the possible existence of hiddenstrange pentaquarks , in which the pair is replaced by . In fact, in the light flavor sector below GeV, understanding the nature of the resonance with spin parity is very challenging Klempt:2007cp (); Crede:2013sze (). One peculiar property of the is that it couples strongly to the channels with strangeness, such as the and , which is difficult to understand in the classical threeconstituentquark models. This finds an explanation within the chiral unitary approach in the work of Ref. Inoue:2001ip (). The strange decay properties of the resonance can also be easily understood by considering large fivequark components in it Liu:2005pm (); Helminen:2000jb (); Zou:2007mk (); An:2008xk (). Within this pentaquark picture, the resonance could be the lowest orbitally excited state with a large admixture of pentaquark component. This makes the heavier than the and also gives a natural explanation of its large couplings to the channels with strangeness Zou:2010tc (). In a very recent quark model study Gao:2017hya (), a state with a mass varying from to MeV is obtained, and its main component is . This state could correspond to the resonance which has only an overall twostar status according to the review by the Particle Data Group (PDG) Olive:2016xmw (). Its existence is supported by the analysis of the new photoproduction data Kashevarov:2017kqb (); Collins:2017sgu (), which finds that the with is crucial in order to describe the cusp observed in the photoproduction at around MeV as well as the fast nearthreshold rise of the total cross section of the reaction Kashevarov:2017kqb (). In Refs. Kashevarov:2017kqb (); Collins:2017sgu (), it was also pointed out that the has strong couplings to both the and channels.
At around GeV, a bound state is predicted in several models Gao:2000az (); Huang:2005gw (); Gao:2017hya (). Such a state can be viewed as a pentaquark. In Ref. Gao:2017hya (), a state dominated by the component is obtained with a mass varying from to MeV. Independently, a resonance ^{1}^{1}1In the editions of the PDG review before 2012, all the evidence for a state with a mass above MeV was filed under a twostar . There is now evidence Anisovich:2011fc () of two states in this region, and the PDG has associate the older data (according to masses) to two states: a threestar and a twostar Olive:2016xmw (). with a mass about GeV is proposed to explain the experimental results Kohri:2009xe (); Kiswandhi:2010ub (); Xie:2010yk (); Kiswandhi:2011cq (); Kim:2011rm (); He:2012ud (); Xie:2013mua (); Seraydaryan:2013ija (); Dey:2014tfa (); Dey:2014npa (); Kiswandhi:2016cav () on the associated strangeness production reactions , , and . The forwarddirection enhancement at around GeV in the reaction can be also reproduced by including a special correlated fivequark configuration of a colorantitriplet () diquark and a colortriplet , which subsequently hadronize into the and proton Lebed:2015fpa (); Lebed:2015dca (). However, it is pointed that such a fivequark configuration is not literally a resonant pentaquark state Lebed:2015fpa (); Lebed:2015dca (). In Ref. He:2017aps (), it is proposed that the states and in the photoproduction are hadronic molecular states from the and interactions, respectively, and they can be regarded as the hiddenstrange partners of the LHCb pentaquarks.
The structures were produced in the process . Analogously, one may study the possible states in the singly Cabibbo suppressed process . As pointed out in Ref. Lebed:2015dca (), the and are entirely comparable if one substitutes . One important difference between the two processes is that the former has a much smaller phase space—the is above the threebody threshold only by 193 MeV. Such a small phase space for the decay restricts that only the neutral pion is possible in the final state, for a hadronic decay, if we want to produce in addition a pair. Because of the small phase space and the weak interaction^{2}^{2}2The interaction should be very weak for two reasons. Firstly, the pion and the meson do not have the same quark flavors, which leads to an Okubo–Zweig–Iizuka (OZI) suppression. Secondly, the small phase space means that the pion is soft, and the interaction between a soft pion and matter fields is weak because of the spontaneous breaking of chiral symmetry in quantum chromodynamics., no other resonances except for the possible contribute to the process. The first experimental measurement of the process has been reported by the Belle Collaboration Pal:2017ypp (). Very recently, the Belle Collaboration reported their searching for the decay of , and no significant signal was observed with an upper limit on the branching fraction of at a confidence level Pal:2017ypp ().
In this paper, we will show that the also receives a contribution from triangle singularities close to the physical region. A triangle singularity appears on the physical boundary in a particular situation when all the intermediate states are on shell, and all the particles move along the same direction (parallel or antiparallel) such that the interactions at all three vertices can happen as classical processes Coleman:1965xm (). Such a physical picture can be easily seen following the analysis of Ref. Bayar:2016ftu (). In addition to the works related to the structures mentioned above, the role played by triangle singularities has been broadly investigated recently in the literature Wu:2011yx (); Aceti:2012dj (); Wu:2012pg (); Liu:2014spa (); Szczepaniak:2015eza (); Ketzer:2015tqa (); Liu:2015taa (); Szczepaniak:2015hya (); Aceti:2016yeb (); Yang:2016sws (); Wang:2016dtb (); Xie:2016lvs (); Debastiani:2016xgg (); Roca:2017bvy (); Debastiani:2017dlz (); Samart:2017scf (); Sakai:2017hpg (). Along this line, we will calculate the triangle singularity contribution to the decay, where the decays into , the () decays to the (or ) and the (or ) rescatter into in the final state, see Fig. 1. In addition to the effects of the triangle mechanism, we consider also the final state interaction (FSI) of and . Were there a resonance, it must couple to both the and and thus may be manifest in the Dalitz plot or in the invariant mass distribution. Yet, because of the small phase space and depending on the mass and width of such a state, it could be difficult to search for it. As will be shown in this paper, on one hand the triangle singularity contribution can enhance the production of such a resonance, on the other hand it makes the identification of the signal more difficult if its mass is around 2.02 GeV.
Ii Formalism
The decay can proceed through the triangle diagrams depicted in Fig. 1^{3}^{3}3Replacing the by the leads to vanishing contribution because the cannot couple to the .. Given the masses of the initial state , the neutral pion in the final state and two of the intermediate states, for example the and , the region for the mass in order to produce a triangle singularity at the physical boundary, i.e., in the physical region^{4}^{4}4The triangle singularity of course cannot be exactly in the physical region since otherwise one would get a logarithmically divergent amplitude. It is shifted into the complex plane because of the finite decay width of at least one of the intermediate states. The amplitude in the physical region is well defined without any divergence. can be worked out Guo:2015umn (); Liu:2015taa (). Using the central values for all of the mentioned hadron masses, the region can be obtained as MeV for diagram (A), while the measured mass of the , MeV, is 4 MeV below. The region is MeV for the exchange in diagram (B), while the measured mass of the MeV almost reaches the lower bound. In this case, the triangle singularities still have sizeable influence on the physical decay amplitude. The lower bound of that region means that the triangle singularity in the invariant mass is located exactly at the twobody threshold of the two particles which rescatter into the . Thus, one expects that the triangle singularity induced effects would be mainly because of the exchange diagram which can produce an enhancement around the threshold at 2007 MeV. In the following, we give the details of the calculation which shows explicitly the enhancement around that energy.
The decay of can proceed by exchange diagram Chau:1995gk (); Lu:2009cm (); Cheng:2010vk (): , and the are hadronized, together with a pair with the vacuum quantum numbers, into the .
The evaluation of the diagram in in Fig. 1 requires us first to provide an expression for the vertex. Because the threshold ( MeV) is very close to the mass of , we consider only the wave coupling. Then we can write
(1) 
where is the isospin factor with for the and for the , and is an effective coupling constant which can be obtained, in general, from the branching ratio of .
The decays of and are in waves, then we can easily write with symmetry
(2)  
(3) 
with obtained from the total decay width MeV and the branching fraction .
After the production of the and , they rescatter into the in the final state, as shown in Fig. 1. The total decay amplitude for the processes shown in Fig. 1 can be written as
(4)  
where we have defined and , and and are matrix elements for the rescattering processes, which will be discussed in the next section. We notice that the mass threshold is close to the mass of and the range of the invariant mass for the decay of interest, MeV, allows us to make nonrelativistic apporoximation for all the involved baryons and vector mesons. Therefore, we can consider only waves for the rescattering. Furthermore, we can make the approximation
(5) 
where the sum runs over the polarizations of the meson.
After performing the contour integration over the temporal component in Eq. (6), in the same way as shown in Refs. Bayar:2016ftu (); Aceti:2015zva (), and including the finite widths of the and resonances, we get
(6)  
where , , , , and with or . Because the wave vertices attached to the initial state and the final state do not introduce any momentum dependence into the loop amplitude, and the wave pionic vertices result in a factor of the pion momentum, the above loop integral is ultraviolet convergent.
The invariant mass mass distribution for the decay then reads
(7) 
where is the momentum in the rest frame of the , and is the momentum in the centerofmass frame of the system. They are given by
with MeV, MeV, and MeV. Finally, the partial decay width of the decay is obtained by integrating Eq. (7) over ,
(8) 
Iii Numerical results
So far we have not specified the input for the rescattering matrix elements. In principle, because of the very small phase space and the closeness of the thresholds, all of the involved hadrons, , or , and proton, can be treated nonrelativistically. Thus, one may construct a nonrelativistic effective field theory describing the interaction between vector mesons and baryons with the leading order defined by a few constant contact terms. However, it does not make much sense doing it in that manner because of the lack of experimental information. We will thus take the model of Ref. Oset:2009vf () where the interaction of the vector mesons with the SU(3) octet baryons is studied in the local hidden gauge formalism using a coupledchannel unitary approach. In that model, a degenerate pair of resonances with and which couple strongly to , and is obtained, and the pole is at MeV Oset:2009vf (). They can be regarded as the states. The prediction in this model was updated in light of the data Ewald:2011gw () in Ref. Ramos:2013wua () to get resonance parameters with a mass of MeV and a width of MeV. However, this state only shows up in the transitions involving the channel. One may regard that model as providing a special set of parameters for the nonrelativistic effective field theory mentioned above. By adjusting the interaction strengths, one can in principle investigate the possibility of with other masses and as well as the possibility without any , i.e. no pole around the threshold.
Here we present the numerical results for the invariant mass distribution for three different cases, which are denoted as Model I, II and III, in Fig. 2. Model I represents the calculation of the triangle diagrams in Fig. 1 with the FSI taken from Ref. Ramos:2013wua (), which includes the contribution of a state with properties specified in that model, which could be well different in other models and in reality, see above. Model II is different from Model I by modelling the FSI by a constant, and it thus represents the case without any resonance. For comparison, we show the phase space without any special dynamics as Model III. To be more explicit, for these three cases, the total decay amplitudes (, II, and III) are given by
(9)  
(10) 
where is the amplitude shown in Eq. (6), and are normalization constants to be adjusted to match the measured event distribution, MeV is the pion decay constant, and are the energies of the and mesons in the centerofmass frame. Here we take MeV and MeV, which are obtained at the mass threshold.
In Fig. 2, the solid, dashed, and dotted curves represent the results of Model I, II, and III, respectively. The parameter of Model II has been adjusted to the strength of the experimental data reported by the Belle Collaboration Pal:2017ypp () at its peak around MeV. The results of Model I and III are normalized such as to have the same integrated partial width as Model II. Model II clearly shows a peak structure around 2.02 GeV. The origin of this peak is the triangle diagrams, in particular the exchange in Fig. 1 (B) which has a triangle singularity close to the threshold ( GeV). The width of this peak is comparable with the width of the , which is about 50 MeV. This is a quite natural consequence as the invariant mass is the same as the invariant mass so that its distribution inherits the width of the . Were the width much smaller, one would get a much narrower peak. For Model I, one might think that there should be also a bump structure around MeV which is the mass of the generated resonance in the vector meson–baryon interaction model we are using Ramos:2013wua (). However, the triangle diagram involving the transition is the predominant contribution in the present case because its triangle singularity is closer to the physical region, while the resonance peak only shows up in the channels involving the .
Here, the FSI results in a nearthreshold enhancement, see Fig. 3 where the kink in the solid line is located at the threshold. It could be that in other models the resonance couples to these vector meson–baryon channels in a different pattern so as to show up as a nearthreshold peak in the invariant mass distribution. Identifying such an enhancement in experiments is difficult as it requires the data to have a high statistics. In particular, it becomes much more difficult if the mass is close to the threshold because of the presence of kinematic singularities there. In any case, the phase space shown as Model III is very different from both Model I and Model II. Despite the low statistics of the current Belle data, the curve of Model II, whose shape is completely fixed, has a remarkable agreement with the data. In particular, the data seem to indeed have a peak around threshold. More data are welcome to clarify the situation.
Using the value of estimated in Appendix A, we can get an estimate of the branching fraction of the threebody decay by integrating over the invariant mass distribution. For Model II with , the result is
(11) 
which is of the same order as the Belle upper limit Pal:2017ypp (). Thus, although not all of the contributions to the decay are from this mechanism we expect the actual branching fraction is of this order.
Iv Summary
We study the decay by considering a triangle singularity mechanism. The decay was proposed to be a channel to search for the hiddenstrange partner of the states. The mechanism is such that the decays into the , the subsequently decays into the (or ), and the then interacts with the (or ) to produce the in the final state. In the FSI, we consider cases with and without a state. For the case with the , we take the model of Refs. Oset:2009vf (); Ramos:2013wua () which produces a resonance at around 2 GeV. The triangle singularities considered in this paper are close to the physical region, and can produce a peak at around 2.02 GeV with a width similar to that of the resonance. The obtained invariant mass distribution agrees with the existing Belle data. Were there a state, it could distort the distribution. However, it is difficult to be identified in the decay under study because of the small phase space and the presence of triangle singularities. We look forward to more data from the BESIII, BelleII and LHCb experiments in the future, which will be decisive to illuminate the role played by triangle singularities in this decay.
Acknowledgments
We would like to thank Eulogio Oset and Wei Wang for useful discussions. This work is partly supported by the National Natural Science Foundation of China (NSFC) under Grant Nos. 11475227 and 11647601, by the DFG and NSFC through founds provided to the SinoGerman CRC 110 “Symmetries and the Emergence of Structure in QCD” (NSFC Grant No. 11621131001), by the CAS Key Research Program of Frontier Sciences (Grant No. QYZDBSSWSYS013), by the Youth Innovation Promotion Association of CAS (Grant No. 2016367), and by the Thousand Talents Plan for Young Professionals.
Appendix A Estimate of the coupling constant
The branching fraction for has not been estimated so far. Yet, an upper limit has been reported as Olive:2016xmw (). Because the and decay dominantly into the and , we thus take as an orderofmagnitude estimate for to estimate the coupling constant using the following decay width formula
(12) 
with
Using the measured masses MeV, MeV, MeV and the total decay width of MeV, we get
(13) 
99
 (1)

R. Aaij et al. [LHCb Collaboration], Phys. Rev. Lett. 115, 072001 (2015) [arXiv:1507.03414 [hepex]].
 (2)
 (3)
 (4)
 (5)
 (6)
 (7)
 (8)
 (9)
 (10)
 (11)
 (12)
 (13)
 (14)
 (15)
 (16)
 (17)
 (18)
 (19)
 (20)
 (21)
 (22)
 (23)
 (24)
 (25)
 (26)
 (27)
 (28)
 (29)
 (30)
 (31)
 (32)
 (33)
 (34)
 (35)
 (36)
 (37)
 (38)
 (39)
 (40)
 (41)
 (42)
 (43)
 (44)
 (45)
 (46)
 (47)
 (48)
 (49)
 (50)
 (51)
 (52)
 (53)
 (54)
 (55)
 (56)
 (57)
 (58)
 (59)
 (60)
 (61)
 (62)
 (63)
 (64)
 (65)
 (66)
 (67)
 (68)
 (69)
 (70)
 (71)
 (72)
H. X. Chen, W. Chen, X. Liu and S. L. Zhu, Phys. Rept. 639, 1 (2016) [arXiv:1601.02092 [hepph]].
F.K. Guo, C. Hanhart, U.G. Meißner, Q. Wang, Q. Zhao and B.S. Zou, arXiv:1705.00141 [hepph].
J.J. Wu, R. Molina, E. Oset and B.S. Zou, Phys. Rev. Lett. 105, 232001 (2010) [arXiv:1007.0573 [nuclth]].
J.J. Wu, R. Molina, E. Oset and B.S. Zou, Phys. Rev. C 84, 015202 (2011) [arXiv:1011.2399 [nuclth]].
W. L. Wang, F. Huang, Z. Y. Zhang and B. S. Zou, Phys. Rev. C 84, 015203 (2011) [arXiv:1101.0453 [nuclth]].
Z.C. Yang, Z.F. Sun, J. He, X. Liu and S.L. Zhu, Chin. Phys. C 36, 6 (2012) [arXiv:1105.2901 [hepph]].
C.W. Xiao, J. Nieves and E. Oset, Phys. Rev. D 88, 056012 (2013) [arXiv:1304.5368 [hepph]].
F.K. Guo, U.G. Meißner, W. Wang and Z. Yang, Phys. Rev. D 92, 071502 (2015) [arXiv:1507.04950 [hepph]].
X.H. Liu, Q. Wang and Q. Zhao, Phys. Lett. B 757, 231 (2016) [arXiv:1507.05359 [hepph]].
F.K. Guo, U.G. Meißner, J. Nieves and Z. Yang, Eur. Phys. J. A 52, 318 (2016) [arXiv:1605.05113 [hepph]].
M. Bayar, F. Aceti, F.K. Guo and E. Oset, Phys. Rev. D 94, 074039 (2016) [arXiv:1609.04133 [hepph]].
N. P. Jurik, Observation of resonances consistent with pentaquark states in decays, CERNTHESIS2016086.
Q. Wang, X.H. Liu and Q. Zhao, Phys. Rev. D 92, 034022 (2015) [arXiv:1508.00339 [hepph]].
V. Kubarovsky and M. B. Voloshin, arXiv:1609.00050 [hepph].
Z. E. Meziani et al., arXiv:1609.00676 [hepex].
A. N. Hiller Blin, C. FernándezRamírez, A. Jackura, V. Mathieu, V. I. Mokeev, A. Pilloni and A. P. Szczepaniak, Phys. Rev. D 94, 034002 (2016) [arXiv:1606.08912 [hepph]].
E. Klempt and A. Zaitsev, Phys. Rept. 454, 1 (2007) [arXiv:0708.4016 [hepph]].
V. Crede and W. Roberts, Rept. Prog. Phys. 76, 076301 (2013) [arXiv:1302.7299 [nuclex]].
T. Inoue, E. Oset and M. J. Vicente Vacas, Phys. Rev. C 65, 035204 (2002) [hepph/0110333].
B.C. Liu and B.S. Zou, Phys. Rev. Lett. 96, 042002 (2006) [nuclth/0503069].
C. Helminen and D. O. Riska, Nucl. Phys. A 699, 624 (2002) [nuclth/0011071].
B.S. Zou, Eur. Phys. J. A 35, 325 (2008) [arXiv:0711.4860 [nuclth]].
C.S. An and B.S. Zou, Eur. Phys. J. A 39, 195 (2009) [arXiv:0802.3996 [nuclth]].
B.S. Zou, Nucl. Phys. A 835, 199 (2010) [arXiv:1001.1084 [nuclth]].
H. Gao, H. Huang, T. Liu, J. Ping, F. Wang and Z. Zhao, Phys. Rev. C 95, 055202 (2017) [arXiv:1701.03210 [hepph]].
C. Patrignani et al. [Particle Data Group], Chin. Phys. C 40, 100001 (2016).
V. L. Kashevarov et al., Phys. Rev. Lett. 118, 212001 (2017) [arXiv:1701.04809 [nuclex]].
P. Collins et al., Phys. Lett. B 771, 213 (2017) [arXiv:1703.00433 [nuclex]].
H. Gao, T. S. H. Lee and V. Marinov, Phys. Rev. C 63, 022201 (2001) [nuclth/0010042].
F. Huang, Z.Y. Zhang and Y.W. Yu, Phys. Rev. C 73, 025207 (2006) [nuclth/0512079].
A. V. Anisovich, R. Beck, E. Klempt, V. A. Nikonov, A. V. Sarantsev and U. Thoma, Eur. Phys. J. A 48, 15 (2012) [arXiv:1112.4937 [hepph]].
H. Kohri et al. [LEPS Collaboration], Phys. Rev. Lett. 104, 172001 (2010) [arXiv:0906.0197 [hepex]].
A. Kiswandhi, J.J. Xie and S. N. Yang, Phys. Lett. B 691, 214 (2010) [arXiv:1005.2105 [hepph]].
J.J. Xie and J. Nieves, Phys. Rev. C 82, 045205 (2010) [arXiv:1007.3141 [nuclth]].
A. Kiswandhi and S. N. Yang, Phys. Rev. C 86, 015203 (2012) Erratum: [Phys. Rev. C 86, 019904 (2012)] [arXiv:1112.6105 [nuclth]].
S. H. Kim, S. i. Nam, Y. Oh and H. C. Kim, Phys. Rev. D 84, 114023 (2011) [arXiv:1110.6515 [hepph]].
J. He and X.R. Chen, Phys. Rev. C 86, 035204 (2012) [arXiv:1209.4138 [nuclth]].
J.J. Xie, E. Wang and J. Nieves, Phys. Rev. C 89, 015203 (2014) [arXiv:1309.7135 [nuclth]].
H. Seraydaryan et al. [CLAS Collaboration], Phys. Rev. C 89, 055206 (2014) [arXiv:1308.1363 [hepex]].
B. Dey et al. [CLAS Collaboration], Phys. Rev. C 89, 055208 (2014) Addendum: [Phys. Rev. C 90, 019901 (2014)] [arXiv:1403.2110 [nuclex]].
B. Dey, arXiv:1403.3730 [hepex].
A. Kiswandhi, S. N. Yang and Y.B. Dong, Phys. Rev. C 94, 015202 (2016) [arXiv:1604.01555 [nuclth]].
R. F. Lebed, Phys. Rev. D 92, 114006 (2015) [arXiv:1510.01412 [hepph]].
R. F. Lebed, Phys. Rev. D 92, 114030 (2015) [arXiv:1510.06648 [hepph]].
J. He, Phys. Rev. D 95, 074031 (2017) [arXiv:1701.03738 [hepph]].
B. Pal et al. [Belle Collaboration], arXiv:1707.00089 [hepex].
S. Coleman and R. E. Norton, Nuovo Cim. 38, 438 (1965).
J.J. Wu, X.H. Liu, Q. Zhao and B.S. Zou, Phys. Rev. Lett. 108, 081803 (2012) [arXiv:1108.3772 [hepph]].
F. Aceti, W.H. Liang, E. Oset, J.J. Wu and B.S. Zou, Phys. Rev. D 86, 114007 (2012) [arXiv:1209.6507 [hepph]].
X.G. Wu, J.J. Wu, Q. Zhao and B.S. Zou, Phys. Rev. D 87, 014023 (2013) [arXiv:1211.2148 [hepph]].
X. H. Liu, Phys. Rev. D 90, 074004 (2014) [arXiv:1403.2818 [hepph]].
A. P. Szczepaniak, Phys. Lett. B 747, 410 (2015) [arXiv:1501.01691 [hepph]].
M. Mikhasenko, B. Ketzer and A. Sarantsev, Phys. Rev. D 91, 094015 (2015) [arXiv:1501.07023 [hepph]].
X.H. Liu, M. Oka and Q. Zhao, Phys. Lett. B 753, 297 (2016) [arXiv:1507.01674 [hepph]].
A. P. Szczepaniak, Phys. Lett. B 757, 61 (2016) [arXiv:1510.01789 [hepph]].
F. Aceti, L.R. Dai and E. Oset, Phys. Rev. D 94, 096015 (2016) [arXiv:1606.06893 [hepph]].
Z. Yang, Q. Wang and U.G. Meißner, Phys. Lett. B 767, 470 (2017) [arXiv:1609.08807 [hepph]].
E. Wang, J. J. Xie, W. H. Liang, F.K. Guo and E. Oset, Phys. Rev. C 95, 015205 (2017) [arXiv:1610.07117 [hepph]].
J.J. Xie, L.S. Geng and E. Oset, Phys. Rev. D 95, 034004 (2017) [arXiv:1610.09592 [hepph]].
V. R. Debastiani, F. Aceti, W.H. Liang and E. Oset, Phys. Rev. D 95, 034015 (2017) [arXiv:1611.05383 [hepph]].
L. Roca and E. Oset, Phys. Rev. C 95, 065211 (2017) [arXiv:1702.07220 [hepph]].
V. R. Debastiani, S. Sakai and E. Oset, Phys. Rev. C 96, 025201 (2017) [arXiv:1703.01254 [hepph]].
D. Samart, W.H. Liang and E. Oset, arXiv:1703.09872 [hepph].
S. Sakai, E. Oset and A. Ramos, arXiv:1705.03694 [hepph].
L.L. Chau, H.Y. Cheng and B. Tseng, Phys. Rev. D 54, 2132 (1996) [hepph/9508382].
C.D. Lu, Y.M. Wang, H. Zou, A. Ali and G. Kramer, Phys. Rev. D 80, 034011 (2009) [arXiv:0906.1479 [hepph]].
H.Y. Cheng and C.W. Chiang, Phys. Rev. D 81, 074031 (2010) [arXiv:1002.2466 [hepph]].
F. Aceti, J. M. Dias and E. Oset, Eur. Phys. J. A 51, 48 (2015) [arXiv:1501.06505 [hepph]].
E. Oset and A. Ramos, Eur. Phys. J. A 44, 445 (2010) [arXiv:0905.0973 [hepph]].
R. Ewald et al. [CBELSA/TAPS Collaboration], Phys. Lett. B 713, 180 (2012) [arXiv:1112.0811 [nuclex]].
A. Ramos and E. Oset, Phys. Lett. B 727, 287 (2013) [arXiv:1304.7975 [nuclth]].