arxiv: v1 [quant-ph] 9 Oct 2008

Bound states of the Klein-Gordon equation for vector and scalar general Hulthén-type potentials in D-dimension Sameer M. Ikhdair 1, 1 Department of Physics, Near East University, Nicosia, North Cyprus, Turkey Dated: May 30, 018) arxiv:0810.1590v1 [quant-ph] 9 Oct 008 Abstract We solve the Klein-Gordon equation in any D-dimension for the scalar and vector general Hulthén-type potentials with any l by using an approximation scheme for the centrifugal potential. Nikiforov-Uvarov method is used in the calculations. We obtain the bound state energy eigenvalues and the corresponding eigenfunctions of spin-zero particles in terms of Jacobi polynomials. The eigenfunctions are physical and the energy eigenvalues are in good agreement with those results obtained by other methods for D = 1 and 3 dimensions. Our results are valid for q = 1 value when l 0 and for any q value when l = 0 and D = 1 or 3. The s-wave l = 0) binding energies for a particle of rest mass m 0 = 1 are calculated for the three lower-lying states n = 0,1,) using pure vector and pure scalar potentials. Keywords: Bound states, Klein-Gordon equation, Hulthén potential, Nikiforov-Uvarov method, Approximation schemes. PACS numbers: 03.65.-w; 03.65.Fd; 03.65.Ge E-mail: sikhdair@neu.edu.tr 1

I. INTRODUCTION In nuclear and high energy physics [1,], it is necessary to obtain the exact bound state energy spectrum of non-relativistic or relativistic wave equations for various potentials. In relativistic quantum mechanics, the Klein-Gordon and Dirac wave equations are most frequently used to understand the dynamics of a particle in large class of potentials. Therefore, much works have been done to solve these relativistic wave equations for various radial and angular potentials. Unfortunately, some quantum mechanical equations with these potentialscanbesolvedanalyticallyonlyfors-waveequationi.e., l = 0states). Someotherworks have been done to study the l 0 bound state solutions of the Klein-Gordon equation with pure scalar, pure vector and mixed vector and scalar potentials in various dimensional space using different methods. The scalar coupling constant is almost taken to be equal to the vector coupling constant i.e., S 0 = V 0 ) [-14]. For this case, the Klein-Gordon equation reduces to a Schrödinger-like equation and thereby the bound state solutions are easily obtained by using the well-known methods developed in nonrelativistic quantum mechanics []. On the other hand, it has been shown that when S 0 V 0, there exist real bound state solutions i.e., the radial wave function must satisfy the boundary condition that it becomes zero when r and also finite at r = 0). However, there are very few exact solvable cases [15-19]. The bound state solutions for the last case is obtained for the s-wave Klein-Gordon equation with the exponential [15] and some other classes of potentials [16]. Chen et al. [19] used the exponential function transformation approach along with an approximation for the centrifugal potential to find bound-state solutions of the Klein-Gordon equation with the vector and scalar Hulthén potentials for different l. The D-dimensional Schrödinger and Klein-Gordon wave equations have been solved for various types of angular and radial dependent potentials using the Nikiforov-Uvarov N-U) method [0-5]. The exact bound state solutions for the Schrödinger, Klein-Gordon and Salpeter wave equations, in one-dimension 1D), with the generalized Woods-Saxon potential are also obtained [6-9]. These solutions of the D-dimensional Klein-Gordon equation describing a spin-zero particle are obtained for the case V 0 = S 0 ring-shaped Kratzer-type and pseudoharmonic potentials [1,]. Furthermore, analytic solution of radial Schrödinger equation has been given for the l 0 general Hulthén potential within an approximation to the centrifugal barrier term [5,9]. Simsek and Egrifes [30] have investigated the reality of

exact bound states of the 1D Klein-Gordon equation with complex and/or PT-symmetric non-hermitian exponential-type like general Hulthén potential. Berkdemir [31] has also applied the method to solve the Klein-Gordon equation of a spin-zero particle for Vr) = Sr) Kratzer-type potentials. Qiang et al. [3] have given an approximate analytic solution of the l 0 Klein-Gordon equation for the general Hulthén-type potential within the approximation given in [9]. Saad [33] has presented an approximate solution of the l 0 bound states of the Klein-Gordon equation in D-dimensions with the general Hulthén-type potentials following the model used in [9,3]. Chen et al. [34] have employed two semiclassical methods to determine the bound state energy spectrum of Klein-Gordon equation when S 0 V 0. Dong et al [35] have also obtained the bound state solutions of the Schrödinger equation for an exponential-type potential following Refs. [9,3]. Very recently, we have obtained the bound-state solutions for the D-dimensional Schrödinger equation with the Manning- Rosen potential [36] which can be reduced to the Hulthén potential and the exponential-type hyperbolic potential [37] using a novel approximation to the centrifugal term [9,36,37]. The purpose of the present paper is to study the relativistic characteristics of the scalar and vector general Hulthén-type potentials written, respectively, as follows [3,33]: Vr) = V 0e r/r 0 1 qe r/r 0, Sr) = S 0e r/r0 1 qe r/r 0, r 0 = α 1, q 0 1) where V 0 and S 0 represent the coupling constants of the vector and scalar general Hulthéntype potentials, respectively, α is the screening range parameter, r 0 represents the spatial range and q is the deformation parameter [6-8]. This potential has been widely used in a number of areas of physics. In atomic physics, V 0 = Ze /r 0 where Z is the atomic number. Equation 1) behaves like a Coulomb potential V C r) = Ze /r when r r 0, but decreases exponentially in the asymptotic region when r r 0 ), so its capacity for bound state is smaller than the Coulomb potential. The Klein-Gordon equation has been solved analytically in [30] only for the s-wave. This equation cannot be solved analytically for l 0 because of the centrifugal term, 1/r. Therefore, for l 0, we must use an approximation for the centrifugal term similar to other authors[9,38-43] to find approximate analytical solutions of the D-dimensional Klein-Gordon equation with vector and scalar general Hulthén-type potentials. However, the resulting analytic solutions are only valid for q = 1 in the l 0 case and for any q value in the l = 0 and D = 1,3. This paper is organized as follows: In Section II, we will derive l 0 state solutions 3

