In this paper, we demonstrated that the multiple turning point problems within the framework of the Wentzel-Kramers-Brillouin (WKB) approximation method can be reduced to two turning point one for a nonsymmetric potential function by using an appropriate Pekeris-type approximation scheme. We solved the Schrödinger equation with the Killingbeck potential plus an inversely quadratic potential (KPIQP) function. The special cases of the modeled potential are discussed. We obtained the energy eigenvalues and the mass spectra of the heavy QQ¯ and heavy-light Qq¯ mesons systems. The results in this present work are in good agreement with the results obtained by other analytical methods and available experimental data in the literature.
1. Introduction
In the Quark model, a meson is a hadronic subparticle which consists of the quark and its antiquark, and the reduced mass is dominated by the light quark mass [1]. The meson system is mediated by the strong interactions which are described by the theory of quantum chromo-dynamics (QCD) [2]. The heavier mesons also known as quarkonia QQ¯ are the constituents of the heavy quarks such as the bottom (b) and charmed (c) and are considered as nonrelativistic bound systems described by the Schrodinger wave equation (SE) [1, 2]. However, the case is different for the light (uds) quark interaction systems qq¯ which are described by relativistic equations [1, 2]. Also, the heavy-light meson Qq¯ system bound states have been studied using both the relativistic and nonrelavistic quark models [1, 2]. The bound state solutions to the wave equations under the quark-antiquark interaction potential functions such as the Killingbeck or the Cornell potentials have attracted much research interest in atomic and high energy physics [3–26]. The Killingbeck potential comprises the sum of the Cornell plus the Harmonic oscillator potentials, while the Cornell potential is the sum of the Coulomb plus linear potentials. The Cornell potential and its extended forms have been extensively solved with the SE [3–22], semirelativistic equation [23], and also with the relativistic equation [24–26]. The solution of the wave equation with some potentials are exactly solvable for l=0, whereas other potentials are insoluble and nontrivial for any arbitrary angular momentum quantum number l≠0. In this case, numerical methods and approximate analytical techniques are required to obtain the solutions of the quantum system of choice [27, 28]. Owing to the nontrivial mathematical properties of the quark-antiquark interaction potentials, several analytical techniques and approximate methods have been used to calculate the energy and mass spectra. The standard methods used in the literature are the Nikoforov Uvarov (NU) method [12–17, 29], asymptotic iterative method (AIM) [30], the Laplace transformation method (LTM) [18], artificial neural network method (ANN) [19], and the analytical exact iterative method (AEIM) [11]. The authors in [12–17] obtained the mass spectrum of the quark-antiquark interaction system using an appropriate Pekeris-type approximation scheme to deal with the orbital centrifugal energy barrier. In this present work, motivated by their approach, the same Pekeris-type approximation scheme is extended to the WKB formalism for the first time to the best of our knowledge. It is well known that the multiple turning points problems within the framework of the WKB approximation method can be reduced to the standard two turning points problem for symmetric potential functions, for example, the isotropic harmonic oscillator and the molecular pseudoharmonic potential with the proper coordinate transformation [31–33]. This approach allows us to obtain the exact bound state solutions of the SE for any arbitrary l quantum numbers in the physical axis 0<r1<r2, where r1,r2 are the classical turning points.
In this present work, we studied the nonrelativistic quark model under the interaction of the nonsymmetric KPIQP. The main focus of this paper is to obtain the mass spectra of the heavy meson such as charmonium cc¯, bottomonium bb¯, bottom-charmed bc¯, and the heavy-light meson such as the charmed-strange cs¯ system via the WKB approximation method.
The KPIQP has the form [15, 30] of
(1)Vr=Ar2+Br−Cr+Dr2,where A,B,C, and D are constant potential parameters.
If we set A=B=D=0,, the KPIQP reduces to the Coulomb potential used in the description of the hydrogenic atom.
(2)Vr=−Cr
The KPIQP reduces to the Cornell potential if we set the constants A=D=0(3)Vr=Br−Cr,B,C>0.where C is a coupling constant and B is a linear confinement parameter.
The SE spectrum generated by the Cornell potential can be used in the investigation of the masses and decay widths of charmonium states [8, 19]. The coulombic term arises from one gluon exchange between the quark and its antiquark and dominates at short distances [25]. While the linear term, which is supported by lattice QCD measurements, dominates at large distances [19]. The addition of an inversely quadratic potential (IQP) and the harmonic oscillator (HO) term to the Cornell potential improves the behavior of the potential in the region r→0. Furthermore, it leads to improved results as compared to the Cornell potential [30].
The remaining part of the paper is organized as follows. In section two, the synopsis of the WKB approximation formalism is presented. Section three contains the analytical solution of the bound states of the SE generated by the KPIQP. The special cases of the obtained energy eigenvalue are discussed. In section four, we present the numerical results of the masses of charmonium cc¯, bottomonium bb¯, bottom-charmed bc¯, and charmed-strange cs¯ mesons. Furthermore, we compared the WKB mass spectra with the ones obtained by other analytical methods and available experimental data. The paper is concluded in section five.
2. WKB Approximation Formalism
The WKB approximation method is an effective tool initially proposed to find approximate solutions to the one dimensional time independent SE in the limiting case of large radial quantum numbers. It is used to obtain the finite wave function and energy eigenvalues of potentials of interest [20, 31–43]. The method can be used to study quantum tunneling rates in a potential barrier, resonance behavior in a continuum, the exponential decay of an unstable system [43], and the quasi-normal nodes of a black hole and quantum cosmology [44]. The method fails at the classical turning point where the momentum vanishes. This failure is one of the major problem associated with the method. The difficulty can be circumvented using the connection formula [42, 43]. The method accuracy varies markedly for the ground and other low lying states depending on the potential function [39]. Also, the leading order WKB approximation scheme does not yield an exact eigenvalue of the radial SE [31]. To overcome this problem, the orbital centrifugal barrier term ll+1 in the radial SE has to be replaced with the term l+1/22. This modification is known as the Langer correction [40]. Sergeenko [38] stated that the Langer correction regularizes the WKB wave function at the origin and ensures the correct asymptotic behaviour at large radial quantum numbers. Also, the centrifugal barrier contribution of the effective potential does not vanish for the s-wave case l=0 which makes the SE nontrivial for some potential functions [31, 45], hence the need to use an appropriate approximation scheme.
The three-dimensional time-independent SE with a reduced mass μ and wave-function ψr,θ,ϕ is given as
(4)ℏ22μ1r2∂∂rr2∂∂r+1r2sinθ∂∂θsinθ∂∂θ+1r2sin2θ∂2∂ϕ2ψr,θ,ϕ+Vrψr,θ,ϕ=Eψr,θ,ϕ.
With the help of the method of separation of variables, we can obtain the radial SE by using the transformation ψr,θ,ϕ=RrYθ,ϕ/r(5)d2Rrdr2+2μℏ2E−Vr−l+1/22ℏ22μr2Rr=0where the effective potential is given as
(6)Veffr=Vr+l+1/22ℏ22μr2
We can rewrite Eq. (5) as
(7)−iℏddr2Rr=2μE−Vrl+1/22ℏ22μr2Rr.
From Eq. (7), the classical momentum is given as
(8)Pr=2μE−Vr−l+1/22ℏ22μr21/2
The standard WKB quantization condition [38, 39] for two turning points (r1,r2) problem is given as
(9)∫r1r2Prdr=πℏn+12⋅r1<r<r2n=0,1,2⋯
The turning points are gotten from Eq. (8) by setting Pr=0.
The semiclassical wave function in the leading ℏ approximation has the form
(10)ψWKBr=NPrexp±iℏ∫Prdr.
3. Solution of the Radial Schrodinger Equation
In this section, we obtain the bound state solution of the SE by substituting Eq. (1) into the WKB standard quantization condition given by Eq. (9)
(11)∫r1r22μE−Ar2−Br+Cr−Dr2−l+1/22ℏ22μr21/2dr=πℏn+12⋅r1<r<r2,n=0,1,2⋯
In this case, the classical momentum is given as
(12)Pr=2μE−Ar2−Br+Cr−Dr2−l+1/22ℏ22μr21/2
Equation (12) will produce four turning points r1,r2,r3 and r4 of which two of the points must lie in the real axis r>0 where Pr=0. The multiple turning point problem with the potential in Eq (1) is beset with some mathematical challenges. (i) The turning points of the polynomial equation in Eq. (12) would have to be determined by some means of algebra. (ii) The evaluation of the integral in Eq. (11) is not so easy to solve analytically even without the centrifugal barrier term.
To obtain the analytical solution of the WKB quantization integral in Eq. (11), we used the Pekeris-type approximation scheme by setting r=1/y.
Changing the variable from r to y, we obtained
(13)−2μ∫y1y2dyy2E−Ay2−By+Cy−Dy2−l+1/2ℏ2y22μ=πℏn+12,y1<y<y2,n=0,1,2
Next, we expanded the terms A/y2 and B/y in power series form to the second-order around r0δ=1/r0 which is assumed to be the characteristic radius of a meson [12–17].
If we let x=y−δ and expand around x=0, we obtained
(14)Ay2=Ax+δ2=Aδ21+xδ2≈Aδ21−2xδ+3x2δ2≈A6δ2−8yδ3+3y2δ4.
In the same vein, we obtained the expansion for B/y as
(15)By≈B3δ−3yδ2+y2δ3.
On substituting Eqs. (14) and (15) into (13), we obtained
(16)−2μ∫y1y2dyy2−T+My−Ny2=πℏn+12.y1<y<y2,n=0,1,2⋯.,where
(17)N=D+3Aδ4+Bδ3+l+1/22ℏ22μ(18)M=C+8Aδ3+3Bδ2(19)−T=E−6Aδ2−3Bδ
Equation (16) can further be simplified as
(20)−2μN∫y1y2dyy2−C+By−y2=πℏn+12⋅y1<y<y2,n=0,1,2⋯.,where
(21)B=MN(22)C=TN
Furthermore, we can write Eq. (20) in a more compact or standard form:
(23)−2μN∫y1y2dyy2y2−yy−y1=πℏn+12.0<y1<y2
Where y1,y2 are real turning points obtained by solving the quadratic equation in the squared root of the integrand of Eq. (20).
(24)y1=B−B2−4C2(25)y2=B+B2−4C2
Performing the integration of (23), we obtained the solution
(26)∫y1y2dyy2y2−yy−y1=πy2y1−1/2y2+y1y2y1,0<y1<y2
By comparing Eqs. (23) and (26), implies that
(27)2μN1/2y2+y1−y2y1y2y1=ℏn+12.
Furthermore, with the help of Eqs. (24) and (25), we can write the sum and products of the turning points as
(28)y2+y1=B=MN(29)y2y1=C=TNsubstituting Eqs. (28) and (29) into (27), we obtained the expression
(30)T=μ2Mℏn+1/2+2μN2
Finally, using the respective notations T,M, and N given in Eqs. ((17)–(19)), we obtained the energy eigenvalue expression for the KPIQP as
(31)Enl=−2μℏ28A/δ3+C+3B/δ22n+1+8μ/ℏ23A/δ4+B/δ3+D+l+1/22ℏ2/2μ2+6Aδ2+3Bδ
4. Discussion
In this section, we present some special cases of the energy eigenvalues of the KPIQP and also obtain the mass spectra of quarkonia and the heavy-light meson system.
If we set A=B=D=0 and C=Ze2 in Eq. (31), we immediately retrieved the Coulomb’s energy eigenvalue equation given as
(32)EnlCP=−μZ2e42ℏ2np2,np=1,2,3where Z is the atomic number, e is the electronic charge, and np is a principal quantum number with the notation np=n+l+1.
Setting A=B=0,D≠0,C≠0, Eq. (1) reduces to the molecular Kratzer-type [46–59] potential and the energy eigenvalue equation has the form of
(33)Enl=−2μℏ2C2n+1+8μ/ℏ2D+l+1/22ℏ2/2μ2
If we let A=D=0, then Eq. (31) reduces to the approximate energy expression of the Cornell potential,
(34)Enl=−2μqq¯ℏ2C+3B/δ22n+1+8μ/ℏ2B/δ3+l+1/22ℏ2/2μ2+3Bδ
In order to compute the mass spectra of meson systems, we use the meson mass relation [12]
(35)Mnl=mq+mq¯+Enl
Where mq and mq¯ are the respective quark and antiquark masses.
By substituting Eq. (34) into (35), the mass spectrum of the meson systems for any arbitrary radial and angular momentum quantum numbers becomes
(36)Mnl=mq+mq¯−2μqq¯ℏ2C+3B/δ22n+1+8μ/ℏ2B/δ3+l+1/22ℏ2/2μ2+3Bδ
Where mqq¯ is the reduced mass of the quark-antiquark systems given by the relation
(37)μqq¯=mqmq¯mq+mq¯
With the numerical support of the MAPLE package, we generated the masses of QQ¯ and Qq¯ mesons systems by fitting Eq. (36) with experimental data given in Tables 1–4. For the charmonium meson, we substituted the experimental data for the 1S, 2S, 3S states into Eq. (36) and obtained the free parameters B,C and δ by solving three algebraic equations. For the bottomonium meson, we substituted the value of B obtained from the charmonium fit, including the substitution of the experimental values for the 1S and 2S states into Eq. (36) to obtain the parameters C and δ. Also, for the bottom-charmed meson, the values of B,C gotten from the charmonium fit with the 1S state experimental data are substituted into Eq. (36) to obtain the parameter δ. Finally, the free parameters for the charmed-strange meson are obtained by inserting B, and the experimental data for the 1S and 2S states into Eq. (36). This enabled us to determine the C and δ parameters from two algebraic equations. The masses for other excited states are obtained and compared with the results obtained by other analytical methods including available experimental data. Generally, the mass spectra increase as the adjacent quantum levels increase. These trends are in consonance with the results reported in the literature [11, 15, 18, 19, 30], including available experimental data [60–62]. Also, the results for the charmed-strange meson tabulated in Table 4 are in excellent agreement with the ones obtained in Refs. [15, 30] and experimental data such as the excited 1D-state [62]. Furthermore, in Table 3, the mass spectrum of the bottom-charmed meson is very close to the ones obtained with the artificial neural network method [19] and also to the 2S state experimental data [60] indicating an improvement compared to the other methods.
Charmonium cc¯ spectra in GeVmc=1.209GeV,B=0.202GeV2,C=1.213,δ=0.236GeV.
State
Present work
AEIM [11]
NU [15]
AIM [30]
LTM [18]
ANN [19]
Exp. [60]
1S
3.098
3.095481
3.095
3.096
3.0963
3.098
3.097
2S
3.689
3.567354
3.685
3.686
3.5681
3.688
3.686
3S
4.041
4.039226
4.040
4.275
4.0400
4.029
4.039
4S
4.266
4.511098
4.262
4.865
4.5119
1P
3.262
3.567735
3.258
3.214
3.5687
3.516
3.511
2P
3.784
4.039607
3.779
3.773
4.0406
3.925
3.927
1D
3.515
4.039683
3.510
3.412
4.0407
3.779
3.770
Bottomonium bb¯ spectra in GeVmb=4.823GeV,B=0.202GeV2,C=1.664,δ=0.361GeV.
State
Present work
AEIM[11]
NU [15]
AIM [30]
LTM [18]
ANN [19]
Exp. [60]
1S
9.461
9.74473
9.460
9.460
9.745
9.460
9.460
2S
10.023
10.02315
10.022
10.023
10.023
10.026
10.023
3S
10.365
10.30158
10.360
10.585
10.3016
10.354
10.355
4S
10.588
10.58000
10.580
11.148
10.580
10.572
10.579
1P
9.608
10.02406
9.609
9.492
10.0246
9.891
9.899
2P
10.110
10.30248
10.109
10.038
10.3029
10.258
10.260
1D
9.841
10.30266
9.846
9.551
10.3032
10.156
10.164
Bottom-charmed bc mass spectra in GeVmc=1.209GeV,mb=4.823GeV2,B=0.202GeV2,C=1.213,δ=0.371GeV.
State
Present work
AEIM [11]
NU [15]
AIM [30]
LTM [18]
ANN [19]
Exp. [60]
1S
6.274
6.277473
6.277
6.277
6.2770
6.274
6.275
2S
6.845
7.037641
7.383
6.814
7.0372
6.839
6.842
3S
7.125
7.797808
7.206
7.351
7.7973
7.245
4S
7.283
7.038623
7.889
7.522
1P
6.519
7.798791
7.042
6.340
7.0381
6.743
2P
6.959
6.663
6.851
7.7983
7.187
1D
6.813
6.452
7.046
Mass spectra of charmed-strange cs meson in GeVmc=1.628GeV,ms=0.419GeV,B=0.202GeV2,C=2.046,δ=0.561GeV.
State
Present work
NU [15]
AIM [30]
Exp.
1S
1.969
1.968
2.512
1.96849±0.00033 [61]
2S
2.709
2.709
2.709
2.709 [62]
3S
2.913
2.932
2.906
4S
2.998
3.102
1P
2.601
2.565
2.649
2P
2.876
2.860
1D
2.862
2.857
2.859
2.859 [62]
5. Conclusion
The energy eigenvalue equation of the SE with the KPIQP has been obtained with elegance. We demonstrated that the multiple turning point problems can be truncated to two turning points one for a nonsymmetric potential function by using an appropriate Pekeris-type approximation scheme to the deal with the orbital centrifugal barrier energy. This approximation scheme allows us to solve the WKB integral analytically for any angular momentum quantum number. We computed the masses of mesons systems by using the Cornell potential as a model. The approximate mass spectra of the meson systems obtained in this present work are in excellent agreement with the results obtained by the other analytical methods and thus reinforce the exactness of the leading order WKB approximation method. Finally, our results may be useful to future experimental works and also can be applied in the study of molecular structures and interactions of diatomic molecules.
Data Availability
The results in this present work were numerically obtained from the theoretical calculations without taking any data. Therefore, no data is used in our manuscript.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
It is a pleasure for us to thank the kind referees for their many useful comments and suggestions which greatly helped us in making improvements to this paper.
MutukH.Mass Spectra and Decay Constants of Heavy-Light Mesons: A Case Study of QCD Sum Rules and Quark Model20182018810.1155/2018/80956532-s2.0-850564980588095653SergeenkoM. N.Light and Heavy Mesons in The Complex Mass Scheme2019https://arxiv.org/abs/1909.10511SyahroniE.SuparmiA.CariC.The determination of the spectrum energy on the model of DNA-protein interactions using WKB approximation method2017795, article 01202710.1088/1742-6596/795/1/0120272-s2.0-85014098852CariC.SuparmiA.YuniantoM.PratiwiB. N.Supersymmetric approach for Killingbeck radial potential plus noncentral potential in Schrodinger equation2016776, article 01209210.1088/1742-6596/776/1/0120922-s2.0-85012979661RoychoudhuryR. K.VarshniY. P.Shifted 1/N expansion and exact solutions for the potential Vr=−Z/r+gr+λr2198821133025303410.1088/0305-4470/21/13/0252-s2.0-0041197164Abu-ShadyM.Heavy Quarkonia and Mesons in the Cornell potential with harmonic oscillator potential in the N-dimensional Schrödinger equation2016221620CiftciH.KisogluH. F.Nonrelativistic Arbitrary l-States of Quarkonium through Asymptotic Iteration Method201820187454970510.1155/2018/45497052-s2.0-85048596690HallR. L.SaadN.Schrödinger spectrum generated by the Cornell potential2015131838910.1515/phys-2015-00122-s2.0-84923698759ChungH. S.LeeJ.KangD.Cornell potential parameters for s-wave heavy quarkonia20085241151115410.3938/jkps.52.11512-s2.0-43149085790VegaA.FloresJ.Heavy quarkonium properties from Cornell potential using variational method and supersymmetric quantum mechanics2016875, article 7310.1007/s12043-016-1278-72-s2.0-84995422341KhokhaE. M.Abu-ShadyM.Abdel-KarimT. A.Quarkonium masses in the N-dimensional space using the analytical exact iteration method2016228692KuchinS. M.MaksimenkoN. V.Theoretical estimation of spin- average mass spectra of heavy quarkonia and Bc Mesons201313295298Al-JamelA. F.WidyanH.Heavy quarkonium mass spectra in a Coulomb field plus quadratic potential using Nikiforov-Uvarov method20124310.5539/apr.v4n3p94AhmadovA. I.AydinC.UzunO.Bound state solution of the Schrödinger equation at finite temperature20191194, article 01200110.1088/1742-6596/1194/1/0120012-s2.0-85065557616Abu-ShadyM.Abdel-KarimT. A.Ezz-AlarabS. Y.Masses and thermodynamic properties of heavy mesons in the non-relativistic quark model using the Nikiforov–Uvarov method201927110.1186/s42787-019-0014-0MansourH.GamalA.Bound State of Heavy Quarks Using a General Polynomial Potential20182018710.1155/2018/72696572-s2.0-850596814427269657OmugbeE.Non-relativistic eigensolutions of molecular and heavy quarkonia interacting potentials via the Nikiforov Uvarov method202010.1139/cjp-2020-0039Abu-ShadyM.Abdel-KarimT. A.KhokhaE. M.Exact Solution of the N-dimensional Radial Schrödinger Equation via Laplace Transformation Method with the Generalized Cornell Potential201821MutukH.Cornell potential: A Neural Network approach20192019910.1155/2019/31053732-s2.0-850628644733105373FaustovR. N.GalkinV. O.TatarintsevA. V.VshivtsevA. S.Algebraic approach to the spectral problem for the Schrödinger equation with power potentials201215220922610.1142/s0217751x000000942-s2.0-0034688118KangD.WonE.Precise numerical solutions of potential problems using the Crank–Nicolson method200822752970297610.1016/j.jcp.2007.11.0282-s2.0-38549105519SeetharamanM.RaghavanS.VasanS. S.Analytic WKB energy expression for the linear plus coulomb potential198316345546210.1088/0305-4470/16/3/0052-s2.0-0039885934HassanabadiS.RajabiA. A.ZarrinkamarS.Cornell and Kratzer potentials within the semirelativistic treatment20122710, article 125005710.1142/S02177323125005752-s2.0-84859452105RubishV. V.LazurV. Y.ReityO. K.ChalupkaS.SalakM.The WKB method for the Dirac equation with the vector and scalar potentials200454989791910.1023/B:CJOP.0000042643.47049.982-s2.0-4844229622IkhdairS. M.Relativistic bound states of spinless particle by the Cornell potential model in external fields201320131010.1155/2013/4916482-s2.0-84886683544491648Abu-ShadyM.Analytic solution of Dirac equation for extended Cornell potential using the Nikiforov-Uvarov method2015111619DaviesP. C. W.BettsD. S.19942ndLondonChapman and Hall10.1007/978-1-4899-2999-0HazarikaB. J.ChoudhuryD. K.Isgur–Wise function in a QCD-inspired potential model with WKB approximation201788310.1007/s12043-016-1357-92-s2.0-85016577625NikiforovA. F.UvarovV. B.1988BaselBirkhauser10.1007/978-1-4757-1595-8RaniR.BhardwajS. B.ChandF.Mass Spectra of Heavy and Light Mesons Using Asymptotic Iteration Method201870217910.1088/0253-6102/70/2/1792-s2.0-85051227668SergeenkoM. N.Zeroth WKB approximation in quantum mechanics2002https://arxiv.org/abs/quant-ph/0206179OmugbeE.OsafileO. E.OkonI. B.WKB Energy Expression for the Radial Schrödinger Equation with a Generalized Pseudoharmonic Potential202082132010.9734/ajopacs/2020/v8i230112SergeenkoM. N.Quasiclassical analysis of three-dimensional Schrödinger's equation and its solution20111528310010.1142/s0217732300000104ItaB. I.LouisH.AkakuruO. U.Nzeata-IbeN. A.IkeubaA. I.MaguT. O.AmosP. I.EdetC. O.Approximate solution to the Schrodinger equation with Manning-Rosen plus a class of Yukawa potential via the WKBJ approximation method201845323333BenedictI. I.LouisH.NelsonN. I.IkeubaA.OziomaA. U.ThomasM. O.PigwehA. I.MichaelM. O.Approximate l-States solutions to the Schrodinger equation with Manning-Rosen plus Hellmann potential via WKB approximation scheme2018191374510.4038/sljp.v19i1.8050LouisH.ItaB. I.AmosP. I.NzeataI. N.JosephI.IvanO.MaguT. O.WKB solution for inversely quadratic Yukawa potential plus inversely Hellmann potential201724109112SergeenkoM. N.Quantum fluctuations of the angular momentum and energy of the ground state2011131333710.1142/s02177323980000612-s2.0-0008108733SergeenkoM. N.Semiclassical wave equation and exactness of the WKB method19965363798380410.1103/PhysRevA.53.37982-s2.0-0008111061HruskaM.KeungW. Y.SukhatmeU.Accuracy of semiclassical methods for shape-invariant potentials19975553345335010.1103/PhysRevA.55.33452-s2.0-0000416368LangerR. E.On the connection formulas and the solutions of the wave equation193751866967610.1103/PhysRev.51.6692-s2.0-36149009714SukhatmeU.PagnamentaA.Finite eigenfunctions in the WKB approximation1991591094494710.1119/1.16678GriffithsD. J.1995Upper Saddle River, NJ USAPrentice Hall Inc.MerzbacherE.19983rdNew YorkJohn Wiley and Sons Inc.LuF.LvB.WangP.YangH.WKB approximation for a deformed Schrodinger-like equation and its applications to quasinormal modes of black holes and quantum cosmology201893750253210.1016/j.nuclphysb.2018.06.0022-s2.0-85056185539OmugbeE.Non-relativistic Energy Spectrum of the Deng-Fan Oscillator via the WKB Approximation Method202081263610.9734/ajopacs/2020/v8i130107AkcayH.SeverR.Analytical solutions of Schrödinger equation for the diatomic molecular potentials with any angular momentum20125071973198710.1007/s10910-012-0013-y2-s2.0-84864288882IkhdairS. M.SeverR.Exact quantization rule to the Kratzer-type potentials: an application to the diatomic molecules20094541137115210.1007/s10910-008-9438-82-s2.0-63049123888IkhdairS.SeverR.Exact solutions of the radial Schrödinger equation for some physical potentials20075451652710.2478/s11534-007-0022-92-s2.0-36749087717OyewumiK. J.Realization of the spectrum generating algebra for the generalized Kratzer potentials2010https://arxiv.org/abs/1008.1515v1RaniR.BhardwajS. B.ChandF.Bound state solutions to the Schrödinger equation for some diatomic molecules201891410.1007/s12043-018-1622-12-s2.0-85051186367BerkdemirC.PahlavaniM. R.Application of the Nikiforov-Uvarov Method in Quantum Mechanics2012Intech Open10.5772/33510ItaB. I.LouisH.MichaelO. I.NelsonN. I.Solutions of the Schrödinger Equation for the Superposed Screened Coulomb plus Kratzer Fues Potential Using the WKB Approximation Method201854276447452OyewumiK. J.Analytical solutions of the Kratzer-Fues potential in an arbitrary number of dimensions2005181758410.1007/s10702-005-2481-92-s2.0-15244342068MolskiM.KonarskiJ.Modified Kratzer-Fues formula for rotation-vibration energy of diatomic molecules199282692793610.12693/APhysPolA.82.927TezcanC.SeverR.A General Approach for the Exact Solution of the Schrodinger Equation2009https://arxiv.org/abs/0807.2304v2HajigeorgiouP. G.Exact analytical expressions for diatomic rotational and centrifugal distortion constants for a Kratzer–Fues oscillator2006235111111610.1016/j.jms.2005.10.0102-s2.0-32044448129HallR. L.SaadN.Smooth transformations of Kratzer’s potential inNdimensions199810982983298610.1063/1.4768892-s2.0-0009529105BayrakO.BoztosunI.CiftciH.Exact analytical solutions to the Kratzer potential by the asymptotic iteration method2007107354054410.1002/qua.211412-s2.0-33846268439IkhdairS.SeverR.On solutions of the Schrödinger equation for some molecular potentials: wave function ansatz20086369770310.2478/s11534-008-0060-y2-s2.0-48349118208PatrignaniC.Particle Data GroupReview of particle physics20164010, article 10000110.1088/1674-1137/40/10/1000012-s2.0-85049178969OliveK. A.Particle Data GroupReview of Particle Physics2014389, article 09000110.1088/1674-1137/38/9/0900012-s2.0-85014877158GodfreyS.MoatsK.2015https://arxiv.org/abs/1409.0874v3