www.ccsenet.org/apr Applied Physics Research Vol. 2, No.

2; November 2010

Approximate Solution of the Schrödinger Equation with Rosen-Morse

Potential Including the Centrifugal Term

Akpan N. Ikot
Department of Physics, University of Uyo, Nigeria
E-mail: [email protected]

Louis E. Akpabio
Department of Physics, University of Uyo, Nigeria
E-mail: [email protected]

Abstract
We derive approximate analytical solutions of the Schrödinger equation with Rosen-Morse potential via the
Nikiforov-Uvarov method. The bound state energy eigenvalues are given in a closed form and the corresponding
eigenfunctions are obtained in terms of the generalize Jacobi Polynomials and hypergeometrical function.
Keywords: Nikiforov-Uvarov methods, Schrödinger equation, Rosen-Morse potential
1. Introduction
The exact analytical solution of the Schrödinger equations are only possible with the angular momentum l=0 for
some potential models. However, when l  0 , the Schrödinger equation can only be solve approximately for
different suitable approximation scheme. One of such approximations are the convectional approximation
scheme suggested by Greene and Aldrich and the new improved approximation scheme.
Many authors have investigated the bound state solution of Schrödinger approximately with exponential-line
potential such as Eckart Potential, Hulthen potential and Manning Rosen Potential.
Presently, an alternative method for solving Schrödinger equation has been proposed. The method called
Nikiforov-Uvarov (NU) is based on solving the second order linear differential equations by reducing to a
generalize equation of hypergeometric type. This method had been used to solve the Schrödinger, Klein-Gordon,
Dirac and Duffin-Kemner Petiau equations with exponential-like potentials such as Woods-Saxon, Huthen and
Poschl-Teller, and generalize Posch-Teller.
In this paper, we attempt to solve approximately the Schrödinger equation with Rosen-Morse potential including
the centrifugal term using the NU-Method. This potential plays a vital role in atomic, molecular and chemical
physics since it can be used to describe molecular vibrations and to obtain the energy spectra of linear and
nonlinear systems. The Rosen-Morse potential is one of the candidate potential used in describing the eigen
function and energy spectra of quark confinement in quantum chromodynamics(QCD). The QCD prescribes the
baryons to be constituted of three quarks in a colour singlet state. Bijker e tal exploits for the description of
baryonic systems algebraic methods developed for the purposes of triatomic molecules. However an interesting
dynamics limit of the three quark system is the one where two of the quarks are squeezed to one independent
entity, a di-quark (QQ), while the third quark (Q) remains spectator. In that limit which corresponds to U
(7)U(3)U(4), one can exploit the chain of reducing U(4) down to O(2). As a principle challenge still
remaining unanswered is finding a suitable quark potential that explain the above scenario. The Rosen-morse
potential with trigonometric function was use in Ref. [C. S. Jia, 2003] to explain the QCD quark-gluon dynamics.
In their report they show that the shape of the potential capture the essentials of QCD quark-gluon dynamics
where one gluon exchange gives rise to an effective coulomb-like potentials while the self-gluon interactions
produce a linear confinement potential as established by lattice calculation of hadrons properties. It was also
observed that the infinite wall shape of the trigonometric Rosen-Morse potential provides the regime suited for
trapped but asymptotically free quarks. With this behaviour of this potential, one is tempted to conclude that this
potential may be a good candidate for a quark potential of QCD traits.

202 ISSN 1916-9639 E-ISSN 1916-9647

www.ccsenet.org/apr Applied Physics Research Vol. 2, No. 2; November 2010

The organization of the paper is as follows. In section two we review the Nikiforov-Uvarov method. We find the
bound solution of the Schrödinger with Rosen-Morse potential in section three. A brief conclusion is presented
in section four.
2. Review of Nikiforov-Uvarov Method
The NU method is based on the solution of a generalized second order differential linear equation with special
orthogonal functions. Therefore, a non-relativistic Schrödinger equation can be solved exactly using this method.
For any given real or complex potential, the Schrödinger equation reduces to a generalized equation of
hypergeometric type with an appropriate, s = s(x) co-ordinate transformation. Thus, it can be written as follows:
 (s)  (s)
  ( s )   ( s )  2  (s)  0, (1 )
 (s)  (s)
