This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Spectroscopic properties of Δ\Delta Baryons

Chandni Menapara    Zalak Shah [email protected]    and Ajay Kumar Rai [email protected]    [email protected] Department of Applied Physics, Sardar Vallabhbhai National Institute of Technology, Surat-395007, Gujarat, India
Abstract

The resonance state of Δ\Delta baryon existing in four isospin (I=32I=\frac{3}{2}) states, has been studied using Hypercentral Constituent Quark Model (hCQM) with a simple linear potential with added first order correction. The calculated data range for 1S-5S, 1P-5P, 1D-4D and 1F-2F with possible spin-parity assignments of all the states. The magnetic moments have also been obtained for all four configuration. The NπN\pi decay channel width has been calculated for few states. The linear nature of the data has been verified through Regge trajectories.

keywords:
Mass spectra, light baryon, magnetic moment, regge trajectory
00footnotetext: ©\scriptstyle\copyright2013 Chinese Physical Society and the Institute of High Energy Physics of the Chinese Academy of Sciences and the Institute of Modern Physics of the Chinese Academy of Sciences and IOP Publishing Ltd

1 Introduction

Hadron spectroscopy is a tool to reveal the dynamics of the quark interactions within the composite systems like baryons, mesons, and exotics. The phenomenological approach of hadron spectroscopy is about using the potential to establish the resonance mass of higher radial and orbital states of a hadron. Also, the various possible decays of a resonance state help in identifying the short-lived hadrons and even missing excited states. A number of resonance states of light and heavy hadrons have been provided by Particle Data Group [2].

The specific target, here, is the study of Δ\Delta baryon, a member of Baryon decuplet (JP=32+J^{P}=\frac{3}{2}^{+}) which is composed of light quarks u and d. Inspite of bieng the same composition of Nucleons The possible four combinations of the symmetric wave function gives four Δ\Delta particle with isospin I=32I=\frac{3}{2} as Δ++\Delta^{++} (uuu, I3=32I_{3}=\frac{3}{2}), Δ+\Delta^{+} (uud, I3=12I_{3}=\frac{1}{2}), Δ0\Delta^{0} (udd, I3=12I_{3}=-\frac{1}{2}) and Δ\Delta^{-} (ddd, I3=32I_{3}=-\frac{3}{2}). The present work is motivated by the fact that heavy quark systems decay into light quark systems through various decay channels and most of the matter is composed of these light quark systems. Δ\Delta(1232) has been observed experimentally in pion-nucleon decays for quite long [3, 4], recent studies have still been exploring the new properties at HADES-GSI [5]. Δ\Deltas, likely an excited state of nucleon (N) with ground state 939 MeV have been extensively studied through photoproduction decays by ELSA[6]. However, the symmetric flavour wavefunction of Δ\Delta differs from mixed symmetry wavefunction of nucleons. Thus, the revealing every known and unknown property of Δ\Delta baryon has always been a matter of interest as discussed in many review articles [7, 8, 9, 10]. The Δ\Delta resonances shall also be focused on at upcoming experimental facilities at PANDA-GSI [11, 12].

Phenomenological and theoretical models for light baryon studies have been developed and modified over time. The light baryon resonances have been explored through well-known Isgur-Karl model basically applied for P-wave states [13] as well as modified with relativised approach [14], Goldstone-boson exchange model due to spontaneous chiral symmetry breaking [15, 16], quark-diquark system along with Gursey-Radicati exchange intercation [17, 18], and semi-relativistic model with SU(6)-invariant and SU(6)-violating terms [19]. Lately, varied approaches based on QCD SUM Rules [20], basis light-front model [21] and light-front relativistic [22], Lattice QCD [23] and covariant Faddeev approach [24] and others based on n and JPJ^{P} values and the respective trajectories against square of mass of a given state. [25, 26]. The spectrum of octet and decuplet light baryons has also been studied in a relativistic approach using instanton induced quark forces [27].

In this paper, a non-relativistic hypercentral Contituent Quark Model (hCQM) has been employed to obtain resonance masses of radial and orbital states of Δ\Delta baryon [28, 29, 30]. The potential term consists of two parts: a Coulomb-like term and a Confinement term. A similar methodology has been employed for heavy baryons using different potential such as screened potential [31], linear[32, 33], etc.

