Analyzing $\Xi(1620)$ in the molecule picture in the Bethe-Salpeter equation approach

In this work, we assume that the observed state $\Xi(1620)$ is a $s$-wave $\Lambda\bar{K}$ or $\Sigma\bar{K}$ bound state. Based on this molecule picture, we establish the Bethe-Salpeter equations for $\Xi(1620)$ in the ladder and instantaneous approximations. We solve the Bethe-Salpeter equations for the $\Lambda\bar{K}$ and $\Sigma\bar{K}$ systems numerically and find that the $\Xi(1620)$ can be explained as $\Lambda\bar{K}$ and $\Sigma\bar{K}$ bound states with $J^P=1/2^-$, respectively. Then we calculate the decay widths of $\Xi(1620)\rightarrow\Xi\pi$ in these two different molecule pictures systems, respectively.


I. INTRODUCTION
The advent of the LHCb, Belle, BESIII, and other facilities and their unexpectedly successful contributions to hadron physics have stimulated of hadron studies.With the observations of some states which do not agree well with the theoretical predictions in the constituent quark model (like Λ(1405), Ξ(1620), X, Y , Z states, and pentaquark states (P c (4380) and P c (4450) states)) [1], it is important to study these unusual states, both to probe the limitations of the quark model and to discover the unrevealed aspects of the quantum chromodynamics (QCD) description of structures of hadron resonances.
Up to now lots of nucleons and S = ±1 hyperon resonances have been discovered and their quantum numbers have also been measured.In the charmed baryon sector, there have been also lots of significant progresses made in the experimental studies by the LHCb, Belle, BESIII and other collaborations.For Ξ states, only the spin-parity quantum numbers of the ground octet state Ξ(1320), the decuplet state Ξ(1530), and the excited state Ξ(1820) have been determined, but for other known Ξ resonances, their spin-parity numbers are incomplete.Foe example, the Ξ(1690) and Ξ(1620) states are cataloged in the Particle Data Group (PDG) with only one-star and three-star [1], respectively.If Ξ(1620) has J P as 1/2 − , it will be similar to the Λ(1405) state, which has been postulated as a meson-baryon molecular state or a pentaquark candidate [2].Determining the masses and quantum numbers of the Ξ resonances is vital for us to understand their structures.
The Ξ(1620) was observed through the Ξ(1620) → Ξπ decay in the 1970's [3,4].Although the mass and the width measurements in the two experiments are consistent, they both have large statistical uncertainties.Recently, the Belle Collaboration reported the observation of Ξ(1620) via its decay to [5].The mass and width are measured to be 1610.4± 6.0(stat) +5.9 −3.5 (syst) MeV and 59.9 ± 4.8(stat) +2.8 −3.0 (syst) MeV, respectively.On the theoretical side, one was shown that it is very difficult to accommodate the Ξ(1620) in the quark models [6,7].On the other hand, the meson-baryon scattering in the strangeness S = −2 sector was also studied in different unitary coupled-channel approaches constrained by QCD chiral symmetry [8][9][10].In all these chiral unitary approaches, the Ξ(1620) is dynamically generated with a relatively large decay width, and couples strongly to the Ξπ and Λ K channels but very weakly to Σ K and Ξη.In addition, the poles of Ξ(1620) are below the threshold of Λ K.
The purpose of this paper is to study the possibilities that the Ξ(1620) is a Λ K or Σ K bound state with quantum numbers J P = 1/2 − in the Bethe-Salpeter equation approach.We will also calculate the decay widths of Ξ(1620) → Ξπ in these two pictures.The Bethe-Salpeter equation is a formally exact equation to describe the relativistic bound state [11,12], and has been applied in many theoretical studies concerning heavy mesons and heavy baryons [13][14][15][16][17][18].In this paper, we will study the s-wave baryon-meson molecular bound state with the kernel introduced by the vector meson exchange interactions.This paper is organized as follows.In the next section, we will briefly review the Bethe-Salpeter equation for the bound state of a meson and a baryon.In Sec.III, we will discuss the normalization condition of the Bethe-Salpeter wave function.In Sec.V, the decay of Ξ(1620) → Ξπ will be calculated.
The numerical results will be presented in Sec.IV.In the last section, we will give a summary.

