# Unified Treatment of a Class of Spherically Symmetric Potentials: Quasi-Exact Solution

Authors: H. Panahi; M. Baradaran
Journal: Advances in High Energy Physics
 (2016)

---

## Abstract

We investigate the Schrödinger equation for a class of spherically symmetric potentials in a simple and unified manner using the Lie algebraic approach within the framework of quasi-exact solvability. We illustrate that all models give rise to the same basic differential equation, which is expressible as an element of the universal enveloping algebra ofs
l
(
2
). Then, we obtain the general exact solutions of the problem by employing the representation theory of s
l
(
2
) Lie algebra.


---

## Body

## 1. Introduction

From the viewpoint of solvability, the spectral problems are divided into two main classes: the exactly solvable (ES) models and exactly nonsolvable models. A quantum model is called ES if, for all its energy levels and corresponding wavefunctions, explicit expressions can be determined algebraically. These models are distinguished by the fact that there is a natural basis in the Hilbert space in which the infinite-dimensional Hamiltonian can be diagonalized with the help of algebraic methods. In the literature, considerable efforts have been devoted to obtaining the exact solutions of the relativistic and nonrelativistic equations using different methods and techniques [1–5]. In contrast, the exactly nonsolvable models are the spectral problems whose infinite-dimensional Hamiltonian cannot be diagonalized algebraically. Unfortunately, there are only a small number of potentials for which the Schrödinger equation can be solved exactly, like the harmonic oscillator [6, 7], Pöschl-Teller [8, 9], Coulomb [10–12], Morse [13, 14], Rosen-Morse [15], Manning-Rosen [16, 17], Tietz [18], and so forth [19–22]. In 1980s, the series of papers by Shifman, Ushveridze, and Turbiner was devoted to the introduction of an intermediate class between the ES and the exactly nonsolvable models for which a certain finite number of eigenvalues and eigenfunctions, but not the whole spectrum, can be calculated exactly by algebraic methods. They were called quasi-exactly solvable (QES) [23–26]. These models are distinguished by the fact that the Hamiltonian is expressible as a quadratic combination of the generators of a finite-dimensional Lie algebra of first-order differential operators preserving a finite-dimensional subspace of functions and thereby can be represented as a block-diagonal matrix with at least one finite block. Thus, the problem reduces to diagonalizing this block and computing the corresponding eigenvalues and eigenfunctions, which always can be done. In this paper, using the Lie algebraic approach, we present a simple unified derivation and exact solution of the Schrödinger equation for a class of four spherically symmetric potentials within the framework of quasi-exact solvability. We demonstrate that all four cases are reducible to the same basic differential equation which can be solved exactly due to the existence of a hidden s
l
(
2
) symmetry. Then, with the aid of the representation theory of s
l
(
2
), the general exact solution to the basic equation is determined.This paper is organized as follows: in Section2, we briefly review the Lie algebraic approach of quasi-exact solvability. In Section 3, we introduce the four systems and transform the corresponding equations into the same form that is suitable for s
l
(
2
) algebraization and can be expressed as an element of the universal enveloping algebra of s
l
(
2
). Then, using the representation theory of s
l
(
2
), we obtain the general exact solution to the basic equation in Section 4. Also, the closed-form expressions for the eigenvalues and eigenfunctions as well as the allowed parameters of potentials are given for each of the systems. We end with conclusions in Section 5.


## 2. Quasi-Exact Solvability through the sl(2) Algebraization

The general problem of quantum mechanics is to solve the Schrödinger equationH
Ψ
=
E
Ψ, where the wavefunction Ψ belongs to the space of square integrable functions L
2
(
R
). An eigenvalue equation is called QES if there is a nontrivial finite-dimensional subspace U of L
2
(
R
) which is invariant under H: that is, H
U
⊆
U [25, 26]. According to [25], any one-dimensional QES differential equation possesses a hidden Lie algebra s
l
(
2
) which is the only Lie algebra of first-order differential operators that possesses a finite-dimensional representation. This implies that they can be rewritten in terms of s
l
(
2
) generators with differential operators [23–25]:(1)
J
n
+
=
-
x
2
d
d
x
+
n
x
,
J
n
0
=
x
d
d
x
-
n
2
,
J
n
-
=
d
d
x
,which obey the commutation relations (2)
J
n
+
,
J
n
-
=
2
J
n
0
,
J
n
0
,
J
n
±
=
±
J
n
±
,and leave invariant the (n
+
1)-dimensional vector space of polynomials:(3)
P
n
+
1
x
=
1
,
x
,
x
2
,
…
,
x
n
.With these properties, it is easy to verify that the most general second-order differential operator in the enveloping algebra of s
l
(
2
) takes the form(4)
H
=
∑
a
,
b
=
0
,
±
C
a
b
J
a
J
b
+
∑
a
=
0
,
±
C
a
J
a
+
C
,where C
a
b, C
a, and C are real constants. On the other hand, this operator as an ordinary differential equation has the general form(5)
H
ϕ
x
=
0
,
H
=
-
P
4
x
d
2
d
x
2
+
P
3
x
d
d
x
+
P
2
x
,where P
j are the polynomials of degree j. This operator can be turned into the Schrödinger-like operator:(6)
H
~
=
e
-
A
z
H
e
A
z
=
-
1
2
d
2
d
z
2
+
A
′
2
-
A
′
′
+
P
2
x
z
,through the following change of variable and gauge transformation:(7)
z
=
±
∫
d
x
P
4
,
ϕ
x
=
e
-
∫
P
3
/
P
4
d
x
+
log
⁡
z
′
ψ
z
.In the next section, using the method given above, we show that exact solutions of the Schrödinger equation for the four models can be simply obtained in a unified treatment.


## 3. The Four Models and the Corresponding Differential Equations

In this section, we introduce the models which will be the object of study in this paper. For each case, we illustrate that the corresponding Schrödinger equation is reducible to the basic differential equation of second order which is QES due to the existence of a hiddens
l
(
2
) algebraic structure.


### 3.1. Nonpolynomial Potential