The paper is organized as follows: after introduction, theoretical framework has been discussed. The third section incorporates the results and discussion of the mass spectra. Sections four, five and six deal with baryon magnetic moment, Regge trajectory and decay widths respectively. Finally, conclusions are drawn in the last section.

2 Hypercentral Constituent Quark Model (hCQM)

Hadron spectroscopy is useful for better understanding of hadron as a bound state of quarks and gluons as well as the spectrum and internal structure of excited baryons. This is a key to strong interactions in the region of quark confinement. The system becomes complex and difficult to deal considering all the interaction of quark-quark, quark-gluon and gluon-gluon. This is the reason for using constituent quark mass incorporating all the other effects in the form of some parameters.
A Constituent Quark Model is a modelization of a baryon as a system of three quarks or anti-quarks bound by some kind of confining interaction. An effective way to study three body systems is through consideration of Jacobi coordinates as

ρ=12(𝐫𝟏𝐫𝟐)\displaystyle{\bf\rho}=\frac{1}{\sqrt{2}}({\bf r_{1}}-{\bf r_{2}}) (1a)
λ=(𝐦𝟏𝐫𝟏+𝐦𝟐𝐫𝟐(𝐦𝟏+𝐦𝟐)𝐫𝟑)𝐦𝟏𝟐+𝐦𝟐𝟐+(𝐦𝟏+𝐦𝟐)𝟐\displaystyle{\bf{\bf\lambda}=\frac{(m_{1}{\bf r_{1}}+m_{2}{\bf r_{2}}-(m_{1}+m_{2}){\bf r_{3}})}{\sqrt{m_{1}^{2}+m_{2}^{2}+(m_{1}+m_{2})^{2}}}} (1b)
x=ρ2+λ𝟐;ξ=arctan(ρλ)x=\sqrt{\rho^{2}}+{\bf\lambda^{2}};\;\;\xi=arctan(\frac{\rho}{\lambda}) (2)

where x is hyperradius and ξ\xi is hyperangle.
The Hamiltonion of the system is expressed as

H=P22m+V0(x)+1mxV1(x)+VSD(x)H=\frac{P^{2}}{2m}+V^{0}(x)+\frac{1}{m_{x}}V^{1}(x)+V_{SD}(x) (3)

where m=2mρmλmρ+mλm=\frac{2m_{\rho}m_{\lambda}}{m_{\rho}+m_{\lambda}} being the reduced mass.

The dynamics are considered in the wave-function ψ(x)\psi(x) which is the solution of hyper-radial equation

[d2dx2+5xddxγ(γ+4)x2]ψ(x)=2m[EV(x)]ψ(x)\left[\frac{d^{2}}{dx^{2}}+\frac{5}{x}\frac{d}{dx}-\frac{\gamma(\gamma+4)}{x^{2}}\right]\psi(x)=-2m[E-V(x)]\psi(x) (4)

The potential incorporated solely depends on hyperradius x of the system and not on hyperangle [34].

V0(x)=τx+αxνV^{0}(x)=-\frac{\tau}{x}+\alpha x^{\nu} (5)

V(x) consists of Coulomb-like term and a confining term which is taken to be linear with power index ν\nu=1. Another part of potential form is the first order correction term with 1mx=(1mρ1mλ)\frac{1}{m_{x}}=(\frac{1}{m_{\rho}}-\frac{1}{m_{\lambda}}).

V1(x)=CFCAαs24x2V^{1}(x)=-C_{F}C_{A}\frac{\alpha_{s}^{2}}{4x^{2}} (6)

where CFC_{F} and CAC_{A} are Casimir elements of fundamental and adjoint representation. αs\alpha_{s} is the running coupling constant.

Alongwith zeroth and first order correction term in hypercentral approximation, spin-dependent term VSDV_{SD}(x) is also incorporated to sharply distinguish the degenerate states[35].

VSD(x)=VSS(x)(𝐒ρ𝐒λ)+VγS(x)(γ𝐒)+VT×[S23(𝐒𝐱)(𝐒𝐱)x2]\begin{split}V_{SD}(x)=V_{SS}(x)({\bf S_{\rho}\cdot S_{\lambda}})+V_{\gamma S}(x)({\bf\gamma\cdot S})\\ +V_{T}\times[S^{2}-\frac{3({\bf S\cdot x})({\bf S\cdot x})}{x^{2}}]\end{split} (7)