II. THE BETHE-SALPETER FORMALISM FOR THE Ξ(1620)
In this section, we will review the general formalism of the Bethe-Salpeter equation and derive the Bethe-Salpeter equation for the system composed of a baryon (Λ or Σ) and a pseudoscalar meson ( K).
Then we will derive the normalization condition for the Bethe-Salpeter wave function in the next section.
Firstly, we define the Bethe-Salpeter wave function for the bound state |P of a baryon (Λ or Σ) and a pseudoscalar meson ( K) as the following: where ψ(x 1 ) and φ(x 2 ) are the field operators of the baryon (Λ or Σ) and pseudoscalar meson ( K) at space coordinates x 1 and x 2 , respectively, P denotes the total momentum of the bound state with mass M and velocity v.In momentum space, the Bethe-Salpeter wave function can be defined as where p represents the relative momentum of the two constituents.
The Bethe-Salpeter equation for the bound state can be written in the following form: where S ψ (p 1 ) and S φ(p 2 ) are the propagators of the baryon (Λ or Σ) and the pseudoscalar meson ( K), respectively, and K(P, p, q) is the kernel which contains two-particle-irreducible diagrams.For convenience, we define p l (= p • v) and p µ t (= p µ − p l v µ ) to be the longitudinal and transverse projections of the relative momentum (p) along the bound state momentum (P ).Then, the propagator of Λ ( or Σ) has the form and the propagator of the K meson can be expressed as where ω 1(2) = m 2 1(2) + p 2 t (in which we have defined is the spinor of the bound state with helicity s), χ P (p) can be written as where g i (i = 1, • • •, 5) are Lorentz-scalar functions.Furthermore, each term in the expansion of χ P (p) transforms exactly in the way that χ P (p) transforms under P -parity and Lorentz transformations, which can help us simplify the form of χ P (p), it is easy to prove that χ P (p) can be simplified as in which f 1 (p) and f 2 (p) are two independent Lorentz-scalar function of p.
As discussed in the introduction, we will study the s-wave bound state of the Λ K and Σ K systems.The isospin field doublets ψ = ψ 0 , ψ − T and φ = −φ + , φ 0 T have the following expansions in momentum space: The isospin quantum number of Ξ(1620) is 1/2, so the flavor wave function of Λ K and Σ K systems can be written as Projecting the bound states on the field operators ψ 1 (x), ψ 2 (x), φ 1 (x), and φ 2 (x), then we have where χ I P is the common Bethe-Salpeter wave function for the bound state with isospin for the Λ K system, and the isospin coefficients are for the Σ K system.
Then the Bethe-Salpeter equation for the bound state can be written as where i(j) and l(k) refer to the components of the ψ(φ) field doublets.Then, the Bethe-Salpeter equation for the I = 1/2 Λ K molecule can be written as and for the I = 1/2 Σ K molecule the Bethe-Salpeter equation can be write as In the Bethe-Salpeter equation approach, the interactions between Λ and K mesons are due to the light vector-meson (ω and φ) exchanges.There is no ρ-exchange contribution, because of the isospin conservation.For the Σ K interaction we will consider the exchanges of vector mesons ρ, ω and φ.The pseudoscalar meson exchanges are forbidden because the K meson is also a pseudoscalar meson.The Lagrangians for the vertices of the strange K meson and one-strange baryon with vector mesons are [19,20]: where c.c. is the complex conjugate of the first term, τ is the Pauli spin matrix.The coupling constants are constrained by the SU (3) symmetry, g KKρ = g KKω = g ρππ /2 and g KKφ = g ρππ / √ 2. The ρππ coupling is determined by , where we take the value α = 1.15 based on the ω coupling constant given in Ref. [21].g ΣΣρ = g ΣΣω = 2αg N N ρ and g ΣΣφ = − √ 2(2α − 1)g N N ρ .g N N ρ is chosen as g ρππ /2 as in Ref. [19,20], Under SU (3) symmetry, the κ BBV (B = Λ, Σ) can be obtained with the relations where f BBV is defined as f BBρ = g BBρ κ BBρ , and κ BBρ = 6.1 and f N N ω = 0 [19].
From the above observations, at the tree level, in the t-channel we have the following kernel for the Bethe-Salpeter equation in the so-called ladder approximation: where m V represents the mass of the exchanged vector meson (ρ, ω and φ), and c I is the isospin coefficient: 1, 1 for ρ, ω, φ mesons, respectively.
In order to describe the phenomena in the real world, we should include a form factor at each interacting vertex of hadrons to include the finite-size effects of these hadrons.For the meson-exchange case, the form factor is assumed to take the following form [22]: where Λ, m and k represent the cutoff parameter, the mass of the exchanged meson and the momentum of the exchanged meson, respectively.
After reducing the above coupled integral equations for f1 (p t ) and f2 (p t ) to one dimensional integral equations, we obtain the following equations: Substituting Eq. ( 3) into above equations, then, one can derive the parametric forms of α P (p) and f1 (q t ). (35) After substituting Eqs. ( 34) and (35) into Eq.( 32), we have

