Investigation of Baryons in the Hypercentral Quark Model

In the present study, we consider baryons as three-body bound systems according to the hypercentral constituent quark model in configuration space and solve the three-body Klein-Gordon equation. Then, we analyze perturbative spin-dependent and isospindependent interaction effects. To find the analytical solution, we use screened potential and calculate the eigenfunctions and eigenvalues of some baryons. We consider exclusive semileptonic decays of bottom and charm baryons and apply the differential decay width with the Isgur-Wise function and arrive at the rates for some semileptonic baryon decays. The results prove more enhanced compared to recent works and comply well with the experimental data.


Introduction
Exact solutions have been of paramount importance in the study of relativistic and nonrelativistic quantum systems as they can yield useful information on these systems. Nevertheless, there are so few exact solutions that approximate methods are often used to investigate quantum systems. For instance, Schrödinger, Klein-Gordon, and Bethe-Salpeter equations have been solved through proper approximation [1][2][3] to the centrifugal term for a number of potentials like Eckart [4], Manning-Rosen [5], Pusch-Terrel [6], Hulthen [7], and Wei et al. [8] employed an approximation to arrive at the approximative analytical l-wave scattering solutions of the Schrödinger equation. Their findings closely comply with other studies as for short-range potentials but differ as for long-range ones. They also solved the Klein-Gordon equation by the unequal scalar and vector Hulthen potential via a creative approximation to the centrifugal term [9].
In this study, the baryons are treated as spinindependent three-quark systems and based on hyperspherical approach; a nonrelativistic quark model is presented for them. Considering the significance of the analytical solution for the qualitative analysis of the spectrum, we present analytical formalisms for wave functions and energy levels by using the screened potential. In the 1980s, a screening mul-tiplicative factor in the potential was presented [10,11]. Accordingly for larger distances, the screened linear term in the potential acts as the dominant term, so the most excited states are supposed as confinement states, which are determined by considering merely the screened linear interaction. Also for short distances, due to the small linear term, the Coulomb term is dominant. Therefore, screened potential accounts for the behavior of quarks according to the QCD rules [12]. Finally, we make use of spin-isospin and flavor-dependent interactions to identify the baryon energy spectrum more accurately.
In the next section of this paper, we examine the heavyto-light semileptonic decays of heavy baryons. These decays are utilized to determine the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. For this purpose, we use the Isgur-Wise function (IWF). The IWF has provided physicists with remarkable results for transitions of bottom baryons into charm baryons [13]. The calculation of the IWF yields useful information on the decay width, branching ratio, and CKM matrix [10,14,15]. Ivanov et al. studied the semileptonic decays of bottom and charm baryons within the relativistic quark model [16]. Recent measurements and the pertinent theoretical approaches contribute greatly to understanding the decays of bottom baryons. A number of scholars have made attempts to study the structure of heavy-baryon transitions via IWF formalism. Ke et al. [17] investigated baryon weak decays in the light-front quark model. Using the heavy diquark model for baryons, Viet described the matrix elements of weak currents of baryons through universal IWF [18]. Ebert et al. [19] used the relativistic quark model to calculate the semileptonic decay rates of heavy baryons. Cardarelli and Simula explored the corresponding form factors of heavy baryons in a lightfront constituent quark model [20]. Körner et al. also studied the heavy baryon transitions and analyzed the decay rates, polarization effects, lifetimes and angular decay distributions [21].
In our previous works, we solved the Schrödinger equation for different potentials for few-quark systems [22][23][24][25]. In a more recent work, we investigated the relativistic Klein-Gordon equation analytically for the Deng-Fan potential and Hulthen plus Eckart potential under the equal vector and scalar potential conditions [26]. In this paper, we study our three-body system and make use of the hyperspherical coordinates and calculate the parameters of the IWF for some baryons. The purpose of this study is to investigate the heavy-baryon decays in hyperspherical coordinates. To this end, we use the hyperradial Klein-Gordon equation in configuration space in Section 2. Afterwards, in Section 3, we solve the Klein-Gordon equation via screened potential. Then, by analyzing the IWF corresponding to hyperspherical coordinates, we study the heavy-baryon decays and calculate the decay width of these particles in Section 4. And finally, the relevant conclusions are presented in Section 5.