First, we consider the nonpolynomial oscillator defined as [27](8)
V
r
=
r
2
+
α
r
2
1
+
β
r
2
,where β
>
0 and α is a real constant. This potential appears in various branches of physics such as the zero-dimensional quantum field theory with nonlinear Lagrangian [28, 29], quantum mechanics [30, 31], laser physics [32, 33], and so forth. This potential has been studied by a variety of methods including the analytic continued fractions [34], the supersymmetric quantum mechanics [35], the 1
/
N expansion method [27], the wavefunction ansatz method [12], and the Bethe ansatz method [36]. In atomic units (m
=
ħ
=
c
=
1), the radial Schrödinger equation with potential (8) is (9)
-
1
2
d
2
d
r
2
+
l
l
+
1
2
r
2
+
r
2
+
α
r
2
1
+
β
r
2
-
E
ψ
r
=
0
.Using the change of variable z
=
-
β
r
2 and making the gauge transformation(10)
ψ
z
=
1
-
z
z
l
+
1
/
2
e
-
z
/
β
2
φ
z
,which preserves the asymptotic behaviour of the wavefunction at the origin and infinity, (9) becomes(11)
z
z
-
1
d
2
d
z
2
+
-
2
β
2
z
2
+
2
β
2
+
l
+
7
2
z
-
l
+
3
2
d
d
z
+
-
E
2
β
+
α
2
β
2
-
1
β
2
l
+
7
2
z
+
E
2
β
+
l
2
+
1
+
1
β
2
l
+
3
2
φ
z
=
0
.


### 3.2. Screened Coulomb Potential

Here, we consider the screened Coulomb potential defined by [37](12)
V
r
=
γ
r
+
δ
r
+
κ
,
γ
<
-
δ
.The corresponding radial Schrödinger equation is given by(13)
-
1
2
d
2
d
r
2
+
l
l
+
1
2
r
2
+
γ
r
+
δ
r
+
κ
ψ
r
=
E
ψ
r
.Several methods and techniques for solving this problem can be found in [12, 36–39]. Applying the transformation(14)
ψ
r
=
r
+
κ
r
l
+
1
e
-
-
2
E
r
+
κ
φ
rand also replacing the variable r by z, we obtain(15)
z
z
+
κ
d
2
d
z
2
+
2
-
-
2
E
z
2
+
-
-
2
E
κ
+
l
+
2
z
+
l
+
2
κ
d
d
z
+
2
-
δ
-
γ
-
-
2
E
l
+
2
z
-
2
κ
l
+
1
-
2
E
+
κ
γ
-
l
-
1
φ
z
=
0
.


### 3.3. Singular Integer Power Potential

The problem of the singular power potentials has been widely carried out in various branches of physics, in both classical and quantum mechanics [40–43]. Here, we consider the singular integer power potential as [44, 45](16)
V
r
=
λ
r
+
μ
r
2
+
ξ
r
3
+
τ
r
4
,with the corresponding Schrödinger equation given by (17)
-
1
2
d
2
d
r
2
+
l
l
+
1
2
r
2
+
λ
r
+
μ
r
2
+
ξ
r
3
+
τ
r
4
ψ
r
=
E
ψ
r
.From the asymptotic behaviour of the wavefunction, we consider the following transformation:(18)
ψ
r
=
r
1
+
ξ
/
2
τ
e
-
-
2
E
r
+
2
τ
/
r
φ
r
.Substituting this into (17) and also replacing r by z, we get(19)
z
2
d
2
d
z
2
+
2
-
-
2
E
z
2
+
1
+
ξ
2
τ
z
+
2
τ
d
d
z
-
2
-
2
E
1
+
ξ
2
τ
+
λ
z
-
2
μ
-
l
l
+
1
-
2
-
4
τ
E
+
ξ
2
τ
+
ξ
2
2
τ
φ
z
=
0
.


### 3.4. Singular Anharmonic Potential

Here, following [46], we consider the potential(20)
V
r
=
ω
r
2
+
ε
r
2
+
σ
r
4
+
χ
r
6
,with the radial Schrödinger equation: (21)
-
1
2
d
2
d
r
2
+
l
l
+
1
2
r
2
+
ω
r
2
+
ε
r
2
+
σ
r
4
+
χ
r
6
ψ
r
=
E
ψ
r
.Similar to the previous cases, we extract the asymptotic behaviour of the wavefunction by making the following transformations:(22)
z
=
r
2
,
ψ
r
=
r
3
/
2
+
σ
/
2
χ
e
-
ω
/
2
r
2
+
χ
/
2
1
/
r
2
φ
r
.After substituting this into (21), we obtain(23)
z
2
d
2
d
z
2
+
-
2
ω
z
2
+
2
+
σ
2
χ
z
+
2
χ
d
d
z
+
E
2
-
ω
2
2
+
σ
2
χ
z
-
1
4
l
′
l
′
+
1
+
2
4
ω
χ
-
σ
2
2
χ
-
2
σ
2
χ
-
3
4
φ
z
=
0
,where(24)
l
′
=
-
1
+
4
l
2
+
4
l
+
8
ε
+
1
2
.


## 3.1. Nonpolynomial Potential

First, we consider the nonpolynomial oscillator defined as [27](8)
V
r
=
r
2
+
α
r
2
1
+
β
r
2
,where β
>
0 and α is a real constant. This potential appears in various branches of physics such as the zero-dimensional quantum field theory with nonlinear Lagrangian [28, 29], quantum mechanics [30, 31], laser physics [32, 33], and so forth. This potential has been studied by a variety of methods including the analytic continued fractions [34], the supersymmetric quantum mechanics [35], the 1
/
N expansion method [27], the wavefunction ansatz method [12], and the Bethe ansatz method [36]. In atomic units (m
=
ħ
=
c
=
1), the radial Schrödinger equation with potential (8) is (9)
-
1
2
d
2
d
r
2
+
l
l
+
1
2
r
2
+
r
2
+
α
r
2
1
+
β
r
2
-
E
ψ
r
=
0
.Using the change of variable z
=
-
β
r
2 and making the gauge transformation(10)
ψ
z
=
1
-
z
z
l
+
1
/
2
e
-
z
/
β
2
φ
z
,which preserves the asymptotic behaviour of the wavefunction at the origin and infinity, (9) becomes(11)
z
z
-
1
d
2
d
z
2
+
-
2
β
2
z
2
+
2
β
2
+
l
+
7
2
z
-
l
+
3
2
d
d
z
+
-
E
2
β
+
α
2
β
2
-
1
β
2
l
+
7
2
z
+
E
2
β
+
l
2
+
1
+
1
β
2
l
+
3
2
φ
z
=
0
.


## 3.2. Screened Coulomb Potential

