﻿ Application of the Frobenius Method to the Schr？dinger Equation for a Spherically Symmetric Hyperbolic Potential

Open Access Library Journal
Vol.05 No.10(2018), Article ID:88136,15 pages
10.4236/oalib.1104950

Application of the Fröbenius Method to the Schrödinger Equation for a Spherically Symmetric Hyperbolic Potential

Hippolyte Nyengeri1, Rénovat Nizigiyima1, Eugène Ndenzako1, Félix Bigirimana2, Dieudonné Niyonkuru3, Audace Girukwishaka1

1Department of Physics, Faculty of Science, University of Burundi, Bujumbura, Burundi

2Petit Séminaire de Mika, Rutana, Burundi

3Faculty of Science and Technology, East African Leadership University, Muyinga, Burundi    Received: September 29, 2018; Accepted: October 26, 2018; Published: October 29, 2018

ABSTRACT

In this paper, an efficient technique for computing the bound state energies and wave functions of the Schrödinger Equation (SE) associated with a new class of spherically symmetric hyperbolic potentials is developed. This technique is based on a recent approximation scheme for the orbital centrifugal term and on the use of the Fröbenius method (FM). The bound state eigenvalues are given as zeros of calculable functions. The corresponding eigenfunctions can be obtained by substituting the calculated energies into the recurrence relations for the expanding coefficients of the Fröbenius series representing the solution. The excellent performance of this technique is illustrated through numerical results for some special cases like Pöschl-Teller potential (PTP), Manning-Rosen potential (MRP) and Pöschl-Teller polynomial potential (PTPP), with an application to the Gaussian potential well (GPW). Comparison with other methods is presented. Our results agree noticeably with the previously reported ones.

Subject Areas:

Quantum Mechanics

Keywords:

Fröbenius Method, Schrödinger Equation, Hyperbolic Potential, Bound States, Wave Functions 1. Introduction

It is well known that the exact solution of the Schrödinger Equation (SE) plays a vital role in quantum mechanics, and solving this equation is still an interesting work in the existing literature  -  . Generally, the SE is very difficult to solve for most physical central potentials    . It is for this reason that approximation and numerical methods are frequently used to arrive at the solution   -  .

Recently, the study of hyperbolic and exponential-type potentials has attracted a lot of interest by many authors due to their wide range of applications in physics    -    -  . However, the bound state solutions can only be obtained analytically for some of these potentials among which Pöschl-Teller     , Hulthen  , Manning-Rosen     and Rosen-Morse     potentials. Moreover, when, in three-dimensional space, the angular momentum quantum number is non zero, one can only solve the SE approximately using a suitable approximation scheme for the centrifugal term  . Such an approximation scheme includes the conventional approximation scheme proposed by Greene and Aldrich   ; the improved approximation scheme    suggested by Jia, Chen and Cui; the elegant approximation scheme   ; the approximation scheme proposed by Wei and Dong  and a recent approximation scheme suggested by Chen, Lou and You , and given by  :

$\frac{1}{{r}^{2}}\approx {\phi }^{2}\left[\frac{S}{{\mathrm{sinh}}^{2}\left(\phi r\right)}+\frac{t}{{\mathrm{cosh}}^{2}\left(\phi r\right)}\right]$ (1)

where $\phi$ , $S$ and $t$ are adjustable dimensionless parameters.

In this work, we use the approximation scheme (1) to deal with the centrifugal term and solve, by means of the Fröbenius method (FM), the SE with the class of hyperbolic potentials defined by