The Three-Body Problem
Calculations of mass and energy of some heavy baryons have been carried out through a simple semirelativistic quark model. Assuming a dominant SU(6)-invariant spectrum in the nonrelativistic hypercentral constituent quark model which led to the apparent success in describing the mass spectra and other observables motivates us to introduce a spin-independent relativistic description to approximate the problem and then deal the SU(6)-breaking interactions perturbatively.
In this study, we investigate baryons using a semirelativistic model based on the hypercentral quark constituent model. Due to the significance of the relativistic analysis of particles behavior within bound systems, we study baryons as three-quark systems while taking relativistic effects into account. Totally relativistic study of baryons demands solving the Dirac equation for each quark (while considering interquark interactions). Totally relativistic solution to a three-body problem is difficult. Thus, we have to apply a few simplifications to solve this problem. As the first step to simplify the particles' spin and isospin, we introduce spin and isospin effects perturbatively. The reason for the assumption is the experimental baryon spectra has a SU(6) symmetry for baryon resonances. Therefore, for the SU(6)invariance part, we can considerations as a three-quark system of zero-spin and use a separate relativistic description of spin for the baryon spectrum. The exact analytical formulas for energy levels and wave functions have been obtained, which can be used for a parametric analysis of the problem.
The three-body Klein-Gordon equation with equal scalar and vector potentials in the natural unit h = c = 1 is of this form [9]: where jψi denotes the total wave function of the three-body system and 3ε is its total energy. V ij stands for the potential of interaction between the i-th and j-th particles. We solve the three-body spherical coordinates system and express the spatial degrees of freedom in terms of Jacobian momentums. Jacobi coordinates to define three-body basis states are given by [27,28] Here r i (i = 1, 2, 3) denote coordinate of the i-th constituent quark and m is the constituent quark mass. In the numerical calculations we take the values of the quark masses as m s = 0:483 GeV, m u = m d = 0:336 GeV, m b = 4:950 GeV, and m c = 1:550 GeV. In hyperspherical coordinates, the hyperradius and hyperangle are defined as in Figure 1.
In addition to x and φ, four polar angles Ω ρl and Ω λl show the orientation of Jacobian vectors according to reference experimental frameworks. For Jacobian momentums, In hyperspherical coordinates in the D-dimensional space, the Laplace operator for N particles generally becomes The term L 2 ðΩÞ/x 2 , as a generalization of the centrifugal barrier for the six-dimensional space, involves the angular coordinates Ω ρ and Ω λ and the hyperangle ξ. In a three-2 Advances in High Energy Physics body system, where λ is the polynomial order and λ = 2n + l ρ + l λ . l ρ ðl λ Þ where n is a positive integer and l ρ and l λ denote the quantum numbers of angular momentums related to λ and ρ coordinates. According to the hypercentral hypothesis, the interacting potential in equation (1) can be expressed as a function of hyperradius x only: The potential VðxÞ is known as hypercentral, which is invariant for any rotation in D-dimensional space spanned by the coordinates and OðDÞ symmetry. Depending on the hyperradius x implies that the potential has a three-body character as the dependence on the single-pair coordinates cannot be disentangled from the other particles. Accordingly, the following three-body Klein-Gordon equation, can be reduced to a six-dimensional time-dependent Schrodinger-like equation: where

Screened Potential
Screened potential is of the following form [12]: The hat over the parameters in the potential distinguishes them from the nonscreened case; σ corresponds to the string tension of the confinement in screened case. α s is a parameter representing the Coulomb strength, and μ stands for a screening parameter. By using the following approximation to the first order, This potential turns out as By using variable changes s = ð1/xÞ, β = 1/x 0 , and z = s − β, 3 Advances in High Energy Physics By simplifying the above equation, we arrive at x 0 is supposed as a characteristic radius of the baryon. The screened potential appears as which can be reformulated as where A, B, and C appear as the following.
By using Nikiforov-Uvarov (NU) method (see appendix), we calculate the values for constant parameters α i s, as we calculated these parameters in [26]. They are shown in Table 1. By introducing these values into equation (18), By introducing the values of A, B, and C into equation (17), And new variables are defined as With the above variable changes and simple calculations employed, the energy eigenvalues will be As α 3 = 0, then to calculate the wave function, the Logger polynomial is used. Therefore, the wave function will be Here, γ, the degree of the polynomial, is called the grand orbital quantum number and is equal to 2n + l ρ + l γ and υ = 0, 1, 2, ⋯ according to the hypercentral hypothesis. Also, g = ð2/3Þα s ðm + εÞ/ðν +γ + 5/2Þ and N ν,γ is a normalization coefficient. The baryon masses are obtained via Similarly, we can repeat this iteration for other modes; we obtain the energy and the corresponding wave function for any arbitrary υ. We have plotted the distribution of hyperradial wave function density for Ω − , Ξ 0 , Ξ + in Figure 2. As expected, the wave function peak height grows larger as the baryon mass increases. By obtaining energy, we calculate the mass of baryon according to the following equation: where E 0 is calculated for the ground state. H int is the perturbation Hamiltonian and depends on the spin, isospin, and flavor degrees of freedom. Table 2 displays the fitted values for calculating the energy eigenvalue.
Advances in High Energy Physics Table 3 presents the results obtained for spinless ground state of some baryons. These results are consistent with particle data group [29]. Even though the terms of tensor and spin-orbit are small, to achieve higher accuracy, these terms are considered as perturbation potentials and the perturbation energy levels are calculated. The perturbation method is helpful when the eigenstates cannot be obtained easily.
To obtain the first-order correction of the energy eigenstates, we calculate the expected value of the perturbed potential while the system is in the unperturbed state.
where H int ðrÞ = H S 1 ,S 2 ðrÞ + H I 1 ,I 2 ðrÞ + H S,I ðrÞ. This interaction for a particle by total spin S and isospin T is of the following form: where A S = 35:46, A I = 33:49, A SI = −3:94 (all in MeV·fm) and S is spin and T is isospin of the baryon. From equation (28), it can be deduced that H int in baryons is nonzero. Only for particles with S = T = 1:5 and S = T = 4:4 is this perturbative energy of order 1 zero, and such particles do not physically exist. Table 4 presents our results as compared to the results from the nonrelativistic method [30] and lattice QCD calculations [31,32]. As observed, there is a good agreement between the two sets of results. References [33,34] also determine the spectra of c quark and b quark baryons, and they use the color coulomb plus power and the hyper-Coulomb plus linear potentials, respectively, to solve six-dimensional Schrödinger equation numerically using the Mathematica method.