where VSS(x)V_{SS}(x), VγS(x)V_{\gamma S}(x) and VT(x)V_{T}(x) are spin-spin, spin-orbit and tensor terms respectively.

The quark masses are taken as mu=md=0.290m_{u}=m_{d}=0.290 GeV. The numerical solution of the six-dimensional Schrodinger equation has been performed using Mathematica Notebook [36].

3 Results and Discussion

Based on the model and potential term discussed in the above section, the resonance masses from 1S-5S, 1P-5P, 1D-4D and 1F-2F with allowed spin-parity assignments have been computed in Table 3. In addition, the present results are compared with different models inspired resonance masses for available states.

\tabcaption

Resonance masses of Δ\Delta baryons (in MeV). State JPJ^{P} Present model PDG[2] Status [16] [17] [18] [19] [28] [26] [13] [14] [25] 1S 32+\frac{3}{2}^{+} 1232 1230-1234 **** 1232 1235 1247 1231 1232 1232 1232 1230 1232 2S 32+\frac{3}{2}^{+} 1611 1500-1640 **** 1659.1 1714 1689 1658 1727 1625 1600 3S 32+\frac{3}{2}^{+} 1934 1870-1970 *** 2090.2 1914 1921 1920 4S 32+\frac{3}{2}^{+} 2256 - - 5S 32+\frac{3}{2}^{+} 2579 - - 1P 12\frac{1}{2}^{-} 1609 1590-1630 **** 1667.2 1673 1830 1737 1573 1645 1685 1555 1P 32\frac{3}{2}^{-} 1593 1690-1730 **** 1667.2 1673 1830 1737 1573 1720 1685 1620 1P 52\frac{5}{2}^{-} 1550 - - 2P 12\frac{1}{2}^{-} 1956 1840-1920 *** 2003 1910 1910 1900 2P 32\frac{3}{2}^{-} 1919 1940-2060 ** 1910 1940 2P 52\frac{5}{2}^{-} 1871 1900-2000 *** 2003 1910 1908 1945 3P 12\frac{1}{2}^{-} 2280 - * 2150 3P 32\frac{3}{2}^{-} 2242 - - 3P 52\frac{5}{2}^{-} 2193 - - 4P 12\frac{1}{2}^{-} 2602 - - 4P 32\frac{3}{2}^{-} 2565 - - 4P 52\frac{5}{2}^{-} 2515 - - 5P 12\frac{1}{2}^{-} 2926 - - 5P 32\frac{3}{2}^{-} 2888 - - 5P 52\frac{5}{2}^{-} 2836 - - 1D 12+\frac{1}{2}^{+} 1905 1850-1950 **** 1873.5 1930 1827 1891 1953 1895 1910 1D 32+\frac{3}{2}^{+} 1868 1870-1970 *** 1930 2042 1935 1920 1D 52+\frac{5}{2}^{+} 1818 1855-1910 **** 1873.5 1930 2042 1891 1901 1895 1905 1D 72+\frac{7}{2}^{+} 1756 1915-1950 **** 1873.5 1930 2042 1891 1955 1950 1950 2D 12+\frac{1}{2}^{+} 2227 - - 2D 32+\frac{3}{2}^{+} 2190 - - 2D 52+\frac{5}{2}^{+} 2140 - ** 2200 2D 72+\frac{7}{2}^{+} 2078 - - 3D 12+\frac{1}{2}^{+} 2556 - - 3D 32+\frac{3}{2}^{+} 2516 - - 3D 52+\frac{5}{2}^{+} 2463 - - 3D 72+\frac{7}{2}^{+} 2397 - - 4D 12+\frac{1}{2}^{+} 2874 - - 4D 32+\frac{3}{2}^{+} 2835 - - 4D 52+\frac{5}{2}^{+} 2784 - - 4D 72+\frac{7}{2}^{+} 2720 - - 1F 32\frac{3}{2}^{-} 2165 - - 1F 52\frac{5}{2}^{-} 2108 - - 1F 72\frac{7}{2}^{-} 2037 2150-2250 *** 2200 1F 92\frac{9}{2}^{-} 1952 - - 2F 32\frac{3}{2}^{-} 2486 - - 2F 52\frac{5}{2}^{-} 2430 - * 2350 2F 72\frac{7}{2}^{-} 2359 - - 2F 92\frac{9}{2}^{-} 2274 - ** 2400