within the approximation given in [9] for the D-dimensional Klein-Gordon equation describing spin-zero particle with the vector and scalar general Hulthén-type potentials. We further separate the wave equation into radial and angular parts. Section III is devoted to a brief description of the N-U method. In Section IV, we present general solutions to the radial and angular equations in D-dimensions. We also give discussions of the solution in various dimensions and the constraints for obtaining the real bound state solutions. Further, numerical values for the three lower-lying states when pure vector and pure scalar potentials are used with coupling constants V 0 = 0.5 and S 0 = 0.5 and unity mass particle. Finally, the relevant conclusions are given in Section V. II. THE KLEIN-GORDON EQUATION WITH SCALAR AND VECTOR POTEN- TIALS The time-independent Klein-Gordon equation in any arbitrary D-dimension) with scalar Sr) and vector Vr) potentials, r = r,describing spin-zero particle of rest mass, m 0, can be written as in the relativistic natural units h = c = 1) [,1,] { D +[E nl Vr)] [m 0 +Sr)] } ψ l D 1=l) l 1 l D x) = 0, D = where E nl denotes the relativistic energy and D D x j=1 j, ) denotes the D-dimensional Laplacian. Usually, it is required that S 0 V 0 and m 0 > E nl for the existence of real bound states [15-17]. Furthermore, the x is a D-dimensional position vector, with the unit vector along x is usually denoted by x = x/r, is the hyperspherical Cartesian components x 1,x,,x D given as follows cf. Refs. [0-5] and the references therein): x 1 = rcosθ 1 sinθ sinθ D 1, x = rsinθ 1 sinθ sinθ D 1, x 3 = rcosθ sinθ 3 sinθ D 1,. x j = rcosθ j 1 sinθ j sinθ D 1, 3 j D 1,. 4

x D 1 = rcosθ D sinθ D 1, D x D = rcosθ D 1, x j = r,d. 3) j=1 We have x 1 = rcosθ 1, x = rsinθ 1, θ 1 = ϕ, for D = and x 1 = rcosθ 1, x = rsinθ 1 sinθ, x 3 = rcosθ, θ = θ, for D = 3. The volume element of the configuration space is given by D D 1 dx j = r D 1 d 3 rdω, dω = sinθ j ) j 1 dθ j, 4) j=1 where r [0, ), θ 1 [0,π] and θ j [0,π], j [,D 1]. The wave function ψ l D 1=l) l 1 l D x) with a given angular momentum l can be decomposed as a product of a radial wave function R l r) and the normalized hyper-spherical harmonics Y l) l 1 l D θ 1,θ,,θ D 1 ) as j=1 ψ l D 1=l) l 1 l D x) = R l r)y l) l 1 l D θ 1,θ,,θ D 1 ), 5) where Y l) l 1 l D θ 1,θ,,θ D 1 ) = Hθ 1 )Hθ,,θ D )Hθ D 1 ), 6) which is the simultaneous eigenfunction of L j : L jy l) l 1 l D θ 1,θ,,θ D 1 ) = l j l j +j 1)Y l) l 1 l D θ 1,θ,,θ D 1 ), l = 0,1,,l k = 0,1,,l k+1, j [1,D 1], k [,D ], l 1 = l, l +1,,l 1,l, L D 1 Y l) l 1 l D θ 1,θ,,θ D 1 ) = ll+d )Y l) l 1 l D θ 1,θ,,θ D 1 ). 7) The angular momentum operators L j are defined as k+1 L k = a<b= L 1 =, θ1 L ab = 1 sin k 1 θ sin k 1 k θ k θ k θ k L ab = i ) [ x a x b x b x a + L k 1 sin, k D 1, θ k ]. 8) Employing the method of separation of variables and substituting Eqs. 5)-7) into Eq. ), we obtain the following Schrödinger-like equation: R l r)+ D 1) R r lr)+ { ll+d ) r +[E nl Vr)] [m 0 +Sr)] 5 } R l r) = 0, 9)

where ll + D )/r is known as the centrifugal term. Furthermore, using R l r) = r D 1)/ gr), it is straightforward to find out the radial wave equation: { g r)+ [E nl Vr)] [m 0 +Sr)] D +l 1)D+l 3) 4r } gr) = 0. 10) Obviously, Eq. 10) can be analytically solved [17,7,30] only for l = 0 s-wave) case. For l 0 case, we have to use an approximation for the centrifugal term similar to the non-relativistic cases which is valid for q = 1 value [9,3,33,38-43]: 1 r α e αr 1 qe αr ). 11) We follow the model used by Qiang et al. [3] and rewrite Eq. 10) for the scalar and vector general Hulthén-type potentials 1) as, g r)+ [E nl m0 + m ] 0S 0 +E nl V 0 )e αr S 0 V0 )e αr D +l 1)D+l 3) 1 qe αr 1 qe αr ) 4r gr) = 0. 1) On the other hand, the other angular equations, obtained from the separation procedures, are [1,]: [ ) 1 d sin j 1 sin j 1 d θ j +l j l j +j 1) l ] j 1l j 1 +j ) θ j dθ j dθ j sin Hθ,,θ D ) = 0, θ j [ 1 sin D θ D 1 d dθ D 1 sin D θ D 1 j [,D ], 13) ) ] d +ll+d ) L D sin Hθ D 1 ) = 0, 14) θ D 1 dθ D 1 d Hθ 1 ) +l dθ 1Hθ 1 1 ) = 0. 15) It is well-known that the solution of Eq. 15) is H l1 θ 1 ) = 1 π exp±il 1 θ 1 ), l 1 = 0,1,,. 16) Hence, the above Eqs. 1)-14), have to be solved by using N-U method [0-31,36] which is reviewed briefly in the following Section. 6