Here, we consider the screened Coulomb potential defined by [37](12)
V
r
=
γ
r
+
δ
r
+
κ
,
γ
<
-
δ
.The corresponding radial Schrödinger equation is given by(13)
-
1
2
d
2
d
r
2
+
l
l
+
1
2
r
2
+
γ
r
+
δ
r
+
κ
ψ
r
=
E
ψ
r
.Several methods and techniques for solving this problem can be found in [12, 36–39]. Applying the transformation(14)
ψ
r
=
r
+
κ
r
l
+
1
e
-
-
2
E
r
+
κ
φ
rand also replacing the variable r by z, we obtain(15)
z
z
+
κ
d
2
d
z
2
+
2
-
-
2
E
z
2
+
-
-
2
E
κ
+
l
+
2
z
+
l
+
2
κ
d
d
z
+
2
-
δ
-
γ
-
-
2
E
l
+
2
z
-
2
κ
l
+
1
-
2
E
+
κ
γ
-
l
-
1
φ
z
=
0
.


## 3.3. Singular Integer Power Potential

The problem of the singular power potentials has been widely carried out in various branches of physics, in both classical and quantum mechanics [40–43]. Here, we consider the singular integer power potential as [44, 45](16)
V
r
=
λ
r
+
μ
r
2
+
ξ
r
3
+
τ
r
4
,with the corresponding Schrödinger equation given by (17)
-
1
2
d
2
d
r
2
+
l
l
+
1
2
r
2
+
λ
r
+
μ
r
2
+
ξ
r
3
+
τ
r
4
ψ
r
=
E
ψ
r
.From the asymptotic behaviour of the wavefunction, we consider the following transformation:(18)
ψ
r
=
r
1
+
ξ
/
2
τ
e
-
-
2
E
r
+
2
τ
/
r
φ
r
.Substituting this into (17) and also replacing r by z, we get(19)
z
2
d
2
d
z
2
+
2
-
-
2
E
z
2
+
1
+
ξ
2
τ
z
+
2
τ
d
d
z
-
2
-
2
E
1
+
ξ
2
τ
+
λ
z
-
2
μ
-
l
l
+
1
-
2
-
4
τ
E
+
ξ
2
τ
+
ξ
2
2
τ
φ
z
=
0
.


## 3.4. Singular Anharmonic Potential

Here, following [46], we consider the potential(20)
V
r
=
ω
r
2
+
ε
r
2
+
σ
r
4
+
χ
r
6
,with the radial Schrödinger equation: (21)
-
1
2
d
2
d
r
2
+
l
l
+
1
2
r
2
+
ω
r
2
+
ε
r
2
+
σ
r
4
+
χ
r
6
ψ
r
=
E
ψ
r
.Similar to the previous cases, we extract the asymptotic behaviour of the wavefunction by making the following transformations:(22)
z
=
r
2
,
ψ
r
=
r
3
/
2
+
σ
/
2
χ
e
-
ω
/
2
r
2
+
χ
/
2
1
/
r
2
φ
r
.After substituting this into (21), we obtain(23)
z
2
d
2
d
z
2
+
-
2
ω
z
2
+
2
+
σ
2
χ
z
+
2
χ
d
d
z
+
E
2
-
ω
2
2
+
σ
2
χ
z
-
1
4
l
′
l
′
+
1
+
2
4
ω
χ
-
σ
2
2
χ
-
2
σ
2
χ
-
3
4
φ
z
=
0
,where(24)
l
′
=
-
1
+
4
l
2
+
4
l
+
8
ε
+
1
2
.


## 4. Solutions of the Basic Differential Equation for the Four Models

In the previous section, we have shown that our QES models, after the appropriate transformations, are expressible as second-order differential equations (11), (15), (19), and (23), respectively. These equations have the same basic structure:(25)
H
ϕ
z
=
0
,
H
=
z
z
-
a
d
2
d
z
2
+
b
2
z
2
+
b
1
z
+
b
0
d
d
z
+
c
1
z
+
c
0
,where a, c
0, c
1, b
0, b
1, and b
2 are real constants. Here, we intend to solve this equation using the Lie algebraic approach within the representation theory of s
l
(
2
). More precisely, from (4), the general form of a one-dimensional QES differential equation is as follows [25]:(26)
H
=
C
+
+
J
n
+
J
n
+
+
C
+
0
J
n
+
J
n
0
+
C
+
-
J
n
+
J
n
-
+
C
0
-
J
n
0
J
n
-
+
C
-
-
J
n
-
J
n
-
+
C
+
J
n
+
+
C
0
J
n
0
+
C
-
J
n
-
+
C
,which clearly preserves the (n
+
1)-dimensional representation space of the s
l
(
2
) algebra as(27)
ϕ
n
z
=
∑
m
=
0
n
p
m
z
m
,
n
=
0,1
,
2
,
…
.Comparing (25) with (26), it is seen that the differential operator H can be expressed as a special case of the general form (26) as(28)
H
=
-
J
n
+
J
n
-
-
a
J
n
0
J
n
-
-
b
2
J
n
+
+
n
+
b
1
J
n
0
+
b
0
-
n
2
a
J
n
-
+
n
2
2
+
n
2
b
1
+
c
0
,if the following constraint on the coefficients holds:(29)
c
1
=
-
n
b
2
.Hence, we have shown that the differential operator H is an element of the universal enveloping algebra of s
l
(
2
) and thereby we can use the representation theory of s
l
(
2
) to determine the solutions of the problem. In the (n
+
1)-dimensional space ϕ
n
(
z
), the operators J
n
+, J
n
0, and J
n
- can be represented by the n
+
1
×
n
+
1 matrices in the basis 1
,
z
,
z
2
,
…
,
z
n
+
1: (30)
J
n
+
=
0
1
0
⋯
0
0
0
0
2
⋯
0
0
⋮
⋮
⋯
3
⋮
⋮
⋮
⋮
⋯
⋯
n
-
1
0
0
0
⋯
⋯
0
n
0
0
⋯
⋯
0
0
,
J
n
-
=
0
0
0
⋯
0
0
n
0
0
⋯
0
0
⋮
n
-
1
⋯
⋯
⋮
⋮
⋮
⋮
⋯
⋯
0
0
0
0
⋯
2
0
0
0
0
⋯
⋯
1
0
,
J
n
0
=
n
2
0
⋯
⋯
0
0
0
n
2
-
1
⋯
⋯
0
0
⋮
⋮
⋯
⋯
⋮
⋮
0
0
⋯
⋯
1
-
n
2
0
0
0
⋯
⋯
0
-
n
2
.Substituting (30) into (28) results in a matrix equation whose nontrivial solution exists if the following condition is fulfilled (Cramer’s rule):(31)
c
0
b
0
0
0
0
-
n
b
2
c
0
+
b
1
2
b
0
-
2
a
0
0
0
-
n
-
1
b
2
c
0
+
2
b
1
+
2
⋱
0
0
⋱
⋱
⋱
⋮
0
0
0
0
-
2
b
2
0
c
0
+
n
-
1
b
1
+
n
-
2
-
b
2
n
b
0
-
n
-
1
a
c
0
+
n
b
1
+
n
-
1
=
0
.This equation yields a (n
+
1)-set of important constraints on the potential parameters that vary for different n. These constraints together with (29) yield the exact solutions of the systems. It can be easily shown that the expansion coefficients p
m’s in (27) satisfy the three-term recursion relation:(32)
p
m
+
1
=
2
b
2
p
m
-
1
-
c
0
+
m
b
1
+
m
-
1
p
m
m
+
1
b
0
-
m
a
,with the boundary conditions p
n
+
1
=
0 and p
-
1
=
0. Therefore, we have succeeded in obtaining the exact expressions for the energies, wavefunctions, and the allowed values of the potential parameters for the n
+
1 first states algebraically. The main advantage of our algebraic method is that we can quickly obtain the general solutions of the systems for any arbitrary n from (27), (29), and (31) without the cumbersome numerical and analytical procedures usually involved in obtaining the solutions for higher states. In the following, we apply the above results to obtain explicit solutions for each of the four systems.