The four star status assigned by Particle Data Group(PDG) ensures the certainty of its existence with quite known properties. The radial states comprise of JP=32+J^{P}=\frac{3}{2}^{+}, the 2S(1600) predicted as 1611 MeV differs by 11 MeV from Ref.[25], 14 MeV from Ref.[26] and nearly 47 MeV Refs.[16, 19]. Similarly the 3S(1920) as 1934 falls within PDG range and differs only by 1-14 MeV from some references.

The first orbital excited state 1P(1620) with 12\frac{1}{2}^{-} is well within the range of PDG and differs by 36 MeV from Ref[26]. However, the 1P(1700) state predicted with 1593 (32\frac{3}{2}^{-}) is underpredicted by 97 MeV from lower range of experimental data. The three star states of 2P with spin-parity assignment 12\frac{1}{2}^{-} and 52\frac{5}{2}^{-} are over- and under-predicted compared to PDG ranges.

The four star designated second orbital state 1D with 12+\frac{1}{2}^{+} is obtained as 1905 MeV differs by 5 MeV from ref[25] and 25 MeV from ref[17]. The two states with (52+,72+\frac{5}{2}^{+},\frac{7}{2}^{+}) have been predicted to be quite low compared to known data as well as other references. Based on current results, Δ(1920)\Delta(1920) state from PDG might be assigned to 1D(32+\frac{3}{2}^{+}) or 3S(32+\frac{3}{2}^{+}) based on the comparison shown in the table.

The predicted 1F 72\frac{7}{2}^{-} state as 2037 MeV is 113 MeV less than lower limit of PDG-range. However, present study has attempted to predict many unknown states too which are least explored by any other models and experiments.

4 Baryon Magnetic Moment

Baryon magnetic moment places a crucial role in providing information regarding the structure and shape of baryon [37]. Magnetic moment of Δ++\Delta^{++} has been precisely measured through pion bremstrahlung analysis [38]. Theoretically, magentic moment of JP=32+J^{P}=\frac{3}{2}^{+} decuplet baryons have been calculated through different approaches which include Quark Model and QCD Sum Rules [39], Chiral Quark Model [40, 41], color dielectric model [42]. However, all of the calculations do not account for complicated effects due to valence quarks, pion cloud, exchange currents, constituent quark mass, etc and thus neglected. In the present study, effective quark mass has been considered to obtain magnetic moment of all four Δ\Delta isospin states. Baryon magentic moment is expressed as [43]

μB=qϕsf|μqz|ϕsf\mu_{B}=\sum_{q}\left\langle\phi_{sf}|\mu_{qz}|\phi_{sf}\right\rangle (8)

where ϕsf\phi_{sf} is the spin-flavour wave function.

μqz=eq2mqeffσqz\mu_{qz}=\frac{e_{q}}{2m_{q}^{eff}}\sigma_{qz} (9)

The effective mass of quark mqeffm_{q}^{eff} would be different from the model based mass as within baryon, mass may vary due to interaction among quarks.

mqeff=mq(1+Hqmq)m_{q}^{eff}=m_{q}(1+\frac{\langle H\rangle}{\sum_{q}m_{q}}) (10)

where H=E+Vspin\langle H\rangle=E+\langle V_{spin}\rangle[43]. A similar study for NN^{*} has been done by Zalak Shah et al [34].

The result in Table 4 shows that the Δ\Delta magnetic moments obtained from present work are in quite accordance with experimental results. The Ref. [40] has compared magnetic moment using different sets of data; so based on that non-relativistic quark model Δ++\Delta^{++} magnetic moment is 5.43μN\mu_{N} differing by 0.86 μN\mu_{N}. The magnetic moment for Δ+\Delta^{+} and Δ\Delta^{-} is 2.72μN\mu_{N} and -2.72μN\mu_{N} respectively which differs by 0.48μN\mu_{N}.

\tabcaption

Magnetic Moment of Δ\Delta(1232) isospin state State Wave-function μ\mu exp(PDG)[38] Δ++\Delta^{++}00 3μu\mu_{u}0 4.568μN\mu_{N}0 4.52μN\mu_{N} Δ+\Delta^{+}00 2μu+μd\mu_{u}+\mu_{d}0 2.28μN\mu_{N}0 2.7μN\mu_{N} Δ0\Delta^{0}00 2μd+μu\mu_{d}+\mu_{u}0 00 - Δ\Delta^{-}00 3μd\mu_{d}0 -2.28μN\mu_{N}0 -