III. NIKIFOROV-UVAROV METHOD The N-U method is breifly outlined here and the details can be found in [0-31,36,37,44]. N-U method is proposed to solve the second-order differential equation of the hypergeometric type: ψ ns)+ τs) σs) ψ ns)+ σs) σ s) ψ ns) = 0, 17) where σs) and σs) are polynomials, at most of second-degree, and τs) is a first-degree polynomial. Using a wave function, ψ n s), of the simple form reduces Eq. 17) into an equation of a hypergeometric type ψ n s) = φ n s)y n s), 18) σs)y ns)+τs)y ns)+λy n s) = 0, 19) where σs) = πs) φs) φ s), 0) and λ is a parameter defined as τs) = τs)+πs), τ s) < 0, 1) λ = λ n = nτ s) nn 1) σ s), n = 0,1,,. ) The polynomial τs) with the parameter s and prime factors show the differentials at first degreebenegative. Theotherparty n s)ofthewavefunctioneq. 18)isthehypergeometrictype function whose polynomial solutions are given by Rodrigues relation y n s) = B n d n ρs) ds n [σn s)ρs)], 3) where B n is the normalization constant and the weight function ρs) can be found by [44] ω s) τs) ωs) = 0, ωs) = σs)ρs). 4) σs) The function πs) and the parameter λ are defined as σ ) πs) = σ s) τs) s) τs) ± σs)+kσs), 5) 7

λ = k +π s), 6) where πs) has to be a polynomial of degree at most one. To obtain k, the expression under the square root sign in Eq. 5) can be arranged to be the square of a polynomial of first degree [44]. This is possible only if its discriminant is zero. Finally, the energy eigenvalues are obtained by comparing Eqs. ) and 6). IV. SOLUTIONS OF THE RADIAL AND ANGLE-DEPENDENT EQUATIONS A. The D-dimensional angular equations At the beginning, we rewrite Eqs. 13) and 14) representing the angular wave equations in the following simple forms [1,] : d Hθ j ) dθ j +j 1)ctgθ j dhθ j ) dθ j + Λ j Λ j 1 sin θ j ) Hθ j ) = 0, j [,D ],D > 3, 7) [ d Hθ D 1 ) dhθ D 1 ) +D )ctgθ dθd 1 D 1 + ll +D ) Λ ] D dθ D 1 sin Hθ D 1 ) = 0, 8) θ D 1 where Λ p = l p l p +p 1), p = j 1, j with l p = 0,1,, are the angular quantum numbers in D-dimensions. Employing s = cosθ j, we transform Eq. 7) to the associated-legendre equation d Hs) ds js dhs) + Λ j Λ j 1 Λ j s Hs) = 0, j [,D ],D > 3. 9) 1 s ds 1 s ) By comparing Eqs. 9) and 17), the corresponding polynomials are obtained τs) = js, σs) = 1 s, σs) = Λ j s +Λ j Λ j 1. 30) Inserting the above expressions into Eq. 5) and taking σ s) = s, one obtains the following function: j ) j ) πs) = s± [ +Λ j k] s +k Λ j +Λ j 1. 31) Following the method, the polynomial πs) is found to have the following four possible values: πs) = j + Λ j 1 )s for k 1 = Λ j Λ j 1, j Λ j 1 )s for k 1 = Λ j Λ j 1, j ) s+ Λ j 1 for k = Λ j + ) j, j ) s Λ j 1 for k = Λ j + ) j, 8 3)

where Λ p = l p +p 1)/,with p = j 1, j and j [,D ], D > 3. Imposing the condition τ s) < 0, where τ = τ l j 1 ), we may select the following physical solutions: k 1 = Λ j Λ j 1 and πs) = l j 1 s, j [,D ],D > 3, 33) which yields τs) = 1+ Λ j 1 )s. 34) Using Eqs. ) and 6), give λ nj = n j 1+ Λ j 1 )+n j n j 1), 35) λ = Λ j Λ j 1 l j 1, 36) and setting λ = λ nj, we obtain n j = l j l j 1. 37) Substituting Eqs. 30), 33) and 34) into Eqs. 0)-1) and 3)-4), give φs) = 1 s ) l j 1 /, ρs) = 1 s ) l j 1 +j )/. 38) Further, substituting the weight function ρs) given in Eq. 38) into the Rodrigues relation 3) gives y nj s) = A nj 1 s ) e Λ j 1 d n j ds n j 1 s ) n j + e Λ j 1, 39) where A nj is the normalization factor. Finally, the angular wave functions are H nj θ j ) = N nj sinθ j ) l j n j P l j n j +j )/,l j n j +j )/) n j cosθ j ), 40) where the quantum numbers n j are defined in Eq. 37) and the normalization factor is l j +j 1)n j! N nj =, j [,D ], D > 3. 41) Γl j +l j 1 +j ) Likewise, in solving Eq. 8), we introduce a new variable s = cosθ D 1. Thus, we can also rearrange it as the universal associated-legendre differential equation d Hs) ds D 1)s dhs) + ν1 s ) Λ D Hs) = 0, 4) 1 s ds 1 s ) where ν = ll+d ). 43) 9

