Received 21 December 2015; accepted 19 January 2016; published 22 January 2016

ABSTRACT

Using the asymptotic iteration method, we obtain the S-wave solution for a short-range three-pa- rameter central potential with 1/r singularity and with a non-orbital barrier. To the best of our knowledge, this is the first attempt at calculating the energy spectrum for this potential, which was introduced by H. Bahlouli and A. D. Alhaidari and for which they obtained the “potential parameter spectrum”. Our results are also independently verified using a direct method of diagonalizing the Hamiltonian matrix in the J-matrix basis.

Keywords:

Short-Range Three-Parameter Central Potential, Asymptotic Iteration Method, Potential Parameter Spectrum Method, J-Matrix, Diagonalizing Method, Energy Eigenvalues

1. Introduction

In previous work [1] , we used the Asymptotic Iteration Method (AIM) to find the energy spectrum for the hyperbolic single wave potential, which was introduced by Bahlouli and Alhaidari [2] [3] . In the present work, we apply the same technique in [1] to the short-range three-parameter central potential, which was introduced by the same authors [3] [4] as

(1)

where is the potential strength and the range parameter l is positive with an inverse length units. The dimensionless parameter g is in the open range. This potential is short-range with 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 [5] with non-zero angular momentum (see Figure 1). However, 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 [3] [4] . In [4] , 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 [2] - [4] 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 [6] 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 [6] .

The calculation of energy eigenvalues is one of the basic problems of elementary quantum mechanics. Many techniques have been suggested to obtain the energy spectrum associated with a given potential. We will not try to give an overview here but one may consult [7] - [15] and references therein for the description of the AIM, which has been used successfully in the past. In this work, we will use this method to solve the time-independent S-wave (zero angular momentum) Schrödinger equation with potential (1) in order to find the energy spectrum.

The paper has the following structures. In Section 2, we briefly present an overview of the AIM and apply it to obtain the energy spectrum for potential (1). In Section 3, we find the eigen-energies where we make a comparison between the results obtained by the AIM and those obtained by the potential parameter spectrum method (PPSM). Additionally, we confirm our results independently by diagonalizing the Hamiltonian matrix, which is calculated in the J-matrix basis [16] [17] .

2. Basics of the AIM and Its Application to Potential (1)

In this section, we present the basic concepts of the AIM. For more details, we refer the reader to [7] - [15] . The method could be used to solve a second-order homogeneous linear differential equations of the form

(2)

where and are continuously differentiable functions over the defined interval of the coordinate x and the prime stands for the derivative with respect to x. According to the asymptotic aspect of the method and for sufficiently large n we have

(3)

where the function is independent of n and where

(4)

The general solution of Equation (2) is obtained as

(5)

The energy eigenvalues, E, of the Schrödinger equation corresponding to Equation (2) are obtained using the iteration terminating condition by means of Equation (3) that reads

(6)

If the problem is analytically solvable (exact solution) then. In this case, the AIM gives the eigenvalues and eigenfunctions in explicit algebraic form. However, in this case there are a limited number of potentials and the condition is satisfied at every point x in the defined interval. That is, is independent of x and the eigenvalues are determined from the n roots of the condition. In case of approximation, depends on both x and E. So, we have to determine a proper initial point value for solving with respect to E. Ideally, the solution should be independent of the choice of. However, an acceptable range of stability of the calculated energy eigenvalue may be found near the minimum value of the potential or the maximum value of the asymptotic wave function [14] [15] . Next, we apply the method to the problem with the potential function (1).

Inserting potential (1) in the time-independent radial Schrödinger’s equation results in the following second order linear differential equation

(7)

where is the angular momentum quantum number and we have adopted the atomic units. Defining the new variable, whose range is between −1 and +1, transforms this equation for S-wave () into the desired Equation (2) where we can apply the AIM with

(8)

(9)

Using these seed functions together with the quantization condition (6) in the iteration, we can calculate the energy eigenvalues. Since the problem is not exactly solvable, we have to select a proper. In this work, we observed that the best initial point is, which corresponds to the middle of the x interval. Therefore, at the end of the iterations when a stable result is reached, the energy spectrum is obtained by substituting.

3. Results and Discussion

