In quantum mechanics, the most basic problem is to solve the Schrödinger equation for the energy eigenvalues En and the associated energy Eigen functions.
There are a number of important cases for which the stationary Schrödinger equation can be solved analytically. However, analytic solutions are possible only for a few simple quantum systems such as the hydrogen atom, the harmonic oscillator and others   . In most cases, many quantum systems can be treated only by approximation methods. Among such approximation methods include Pekeris approximation  , Semi-classical method  , and asymptotic iteration method  -  .
Recently, the study of exponential-type potentials has attracted much attention from many authors    . These potentials include the Hulthén potential  , the multi parameter exponential-type potentials  , the Manning-Rosen potential  and the Eckart-type potential  . It should be mentioned that most contributions appearing in the literature are concerned with the s-wave case.
In this work, we study the arbitrary l-state solutions of the Schrödinger equation with a short rang three parameters central potential which was introduced by H. Bahlouli and A. D. Alhaidari  
where is the potential strength, and the range parameter λ is positive with an inverse unit length. The dimensionless parameter γ is in the open range . This potential is short-range with 1/r singularity at the origin. It is also interesting to note that, at short distance and with , there is a clear resemblance of this potential with to the attractive Coulomb potential  with non-zero angular momentum. The potential valley here is not due to the centrifugal force attributed to the angular momentum. Moreover, it does not have the long-range behavior of the Coulomb potential   . In  , Alhaidari argued that in contrast to the Coulomb potential, the number of bound states for this potential is finite and that it could be used as a more appropriate model for the description of an electron interacting with an extended molecule whose electron cloud is congregated near the center of the molecule. The authors of    found the “potential parameter spectrum” (PPS) for the hyperbolic single wave potential and for potential (1). The concept of a PPS was introduced for the first time in the solution of the wave equation in  where for a given energy, the problem becomes exactly solvable for a discrete set (finite or infinite) of values of the potential parameters. If the map that associates the parameter spectrum with the energy is invertible, then in principle one could obtain the energy spectrum for a given choice of potential parameters  .
In a previous article  , we used the Asymptotic Iteration Method (AIM) to find the energy spectrum for the s-wave (zero angular momentum) Schrodinger Equation with potential in (1) which was introduced by Bahlouli and Alhaidari   . In the present work, we apply the same technique in  -  to the same potential in (1), and we will find the eigenvalues for the time-independent radial Schrödinger equation for any angular momentum l.
The paper has the following structures. In Section 2, we briefly present an overview of the AIM which introduced to find the solutions for the second-order differential equation. In Section 3, change of variables and approximation scheme has been done which allows as transforming Schrödinger equation to another form in order to apply the method to solve the equation with the short rang three parameters central potential. In Section 4, our numerical calculation results have been presented for the eigenvalues.
2. Concept of Asymptotic Iteration Method
In this section we shall outline the general procedure of the AIM for determining the eigenvalues differential equation. The AIM has been proposed and used to solve the homogenous linear second-order differential equation of the form  -  ,
where is a function of x, are the first and second derivatives with respect to x; are arbitrary functions in and . If we differentiate Equation (2) with respect to x, we find that;
And we continue with the same way for higher derivatives and where is the iteration number, and we get,
Iterations will continue until we get the convergence condition of the AIM  .
Equation (5) is called the recurrence relation of Equation (2). Taking the ratio of the and derivatives we get;
Assuming that for sufficiently large k;
holds, which is the “asymptotic aspect” of the method, Equation (6) reduces to
Substituting from Equation (5) and then using in the right hand side of Equation (8) we obtain;
in which is the integration constant. Inserting Equation (9) into Equation (4) and solving for , we obtain the general solution of Equation (2) as
The energy eigenvalues can be determined by the quantization condition given by the termination condition in Equation (7). Thus, one can write the quantization condition combined with Equation (5) as
If the eigenvalues problem is an exactly solvable problem, an explicit expression for energy can directly be obtained from the roots of this equation that depends only on the eigenvalues E. In this case, for any given quantum number n, the energy eigenvalues can be calculated from the roots of the Equation (11) at some suitable point. The chosen value of is arbitrary in principle and can be critical only to the speed of the convergence of the method. This starting value may be determined generally as the minimum value of the potential or the maximum value of the asymptotic wave function  .
3. Solution and Energy Spectrum
For a given energy E and angular momentum l, time-independent radial Schrödinger equation for the reduced radial wave function, , can be written as
where the potential given in equation (1); and we have adopted the atomic units .
In order to solve Equation (12) for , we need to apply the approximation scheme to the centrifugal term given by 
Substituting Equation (13) into Equation (12), we obtain the following equation
Now introducing the change in variable whose range is between −1 and +1, transforms Equation (14) into the desired Equation (1) where we can apply the AIM with
In order to obtain the energy eigenvalues from equation the (14), using Equation (5), we obtain the , and in terms of , and . Then, using the quantization condition of the method given by Equation (11), we obtain the energy eigenvalues. This straightforward application of AIM gives us the energy eigenvalues. We have observed that the energy eigenvalues converge within a reasonable number of iterations. This result agrees with the principle of AIM; as the number of iteration increases, the method should converge and should not oscillate.
4. Results and Discussion
The arbitrary l-state solutions of the Schrödinger equation with short rang three parameters central potential have been obtained. To show the accuracy of our results, we have calculated the eigenvalues numerically for many l-states and found that the results obtained by Equation (11) are in good agreement with those obtained by other methods.
Our numerical results for the short rang three parameters central potential are listed in Tables 1-3. They show a comparison between the CSM and AIM eigen states. The results are in full agreement in the ground states, and agree very well for the higher states.
We used the data in Table 4 to draw the relation between the angular momentum l and the corresponding eigen states. Figure 1 shows that the ground states E0 and E1 are very deeply bounded. Figure 2 shows the same behavior for
Table 1. A comparison of the energy eigenvalues En of the potential (1) obtained here by the AIM and compared to those obtained by the complex scaling Method (CSM) in  . We took , , and for various values of the angular momentum l.
Table 2. A comparison of the energy eigenvalues En of the potential (1) obtained here by the AIM and compared to those obtained by the complex scaling Method (CSM) in  . We took , , and for various values of the angular momentum l.
Table 3. A comparison of the energy eigenvalues En of the potential in equation (1) obtained here by the AIM and compared to those obtained by the complex scaling Method (CSM) in  . We took , , for various values of the angular momentum l.
Figure 1. The energy eigenvalues as function of angular momentum l for the gound (Eo) and first (E1) excited states.
Table 4. The energy eigenvalues for different angular momentum l; calculated by AIM. The values of the potential parameters are , , .
Figure 2. The energy eigenvalues as function of angular momentum l for the second (E2) and third (E3) excited states.
higher eigen states. As the angular momentum increases the energy eigen states increased until l is greater than 2, we reached like a continuous state or we can say that the potential is not affecting any more.
We used the AIM as a powerful engine to solve the Schrödinger equation for a short range three parameters potential with any angular momentum values. Our results are considered to be excellent compared to the CSM. We also introduced a new idea by studying the change of angular momentum with energy eigenvalues and clarifying the relation by Figure 1 and Figure 2. This relation was not being able to reach by other methods especially for higher l values.
 Ciftci, H., Hall, R.L. and Saad, N. (2003) Asymptotic Iteration Method for Eigen value Problems. Journal of Physics A: Mathematical and Theoretical, 36, Article ID: 11807.
 Ozer, O. and Roy, P. (2009) The Asymptotic Iteration Method Applied to Certain Quasi Normal Modes and Non Hermitian Systems. Central European Journal of Physics, 7, 747-752.
 Soylu, A., Bayrak, O. and Boztosun, I. (2007) An Approximate Solution of Dirac-Hulthén Problem with Pseudo Spin and Spin Symmetry for Any k state. Journal of Mathematical Physics, 48, 082302.
 Sous, A.J. and El-Kawni, M.I. (2009) General Eigenvalue Problems with Unbounded Potential from Below. International Journal of Modern Physics A, 24, 4169.
 Barakat, T. (2005) The Asymptotic Iteration Method for the Eigenenergies of the Anharmonic Oscillator Potential V(x) = Ax2a + Bx2. Physics Letters A, 344, 411-417.
 Arda, A. and Sever, R. (2009) Exact Solutions of Effective Mass Dirac Equation with Non-PT-Symmetric and Non-Hermitian Exponential-Type Potentials. Chinese Physics Letters, 26, 90305-090305.
 Soylu, A., Bayrak, O. and Boztosun, I. (2008) Exact Solutions of Klein-Gordon Equation with Scalar and Vector Rosen-Morse-Type Potentials. Chinese Physics Letters, 25, 2754-2757.
 Taskin, F. and Kocak, G. (2010) Approximate Solutions of Schrodinger Equation for Eckart Potential with Centrifugal Term. Chinese Physics B, 19, 090314-6.
 Sous, A.J. and Alhaidari, A.D. (2016) Energy Spectrum for a Short-Range 1/r Singular Potential with a Non-Orbital Barrier Using the Asymptotic Iteration Method. Journal of Applied Mathematics and Physics, 4, 79-85
 Alhaidari, A.D. (2007) Representation Reduction and Solution Space Contraction in Quasiexactly Solvable Systems. Journal of Physics A: Mathematical and Theoretical, 40, 6305.