5 Regge Trajectory

An important property concluded from baryon spectrum is the plot of J, total angular momentum against M2M^{2} as well as principle quantum number n against M2M^{2}. These lines are so far observed to be linear and non-intersecting for light baryon spectrum [44]. These plots provide a confirmation between experimental and theoretical predicted masses of excited state with their respective quantum numbers[45]. This holds true for positive and negative parity states as well. Regge trajectories have been widely employed in heavy hadron studies too [46, 47]. The equations are as follows

J=αM2+α0\displaystyle J=\alpha M^{2}+\alpha_{0} (11a)
n=βM2+β1\displaystyle n=\beta M^{2}+\beta_{1} (11b)
Refer to caption
Figure 1: (n,M2M^{2}) Regge trajectory for Δ\Delta states
Refer to caption
Figure 2: (J,M2M^{2}) Regge trajectory for Δ\Delta states

The trajectory in Fig.1 based on equation (10b) shows that calculated data has been in good agreement with the nature of experimental data as all the calculated resonance squared mass fall on linear curve. Also, few individual experimental points are marked in the graph agree with the total angular momentum and spin configuration assigned in the calculated data.
The plot of total angular momentum quantum number J with natural parity P against the squared mass is represented in Fig. 2 also follows the linear curve.

6 Decay Widths

The observations of decays of baryon resonances afford a valuable guidance in assigning the resonances their correct places in various symmetry schemes. The correct isotopic spin assignment is likely to be implied by the experimental branching ratio into different charge states of particles produced by the decay, while experimental decay widths provide a means of extracting phenomenological coupling constants.

The chiral quark model, in which constituent quarks couple directly to mesons, is known to describe the properties of the ground state octet and decuplet baryons quite well [48].

The prominant decay channel for Nucleons including Δ\Delta has been observed to be NN^{*} and pion, depending on the charge of respective parent [49]. The transition couplings of vector mesons has been obtained along with other constants by Riska et al. [50]. In the present work, the constants and decay widths provided by Particle Data group has been employed to establish the decay width of some well-established resonance mass. For Δ(1600)\Delta(1600) decay to NπN\pi,

Γ=13f24πE+mNmk3mπ2\Gamma=\frac{1}{3}\frac{f^{2}}{4\pi}\frac{E^{{}^{\prime}}+m_{N}}{m^{*}}\frac{k^{3}}{m_{\pi}^{2}} (12)

where, EE^{{}^{\prime}} is the energy of the final nucleon and k is pion momentum.

E=m2mπ2+mN22mE^{{}^{\prime}}=\frac{m^{*2}-m_{\pi}^{2}+m_{N}^{2}}{2m^{*}} (13)
k=[m2(mN+mπ)2][m2(mNmπ)2]2mk=\frac{\sqrt{[m^{*2}-(m_{N}+m_{\pi})^{2}][m^{*2}-(m_{N}-m_{\pi})^{2}]}}{2m^{*}} (14)

Here mm^{*} is resonance mass calculated using above model, mNm_{N} is nucleon mass 939 MeV and mπm_{\pi} is pion mass 139 MeV. Using m=1611m^{*}=1611 and f=0.51f=0.51, Γ=24.8%\Gamma=24.8\% which is well within the PDG range 824%8-24\%.
For Δ(1620)\Delta(1620) decaying to NπN\pi,

Γ=f24πE+mNmkmπ2(mmN)2\Gamma=\frac{f^{2}}{4\pi}\frac{E^{{}^{\prime}}+m_{N}}{m^{*}}\frac{k}{m_{\pi}^{2}}(m^{*}-m_{N})^{2} (15)

m=1609m^{*}=1609, f=0.34f=0.34, Γ=92%\Gamma=92\% whereas PDG range is 2535%25-35\%.

For Δ(1700)\Delta(1700) decaying to NπN\pi,

Γ=13f24πEmNmk3mπ2\Gamma=\frac{1}{3}\frac{f^{2}}{4\pi}\frac{E^{{}^{\prime}}-m_{N}}{m^{*}}\frac{k^{3}}{m_{\pi}^{2}} (16)

m=1593m^{*}=1593, f=1.31f=1.31, Γ=14.83%\Gamma=14.83\% whereas PDG range is 1020%10-20\%.