### 4.1. Nonpolynomial Potential

In this case, from (11) with (25), we get (33)
a
=
1
,
b
2
=
-
2
β
2
,
b
1
=
2
β
2
+
l
+
7
2
,
b
0
=
-
l
+
3
2
,
c
1
=
-
E
2
β
+
α
2
β
2
-
1
β
2
l
+
7
2
,
c
0
=
E
2
β
+
l
2
+
1
+
1
β
2
l
+
3
2
.Then by (29), the closed form of the energy of the system is obtained as(34)
E
n
=
α
β
-
2
2
n
+
l
+
7
2
,which together with (31) yields the exact solutions of the system. Also, from (10) and (27), the wavefunction of the model is obtained as(35)
ψ
n
z
=
1
-
z
z
l
+
1
/
2
e
-
z
/
β
2
∑
m
=
0
n
p
m
z
m
,
n
=
0,1
,
2
,
…
,where the expansion coefficients p
m’s satisfy the recursion relation:(36)
p
m
+
1
=
-
4
/
β
2
p
m
-
1
-
E
/
2
β
+
l
/
2
+
1
+
1
/
β
2
l
+
3
/
2
+
m
2
/
β
2
+
l
+
7
/
2
+
m
-
1
p
m
-
m
+
1
l
+
3
/
2
+
m
.The results obtained for the first three states of this model are displayed in Table 1.Table 1
Exact solutions of the Schrödinger equation with nonpolynomial potential for the ground, first, and second excited states.
n
Energy
Relation between potential parameters andl
The radial wavefunction
0
E
0
=
α
β
-
2
l
+
7
2
α
2
β
2
-
2
β
+
l
2
+
1
=
0
ψ
0
(
r
)
=
1
+
β
r
2
-
r
β
l
+
1
e
-
(
r
2
/
2
)
p
0
1
E
1
=
α
β
-
2
l
+
11
2
8
l
2
+
46
l
+
53
β
+
4
E
1
l
+
10
E
1
2
+
6
l
2
+
30
l
+
36
β
2
+
8
l
+
22
β
α
β
-
2
l
+
11
2
+
2
E
1
2
+
4
l
2
+
20
l
+
21
=
0
ψ
1
r
=
1
+
β
r
2
-
r
β
l
+
1
e
-
r
2
/
2
p
0
-
p
1
β
r
2,
p
1
=
2
2
l
+
3
E
1
2
β
+
l
2
+
1
+
1
β
2
l
+
3
2
p
0
2
E
2
=
α
β
-
2
l
+
19
2
1
64
2
2
l
+
2
l
β
+
2
E
2
+
3
2
+
4
β
3
+
1
32
6
l
β
+
25
β
+
6
2
2
2
l
+
2
l
β
+
2
E
2
+
3
2
+
4
β
2
+
2
+
4
2
β
l
+
3
+
β
2
l
+
9
2
l
+
7
2
2
l
+
3
2
+
E
2
+
2
+
l
β
+
4
l
+
6
2
β
2
9
2
+
2
β
+
l
-
1
=
0
ψ
2
r
=
1
+
β
r
2
-
r
β
l
+
1
e
-
r
2
/
2
p
0
-
p
1
β
r
2
+
p
2
β
2
r
4,
p
2
=
1
2
l
+
5
4
β
2
p
0
+
E
2
2
β
+
3
l
2
+
1
β
2
l
+
7
2
+
9
2
p
1,
p
1
=
2
2
l
+
3
E
2
2
β
+
l
2
+
1
+
1
β
2
l
+
3
2
p
0


### 4.2. Screened Coulomb Potential

In this case, comparing (15) with (25), we obtain (37)
a
=
-
κ
,
b
2
=
-
2
-
2
E
,
b
1
=
2
-
-
2
E
κ
+
l
+
2
,
b
0
=
2
l
+
4
κ
,
c
1
=
2
-
δ
-
γ
-
-
2
E
l
+
2
,
c
0
=
-
2
κ
l
+
1
-
2
E
+
κ
γ
-
l
-
1
.Then by (29) and (14), we obtain the following relations for the energy eigenvalues and the corresponding wavefunctions:(38)
E
n
=
-
1
2
δ
+
γ
n
+
l
+
2
2
,
ψ
n
r
=
r
+
κ
r
l
+
1
e
-
-
2
E
r
+
κ
∑
m
=
0
n
p
m
r
m
,
n
=
0,1
,
2
,
…
,which together with the determinant relation (31) give the exact solutions of this system. Also, the expansion coefficients p
m’s obey the three-term recursion relation:(39)
p
m
+
1
=
-
4
-
2
E
p
m
-
1
-
-
2
κ
l
+
1
-
2
E
+
κ
γ
-
l
-
1
+
m
2
-
-
2
E
κ
+
l
+
2
+
m
-
1
p
m
m
+
1
2
l
+
4
+
m
κ
.The results for the ground, first, and second excited states of this model are reported in Table 2.Table 2
Exact solutions of the Schrödinger equation with screened Coulomb potential for the ground, first, and second excited states.
n
Energy
Relation between potential parameters andl
The radial wavefunction
0
E
0
=
-
1
2
δ
+
γ
l
+
2
2
κ
l
+
1
-
2
E
0
+
κ
γ
-
l
-
1
=
0
ψ
0
(
r
)
=
(
r
+
κ
)
r
l
+
1
e
-
-
2
E
0
(
r
+
κ
)
p
0
1
E
1
=
-
1
2
δ
+
γ
l
+
3
2
l
+
2
κ
-
2
E
1
+
-
-
2
E
1
κ
+
l
+
2
κ
l
+
1
-
2
E
1
+
κ
γ
-
l
-
1
-
-
κ
l
+
1
-
2
E
1
-
κ
γ
+
l
+
1
2
=
0
ψ
1
(
r
)
=
(
r
+
κ
)
r
l
+
1
e
-
-
2
E
1
(
r
+
κ
)
(
p
0
+
p
1
r
),
p
1
=
2
(
κ
(
l
+
1
)
-
2
E
1
+
κ
γ
-
l
-
1
)
2
l
+
4
κ
p
0
2
E
2
=
-
1
2
δ
+
γ
l
+
4
2
-
κ
l
+
1
-
2
E
2
-
κ
γ
+
l
+
1
3
+
3
κ
l
+
1
-
2
E
2
+
κ
γ
-
l
-
1
2
-
-
2
E
2
κ
+
l
+
7
3
-
2
κ
l
+
1
-
2
E
2
+
κ
γ
-
l
-
1
-
2
E
2
κ
2
+
-
4
κ
l
-
9
κ
-
2
E
2
+
9
l
2
+
5
+
l
2
-
2
κ
l
+
2
-
2
E
2
5
-
2
-
2
E
2
κ
+
2
l
=
0
ψ
2
(
r
)
=
(
r
+
κ
)
r
l
+
1
e
-
-
2
E
2
(
r
+
κ
)
(
p
0
+
p
1
r
+
p
2
r
2
), ⁢
p
2
=
2
-
2
E
2
p
0
-
(
l
+
2
)
κ
-
2
E
2
+
κ
γ
-
2
l
-
3
p
1
-
2
l
-
5
κ,
p
1
=
2
(
κ
(
l
+
1
)
-
2
E
2
+
κ
γ
-
l
-
1
)
2
l
+
4
κ
p
0