Equation 4) has been recently solved in D and 3D by the N-U method in [45,46]. By comparing Eqs. 4) and 17), the corresponding polynomials are obtained τs) = D 1)s, σs) = 1 s, σs) = νs +ν Λ D. 44) Inserting the above expressions into Eq. 5) and taking σ s) = s, one obtains the following function: D 3) D ) πs) = s± [ 3 +ν k] s +k ν +Λ D, 45) which gives the possible solutions: D 3 +m ) s for k 1 = ν Λ D, D 3 m ) s for k 1 = ν Λ D, πs) = D 3) s+m for k = ν + ) D 3 46), D 3) s m for k = ν + ) D 3, where m = l D + D 3. Imposing the condition τ s) < 0, where τ = τ l D ), one may select: k 1 = ν Λ D and πs) = l D s, 47) which yields τs) = 1+m )s. 48) The following expressions for λ are obtained, respectively, λ = λ nd 1 = n D 1 1+m )+n D 1 n D 1 1), 49) λ = ll+d ) l D l D +D ), 50) and from λ = λ nd 1, we obtain the angular momentum quantum number: l = l D +n D 1. 51) where n D 1 = n and l 1 = m in 3D space [46]. Using Eqs. 0)-1) and 3)-4), we obtain φs) = 1 s ) l D / ), ρs) = 1 s m. 5) The Rodrigues relation 3) gives the following wave functions: y nd 1 s) = B nd 1 1 s ) m d n D 1 ds n D 1 1 s ) n D 1 +m, 53) 10

where B nd 1 is the normaliation factor. Finally, the angular wave functions are where the normalization factor is H nd 1 θ D 1 ) = N nd 1 sinθ D 1 ) l D P m,m ) n D 1 cosθ D 1 ), 54) N nd 1 = with m is defined after Eq. 46). n D 1 +m +1)n D 1!, 55) Γn D 1 +m ) B. Bound States of the D-dimensional Radial Equation The D-dimensional Klein-Gordon radial energy eigenvalue equation with the vector and scalar general Hulthén-type potentials can be rewritten as g r)+ [E nl m0 + m 0S 0 +E nl V 0 )e αr 1 qe αr 1 D +l 1)D 4 +l 3)α e αr +S0 V 1 qe αr ) 0 )e αr ] gr), g0) = 0, 56) On account of the wave function gr) satisfying the standard bound state condition, which is when r, the wave function gr) 0. Equation 56) can be further transformed by using a new variable s = qe αr r [0, ), s [q,0]), g s)+ 1 s) 1 [ s1 s) g s)+ ε [s1 s)] nl +β 1 γ +ε nl )s β 1 +β +ε nl )s] gs) = 0, where m ε nl = 0 Enl α, β 1 = m 0S 0 +E nl V 0 ) α q, β = S 0 V 0, γ = α q 57) D +l 1)D+l 3). 4q 58) For bound states, E nl m 0, ε nl 0 [17,47,48]. Comparing Eqs. 57) and 17), the corresponding polynomials are obtained: τs) = 1 s,σs) = s1 s), σr) = ε nl +β 1 γ +ε nl )s β 1 +β +ε nl )s. 59) The substitution of Eq. 59) into Eq. 5) and taking σ s) = 1 s, give the polynomial: πs) = s ± 1 [1+4β 1 +β +ε nl k)]s +4k β 1 +γ ε nl )s+4ε nl. 60) 11

If the expression under the square root in Eq. 60) is set equal to zero and solved for k, we obtain: k = β 1 γ ±ε nl 1+4β +γ). 61) In view of that, we arrive at the following four possible functions of πs) : [ ] s +ε nl ε nl 1 1+4β +γ) s for k 1 = β 1 γ +ε nl 1+4β +γ), [ ] s πs) = ε nl + ε nl 1 1+4β +γ) s for k 1 = β 1 γ +ε nl 1+4β +γ), [ ] s +ε nl ε nl + 1 1+4β +γ) s for k = β 1 γ ε nl 1+4β +γ), [ ] s ε nl + ε nl + 1 1+4β +γ) s for k = β 1 γ ε nl 1+4β +γ). Thecorrectvalueofπs)ischosensuchthatthefunctionτs)inEq. 1)must havenegative derivative [44]. So we can select the physical values to be k = β 1 γ ε nl 1+4β +γ) and πs) = s +ε nl 6) [ ε nl + 1 ] 1+4β +γ) s, 63) which yield [ τs) = 1+ε nl 1+ε nl + 1 ] 1+4β +γ) s, [ τ s) = 1+ε nl + 1 ] 1+4β +γ) < 0. 64) Using Eqs. ) and 6), the following expressions for λ are obtained, respectively, λ = λ n = n + [ 1+ε nl + ] 1+4β +γ) n, n = 0,1,, ), 65) λ = β 1 γ 1 [ 1+ε nl) 1+ ] 1+4β +γ), 66) where n is the radial quantum number. Solving the last two equations and using β = δ ± δ ± γ, give ε D) nl = β 1 γ n ) n+1)δ ± n+δ ± ) δ ± = 1 = qm 0S 0 +E nl V 0 )+S0 V0 q α n+δ ± ) 1± a ), a = q = 4q[β 1 n n+1)δ ± ] D+l 1)D +l 3) 8qn+δ ± ) n+δ ±, n = 0,1,, ), q + 4S 0 V 0 ) α +qd +l 1)D+l 3), l = 0,1,, ), 67) 1

