The Asymptotic Iteration Method for the Eigenenergies of the a Novel Hyperbolic Single Wave Potential

By using the asymptotic iteration method, we have calculated numerically the eigenvalues n E of the hyperbolic single wave potential ( ) ( ) ( ) x V x V x λ γ λ 0 2 tanh cosh + = which is introduced by H. Bahlouli, and A. D. Alhaidari. They found a new approach (the “potential parameter” approach) which has been adopted for this eigenvalues problem. For a fixed energy, the problem is solvable for a set of values of the potential parameters (the “parameter spectrum”). This paper will introduce a related work to complete the goal of finding the eigenvalues, the Schrödinger equation with hyperbolic single wave potential is solved by using asymptotic iteration method. It is found that asymptotically this method gives accurate results for arbitrary parameters 0 V , γ , and λ .


Introduction
The numerical solutions of the Schrödinger equation have taken a great deal of interest in the description of nuclear motion in molecules in the Born-Oppenheimer approximation or in the description of atoms and molecules in self-consistent field approximations.In most cases the use of numerical methods is necessary, where analytical solutions are not possible or difficult to find and the number of exactly solvable problems is very limited [1].Over the past decade, several numerical methods have been developed for the solution of the Schrödinger equation, such as the Nikiforov-Uvarov (NU) method, Runge-Kutta methods, Matrix Diagonaliza-tion Method (MDM), Discrete Variable Method (DVM), Spectral Method [2]- [6], and the asymptotic iteration method (A.I.M).
Ever since the inception of empirical Morse potential about 85 years ago, a vast number of potential functions have been reported for molecules, with varying degrees of flexibility and accuracy [7].In this work we are interested to calculate the eigenvalues of hyperbolic single wave potential introduced by H. Bahlouli and A. D. Alhaidari [8] [9] some years ago.The new solvable hyperbolic single wave potential introduced given as where 0 V and γ are real potential parameters and λ is a length scale that determines the range of the po- tential.The first term of this potential is completely new and cannot be predicted by any of the classical approaches based on diagonal representations.However, the second part is a special case of the hyperbolic Pöschl-Teller potential that has an exact conventional (diagonal representation) solution.Physically, the most interesting situation is when the parameter γ lies between 1 − and 1 + in which case the shape of this potential becomes a hyperbolic single wave about the origin [8] [9].The potential has two local extrema at ( ) They have dealt with the case 1 γ < which carries a richer structure for a given potential.
In [8] they have dealt with case 1 γ < which carries a richer structure for a given potential.The configura- tion of this potential allows for both resonances and bound states, which is contrary to the case 1 γ ≥ where either bound or resonance energy states are allowed.Nonetheless, their approach can still handle the case 1 γ ≥ as long as 0 0 V γ < .In this work, we apply the asymptotic iteration method, in order to find eigenvalues of hyperbolic single wave potential.
This paper is divided into the following sections.In Section 2, we describe the asymptotic iteration method which introduced to find the solutions for the second-order differential equation.In Section 3, change of variables has been done which allows as transforming Schrödinger equation to another form in order to apply the method to solve the equation with hyperbolic single wave potential.In Section 4, our numerical calculation results have been presented for the eigenvalues of the hyperbolic single wave potential.

Solution Methodology
Recently a technique called the asymptotic iteration method (AIM) has been introduced [10] to obtain eigenvalues of second-order homogeneous differential equations.In the case of the Schrödinger equation the AIM was found to reproduce the energy spectrum exactly for most exactly solvable potentials [11]- [13], while for nonexactly solvable potentials it produces very good results [14]- [18].
The first step in applying this method to solve Schrödinger-type equations is to transform these equations, with the aid of appropriate asymptotic forms, to second-order homogeneous linear differential equations of the general form where the primes of ( ) y x denote derivatives with respect to x, ( ) . The solution of (3) follows from the symmetric structure of the right-hand side of Equation (3).Therefore, if we differentiate (3) with respect to x, we get where By differentiating Equation (4), we get for which In general, the ( ) and ( ) Respectively, where The ratio of the ( ) and ( ) derivatives can be expressed as For sufficiently large j, we can introduce the "asymptotic" aspect of the method, that is, if Thus, Equation ( 13) can be reduced to Which yields where 1 C is the integration constant, and the right-hand side of Equation ( 13) follows from Equation ( 9) and the definition of ( ) . Substituting Equation (13) into Equation ( 8) we obtain a first-order differential equation Which, in turn, yields the general solution to Equation (3)