IV. NUMERICAL RESULTS FOR THE BETHE-SALPETER WAVE FUNCTIONS
In this part, we will solve the Bethe-Salpeter equation numerically and try to search for possible solutions of the Λ K and Σ K bound states.It can be seen from Eq. ( 25) that there is only one free parameter in our model, the cutoff Λ, which contains the information about the nonpoint interactions due to the structure of hadrons at the interaction vertices.Although the value of Λ cannot be exactly determined and depends on the specific process, it should be typically the scale of low energy physics, which is about 1 GeV.In this work, we treat the cutoff in the form factors as a parameter varying in a much wider range 0.8 − 4.8 GeV.
To find out the possible molecule bound states, one only needs to solve the homogeneous Bethe-Salpeter equations.One numerical solution of the homogeneous Bethe-Salpeter equation corresponds to a possible bound state.The integration region in each integral will be discretized into n pieces, with n being sufficiently large.In this way, the integral equation will be converted into an n×n nmatrix equation, and the scalar wave functions of each equation will now be regarded as an n-dimensional vector.Then, the two coupled integral equations can be illustrated as where f1( 2) is an n-dimensional vector, and where ǫ is a parameter introduced to avoid divergence in numerical calculations, w and y are parameters used in controlling the slope of wave functions and finding the proper solutions for these functions.Then one can obtain the numerical results of the Bethe-Salpeter wave functions by requiring the eigenvalue of the eigenvalue equation to be 1.
In the rest frame, we define p ′ 1 = (E ′ 1 , p ′ ) and p ′ 2 = (E ′ 2 , −p ′ ) (p ′ is three-momentum) to be the momenta of Ξ and π, respectively.The masses of Ξ and π are m ′ 1 and m ′ 2 , respectively.According to the kinematics in the rest frame of the two-body decay, one has and where |p ′ | is the norm of the 3-momentum of either particle in the final state in the rest frame of the initial bound state and M is the Lorentz-invariant decay amplitude of the process.
According to the above interactions, the decay Ξ(1620) → Ξ − π + is shown in Fig. 3.We can write down the amplitude as where q ′ = λ 2 q ′ 1 − λ 1 q ′ 2 is not the relative momentum of particles in the final state, λ 1 and λ 2 are defined as λ i = m i /(m 1 + m 2 ), and m 1 and m 2 are the masses of the component particles of the bound states but not of the final state.
Then, we apply the numerical solution of the Bethe-Salpeter amplitude to calculate the decay width of Ξ(1620) → Ξπ.The decay widths are 36.94MeV and 9.35 MeV for the Λ K and Σ K bound stats, respectively.

VI. SUMMARY
In this paper, we applied the Bethe-Salpeter equations to study the possibilities that the is Ξ(1620) is s-wave Λ K or Σ K bound states with the quantum numbers J P = 1/2 − .Considering the interaction kernels based on ω and φ mesons exchange diagrams for the Λ K system and ρ, ω, and φ mesons exchange diagrams for the Σ K system, we established the Bethe-Salpeter equations in the ladder and instantaneous approximations.Because the constituent particles and the exchanged particles in the Λ K and Σ K systems are not pointlike, we introduced a form factor including a cutoff Λ which reflects the effects of the structure of these particles.Since Λ is controlled by nonperturbative QCD and cannot be determined at present, we let it vary in a reasonable range within which we examined whether Λ K and Σ K bound states could be the Ξ(1620) state by solving the Bethe-Salpeter equations.From our calculations, we found that the Ξ(1620) state can be treated as the Λ K and Σ K bound states when Λ = 1632 MeV and 1356 MeV, respectively.
Then, we applied the numerical solutions of the Bethe-Salpeter wave functions to calculate the decay widths of Ξ(1620) → Ξπ which are induced by K * exchange meson.We obtained that the decay widths are 36.94MeV and 9.35 MeV for the Λ K and Σ K bound states, respectively.Comparing the magnitides of these two decay widths, it is obvious that the Ξ(1620) has a larger contribution from the Λ K system than the Σ K system.The same conclusion was also found in Ref. [8] from the chiral perturbation theory.
Clearly, more theoretial and experimental efforts will be needed to fully understand the nature of the one-star Ξ(1620).
m 1 and m 2 are the masses of Λ(Σ) and K mesons In general, for a baryon and a pseudoscalar meson bound state, considering

FIG. 2 :
FIG.2: Numerical result for the Bethe-Salpeter wave functions in the Σ K system.
which corresponds to the matrix labeled by p t and q t in each integral equation.Generally, |p t | (and |q t |) varies from 0 to +∞.Here, |p t | (and |q t |) will be transformed into a new variable t that varies from −1 to 1 based on the Gaussian integration method,