### 4.3. Singular Integer Power Potential

In this case, from (19) and (25), we have(40)
a
=
0
,
b
2
=
-
2
-
2
E
,
b
1
=
2
1
+
ξ
2
τ
,
b
0
=
2
2
τ
,
c
1
=
-
2
-
2
E
1
+
ξ
2
τ
+
λ
,
c
0
=
-
2
μ
-
l
l
+
1
-
2
-
4
τ
E
+
ξ
2
τ
+
ξ
2
2
τ
.Then from (29) and (18), we obtain the following relations for energy and wavefunction: (41)
E
n
=
-
λ
2
2
n
+
ξ
2
τ
+
1
-
2
,
ψ
n
r
=
r
1
+
ξ
/
2
τ
e
-
-
2
E
r
+
2
τ
/
r
∑
m
=
0
n
p
m
r
m
,
n
=
0,1
,
2
,
…
,where the expansion coefficients p
m’s satisfy the following recursion relation:(42)
p
m
+
1
=
-
4
-
2
E
p
m
-
1
-
-
2
μ
-
l
l
+
1
-
2
-
4
τ
E
+
ξ
/
2
τ
+
ξ
2
/
2
τ
+
m
2
1
+
ξ
/
2
τ
+
m
-
1
p
m
2
m
+
1
2
τ
.These relations together with (31) yield the exact solutions of the system. The results determined for the ground, first, and second excited states are displayed in Table 3.Table 3
Exact solutions of the Schrödinger equation with singular integer power potential for the ground, first, and second excited states.
n
Energy
Relation between potential parameters andl
The radial wavefunction
0
E
0
=
-
λ
2
2
ξ
2
τ
+
1
2
-
2
μ
-
l
l
+
1
-
2
-
4
τ
E
0
+
ξ
2
τ
+
ξ
2
2
τ
=
0
ψ
0
(
r
)
=
r
1
+
ξ
/
2
τ
e
-
-
2
E
0
r
+
2
τ
/
r
p
0
1
E
1
=
-
λ
2
2
ξ
2
τ
+
2
2
-
4
2
τ
-
2
E
1
+
1
2
τ
ξ
2
τ
-
2
l
l
+
1
τ
-
8
-
τ
E
1
τ
-
4
τ
μ
+
ξ
2
2
+
ξ
2
τ
+
-
2
μ
-
l
l
+
1
-
4
-
τ
E
1
+
1
2
ξ
2
τ
+
1
2
ξ
2
τ
2
=
0
ψ
1
(
r
)
=
r
1
+
ξ
/
2
τ
e
-
-
2
E
1
r
+
2
τ
/
r
(
p
0
+
p
1
r
),
p
1
=
-
1
2
2
τ
-
2
μ
-
l
(
l
+
1
)
-
2
-
4
τ
E
1
+
ξ
2
τ
+
ξ
2
2
τ
p
0
2
E
2
=
-
λ
2
2
ξ
2
τ
+
3
2
-
2
μ
-
l
l
+
1
-
4
-
τ
E
2
+
ξ
2
2
τ
+
1
2
ξ
2
τ
3
+
3
4
τ
2
ξ
2
τ
+
8
3
ξ
2
τ
-
2
l
l
+
1
τ
-
8
-
τ
E
2
τ
-
4
τ
μ
+
ξ
2
2
+
2
ξ
2
τ
2
+
5
ξ
τ
2
τ
-
16
-
E
2
τ
+
6
τ
ξ
2
τ
-
2
l
l
+
1
τ
-
8
τ
-
τ
E
2
-
4
τ
μ
+
ξ
2
-
16
2
τ
-
2
E
2
3
+
ξ
2
τ
=
0
ψ
2
(
r
)
=
r
1
+
ξ
/
2
τ
e
-
-
2
E
2
r
+
2
τ
/
r
(
p
0
+
p
1
r
+
p
2
r
2
),
p
2
=
-
-
2
E
2
p
0
+
1
/
4
2
μ
+
l
(
l
+
1
)
+
2
-
4
τ
E
2
-
3
ξ
/
2
τ
-
ξ
2
/
2
τ
-
2
p
1
2
τ,
p
1
=
-
1
2
2
τ
-
2
μ
-
l
(
l
+
1
)
-
2
-
4
τ
E
2
+
ξ
2
τ
+
ξ
2
2
τ
p
0


### 4.4. Singular Anharmonic Potential