with δ = δ + for q > 0 and δ = δ for q < 0. Furthermore, from setting ε D) nl = ε nl and using Eqs. 58) and 67), we obtain the energy equation m 0 E nl = qm 0S 0 +E nl V 0 )+S 0 V 0 q αn+δ) αn+δ), 68) or equivalently the explicit expression of the energy eigenvalues E D) nl = η ) nlv 0 m 0 ηnl ±κ nl,16q m q 4V0 +κ 0 nl 4q η nl 4V 0 +κnl), η nl = 4V 0 S 0 )+κ nl 8qm 0S 0, 4V0 +κ nl κ nl = qαn+1)± q α +4S0 V 0 )+qα D +l 1)D+l 3), n = 0,1,,3,, l = 0,1,, ), 69) and q α +qα D+l ) +4S 0 qα +4V 0, 70) is a constraint over the potential parameters. In the case of pure vector potential S 0 = 0,V 0 0), the energy equation 69) reduces to E D) nl = V 0 q ±κ m 0 nl 4V0 +κ nl κ nl = qαn+1)± ) 1, 16q m 0 4q 4V 0 +κ nl, q α +qα D +l 1)D+l 3) 4V 0, D 1, n = 0,1,,3,, l = 0,1,, ), 71) with the constraint q α +qα D+l ) qα +4V 0, 7) forrealboundstates. Fromtheaboveresult, itisnotdifficulttoconcludethatthetwoenergy solutions are valid for the particle and the second one corresponds to the anti-particle energy. The restriction that gives the critical coupling value leads to the result n 1 ) qα 4q qα m 0 V0 [q +D +l 1)D+l 3)] V 0 1 4, 73) i.e., there are only finitly many eigenvalues. In order that at least one level might exist, it is necessary that the inequality: qα+ qα [q +D +l 1)D+l 3)] 4V0 4q m 0 V 0, 74) 13

is fulfilled. As can be seen from Eq. 73), there are at most only two lower-lying states n = 0,1)for theklein-gordonparticleof mass unity when theparameter α = 1 andq = ±1 for any arbitrary value of V 0. Therefore, we have the inequality: D ) n ± +l 4 V0 V0 1. 75) For instance, if one selects V 0 = D+l )/, then the real bound states must be restricted 16 D ) by n +l ) 1 /. Having solved the D-dimensional Klein-Gordon equation for scalar and vector general Hulthén-type potentials, we should make some remarks. i) For s-wave l = 0), the exact energy eigenvalues of the 1D Klein-Gordon equation becomes E n = η nv 0 q ± κ n 4q4V0 +κ n ) [κ n +4V0 S0)] [ S 0 +4qm 0 ) ] κ n 4V0, 76) where κ n = qαn+1)+ q α +4S 0 V 0 ), η n = 4V 0 S 0 )+κ n 8qm 0S 0. 77) 4V0 +κ n In order that at least one level might exist, it is necessary that the inequality 16q m 0 η n ) 4V 0 +κ n, q α +4S0 4V 0, 78) is fulfilled. In the case of pure vector potential S 0 = 0,V 0 0), the energy spectrum E n = V [ ] 0 q ± qαn+1)+ q α 4V0 m 0 with [ 4V0 + qαn+1)+ q α 4V 0 ] 1 16q, 79) [ 16q m 0 4V0 + qαn+1)+ q α 4V0], qα V 0. 80) We notice that the result given in Eq. 79) is identical to Eq. 31) of Ref. [30]. There are only two lower-lying states n = 0,1) for the Klein-Gordon particle of a rest mass m 0 = 1.0 and α = 1.0. For example, one may calculate the ground state energy for the coupling strength V 0 = qα/ as E 0 = V 0 q [ 1± 14 m 0 q V 0 1 ]. 81)

Further, in the case of pure scalar potential V 0 = 0,S 0 0), the energy spectrum [ 1 E n = ± [ ] 4q qαn+1)+ qαn+1)+ q α +4S0 S 0 +4m 0 q) q α +4S 0] 4S 0 [ qαn+1)+ q α +4S0]. 8) In the above equation, all bound states appear in pairs, with energies ±E n. Since the Klein- Gordon equation is independent of the sign of E n for scalar potentials, the wavefunctions become the same for both energy values. We notice that Eq. 8) is identical to Eq. 4) of Ref. [47] obtained by N-U method and to Eq. 0) of Ref. [48] obtained by supersymmetric method. If the range parameter α is chosen to be α = 1/λ c, where λ c = h/m 0 c = 1/m 0 denotes the Compton wavelength of the Klein-Gordon particle. It can be seen easily that while S 0 0 in ground state n = 0), all energy eigenvalues tend to the value E 0 0.866 m 0. ii) For D = 3, the mixed scalar and vector Hulthén potentials q = 1,l 0), the energy eigenvalues are E nl = η [ nlv 0 ± αn+1)+ ] m α 0 +4B 4V0 +[ αn+1)+ α +4B ] where ηnl 4 ), 83) η nl = 4V 0 S 0 )+[ αn+1)+ α +4B ] 8m0 S 0 4V0 + [ αn+1)+ α +4B ], B = S0 V 0 +α ll +1). 84) In order that at least one level might exist, it is necessary that the inequality { [ 16m 0 ηnl 4V0 + αn+1)+ ] } α +4B, α +4B 0, 85) is fulfilled. In the case of pure vector S 0 = 0,V 0 0),the energy eigenvalues are E nl = V [ ] 0 ± αn+1)+ α l+1) 4V0 m 0 [ 4V0 + αn+1)+ ] 1 α l+1) 4V0 16, 86) with the following constraint over the potential parameters: 16m 0 4V 0 + [ ] αn+1)+ α l+1) 4V0, 15