where (s) and  (s) are polynomials at most second degree and  (s ) is a first degree polynomials.
In order to find the exact solution to equation (1), we set the wave function as

 ( s )   ( s )  ( s ), ( 2)

and on substituting equation (2) into equation (1) reduces equation (1) into hypergeometric type,

 ( s )  ( s )   ( s )  ( s )   ( s )  0, (3)

where the wave function (s) is defined as the logarithmic derivative [G. Szego, 1939]
 ( s )  ( s )
 , (4)
 (s)  (s)
where (s) is at most first degree polynomials.
Also, the hypergeometric type function (s) in equation (3) for a fixed n is given by the Rodriques relation
Bn d n
 n (s)   (s)  (s), (5)
 ( s ) ds n
Where Bn is the normalization constant and the weight function (s) must satisfy the condition [N. Cottas, 2004]

d
 ( s )  ( s )    ( s )  ( s ), (6 )
ds

with

 ( s)   ( s)  2 ( s). (7)
To accomplish the conditions imposed on the weight function (s), it is necessary that the classical orthogonal
polynomials (s) be equal to zero to some point of an interval (a, b) and its derivative at this interval at (s) > 0
will be negative, that is
d ( s)
 0. (8)
ds
Therefore, the function  (s ) and the parameter  required for the NU-method are defined as follows:
2
 1    
 (s)        k (9 )
2  2 

  k   (s). (10)
The k-values in the square-root of equation (9) is possible to evaluate if the expression under the square root
must be square of polynomials. This is possible if its discriminate is zero. Thus, the new eigen-value equation for
the Schrödinger equation becomes,
nd n(n  1) d 2
  n    , n  0,1,2,... (11)
ds 2 ds 2

Published by Canadian Center of Science and Education 203

www.ccsenet.org/apr Applied Physics Research Vol. 2, No. 2; November 2010

where (s) is as defined in equation (7) and on comparing equations (10) and (11), we obtain the energy eigen
values.
3. Solutions of the Schrödinger Equation with the Rosen-Morse Potential
The generalized Rosen-Morse potential is given by

v( x)  v1 sech12 (r)  v2 tan(r), (12)

displayed in fig.1.The shape of this potential also captures the essentials behaviour of the QCD quark-gluon
dynamics. The hyperbolic functions are defined as
e r  e  r e r  e r
sinh r  , cosh r 
2 2
sinh  r
tanh  r . (13 )
cosh  r
The radial Schrödinger equation with Rosen-Morse potential is

 2  d 2 2 d l (l  1) 
     nl (r )
2  dr 2 r dr  2 r 2 


 v 1 sec h 2 ( r )  v 2 tanh   nl ( r )  E nl  nl ( r ), (14 )

where l is the angular momentum number,  is the reduced mass and l is the centrifugal term.
r2
Now introducing a new variable for the wave function as
R(r )
 nl (r )  , (15)
r
reduces equation (14) into

d 2 R (r ) 2    2 l (l  1) 
  E  v 1 sec h 2
( r )  v 2 tanh(  r )   R (r )  0 (16 )
dr 2 2  2 r 2 

Due to the centrifugal term in equation (16), this equation cannot be solved analytically when the angular
momentum quantum number is not equal to zero.
Fortunately, approximate analytical solutions of with l  0 has been devised to deal with centrifugal term.
When r  1 , we use a new improved approximation scheme for the centrifugal term as [Y, Xu, 2010; C. S. Jia,
2009]

1  e  2r  2 1 
 4 2
 c    4 c 0  , (17)
r2 
0

1  e  2r 2
  4 sinh r 
2

where the parameter c  1 is a dimensional constant.

0
12

The convectional approximate scheme is obtain by setting c0=0 in equation (17). In order to apply the
NU-method, we substitute equation (17) into equation (16) and using a new variable of the form.

s  tanh(r ), (18)
we obtain the following,