Formulation of the Problem
The time-independent Schrödinger equation, in the field of potential ( ) V x can be written as where ( ) is the hyperbolic single wave potential in Equation (1) Substituting Equation (1) into Schrödinger equation ( 16) we obtain Before applying AIM to this problem, we have to transform the modified Schrödinger equation ( 17) to an amenable form for AIM.By using the change of variables, we can simplify Equation ( 17 AIM can be applied to the final form of the Schrödinger equation ( 19).
In order to calculate the eigenvalues n E , 0,1, 2,3, n =  we need to iterate the expansion In each iteration, the expansion will depend on two variables E and r.The calculated eigenvalues n E by the mean of the condition ( ) 0 r ϕ = should depend on the choice of r.Nevertheless, the choice of r is observed to be critical only to the speed of the convergence of the eigenvalues n E , as well as for the stability of the process.In our work, it is observed that the best starting value for r must be 0 r = .Therefore, at the end of the iterations we put 0 r = .

Results and Discussions
In this section, the paper shows the results of applying and using of AIM numerically, also the results of the calculations show that the AIM is powerful and the spectrum of energy for hyperbolic single wave potential can be obtained for all energy states.The resulting energy levels of the AIM are not limited by the magnitude of the parameters 0 , V γ and λ of the potential.Meanwhile, in the paper which has been done by H. Bahlouli, and A. D. Alhaidari found out that it is not so easy to obtain the energy spectrum for the potential by using the method which has been used to obtain the potential parameter spectrum, since it has to be inverted itself to obtain the energy spectrum, on the other hand, the inversion process is very hard to compute numerically the energy spectrum.
This paper aims to compute the energy spectrum which has not been computed in [8] [9].To the best of our knowledge, this paper is the first study to compute the eigenenergies for the hyperbolic single wave potential.Before presenting the results, test where made to find out whether the AIM is suitable for finding the energy spectrum for hyperbolic single wave potential., they were in a good agreement with them.The conclusion to be drawn from this behavior which taken from Table 1, is that for the application of the method, the number of iterations depends on the degree of convergence required.In order to get accurate result for higher excited states, one needs to increase the number of iterations j for ( ) r ϕ until the convergence for n E take place.Table 2 a comparison between AIM results and numerical results obtained through a private communication with the authors of [8] [9] are made for the potential (1) with 0 5 V = − , 1, 1 m λ = = , and different values of the parameter γ , and it is found that the results obtained by AIM are in good agreement with the other.In conclusion, AIM was used to obtain the energy spectrum hyperbolic single wave potential.One can note that the method gives the eigenvalues if one can transform the Schrödinger equation into a form of ( ) ( ) ( ) ( ) ( ) 0 0 y x k x y x z x y x ′′ ′ = + .The obtained numerical eigenenergies are in good agreement with the calculation results.
convenient change of variable, a straightforward calculation shows that Equation (17) becomes

Table 1 ,
elucidate the rate of convergence of the AIM for 0

Table 1 .
The rate of convergence of the AIM for the eigenenergies of the lowest states E 0 , E 1 with hyperbolic single wave potential in Equation (1) with V 0 = 5, λ = 0.2, γ = −0.4,and m = 1, computed for a number of iterations j.

Table 2 .
A comparison of the energy eigenvalues E n of the potential (1) with λ = 1, V 0 −5, m 1 and for different parameter γ.