Frobenius Series Solutions of the Schrodinger Equation with Various Types of Symmetric Hyperbolic Potentials in One Dimension

The Schrodinger equation (SE) for a certain class of symmetric hyperbolic potentials is solved with the aid of the Frobenius method (FM). The bound state energies are given as zeros of a calculable function. The calculated bound state energies are successively substituted into the recurrence relations for the expanding coefficients of the Frobenius series representing even and odd solutions in order to obtain wave functions associated with even and odd bound states. As illustrative examples, we consider the hyperbolic Poschl-Teller potential (HPTP) which is an exactly solvable potential, the Manning potential (MP) and a model of the Gaussian potential well (GPW). In each example, the bound state energies obtained by means of the FM are presented and compared with the exact results or the literature ones. In the case of the HPTP, we also make a comparison between exact bound state wave functions and the eigenfunctions obtained by means of the present approach. We find that our results are in good agreement with those given by other methods considered in this work, and that our class of potentials can be a perfect candidate to model the GPW.

Share and Cite:

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. doi: 10.4236/oalib.1104728. 1. Introduction

The exact solution of the Schrödinger equation (SE) can be obtained only for a few particular forms of potentials. In other cases, one has to appeal to approximations or numerical techniques.

Many approximation methods have been developed for solving problems in one-dimensional space, among which the power series method and its extension were known as the Fröbenius method (FM)  .

The power series method is the standard method for solving linear ordinary differential equations with variable coefficients. It gives solutions in the form of a power series expansion about any desired point $x={x}_{0}$ . The FM generalizes the power series one. It gives solutions in the form of a power series expansion, multiplied by a logarithm term $\mathrm{ln}|x-{x}_{0}|$ or a fractional power ${\left(x-{x}_{0}\right)}^{r}$  .

In this paper, we show that highly accurate solutions to the SE can be determined for various types of symmetric hyperbolic potentials by the use of the FM. The potentials we consider here belong to a class of hyperbolic potentials given by

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

