^{1}

^{1}

^{1}

^{1}

^{1}

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.

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 ln | x − x 0 | or a fractional power ( x − x 0 ) 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 ( X ) = ∑ j = 1 K B j sech 2 j ( α X ) (1)

where K ∈ ℕ 0 , α ∈ ] 0 , + ∞ [ and coefficients B j , j = 1 , 2 , ⋯ , K are not simultaneously zero. This potential can be reduced to the hyperbolic Pöschel-Teller potential (HPTP) [

V p , q ( X ) = − V 0 sinh 2 q ( α X ) cosh 2 p ( α X ) (2)

where V 0 , α ∈ ] 0 , + ∞ [ and q = 1 , 2 , ⋯ , 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 ( X ) = − exp ( − X 2 ) (3)

by suitably choosing the values of α , K and B j , j = 1 , 2 , ⋯ , 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 − exp ( − X 2 ) , with an application to the Gaussian potential well (GPW). The conclusion is given in Section 5.

The time-independent Schrödinger equation for a particle of mass m that moves in one dimension ( − ∞ < X < + ∞ ) under the effect of the potential (1) is

− ℏ 2 2 m d 2 ψ ( X ) d X 2 + [ ∑ j = 1 K B j sech 2 j ( α X ) ] ψ ( X ) = E ψ ( X ) (4)

or

d 2 ψ ( x ) d x 2 + [ ε − ∑ j = 1 K b j sech 2 j ( x ) ] ψ ( x ) = 0 (5)

where

x = α X , ε = 2 m ℏ 2 α 2 E , b j = 2 m ℏ 2 α 2 B j . (6)

Rewriting Equation (5) with a new variable of the form ξ = sech 2 ( x ) , such that the domain − ∞ < x < + ∞ maps to 0 < ξ < 1 and

d 2 ψ ( x ) d x 2 = 4 ξ 2 ( 1 − ξ ) d 2 ψ d ξ 2 + ( 4 ξ − 6 ξ 2 ) d ψ d ξ , (7)

we obtain

ξ 2 ( 1 − ξ 2 ) d 2 ψ d ξ 2 + ξ ( 1 − 3 2 ξ ) d ψ d ξ + 1 4 ( ε − ∑ j = 1 K b j ξ j ) ψ ( ξ ) = 0. (8)

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

ψ ( ξ ) = ξ δ ∑ n = 0 ∞ a n ξ n , a 0 ≠ 0 (9)

or

ψ ( ξ ) = ( ξ − 1 ) δ ∑ n = 0 ∞ a n ( ξ − 1 ) n , a 0 ≠ 0 . (10)

Let us seek ψ ( ξ ) in the form (10). We first make the substitution η = ξ − 1 . The resulting differential equation, i.e.

( η + 2 η 2 + η 3 ) d 2 ψ ( η ) d η 2 + ( 1 2 + 2 η + 3 2 η 2 ) d ψ ( η ) d η − 1 4 ( ε − ∑ j = 0 K c j η j ) ψ ( η ) = 0 , (11)

can be solved near η = 0 by writing

ψ ( η ) = η δ ∑ n = 0 ∞ a n η 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 = ∑ j = 1 K b j , (13)

c i = ∑ j = i K ( j i ) b j , i = 1 , 2 , ⋯ , K (14)

where ( j i ) 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:

∑ n = 0 ∞ ( n + δ ) ( n + δ − 1 2 ) a n η n + δ − 1 + ∑ n = 0 ∞ [ 2 ( n + δ ) 2 − ε 4 + c 0 4 ] a n η n + δ + ∑ n = 0 ∞ [ 3 2 ( n + δ ) + c 1 4 ] a n η n + δ + 1 + c 2 4 ∑ n = 0 ∞ a n η n + δ + 2 + ⋯ + c K 4 ∑ n = 0 ∞ a n η n + δ + K = 0. (15)

The smallest power is x δ − 1 . It appears in the first series. By equating its coefficient to zero, we obtain the indicial equation

δ ( δ − 1 2 ) = 0 (16)

whose solutions are δ 1 = 0 and δ 2 = 1 / 2 . We can therefore conclude that the two solutions of the SE obtained as generalized series, one with δ = δ 1 and the other with δ = δ 2 , are linearly independent. The value of δ determines the behavior of ψ ( η ) for η → − 1 ( X → ± ∞ ) and η = 0 ( X = 0 ) . Clearly, only δ = δ 2 is acceptable for odd states since in this case ψ ( X = 0 ) = 0 . On the other hand, δ = δ 1 is suitable for even states.

To sum up:

・ Odd solutions contain only the series with δ = δ 2 = 1 / 2 , and in the following will be denoted by

ψ 0 ( η , ε ) = η δ 2 ∑ n = 0 ∞ a n η n (17)

where the dependence on the eigenvalue is explicitly marked. The coefficients a n are given by the ( K + 2 ) -term recurrence relation

a n + 1 = − 1 ( n + 3 / 2 ) ( n + 1 ) { [ 2 ( n + 1 2 ) 2 − ε 4 + c 0 4 ] a n + [ n ( n − 1 2 ) + c 1 4 ] a n − 1 + 1 4 ∑ j = 2 K c j a n − j } (18)

with the understanding that a − K = a − K + 1 = ⋯ = a − 1 = 0 .

・ An even solution can be written as a power series

ψ e ( η , ε ) = ∑ n = 0 ∞ a n η n (19)

where the coefficients a n are given by the ( K + 2 ) -term recurrence relation

a n + 1 = − 1 ( n + 1 ) ( n + 1 2 ) { ( 2 n 2 − ε 4 + c 0 4 ) a n + [ ( n − 1 ) ( n − 1 2 ) + c 1 4 ] a n − 1 + 1 4 ∑ j = 2 K c j a n − j } (20)

The function ψ o ( η , ε ) has an important property that ψ o ( η = 0 , ε ) = 0 , which can be regarded as a boundary condition for odd wave functions. For even wave functions, we choose ψ e ( η = 0 , ε ) = a 0 as boundary condition at this point. These conditions are sufficient to ensure that ψ o ( η , ε ) and ψ e ( η , ε ) are linearly independent.

If ψ o ( η , ε ) (resp. ψ e ( η , ε ) ) is associated with a bound state, it has another important property, i.e. ψ o ( η = − 1 , ε ) = 0 (resp. ψ e ( η = − 1 , ε ) = 0 ), which can be considered as the boundary condition for bound state wave functions at x = ± ∞ . It means that for arbitrary ε , the wave functions ψ o ( η , ε ) and ψ e ( η , ε ) 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 ψ o ( η = − 1 , ε ) (resp. ψ e ( η = − 1 , ε ) ) 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 ε whose zeros correspond to bound state values of ε . If for two neighboring values of ε the wave function takes values of different signs, we can deduce that one bound state reduced energy ( ε ) lies between these two values. The bound state reduced energy in question can be computed numerically by means of the dichotomy method [

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 ψ o ( x , ε ) and ψ e ( x , ε ) in the forms

ψ o ( X , ε ) = N o ( ε ) ( 1 − sech 2 ( α X ) ) 1 / 2 ∑ n = 0 N ( − 1 ) n a n o d d tanh 2 n ( α X ) (21)

and

ψ e ( X , ε ) = N e ( ε ) ∑ n = 0 N ( − 1 ) n a n e v e n tanh 2 n ( α X ) (22)

wherewe have used the fact that ( sech 2 ( α X ) − 1 ) n = ( − 1 ) n ( 1 − sech 2 ( α X ) ) n and ( 1 − sech 2 ( α X ) ) = tanh 2 ( α X ) . Note that the parity of the wave functions

areexplicitly marked in the a j coefficients, and that N o ( ε ) and N e ( ε ) design the normalization constants for odd and even states respectively. These constants are obtained from the equations

∫ − ∞ + ∞ [ ψ o ( x , ε ) ] 2 d x = 1 (23)

and

∫ − ∞ + ∞ [ ψ e ( x , ε ) ] 2 d x = 1 . (24)

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

∑ n = 0 N ( − 1 ) n a n o d d ≈ 0 , ∑ n = 0 N ( − 1 ) n a n e v e n ≈ 0 . (25)

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

ψ o ( X , ε ) = N o ( ε ) ( 1 − sech 2 ( α X ) ) 1 / 2 ∑ n = 1 N ( − 1 ) n + 1 a n o d d ( 1 − tanh 2 n ( α X ) ) (26)

and

ψ e ( X , ε ) = N o ( ε ) ∑ n = 1 N ( − 1 ) n + 1 a n e v e n ( 1 − tanh 2 n ( α X ) ) (27)

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

( ∑ i = 1 m x i ) ( ∑ j = 1 n y i ) = ∑ i = 1 m ∑ j = 1 n x i y j (28)

and

1 − a n = ( 1 − a ) ∑ i = 0 n − 1 a i , (29)

we obtain:

( N o ( ε ) ) 2 ∑ j = 1 N ∑ p = 1 N ( − 1 ) j + p a j o d d a p o d d ∫ − ∞ + ∞ f ( X ) d X = 1 , (30)

( N e ( ε ) ) 2 ∑ j = 1 N ∑ p = 1 N ( − 1 ) j + p a j e v e n a p e v e n ∫ − ∞ + ∞ g ( X ) d X = 1 (31)

with

f ( X ) = ( 1 − tanh 2 ( α x ) ) ∑ q = 0 j − 1 tanh 2 q + 2 ( α x ) ( 1 − tanh 2 p ( α x ) ) , (32)

g ( X ) = ( 1 − tanh 2 ( α X ) ) ∑ q = 0 j − 1 tanh 2 q ( α X ) ( 1 − tanh 2 p ( α X ) ) . (33)

Changing variables from x to u = tanh ( α x ) and using the fact that

∫ − 1 1 u 2 n d u = 2 ∫ 0 1 u 2 n d u = 2 2 n + 1 , (34)

we find

2 α ( N o ( ε ) ) 2 ∑ j = 1 N ∑ p = 1 N ( − 1 ) j + p a j o d d a p o d d ( ∑ q = 0 j − 1 I q + 1 − ∑ q = 0 j − 1 I q + p + 1 ) = 1 , (35)

2 α ( N e ( ε ) ) 2 ∑ j = 1 N ∑ p = 1 N ( − 1 ) j + p a j e v e n a p e v e n ( ∑ q = 0 j − 1 I q − ∑ q = 0 j − 1 I q + p ) = 1 (36)

where I q = ∫ 0 1 u 2 q d u = ( 2 q + 1 ) − 1 .

N o ( ε ) = α 2 [ ∑ j = 1 N ∑ p = 1 N ( − 1 ) j + p a j o d d a p o d d ( ∑ q = 0 j − 1 I q + 1 − ∑ q = 0 j − 1 I q + p + 1 ) ] − 1 / 2 (37)

N e ( ε ) = α 2 [ ∑ j = 1 N ∑ p = 1 N ( − 1 ) j + p a j e v e n a p e v e n ( ∑ q = 0 j − 1 I q − ∑ q = 0 j − 1 I q + p ) ] − 1 / 2 . (38)

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

V ( X ) = − V 0 sech 2 ( α X ) , V 0 , α > 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 = − ℏ 2 α 2 2 m [ ( n + 1 2 − ν + 1 4 ) 2 ] , n = 0 , 1 , ⋯ ; n < ν + 1 4 − 1 2 (40)

with ν = 2 m V 0 / ( ℏ 2 α 2 ) . The associated wave functions can be written as

ψ n ( x ) = N n sech β n ( α x ) F 2 1 ( − n , 4 ν + 1 − n ; 1 + β n ; 1 2 ( 1 − tanh ( α x ) ) ) (41)

where

β n = ν + 1 4 − 1 2 − n . (42)

The normalization constants N n are obtained from the equation ∫ − ∞ + ∞ [ ψ n ( x ) ] 2 d x = 1 and can be expressed as follows:

N n = [ ∑ j = 0 p = 0 n D ( n , j ) D ( n , p ) 4 β n 2 α Γ ( β n + j + p ) Γ ( β n ) Γ ( 2 β n + j + p ) ] − 1 / 2 (43)

with

D ( n , j ) = ( − n ) j ( 4 ν + 1 − n ) j ( 1 + β n ) j j ! (44)

where ( a ) j is the Pochhammer symbol, i.e.,

( a ) j = a ( a + 1 ) ⋯ ( a + j − 1 ) = Γ ( a + j ) Γ ( a ) , ( a ) 0 = 1 . (45)

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 [

V h c = 1 k ρ 2 [ − β 2 ( β 2 + 1 2 ) sech 2 X 2 ρ − D ( sech 2 X 2 ρ − sech 4 X 2 ρ ) ] cm − 1 (46)

where k = 8 π 2 c μ / h , μ 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 ρ , β and D are arbitrary constants. The SE for the vibration is

− ℏ 2 2 μ d 2 ψ d X 2 + h c [ − E − β 2 ( β 2 + 1 2 ) + D k ρ 2 sech 2 X 2 ρ + D k ρ 2 sech 4 X 2 ρ ] ψ = 0 (47)

n | E n e x a c t | E n F M (N = 500) | E n F M (N = 1000)_{ } | E n F M (N = 2000)_{ } |
---|---|---|---|---|

0 | −1.4112 | −1.41120000000000 | −1.41120000000000 | −1.41120000000000 |

1 | −1.0368 | −1.03680000000000 | −1.03680000000000 | −1.03680000000000 |

2 | −0.7200 | −0.71999999999792 | −0.71999999999994 | −0.72000000000000 |

3 | −0.4608 | −0.46080000000000 | −0.46080000000000 | −0.46080000000000 |

4 | −0.2592 | −0.25920020140700 | −0.25920002534412 | −0.25920000317865 |

5 | −0.1152 | −0.11520000000000 | −0.11520000000000 | −0.11520000000000 |

6 | −0.0288 | −0.02711385566710 | −0.02797214839035 | −0.02839059824400 |

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

ε = 4 k ρ 2 E , K = 2 , b 1 = − [ β ( β + 1 ) + 4 D ] and b 2 = 4 D . (48)

Following Manning [

^{−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

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

V G ( X ) = − V 0 exp ( − X 2 / X 0 2 ) , ( V 0 , X 0 > 0 ) (49)

i.e.

− ℏ 2 2 m d 2 ψ ( X ) d X 2 − [ E + V 0 exp ( − X 2 / X 0 2 ) ] ψ ( X ) = 0 . (50)

n | ε ˜ n (WM) | ε ˜ n (FM) | 2 k ρ 2 ( ε ˜ n − ε ˜ 0 ) (FM) | Ref. [ |
---|---|---|---|---|

0 | −643.845703 | −643.84571 | 0 | 0 |

1 | −643.83502 | −643.83499 | 0.742728866 | 0.83 |

2 | −630.4740 | −630.47403 | 926.9789526 | 935 |

3 | −630.041077 | −630.04109 | 956.9925447 | 961 |

N | ε ˜ n (WM) | ε ˜ n (FM) | 2 k ρ 2 ( ε ˜ n − ε ˜ 0 ) (FM) | Ref. [ |
---|---|---|---|---|

0 | −1096.739532 | −1096.739503 | 0 | 0 |

1 | −1096.738464 | −1096.738420 | 0.0442 | <0.2 |

2 | −1078.470184 | −1078.470125 | 745.8411 | 746 |

3 | −1078.404724 | −1078.404795 | 748.5082 | 748.5 |

4 | −1062.958557 | −1062.958552 | 1379.096 | 1379 |

5 | −1061.637298 | −1061.637339 | 1433.034 | 1434 |

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 − ∞ < X < + ∞ , whereas the finite square well goes to zero discontinuously. Furthermore, the Gaussian well of the form (49) has been used by Bardsley and Comella [

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 [

n b s t = ⌊ 2 X 0 V 0 / π + 1 2 ⌋ + 1 (51)

where ⌊ ⌋ 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

− ℏ 2 2 m X 0 2 d 2 ψ ( χ ) d χ 2 − [ E + V 0 exp ( − χ 2 ) ] ψ ( χ ) = 0 (52)

by changing variables from X to χ = X / X 0 .

The approximation solution we propose here relies on approximating the Gaussian function − exp ( − ξ 2 ) by a potential of type (1). More precisely, we write

− exp ( − χ 2 ) ≈ ∑ j = 1 K B j sech 2 j ( α χ ) (53)

where α , B 1 , B 2 , ⋯ , 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 [ χ min , χ max ] and evaluate the function − exp ( − χ 2 ) at χ k = χ min + k h (for k = 1 , 2 , ⋯ , p ) where p is the number of mesh points, h is the step size, and χ min and χ max are the minimum and maximum values of the variable χ . We finally use the maple 18 Fit function to calculate the B j coefficients. With χ max = 6 = − χ min and h = 0.01 , we find:

{ 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 (54)

Using these parameters, we plot the variations of f G a p p r o x ( χ ) ≡ ∑ j = 1 8 B j sech 2 j ( α χ ) and f G e x a c t ( χ ) ≡ − exp ( − χ 2 ) with χ . We also consider the difference between the two functions, i.e. Δ f G ( χ ) = f G e x a c t ( χ ) − f G a p p r o x ( χ ) .

It is obvious from ^{−6} and hence very small.

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

ε = 2 m X 0 2 ℏ 2 α 2 E end b j = 2 m X 0 2 ℏ 2 α 2 V 0 B j , j = 1 , 2 , ⋯ , K . (55)

In

ϕ j M ( X ) = α ˜ M ( j − M ) ! ( j + M ) ! P j M ( tanh ( α ˜ X ) )

n | E n d i a g | E n F M (N = 300) | E n F M (N = 500) | E n F M (N = 1000) |
---|---|---|---|---|

0 | −2.156866939625 | −2.156871437809 | −2.156868537378 | −2.156868558806 |

1 | −1.513145173222 | −1.513127087344 | −1.513127087344 | −1.513139193491 |

2 | −0.958808621450 | −0.958794700980 | −0.958794700980 | −0.958794700980 |

3 | −0.504841844741 | −0.504830341476 | −0.504830341476 | −0.504830341476 |

4 | −0.170334877363 | −0.170345579080 | −0.170345579080 | −0.170345579080 |

5 | −0.002383210220 | −0.002215413790 | −0.002247696848 | −0.002312262964 |

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.

We have shown that the application of the Fröbenius method to the class of symmetric hyperbolic potentials of the form V ( X ) = ∑ j = 1 K B j sech 2 j ( α X ) 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 ( X ) 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 ( X ) , 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 ( X ) = − V 0 exp ( − X 2 ) by suitably choosing the values of α , K and the B j coefficients.

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

Nyengeri, H., Simbizi, R., Girukwishaka, A., Nizigiyimana, R. and Ndenzako, E. (2018) Fröbenius Series Solutions of the Schrödinger Equation with Various Types of Symmetric Hyperbolic Potentials in One Dimension. Open Access Library Journal, 5: e4728. https://doi.org/10.4236/oalib.1104728