Structure, formation and decay of $\bar{K}NN$ system by Faddeev-AGS calculations
Sajjad Marri, S. Z. Kalantari, J. Esmaili

TL;DR
This paper uses Faddeev-AGS calculations to analyze the $ar{K}NN$ system, demonstrating the potential to observe the $K^{-}pp$ quasi-bound state and branch points in $ar{K}N- ext{pi}\Sigma$ interactions through mass spectra.
Contribution
It introduces a detailed Faddeev-AGS approach with various potentials to identify signatures of the $K^{-}pp$ state and branch points in experimental observables.
Findings
Signature of $K^{-}pp$ quasi-bound state can be observed in mass spectra.
Trace of branch points can be detected in the observables.
Different interaction models consistently show these features.
Abstract
The Faddeev AGS equations are solved for coupled-channels system with quantum numbers and . Using separable potentials for interaction, we have calculated the transition probability for the reaction. The possibility to observe the trace of quasi-bound state in the mass spectra was studied. Various types of chiral based and phenomenological potentials are used to describe the interaction. It was shown that not only we can see the signature of the quasi-bound state in the mass spectra, but also, one can see the trace of branch points in the observables.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Structure, formation and decay of system by Faddeev-AGS calculations
S. Marri
Department of Physics, Isfahan University of Technology, Isfahan 84156-83111, Iran
S. Z. Kalantari
Department of Physics, Isfahan University of Technology, Isfahan 84156-83111, Iran
J. Esmaili
Department of Physics, Faculty of Basic Sciences, Shahrekord University, Shahrekord, 115, Iran
Abstract
The Faddeev AGS equations are solved for coupled-channels system with quantum numbers and . Using separable potentials for interaction, we have calculated the transition probability for the reaction. The possibility to observe the trace of quasi-bound state in the mass spectra was studied. Various types of chiral based and phenomenological potentials are used to describe the interaction. It was shown that not only we can see the signature of the quasi-bound state in the mass spectra, but also, one can see the trace of branch points in the observables.
pacs:
13.75.Jz, 14.20.Pt, 21.85.+d, 25.80.Nv
I Introduction
The wish to define a precise interaction model for interaction is a basic goal in strangeness nuclear physics. For the past two decades, an enormous amount of efforts has been made to study the structure of dense kaonic nuclear clusters akaishi ; yamazaki1 ; dot1 ; dot2 . An important kaonic cluster is the system, which is a highly controversies issue in studying the kaonic systems. Many theoretical calculations were performed, focusing on the system shev1 ; shev2 ; ikeda1 ; ikeda2 ; dote1 ; dote2 ; ikeda3 ; kh1 ; kh2 ; kh3 ; maeda . All few-body calculations have shown that the is bound, but with some variation in the values of the extracted pole energy.
If system is indeed bound, then the remaining question is whether this state is sufficiently narrow to allow observation and identification. Due to strong absorption of antikaon by the nucleus, the quasi-bound state in system can have a large width. Thus, this may provide difficulties for direct experimental observation of kaonic bound states in nuclei koike . Many experimental efforts have been also performed to explore the pole structure of the system. An exclusive analysis of the reaction at 2.85 GeV yamazaki2 indicated a large peak both in the invariant-mass and missing-mass spectra, which had been predicted in a theoretical works yamazaki3 ; yamazaki4 . The observed peak corresponds to the binding energy of about 103 MeV and the width is given as 118 MeV. The quasi-bound state can be produced in kaon-induced reactions on light nuclei such as and deuteron, and the signal of the resonance may be observed in the mass spectra of the final particles. The investigation for the quasi-bound state have been further explored through incident reaction by E27 experiment at J-PARC e27 . The experiment also revealed a distinct peak in the missing-mass spectrum, nearly at the same mass and width as the DISTO peak X. The investigation for the quasi-bound state could be reached through the reaction (see Fig. 1). This reaction was performed as an E15 experiment at J-PARC e15 . The E15 group suggests a broad quasi-bound state structure at 15 MeV just below threshold e15 .
The kaon-induced reactions have been studied by Koike-Harada harada and Yamagata et al. yamagata using the optical potential method. In Ref. ohnishi the has been studied by using Faddeev approach. In this calculation they employed chiral based potentials for the -wave interaction. Within their model, they have found a clear signal corresponding to the strange-dibaryon resonances in the Faddeev scattering amplitudes and the transition probabilities. The reaction was also studied by Sekihara et al., sek1 ; sek2 ; sek3 to investigate the origin of observed peak close to the threshold in first run of E15 experiment at J-PARC e15 . Two scenarios were considered to produce the peak. In the first scenario, the (1405) resonance can be generated but it does not correlate with , and the uncorrelated system subsequently decays into and in the other one, the quasi-bound state should be generated and decays into . From the calculation of the invariant mass spectrum, the experimental signal was reproduced in the second scenario and they definitely discarded the scenario that the (1405) does not correlate with .
This study is devoted to study the pole structure of the three-body system. We study how well the signal quasi-bound state can be observed in the mass spectra resulting from reaction under consideration. We performed few-body calculations for the system by using coupled-channels Faddeev AGS equations. The transition probabilities for the reaction are calculated. With this method, we investigated the behavior of the transition probability for reaction. Several chiral based and phenomenological potentials are used shev3 ; shev4 ; ohnishi to investigate the sensitivity of the the three-body observables on two-body inputs.
The paper is organized as follows: in Sect. II, we will explain the formalism used for the coupled channel three-body system and give a brief description of the transition probability formula for break-up reactions. Sect. III is devoted to the two-body inputs of the calculations. The computed transition probabilities are presented in Sect. IV and in Section V, we give conclusions.
II Three particle system
The calculation of the three body system is based on the Faddeev treatment alt . Separable potentials were used for describing the two-body interactions
[TABLE]
where is the form factor of the interacting two-body subsystem , with relative momentum and isospin . Here, is the strength parameter of the interaction. To take the coupling directly into account, the potentials are further labeled with the values. The two-body -matrices in separable form can be given by
[TABLE]
where and are the total energy of the system and the energy of the spectator particle in channel, respectively.
[TABLE]
the quantity , is the reduced mass, when particle in channel is spectator. The operator is also the usual two-body propagator. Using separable potential for two-body interactions, the three-body Faddeev equations shev2 in the AGS take the form
[TABLE]
After partial wave decomposition and where we assumed that only s-wave contribution will be significant in our calculations, we get the following equations
[TABLE]
Here, the operators are the transition amplitudes between Faddeev channels and particle channels shev2 . The operators are the corresponding Born terms. The inputs for the AGS system of equations (5) are two-body -matrices, embedded in the three-body Hilbert space. Faddeev partition indices are used to define the interacting pair and also the spectator particle. The Faddeev equations are modified shev2 ; ikeda1 to take the coupling directly into account. Thus, in addition to the Faddeev indices the particle indices () are also added for each state () revei .
[TABLE]
Since the total isospin of the system is , therefore, depending on the spin of the two baryons, we should treat or system. The total spin of the system remains unchanged. Therefore, the baryon spins do not enter explicitly and the operators will be labeled by isospin indices. In the system the spin component is antisymmetric, so all operators in isospin base should be symmetric.
In present Calculations, we used the quasi-particle approach to solve the Faddeev equations for bound state problem. The most important part of the quasi-particle method is the separable expansion of the scattering amplitudes in the two- and three-body systems grass ; nadro ; fonce . To find the pole position, the separable representation must be defined for the three-body amplitudes and driving terms. For this purpose, we used the Hilbert-Schmidt expansion method nadro ; fix .
[TABLE]
the form factors are taken as the eigenfunctions of the kernel of the equation (5). The separable form of the Faddeev transition amplitudes is given by
[TABLE]
where the functions obey the equation
[TABLE]
Applying Hilbert-Schmidt expansion method fix to the Faddeev equations of , the following homogeneous integral equations for are obtained
[TABLE]
To solve the homogeneous system, we should search for a complex energy at which one of the eigenvalues () of the kernel matrix becomes equal to one. We must work on the physical and unphysical energy sheet of the and channels, respectively.
Another purpose of this work is to study the possible signature of the quasi-bound state in the mass spectra from the reaction . The break-up amplitude for this reaction in terms of the Faddeev transition amplitudes can be given by afnan .
[TABLE]
where is the three-body energy. To find the two-body energy, we should reduce it by spectator particle energy . Here, is the relative momentum between the interacting pair (). The quantities are Faddeev amplitudes, which are derived from Faddeev equation (5). Since the interaction is neglected, in this equation the Faddeev transition amplitudes corresponding to system are missing. Using Eq.(10), we define the transition probability of as follows,
[TABLE]
where is given by
[TABLE]
III TWO-BODY INPUT
In this section, we give a brief survey on two-body interactions, which are the central inputs for the present few-body calculations. The interaction is dominated by the -wave (1405) resonance. Therefore, the orbital angular momentum for interaction is taken to be zero. The and interactions were also taken in state. Since, the interaction in subsystem is dominated by the -wave component. Thus, in our few-body calculations the interaction is neglected. All separable potentials in momentum representation have the form of Eq. (1).
III.1 coupled-channel system
During the past two decades, different phenomenological dalit ; akaishi ; shev3 and chiral based chi1 ; chi2 ; chi3 ; chi4 ; chi5 ; chi6 ; chi7 potentials are constructed to describe the interaction. The phenomenological models of interaction consider the resonance as a quasi-bound in system embedded in the continuum. The chiral SU(3) dynamics has also turned out to be a successful approach to describe the interaction and the (1405) resonance chi1 ; chi2 ; chi3 ; chi4 ; chi5 ; chi6 ; chi7 . At and above threshold, the phenomenological and the chiral SU(3) models of interaction produce comparable results, while for subthreshold energies their results are different. The phenomenological interactions are constructed to reproduce a single pole nature for the (1405) resonance as a quasi-bound state of the system around 1405 MeV. The coupled-channels amplitude resulting from chiral SU(3) based potentials has two poles, one of them is located around 1420 MeV chi2 . while the other pole with large width is located above the threshold. Therefore, the chiral based potentials produce a binding energy of about 15 MeV for system, which is about half the binding produced with the purely phenomenological models of interaction.
We used different models to describe the s-wave interaction, which is the most important interaction in the present three-body calculations with and coupled-channels. We considered four types of phenomenological potentials. They reproduce the one- and two-pole structure for resonance and their parameters are given in Refs. shev3 ; shev4 . The parameters are adjusted to reproduce all existing data on low-energy interaction. The potentials in Ref. shev4 are adjusted to reproduce the experimental results of the SIDDHARTA experiment bazzi . Depending on a pole structure of the , we refer these potentials as “SIDD, one-pole”and “SIDD, two-pole”potential. The parameters of the potentials in Ref. shev3 are adjusted to reproduce the experimental results of the KEK experiment kek1 ; kek2 . Depending on a pole structure of the , we refer these potentials as “KEK, one-pole”and “KEK, two-pole”potential.
Plus the phenomenological potentials, we also used two chiral based potentials, which are given in Ref. ohnishi . These chiral based potentials reproduce the elastic and inelastic cross sections for the reaction as well as the mass spectra.
III.2 interaction
In order to investigate the dependence of the invariant mass on nucleon-nucleon interaction models, we used two different potentials for interaction. The first one is a two-term separable potential gal
[TABLE]
where is negative to take into account the short range repulsion part of the interaction. The parameters of this potential are given in Ref. gal .
We also used one-term PEST potential from Ref. pest . The strength parameter of the PEST potential is equal to one and the form-factor is defined by
[TABLE]
where the parameters of the potential are given in Ref. pest . The PEST potential is not repulsive at short distances, but at low energies its phase shifts are close to the rank-two potential.
For the -wave interaction, we follow the form given in Ref. tores ,
[TABLE]
where the form factor defined by . The coupling constants, , and the range parameters are given in Ref. tores .
IV RESULTS AND DISCUSSION
Solutions of the Faddeev equations corresponding to bound states and resonance poles in the channel of the three-body system were found by applying search procedures described in Sec. II. In Table 1 and 2, the results of the present work for three-body quasi-bound state are presented. The sensitivity of the pole position to the interaction is investigated by using different potential models. In Table 1 the pole position of the quasi-bound states in the systems is calculated for phenomenological models of the interaction and in Table 2, we calculated the pole energies for energy-dependent and energy-independent chiral potentials.
The position of a quasi-bound state in the three-body problem is usually defined by solving the homogeneous integral equations (9) which comes from the separable expansion of the Faddeev amplitudes. To find the resonance energy of the system using these equations, one should search for a complex energy at which one of the eigenvalues of the kernel matrix becomes equal to one. Therefore, as one can see from Eqs. 7 and 8, the Faddeev amplitudes will have a pole at this energy.
To examine the efficiency of the separable expansion method, we used another way to find the pole position(s) without using the integration in the complex momentum plane. The signal of the bound state would be seen in the Faddeev transition amplitudes. In the present work, we studied how the signature of the quasi-bound state shows up in the three-body scattering amplitudes by using coupled-channel Faddeev AGS equations. To achieve this goal, we must solve the inhomogeneous integral equations for the amplitudes defined in Eq. (5).
Fig. 2 shows the calculated three-body scattering amplitude whose initial and final states are . The off-shell momenta and are equal, 150 MeV/c and the real and imaginary part of the three-body energy, , change from -100 MeV to 0 MeV. We used one-pole (left) and two-pole (right) version of the KEK potential for describing the interaction. Since the input energy of the AGS equations is complex the moving singularities which are caused by the open channel , will not appear in the three-body amplitudes. The calculated resonance energies of the system by this method, have presented in Table 3. Comparing these results with those in Table 1, one can see that both results are in good agreement with each other.
When at least one of the intermediate particles is unstable, plus the signal of the resonance states, one can see a branch point in the complex plane. In Fig. 2, plus the signature of pole position, we can see the branch points i.e., a threshold opening associated with the pole, situated at , sum of nucleon mass and pole position. In the second row of Fig. 2, we have shown the branch points for intermediate state. The branch point is clearly visible, together with the cut that in this picture is chosen in the positive Re direction. In the second row, to make the branch points more visible the imaginary part of the three-body energy, , was chosen to be between -50 MeV and -12 MeV and the real part change from 2347 MeV to 2375 MeV. We used one-pole (left) and two-pole (right) version of the KEK potential to extract the branch points.
We investigated the dependence of the two- and three-body pole energy trajectories on the magnitude , when the strength parameter is increased from its physical value. Let stand for an enhancement factor of strength of the interaction:
[TABLE]
We calculated the pole trajectory for one- and two-pole version of the KEK potential. The behavior of the two-body and three-body pole energy trajectories are quite different at threshold. The pole energies obtained for three-body system are shown in Fig. 3 (B). The blue dashed and black solid curves in Fig. 3 (A) correspond to the two-body results. The numbers attached to the circles and squares give the corresponding values of the enhancement factor . As increases, the binding energy of the system increases for both the two- and three-body systems. In the two-body calculations of two-pole potential, the imaginary part of the resonance energy becomes smaller as the binding energy increases and for 1.33 (at threshold) the resonance almost becomes a bound state in channel. In contrary, in the three-body system the resonance energy will have a non zero imaginary part at the threshold as grows, since the channel is included effectively.
IV.1 Calculation of the transition probability
The calculated resonance energies that have presented in Table 1 and 2, give only pole positions of the system. However, we know that these results are not a quantity that can be directly measured in any experiments. To examine the existence of the quasi-bound state in system by experiments, one has to calculate the cross sections of production reactions. We can use the calculated results in Table 1 and 2 and also in Fig. 3 as guideline to study these reactions. As it was said in Sect. I, the quasi-bound state can be produced through kaon-induced reactions on light nuclei such as and deuteron. The trace of the resonances would be seen in the mass spectra of the final particles. In the present calculations, we studied how good the signature of the system shows up in the observables of the three-body reactions by using coupled-channel Faddeev equations in the AGS form. To achieve this goal, we must solve the coupled integral equations for the amplitudes defined in Eq. (5), and then construct the breakup amplitudes defined in Eq. (10). Since the kernel of AGS equations has the standard moving singularities that are caused by the opened channel and are encountered in any three-body breakup problem, we have followed the same procedure implemented in Refs. poin1 ; poin2 . Using the so called “point-method ”, we computed the cross section of reaction and studied the behavior of mass spectra. The transition probabilities for phenomenological potentials are depicted in Figs. 4 and 5. In Fig. 6 the three-body calculations are performed by chiral based potentials for interaction and PEST potential for interaction.
In Fig. 4, we calculated the mass spectra using one-pole (A) and two-pole (B) version of KEK potentials for interaction given in Ref. shev3 . To investigate the energy dependence of the transition probability, we calculated for MeV/c. We investigated the dependence of mass spectra on two-body interactions, necessary for the description of the system. Therefore, in Fig. 5, we calculated the mass spectra using the one- and two-pole version of the SIDDHARTA potential for interaction given in Ref. shev3 . The results suggest that a distinct peak of bound kaonic states should be observed, regardless of the momentum value and the class of the interaction. In the calculated mass spectra for the two-pole model of the KEK and SIDDHARTA potentials, the second pole of resonance with its large width, does not affect the invariant mass. As one can see from Figs. 4 and 5, all potential models will produce the mass spectra with the similar behavior and two distinct bump structures can be seen in the invariant mass.
The results of the full coupled-channel calculations of the scattering using two versions of the energy-dependent and energy-independent potentials derived based on chiral SU(3) dynamic and non relativistic kinematics are shown in Fig. 6. It is seen, that the three-body results corresponding to each version of interaction differ sufficiently. Therefore, in principle, it would be possible to favor one version of the potential by comparing with experimental results. Within this model, we have found two bump structures appearing in the transition probabilities in the energy region around the pole position and . As it was said before, the second bump which is situated at is actually originated from a branch point in the complex plane (see Fig.2), i.e., a threshold opening associated with the pole.
To show that, these bumps are really corresponding to the quasi-bound state in the system and pole and are not caused by threshold effects. Let us investigate these bump structures in the mass spectra and clarify the origin of these bumps. In Fig. 7, we calculated the mass spectra using one- and two-pole version of the KEK and also energy-dependent chiral potentials for interaction when the magnitude of the strength parameter is increased from its physical value. We calculated the mass spectra for three values of the enhancement factor 1.0, 1.05 and 1.10. As it was shown in Fig. 3 and Table 4, when we increase the parameter, the binding energy of the system will increase for both the two- and three-body systems and the pole energies will go toward the and threshold, respectively. Comparing the results of the mass spectra with those presented in Table 4 and Fig. 3, one can see that the bump structures in the mass spectra and the quasi-bound states in Table 4 will locate at the same energy and have the same movement. Therefore, one can say that the first bump structure should be corresponding to a quasi-bound state in system and the second bump structure is derived from a branch point in the complex plane.
In order to compare the present results with those in Ref. ohnishi , we calculated for MeV/c using the same and a two-term type potential gal with a repulsive core and an intermediate-range attraction is used to describe the nucleon-nucleon interaction. The invariant mass obtained with the two-terms are shown in Fig. 8. Energy-dependent set of potential was used together.
In contrast to the results of Ref. ohnishi , our results show that, plus the bump related to the quasi-bound state in systems, a typical bump structure manifests itself in the invariant mass at the energy related to the quasi-bound state in system. However, this bump structure in the observables dose not derives from a resonance pole and the origin of this structure is the branch points. The behavior of the mass spectra is similar to the extracted results for the phenomenological potentials. The difference between the present results and those by Ohnishi et al., can be important. In our results, we have two bump structure close to each other in the mass spectra, one is related to the quasi-bound state in system and the other corresponding to branch points which originates from intermediate . Thus, this effect should be taken into account in theoretical interpretation of the experimental results by E15.
IV.2 Averaged transition probability
In Subsection IV.1, we calculated the transition probabilities for four discrete values of momentum, but, in actual situation the momentum can occupy any value over a continuous range. To include all these momenta into consideration, in this subsection, we calculated the averaged transition probability , which is given by
[TABLE]
where the function can be defined by
[TABLE]
To define the wave function of the three-body system, we used the so called ”exact optical” potential shev3 . Therefore the channel has not included directly into account and we can drop the particle channel indices. Using Faddeev equations in 9, we can calculate three Faddeev components, , which are given by
[TABLE]
where is the reduced mass of the interacting pairs and also and are the Jacobi momenta of the spectator particle and interacting particles, respectively. The three-body wave function can be defined by
[TABLE]
where is the normalized wave function of the system, which is defined as a sum of the above components. Figure 9 (up) shows the momentum distribution of the spectator nucleon, , in system for various models of the potential and Figure 9 (down) shows the same distributions but multiplied by .
[H]
Fig. 10 shows the calculated mass spectra using one- and two-pole version of KEK (black curves) and SIDD potential (blue curves) for interaction. As one can see, the mass spectra around the threshold are affected by the kinematical effects and the peaks corresponding to the branch point and quasi-bound states are not as clear as in Figs. 4 and 5. According to the Figs. 4, 5 and 9, these changes were expected, because the threshold effects are stronger for low values of the momentum and the weight of them in the momentum distribution are larger than high momentum.
In general, Faddeev equations need as input a potential that describes the interaction between two individual particles. It is also possible to introduce a term in the equation in order to take three-body forces into account. Although, we think, that while our information about the two-body interaction is not completed, the inclusion of three-body forces can not be necessary. One can also investigate the dependence of the mass spectra on the two-body local potentials chi7 . The three-body theory of reactions can also be formulated for local potentials on the basis of the Faddeev equations. The work in this direction is underway and will be reported elsewhere.
V CONCLUSION
In summary, in this work exact Faddeev-type calculations of system were performed to define the binding energy and width of system. The efficiency of the so called HSE method was investigated. We have calculated the transition probability (11) for reaction in the energy region between the and thresholds. We have examined how the signature of the quasi-bound state in the three-body system manifests itself in the transition probabilities on the real energy axis. To investigate the dependence of the resulting transition probabilities on models of interaction, several versions of potentials, which can produce different structures for (1405) resonance, were used. Within this model, we have found a bump produced by system appearing in the transition probabilities in the energy region around the pole position. We found, that we can find a distinct peak in the mass spectrum for momentum MeV/c. In the present calculations, we also found that the shape and position of the peaks in the transition probability are independent of the momentum of the initial channel. Therefore, this fact implies that the bumps are corresponding to the (1405) and quasi-bound states. Since, the nucleon in the initial state covers a continuous, we should include the effect of all momenta in transition amplitude. We calculated the averaged transition probability . Furthermore, we have shown that not only we can see the signature of the quasi-bound state, but also, we can see the effect of the branch points in the Faddeev amplitudes (complex plane) and transition probabilities which are resulting from (1405) pole. we have shown that the bump structures related to the the branch points can affect the peak corresponding to the quasi-bound state. Thus, in the mesonic decay channel, we should consider the effect of the branch points and this reaction would also be helpful to reveal the dynamical origin of (1405) resonance.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1(1) Y. Akaishi and T. Yamazaki, Phys. Rev. C, 2002, 65 : 044005.
- 2(2) T. Yamazaki and Y. Akaishi, Phys. Lett. B, 2002, 535 : 70.
- 3(3) A. Dote, H. Horiuchi, Y. Akaishi and T. Yamazaki, Phys. Lett. B, 2004, 590 : 51.
- 4(4) A. Dote, H. Horiuchi, Y. Akaishi and T. Yamazaki, Phys. Rev. C, 2004, 70 : 044313.
- 5(5) N. V. Shevchenko, A. Gal and J. Mares, Phys. Rev. Lett., 2007, 98 : 082301.
- 6(6) N. V. Shevchenko, A. Gal, J. Mares and J. Revai, Phys. Rev. C, 2007, 76 : 044004.
- 7(7) Y. Ikeda and T. Sato, Phys. Rev. C, 2007, 76 : 035203.
- 8(8) Y. Ikeda and T. Sato, Phys. Rev. C, 2009, 79 : 035201.