7 Conclusion

In the present work, Δ\Delta resonance masses have been calculated using Hypercentral Constituent Quark Model employed with linear potential. Also, the first order correction has been included. All the masses upto 2F states have been compared with available experimental data as well as different theoretical and phenomenological models in Table 3. Therein Δ(1232)\Delta(1232), Δ(1600)\Delta(1600), Δ(1620)\Delta(1620), Δ(1700)\Delta(1700), Δ(1905)\Delta(1905), Δ(1910)\Delta(1910) and Δ(1950)\Delta(1950) - four star states; Δ(1900)\Delta(1900), Δ(1920)\Delta(1920), Δ(1930)\Delta(1930) and Δ(2200)\Delta(2200) - three star states and other fairly established states have been predicted.

It is evident that radial excited states as well as orbital excited states with lower spin state agree to a considerable level with PDG-range and few of the models vividly discussed in section 3. However, higher spin states of orbital excited states are mostly under predicted compared to experimental range.

The Regge trajectories have been plotted with principle quantum number n and angular momentum J against square of resonance mass. Fig. 1 shows that Regge trajectories are linear but not exactly parallel. However, experimental points are not very far from the respective lines. Fig. 2 resolves that the spin-parity assignment for orbital excited states also follow the linear relation.

The baryon magentic moment has been calculated for all four isospin state of Δ\Delta as described in Table 4, however values of two isospin states are not obtained experimentally so far. The Δ++\Delta^{++} magnetic moment is almost similar to PDG value and Δ+\Delta^{+} magnetic moment differs by 0.58μN\mu_{N} from that of PDG.

Finally, decay widths have been obtained for strong decay through NπN\pi channel for three states Δ(1600)\Delta(1600), Δ(1620)\Delta(1620) and Δ(1700)\Delta(1700) using the nucleon to vector meson transition couplings. For Δ(1600)\Delta(1600) and Δ(1700)\Delta(1700) decay width are well withing the range but Δ(1620)\Delta(1620) decay width predicted is higher than the experimental range.

Thus, present work has effectively explored the known and unknown properties of Δ\Delta baryon in a similar approach of earlier N* spectroscopy[34]. The accomplishments and shortcomings from this study is expected to inspire for improvements and further exploring other light baryons in addition to experimental facilities PANDA-GSI[11, 12].

Acknowledgement

One of the authors, Ms. Chandni Menapara would like to acknowledge the support from Department of Science and Technology (DST) under INSPIRE-FELLOWSHIP scheme.