l +1)α V 0. 87) A preliminary analysis about the possibility of obtaining real bound states for the Klein- Gordon particle of mass unity is done for various dimensions as follows: when D = 3, we have only two lower-lying real bound states n = 0,1) for the angular quantum number l = 0, one lower-lying state n = 0) for l = 1 and no states for l. When D = 4, we have two lower-lying states for l = 0 and no states for l 1. Further, when D = 5, we have one lower-lying state for l = 0 and no states for l 1. Finally, when D 6, there are no real bound states for l 0, however, only scattering states will be possible. iii) When D = 3 and l = 0, the centrifugal term D+l 1)D+l 3) 4r = 0, and the approximation term D+l 1)D+l 3)α e αr = 0, too. Thus, letting l = 0 and D = 3 in Eq. 41 qe αr ) 10), it reduces to the exact spectrum formula and normalized radial eigenfunctions of the Klein-Gordon equation for vector and scalar general Hulthén-type potentials: which gives 16q m 0 η n m 0 E n = qr 0m 0 S 0 +E nl V 0 )+r 0 S0 V 0 ) n+δ, q n+δ) r 0 δ = 1 [ 1+ 1 ] q q +4r0 S 0 V 0 ), n = 0,1,,3, ) 88) E n = η ) nv 0 q + qαn+1)+ q α +4S0 V0 ) m 0 4V0 + qαn+1)+ ) q α +4S0 V 0 ) η n = [ 4V 0 + ) ηn, n = 0,1,,3, ) 4q ) ] qαn+1)+ q α +4S0 V0 ), q α +4S0 V0 ) 0, 4V0 S0)+ qαn+1)+ q α +4S0 V0 )) 8qm0 S 0 4V0 + qαn+1)+ ). 89) q α +4S0 V 0 ) iv) In the case S 0 = V 0 Hulthén potential, Eqs. 67) and 68) can be reduced to the relativistic energy equation: m 0 E R = r 0V 0 m 0 +E R ) n+δ, δ = n+δ) r 0 which is Eq. ) of Ref. [19]. D +l 1), n,l = 0,1,, ). 90) 16

v) We discuss non-relativistic limit of the energy equation 90). When V 0 = S 0, Eq. 10) reduces to a Schrödinger-like equation for the potential Vr). In other words, the non-relativistic limit is the Schrödinger equation for the potential V 0 e r/r 0 / [ 1 e r/r 0]. Hence, using the transformation m 0 + E R m 0 and m 0 E R E NR [], we obtain the non-relativistic energy equation: E NR = 1 8m 0 α [ 4m0 V 0 α n+δ) n+δ) ], α = r 1 0 91) which is Eq. 3) of Ref. [19] with δ is given in Eq. 90). It is worthwhile to remark that Eq. 91) is identical to Eq. 59) of Ref. [30] when the potential is Vr) and α becomes pure imaginary, i.e., α iα. Thus, in the weak coupling condition, [n+δ)/m 0 r 0 ] 1, [V 0 r 0 /n+δ)] 1, expanding the energy equation 90), retaining only the term containing the power of 1/m 0 r 0 ) and r 0 V 0 ) 4, we have the relativistic energy ) 4 V0 r 0 E R E NR +m 0 +4m 0, 9) qn+δ) which iseq. 4)ofRef. [19]with δ isgiven ineq. 90). Thefirst termisthenon-relativistic energy and third term is the relativistic approximation to energy. vi) For a more specific case where q = 1, the usual Hulthén potential is reduced to the shifted usual Woods-Saxon WS) potential Vr) = V 0 + V 0 1+e, Sr) = S S 0 r/r 0 + 0 1+e, r r/r 0 = α 1, 93) 0 and hence the energy eigenvalues, Eq. 69), for the general WS-type potentials, i.e., q q and then q = 1, become where E D) nl = ξ nlv 0 ± [ ] α +4S0 V 0 ) α D +l ) αn+1) ξnl ) m 0 [ α ], 4V0 + +4S0 V0 ) α D +l ) 4 αn+1) { [ ] } 16m 0 ξ nl 4V0 + α +4S0 V 0 ) α D+l ) αn+1), 94) ξ nl = [ 4V0 S0)+ α +4S0 V0 ) α D +l ) αn+1)] +8m0 S 0 4V 0 + [ α +4S 0 V 0 ) α D+l ) αn+1) ], 95) 17

and the inequality α +4S 0 4V 0 +α D +l ), 96) must be fulfilled. In the pure vector potential S 0 = 0,V 0 0), the energy eigenvalues are [ ] α α D +l ) 4V0 αn+1) E D) nl = V 0 ± m 0 4V 0 + [ α α D+l ) 4V 16m 0 4V 0 + 1 16, ] 0 αn+1) [ α α D +l ) 4V0 αn+1)], 97) where α 4V 0 +α D+l ). 98) From the above equation, for any given α, the Klein-Gordon equation with the usual shifted WS potential has negartive eigenvalues, i.e, E nl < 0. The given restriction in 99) imposes the critical coupling value and thus leads to the result n 1 α ) α 4 [ D +l ) ] V0 4m 0 V0 1, 99) i.e., there are only finitly many eigenvalues. Further, it is necessary that the inequality: must be fulfilled. α [ D +l ) ] 4V0 α+ 4m 0 V0, 100) Now, let us find the wave function y n s), which is the polynomial solution of hypergeometric-type equation, we multiply Eq. 19) by the weight function ρs) so that it can be rewritten in self-adjoint form [36] [ωs)y ns)] +λρs)y n s) = 0. 101) The weight function ρs) which satisfies Eqs. 4) and 10) has the form ρs) = s ε nl 1+4β 1 s) +γ), 10) and consequently from the Rodrigues relation 3), we obtain y nl s) = B nl s ε nl 1 s) [ 1+4β +γ) dn s n+ε nl 1 s) n+ ] 1+4β +γ) ds n 18