Isgur-Wise Function and Semileptonic Decays of Baryons
The IWF may be parameterized in different ways [17,35,36]. In the approximation of infinite heavy quark masses, the form factors F i , G i ði = 1, 2, 3Þ define the semileptonic transition of Λ b related to a particular universal IWF. These form factors appear as However, all the models similarly do this at the zero recoil point as Taylor's series as where ρ is the magnitude of the slope and c is the curvature (convexity parameter) of the Isgur-Wise function at ω = 1. The higher terms in the IWF may be negligible since slope and curvature are the most dominant parameters. The Isgur-Wise function for the weak decay of heavy baryons transition in the HCQM can be written as an overlap integral of the baryon wave functions and is of the following form: Generally, the overlap integral which involves the final and the initial wave function is used to calculate transition matrix elements. In the above equation, only jψ 0,0 ðxÞj 2 comes into the picture instead of the overlap integral of the final and the initial state. This is because we have where p 2 = 2m 2 ðω − 1Þ. p 2 stands for the square of virtual momentum transfer. Then, x Figure 3: Differential decay width for c baryons.    (30) and (33), we arrive at the slope (ρ 2 ) and curvature (c) parameters as The Isgur-Wise function mentioned by equation (35) depends on the product of two terms: the first is the square of the modulus of the wave function (jψ ν,γ ðxÞj 2 ) and the second is cos ðpxÞ. The value of the cosine term appeared in the Isgur-Wise function becomes 1, when we put ω = 1. And the remaining term (jψ ν,γ ðxÞj 2 ) gives ξðω = 1Þ, while integrating it for the ground state ðν = 0, γ = 0Þ wave function at the zero recoil point. The obtained slope and curvature for Ξ + c , and Ω c baryons appear in Table 5. The IWF presented in this section is the key to calculating the decay width for heavy baryon transitions. Other works also have employed this method, but they have used different interaction potentials. For instance, [37] uses Cornell potential, [38] uses a Killingbeck potential. Reference [39] makes use of the Hulthen plus linear to describe baryon interactions. However, our findings are well consistent with theirs. In Figure 3, we illustrate the behavior of dΓ/Adω in terms of ω for a few baryons in two graphs. For example, for semileptonic decay Λ b ⟶ Λ c transition, we have calculated differential decay width by where where η = m Λ b /m Λ c and BðΛ c ⟶ abÞ is the branching ratio for the decay Λ c ⟶ að1/2 + Þ + bð0 − Þ and |V cb | stands for the Kobayashi-Maskawa matrix element. Also m 1 and m 2 refer to the masses of Λ b and Λ c , respectively. To calculate the total decay width, we have integrated equation (36) over the ω as The maximum of the ω values is calculated from Figures 3 and 4 depict the variations of the differential decay width for a few heavy-baryon decays. As expected, the heavier the baryon, the larger the parameters ρ 2 and c, and the more quickly dΓ/Adω tends to zero. The semileptonic decay of Λ 0 b ⟶ Λ + c transitions is shown in Figure 5. In this transition, b quark converts to c quark and a virtual W [40], which becomes a charged lepton and neutrino. Using the decay rates of baryons in Table 6 and integrating equation (36), we obtain function Γ and calculate its error. To identify the error of equation (38), we applied the natural logarithm to both sides of equation (38) and then derived from it. By considering each term's error, we determined the absolute error of Γ. Also from the obtained decay width and decay lifetime of each baryon adopted from PDG, we found the branching ratio values for semileptonic decay. The values obtained for Λ 0 b and Λ + c comply with PDG values. The branching ratio for other decay channels is not given in PDG, and we indeed propose branching ratio values and present them in the last column of

Conclusions
In this study, we formulated a simple three-body quark model for some baryons. The present model, therefore, has a main SU(6)-invariant part and a perturbative SU(6)-breaking interaction for the splitting within the multiplets. Such a structure may also emerge from lattice QCD calculations. Thus, first we switched off spin and isospin degrees of freedom for the SU(6)-invariant part of the spectrum and considered the baryon as a spin-independent three-body bound system and suggested a semirelativistic quark model for baryons based on the spin-independent three-particle Klein-Gordon equation. We calculated the observable properties of semileptonic decays of bottom and charm baryons, i.e., Isgur-Wise functions, decay rates, and distributions too. Finally, we found the branching ratio values for semileptonic decay. The computed results are compared with other theoretical predictions as well as with the available experimental observations. It can pave the way for further studies, including double heavy quark baryons and their decays. Such stud-ies are important as the current interest in heavy-flavor hadron physics is growing.