In this case, from (23) and (25), we get (43)
a
=
0
,
b
2
=
-
2
ω
,
b
1
=
2
+
σ
2
χ
,
b
0
=
2
χ
,
c
1
=
E
2
-
ω
2
2
+
σ
2
χ
,
c
0
=
-
1
4
l
′
l
′
+
1
+
2
4
ω
χ
-
σ
2
2
χ
-
2
σ
2
χ
-
3
4
.Then from (29) and (22), we obtain(44)
E
n
=
2
ω
2
n
+
σ
2
χ
+
2
,
ψ
r
=
r
3
/
2
+
σ
/
2
χ
e
-
ω
/
2
r
2
+
χ
/
2
1
/
r
2
∑
m
=
0
n
p
m
r
m
,
n
=
0,1
,
2
,
…
,for the energy and wavefunction, respectively, which together with (31) yield the exact solutions of the system. Also, the coefficients p
m’s satisfy the recursion relation: (45)
p
m
+
1
=
-
2
2
ω
p
m
-
1
+
1
/
4
l
′
l
′
+
1
+
2
4
ω
χ
-
σ
2
/
2
χ
-
2
σ
/
2
χ
-
3
/
4
-
m
2
+
σ
/
2
χ
+
m
-
1
p
m
m
+
1
2
χ
.Solutions of the ground, first, and second excited states of this system are reported in Table 4.Table 4
Exact solutions of the Schrödinger equation with singular anharmonic potential for the ground, first, and second excited states.
n
Energy
Relation between potential parameters andl
The radial wavefunction
0
E
0
=
2
ω
σ
2
χ
+
2
l
′
(
l
′
+
1
)
+
2
4
ω
χ
-
σ
2
2
χ
-
2
σ
2
χ
-
3
4
=
0
ψ
0
r
=
r
3
/
2
+
σ
/
2
χ
e
-
ω
/
2
r
2
+
χ
/
2
1
/
r
2
p
0,
1
E
1
=
2
ω
σ
2
χ
+
4
-
2
ω
χ
+
2
+
σ
2
χ
-
l
′
l
′
+
1
4
-
ω
χ
+
σ
2
8
χ
+
σ
2
2
χ
+
3
16
+
-
l
′
l
′
+
1
4
-
ω
χ
+
σ
2
8
χ
+
σ
2
2
χ
+
3
16
2
=
0
ψ
1
(
r
)
=
r
3
/
2
+
σ
/
2
χ
e
-
ω
/
2
r
2
+
χ
/
2
1
/
r
2
(
p
0
+
p
1
r
2
),
p
1
=
1
4
2
χ
l
′
(
l
′
+
1
)
+
2
4
ω
χ
-
σ
2
2
χ
-
2
σ
2
χ
-
3
4
p
0
2
E
2
=
2
ω
σ
2
χ
+
6
-
l
′
l
′
+
1
4
-
ω
χ
+
σ
2
8
χ
+
σ
2
2
χ
3
16
3
+
8
+
3
σ
2
χ
-
l
′
l
′
+
1
4
-
ω
χ
+
σ
2
8
χ
+
σ
2
2
χ
+
3
16
2
+
4
-
8
ω
χ
+
2
2
+
σ
2
χ
2
+
σ
2
χ
σ
2
8
χ
+
σ
2
2
χ
+
3
16
-
l
′
l
′
+
1
4
-
ω
χ
-
8
ω
χ
3
+
σ
2
χ
=
0
ψ
2
(
r
)
=
r
3
/
2
+
σ
/
2
χ
e
-
ω
/
2
r
2
+
χ
/
2
(
1
/
r
2
)
p
0
+
p
1
r
2
+
p
2
r
4,
p
2
=
-
8
2
ω
p
0
+
l
′
l
′
+
1
+
2
4
ω
χ
-
σ
2
/
2
χ
-
6
σ
/
2
χ
-
35
/
4
p
1
8
2
χ


## 4.1. Nonpolynomial Potential

In this case, from (11) with (25), we get (33)
a
=
1
,
b
2
=
-
2
β
2
,
b
1
=
2
β
2
+
l
+
7
2
,
b
0
=
-
l
+
3
2
,
c
1
=
-
E
2
β
+
α
2
β
2
-
1
β
2
l
+
7
2
,
c
0
=
E
2
β
+
l
2
+
1
+
1
β
2
l
+
3
2
.Then by (29), the closed form of the energy of the system is obtained as(34)
E
n
=
α
β
-
2
2
n
+
l
+
7
2
,which together with (31) yields the exact solutions of the system. Also, from (10) and (27), the wavefunction of the model is obtained as(35)
ψ
n
z
=
1
-
z
z
l
+
1
/
2
e
-
z
/
β
2
∑
m
=
0
n
p
m
z
m
,
n
=
0,1
,
2
,
…
,where the expansion coefficients p
m’s satisfy the recursion relation:(36)
p
m
+
1
=
-
4
/
β
2
p
m
-
1
-
E
/
2
β
+
l
/
2
+
1
+
1
/
β
2
l
+
3
/
2
+
m
2
/
β
2
+
l
+
7
/
2
+
m
-
1
p
m
-
m
+
1
l
+
3
/
2
+
m
.The results obtained for the first three states of this model are displayed in Table 1.Table 1
Exact solutions of the Schrödinger equation with nonpolynomial potential for the ground, first, and second excited states.
n
Energy
Relation between potential parameters andl
The radial wavefunction
0
E
0
=
α
β
-
2
l
+
7
2
α
2
β
2
-
2
β
+
l
2
+
1
=
0
ψ
0
(
r
)
=
1
+
β
r
2
-
r
β
l
+
1
e
-
(
r
2
/
2
)
p
0
1
E
1
=
α
β
-
2
l
+
11
2
8
l
2
+
46
l
+
53
β
+
4
E
1
l
+
10
E
1
2
+
6
l
2
+
30
l
+
36
β
2
+
8
l
+
22
β
α
β
-
2
l
+
11
2
+
2
E
1
2
+
4
l
2
+
20
l
+
21
=
0
ψ
1
r
=
1
+
β
r
2
-
r
β
l
+
1
e
-
r
2
/
2
p
0
-
p
1
β
r
2,
p
1
=
2
2
l
+
3
E
1
2
β
+
l
2
+
1
+
1
β
2
l
+
3
2
p
0
2
E
2
=
α
β
-
2
l
+
19
2
1
64
2
2
l
+
2
l
β
+
2
E
2
+
3
2
+
4
β
3
+
1
32
6
l
β
+
25
β
+
6
2
2
2
l
+
2
l
β
+
2
E
2
+
3
2
+
4
β
2
+
2
+
4
2
β
l
+
3
+
β
2
l
+
9
2
l
+
7
2
2
l
+
3
2
+
E
2
+
2
+
l
β
+
4
l
+
6
2
β
2
9
2
+
2
β
+
l
-
1
=
0
ψ
2
r
=
1
+
β
r
2
-
r
β
l
+
1
e
-
r
2
/
2
p
0
-
p
1
β
r
2
+
p
2
β
2
r
4,
p
2
=
1
2
l
+
5
4
β
2
p
0
+
E
2
2
β
+
3
l
2
+
1
β
2
l
+
7
2
+
9
2
p
1,
p
1
=
2
2
l
+
3
E
2
2
β
+
l
2
+
1
+
1
β
2
l
+
3
2
p
0