204 ISSN 1916-9639 E-ISSN 1916-9647

www.ccsenet.org/apr Applied Physics Research Vol. 2, No. 2; November 2010

d 2

2s
  2
 
  2 s 2  r 2 s   2   2   2  2  0,  (19a)
ds 2

1 s2 2

where

2E 2v
2  , E  0,  2  2 12 , (19b)

2 2

2v 2
2  ,  2  4l (l  1)c0 , 2  l (l  1)
 2 2
are dimensionless parameters.
On comparing equation (19) with equation (1), we get

 ( s)  2s,  ( s)  1  s 2 ,

 ( s )    2   2 s 2   2 s   2   2   2   2 . ( 20 )

Substituting equations (20) into equation (9), we get

 ( s )     2   2  k s 2   2 s   2   2   2   2  k  ( 21)
The expression in the square root must be square of polynomials in respect of the NU-method. Therefore, we can
determine the constant k by the condition that the discriminate of the square root be zero, that is

K  
2  2
 2  2  2
 2
   2
 2

2
 4
, ( 22 )
2
In view of equation (22), we find new possible functions of  (s) for each k as


 2

 2  2 s  2  2  2

for , k
2 
 2 2   2   2 
2
  2

 2  4
 2
 ( s ) z  12 
  2
 2  2 
s   2   2   2 , for
( 23)


 for , k
2  2
 2 2   2   2    2
 2 2
 4
 2

For the polynomial of     2 ( s ) , which has a negative derivative, we get

k
2 2
 2   2   2    2
 2 
2
 4
, (24)
2

 ( s )   12  2
 2  2 s   1
2  2  2  2 (25)

Now using   k   (s) together with the values k and  ( s ),  , we obtain


 ( s)  2s  12  2   2   2 s  12  2   2   2  (26)


2 2
 22   2   2    2
 2 
2
 4

1
 2  2  2 (27)
2 2

Published by Canadian Center of Science and Education 205

www.ccsenet.org/apr Applied Physics Research Vol. 2, No. 2; November 2010

Also from the definition of  n , we find

 n  2n  n 2  2   2   2  n(n  1). ( 28)

Comparing equations (27) and (28), we obtain approximately the energy spectrum of the radial Schrödinger
equation with Rosen-Morse potential as

  2

  2  v4   2

 2  2 n  2  2  2  
 
  2   2  2 2  22  2n(n  1)  4n, (29)

or simplifying equation (29) yields

 2
 2   22
2 2
 22  2n(n  1)  4n 
 3(n  1)   2   2    2  3(n  1) 2 


 2 2  22  2n(n  1)  4n 2
(30)

Solving equation (30) explicitly, we obtain the energy spectrum as

 2 2  4v1
E  2 2  2(  1)1  2c0 
2   

 2n(n  1)  3 2 n  1  1
2  4n 

 4 v
  2 12  2[  1]  2n(n  1)  4n  3 2 (n  1)  1 2
 
1

  2 
2
2v  4  v1 
 2 22 1  3 ( n  1)       2  (   1 )  2 ( n  1 )  4 n   , ( 31 )
     
2 2
  

where we have use equation (19b) in obtaining (31) explicitly.
From equation (6), we obtain the weight function  (s) as

 ( s )  1  s 
1  2  2  2
1  s   1  2   2  

2  1
 4 , ( 32 )
4
and on substituting equation (32) into the Rodriques relation of equation (5), we get,

 1 4  1  2  2  r  1 
4  1  2  2  r 
 ( s )  Bn (1  s )  
(1  s )  

dn 1  s 2 1  s  14 1  2  2  r 
1  s  
  2  2  r 
4  1
1
x   
ds n 

   , 
 Bn Pn ( s) , (33)

206 ISSN 1916-9639 E-ISSN 1916-9647

www.ccsenet.org/apr Applied Physics Research Vol. 2, No. 2; November 2010

 2  2  r 2
where   14 ,   1
4 and Pn is the Jacobi polynomials.
The remaining part of the wave function in equation (4) is found as

 ( s )  1  s 2    v  1  s 2    v  (34)