References

  • [1]
  • [2] P.A. Zyla et al., (Particle Data Group), Prog. Theor. Exp. Phys. 2020, 083C01 (2020)
  • [3] A. Barnicha, G. Lopez Castro and J. Pestieau, Nucl. Phys. A 597, 623-635 (1996)
  • [4] E. Pedroni et al., Nucl. Phys. A 300 321-347 (1978), DOI:10.1016/0375-9474(78)96136-5
  • [5] J. Adamczewski-Musch et al. (HADES Collaboration) Phys. Rev. C 95, 065205 (2017)
  • [6] U. Thoma et al. (CB-ELSA Collaboration), (2007), arXiv: 0707.3592v4 [hep-ph]
  • [7] V. Crede and W. Roberts, Rept. Prog. Phys., 76: 076301 (2013), arXiv:1302.7299 [hep-ph]
  • [8] M. M. Giannini, E. Santopinto, Chin. J. Phys., 53, 020301 (2015)
  • [9] N. P. Samios, M. Goldberg, and B. T. Meadows, Rev. Mod. Phys., 46, 49 (1974)
  • [10] A. Valcarce, H. Garcilazo, F. Fernandez, and P. Gonzalez, Rept. Prog. Phys., 68, 965 (2005), arXiv:hep-ph/0502173 [hep-ph]
  • [11] G.Barucca et al., (PANDA Collaboration), Eur. Phys. J. A 55, 42 (2019)
  • [12] B. Singh et al., (PANDA Collaboration), Phys. Rev. D 95, 032003 (2017)
  • [13] N. Isgur and G. Karl, Phys. Rev. D. 18, 4187 (1978)
  • [14] S. Capstick and N. Isgur, Phys. Rev. D 34, 2809 (1986)
  • [15] L. Ya. Glozman, Z. Papp and W. Plessas, Phys. Lett. B 381, 311 (1996)
  • [16] Z. Ghalenovi and M. Moazzen, Eur. Phys. J. Plus 132, 354 (2017)
  • [17] E. Santopinto, Phys. Rev. C 72, 022201(R) (2005)
  • [18] E. Santopinto and J. Ferretti, Phys. Rev. C 92, 025202 (2015)
  • [19] M. Aslanzadeh and A. A. Rajabi, Int.J.Mod.Phys. E26, 1750042 (2017)
  • [20] K. Azizi, Y. Sarac and H. Sundu, J.Phys.G 47, 9, 095005 (2020) arXiv:1909.03323v2 [hep-ph]
  • [21] J.P.Vary et al., Few Body Syst. 59, 56 (2018), arXiv:1804.07865 [nucl-th]
  • [22] I. G. Aznauryan and V. D. Burkert, Phys. Rev. C 85, 055202 (2012), arXiv:1201.5759 [hep-ph]
  • [23] R. G. Edwards, J. J. Dudek, D. G. Richards, and S. J. Wallace Phys. Rev. D 84, 074508 (2011)
  • [24] H. Sanchis-Alepuz and C. S. Fischer, Phys. Rev. D 90, no. 9, 096001 (2014)
  • [25] Y. Chen, B.Q. Ma, Chin. Phys. Lett. 25, 3920 (2008)
  • [26] E. Klempt, Phys. Rev. C 66, 058201 (2002)
  • [27] U. Loring, B. C. Metsch and H. R. Petry, Eur. Phys. J. A 10, 395-446 (2001)
  • [28] M. M. Giannini, E. Santopinto and A. Vassallo, Eur. Phys. J. A 12, 447-452 (2001)
  • [29] C. Menapara, Z. Shah and A. K. Rai, AIP Conf.Proc. 2220, 140014 (2020)
  • [30] C. Menapara and A. K. Rai, DAE Symp.Nucl.Phys. 64 673-674 (2020)
  • [31] K. Gandhi and A. K. Rai, Eur. Phys. J. Plus, 135, 213 (2020)
  • [32] K. Gandhi, Z.Shah and A. K. Rai, Int. J. Ther. Phys. 59, 1129 (2020)
  • [33] Z. Shah and A. K. Rai, Few-Body Syst, 59, 112 (2018)
  • [34] Z. Shah, K. Gandhi and A. K. Rai, Chin. Phys. C 43, 034102 (2019), arXiv:1812.04858 [hep-ph]
  • [35] M. B. Voloshin, Prog. Part. Nucl. Phys. 61, 455-511 (2008)
  • [36] W. Lucha, F. schoberls, Int. J. Modern Phys. C. 10 607 (1997)
  • [37] A. J. Buchmann, Few-Body Syst. 59, 145 (2018), arXiv:1902.10166v1 [hep-ph]
  • [38] A. Bossard et al., Phys. Rev. D 44 (1991)
  • [39] T. M. Aliev and V. S. Zamiralov, Adv. High Energy Phys. 2015,406875 (2015)
  • [40] H. Daliya and M. Gupta, Phys. Rev.D 67, 114015 (2003)
  • [41] H. Daliya and M. Gupta, Int. J. Mod. Phys. A 19, 5027-5042 (2004)
  • [42] S. Sahu, Rev. Mex. Fis. 48, 48-51 (2002)
  • [43] B. Patel, A. Majethiya and P. C. Vinodkumar, Pramana 72, 679-688 (2009)
  • [44] J. A. Silva Castro et al., Phys. Rev. D. 99, 034003 (2019)
  • [45] P. Masjuan and E. R. Arriola, Phys. Rev. D 96, 054006 (2017)
  • [46] Z. Shah, K. Thakkar and A. K. Rai, Eur. Phys. J. C 76, 530 (2016)
  • [47] Z. Shah and A. K. Rai, Eur. Phys. J. A 53, 195 (2017)
  • [48] L.Ya. Glozman, D.O. Riska, Phys. Rep. 268 263-303 (1996)
  • [49] R. Bijker, J. Ferretti, G. Galata, H. Garcia-Tecocoatzi and E. Santopinto, Phys. Rev. D. 94, 074040 (2016)
  • [50] D. O. Riska and G. E. Brown, Nucl. Phys. A. 679, 577-596 (2001)