where $K\in {ℕ}_{0}$ , $\alpha \in \right]0,+\infty \left[$ and coefficients ${B}_{j}$ , $j=1,2,\cdots ,K$ are not simultaneously zero. This potential can be reduced to the hyperbolic Pöschel-Teller potential (HPTP)    which occurs in the study of solitons   and the Manning potential (MP)  which can be used to study polyatomic molecular vibrations   . We can verify that it also contains the family of potentials ${V}_{p,q}\left(X\right)$ defined by

${V}_{p,q}\left(X\right)=-{V}_{0}\frac{{\mathrm{sinh}}^{2q}\left(\alpha X\right)}{{\mathrm{cosh}}^{2p}\left(\alpha X\right)}$ (2)

where ${V}_{0},\alpha \in \right]0,+\infty \left[$ and $q=1,2,\cdots ,p$ . Furthermore, the class of potentials

under consideration has another interesting aspect that it can be used to model the Gaussian function of the form

${f}_{G}\left(X\right)=-\mathrm{exp}\left(-{X}^{2}\right)$ (3)

by suitably choosing the values of $\alpha ,K$ and ${B}_{j}$ , $j=1,2,\cdots ,K$ .

Based on the foregoing, it can be said that the class of potentials (1) is composed of many types of hyperbolic potentials and several approximations of the Gaussian function. It therefore has the ability to explain various physical phenomena. Consequently, it is of considerable interest to find general solutions of its associated SE. In this work, we aim to present Fröbenius series solutions of this equation and discuss several members of the family of potentials under consideration.

The rest of this paper is organized as follows. In section 2, we convert the SE into a dimensionless eigenvalue problem which can be solved by means of the FM, and give recurrence relations for the expanding coefficients of the generalized power series for the wave functions. We also show how to determine bound state energies and wave functions. Furthermore, we derive the normalization constants for odd and even bound state eigenfunctions. In section 3, we present and comment results for two special cases, i.e., the case of the HPTP and that of the MP. In section 4, we apply the FM to a potential of type (1) modeling the function $-\mathrm{exp}\left(-{X}^{2}\right)$ , with an application to the Gaussian potential well (GPW). The conclusion is given in Section 5.

2. Derivation of Even and Odd Series Solutions

2.1. Dimensionless Schrödinger Equation

The time-independent Schrödinger equation for a particle of mass m that moves in one dimension $\left(-\infty under the effect of the potential (1) is

$-\frac{{\hslash }^{2}}{2m}\frac{{\text{d}}^{2}\psi \left(X\right)}{\text{d}{X}^{2}}+\left[\underset{j=1}{\overset{K}{\sum }}{B}_{j}{\text{sech}}^{2j}\left(\alpha X\right)\right]\psi \left(X\right)=E\psi \left(X\right)$ (4)

or

$\frac{{\text{d}}^{2}\psi \left(x\right)}{\text{d}{x}^{2}}+\left[\epsilon -\underset{j=1}{\overset{K}{\sum }}{b}_{j}{\text{sech}}^{2j}\left(x\right)\right]\psi \left(x\right)=0$ (5)

where

$x=\alpha X,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\epsilon =\frac{2m}{{\hslash }^{2}{\alpha }^{2}}E,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{b}_{j}=\frac{2m}{{\hslash }^{2}{\alpha }^{2}}{B}_{j}$ . (6)

Rewriting Equation (5) with a new variable of the form $\xi ={\text{sech}}^{2}\left(x\right)$ , such that the domain $-\infty maps to $0<\xi <1$ and

$\frac{{\text{d}}^{2}\psi \left(x\right)}{\text{d}{x}^{2}}=4{\xi }^{2}\left(1-\xi \right)\frac{{\text{d}}^{2}\psi }{\text{d}{\xi }^{2}}+\left(4\xi -6{\xi }^{2}\right)\frac{\text{d}\psi }{\text{d}\xi },$ (7)

we obtain

${\xi }^{2}\left(1-{\xi }^{2}\right)\frac{{\text{d}}^{2}\psi }{\text{d}{\xi }^{2}}+\xi \left(1-\frac{3}{2}\xi \right)\frac{\text{d}\psi }{\text{d}\xi }+\frac{1}{4}\left(\epsilon -\underset{j=1}{\overset{K}{\sum }}{b}_{j}{\xi }^{j}\right)\psi \left(\xi \right)=0.$ (8)

2.2. Expansion around a Regular-Singular Point

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

$\psi \left(\xi \right)={\xi }^{\delta }\underset{n=0}{\overset{\infty }{\sum }}{a}_{n}{\xi }^{n}$ , ${a}_{0}\ne 0$ (9)

or

$\psi \left(\xi \right)={\left(\xi -1\right)}^{\delta }\underset{n=0}{\overset{\infty }{\sum }}{a}_{n}{\left(\xi -1\right)}^{n}$ , ${a}_{0}\ne 0$ . (10)

Let us seek $\psi \left(\xi \right)$ in the form (10). We first make the substitution $\eta =\xi -1$ . The resulting differential equation, i.e.

$\left(\eta +2{\eta }^{2}+{\eta }^{3}\right)\frac{{\text{d}}^{2}\psi \left(\eta \right)}{\text{d}{\eta }^{2}}+\left(\frac{1}{2}+2\eta +\frac{3}{2}{\eta }^{2}\right)\frac{\text{d}\psi \left(\eta \right)}{\text{d}\eta }-\frac{1}{4}\left(\epsilon -\underset{j=0}{\overset{K}{\sum }}{c}_{j}{\eta }^{j}\right)\psi \left(\eta \right)=0,$ (11)

can be solved near $\eta =0$ by writing

$\psi \left(\eta \right)={\eta }^{\delta }\underset{n=0}{\overset{\infty }{\sum }}{a}_{n}{\eta }^{n}.$ (12)

It is worth noting that the ${c}_{j}$ coefficients that appear in Equation (11) are calculated from the two following expressions:

${c}_{0}=\underset{j=1}{\overset{K}{\sum }}{b}_{j},$ (13)

${c}_{i}=\underset{j=i}{\overset{K}{\sum }}\left(\begin{array}{l}j\\ i\end{array}\right){b}_{j}$ , $i=1,2,\cdots ,K$ (14)

where $\left(\begin{array}{c}j\\ i\end{array}\right)$ denotes the number of combinations of j objects taken i at a time.

Putting Equation (12) and its two first derivatives into Equation (11) we get, after some reduction:

$\begin{array}{l}\underset{n=0}{\overset{\infty }{\sum }}\left(n+\delta \right)\left(n+\delta -\frac{1}{2}\right){a}_{n}{\eta }^{n+\delta -1}+\underset{n=0}{\overset{\infty }{\sum }}\left[2{\left(n+\delta \right)}^{2}-\frac{\epsilon }{4}+\frac{{c}_{0}}{4}\right]{a}_{n}{\eta }^{n+\delta }\\ +\underset{n=0}{\overset{\infty }{\sum }}\left[\frac{3}{2}\left(n+\delta \right)+\frac{{c}_{1}}{4}\right]{a}_{n}{\eta }^{n+\delta +1}+\frac{{c}_{2}}{4}\underset{n=0}{\overset{\infty }{\sum }}{a}_{n}{\eta }^{n+\delta +2}\\ +\cdots +\frac{{c}_{K}}{4}\underset{n=0}{\overset{\infty }{\sum }}{a}_{n}{\eta }^{n+\delta +K}=0.\end{array}$ (15)

The smallest power is ${x}^{\delta -1}$ . It appears in the first series. By equating its coefficient to zero, we obtain the indicial equation

$\delta \left(\delta -\frac{1}{2}\right)=0$ (16)

whose solutions are ${\delta }_{1}=0$ and ${\delta }_{2}=1/2$ . We can therefore conclude that the two solutions of the 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$ $\left(X\to ±\infty \right)$ and $\eta =0$ $\left(X=0\right)$ . Clearly, only $\delta ={\delta }_{2}$ is acceptable for odd states since in this case $\psi \left(X=0\right)=0$ . On the other hand, $\delta ={\delta }_{1}$ is suitable for even states.

To sum up:

・ Odd solutions contain only the series with $\delta ={\delta }_{2}=1/2$ , and in the following will be denoted by

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

where the dependence on the eigenvalue is explicitly marked. The coefficients ${a}_{n}$ are given by the $\left(K+2\right)$ -term recurrence relation

$\begin{array}{c}{a}_{n+1}=\frac{-1}{\left(n+3/2\right)\left(n+1\right)}\left\{\left[2{\left(n+\frac{1}{2}\right)}^{2}-\frac{\epsilon }{4}+\frac{{c}_{0}}{4}\right]{a}_{n}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\left[n\left(n-\frac{1}{2}\right)+\frac{{c}_{1}}{4}\right]{a}_{n-1}+\frac{1}{4}\underset{j=2}{\overset{K}{\sum }}{c}_{j}{a}_{n-j}\right\}\end{array}$ (18)

with the understanding that ${a}_{-K}={a}_{-K+1}=\cdots ={a}_{-1}=0$ .

・ An even solution can be written as a power series

${\psi }_{e}\left(\eta ,\epsilon \right)=\underset{n=0}{\overset{\infty }{\sum }}{a}_{n}{\eta }^{n}$ (19)

where the coefficients ${a}_{n}$ are given by the $\left(K+2\right)$ -term recurrence relation

$\begin{array}{c}{a}_{n+1}=\frac{-1}{\left(n+1\right)\left(n+\frac{1}{2}\right)}\left\{\left(2{n}^{2}-\frac{\epsilon }{4}+\frac{{c}_{0}}{4}\right){a}_{n}+\left[\left(n-1\right)\left(n-\frac{1}{2}\right)+\frac{{c}_{1}}{4}\right]{a}_{n-1}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{1}{4}\underset{j=2}{\overset{K}{\sum }}{c}_{j}{a}_{n-j}\right\}\end{array}$ (20)

2.3. Bound State Energies and Wave Functions

The function ${\psi }_{o}\left(\eta ,\epsilon \right)$ has an important property that ${\psi }_{o}\left(\eta =0,\epsilon \right)=0$ , which can be regarded as a boundary condition for odd wave functions. For even wave functions, we choose ${\psi }_{e}\left(\eta =0,\epsilon \right)={a}_{0}$ as boundary condition at this point. These conditions are sufficient to ensure that ${\psi }_{o}\left(\eta ,\epsilon \right)$ and ${\psi }_{e}\left(\eta ,\epsilon \right)$ are linearly independent.

If ${\psi }_{o}\left(\eta ,\epsilon \right)$ (resp. ${\psi }_{e}\left(\eta ,\epsilon \right)$ ) is associated with a bound state, it has another important property, i.e. ${\psi }_{o}\left(\eta =-1,\epsilon \right)=0$ (resp. ${\psi }_{e}\left(\eta =-1,\epsilon \right)=0$ ), which can be considered as the boundary condition for bound state wave functions at $x=±\infty$ . It means that for arbitrary $\epsilon$ , the wave functions ${\psi }_{o}\left(\eta ,\epsilon \right)$ and ${\psi }_{e}\left(\eta ,\epsilon \right)$ obtained from Equations (17) and (19) are not square integrable. In other words, this condition determines the energies of the discrete spectrum. Substituting the calculated bound state energy to the recursion relation (18) or (20), the coefficients of the generalized series (17) or the power series (19) can be successively determined in order to obtain the wave function as a sum of the series.

We have to emphasize that in practice, the function ${\psi }_{o}\left(\eta =-1,\epsilon \right)$ (resp. ${\psi }_{e}\left(\eta =-1,\epsilon \right)$ ) has to be approximated by truncating the series in Equation (17) (resp. (19)) at suitably high order N. The truncated function is a polynomial of degree N in the variable $\epsilon$ whose zeros correspond to bound state values of $\epsilon$ . If for two neighboring values of $\epsilon$ the wave function takes values of different signs, we can deduce that one bound state reduced energy ( $\epsilon$ ) lies between these two values. The bound state reduced energy in question can be computed numerically by means of the dichotomy method  with an arbitrary chosen precision.

2.4. Bound State Wave Function Normalization Constants

In order to facilitate the calculation of odd and even wave function normalization constants, we consider only the first $N+1$ terms in the generalized series (17) and the power series (19), N being sufficiently large. More precisely, we write the wave functions ${\psi }_{o}\left(x,\epsilon \right)$ and ${\psi }_{e}\left(x,\epsilon \right)$ in the forms

${\psi }_{o}\left(X,\epsilon \right)={N}_{o}\left(\epsilon \right){\left(1-{\text{sech}}^{2}\left(\alpha X\right)\right)}^{1/2}\underset{n=0}{\overset{N}{\sum }}{\left(-1\right)}^{n}{a}_{n}^{odd}{\mathrm{tanh}}^{2n}\left(\alpha X\right)$ (21)

and

${\psi }_{e}\left(X,\epsilon \right)={N}_{e}\left(\epsilon \right)\underset{n=0}{\overset{N}{\sum }}{\left(-1\right)}^{n}{a}_{n}^{even}{\mathrm{tanh}}^{2n}\left(\alpha X\right)$ (22)

wherewe have used the fact that ${\left({\text{sech}}^{2}\left(\alpha X\right)-1\right)}^{n}={\left(-1\right)}^{n}{\left(1-{\text{sech}}^{2}\left(\alpha X\right)\right)}^{n}$ and $\left(1-{\text{sech}}^{2}\left(\alpha X\right)\right)={\mathrm{tanh}}^{2}\left(\alpha X\right)$ . Note that the parity of the wave functions

areexplicitly marked in the ${a}_{j}$ coefficients, and that ${N}_{o}\left(\epsilon \right)$ and ${N}_{e}\left(\epsilon \right)$ design the normalization constants for odd and even states respectively. These constants are obtained from the equations

$\underset{-\infty }{\overset{+\infty }{\int }}{\left[{\psi }_{o}\left(x,\epsilon \right)\right]}^{2}\text{d}x=1$ (23)

and

$\underset{-\infty }{\overset{+\infty }{\int }}{\left[{\psi }_{e}\left(x,\epsilon \right)\right]}^{2}\text{d}x=1$ . (24)

In view of the boundary condition for bound state wave functions at $x=±\infty$ , we have the following approximate equations:

$\underset{n=0}{\overset{N}{\sum }}{\left(-1\right)}^{n}{a}_{n}^{odd}\approx 0$ , $\underset{n=0}{\overset{N}{\sum }}{\left(-1\right)}^{n}{a}_{n}^{even}\approx 0$ . (25)

Using these approximations, we can rewrite Equations (21) and (22) as

${\psi }_{o}\left(X,\epsilon \right)={N}_{o}\left(\epsilon \right){\left(1-{\text{sech}}^{2}\left(\alpha X\right)\right)}^{1/2}\underset{n=1}{\overset{N}{\sum }}{\left(-1\right)}^{n+1}{a}_{n}^{odd}\left(1-{\mathrm{tanh}}^{2n}\left(\alpha X\right)\right)$ (26)

and

${\psi }_{e}\left(X,\epsilon \right)={N}_{o}\left(\epsilon \right)\underset{n=1}{\overset{N}{\sum }}{\left(-1\right)}^{n+1}{a}_{n}^{even}\left(1-{\mathrm{tanh}}^{2n}\left(\alpha X\right)\right)$ (27)

By inserting Equation (26) (resp. (27)) in the left hand side of Equation (23) (resp. (24)) and then using the remarkable identities

$\left(\underset{i=1}{\overset{m}{\sum }}{x}_{i}\right)\left(\underset{j=1}{\overset{n}{\sum }}{y}_{i}\right)=\underset{i=1}{\overset{m}{\sum }}\underset{j=1}{\overset{n}{\sum }}{x}_{i}{y}_{j}$ (28)

and

$1-{a}^{n}=\left(1-a\right)\underset{i=0}{\overset{n-1}{\sum }}{a}^{i},$ (29)

we obtain:

${\left({N}_{o}\left(\epsilon \right)\right)}^{2}\underset{j=1}{\overset{N}{\sum }}\underset{p=1}{\overset{N}{\sum }}{\left(-1\right)}^{j+p}{a}_{j}^{odd}{a}_{p}^{odd}\underset{-\infty }{\overset{+\infty }{\int }}f\left(X\right)\text{d}X=1,$ (30)

${\left({N}_{e}\left(\epsilon \right)\right)}^{2}\underset{j=1}{\overset{N}{\sum }}\underset{p=1}{\overset{N}{\sum }}{\left(-1\right)}^{j+p}{a}_{j}^{even}{a}_{p}^{even}\underset{-\infty }{\overset{+\infty }{\int }}g\left(X\right)\text{d}X=1$ (31)

with

$f\left(X\right)=\left(1-{\mathrm{tanh}}^{2}\left(\alpha x\right)\right)\underset{q=0}{\overset{j-1}{\sum }}{\mathrm{tanh}}^{2q+2}\left(\alpha x\right)\left(1-{\mathrm{tanh}}^{2p}\left(\alpha x\right)\right),$ (32)

$g\left(X\right)=\left(1-{\mathrm{tanh}}^{2}\left(\alpha X\right)\right)\underset{q=0}{\overset{j-1}{\sum }}{\mathrm{tanh}}^{2q}\left(\alpha X\right)\left(1-{\mathrm{tanh}}^{2p}\left(\alpha X\right)\right).$ (33)

Changing variables from x to $u=\mathrm{tanh}\left(\alpha x\right)$ and using the fact that

$\underset{-1}{\overset{1}{\int }}{u}^{2n}\text{d}u=2\underset{0}{\overset{1}{\int }}{u}^{2n}\text{d}u=\frac{2}{2n+1}$ , (34)

we find

$\frac{2}{\alpha }{\left({N}_{o}\left(\epsilon \right)\right)}^{2}\underset{j=1}{\overset{N}{\sum }}\underset{p=1}{\overset{N}{\sum }}{\left(-1\right)}^{j+p}{a}_{j}^{odd}{a}_{p}^{odd}\left(\underset{q=0}{\overset{j-1}{\sum }}{I}_{q+1}-\underset{q=0}{\overset{j-1}{\sum }}{I}_{q+p+1}\right)=1$ , (35)

$\frac{2}{\alpha }{\left({N}_{e}\left(\epsilon \right)\right)}^{2}\underset{j=1}{\overset{N}{\sum }}\underset{p=1}{\overset{N}{\sum }}{\left(-1\right)}^{j+p}{a}_{j}^{even}{a}_{p}^{even}\left(\underset{q=0}{\overset{j-1}{\sum }}{I}_{q}-\underset{q=0}{\overset{j-1}{\sum }}{I}_{q+p}\right)=1$ (36)

where ${I}_{q}=\underset{0}{\overset{1}{\int }}{u}^{2q}\text{d}u={\left(2q+1\right)}^{-1}$ .

It follows immediately from Equations (35) and (36) that

${N}_{o}\left(\epsilon \right)=\sqrt{\frac{\alpha }{2}}{\left[\underset{j=1}{\overset{N}{\sum }}\underset{p=1}{\overset{N}{\sum }}{\left(-1\right)}^{j+p}{a}_{j}^{odd}{a}_{p}^{odd}\left(\underset{q=0}{\overset{j-1}{\sum }}{I}_{q+1}-\underset{q=0}{\overset{j-1}{\sum }}{I}_{q+p+1}\right)\right]}^{-1/2}$ (37)

and

${N}_{e}\left(\epsilon \right)=\sqrt{\frac{\alpha }{2}}{\left[\underset{j=1}{\overset{N}{\sum }}\underset{p=1}{\overset{N}{\sum }}{\left(-1\right)}^{j+p}{a}_{j}^{even}{a}_{p}^{even}\left(\underset{q=0}{\overset{j-1}{\sum }}{I}_{q}-\underset{q=0}{\overset{j-1}{\sum }}{I}_{q+p}\right)\right]}^{-1/2}$ . (38)

3. Numerical Results for Two Special Cases

3.1. Hyperbolic Pöschl-Teller Potential

In order to test the reliability and efficiency of the FM for the class of potentials (1), we here consider the HPTP

$V\left(X\right)=-{V}_{0}{\text{sech}}^{2}\left(\alpha X\right)$ , ${V}_{0},\alpha >0$ (39)

which corresponds to the special case where $K=1$ and ${B}_{1}=-{V}_{0}$ . This potential is exactly solvable. The bound state energies are given by   :

${E}_{n}=-\frac{{\hslash }^{2}{\alpha }^{2}}{2m}\left[{\left(n+\frac{1}{2}-\sqrt{\nu +\frac{1}{4}}\right)}^{2}\right],\text{\hspace{0.17em}}n=0,1,\cdots ;\text{\hspace{0.17em}}n<\sqrt{\nu +\frac{1}{4}}-\frac{1}{2}$ (40)

with $\nu =2m{V}_{0}/\left({\hslash }^{2}{\alpha }^{2}\right)$ . The associated wave functions can be written as

${\psi }_{n}\left(x\right)={N}_{n}{\text{sech}}^{{\beta }_{n}}\left(\alpha x\right){}_{2}F{}_{1}\left(-n,\sqrt{4\nu +1}-n;1+{\beta }_{n};\frac{1}{2}\left(1-\mathrm{tanh}\left(\alpha x\right)\right)\right)$ (41)

where

${\beta }_{n}=\sqrt{\nu +\frac{1}{4}}-\frac{1}{2}-n$ . (42)

The normalization constants ${N}_{n}$ are obtained from the equation $\underset{-\infty }{\overset{+\infty }{\int }}{\left[{\psi }_{n}\left(x\right)\right]}^{2}\text{d}x=1$ and can be expressed as follows:

${N}_{n}={\left[\underset{\begin{array}{l}j=0\\ p=0\end{array}}{\overset{n}{\sum }}D\left(n,j\right)D\left(n,p\right)\frac{{4}^{{\beta }_{n}}}{2\alpha }\frac{\Gamma \left({\beta }_{n}+j+p\right)\Gamma \left({\beta }_{n}\right)}{\Gamma \left(2{\beta }_{n}+j+p\right)}\right]}^{-1/2}$ (43)

with

$D\left(n,j\right)=\frac{{\left(-n\right)}_{j}{\left(\sqrt{4\nu +1}-n\right)}_{j}}{{\left(1+{\beta }_{n}\right)}_{j}j!}$ (44)

where ${\left(a\right)}_{j}$ is the Pochhammer symbol, i.e.,

${\left(a\right)}_{j}=a\left(a+1\right)\cdots \left(a+j-1\right)=\frac{\Gamma \left(a+j\right)}{\Gamma \left(a\right)}$ , ${\left(a\right)}_{0}=1$ . (45)

Table 1 displays bound state energies for $m={m}_{e}$ , ${V}_{0}=-1.6128$ and $\alpha =0.24$ in atomic units ( $\hslash =1$ and ${m}_{e}=1$ ), where ${m}_{e}$ designs the electron mass. The exact energies ( ${E}_{n}^{exact}$ ) are reported in the second column; those calculated by means of the FM are given in the third, fourth and fifth columns for three values of N ( ${N}_{1}=500$ , ${N}_{2}=1000$ and ${N}_{3}=2000$ ). We have chosen ${a}_{0}=1$ . As expected, the FM leads to more accurate results as the truncation order N goes to $\infty$ . Note that our calculations have been performed with quadruple precision.

Figure 1 shows four normalized bound state wave functions as functions of X. We remark that there is a good agreement between analytical results and those obtained by the use of the FM, the number of terms in the series (17) and (19) being ${N}_{3}$ .

3.2. Manning Potential

As an interesting application of the FM to a physical problem, we consider the MP that has been used in discussing vibrational energies of polyatomic molecules such as ammonia  , formamide  , cyanamide  and their deuterated species. In order to study one of the vibrational normal modes of $N{H}_{3}$ and $N{D}_{3}$ , Manning  proposed a potential defined by the expression

$\frac{V}{hc}=\frac{1}{k{\rho }^{2}}\left[-\frac{\beta }{2}\left(\frac{\beta }{2}+\frac{1}{2}\right){\text{sech}}^{2}\frac{X}{2\rho }-D\left({\text{sech}}^{2}\frac{X}{2\rho }-{\text{sech}}^{4}\frac{X}{2\rho }\right)\right]{\text{cm}}^{-1}$ (46)

where $k=8{\text{π}}^{2}c\mu /h$ , $\mu$ is the reduced mass of the molecule for the type of vibration considered; X is the distance from the nitrogen atom to the plane of the hydrogen atoms; and $\rho$ , $\beta$ and D are arbitrary constants. The SE for the vibration is

$-\frac{{\hslash }^{2}}{2\mu }\frac{{\text{d}}^{2}\psi }{\text{d}{X}^{2}}+hc\left[-E-\frac{\frac{\beta }{2}\left(\frac{\beta }{2}+\frac{1}{2}\right)+D}{k{\rho }^{2}}{\text{sech}}^{2}\frac{X}{2\rho }+\frac{D}{k{\rho }^{2}}{\text{sech}}^{4}\frac{X}{2\rho }\right]\psi =0$ (47)

Figure 1. The normalized wave functions of the ground state, the first, fifth and sixth excited states contained within the HPTP described by the parameters ${V}_{0}=1.6228$ atomic units and $\alpha =0.24$ atomic units.

Table 1. Bound state energies (in atomic units) for the one-dimensional HPTP with ${V}_{0}=-1.6128$ and $\alpha =0.24$ .

where E is the energy in cm−1. If we define $x=X/2\rho$ , we obtain the dimensionless Schrödinger Equation (5) with

$\epsilon =4k{\rho }^{2}E$ , $K=2$ , ${b}_{1}=-\left[\beta \left(\beta +1\right)+4D\right]$ and ${b}_{2}=4D$ . (48)

Following Manning  , we choose $\beta =70$ (91.4), $D=1920$ (3261) and $k{\rho }^{2}=2.885×{10}^{-2}$ $\left(4.899×{10}^{-2}\right)$ for $N{H}_{3}$ $\left(N{D}_{3}\right)$ .

Figure 2 shows the MP in cm−1 for both sets of parameters, as a function of x. This potential exhibits two minima and a barrier with maximum at x = 0. Although the barriers appear to be rather low, there are several energies between a minimum of the potential and the maximum of the barrier. We have calculated these energies by means of the FM. Our results are reported in Table 2 and Table 3, and compared with Manning’s ones  and those calculated by the use of the Wronskian method (WM)  . The remarkable agreement between the results of the two completely independent approaches shows that the present calculation is sufficiently accurate. Note that ${\stackrel{˜}{\epsilon }}_{n}$ is defined by ${\stackrel{˜}{\epsilon }}_{n}=k{\rho }^{2}{E}_{n}/2=\frac{{\epsilon }_{n}}{8}$ .

4. Application to the Gaussian Potential Well

As another interesting application of the FM, we consider the one-dimensional SE for the attractive Gaussian potential

${V}_{G}\left(X\right)=-{V}_{0}\mathrm{exp}\left(-{X}^{2}/{X}_{0}^{2}\right)$ , $\left({V}_{0},{X}_{0}>0\right)$ (49)

i.e.

$-\frac{{\hslash }^{2}}{2m}\frac{{\text{d}}^{2}\psi \left(X\right)}{\text{d}{X}^{2}}-\left[E+{V}_{0}\mathrm{exp}\left(-{X}^{2}/{X}_{0}^{2}\right)\right]\psi \left(X\right)=0$ . (50)

Figure 2. Manning potential for NH3 and ND3.

Table 2. Energies of the Manning potential for NH3.

Table 3. Energies of the Manning potential for ND3.

The GPW is not a long range potential, because it falls off faster then $1/{X}^{2}$ . However, it has a crucial advantage over the widely used “finite square well” potential in that it is continuous throughout the entire range $-\infty , whereas the finite square well goes to zero discontinuously. Furthermore, the Gaussian well of the form (49) has been used by Bardsley and Comella  for describing model neutral atoms with a few bound states.

It is worth noting that by an estimation based on the WKB method, the number of bound states supported by the potential (49) is approximated by 

$nbst=⌊2{X}_{0}\sqrt{{V}_{0}/\text{π}}+\frac{1}{2}⌋+1$ (51)

where $⌊\text{ }⌋$ represents the floor, which for positive numbers is simply the integer

part. But so far, to the best of our knowledge, exact analytical forms of the bound state energies and wave functions of this potential have not been reported. Therefore, it is of considerable interest to find a good approximation for the solution of the SE (50) which can be rewritten as

$-\frac{{\hslash }^{2}}{2m{X}_{0}^{2}}\frac{{\text{d}}^{2}\psi \left(\chi \right)}{\text{d}{\chi }^{2}}-\left[E+{V}_{0}\mathrm{exp}\left(-{\chi }^{2}\right)\right]\psi \left(\chi \right)=0$ (52)

by changing variables from X to $\chi =X/{X}_{0}$ .

The approximation solution we propose here relies on approximating the Gaussian function $-\mathrm{exp}\left(-{\xi }^{2}\right)$ by a potential of type (1). More precisely, we write

$-\mathrm{exp}\left(-{\chi }^{2}\right)\approx \underset{j=1}{\overset{K}{\sum }}{B}_{j}{\text{sech}}^{2j}\left(\alpha \chi \right)$ (53)

where $\alpha ,{B}_{1},{B}_{2},\cdots ,{B}_{K}$ are adjustable parameters and K is a positive integer.

In order to show that the above approximation for the Gaussian function is suitable, we choose K = 8 and impose α = 0.35. We then discretize the problem in the interval $\left[{\chi }_{\mathrm{min}},{\chi }_{\mathrm{max}}\right]$ and evaluate the function $-\mathrm{exp}\left(-{\chi }^{2}\right)$ at ${\chi }_{k}={\chi }_{\mathrm{min}}+kh$ (for $k=1,2,\cdots ,p$ ) where p is the number of mesh points, h is the step size, and ${\chi }_{\mathrm{min}}$ and ${\chi }_{\mathrm{max}}$ are the minimum and maximum values of the variable $\chi$ . We finally use the maple 18 Fit function to calculate the ${B}_{j}$ coefficients. With ${\chi }_{\mathrm{max}}=6=-{\chi }_{\mathrm{min}}$ and $h=0.01$ , we find:

$\left\{\begin{array}{l}{B}_{1}=-0.00056186637277292680096387645307\\ {B}_{2}=0.023603992393495810256516883997296\\ {B}_{3}=-0.3203936279047603522539856166491\\ {B}_{4}=2.031520751522174095847244814245208\\ {B}_{5}=-6.85096256835970745670929664345292\\ {B}_{6}=12.49340526371948028270841849878273\\ {B}_{7}=-10.9323757769502446942988487237498\\ {B}_{8}=2.555762222221351947020248045060859\end{array}$ (54)

Using these parameters, we plot the variations of ${f}_{G}^{approx}\left(\chi \right)\equiv \underset{j=1}{\overset{8}{\sum }}{B}_{j}{\text{sech}}^{2j}\left(\alpha \chi \right)$ and ${f}_{G}^{exact}\left(\chi \right)\equiv -\mathrm{exp}\left(-{\chi }^{2}\right)$ with $\chi$ . We also consider the difference between the two functions, i.e. $\Delta {f}_{G}\left(\chi \right)={f}_{G}^{exact}\left(\chi \right)-{f}_{G}^{approx}\left(\chi \right)$ .

It is obvious from Figure 3 that the approximation (53) is a very good one. ${f}_{G}^{exact}\left(\chi \right)$ and ${f}_{G}^{approx}\left(\chi \right)$ seem to coincide because the difference between them is of order 10−6 and hence very small.

By inserting this approximation in Equation (52) and then changing variables from $\chi$ to $x=\alpha \chi$ , we obtain the dimensionless Schrödinger Equation (5) with

$\epsilon =\frac{2m{X}_{0}^{2}}{{\hslash }^{2}{\alpha }^{2}}E$ end ${b}_{j}=\frac{2m{X}_{0}^{2}}{{\hslash }^{2}{\alpha }^{2}}{V}_{0}{B}_{j}$ , $j=1,2,\cdots ,K$ . (55)

In Table 4, we present bound state energies obtained with the FM for ${V}_{0}=2.5$ and ${X}_{0}=3.05$ atomic units. This table contains also numerical results $\left({E}_{n}^{diag}\right)$ obtained with the aid of the exact Hamiltonian diagonalization on the finite-real basis ${\left\{{\varphi }_{j}^{M}\left(X\right)\right\}}_{j=M}^{Ns+M-1}$ where

${\varphi }_{j}^{M}\left(X\right)=\sqrt{\frac{\stackrel{˜}{\alpha }M\left(j-M\right)!}{\left(j+M\right)!}}{P}_{j}^{M}\left(\mathrm{tanh}\left(\stackrel{˜}{\alpha }X\right)\right)$ , $j\ge M$ , $\stackrel{˜}{\alpha }\ge 0$ . (56)

Table 4. Bound state energies of the Gaussian potential well for ${V}_{0}=2.5$ atomic units and ${X}_{0}=3.05$ atomic units.

Figure 3. Plots of the Gaussian function ( ${f}_{G}^{exact}$ ), its approximation ( ${f}_{G}^{approx}$ ) and the difference between them ( $\Delta {f}_{G}$ ).

We have chosen ${N}_{s}=500$ as the number of basis functions, $M=2$ as the order of the associated Lengendre polynomials intervening in Equation (56), and ${a}_{0}={10}^{-40}$ . It is clear that our results are in good agreement with those obtained by means of the exact Hamiltonian diagonalization.

5. Conclusion

We have shown that the application of the Fröbenius method to the class of symmetric hyperbolic potentials of the form $V\left(X\right)=\underset{j=1}{\overset{K}{\sum }}{B}_{j}{\text{sech}}^{2j}\left(\alpha X\right)$ allows an easy determination of the bound state energies and wave functions. The bound state energies associated with these potentials have been obtained numerically as zeros of a function, calculated from its power series representation. Since the potential $V\left(X\right)$ is even, the problem of computing its spectrum has been split into two calculations, dealing respectively with odd and even states. In order to demonstrate the performance of the FM on the potentials $V\left(X\right)$ , we have considered two special cases, namely the case of the HPTP which is an exactly solvable potential, and the case of the MP suggested by several authors to study polyatomic molecular vibrations. The numerical and analytical results presented in this work show that our results are sufficiently accurate and that our family of potentials can be used to model the GPW of the form ${V}_{G}\left(X\right)=-{V}_{0}\mathrm{exp}\left(-{X}^{2}\right)$ by suitably choosing the values of $\alpha$ , K and the ${B}_{j}$ coefficients.

Conflicts of Interest

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

Conflicts of Interest

The authors declare no conflicts of interest.

  Kreyszig, G.E. (2011) Advanced Engineering Mathematics. 10th Edition, John Wiley & Sons Ltd., Amsterdam.  Greenberg, M.D. (1998) Advanced Engineering Mathematics. 2nd Edition, Simon & Schuster Asia Pte, Ltd., Singapore.  Nieto, M.M. (1978) Exact Wave-Function Normalization Constants for the and Poschel-Teller Potentials. Physical Review A, 17, 1273-1283. https://doi.org/10.1103/PhysRevA.17.1273  Flügge, S. (1971) Practical Quantum Mechanics I. Springer-Verlag, Berlin.  Nieto, M.M. (1978) Exact Wave-Function.  Brezon, E., Le Guillon, J.C. and Zinn-Justin, J. (1077) Perturbation Theory of Large Order. I. The Interaction. Physical Review D, 15, 1544-1557. https://doi.org/10.1103/PhysRevD.15.1544  Lamb, G.L. (1980) Elements of Soliton Theory. Jr. Editeur, New York.  Manning, M.F. (1935) Energy Levels of a Symmetrical Double Minima Problem with Applications to the and Moleculs. The Journal of Chemical Physics, 3, 136-138. https://doi.org/10.1063/1.1749619  Townes, C.H. and Schawlow, A.L. (1975) Microwave Spectroscopy. Dover Publications, New York.  Kincaid, D. and Chemey, W. (2002) Numerical Analysis: The Mathematics of Scientific Computing. 3rd Edition, Brooks/Cole, Pacific Groove, CA.  Morse, P.M. and Feshbach, H. (1953) Methods of Theoretical Physics, MacGraw-Holl, 2 Pennsylvania Plaza, New York City.  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  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  Fernández, F.M. (2011) Wronskian Method for Bound States. European Journal of Physics, 32, 723-732. https://doi.org/10.1088/0143-0807/32/3/008  Bardsley, J.N. and Comella, M.J. (1989) Ac Stark Effect for Short-Range Potentials with Intense Electromagnetic Fields. Physical Review A, 39, 2252-2255. https://doi.org/10.1103/PhysRevA.39.2252  Saikat, N. (2010) The Quantum Gaussian Well. American Journal of Physics, 78, 1341-1344. https://doi.org/10.1119/1.3474665 