$\begin{array}{c}V\left(r\right)={C}_{-2}{\mathrm{coth}}^{2}\left(\alpha r\right)+{C}_{-1}\mathrm{coth}\left(\alpha r\right)+\sum _{i=1}^{3}{C}_{i}{\mathrm{tanh}}^{i}\left(\alpha r\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\sum _{j=1}^{K}{B}_{j}{\text{sech}}^{2j}\left(\alpha r\right)\end{array}$ (2)

where ${C}_{-2}$ , ${C}_{-1}$ , ${C}_{1}$ , ${C}_{2}$ , ${C}_{3}$ , ${B}_{1}$ , $\cdots$ , ${B}_{K}\in ℝ$ ; $\alpha \in {ℝ}_{0}$ and $K\in {ℕ}_{0}$ . This class includes many potentials, among which:

• The modified Pöschl-Teller potential (PTP)  which has been used to describe bending molecular vibrations    ;

• The PTP of the form 

${V}_{PT}\left(r\right)=\frac{-A\left(A+\alpha \right)}{{\mathrm{cosh}}^{2}\left(\alpha r\right)}+\frac{B\left(B-\alpha \right)}{{\mathrm{sinh}}^{2}\left(\alpha r\right)}$ (3)

This potential is an important diatomic molecular potential which has wide applications in physics and chemical physics  . It has been used to account for the physics of many systems which includes the excitons, quantum wires, and quantum dots  ;

• The Rosen-Morse potential (RMP)   and the Manning potential  which are generally used to explain vibrations of polyatomic molecular structures      ;

• The Manning-Rosen potential (MRP)   which is an important mathematical model for molecular vibrations and rotations   ;

• The Hartmann potential  . This potential exhibits in general an asymmetric double-well. It is therefore hoped that it can serve as a useful tool in the study of phenomena whose behavior is described by asymmetric double-wells  .

We have to add that the class of potentials under consideration contains all the approximation schemes of the orbital centrifugal term mentioned above. It also can be used to model the Gaussian potential well (GPW) through the expression 

$-\mathrm{exp}\left(-{r}^{2}\right)\approx \sum _{j=1}^{8}{B}_{j}{\text{sech}}^{2j}\left(\alpha r\right)$ (4)

where $\alpha$ , ${B}_{1}$ , $\cdots$ , ${B}_{8}$ are adjustable parameters.

Based on the foregoing, it can be said that the potential under consideration can be reduced to a lot of potentials which are known to play a very important role in many fields of physics such as molecular physics, solid state, and chemical physics. Consequently, it is for considerable interest to find general solution of its associated SE.

The rest of this paper is organized as follows. In Section 2, we apply the approximation (1) to the exact SE, and then convert the resulting equation into a dimensionless eigenvalue problem that can be solved by the use of the FM. In Section 3, we establish the Fröbenius series solution of this problem and show how to determine the bound state spectrum. Section 4 is devoted to three special cases, namely, the MRP, the PTP and the Pöschl-Teller polynomial potential (PTPP), with an application to the GPW in three dimensional space. In each case, the bound state energies obtained by means of the FM are presented and compared with those given by other methods. Finally, we present our conclusions in Section 5.

2. The Approximate Radial Schrödinger Equation

The time-independent SE for a particle of mass M that moves in three-dimensional space under the effect of a spherically symmetric potential is given by

$-\frac{{\hslash }^{2}}{2M}\left[{\Delta }_{r}+V\left(r\right)\right]\Psi \left(\stackrel{\to }{r}\right)=E\Psi \left(\stackrel{\to }{r}\right)$ (5)

Taking $\Psi \left(\stackrel{\to }{r}\right)={r}^{-1}\psi \left(r\right){Y}_{lm}\left(\theta ,\phi \right)$ and considering potential (2), we obtain the radial SE as

$\begin{array}{l}-\frac{{\hslash }^{2}}{2M}\frac{{\text{d}}^{2}\psi \left(r\right)}{\text{d}{r}^{2}}+\left[\frac{l\left(l+1\right)}{2M{r}^{2}}+{C}_{-2}{\mathrm{coth}}^{2}\left(\alpha r\right)+{C}_{-1}\mathrm{coth}\left(\alpha r\right)\\ +\sum _{i=1}^{3}{C}_{i}{\mathrm{tanh}}^{i}\left(\alpha r\right)+\sum _{j=1}^{K}{B}_{j}{\text{sech}}^{2j}\left(\alpha r\right)\right]\psi \left(r\right)=E\psi \left(r\right)\end{array}$ (6)

Using a new variable of the form

$R=br$ , $b\in \right]1,+\infty \left[$ (7)

we obtain the following differential equation:

$\begin{array}{l}-\frac{{\hslash }^{2}{b}^{2}}{2M}\frac{{\text{d}}^{2}\psi \left(R\right)}{\text{d}{R}^{2}}+\left[\frac{l\left(l+1\right){b}^{2}}{2M{R}^{2}}+{C}_{-2}{\mathrm{coth}}^{2}\left(\stackrel{˜}{\alpha }R\right)+{C}_{-1}\mathrm{coth}\left(\stackrel{˜}{\alpha }R\right)\\ +\sum _{i=1}^{3}{C}_{i}{\mathrm{tanh}}^{i}\left(\stackrel{˜}{\alpha }R\right)+\sum _{j=1}^{K}{B}_{j}{\text{sech}}^{2j}\left(\stackrel{˜}{\alpha }R\right)\right]\psi \left(R\right)=E\psi \left(R\right)\end{array}$ (8)

where $\stackrel{˜}{\alpha }=\alpha /b$ . Obviously, Equation (8) cannot be solved exactly when $l\ne 0$ due to the centrifugal term $1/{R}^{2}$ . To reduce this equation into a differential equation that can be solved by means of the FM, we must use a suitable approximation for the centrifugal term. We here consider the approximation scheme (1) proposed by Chen, Lu and You  . Substituting this approximation into Equation (8) and making the change of variable $\xi =\mathrm{tanh}\left(\stackrel{˜}{\alpha }R\right)$ , such that the domain $0 maps to $0<\xi <1$ , we find

$\begin{array}{l}{\xi }^{2}{\left({\xi }^{2}-1\right)}^{2}\frac{{\text{d}}^{2}\psi \left(\xi \right)}{\text{d}{\xi }^{2}}+2{\xi }^{2}\left({\xi }^{2}-1\right)\frac{\text{d}\psi \left(\xi \right)}{\text{d}\xi }\\ +\frac{2M}{{\hslash }^{2}{\alpha }^{2}}\left[-\left({C}_{-2}+\frac{l\left(l+1\right){\alpha }^{2}s}{2M}\right)-{C}_{-1}\xi +\left(E+\frac{l\left(l+1\right){\alpha }^{2}s}{2M}\right){\xi }^{2}\\ -{C}_{1}{\xi }^{3}-{C}_{2}{\xi }^{4}-{C}_{3}{\xi }^{5}-\sum _{j=1}^{K}{\stackrel{˜}{B}}_{j}{\left(1-{\xi }^{2}\right)}^{j}{\xi }^{2}\right]\psi \left(\xi \right)=0\end{array}$ (9)

where the ${\stackrel{˜}{B}}_{j}$ coefficients are such that

${\stackrel{˜}{B}}_{1}=\frac{l\left(l+1\right){\alpha }^{2}t}{2M}+{B}_{1}$ ; ${\stackrel{˜}{B}}_{j}={B}_{j}$ , $j=2,3,\cdots ,K$ (10)

3. Expansion around a Regular Singular Point

It is clear that $\xi =1$ and $\xi =0$ are both regular singular points for the differential Equation (9). The FM can therefore be applied with the wave function represented as generalized power series

$\psi \left(\xi \right)={\xi }^{\delta }\sum _{n=0}^{\infty }{a}_{n}{\xi }^{n}$ (11)

or

$\psi \left(\xi \right)={\left(\xi -1\right)}^{\delta }\sum _{n=0}^{\infty }{a}_{n}{\left(\xi -1\right)}^{n}$ (12)

where ${a}_{0}\ne 0$ . Let us seek $\psi \left(\xi \right)$ in the form (12). We first make the substitution $\eta =\xi -1$ . The resulting differential equation can be written as

$\begin{array}{l}\left(4{\eta }^{2}+12{\eta }^{3}+13{\eta }^{4}+6{\eta }^{5}+{\eta }^{6}\right)\frac{{\text{d}}^{2}\psi \left(\eta \right)}{\text{d}{\eta }^{2}}+\left(4\eta +14{\eta }^{2}+18{\eta }^{3}+10{\eta }^{4}+2{\eta }^{5}\right)\frac{\text{d}\psi \left(\eta \right)}{\text{d}\eta }\\ +\frac{2M}{{\hslash }^{2}{\alpha }^{2}}\left[-\left({C}_{-2}+\frac{l\left(l+1\right){\alpha }^{2}s}{2M}\right)-{C}_{-1}\left(1+\eta \right)+\left(E+\frac{l\left(l+1\right){\alpha }^{2}s}{2M}\right)\left(1+2\eta +{\eta }^{2}\right)\\ -{C}_{1}{\left(1+\eta \right)}^{3}-{C}_{2}{\left(1+\eta \right)}^{4}-{C}_{3}{\left(1+\eta \right)}^{5}-\sum _{j=1}^{K}{\stackrel{^}{B}}_{j}{\eta }^{j}{\left(1+\eta \right)}^{2}\right]\psi \left(\eta \right)=0\end{array}$

(13)

where the ${\stackrel{^}{B}}_{j}$ coefficients are calculated from the two following equations:

${\stackrel{^}{B}}_{2p+1}=\sum _{\begin{array}{c}j=p+1\\ j\le K\end{array}}^{2p+1}{\left(-1\right)}^{j}\left(\begin{array}{c}j\\ 2p+1-j\end{array}\right){2}^{2j-2p-1}{\stackrel{˜}{B}}_{j}$ ; $p=0,1,\cdots ,K-1$ (14)

${\stackrel{^}{B}}_{2p}=\sum _{\begin{array}{c}j=p\\ j\le K\end{array}}^{2p}{\left(-1\right)}^{j}\left(\begin{array}{c}j\\ 2p-j\end{array}\right){2}^{2j-2p}{\stackrel{˜}{B}}_{j}$ ; $p=1,2,\cdots ,K$ (15)

It can be verified that

$\sum _{p=1}^{2K}{\stackrel{^}{B}}_{p}{\eta }^{p}{\left(1+\eta \right)}^{2}=\sum _{p=1}^{2K+2}{\beta }_{p}{\eta }^{p}$ (16)

where

${\beta }_{1}={\stackrel{^}{B}}_{1}$ ; ${\beta }_{2}=2{\stackrel{^}{B}}_{1}+{\stackrel{^}{B}}_{2}$ ; ${\beta }_{j}={\stackrel{^}{B}}_{j-2}+2{\stackrel{^}{B}}_{j-1}+{\stackrel{^}{B}}_{j}$ , $j=3,4,\cdots ,2K$ (17)

${\beta }_{2K+1}={\stackrel{^}{B}}_{2K-1}+2{\stackrel{^}{B}}_{2K}$ ; ${\beta }_{2K+2}={\stackrel{^}{B}}_{2K}$ . (18)

Consequently, Equation (13) can be rewritten as

$\begin{array}{l}\left(4{\eta }^{2}+12{\eta }^{3}+13{\eta }^{4}+6{\eta }^{5}+{\eta }^{6}\right)\frac{{\text{d}}^{2}\psi \left(\eta \right)}{\text{d}{\eta }^{2}}\\ +\left(4\eta +14{\eta }^{2}+18{\eta }^{3}+10{\eta }^{4}+2{\eta }^{5}\right)\frac{\text{d}\psi \left(\eta \right)}{\text{d}\eta }\\ +\left[{d}_{0}+\epsilon +\left({d}_{1}+2\epsilon \right)\eta +\left({d}_{2}+\epsilon \right){\eta }^{2}+\sum _{i=3}^{5}{d}_{i}{\eta }^{i}+\sum _{p=6}^{2K+2}{b}_{p}{\eta }^{p}\right]\psi \left(\eta \right)=0\end{array}$ (19)

where

$\epsilon =\frac{2M}{{\hslash }^{2}{\alpha }^{2}}E$ ; ${d}_{0}=-\frac{2M}{{\hslash }^{2}{\alpha }^{2}}\left[{C}_{-2}+{C}_{-1}+{C}_{1}+{C}_{2}+{C}_{3}\right]$ (20)

${d}_{1}=\frac{2M}{{\hslash }^{2}{\alpha }^{2}}\left[\frac{2l\left(l+1\right){\alpha }^{2}s}{2M}-{C}_{-1}-3{C}_{1}-4{C}_{2}-5{C}_{3}-{\beta }_{1}\right]$ (21)

${d}_{2}=\frac{2M}{{\hslash }^{2}{\alpha }^{2}}\left[\frac{l\left(l+1\right){\alpha }^{2}s}{2M}-3{C}_{1}-6{C}_{2}-10{C}_{3}-{\beta }_{2}\right]$ (22)

${d}_{3}=-\frac{2M}{{\hslash }^{2}{\alpha }^{2}}\left[{C}_{1}+4{C}_{2}+10{C}_{3}+{\beta }_{3}\right]$ ; ${d}_{4}=-\frac{2M}{{\hslash }^{2}{\alpha }^{2}}\left[{C}_{2}+5{C}_{3}+{\beta }_{4}\right]$ (23)

${d}_{5}=-\frac{2M}{{\hslash }^{2}{\alpha }^{2}}\left({C}_{3}+{\beta }_{5}\right)$ ; ${b}_{j}=-\frac{2M}{{\hslash }^{2}{\alpha }^{2}}{\beta }_{j}$ , $j=6,7,\cdots ,2K+2$ (24)

Note that Equation (19) can be solved near $\eta =0$ by writing

$\psi \left(\eta \right)={\eta }^{\delta }\sum _{n=0}^{\infty }{a}_{n}{\eta }^{n}$ (25)

Substitution of Equation (25) into (19) leads to the recurrence relation

$\begin{array}{l}\left[4+{d}_{0}+\epsilon +4\delta \left(2+\delta \right)+8n\left(1+\delta \right)+4{n}^{2}\right]{a}_{n+1}+\left[{d}_{1}+2\epsilon +2\delta \left(1+6\delta \right)\\ +2n\left(1+12\delta \right)+12{n}^{2}\right]{a}_{n}+\left[8+{d}_{2}+\epsilon -\left(21-13\delta \right)\delta -\left(21-26\delta \right)n\\ +13{n}^{2}\right]{a}_{n-1}+\left[16+{d}_{3}-2\delta \left(10-3\delta \right)-4n\left(5-3\delta \right)+6{n}^{2}\right]{a}_{n-2}\\ +\left[6+{d}_{4}-\left(5-\delta \right)\delta -\left(5-2\delta \right)n+{n}^{2}\right]{a}_{n-3}+{d}_{5}{a}_{n-4}+\sum _{p=6}^{2K+2}{b}_{p}{a}_{n-p+1}=0\end{array}$ (26)

for the expanding coefficients, with the understanding that ${a}_{i}=0$ for $i<0$ . Setting $n=-1$ in the above relation, we obtain the indicial equation $\left(4{\delta }^{2}+{d}_{0}+\epsilon \right){a}_{0}=0$ which is solved by ${\delta }_{1}=-\frac{1}{2}\sqrt{-{d}_{0}-\epsilon }$ and ${\delta }_{2}=\frac{1}{2}\sqrt{-{a}_{0}-\epsilon }$ . We can therefore conclude that the two solutions of the radial SE obtained as generalized series, one with $\delta ={\delta }_{1}$ and the other with $\delta ={\delta }_{2}$ , are linearly independent. The value of $\delta$ determines the behavior of $\psi \left(\eta \right)$ for $\eta \to -1$ ( $r\to 0$ ) and $\eta \to 0$ ( $r\to \infty$ ). Clearly, only $\delta >0$ is acceptable, since only in this case $\psi \left(r=0\right)=0$ . Such a solution to the radial SE contains only the series with $\delta ={\delta }_{2}$ , and in the following will be denoted by

$\psi \left(\eta ,\epsilon \right)={\eta }^{{\delta }_{2}}\sum _{n=0}^{\infty }{a}_{n}{\eta }^{n}$ (27)

where the dependence on the energy eigenvalue is explicitly marked. The function $\psi \left(\eta ,\epsilon \right)$ has an important property that $\psi \left(\eta =0,\epsilon \right)=0$ , which can be regarded as a boundary condition for the radial equation at $r=\infty$ . If $\psi \left(\eta ,\epsilon \right)$ is associated with a bound state, it has another important property, i.e., $\psi \left(\eta =-1,\epsilon \right)=0$ , which can be considered as the boundary condition for bound state wave function at $r=0$ . It means that for arbitrary $\epsilon$ , the wave function $\psi \left(r,\epsilon \right)$ obtained from (27) is not square integrable. If, for a particular value of $\epsilon$ , $\psi \left(\eta =-1,\epsilon \right)=0$ , the wave function is square integrable. In other words, this condition determines the energies of the discrete spectrum. We have to emphasize that in practice, the function $\psi \left(\eta ,\epsilon \right)$ has to be approximated by truncating the series in Equation (27) at suitably high order N, which requires high precision computing since some of the expanding coefficients ${a}_{i}$ in Equation (27) are extremely big for high values of i, especially when the ${C}_{i}$ and/or the ${B}_{j}$ parameters are significantly large, and ${a}_{0}$ is not sufficiently small. To deal with this problem, one can appeal to a new software package for arbitrary precision computation, named MPFUN 2015 and developed by David H. Bailey  . In our calculations, the precision level is set to 60 digits and ${a}_{0}$ equals 10−9. The number of terms in the series (27) is 501. We have chosen fortran 90 as the programming language. One has to add that the truncated series is a function of $\epsilon$ whose zeros correspond to bound state values of $\epsilon$ .

Before finishing this section, we have to indicate that when the class of potentials (2) is subject to the condition that ${C}_{-1}={C}_{1}={C}_{3}=0$ , we found it most appropriate to make the change of variable $\xi ={\text{sech}}^{2}\left(\stackrel{˜}{\alpha }R\right)$ instead of using the transformation $\xi =\mathrm{tanh}\left(\stackrel{˜}{\alpha }R\right)$ . The radial SE can in this case be modified as

$\begin{array}{l}{\xi }^{2}\left(1-\xi \right)\frac{{\text{d}}^{2}\psi \left(\xi \right)}{\text{d}{\xi }^{2}}+\xi \left(1-\frac{3}{2}\xi \right)\frac{\text{d}\psi \left(\xi \right)}{\text{d}\xi }\\ +\frac{1}{4}\left\{\epsilon -\left[{c}_{-2}+\frac{{\alpha }^{2}}{{\hslash }^{2}}l\left(l+1\right)s\right]\frac{1}{1-\xi }\\ -\left[{c}_{2}-l\left(l+1\right)\frac{{\alpha }^{2}}{{\hslash }^{2}}s\right]-\sum _{j=1}^{K}{\stackrel{˜}{b}}_{j}{\xi }^{j}\right\}\psi \left(\xi \right)=0\end{array}$ (28)

where

$\epsilon =\frac{2M}{{\hslash }^{2}}E$ , ${c}_{-2}=\frac{2M}{{\hslash }^{2}}{C}_{-2}$ , ${c}_{2}=\frac{2M}{{\hslash }^{2}}{C}_{2}$ , ${\stackrel{˜}{b}}_{j}=\frac{2M}{{\hslash }^{2}}{\stackrel{˜}{B}}_{j}$ (29)

By putting the series (12) and its two first derivatives into Equation (28), and then changing variables to $\eta =\xi -1$ , we obtain the recurrence relation

$\begin{array}{l}\left[\left(n+\delta +1\right)\left(n+\delta +\frac{1}{2}\right)-\frac{{c}_{-2}}{4}-\frac{{\alpha }^{2}}{4{\hslash }^{2}}l\left(l+1\right)s\right]{a}_{n+1}\\ +\left[2{\left(n+\delta \right)}^{2}-\frac{\epsilon -{c}_{2}-{b}_{0}}{4}-\frac{{\alpha }^{2}}{4{\hslash }^{2}}l\left(l+1\right)s\right]{a}_{n}\\ +\left[\left(n+\delta -1\right)\left(n+\delta -\frac{1}{2}\right)+\frac{{b}_{1}}{4}\right]{a}_{n-1}+\frac{{b}_{2}}{4}{a}_{n-2}+\cdots +\frac{{b}_{K}}{4}{a}_{n-K}=0\end{array}$ (30)

for the expanding coefficients. The ${b}_{j}$ coefficients are such that

${b}_{0}=\sum _{j=1}^{K}{\stackrel{˜}{b}}_{j}$ ; ${b}_{j}=\sum _{i=j}^{K}\left(\begin{array}{c}i\\ j\end{array}\right){\stackrel{˜}{b}}_{i}$ , $j=1,2,\cdots ,K$ (31)

The indicial equation, namely ${\delta }^{2}-\frac{1}{2}\delta -\frac{{a}_{-2}}{4}-\frac{{\alpha }^{2}}{4{\hslash }^{2}}l\left(l+1\right)s=0$ has two solutions, i.e.,

${\delta }_{1}=\frac{1/2-\sqrt{1/4+{c}_{-2}-\frac{{\alpha }^{2}}{{\hslash }^{2}}l\left(l+1\right)s}}{2}$ , ${\delta }_{2}=\frac{1/2+\sqrt{1/4+{c}_{-2}-\frac{{\alpha }^{2}}{{\hslash }^{2}}l\left(l+1\right)s}}{2}$ (32)

When $\delta >0$ , the boundary conditions for bound wave functions are satisfied.

4. Numerical Results for Some Special Cases

To test the accuracy of our approximate results, we consider three special cases, namely the case of the MRP   , that of the Second Pöschl-Teller potential (SPTP)  and the case where $V\left(r\right)=\sum _{j=1}^{K}{B}_{j}{\text{sech}}^{2}\left(\alpha r\right)$ , with an application to the Gaussian potential well (GPW)   . For each case, we give a numerical comparison of the energy eigenvalues obtained by using the FM with those obtained with the use of other methods.

4.1. The Manning-Rosen Potential

The MRP is defined by     :

$V\left(r\right)=\frac{{\hslash }^{2}}{2M{b}^{2}}\left[\frac{a\left(a-1\right)}{{\left({\text{e}}^{r/b}-1\right)}^{2}}-\frac{A}{{\text{e}}^{r/b}-1}\right]=\frac{{\hslash }^{2}}{2M{b}^{2}}\left[\frac{a\left(a-1\right){\text{e}}^{-2r/b}}{{\left(1-{\text{e}}^{-r/b}\right)}^{2}}-\frac{A{\text{e}}^{-r/b}}{1-{\text{e}}^{-r/b}}\right]$ (33)

Here, $a$ and the strength parameter A are two dimensionless parameters, while the screening parameter b, having a dimension of length, is related to the range of potential. It can be verified that the class of potentials (2) reduces to MRP by choosing

$K=1$ , ${C}_{-2}=a\left(a-1\right){\hslash }^{2}/\left(8M{b}^{2}\right)$ , ${C}_{-1}=-\left[a\left(a-1\right)+A\right]{\hslash }^{2}/\left(4M{b}^{2}\right)$ ,

${C}_{1}=0$ , ${C}_{2}={B}_{1}=\left[a\left(a-1\right)+2A\right]{\hslash }^{2}/\left(8M{b}^{2}\right)$ and $\alpha =1/\left(2b\right)$ .

The MRP is used as an important mathematical model for molecular vibrations and rotations   . It has considerable applications in several bound state and scattering problems in physics. For $a=0$ or 1, this potential reduces to the familiar short-range Hulthen potential  having useful applications in nuclear, particle, solid-state and chemical physics. There is a relative minimum value at ${r}_{0}=b\mathrm{ln}\left[1+2a\left(a-1\right)/A\right]$ with value $V\left({r}_{0}\right)=-{\hslash }^{2}{A}^{2}/\left[8M{b}^{2}a\left(a-1\right)\right]$ for $a<0$ or $a>1$ and $A>0$ .

Let us indicate that approximate analytical solutions of $l$ -wave SE for the MRP have been presented in  by a proper approximation of the centrifugal term. Furthermore, a purely numerical integrating procedure has been programmed  as well for bound states, invoking the Mathematica package, which offers decent results, especially in the short potential well, i.e., small $l$ and $a$ .

In Table 1, we report our calculated energies for selected $l\ne 0$ states having principal quantum number $n\le 6$ , two different values of parameter $a$ , $A=2b$ and three different values of the screening parameter b. The parameter $b$ has been chosen such that $\stackrel{˜}{\alpha }=0.0001$ , which implies that $S=0.9999999999999999864889$ and $t=0.3333422254140179276$ . Our results are compared to both the previous ones  and the exact results obtained with the aid of the Mathematica package program by Lucha and Schöberl  as shown in Table 1. It is obvious that the present results are in good agreement with the accurate numerical ones and those calculated by Qiang and Dong  .

4.2. The Second Pöschl-Teller Potential (SPTP)

We here consider the PTP given by 

$V\left(r\right)=-\frac{{\hslash }^{2}{\alpha }^{2}}{2M}\left[\frac{\lambda \left(\lambda +1\right)}{{\mathrm{cosh}}^{2}\left(\alpha r\right)}-\frac{\kappa \left(\kappa -1\right)}{{\mathrm{sinh}}^{2}\left(\alpha r\right)}\right]$ (34)

where the parameters $\lambda$ and $\kappa$ are relevant to the potential depth, while $\alpha$ describes the width of the potential. The SPTP can be obtained from the class of potentials (2) by choosing $K=1$ and

${C}_{-2}=\frac{{\hslash }^{2}{\alpha }^{2}}{2M}\kappa \left(\kappa -1\right)=-{C}_{2}$ , ${C}_{-1}={C}_{1}=0$ , ${B}_{1}=-\frac{{\hslash }^{2}{\alpha }^{2}}{2M}\left[\lambda \left(\lambda +1\right)+\kappa \left(\kappa -1\right)\right]$

We show, in Table 2, the energy levels calculated by the FM, those obtained numerically by means of the Mathematica package developed by Lucha and Schöberl  and the approximate analytical ones  , respectively for a number of given principal quantum number n and angular quantum number $l$ with some typical values of parameter $\alpha$ . It is found that our results agree noticeably with numerically exact ones. Some of them are even better than those obtained by Chen and Dong  .

Table 1. Eigenvalues as a function of b for 2p-6g states and for $a=0.75$ and $a=1.5$ ; A = 2b, $S=0.9999999999999999864889$ and $t=0.3333422254140179276$ .

Table 2. Energy eigenvalues in atomic units ( $\hslash =M=1$ ) for $\lambda =20$ and $\kappa =4$ .

4.3. The Pöschl-Teller Polynomial Potential

As another interesting special case of the class of potentials under consideration, we consider the potential of the form

$V\left(r\right)=\sum _{j=1}^{K}{B}_{j}{\text{sech}}^{2j}\left(\alpha r\right)$ (35)

which we call “Pöschl-Teller polynomial potential” since it is a polynomial function in the PTP of the form ${\text{sech}}^{2}\left(\alpha r\right)$ . The PTPP can be obtained from Equation (2) by choosing ${C}_{-2}={C}_{-1}={C}_{1}={C}_{2}={C}_{3}=0$ . It has a very interesting aspect that it can be used to model the Gaussian function ${f}_{G}\left(r\right)=-\mathrm{exp}\left(-r\right)$ by suitably choosing the values of $\alpha$ , K and the ${B}_{j}$ coefficients  .

In order to show that the potential (35) can have useful applications in physics, we study the three-dimensional SE for the attractive Gaussian potential defined by

${V}_{G}\left(r\right)=-\gamma \mathrm{exp}\left(-\lambda {r}^{2}\right)$ (36)

where $\gamma >0$ is the well depth and $\lambda >0$ determines its width. Note that this potential has been introduced to model nucleon-nucleon scattering  and has played an important role in nuclear physics. In semiconductor physics, the potential (36) has been used as a model potential to investigate the electronic structure of the quantum dots with single or more electrons   , impurity  and excitons  .

The radial SE associated with the attractive Gaussian potential (36) can be written as

$-\frac{\lambda {\hslash }^{2}}{2M}\frac{{\text{d}}^{2}\psi \left(r\right)}{\text{d}{r}^{2}}+\left[\frac{l\left(l+1\right)\lambda }{2M{r}^{2}}-\gamma \mathrm{exp}\left(-{r}^{2}\right)\right]\psi \left(r\right)=E\psi \left(r\right)$ (37)

where $r=\sqrt{\lambda }r$ . By using the approximation scheme (4) in Equation (37), we obtain the following differential equation:

$-\frac{\lambda {\hslash }^{2}{b}^{2}}{2M}\frac{{\text{d}}^{2}\psi \left(R\right)}{\text{d}{R}^{2}}+\left[\frac{l\left(l+1\right)\lambda {b}^{2}}{2M{R}^{2}}+\sum _{j=1}^{8}{B}_{j}{\text{sech}}^{2j}\left(\stackrel{˜}{\alpha }R\right)\right]\psi \left(R\right)=E\psi \left(R\right)$ (38)

Here, $\stackrel{˜}{\alpha }=\alpha /b$ , $R=br$ and ${B}_{j}=\gamma {B}_{j}$ , $j=1,\cdots ,8$ . Using the approximation scheme (1) for the centrifugal term $1/{R}^{2}$ , and changing variables from R to $\xi ={\text{sech}}^{2}\left(\stackrel{˜}{\alpha }R\right)$ , we obtain the dimensionless SE (28) with $\epsilon =2ME/\left(\lambda {\hslash }^{2}{\alpha }^{2}\right)$ , ${c}_{-2}={c}_{2}=0$ , ${\stackrel{˜}{b}}_{j}=2M{\stackrel{˜}{B}}_{j}/\left(\lambda {\hslash }^{2}{\alpha }^{2}\right)$ .

Figure 1 shows the variation of potential as a function of the radial position (dashed curve) and the change in discrete electron energy levels ( ${E}_{N,l}$ ) for different values of $N$ and $l$ when $\lambda =1/{a}_{B}^{2}$ and $\gamma =400$ Ryd. Here ${a}_{B}$ and Ryd are the Bohr radius and Rydberg energy, respectively. Note that $N=0,1,2,\cdots$ and $l=0,1,2,\cdots$ are radial and angular quantum numbers, respectively. Our results are compared with those obtained by means of the Köksal’s empirical formula  and the numerically calculated ones by the exact Hamiltonian diagonalization on a finite basis of Coulomb sturmian functions defined by

Figure 1. Bound state energy values of attractive Gaussian potential for different radial and angular quantum numbers.

${S}_{n,l}^{\kappa }\left(r\right)={N}_{n,l}^{\kappa }{r}^{l+1}{\text{e}}^{-\kappa r}{L}_{n-l-1}^{2l+1}\left(2\kappa r\right)$ (39)

where ${L}_{m}^{\alpha }\left(x\right)$ denotes the associated Laguerre polynomial and n the principal quantum number. The normalization constant ${N}_{n,l}^{\kappa }$ , given by

${N}_{n,l}^{\kappa }=\sqrt{\frac{\kappa }{n}}{\left(2\kappa \right)}^{l+1}{\left[\frac{\left(n-l-1\right)!}{\left(n+l\right)!}\right]}^{1/2}$ (40)

is obtained from the normalization condition ${\int }_{0}^{\infty }{\left[{S}_{n,l}^{\kappa }\left(r\right)\right]}^{\ast }{S}_{n,l}^{\kappa }\left(r\right)\text{d}r=1$ . We have chosen ${N}_{s}=500$ as the number of Coulomb sturmian functions and 0.75 as the value of $\kappa$ . For the case of $N=0$ , we also consider results calculated by using a simple analytical expression provided by the variational method  . It is obvious that our results $\left({E}_{N,l}\left(FM\right)\right)$ are much better than those obtained using the Köksal’s empirical formula $\left({E}_{N,l}\left(Koksal\right)\right)$  . They are in good agreement with the numerically calculated ones $\left({E}_{N,l}\left(Diag\right)\right)$ . In the case of $N=0$ , our results are very closer to both the numerically obtained ones, the Köksal’s ones and the variational ones $\left({E}_{N,l}\left(F&G\right)\right)$ . We appreciate that the present results appears to approach the numerical eigenvalues reasonably well even when $N>0$ and $l$ is not small, while the accuracy of the empirical formula decreases noticeably with $l$ and most remarkably with $N$  .

5. Conclusions

In this work, we have applied a recent approximation scheme to the centrifugal term ${r}^{-2}$ to solve the SE with a new class of spherically symmetric hyperbolic potentials in the framework of the FM. First, we have shown that our class of potentials can be reduced to a lot of hyperbolic or exponential potentials which are known to play a very important role in many fields of physics. Next, we have proved that the application of the FM to the family of potentials in question allows an easy determination of the bound state spectrum provided that the centrifugal potential is adequately approximated. Finally, we have calculated the bound state energy eigenvalues for three special cases, namely the MRP, the SPTP and the PTPP, with an application to the GPW. For each case, it is found that our results are in good agreement with those obtained by other methods for short-range potentials, i.e., small $l$ and $\alpha$ .

One has to add that our approach can also be used to determine the approximate bound state wave functions. After substituting the calculated bound state energies to the recursion relation (26) or (30), the coefficients of the generalized series (25) can be successively determined in order to obtain the unnormalized wave functions.

Acknowledgements

Nyengeri Hippolyte would like to thank the Académie de Recherche et d’Enseignement Supérieur-Commission de la Coopération au Développement (ARES-CCD) for having provided him with a desktop and a printer which have greatly facilitated the achievement of this work.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Nyengeri, H., Nizigiyima, R., Ndenzako, E., Bigirimana, F., Niyonkuru, D. and Girukwishaka, A. (2018) Application of the Fröbenius Method to the Schrödinger Equation for a Spherically Symmetric Hyperbolic Potential. Open Access Library Journal, 5: e4950. https://doi.org/10.4236/oalib.1104950

References

1. 1. Dong, S.H. and Garcia-Ravalo, J. (2007) Exact Solutions of the s-Wave Schrodinger Equation with Manning-Rosen Potential. Physica Scripta, 75, 307-309.
https://doi.org/10.1088/0031-8949/75/3/013

2. 2. Qiang, W.C. and Dong, S.H. (2009) The Manning-Rosen Potential Studied by a New Approximation Scheme to the Centrifugal Term. Physica Scripta, 79, Article ID: 045004.
https://doi.org/10.1088/0031-8949/79/04/045004

3. 3. Ikot, A.N., Antia, A.D., Akpan, I.O. and Awoga, O.A. (2013) Bound State Solutions of Schrodinger Equation with Modified Hylleraas plus Exponential Rosen Morse Potential. Revista Mexicana de Fisica, 59, 46-53.

4. 4. Dong, S., Garcia-Ravelo, J. and Dong, S.H. (2007) Analytical Approximations to the l-Wave Solutions of the Schrodinger Equation with an Exponential-Type Potential. Physica Scripta, 76, 393-396.
https://doi.org/10.1088/0031-8949/76/4/019

5. 5. Wen, F.K., Yang, Z.Y., Chong, L., Yang, W.L. and Zhang, Y.Z. (2014) Exact Polynomial Solutions of the Schrodinger Equation with Various Hyperbolic Potentials. Communications in Theoretical Physics, 61, 153-159.
https://doi.org/10.1088/0253-6102/61/2/02

6. 6. Nyengeri, H., Simbizi, R., Girukwishaka, A., Nizigiyimana, R. and Ndenzako, E. (2018) Frobenius Series Solutions of the Schrodinger Equation with Various Types of Symmetric Hyperbolic Potentials in One Dimension. Open Access Library Journal, 5, 1-14.

7. 7. Roy, A.K. (2014) Studies of Bound States Spectra of Manning-Rosen Potential. Modern Physics Letters A, 29, Article ID: 1450042.
https://doi.org/10.1142/S0217732314500424

8. 8. Ikhdair, S.M. (2011) On the Bound-State Solutions on the Manning-Rosen Potential Including Improved Approximation to the Orbital Centrifugal Term. Physica Scripta, 83, Article ID: 015010.
https://doi.org/10.1088/0031-8949/83/01/015010

9. 9. Chen, C.Y., Lu, F.L. and You, Y. (2012) Scattering States of Modified P?schl-Teller potential in D-Dimension. Chinese Physics B, 21, Article ID: 030302.
https://doi.org/10.1088/1674-1056/21/3/030302

10. 10. Hammed, R.H. (2012) Approximate Solution of Schrodinger Equation with Manning-Rosen Potential in Two Dimensions by Using the shifted 1/N Expansion Method. Journal of Basrah Researches (Sci-ences), 38, 51-59.

11. 11. Chaudhuri, R.N. and Mondal, M. (1995) Eigenvalues of Anharmonic Oscillators and the Perturbed Coulomb Problem in N-Dimensional Space. Physical Review A, 52, Article ID: 1850.
https://doi.org/10.1103/PhysRevA.52.1850

12. 12. Miao, C. H. and Qian, S. W. (1997) Variational Sypersymmetric WKB Ap-proximation. Physical Review A, 56, Article ID: 2412.
https://doi.org/10.1103/PhysRevA.56.2412

13. 13. Green, R.L. and Aldrich, C. (1976) Variational Wave Functions for a Screened Coulomb Potential. Physical Review A, 14, Article ID: 2363.
https://doi.org/10.1103/PhysRevA.14.2363

14. 14. Tang, A.Z. and Chan, F.T. (1987) Shifted 1/N Expansion for the Hulthen Potential. Physical Review A, 35, Article ID: 911.
https://doi.org/10.1103/PhysRevA.35.911

15. 15. Dutt, R., Mukherji, U. and Varshni, Y.P. (1986) Shifted Large-N Expansion for the Bound States of the Hellmann Potential. Physical Review A, 34, Article ID: 777.
https://doi.org/10.1103/PhysRevA.34.777

16. 16. Nioto, M.M. (1978) Exact Wave-Function Normalization Constants for the B0tanhz-U0coshz and Poschl-Teller Potentials. Physical Review A, 17, 1273-1284.
https://doi.org/10.1103/PhysRevA.17.1273

17. 17. Onate, C.A. (2015) Approximate Solutions of the Non-Relativistic Schrodinger Equation with Poschl-Teller Potential. Chinese Journal of Physics, 53, Article ID: 060002.

18. 18. Ikot, A.N. (2010) Approximate Solution of the Schrodinger Equation with Rosen-Morse Potential Including the Centrifugal Term. Applied Physics Research, 2, 202-208.
https://doi.org/10.5539/apr.v2n2p202

19. 19. Downing, C.A. (2013) On a Solution of the Schrodinger Equation with Hyperbolic Double-Well Potential. Journal of Mathematical Physics, 54, Article ID: 072101.
https://doi.org/10.1063/1.4811855

20. 20. Hartman, R.R. (2014) Bound States in Hyperbolic Asymmetric Double-Well. Journal of Mathematical Physics, 55, Article ID: 012105.
https://doi.org/10.1063/1.4861938

21. 21. Rosen, N. and Morse, P.M. (1932) On the Vibration of Polyatomic Molecules. Physical Review, 42, 210-217.
https://doi.org/10.1103/PhysRev.42.210

22. 22. Zlatev, S.I. (2018) Pekeris-Type Approximation for the l-Wave in P?schl-Teller Potential. Unpublished.
https://arxiv.org/pdf/1311.5794v1.pdf

23. 23. Jia, C.S., Chen, T. and Cui, L.G. (2009) Approximate Analytical Solutions of the Dirac Equation with the Generalized Poschl-Teller Potential Including the Pseudo-Centrifugal Term. Physics Letters A, 373, 1621-1626.
https://doi.org/10.1016/j.physleta.2009.03.006

24. 24. Ourdi, R., Hassanabadi, S., Rajabi, A.A. and Hasanabadi, H. (2012) Approximate Bound State Solutions of DKP Equation for Any J State in the Presence of Woods-Saxon Potential. Communications in Theoretical Physics, 57, 15-18.
https://doi.org/10.1088/0253-6102/57/1/04

25. 25. Wei, G.F. and Dong, S.H. (2008) Approximately Analytical Solutions of the Manning-Rosen Potential with the Spin-Orbit Coupling Term and Spin Symmetry. Physics Letters A, 373, 49-57.
https://doi.org/10.1016/j.physleta.2008.10.064

26. 26. Flüge, S. (1999) Practical Quantum Mechanics. Springer-Verlag, Berlin.

27. 27. Carrington, T. (1990) Perturbation Theory for Bending Potentials. Molecular Physics, 70, 757-766.
https://doi.org/10.1080/00268979000101331

28. 28. Iachello, F. and Oss, S. (1993) Algebraic Model of Bending Vibrations of Complex Molecules. Chemistry Letters, 205, 285-289.
https://doi.org/10.1016/0009-2614(93)89244-C

29. 29. Iachello, F. and Oss, S. (1993) Vibrational Spectroscopy and Intramolecular Relaxation of Benzene. The Journal of Chemical Physics, 99, 7337-7349.
https://doi.org/10.1063/1.465715

30. 30. Wei, G.F. and Dong, S.H. (2010) A Novel Algebraic Approach to Spin Symmetry for Dirac Equation with Scalar and Vector Second Poschl-Teller Potential. The European Physical Journal A, 43, 185-190.
https://doi.org/10.1140/epja/i2009-10901-8

31. 31. Hasan, Y. and Tomak, M. (2005) Nonlinear Optical Properties of a Poschl-Teller Quantum Well. Physical Review B, 72, Article ID: 115340.
https://doi.org/10.1103/PhysRevB.72.115340

32. 32. Manning, M.F. (1935) Energy Levels of a Symmetrical Double Minima Problem with Applications to the NH3 and ND3 Moleculs. The Journal of Chemical Physics, 3, 136-138.
https://doi.org/10.1063/1.1749619

33. 33. Costain, C.C. and Dowling, J.M. (1960) Microwave Spectrum and Molecular Structure of Formamide. The Journal of Chemical Physics, 32, 158-165.
https://doi.org/10.1063/1.1700891

34. 34. Fletcher, W.H. and Brown, F.B. (1963) Vibrational Spectra and the Inversion Phenomenon in Cyanamide and Deuterated Cyanamide. The Journal of Chemical Physics, 39, 2478-2490.
https://doi.org/10.1063/1.1734051

35. 35. Hilmi, Y., Oktay, A. and Mustafa, S. (2016) Modeling of Diatomic Molecules. Molecular Physics, 114, 3134-3142.
https://doi.org/10.1080/00268976.2016.1220645

36. 36. Manning, F.M. and Rosen, N. (1933) A Potential Function for the Vi-brations of Diatomic Molecules. Physical Review, 44, 953.

37. 37. LeRoy, R.J. and Bernstein, R.B. (1970) Dissociation Energy and Long-Range Potential of Diatomic Molecules from Vibrational Spacings of Higher Levels. The Journal of Chemical Physics, 52, 3869-3879.
https://doi.org/10.1063/1.1673585

38. 38. Cai, J.M., Cai, P.Y. and Inomata, A. (1986) Path-Integral Treatment of the Hulthen Potential. Physical Review A, 34, 4621-4628.
https://doi.org/10.1103/PhysRevA.34.4621

39. 39. Bailey, D.H. (2015) A Thread-Safe Arbitrary Precision Computation Package.
http://www.davidhbailey.com/dhbpapers/mpfun2015.pdf

40. 40. You, Y., Lu, F.L., Sun, D.S., Chen, C.Y. and Dong, S.H. (2013) Solutions of the Second Poschl-Teller Potential Solved by an Improved Scheme to the Centrifugal Term. Few-Body Systems, 54, 2125-2132.
https://doi.org/10.1007/s00601-013-0725-y

41. 41. Koksal, K. (2012) A Simple Analytical Expression for Bound State Energies for an Attractive Gaussian Confining Potential. Physica Scripta, 86, Article ID: 035006.
https://doi.org/10.1088/0031-8949/86/03/035006

42. 42. Fenandez, F.M. and Garcia, J. (2013) Comment on “A Simple Analytical Expression for Bound State Energies for an Attractive Gaussian Confining Potential”. Physica Scripta, 87, Article ID: 027001.
https://doi.org/10.1088/0031-8949/87/02/027001

43. 43. Lucha, W. and Schoberl, F.F. (1999) Solving the Schrodinger Equation for Bound States with Mathematica 3.0. International Journal of Modern Physics C, 10, 607-619.
https://doi.org/10.1142/S0129183199000450

44. 44. Buck, B., Friedrich, H. and Wheatley, C. (1977) Local Potential Models for the Scattering of Complex Nuclei. Nuclear Physics A, 275, 246-268.
https://doi.org/10.1016/0375-9474(77)90287-1

45. 45. Gomez, S. and Romero, R. (2009) Few-Electron Semiconductor Quantum Dots with Gaussian Confinement. Central European Journal of Physics, 7, 12-21.

46. 46. Adamowski, J., Sobkowiez, M., Szafran, B. and Bednarek, S. (2000) Electron Pair in a Gaussian Confining Potential. Physical Review B, 62, 4234-4237.
https://doi.org/10.1103/PhysRevB.62.4234

47. 47. Wenfang and Xie (2010) Optical Properties of an Off-Center Hydrogenic Impurity in a Spherical Quantum Dot with Gaussian Potential. Superlattices and Microstructures, 48, 239-247.
https://doi.org/10.1016/j.spmi.2010.04.015

48. 48. Hours, J., Senellart, P., Peter, E., Cavanna, A. and Bloch, J. (2005) Exciton Radiative Lifetime Controlled by the Lateral Confinement Energy in a Single Quantum Dot. Physical Review B, 71, Article ID: 161306.
https://doi.org/10.1103/PhysRevB.71.161306