Combining the Jacob polynomials of equation (33) and equation (34), we find the s-wave function as
  v    v      ,  
R n ( s )  C n (1  s ) 2 (1  s ) 2 Pn (s) ( 35 )


where Cn is the new normalization constant and it is determined by the condition

Rn ( s)ds 1 .

4. Conclusion
In this research work, we investigated the bound state solutions of the Schrödinger equation with Rosen-Morse
potential for any arbitrary angular momentum number l approximately. With the new improved scheme for the
centrifugal term, we obtained the energy eigenvalues and the unnormalized wave function in terms of the
hypergeometric functions via the Nikiforov-Uvarov method. The shape of this potential captures the essential
traits of the quark-gluon dynamics of QCD and it is a good choice that can be employed in the calculations of
spectroscopic properties of non-strange resonances.
Finally, our results are in consistent with those of refs. [C. S. Jia, 2003; 2002].
Acknowledgement
We acknowledge the financial support of Nandy-Leabio grants No. NL5563.
References
A. Berkedemir, C. Berkdemir and R. Sever. Quant-Ph/0410153.
A. F. Nikiforov, and V. B. Uvaroo. (1988). Special Function of Mathematical Physics, Birkhausev, Basel.
C. B. Compean and M. Kirchbach. (2006). Proceedings of the Mini-Workshop-Progress in Quark Models,
Angular Momentum Dependent Quark Potential of QCD Traits and Dynamical O(4) Symmetry, (Bled,
Slovenia).
C. Berkdemir, A. Berkdemir and R. Sever. (2004). Arxir: nucl-ph/0412021.
C. Berkedemir, A-Berkdemir and R. Saver. (2005). Phys. Rev., C 72, 027001.
C. S. Jia, Y. Li, Y. Sun, J. Y. Liu, L. T. sun. (2003). Phys. Lett. A 311, 115.
C. S. Jia, T. Chen and L. G. Cui. (2009). Phys. Lett. A 373, 1621.
C. S. Jia, X. L. Zeng, L. T. Sun. (2002). Phys. Lett., A 294. 185.
F. Taskin. (2009). Int. J. Theor. Phys., 48, 1142.
G. Szego. (1939). Orthogonal Polynomials, (American Mathematical Society, New York.
H. Egrifes, D. Demirhan and F. Buyukkihc. (1999). Physics Scripta, 59, 90.
H. Egrites, D. Dimerhan and F. Buyukkihc. (1999). Physics scripta, 195.
L. Chetouchi, L. Guechi, A. Lecheheb, T. F. Hamman and A. Messouber. (1996). Physic A, 234, 529.
L. H. Zhang, X. P. Li and C. S. Yia. (2008). Phys. Lett., A 372, 2201.
N. Cottas. (2004). CEJP, 2, 456.
N. Rosen and P. M. Morse. (1932). Phys. Rev., 42, 192.
O. Yesiltas, M. Simsek, R. Sever and C. Tezcan. (2003). Phys. Scripta, T67, 472.
R. Biyker, F. Iachello, A. Leviatan, Yu:F Smirnov. (2001). Phys. Rev. D44, 114005.
R. L. Greene and C. Aldrich. (1976). Phys. Rev., A. 14, 2363.
S. Flugge. (1974). Practical Quantum Mechanics, (Springer, - Verlag, Berlin).
V. H. Badalov, H. I. Almadov and A. I Almadov. (2009). Arxio: Math-ph/0905.273.

Published by Canadian Center of Science and Education 207

www.ccsenet.org/apr Applied Physics Research Vol. 2, No. 2; November 2010

W. C. Qiang and S. H. Dong. (2008). Phys. Lett., A 372, 4789.

Y, Xu, S. He and C. S. Jia. (2010). Phys. Scripta, 81, 045001.

Figure 1. The variation of Rosen-Morse potential V(r) with radius r. The shape of the potential shows a
typical QCD quark-gluon dynamics

208 ISSN 1916-9639 E-ISSN 1916-9647