## 4.2. Screened Coulomb Potential

In this case, comparing (15) with (25), we obtain (37)
a
=
-
κ
,
b
2
=
-
2
-
2
E
,
b
1
=
2
-
-
2
E
κ
+
l
+
2
,
b
0
=
2
l
+
4
κ
,
c
1
=
2
-
δ
-
γ
-
-
2
E
l
+
2
,
c
0
=
-
2
κ
l
+
1
-
2
E
+
κ
γ
-
l
-
1
.Then by (29) and (14), we obtain the following relations for the energy eigenvalues and the corresponding wavefunctions:(38)
E
n
=
-
1
2
δ
+
γ
n
+
l
+
2
2
,
ψ
n
r
=
r
+
κ
r
l
+
1
e
-
-
2
E
r
+
κ
∑
m
=
0
n
p
m
r
m
,
n
=
0,1
,
2
,
…
,which together with the determinant relation (31) give the exact solutions of this system. Also, the expansion coefficients p
m’s obey the three-term recursion relation:(39)
p
m
+
1
=
-
4
-
2
E
p
m
-
1
-
-
2
κ
l
+
1
-
2
E
+
κ
γ
-
l
-
1
+
m
2
-
-
2
E
κ
+
l
+
2
+
m
-
1
p
m
m
+
1
2
l
+
4
+
m
κ
.The results for the ground, first, and second excited states of this model are reported in Table 2.Table 2
Exact solutions of the Schrödinger equation with screened Coulomb potential for the ground, first, and second excited states.
n
Energy
Relation between potential parameters andl
The radial wavefunction
0
E
0
=
-
1
2
δ
+
γ
l
+
2
2
κ
l
+
1
-
2
E
0
+
κ
γ
-
l
-
1
=
0
ψ
0
(
r
)
=
(
r
+
κ
)
r
l
+
1
e
-
-
2
E
0
(
r
+
κ
)
p
0
1
E
1
=
-
1
2
δ
+
γ
l
+
3
2
l
+
2
κ
-
2
E
1
+
-
-
2
E
1
κ
+
l
+
2
κ
l
+
1
-
2
E
1
+
κ
γ
-
l
-
1
-
-
κ
l
+
1
-
2
E
1
-
κ
γ
+
l
+
1
2
=
0
ψ
1
(
r
)
=
(
r
+
κ
)
r
l
+
1
e
-
-
2
E
1
(
r
+
κ
)
(
p
0
+
p
1
r
),
p
1
=
2
(
κ
(
l
+
1
)
-
2
E
1
+
κ
γ
-
l
-
1
)
2
l
+
4
κ
p
0
2
E
2
=
-
1
2
δ
+
γ
l
+
4
2
-
κ
l
+
1
-
2
E
2
-
κ
γ
+
l
+
1
3
+
3
κ
l
+
1
-
2
E
2
+
κ
γ
-
l
-
1
2
-
-
2
E
2
κ
+
l
+
7
3
-
2
κ
l
+
1
-
2
E
2
+
κ
γ
-
l
-
1
-
2
E
2
κ
2
+
-
4
κ
l
-
9
κ
-
2
E
2
+
9
l
2
+
5
+
l
2
-
2
κ
l
+
2
-
2
E
2
5
-
2
-
2
E
2
κ
+
2
l
=
0
ψ
2
(
r
)
=
(
r
+
κ
)
r
l
+
1
e
-
-
2
E
2
(
r
+
κ
)
(
p
0
+
p
1
r
+
p
2
r
2
), ⁢
p
2
=
2
-
2
E
2
p
0
-
(
l
+
2
)
κ
-
2
E
2
+
κ
γ
-
2
l
-
3
p
1
-
2
l
-
5
κ,
p
1
=
2
(
κ
(
l
+
1
)
-
2
E
2
+
κ
γ
-
l
-
1
)
2
l
+
4
κ
p
0


## 4.3. Singular Integer Power Potential

In this case, from (19) and (25), we have(40)
a
=
0
,
b
2
=
-
2
-
2
E
,
b
1
=
2
1
+
ξ
2
τ
,
b
0
=
2
2
τ
,
c
1
=
-
2
-
2
E
1
+
ξ
2
τ
+
λ
,
c
0
=
-
2
μ
-
l
l
+
1
-
2
-
4
τ
E
+
ξ
2
τ
+
ξ
2
2
τ
.Then from (29) and (18), we obtain the following relations for energy and wavefunction: (41)
E
n
=
-
λ
2
2
n
+
ξ
2
τ
+
1
-
2
,
ψ
n
r
=
r
1
+
ξ
/
2
τ
e
-
-
2
E
r
+
2
τ
/
r
∑
m
=
0
n
p
m
r
m
,
n
=
0,1
,
2
,
…
,where the expansion coefficients p
m’s satisfy the following recursion relation:(42)
p
m
+
1
=
-
4
-
2
E
p
m
-
1
-
-
2
μ
-
l
l
+
1
-
2
-
4
τ
E
+
ξ
/
2
τ
+
ξ
2
/
2
τ
+
m
2
1
+
ξ
/
2
τ
+
m
-
1
p
m
2
m
+
1
2
τ
.These relations together with (31) yield the exact solutions of the system. The results determined for the ground, first, and second excited states are displayed in Table 3.Table 3
Exact solutions of the Schrödinger equation with singular integer power potential for the ground, first, and second excited states.
n
Energy
Relation between potential parameters andl
The radial wavefunction
0
E
0
=
-
λ
2
2
ξ
2
τ
+
1
2
-
2
μ
-
l
l
+
1
-
2
-
4
τ
E
0
+
ξ
2
τ
+
ξ
2
2
τ
=
0
ψ
0
(
r
)
=
r
1
+
ξ
/
2
τ
e
-
-
2
E
0
r
+
2
τ
/
r
p
0
1
E
1
=
-
λ
2
2
ξ
2
τ
+
2
2
-
4
2
τ
-
2
E
1
+
1
2
τ
ξ
2
τ
-
2
l
l
+
1
τ
-
8
-
τ
E
1
τ
-
4
τ
μ
+
ξ
2
2
+
ξ
2
τ
+
-
2
μ
-
l
l
+
1
-
4
-
τ
E
1
+
1
2
ξ
2
τ
+
1
2
ξ
2
τ
2
=
0
ψ
1
(
r
)
=
r
1
+
ξ
/
2
τ
e
-
-
2
E
1
r
+
2
τ
/
r
(
p
0
+
p
1
r
),
p
1
=
-
1
2
2
τ
-
2
μ
-
l
(
l
+
1
)
-
2
-
4
τ
E
1
+
ξ
2
τ
+
ξ
2
2
τ
p
0
2
E
2
=
-
λ
2
2
ξ
2
τ
+
3
2
-
2
μ
-
l
l
+
1
-
4
-
τ
E
2
+
ξ
2
2
τ
+
1
2
ξ
2
τ
3
+
3
4
τ
2
ξ
2
τ
+
8
3
ξ
2
τ
-
2
l
l
+
1
τ
-
8
-
τ
E
2
τ
-
4
τ
μ
+
ξ
2
2
+
2
ξ
2
τ
2
+
5
ξ
τ
2
τ
-
16
-
E
2
τ
+
6
τ
ξ
2
τ
-
2
l
l
+
1
τ
-
8
τ
-
τ
E
2
-
4
τ
μ
+
ξ
2
-
16
2
τ
-
2
E
2
3
+
ξ
2
τ
=
0
ψ
2
(
r
)
=
r
1
+
ξ
/
2
τ
e
-
-
2
E
2
r
+
2
τ
/
r
(
p
0
+
p
1
r
+
p
2
r
2
),
p
2
=
-
-
2
E
2
p
0
+
1
/
4
2
μ
+
l
(
l
+
1
)
+
2
-
4
τ
E
2
-
3
ξ
/
2
τ
-
ξ
2
/
2
τ
-
2
p
1
2
τ,
p
1
=
-
1
2
2
τ
-
2
μ
-
l
(
l
+
1
)
-
2
-
4
τ
E
2
+
ξ
2
τ
+
ξ
2
2
τ
p
0