= B nl P ε nl, 1+4β +γ)) n 1 s). 103) On the other hand, inserting the values of σs),πs) and τs) given in Eqs. 59), 63) and 64) into Eq. 0), one can find the other part of the wave function as φs) = s ε nl 1 s) 1 Hence, the wave function in Eq. 18) becomes h gs) = C nl s ε nl 1 s) 1 1+ i 1+4β +γ) h i 1+ 1+4β +γ). 104) P ε nl, 1+4β +γ)) n 1 s) = C nl s εd) nl 1 s) δ P εd) n 1 s), s [q,0). 105) nl,δ 1) Finally, the radial wave functions of the Klein-Gordon equation are obtained as m 0 E nl,a/q) R l r) = N nl r D 1)/ e m 0 E nl r 1 qe r/r 0 ) q+a)/q) P r 0 n 1 qe r/r 0 ), 106) where N nl is the radial normalization factor. Thus, the radial wave functions for the s-wave Klein-Gordon equation with pure vector Hulthén potential in 1D reduces to m 0 E nl,a/q) R n r) = C nl e m 0 E nl r 1 qe r/r 0 ) q+a)/q) P r 0 n 1 qe r/r 0 ), 107) where a = q 4V 0 /α which is identical to Eq. 35) of Ref. [30]. Finally, from Eqs. 5) and 6), the total wave functions for the usual Hulthén potential is ψ l D 1=l) m 0 E nl,δ 1) l 1 l D x) = N nl r D 1)/ e m 0 E nl 1 e r/r 0 ) δ P r 0 n 1 e r/r 0 ) D 1 l j +j 1)n j! exp±il 1 θ 1 ) π Γl j= j +l j 1 +j ) sinθ j) lj nj P l j n j +j )/,l j n j +j )/) n j cosθ j ) n D 1 +m +1)n D 1! sinθ Γn D 1 +m D 1 ) l D P m,m ) n ) D 1 cosθ D 1 ), 108) where δ = 1 ) 1+ D+l ) +4r0 S 0 V 0 ), l = 0,1,, ). 109) V. CONCLUSIONS We have analytically found an approximate bound state energy eigenvalues and their corresponding wave functions of the D-dimensional Klein-Gordon equation for the spinzero particle with the scalar and vector general Hulthén-type potentials using the N-U 19

method. The analytic energy equation and the wave functions expressed in terms of Jacobi polynomials and can be reduced to their well-known 1D and 3D solutions. The relativistic energies E nl in D-dimensions defined explicitly in Eq. 68) is for the particle and antiparticle energies for a given dimension D, quantum numbers n and l and also coupling constants satisfying the given particular constraints. For pure attractive scalar potential, all bound states appear in pairs, with energies ±E nl. Since the Klein-Gordon equation is independent of the sign of E nl for scalar potential, the wave functions become the same for both energy values. When l = 0, the results in this work reduce to exact solution of bound states of s-wave Klein-Gordon equation with vector and scalar general Hulthén-type potentials. A preliminary analysis about the possibility of obtaining finite number states for the Klein-Gordon particle of mass unity shows that the bound state energies are positive for the general Hulthén-type potentials and negative for general Woods-Saxon-type potentials. In Table 1, we have obtained numerical results for the binding energies of ground state for mass unity 1D Klein Gordon equation with pure vector and pure scalar cases. In pure vector potential, for fixed coupling constants V 0 = 0.5 m 0 and various values of α = 0.5,1.0 and.0, most of the binding energies E 0 are decreasing with increasing deformationconstant q. Also the values obtained with pure scalar potential, S 0 = 0.5 m 0, are also given in Table 1 for comparison. Furthermore, in Table, we give the binding energies for the n = 1, α = 0.5,1.0 and n =, α = 0.5. Obviously, we have no bound states for n = 1 when α =.0 and for n = when α = 1.0 and.0. Acknowledgments Work partially supported by the Scientific and Technological Research Council of Turkey TÜBİTAK). 0

[1] T. Y. Wu and W. Y. Pauchy Hwang, Relativistic Quantum Mechanics and Quantum Fields World Scientific, Singapore, 1991). [] W. Greiner, Relativistic Quantum Mechanics: Wave Equations, 3rd edn springer, Berlin, 000). [3] A. D. Alhaidari, Phys. Rev. Lett. 87 001) 10405; 88 00) 189901. [4] G. Chen, Mod. Phys. Lett. A 19 004) 009; J. -Y. Guo, J. Meng and F. -X. Xu, Chin. Phys. Lett. 0 003) 60; A. D. Alhaidari, J. Phys. A: Math. Gen. 34 001) 987; 35 00) 607; M. Şimşek and H. Eğrifes, J. Phys. A: Math. Gen. 37 004) 4379. [5] J. -Y. Guo, X. -Z. Fang and F. -X. Xu, Phys. Rev. A 66 00) 06105; C. Berkdemir, A. Berkdemir and R. Sever, J. Phys. A: Math. Gen. 39 006) 13455. [6] G. Chen, Acta Phys. Sinica 50 001) 1651; Ö. Yeşiltaş, Phys. Scr. 75 007) 41. [7] G. Chen and Z.M. Lou, Acta Phys. Sinica 5 003) 1071. [8] G. Chen, Z. D. Chen and Z. M. Lou, Chin. Phys. 13 004) 79. [9] W. C. Qiang, Chin. Phys. 1 003) 136. [10] W. C. Qiang, Chin. Phys. 13 004) 571. [11] G. Chen, Phys. Lett. A 38 004) 116; Y. F. Diao, L. Z. Yi and C. S. Jia, Phys. Lett. A 33 004) 157. [1] L. Z. Yi et al, Phys. Lett. A 333 004) 1. [13] X. Q. Zhao, C. S. Jia and Q. B. Yang, Phys. Lett. A 337 005) 189. [14] A. D. Alhaidari, H. Bahlouli and A. Al-Hasan, Phys. Lett. A 349 006) 87. [15] G. Chen, Phys. Lett. A 339 005) 300. [16] A. de Souza Dutra and G. Chen, Phys. Lett. A 349 006) 97. [17] F. Dominguez-Adame, Phys. Lett. A 136 1989) 175; L. Chetouani et al., Physics A 34 1996) 59. [18] A. S. de Castro, Phys. Lett. A 338 005) 81. [19] C. Y. Chen, D. S. Sun and F. L. Lu, Phys. Lett. A 370 007) 19. [0] S. M. Ikhdair and R. Sever, Int. J. Mod. Phys. C 19 008) 1. [1] S. M. Ikhdair and R. Sever, Cent. Eur. J. Phys. 6 008) 141. [] S. M. Ikhdair and R. Sever, arxiv:0801.4857; to appear in the Int. J. Mod. Phys. C. 1

