EXACT SOLUTIONS OF THE KLEIN-GORDON EQUATION WITH HYLLERAAS POTENTIAL Akpan N. Ikot +1, Oladunjoye A. Awoga 1 and Benedict I. Ita 2 1 Theoretical Physics Group, Department of Physics, University of Uyo-Nigeria. 2 Quantum Chemistry Group, Department of Chemistry University of Calabar-Nigeria + e-mail: ndemikot2005@yahoo.com Abstract We present the exact solution of the Klein-Gordon with Hylleraas Potential using the Nikiforov- Uvarov method. We obtain explicitly the bound state energy eigenvalues and the corresponding eigen function are also obtained and expressed in terms of Jacobi Polynomials. KEY WORDS: Hylleraas potential, Bound state, energy eigenvalues PACS: 03.65Ge, 03.65 w, 03.65Ca
1. Introduction The study of the internal structure of any quantum mechanical systems such as nuclei, atom, molecules and solids form an integral part of different sets of orthogonal polynomials [1-5]. In nuclear and high energy physics the study of exact solution of the Klein-Gordon equation (KGE) is of high importance for mixed scalar and vector potentials. However, the problem of exact solution of the KGE for a number of special potential has been of great interests in the recent years. Different techniques have been assumed by different authors to obtain the exact or approximate solution of the KGE with some typical potential. These methods include the Standard Method [1] Supersymmetric Quantum Mechanics [2] Factorization Method [3], Asymptotic Iteration Method (AIM) [4], the Nikiforov-Uvarov (NU) Method [5] and others. The potentials investigated with these techniques include, Manning Rosen potential [6-7], Hulthen potential [8-9], Woods-Saxon potential [10-11], pseudoharmonic[12-13],eckart potential [14-15]harmonic potential [16], Kratzer potential [17-18], Mie-type potential [19], ring shape non-spherical oscillator [20], Hartmann potential [21], Scarf potential [22], Poschl-Teller potential [23], five parameters exponential like potential [24] among others. The relativistic Klein-Gordon equations contain two major objects: the four-vector linear momentum operator and the scalar rest mass. This allows one to introduce two types of potential, (i) the four vector potential V(r) and (ii) the space time scalar potentials(r). Alberto et al [25] have shown that when S(r) = V(r) or S(r) = V(r), the KGE and Dirac equation share the same energy spectrum. The choice of S(r) = V(r) = 2V(r) have been shown by Alhaidari et al [26-30] to yield a non relativistic limit for the KGE. However, analytical solutions of the KGE are only possible for the s-wave with the orbital angular momentum quantum number l = 0
for some potential [31-32]. On the other hand, the case l 0, the KGE can only be solve approximately for some potentials using a suitable approximation techniques [33]. The Hylleraas potential is given by [34-35] V(r) = D 1 ()()() (1) ()()() where D is the dissociation energy, with the function s = e () (2) and the intermediates quantities a, b, c are defined by a =, b = c = (3) with k, k and K being the parameters. The Hylleraas potential is a special case of the multiparameter exponential type potential. This potential contains six parameters and it is very complicated and difficult to evaluate. Since Hylleraas [34] introduce this potential no much work have been reported on the bound state solution with this potential. Recently, we investigated the non-relativistic bound state solution with Hylleraas potential [36]. Motivated by the success in obtaining the bound state solution of the Schrödinger equation with the Hylleraas potentials. We are tempted to solve the solution of the KGE equation with the Hylleraas potential for S(r) = V(r). 2. Review of NU method The NU is based on solving the second order linear differential equation by reducing to a generalized equation of hyper-geometric type. This method has been used to solve the Schrödinger, Klein Gordon and Dirac equation for different kind of potentials [28]. The second-order differential equation of the NU method has the form [5].
ψ (r) + () () ψ (s) + () ψ(s) = 0 (4) () where σ(s) and σ(s), are polynomials at most second degree and τ (s) is a first-degree polynomials. In order to find a particular solution to Eq. (4), we use the ansatz for the wave function as ψ(s) = φ(s)χ(s). (5) It reduces Eq. (4) into an equation of hyper-geometric type σ(s)χ (s) + τ(s)χ (s) + λχ(s) = 0 (6) and the other have function φ(s) is defined as a logarithmic derivative () = () () () (7) The solution of the hyper-geometric type function in Eq. (6) are given by the Rodriques relation χ (s) = () σ (s)ρ(s) (8) where B is a normalization constant and the weight function ρ(s) must satisfy the condition (σρ) = τρ (9) The function π(s) and the parameter λ required for the NU method are defined as follows: π(s) = ± σ + kσ (10) λ = k + π (11) In order to find the value of k in eq. (10), then the expression under the square root must be the square of the polynomials. Therefore, the new eigenvalue becomes where λ = λ = nτ () σ (12) τ(s) = τ + 2π(s) (13)
whose derivative is negative. By comparing Eqs.(11) and (12), we obtain the energy eigen values 3. Exact Bound State Solution of KGE In the relativistic atomic units (ћ = c = 1) for a spindle particle, the three dimensional radial s wave of KGE is [37] [ + (M + S(r)) ]Ψ(r, θ, φ) = [E V(r)] Ψ(r, θ, φ) (14) where is the Laplacian and E is the relativistic energy. The common ansatz for the wave function is Ψ(r, θ, φ) = R(r) r Y (θ, φ),where Y (θ, φ) is the spherical harmonic whose solutions is well-known [38]. Thus the radial part of the KGE becomes, () + (E M ) + V (r) S (r) 2EV(r) + MS(r)R(r) = 0 (15) Now considering the case of equal scalar and vector Hylleraas potential, Eq. (15) becomes () + (E M ) 2D (E + M) 1 ()()() R(r) = 0. (16) ()()() Introducing the new variable, using Eq. (2), then it is straightforward to show that, Eq. (16) reduces to the form, + ()() ()()() + () () () [ ε s + β s + γ ]R(r) = 0 (17) where we have used the following physical parameters: ε = () () (18) β = ()()()() () (19) γ = ()()() () (20) where E = E M and V = 2D (E + M). Equation (17) is now amenable to NU method and in order to find the solution of this equation, it is necessary to compare Eq. (4) with Eq. (17). By this comparison, we obtain the following polynomials:
τ (s) = 2(1 + b)(s + c), σ(s) = 2(s + a)(s + c)(1 + b), σ(s) = ε s + β s + γ (21) where β = ε β (22) γ = ε γ (23) and β and γ are defined as respectively. β = ()() () (24) γ = ()() () (25) Substituting these polynomials into Eq. (10), we obtain the π(s) function as π(s) = α (s + a) ± (α + ε + 2α k)s + (ε + ξ + α k)s + ε + ξ + α k, (26) where ξ = 2aα β (27) α = (1 + b), α = 2α (a + c) (28) ξ = α a γ, α = 2α ac. (29) According to the NU method the expression in the square root must be square of the polynomial. Thus, we determine the complex k value as = (Λ + Λ ε ) ± δ(ε + A) A B, (30) where Λ = 2α ξ 4α α 8α ξ (31) Λ = (2α 4α 8α ) (32)
A = (33) B = (34) δ = Λ + 12Λ (35) and Λ and Λ are defined as Λ = α 8α α (36) Λ = ξ 4α ξ (37) respectively. The polynomial of π(s) for each k value is π(s) = α (s + a) ± δ(ε + A) + A Bs + δ(ε + A) A B for k = (Λ + Λ ε ) + δ(ε + A) A B δ(ε + A) + A Bs δ(ε + A) A B for k = (Λ + Λ ε ) δ(ε + A) A B The correct value of π(s) is chosen such that the function τ(s) has a negative derivative. Thus we select the k and π(s) values as (38) k = (Λ + Λ ε ) δ(ε + A) A B (39) π(s) = α (s + a) δ(ε + A) + A Bs + δ(ε + A) + A B (40) which gives τ(s) = 2α [2s + (a + c)] 2 δ(ε + A) + A Bs δ(ε + A) A B (41) and
τ (s) = 2 δ(ε + A) + A B 4α < 0. (42) Using λ = k + π and its counterpart definition, λ = nτ [(n(n 1)/2)]σ given by Eqs. (12 ) and (13 ), the following expressions for λ and λ are obtained λ = (Λ + Λ ε ) U V (43) λ = 2n U + V 2α n(n + 3) (44) where U = δ (ε + A) and V = (A B). Now using Eqs. (43) and (44), we obtain the complicated energy eigenvalues of the Klein- Gordon equation with Hylleraas potential as E = () Λ () + 1 + A B () ± () () X () () (() (45) Equation(45) is rather complicated transcendental equation. However, the energy eigenvalues can be found by setting the Hyllerass parameters in Appendix A into equation (45). Let us turn attention to find the radial wave function for potential. Using τ(s), π(s) and σ(s) in Eqs (7) and (9), we find φ(s) = (s + a) (s + c) (46a) where D = () (),F = () (46b) and
μ = δ(ε + A) + A B, ν = δ(ε + A) A B, (46c) Then from Eq. (8), we obtain ρ(s) = (s + a) (s + c) (47) χ (s) = B (s + a) (s + c) (s + a) (s + c) (48) and the total wave function can be written in the form ψ(r, θ, φ) = N [a + exp(2(1 + K)ωr)] [c + exp(2(1 + K)ωr)] P (,) (r)y (θ, φ)(49) where N is normalization constant and P (,) (r) is the Jacobi polynomials. The wave function obeys the normalization condition, 6. Conclusion r ψ(r, θ, φ) dτ, where dτ is the volume element. In this work, we used the NU method for Hylleraas potential to obtain the exact solution of the Klein-Gordon equation. We obtain the exact energy eigen values equation and the unnormalized radial wave function expressed in terms of the hypergeometric function. Finally, this result will open ways for more research work on the bound state solution of the Klein-Gordon equation with the six parameter Hylleraas potential and this will leads many authors to study the bound state energy spectra of some diatomic molecules with Hylleraas potential [34-35]. Appendix A: Values of formulas in terms of the Hylleraas potential The following are the useful formulas and relation satisfied by the energy eigen values of the Hylleraas potential of Eq. (45). The α, α, α are given in terms of Hylleraas parameters as α = (1 + b) (A. 1) α = 2(1 + b)(a + c) (A. 2) α = 2ac(1 + b) (A. 3)
The expression β is defined as β = ()() () (A. 4) The parameters Λ, Λ, Λ, Λ are Λ = 4(1 + b) (a 14ac + c ) (A. 5) Λ = 4a(1 + b) (2a c + 1) ()()()() () (A. 6) Λ = 4(1 + b)(a + c 2ac 2) (A. 7) Λ = ()() () b(1 + b) + ()() () (A. 8) Other essentials parameters expressed in terms of the Hylleraas parameters are defined as follows: δ = 64(1 + b) [a(1 c) a(8c + 1) + c (1 a) + a + 4] (A. 9) ξ = 2a(1 + b) ()() () (A. 10) ξ = a (1 + b) ()()() () (A. 11) γ = ()()() () (A. 12) Finally, the constants A and B are defined as A = () ()() () () [ ] B = () () ()() () ()() () () [ ] (A. 14) (A. 13)
REFERENCES 1. Y. Xu, S. He and C. S. Jia, Phys. Scr. 81 (2010) 045001. 2. M. R. Setare and Z. Nazari, Acta Polonica B. 40(10) (2009) 2809. 3. G. F. Wei, X. Y. Liu and W. L. Chen., Int. J. Theor. Phys. 48 (2009) 1649. 4. T. Baraket, Int. J. Mod. Phys. A 21 (2006) 4127. 5. A. F. Nikiforov and V. B. Uvarov, Special Functions of Mathematical Physics, (Birkhausen, Basel 1988). 6. F. Taskin, Int. J. Theor. Phys. 48 (2009) 1142. 7. M. F. Manning and N. Rosen, Phys. Rev. 43 (1933). 8. N. Saad. Phys. Scr. 76 (2007) 623. 9. A. N. Ikot, L. E. Akpabio and E. J. Uwah, EJTP, 8(25) (2011) 225. 10. C. S. Jia, P. Gao and X. L. Peny, J. Phys. A: Math. Gen. 39 (2006) 7737. 11. J. Y. Guo and Z. Q. Sheng, Phys. Lett. A338, (2005) 90. 12. G. Chen, Z. D. Chen and Z. M. Lou, Chin. Phys. 13 (2004) 279. 13. K. J. Oyewumi, F. O. Akinpelu and A. D. Agboola, Int. J. Theor. Phys. 47 (2008) 1039. 14. G. F. Wei, S. H. Dong and V. B Bezerra, Int. J. Mod. Phys. A 24 (2009) 161. 15. A. D. Alhaidari, Found. Phys. 40 (2010) 1088. 16. R. K. Sun and Z. Q. Ma, J. Phys. A: Math Gen. 19 (1986) 1739. 17. W. C. Qiang, Chin. Phys. 13 (2004) 575. 18. K.J.Oyewumi,Found.Phys.Lett.18(1)(2005)75 19. S.Ikhdair and R.Sever,Cent.Euro.J.Phys.6(2008)221 20. X. Q. Zhao, C. S. Jia and Q. B. Yang, Phys. Lett. A 337 (2005) 189. 21. H.Hartman and D.Schuch,Int.J.Quant.Chem.18(1980)125.
22. R.De,R.Dutt and U.Sukhatme,J.Maths.Gen.25(1992)L843. 23. G. Poschl and E. Teller, Z. Phys. 83 (1933) 143. 24. A. N. Ikot, L. E. Akpabio and J. A. Obu, Jour. Vect. Relat, 6 (2011)1. 25. P. Alberto, A. S. de Castro and M. Malheiro, Phys. Rev. C 75 (2007) 047303. 26. A. D. Alhaidari, H. Bahlouli and A. Al-Hasan, Phys. Lett. A 349 (2006) 87. 27. A. D. Alhaidari, Phys.Lett. B699 (2011) 309 28. A. D. Alhaidari, J. Phys. A. Math. Gen. 34 (2001) 9827. 29. M.R.Setare and Z.Nazari,Acta Physica Polonica B 41(11)(2010)2459. 30. A. D. Alhaidari, Phys. Rev. Lett. 87 (2001) 210405-1 31. A. Arda, R. Server and C. Tezcan, Phys. Scr. 79 (2009) 5006. 32. A. Arda, and R. Server, Int. J. Theor. Phys. 48 (2009) 945. 33. R. L. Greene and C. Aldrich, Phys. Rev. A14 (1976) 2363. 34. E. A. Hylleraas, J. Chem. Phys. 3 (1938) 595. 35. Y. P. Varshni, Rev. Mod. Phys. 29 (4) (1957) 664. 36. A.N.Ikot etal,bound state solutions of Schrodinger equation with Hyllerass potential Using Nikiforov-Uvarov Method(Submitted for publication) 37. W. Greiner, Relativistic Quantum Mechanics, (Springer, Berlin, 2000). 38. S. H. Dong, Factorization Methods in Quantum Mechanics, (Springer, Dordrecht, 2007).