## 4.4. Singular Anharmonic Potential

In this case, from (23) and (25), we get (43)
a
=
0
,
b
2
=
-
2
ω
,
b
1
=
2
+
σ
2
χ
,
b
0
=
2
χ
,
c
1
=
E
2
-
ω
2
2
+
σ
2
χ
,
c
0
=
-
1
4
l
′
l
′
+
1
+
2
4
ω
χ
-
σ
2
2
χ
-
2
σ
2
χ
-
3
4
.Then from (29) and (22), we obtain(44)
E
n
=
2
ω
2
n
+
σ
2
χ
+
2
,
ψ
r
=
r
3
/
2
+
σ
/
2
χ
e
-
ω
/
2
r
2
+
χ
/
2
1
/
r
2
∑
m
=
0
n
p
m
r
m
,
n
=
0,1
,
2
,
…
,for the energy and wavefunction, respectively, which together with (31) yield the exact solutions of the system. Also, the coefficients p
m’s satisfy the recursion relation: (45)
p
m
+
1
=
-
2
2
ω
p
m
-
1
+
1
/
4
l
′
l
′
+
1
+
2
4
ω
χ
-
σ
2
/
2
χ
-
2
σ
/
2
χ
-
3
/
4
-
m
2
+
σ
/
2
χ
+
m
-
1
p
m
m
+
1
2
χ
.Solutions of the ground, first, and second excited states of this system are reported in Table 4.Table 4
Exact solutions of the Schrödinger equation with singular anharmonic potential for the ground, first, and second excited states.
n
Energy
Relation between potential parameters andl
The radial wavefunction
0
E
0
=
2
ω
σ
2
χ
+
2
l
′
(
l
′
+
1
)
+
2
4
ω
χ
-
σ
2
2
χ
-
2
σ
2
χ
-
3
4
=
0
ψ
0
r
=
r
3
/
2
+
σ
/
2
χ
e
-
ω
/
2
r
2
+
χ
/
2
1
/
r
2
p
0,
1
E
1
=
2
ω
σ
2
χ
+
4
-
2
ω
χ
+
2
+
σ
2
χ
-
l
′
l
′
+
1
4
-
ω
χ
+
σ
2
8
χ
+
σ
2
2
χ
+
3
16
+
-
l
′
l
′
+
1
4
-
ω
χ
+
σ
2
8
χ
+
σ
2
2
χ
+
3
16
2
=
0
ψ
1
(
r
)
=
r
3
/
2
+
σ
/
2
χ
e
-
ω
/
2
r
2
+
χ
/
2
1
/
r
2
(
p
0
+
p
1
r
2
),
p
1
=
1
4
2
χ
l
′
(
l
′
+
1
)
+
2
4
ω
χ
-
σ
2
2
χ
-
2
σ
2
χ
-
3
4
p
0
2
E
2
=
2
ω
σ
2
χ
+
6
-
l
′
l
′
+
1
4
-
ω
χ
+
σ
2
8
χ
+
σ
2
2
χ
3
16
3
+
8
+
3
σ
2
χ
-
l
′
l
′
+
1
4
-
ω
χ
+
σ
2
8
χ
+
σ
2
2
χ
+
3
16
2
+
4
-
8
ω
χ
+
2
2
+
σ
2
χ
2
+
σ
2
χ
σ
2
8
χ
+
σ
2
2
χ
+
3
16
-
l
′
l
′
+
1
4
-
ω
χ
-
8
ω
χ
3
+
σ
2
χ
=
0
ψ
2
(
r
)
=
r
3
/
2
+
σ
/
2
χ
e
-
ω
/
2
r
2
+
χ
/
2
(
1
/
r
2
)
p
0
+
p
1
r
2
+
p
2
r
4,
p
2
=
-
8
2
ω
p
0
+
l
′
l
′
+
1
+
2
4
ω
χ
-
σ
2
/
2
χ
-
6
σ
/
2
χ
-
35
/
4
p
1
8
2
χ


## 5. Conclusions

In this paper, we have studied the Schrödinger equation for a class of spherically symmetric potentials and illustrated that these models can be treated in a simple and unified manner in the Lie algebraic approach. We have shown that all these models give rise to the same basic differential equation, which is expressible as an element of the universal enveloping algebra ofs
l
(
2
). We have then obtained the general exact solutions of the basic equation within the framework of representation theory of s
l
(
2
) Lie algebra. Also, we have reported the explicit expressions for the energy, wavefunction, and the constraint on the potential parameters for each of the systems. The advantage of our algebraic method is that we can quickly obtain the general solutions of the systems for any arbitrary n, without the cumbersome procedures of obtaining the solutions for higher states. This method is found to be computationally much simpler than other methods.



---

*Source: 8710604-2016-07-31.xml*