The calculation of the energy eigenvalues for potential (1) is not as easy as that for the hyperbolic single wave potential, which we have studied in [1] . This is because the present potential has singularity at the origin. In this work, we take the value of the parameter g between zero and one. If g is less than zero or larger than one then the potential is not too interesting since it will not have the potential valley or potential hill (see Figure 1) and will only be similar to the well-known Hulthen or Yukawa potentials [18] . The necessary but not sufficient conditions for the existence of bound states for this potential are as follows:

1) If g is between zero and one, then any value of positive or negative could in principle support bound states.

2) If g is greater than one or less than zero, then the sign of must be the same as that of g.

In Table 1, we compare the energy spectrum obtained using the AIM outlined above and those obtained by the PPSM for, and for several values of g between zero and one. The Table shows good agreement between the two results. However, the PPSM results are in better agreement with the direct method of diagonalization of the Hamiltonian matrix in the J-matrix basis (HDM). In the Appendix, we show how this diagonalization procedure is carried out. Table 2 shows also another good agreement between the results obtained

by the two methods but now g is fixed while is varied.

Acknowledgements

The second author (ADH) would like to dedicate his contribution in this work to his friend, teacher and mentor Prof. Hashim A. Yamani on the occasion of his 70th birthday. Without his continued support and counsel, I could not have utilized his formulation of the J-matrix theory to the advancement of the tridiagonal representation approach, which was instrumental in the solution of the current and many other problems.

Cite this paper

Abdullah J.Sous,Abdulaziz D.Alhaidari, (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,04,79-85. doi: 10.4236/jamp.2016.41012

References

Appendix: Diagonalization of the Hamiltonian Matrix in the J-Matrix Basis

The Hamiltonian of the problem is, where the kinetic energy operator is and V(r) is the potential in (1). Diagonalization of the matrix representation of this Hamiltonian in a given basis is difficult for two reasons. First, it’s due to the singular behavior at the origin, which goes like for T and like for V. Second, it’s the long-range behavior of T, which dictates that we need to deal with infinite dimensional matrices. On the other hand, if we manage to transfer the singularity from V and add it to T by writing, where, then and we can write where the reference Hamiltonian. Consequently, the new potential U(r) is regular everywhere (see Figure A1)

and thus it could be approximated very well by its matrix elements in a finite subset of a square integrable basis. Then, what remains is to have a full account for the matrix representation of the reference Hamiltonian. However, is just the Coulomb Hamiltonian, which is known to have an exact tridiagonal matrix representation in the square integrable Laguerre basis with the following elements [16] [17]

(A1)

where is the associated Laguerre polynomial of order n in z and m is a length scale parameter. The normalization constant is chosen as. This gives an infinite tridiagonal symmetric matrix representation for the reference Hamiltonian whose elements read as follows

(A2)

As noted above, the potential U(r) is easily accounted for by its matrix elements in a finite subset of the basis (A1). Increasing the size of this subset will improve the accuracy of the results. These matrix elements are obtained by evaluating the integral

(A3)

where. The evaluation of this integral is performed using the Gauss quadrature approximation scheme associated with the Laguerre polynomials [19] , which gives

(A4)

for adequately large integer N. and are the N eigenvalues and corresponding normalized eigenvectors of the N ´ N tridiagonal basis overlap matrix whose elements are

(A5)

In an ideal situation where N is infinite the physical results should be independent of the choice of value of the numerical scale parameter m. Nonetheless, for finite calculation we should be able to find a range of values of m within which the results are stable and accurate to the desired number of significant digits. Increasing the size of the representation N will increase this range, which we refer to as the plateau of stability. Thus, this parameter is analogous to the parameter in the AIM noted in Section 2 above. For a given set of physical parameters, we add (A2) to (A4) giving the matrix elements of the total Hamiltonian which could be diagonalized numerically for a proper value of the numerical scale parameter m chosen from within the plateau of stability and for large enough matrix size N.

Table 1 and Table 2 show an excellent agreement of the results obtained by this diagonalization method (HDM) for with the PPSM and shows good to fair agreement with the AIM results. The agreement is less pronounced for the highest energy states. In the HDM, we used a basis size of N = 100.