^{1}

^{*}

^{2}

Recently, the Asymptotic Iteration Method (AIM) was used to calculate the energy spectrum for a short - rang e three parameters central potential which was introduced by H. Bahlouli and A. D. Alhaidari. The S-orbital wave solution of the Schrödinger equation was obtained for different parameters of the potential. In this work , a non-zero angular momentum term is introduced to the problem and the energy eigenvalues were obtained for different potential parameters. Our results show very good agreements compared with other methods such as potential parameter spectrum method (PPSM) and the complex scaling method (CSM).

In quantum mechanics, the most basic problem is to solve the Schrödinger equation for the energy eigenvalues E_{n} 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 [

Recently, the study of exponential-type potentials has attracted much attention from many authors [

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 [

V ( r ) = V 0 e − λ r − γ e λ r − 1 (1)

where V 0 is the potential strength, and the range parameter λ is positive with an inverse unit length. The dimensionless parameter γ is in the open range 0 < γ < 1 . This potential is short-range with 1/r singularity at the origin. It is also interesting to note that, at short distance and with 0 < γ < 1 , there is a clear resemblance of this potential with V 0 > 0 to the attractive Coulomb potential [

In a previous article [

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.

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 [

f ″ ( x ) = λ 0 ( x ) f ′ ( x ) + s 0 ( x ) f ( x ) (2)

where f ( x ) is a function of x, f ′ ( x ) , f ″ ( x ) are the first and second derivatives with respect to x; λ 0 ( x ) , s 0 ( x ) are arbitrary functions in C ∞ ( a , b ) and λ 0 ( x ) ≠ 0 . If we differentiate Equation (2) with respect to x, we find that;

f ‴ ( x ) = λ 1 ( x ) f ′ ( x ) + s 1 ( x ) f ( x ) (3)

where,

λ 1 ( x ) = λ ′ 0 ( x ) + s 0 ( x ) + λ 0 2 ( x ) , s 1 ( x ) = s ′ 0 ( x ) + s 0 ( x ) λ 0 (x)

And we continue with the same way for higher derivatives ( k + 1 ) th and ( k + 2 ) th where k = 1 , 2 , 3 , ⋯ is the iteration number, and we get,

f ( k + 1 ) ( x ) = λ k − 1 ( x ) f ′ ( x ) + s k − 1 ( x ) f ( x ) f ( k + 2 ) ( x ) = λ k ( x ) f ′ ( x ) + s k ( x ) f ( x ) (4)

where

λ k ( x ) = λ ′ k − 1 ( x ) + s k − 1 ( x ) + λ 0 ( x ) λ k − 1 ( x ) , s k ( x ) = s ′ k − 1 ( x ) + s 0 ( x ) λ k − 1 ( x ) (5)

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 ( k + 2 ) th and ( k + 1 ) th derivatives we get;

d d x ln [ f ( k + 1 ) ( x ) ] = f ( k + 2 ) ( x ) f ( k + 1 ) ( x ) = λ k ( x ) [ f ′ ( x ) + s k ( x ) λ k ( x ) f ( x ) ] λ k − 1 ( x ) [ f ′ ( x ) + s k − 1 ( x ) λ k − 1 ( x ) f ( x ) ] (6)

Assuming that for sufficiently large k;

s k ( x ) λ k ( x ) = s k − 1 ( x ) λ k − 1 ( x ) = α ( x ) (7)

holds, which is the “asymptotic aspect” of the method, Equation (6) reduces to

d d x ln [ f ( k + 1 ) ( x ) ] = f ( k + 2 ) ( x ) f ( k + 1 ) ( x ) = λ k ( x ) λ k − 1 ( x ) (8)

Substituting λ k ( x ) from Equation (5) and then using α ( x ) in the right hand side of Equation (8) we obtain;

f ( k + 1 ) ( x ) = C exp ( ∫ λ k ( x ) λ k − 1 ( x ) d x ) = C 1 λ k − 1 ( x ) exp ( ∫ x [ α ( x 1 ) + λ 0 ( x 1 ) ] d x 1 ) (9)

in which C 1 is the integration constant. Inserting Equation (9) into Equation (4) and solving for f ( x ) , we obtain the general solution of Equation (2) as

f ( x ) = exp ( − ∫ x α ( x ) d x 1 ) [ C 2 + C 1 ∫ x exp ( ∫ x 1 [ λ 0 ( x 2 ) + 2 α ( x 2 ) ] d x 2 ) d x 1 ] (10)

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

Δ k ( x ) = λ k ( x ) s k − 1 ( x ) − λ k − 1 ( x ) s k ( x ) (11)

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 x 0 point. The chosen value of x 0 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 [

For a given energy E and angular momentum l, time-independent radial Schrödinger equation for the reduced radial wave function, R ( r ) , can be written as

( − 1 2 d 2 R ( r ) d r 2 + l ( l + 1 ) 2 r 2 R ( r ) ) + V ( r ) R ( r ) = E R ( r ) (12)

where the potential V ( r ) given in equation (1); and we have adopted the atomic units h = m = 1 .

In order to solve Equation (12) for l ≠ 0 , we need to apply the approximation scheme to the centrifugal term given by [

1 r 2 ≈ λ 2 ( 1 12 + e λ r ( e λ r − 1 ) 2 ) (13)

Substituting Equation (13) into Equation (12), we obtain the following equation

( − 1 2 d 2 R ( r ) d r 2 + λ 2 l ( l + 1 ) 2 ( 1 12 + e λ r ( e λ r − 1 ) 2 ) R ( r ) ) + V ( r ) R ( r ) = E R ( r ) (14)

Now introducing the change in variable x = 1 − 2 e − λ r whose range is between −1 and +1, transforms Equation (14) into the desired Equation (1) where we can apply the AIM with

λ 0 ( x ) = 1 1 − x (15)

s 0 ( x ) = 1 12 λ 2 ( x + 1 ) 2 ( x − 1 ) 2 ( 12 V 0 x 3 + η x 2 + ξ x + ν ) (16)

where

η = l ( l + 1 ) λ 2 + 12 ( 2 γ − 1 ) V 0 − 24 E (17)

ξ = − 22 l ( l + 1 ) λ 2 − 12 V 0 − 48 E (18)

And

ν = 25 l ( l + 1 ) λ 2 + ( 12 − 24 γ ) V 0 − 24 E (19)

In order to obtain the energy eigenvalues from equation the (14), using Equation (5), we obtain the λ k ( x ) , and s k ( x ) in terms of λ 0 ( x ) , and s 0 ( x ) . 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.

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 _{0} and E_{1} are very deeply bounded.

l | n | E_{n} by CSM | E_{n} by AIM |
---|---|---|---|

0 | 0 | −1094.42109160 | −1094.42109160 |

1 | −187.97359168 | −187.973591677 | |

2 | −36.02622806 | −36.0262280564 | |

1 | 0 | −185.38841241 | −185.3881967 |

1 | 33.95322592 | −33.95168758 | |

2 | 0 | −29.64332195 | −29.63997992 |

l | n | E_{n} by CSM | E_{n} by AIM |
---|---|---|---|

0 | 0 | −1406.11040577 | −1406.11040577 |

1 | −223.29635015 | −223.296350148 | |

2 | −27.18320883 | −27.1832088304 | |

1 | 0 | −219.66959141 | −219.6694254 |

1 | −24.21918006 | −24.21795546 | |

2 | 0 | −18.01714564 | −18.01447261 |

l | n | E_{n} by CSM | E_{n} by AIM |
---|---|---|---|

0 | 0 | −679.95986643 | −679.959866430 |

1 | −32.96147955 | −32.961479552 | |

1 | 0 | −27.29186980 | −27.29153038 |

l | E_{0} | E_{1} | E_{2} | E_{3} |
---|---|---|---|---|

0 | −2.948040331 × 10^{5} | −7813.547139 | −1748.849560 | −640.0248682 |

1 | −2.859760644 × 10^{5} | −2.428306994 × 10^{5} | −637.9431013 | −265.3244481 |

2 | −633.7598994 | −261.5672734 | −33.72115660 | −4.266342233 |

3 | −255.8552832 | −99.68104017 | −29.72230919 | −1.541952406 |

4 | −92.89079567 | −24.18745837 | 2.003775523 | 2.116254874 |

5 | −16.85605236 | 2.897862067 | 4.531566676 | 6.360648382 |

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

Sous, A.J. and El-Kawni, M.I. (2018) Numerical Solution of the Schrödinger Equation for a Short-Range Potential with Any ℓ Angular Momentum. Journal of Applied Mathematics and Physics, 6, 901-909. https://doi.org/10.4236/jamp.2018.64077