[3] S. M. Ikhdair and R. Sever, Cent. Eur. J. Phys.6 008) 685. [4] S. M. Ikhdair, Chin. J. Phys.46 008) 91. [5] S. M. Ikhdair and R. Sever, arxiv:0801.3518; arxiv:0801.471; S. M. Ikhdair and R. Sever, J. Mol. Struct.-Theochem 806 007) 155; S. M. Ikhdair and R. Sever, J. Mol. Struct.-Theochem 855 008) 13; S. M. Ikhdair and R. Sever, Cent. Eur. J. Phys. 5 007) 516; S. M. Ikhdair and R. Sever, Cent. Eur. J. Phys.6 008) 697. [6] S. M. Ikhdair and R. Sever, Int. J. Theor. Phys. 46 007) 1643; S. M. Ikhdair and R. Sever, Int. J. Theor. Phys. 46 007) 384. [7] S. M. Ikhdair and R. Sever, Ann. Phys. Leipzig) 16 007) 18. [8] S. M. Ikhdair and R. Sever, Int. J. Mod. Phys. E 17 6) 008) 1107. [9] S. M. Ikhdair and R. Sever, J. Math. Chem. 4 007) 461. [30] M. Simsek and H. Egrifes, J. Phys. A: Math. Gen. 37 004) 4379. [31] C. Berkdemir, Am. J. Phys. 75 007) 81. [3] W. C. Qiang, R. S. Zhou and Y. Gao, Phys. Lett. A 371 007) 01. [33] N. Saad, arxiv:0709.4014. [34] G. Chen, Z. D. Chen and P. C. Xuan, Phys. Scr. 74 0069 367. [35] S. Dong, J. Garcia-Ravelo and S. H. Dong, Phys. Scr. 76 007) 393. [36] S. M. Ikhdair and R. Sever, Ann. phys.berlin) 17 008) [11] DOI:10.100/andp.008103); arxiv:0807.085, submitted to J. Math. Chem. [37] Submitted to Ann. Phys. Berlin). [38] B. Gonul, O. Ozer, Y. Cancelik and M. Kocak, Phys. Lett. A 75 000) 38; M. Aktas and R. Sever, J. Mol. Struct. 710 004) 19. [39] R. L. Greene and C. Aldrich, Phys. Rev. A 14 1976) 363. [40] U. Myhrman, J. Phys. A: Math. Gen. 16 1983) 63. [41] A. Bechlert and W. Bühring, J. Phys. B: At. Mol. Opt. Phys. 1 1988) 817. [4] S. W. Qian, B. W. Huang and Z. Y. Gu, New J. Phys. 4 00) 13. [43] O. Bayrak, G. Kocak and I. Boztosun, J. Phys. A: Math. Gen. 39 006) 1151. [44] A. F. Nikiforov and V. B. Uvarov, Special Functions of Mathematical Physics Birkhauser, Bassel, 1988). [45] C. Y. Chen and S.H. Dong, Phys. Lett. A 335 005) 374. [46] Y. F. Cheng and T. Q. Dai, Phys. Scr. 75 007) 74.

[47] H. Egrifes and R. sever, Int. J. Theor. Phys. 46 007) 935. [48] G. Chen, Z. D. Chen and Z. M. Lou, Phys. Lett. A 331 004) 374. 3

TABLE I: Relativistic ground state binding energies, E 0, for a particle of mass, m 0 = 1 given by Eqs. 79) and 8) for pure vector and pure scalar potentials, respectively, as a function of q for various values of α. Pure vector V 0 = 0.5) a Pure Scalar S 0 = 0.5) q α = 0.5 α = 1.0 α =.0 α = 0.5 α = 1.0 α =.0 0.1 0.755639 0.947484 0.946410 0.5 0.911438 0.500000 0.9981 0.996314 0.64597 1.0 0.80971 0.970804 0.50000 0.98645 0.964541 0.480683 1.5 0.991673 0.940945 0.166683 0.998606 0.94146 0.398341.0 0.999840 0.93893 0.15000 0.999903 0.9656 0.34733.5 0.999140 0.913089 0.100000 0.998350 0.916087 0.311871 5.0 0.988667 0.890301 0.050000 0.988537 0.89966 0.136 7.5 0.98893 0.88371 0.033334 0.98996 0.884414 0.181787 10 0.979613 0.878345 0.05000 0.979771 0.879977 0.157607 a For real bound states q α 0.5. 4

TABLE II: Relativistic binding energies of the excited states, E 1 and E, for a particle of mass, m 0 = 1 given by Eqs. 79) and 8) for pure vector and pure scalar potentials, respectively, as a function of q for various values of α. n = 1 a n = n = 1 b n = q α = 0.5 α = 1.0 α = 0.5 α = 0.5 α = 1.0 α = 0.5 0.1-0.995674 0.771938 0.9413 0.5 0.830948-0.9840 0.57183 0.89156 1.0 0.99641 0.3479 0.880588 0.954903 0.4507 0.779514 1.5 0.94940 0.959 0.769589 0.935491 0.381304 0.75337.0 0.98534 0.171407 0.737131 0.9604 0.336188 0.735135.5 0.91607 0.136934 0.719697 0.913600 0.30388 0.73317 5.0 0.89105 0.06834 0.688449 0.898 0.19316 0.695608 7.5 0.8869 0.045546 0.678994 0.884103 0.18055 0.684996 10 0.87855 0.034156 0.674438 0.879801 0.156613 0.679406 a Pure vector case with coupling constant V 0 = 0.5. b Pure scalar case with coupling constant S 0 = 0.5. 5