# Hulthén and Coulomb-Like Potentials as a Tensor Interaction within the Relativistic Symmetries of the Manning-Rosen Potential

Authors: Hadi Tokmehdashi; Ali Akbar Rajabi; Majid Hamzavi
Journal: Advances in High Energy Physics
 (2014)

---

## Abstract

The bound-state solutions of the Dirac equation for the Manning-Rosen potential are presented approximately for arbitrary spin-orbit quantum numberκ with the Hulthén and Coulomb-like potentials as a tensor interaction. The generalized parametric Nikiforov-Uvarov (NU) method is used to obtain energy eigenvalues and corresponding two-component spinors of the two Dirac particles and these are obtained in the closed form by using the framework of the spin symmetry and p-spin symmetry concept. We have also shown that tensor interaction removes degeneracies between spin and p-spin doublets. Some numerical results are also given.


---

## Body

## 1. Introduction

The spin and pseudospin symmetry concepts introduced in nuclear theory [1, 2] have been used to explain the features of deformed nuclei [3] and superdeformation [4] and to establish an effective shell-model coupling scheme [5]. Within the framework of the relativistic mean field theory, Ginocchio [6, 7] has found that a Dirac Hamiltonian with scalar and vector harmonic oscillator potentials in the case of (
V
(
r
)
-
S
(
r
)
=
0
) possesses not only a spin symmetry but also a U
(
3
) symmetry, but a Dirac Hamiltonian in the case of (
V
(
r
)
+
S
(
r
)
=
0
) possesses a pseudospin symmetry and a pseudo-U
(
3
) symmetry. Meng et al. [8] have showed that the pseudospin symmetry is exact under the condition (
d
(
V
(
r
)
+
S
(
r
)
)
/
d
r
=
0
). In addition, Alhaidari et al. [9] have investigated in detail physical interpretation on the three-dimensional Dirac equation in the case of spin symmetry limit (
V
(
r
)
-
S
(
r
)
=
0
) and pseudospin symmetry limit (
V
(
r
)
+
S
(
r
)
=
0
). In recent years, by considering the importance of spin and pseudospin symmetries, some authors have contributed many works in this field. For more review of this, one can read the recent works by Wei and Dong [10–13].The p-spin symmetry refers to a quasidegeneracy of single nucleon doublets with nonrelativistic quantum number(
n
,
l
,
j
=
l
+
1
/
2
) and (
n
-
1
,
l
+
2
,
j
=
l
+
3
/
2
), where n, l, and j are single nucleon radial, orbital, and total angular quantum numbers, respectively [1, 2]. The total angular momentum j
=
l
~
+
s
~ with l
~
=
l
+
1 is a pseudoangular momentum and s
~ is p-spin angular momentum [14–18].Motivated by recent work about the spin symmetry and pseudospin symmetry solutions of the Dirac equation for arbitrary spin-orbit coupling quantum numberκ with the Eckart potential and Pöschl-Teller potential carried out by Jia et al. [19, 20], in this work we attempt to study the spin symmetry and pseudospin symmetry solutions of the Dirac equation for arbitrary quantum number κ with the Manning-Rosen potential:
(1)
V

(
r
)
=
1
K
b
2
[
-
A
e
-
r
/
b
1
-
e
-
r
/
b
+
α
(
α
-
1
)
e
-
2
r
/
b
(
1
-
e
-
r
/
b
)
2
]
,
000000000000000000000000
K
=
2
M
ℏ
2
,
where α and A are two dimensionless parameters and b is the screening parameter which is used for determining the range of potential and has dimension of length [21]. It is shown that this potential remains invariant by mapping α
↔
1
-
α. The s-wave scattering solutions and bound state energy eigenvalues for the Manning-Rosen potential have been calculated by using the function analysis method [22] and path integral approach [23], respectively. Recently, Ikhdair and Sever have also studied the arbitrary l-wave solutions of the D-dimensional Schrödinger equation with the Manning-Rosen potential by using an approximate to the centrifugal term [24].For more review of the Manning-Rosen potential, one can refer to the [25–29].Tensor potentials were introduced into the Dirac equation with the substitutionp
→
→
p
→
-
im
ω
β
·
r
^
U
(
r
) and a spin-orbit coupling is added to the Dirac Hamiltonian [30, 31].Recently, different kinds of potentials (e.g., Coulomb-like, Linear plus Coulomb-like, and Yukawa), as tensor couplings, have been used widely in the studies of nuclear properties. In this regard, see [32–39]. Here we study a tensor potential in the Hulthén (2) and Coulomb-like (3) form as follows:
(2)
U
(
r
)
=
-
U
0
H
e
-
r
/
b
1
-
e
-
r
/
b
,
b
=
δ
-
1
,
U
0
H
=
Z
e
2
δ
,
where U
0
H is the potential depth, δ is the screening range parameter, and b represents the spatial range. If the potential is used for atoms, then U
0
H
=
Z
δ (in the relativistic units ℏ
=
c
=
e
=
1), where Z is identified as the atomic number. The Hulthén potential behaves like the Coulomb potential near the origin (i.e., r
→
0 or r
≪
b) U
C
(
r
)
=
-
Z
e
2
/
r but decreases exponentially in the asymptotic region when r
≫
0, so its capacity for bound states is smaller than the Coulomb potential. The Hulthén potential is widely used for the description of the nucleon-heavy nucleus interactions (see [40] and the references therein). This potential has been applied to a number of areas such as nuclear and particle physics [41–43], atomic physics [44, 45], molecular physics [46, 47], and chemical physics [48].Also, the Coulomb-like potential [49] for the tensor due to a charge Z
a
e interacting with a charge Z
b
e, distributed uniformly over a sphere of radius R
C, is added:
(3)
U
(
r
)
=
-
U
0
C
r
,
U
0
C
=
Z
a
Z
b
e
2
4
π
ε
0
,
r
≥
R
C
,
where R
C
=
7.78 fm is the Coulomb radius and Z
b and Z
b denote the charges of the projectile a and the target nuclei b, respectively.These potentials have been used in several branches of physics and their discrete and continuum states have been studied by a variety of techniques such as the algebraic perturbation calculations which are based upon the dynamical group structure SO(2,1) [50], the formalism of supersymmetric quantum mechanics within the framework of the variational method [51], the supersymmetry and shape invariance property [52], the asymptotic iteration method [53, 54], the NU method [55], and the approach proposed by Biedenharn et al. for the Dirac-Coulomb problem [56, 57]. Otherwise, some authors studied relativistic and nonrelativistic equations with different potentials [58–90]. In this paper, we solve approximately the Dirac equation with the Manning-Rosen potential for the spin-orbit quantum number κ and we obtained the energy eigenvalues and wavefunctions of Dirac particle in the fields of the scalar and vector Manning-Rosen potential including a Hulthén and Coulomb tensor interaction under consideration of spin symmetry and pseudospin symmetry case.The organization of this paper is as follows. In Section2, we briefly introduce the Dirac equation with the Manning-Rosen potential and tensor potentials with arbitrary spin-orbit coupling number κ including tensor interaction under spin and p-spin symmetric limits. Also the generalized parametric NU method is presented in this section. The energy eigenvalue equations and corresponding wave functions for Hulthén and Coulomb tensor interaction in terms of the Jacobi polynomials are obtained in Sections 3 and 4. We also give some remark graphs and numerical results. We end with conclusion in Section 5.


## 2. Dirac Equation Including Tensor Coupling

The Dirac equation for fermionic massive spin-1
/
2 particles moving in attractive scalar S
(
r
), repulsive vector V
(
r
), and tensor U
(
r
) potentials is (in units ℏ
=
c
=
1)
(4)
[
α
→
·
p
→
+
β
(
M
+
S
(
r
)
)
-
i
β
α
→
·
r
^
U
(
r
)
]
ψ
(
r
→
)
=
[
E
-
V
(
r
)
]
ψ
(
r
→
)
,
where E is the relativistic energy of the system, p
→
=
-
i
∇
→ is the three-dimensional momentum operator, and M is the mass of the fermionic particle. Further, α
→ and β are the 4
×
4 Dirac matrices given by
(5)
α
→
=
(
0
σ
→
σ
→
0
)
,
β
=
(
I
0
0
-
I
)
,
where I is 2
×
2 unitary matrix and σ
→ are three-vector spin matrices
(6)
σ
1
=
(
0
1
1
0
)
,
σ
2
=
(
0
-
i
i
0
)
,
σ
3
=
(
1
0
0
-
1
)
.
The total angular momentum operator J
→ and spin-orbit K
=
(
σ
→
·
L
→
+
1
), where L
→ is orbital angular momentum of the spherical nucleons, commute with the Dirac Hamiltonian. The eigenvalues of spin-orbit coupling operator are κ
=
(
j
+
1
/
2
)
>
0 and κ
=
-
(
j
+
1
/
2
)
<
0 for unaligned spin j
=
l
-
1
/
2 and the aligned spin j
=
l
+
1
/
2, respectively. (
H
2
,
K
,
J
2
,
J
z
) can be taken as the complete set of the conservative quantities. Thus, the spinor wave functions can be classified according to their angular momentum j; spin-orbit quantum number κ and the radial quantum number n can be written as follows:
(7)
ψ
n
κ
(
r
→
)
=
(
f
n
κ
(
r
→
)
g
n
κ
(
r
→
)
)
=
(
F
n
κ
(
r
)
r
Y
j
m
l
(
θ
,
φ
)
i
G
n
κ
(
r
)
r
Y
j
m
l
~
(
θ
,
φ
)
)
,
where f
n
κ
(
r
→
) is the upper (large) component and g
n
κ
(
r
→
) is the lower (small) component of the Dirac spinors. Y
j
m
l
(
θ
,
φ
) and Y
j
m
l
~
(
θ
,
φ
) are spin and p-spin spherical harmonics, respectively, and m is the projection of the angular momentum on the z-axis. Substituting (7) into (4) and using the following relations:
(8)
(
σ
→
·
A
→
)
(
σ
→
·
B
→
)
=
A
→
·
B
→
+
i
σ
→
·
(
A
→
×
B
→
)
,
(
σ
→
·
P
→
)
=
σ
→
·
r
^
(
r
^
·
P
→
+
i
σ
→
·
L
→
r
)
together with the following properties:
(9)
(
σ
→
·
L
→
)
Y
j
m
l
~
(
θ
,
ϕ
)
=
(
κ
-
1
)
Y
j
m
l
~
(
θ
,
ϕ
)
,
(
σ
→
·
L
→
)
Y
j
m
l
(
θ
,
ϕ
)
=
-
(
κ
-
1
)
Y
j
m
l
(
θ
,
ϕ
)
,
(
σ
→
·
r
^
)
Y
j
m
l
~
(
θ
,
ϕ
)
=
-
Y
j
m
l
(
θ
,
ϕ
)
,
(
σ
→
·
r
^
)
Y
j
m
l
(
θ
,
ϕ
)
=
-
Y
j
m
l
~
(
θ
,
ϕ
)
,
one obtains two coupled differential equations for upper and lower radial wave functions F
n
κ
(
r
) and G
n
κ
(
r
) as
(10)
(
d
d
r
+
κ
r
-
U
(
r
)
)
F
n
κ
(
r
)
=
(
M
+
E
n
κ
-
Δ
(
r
)
)
G
n
κ
(
r
)
,
(11)
(
d
d
r
-
κ
r
+
U
(
r
)
)
G
n
κ
(
r
)
=
(
M
-
E
n
κ
+
Σ
(
r
)
)
F
n
κ
(
r
)
,
where
(12)
Δ
(
r
)
=
V
(
r
)
-
S
(
r
)
,
Σ
(
r
)
=
V
(
r
)
+
S
(
r
)
.
Eliminating F
n
κ
(
r
) and G
n
κ
(
r
) from (10) and (11), we finally obtain the following two Schrödinger-like differential equations for the upper and lower radial spinor components, respectively:
(13)
[
d
2
d
r
2
-
κ
(
κ
+
1
)
r
2
+
2
κ
r
U
(
r
)
-
d
U
(
r
)
d
r
-
U
2
(
r
)
]
F
n
κ
(
r
)
+
d
Δ
(
r
)
/
d
r
M
+
E
n
κ
-
Δ
(
r
)
(
d
d
r
+
κ
r
-
U
(
r
)
)
F
n
κ
(
r
)
=
[
(
M
+
E
n
κ
-
Δ
(
r
)
)
(
M
-
E
n
κ
+
Σ
(
r
)
)
]
F
n
κ
(
r
)
,
(14)
[
d
2
d
r
2
-
κ
(
κ
-
1
)
r
2
+
2
κ
r
U
(
r
)
+
d
U
(
r
)
d
r
-
U
2
(
r
)
]
G
n
κ
(
r
)
+
d
Σ
(
r
)
/
d
r
M
-
E
n
κ
+
Σ
(
r
)
(
d
d
r
-
κ
r
+
U
(
r
)
)
G
n
κ
(
r
)
=
[
(
M
+
E
n
κ
-
Δ
(
r
)
)
(
M
-
E
n
κ
+
Σ
(
r
)
)
]
G
n
κ
(
r
)
,
where κ
(
κ
-
1
)
=
l
~
(
l
~
+
1
) and κ
(
κ
+
1
)
=
l
(
l
+
1
).The quantum numberκ is related to the quantum numbers for spin symmetry l and p-spin symmetry l
~ as
(15)
κ
=
{
-
(
l
+
1
)
=
-
(
j
+
1
2
)
,
(
s
1
/
2
,
p
3
/
2
,
etc
.
)
j
=
l
+
1
2
,
aligned
spin
(
κ
<
0
)
,
+
l
=
+
(
j
+
1
2
)
,
(
p
1
/
2
,
d
3
/
2
,
etc
.
)
j
=
l
-
1
2
,
unaligned
spin
(
κ
>
0
)
,
and the quasidegenerate doublet structure can be expressed in terms of a p-spin angular momentum s
~
=
1
/
2 and pseudoorbital angular momentum l
~, which can be defined as
(16)
κ
=
{
-
l
~
=
-
(
j
+
1
2
)
,
(
s
1
/
2
,
p
3
/
2
,
etc
.
)
j
=
l
~
-
1
2
,
alinged
p-spin
(
κ
<
0
)
,
+
(
l
~
+
1
)
=
+
(
j
+
1
2
)
,
(
d
3
/
2
,
f
5
/
2
,
etc
.
)
j
=
l
~
+
1
2
,
unaligned
p-spin
(
κ
>
0
)
,
where κ
=
±
1
,
±
2
,
…. For example, (
1
s
1
/
2
,
0
d
3
/
2
) and (
1
p
3
/
2
,
0
f
5
/
2
) can be considered as p-spin doublets.


### 2.1. Spin Symmetry Limit

In this section, we will solve Dirac equation under spin symmetry limit with Manning-Rosen potential and Hulthén potential as a tensor interaction. The exact spin symmetry occurs in Dirac equation whend
Δ
/
d
r
=
d
(
V
V
-
V
S
)
/
d
r
=
0 or Δ
=
C
s
= constant [91–94]. We are taking Σ
(
r
) and tensor U
(
r
) potentials as
(17)
Σ
(
r
)
=
V
0
e
-
δ
r
1
-
e
-
δ
b
+
V
1
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
,
U
(
r
)
=
-
U
0
H
e
-
δ
r
1
-
e
-
δ
r
,
where δ
=
1
/
b, V
0
=
-
A
/
K
b
2, and V
1
=
α
(
α
-
1
)
/
K
b
2.In Figure1, structure of the Manning-Rosen potential for various values of α screening parameter and r is shown.Figure 1
The variation of the Manning-Rosen potential according to various values ofα screening parameter and (δ
=
1
/
b
=
0.025, A
=
2
b, K
=
2).Under this symmetry, from (13), one obtains
(18)
[
d
2
d
r
2
-
κ
(
κ
+
1
)
r
2
-
2
κ
U
0
H
e
-
δ
r
r
(
1
-
e
-
δ
r
)
00
-
U
0
H
δ
(
e
-
δ
r
1
-
e
-
δ
r
+
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
)
00
-
U
0
H
2
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
-
γ
(
V
0
e
-
δ
r
1
-
e
-
δ
b
+
V
1
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
)
00
-
β
2
d
2
d
r
2
]
F
n
κ
(
r
)
=
0
,
where γ
=
M
+
E
n
κ
-
C
s and β
2
=
(
M
-
E
n
κ
)
(
M
+
E
n
κ
-
C
s
). Also, κ
=
l and κ
=
-
l
-
1 for κ
<
0 and κ
>
0, respectively.


### 2.2. p-Spin Symmetry Limit

Within the pseudospin symmetry case,d
[
V
(
r
)
+
S
(
r
)
]
/
d
r
=
d
Σ
(
r
)
/
d
r
=
0 or Σ
(
r
)
=
C
ps
= constant, and p-spin symmetry is exact in the Dirac equation [8, 95–97]. In this part, we consider Δ
(
r
) as
(19)
Δ
(
r
)
=
V
0
e
-
δ
r
1
-
e
-
δ
b
+
V
1
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
and tensor U
(
r
) potential as
(20)
U
(
r
)
=
-
U
0
H
e
-
δ
r
1
-
e
-
δ
r
.
Under this symmetry, from (14), we obtain
(21)
[
d
2
d
r
2
-
κ
(
κ
-
1
)
r
2
-
2
κ
U
0
H
e
-
δ
r
r
(
1
-
e
-
δ
r
)
00
+
U
0
H
δ
(
e
-
δ
r
1
-
e
-
δ
r
+
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
)
00
-
U
0
H
2
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
-
γ
~
(
V
0
e
-
δ
r
1
-
e
-
δ
b
+
V
1
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
)
00
-
β
~
2
d
2
d
r
2
]
G
n
κ
(
r
)
=
0
,
where γ
~
=
E
n
κ
-
M
-
C
ps and β
~
2
=
(
M
+
E
n
κ
)
(
M
-
E
n
κ
+
C
ps
). Also, κ
=
-
l
~ and κ
=
l
~
+
1 for κ
<
0 and κ
>
0, respectively.The Dirac equation for the Manning-Rosen potential given by (18) and (21) cannot be solved analytically. To obtain analytical approximate solutions for the Manning-Rosen type potential with spin-orbit coupling term, we have to use an approximation in order to deal with the spin-orbit coupling term [57] as
(22)
1
r
≈
δ
e
-
δ
r
1
-
e
-
δ
r
.
Or equivalently
(23)
1
r
2
≈
δ
2
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
,
which is only valid for small values of the parameter δ.


### 2.3. Basic Equations of Parametric Nikiforov-Uvarov (NU) Method

The NU method is a powerful mathematical tool that recently has been used widely to solve second-order differential equations [98, 99]. In this section, we briefly describe this method. This technique is based on solving the hypergeometric type second-order differential equations by means of the special orthogonal functions [100]. For a given potential, the Schrödinger or Schrödinger-like equations in spherical coordinates are reduced to a generalized equation of hypergeometric type with an appropriate coordinate transformation r
→
s and then they are solved systematically to find the exact or particular solutions. The main equation which is closely associated with the method is given in the following form [101, 102]:
(24)
ψ
n
′′
(
s
)
+
τ
~
(
s
)
σ
(
s
)
ψ
n
′
(
s
)
+
σ
~
(
s
)
σ
2
(
s
)
ψ
n
(
s
)
=
0
,
where σ
(
s
) and σ
~
(
s
) are polynomials, at most of second degree, and τ
~
(
s
) is a first-degree polynomial. To make the application of the NU method simpler and direct without need to check the validity of solution, we present a shortcut for the method. Hence, firstly we write the general form of the Schrödinger-like equation (24) in a more general form as
(25)
ψ
n
′′
(
s
)
+
(
c
1
-
c
2
s
s
(
1
-
c
3
s
)
)
ψ
n
′
(
s
)
+
(
-
ξ
2
s
2
+
ξ
1
s
-
ξ
0
s
2
(
1
-
c
3
s
)
2
)
ψ
n
(
s
)
=
0
,
satisfying the wave functions
(26)
ψ
n
(
s
)
=
ϕ
(
s
)
y
n
(
s
)
.
Secondly, we compare (27) with its counterpart equation (26) to obtain the following parameter values:
(27)
τ
~
(
s
)
=
c
1
-
c
2
s
,
σ
(
s
)
=
s
(
1
-
c
3
s
)
,
σ
~
(
s
)
=
-
ξ
2
s
2
+
ξ
1
s
-
ξ
0
.
Now, following the NU method [101], we obtain the energy equation [102, 103]
(28)
c
2
n
-
(
2
n
+
1
)
c
5
+
(
2
n
+
1
)
(
c
9
-
c
3
c
8
)
+
n
(
n
-
1
)
c
3
+
c
7
+
2
c
3
c
8
-
2
c
8
c
9
=
0
.
The values for the parametric constants are shown in Table 1.Table 1
The values for the parametric constants of the NU method.
Parametric NU constants
ξ
0
c
1
c
4
=
1
2
(
1
-
c
1
)
{
∫
T
T
s
}
c
7
=
2
c
4
c
5
-
ξ
1
c
10
=
c
1
+
2
c
4
-
2
c
8
-
1
>
-
1
ξ
1
c
2
c
5
=
1
2
(
c
2
-
2
c
3
)
{
∫
T
T
s
}
c
8
=
c
4
2
+
ξ
0
c
11
=
1
-
c
1
-
2
c
4
+
2
c
3
c
9
>
-
1
,
c
3
≠
0
ξ
2
c
3
c
6
=
c
5
2
+
ξ
2
c
9
=
c
3
(
c
7
+
c
3
c
8
)
+
c
6
c
12
=
c
4
-
c
8
>
0
c
13
=
-
c
4
+
1
c
3
(
c
9
-
c
5
)
>
0
,
c
3
≠
0
{
∫
T
T
s
}Wherec
12
>
0, c
13
>
0 and s
∈
[
0
,
1
/
c
3
], c
3
≠
0, the corresponding wave functions are
(29)
ρ
(
s
)
=
s
c
10
(
1
-
c
3
s
)
c
11
,
ϕ
(
s
)
=
s
c
12
(
1
-
c
3
s
)
c
13
,
0000000000000000000
00
00
c
12
>
0
,
c
13
>
0
,
y
n
(
s
)
=
P
n
(
c
10
,
c
11
)
(
1
-
2
c
3
s
)
,
c
10
>
-
1
,
c
11
>
-
1
,
ψ
n
κ
(
s
)
=
N
n
κ
s
c
12
(
1
-
c
3
s
)
c
13
P
n
(
c
10
,
c
11
)
(
1
-
2
c
3
s
)
,
where P
n
(
μ
,
ν
)
(
x
), μ
>
-
1, ν
>
-
1, and x
∈
[
-
1,1
] are Jacobi polynomials.


## 2.1. Spin Symmetry Limit

In this section, we will solve Dirac equation under spin symmetry limit with Manning-Rosen potential and Hulthén potential as a tensor interaction. The exact spin symmetry occurs in Dirac equation whend
Δ
/
d
r
=
d
(
V
V
-
V
S
)
/
d
r
=
0 or Δ
=
C
s
= constant [91–94]. We are taking Σ
(
r
) and tensor U
(
r
) potentials as
(17)
Σ
(
r
)
=
V
0
e
-
δ
r
1
-
e
-
δ
b
+
V
1
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
,
U
(
r
)
=
-
U
0
H
e
-
δ
r
1
-
e
-
δ
r
,
where δ
=
1
/
b, V
0
=
-
A
/
K
b
2, and V
1
=
α
(
α
-
1
)
/
K
b
2.In Figure1, structure of the Manning-Rosen potential for various values of α screening parameter and r is shown.Figure 1
The variation of the Manning-Rosen potential according to various values ofα screening parameter and (δ
=
1
/
b
=
0.025, A
=
2
b, K
=
2).Under this symmetry, from (13), one obtains
(18)
[
d
2
d
r
2
-
κ
(
κ
+
1
)
r
2
-
2
κ
U
0
H
e
-
δ
r
r
(
1
-
e
-
δ
r
)
00
-
U
0
H
δ
(
e
-
δ
r
1
-
e
-
δ
r
+
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
)
00
-
U
0
H
2
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
-
γ
(
V
0
e
-
δ
r
1
-
e
-
δ
b
+
V
1
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
)
00
-
β
2
d
2
d
r
2
]
F
n
κ
(
r
)
=
0
,
where γ
=
M
+
E
n
κ
-
C
s and β
2
=
(
M
-
E
n
κ
)
(
M
+
E
n
κ
-
C
s
). Also, κ
=
l and κ
=
-
l
-
1 for κ
<
0 and κ
>
0, respectively.


## 2.2. p-Spin Symmetry Limit

Within the pseudospin symmetry case,d
[
V
(
r
)
+
S
(
r
)
]
/
d
r
=
d
Σ
(
r
)
/
d
r
=
0 or Σ
(
r
)
=
C
ps
= constant, and p-spin symmetry is exact in the Dirac equation [8, 95–97]. In this part, we consider Δ
(
r
) as
(19)
Δ
(
r
)
=
V
0
e
-
δ
r
1
-
e
-
δ
b
+
V
1
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
and tensor U
(
r
) potential as
(20)
U
(
r
)
=
-
U
0
H
e
-
δ
r
1
-
e
-
δ
r
.
Under this symmetry, from (14), we obtain
(21)
[
d
2
d
r
2
-
κ
(
κ
-
1
)
r
2
-
2
κ
U
0
H
e
-
δ
r
r
(
1
-
e
-
δ
r
)
00
+
U
0
H
δ
(
e
-
δ
r
1
-
e
-
δ
r
+
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
)
00
-
U
0
H
2
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
-
γ
~
(
V
0
e
-
δ
r
1
-
e
-
δ
b
+
V
1
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
)
00
-
β
~
2
d
2
d
r
2
]
G
n
κ
(
r
)
=
0
,
where γ
~
=
E
n
κ
-
M
-
C
ps and β
~
2
=
(
M
+
E
n
κ
)
(
M
-
E
n
κ
+
C
ps
). Also, κ
=
-
l
~ and κ
=
l
~
+
1 for κ
<
0 and κ
>
0, respectively.The Dirac equation for the Manning-Rosen potential given by (18) and (21) cannot be solved analytically. To obtain analytical approximate solutions for the Manning-Rosen type potential with spin-orbit coupling term, we have to use an approximation in order to deal with the spin-orbit coupling term [57] as
(22)
1
r
≈
δ
e
-
δ
r
1
-
e
-
δ
r
.
Or equivalently
(23)
1
r
2
≈
δ
2
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
,
which is only valid for small values of the parameter δ.


## 2.3. Basic Equations of Parametric Nikiforov-Uvarov (NU) Method

The NU method is a powerful mathematical tool that recently has been used widely to solve second-order differential equations [98, 99]. In this section, we briefly describe this method. This technique is based on solving the hypergeometric type second-order differential equations by means of the special orthogonal functions [100]. For a given potential, the Schrödinger or Schrödinger-like equations in spherical coordinates are reduced to a generalized equation of hypergeometric type with an appropriate coordinate transformation r
→
s and then they are solved systematically to find the exact or particular solutions. The main equation which is closely associated with the method is given in the following form [101, 102]:
(24)
ψ
n
′′
(
s
)
+
τ
~
(
s
)
σ
(
s
)
ψ
n
′
(
s
)
+
σ
~
(
s
)
σ
2
(
s
)
ψ
n
(
s
)
=
0
,
where σ
(
s
) and σ
~
(
s
) are polynomials, at most of second degree, and τ
~
(
s
) is a first-degree polynomial. To make the application of the NU method simpler and direct without need to check the validity of solution, we present a shortcut for the method. Hence, firstly we write the general form of the Schrödinger-like equation (24) in a more general form as
(25)
ψ
n
′′
(
s
)
+
(
c
1
-
c
2
s
s
(
1
-
c
3
s
)
)
ψ
n
′
(
s
)
+
(
-
ξ
2
s
2
+
ξ
1
s
-
ξ
0
s
2
(
1
-
c
3
s
)
2
)
ψ
n
(
s
)
=
0
,
satisfying the wave functions
(26)
ψ
n
(
s
)
=
ϕ
(
s
)
y
n
(
s
)
.
Secondly, we compare (27) with its counterpart equation (26) to obtain the following parameter values:
(27)
τ
~
(
s
)
=
c
1
-
c
2
s
,
σ
(
s
)
=
s
(
1
-
c
3
s
)
,
σ
~
(
s
)
=
-
ξ
2
s
2
+
ξ
1
s
-
ξ
0
.
Now, following the NU method [101], we obtain the energy equation [102, 103]
(28)
c
2
n
-
(
2
n
+
1
)
c
5
+
(
2
n
+
1
)
(
c
9
-
c
3
c
8
)
+
n
(
n
-
1
)
c
3
+
c
7
+
2
c
3
c
8
-
2
c
8
c
9
=
0
.
The values for the parametric constants are shown in Table 1.Table 1
The values for the parametric constants of the NU method.
Parametric NU constants
ξ
0
c
1
c
4
=
1
2
(
1
-
c
1
)
{
∫
T
T
s
}
c
7
=
2
c
4
c
5
-
ξ
1
c
10
=
c
1
+
2
c
4
-
2
c
8
-
1
>
-
1
ξ
1
c
2
c
5
=
1
2
(
c
2
-
2
c
3
)
{
∫
T
T
s
}
c
8
=
c
4
2
+
ξ
0
c
11
=
1
-
c
1
-
2
c
4
+
2
c
3
c
9
>
-
1
,
c
3
≠
0
ξ
2
c
3
c
6
=
c
5
2
+
ξ
2
c
9
=
c
3
(
c
7
+
c
3
c
8
)
+
c
6
c
12
=
c
4
-
c
8
>
0
c
13
=
-
c
4
+
1
c
3
(
c
9
-
c
5
)
>
0
,
c
3
≠
0
{
∫
T
T
s
}Wherec
12
>
0, c
13
>
0 and s
∈
[
0
,
1
/
c
3
], c
3
≠
0, the corresponding wave functions are
(29)
ρ
(
s
)
=
s
c
10
(
1
-
c
3
s
)
c
11
,
ϕ
(
s
)
=
s
c
12
(
1
-
c
3
s
)
c
13
,
0000000000000000000
00
00
c
12
>
0
,
c
13
>
0
,
y
n
(
s
)
=
P
n
(
c
10
,
c
11
)
(
1
-
2
c
3
s
)
,
c
10
>
-
1
,
c
11
>
-
1
,
ψ
n
κ
(
s
)
=
N
n
κ
s
c
12
(
1
-
c
3
s
)
c
13
P
n
(
c
10
,
c
11
)
(
1
-
2
c
3
s
)
,
where P
n
(
μ
,
ν
)
(
x
), μ
>
-
1, ν
>
-
1, and x
∈
[
-
1,1
] are Jacobi polynomials.


## 3. Bound States of the Manning-Rosen Potential with Hulthén Tensor Interaction

### 3.1. Spin Symmetry Solution

Substituting (22) and (23) into (18) gives
(30)
[
d
2
d
r
2
-
{
κ
(
κ
+
1
)
δ
2
+
(
2
κ
+
1
)
U
0
H
δ
+
U
0
H
2
+
γ
V
1
δ
2
}
00
×
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
+
{
-
U
0
H
δ
-
γ
V
0
}
00
×
e
-
δ
r
1
-
e
-
δ
b
-
β
2
]
F
n
κ
(
r
)
=
0
.
Defining a new variable of the form s
=
e
-
δ
r and using (30), we obtain
(31)
d
2
F
n
κ
(
s
)
d
s
2
+
1
-
s
s
(
1
-
s
)
d
F
n
κ
(
s
)
d
s
+
1
δ
2
s
2
(
1
-
s
)
2
×
[
-
(
A
+
B
+
β
2
)
s
2
+
(
B
+
2
β
2
)
s
-
β
2
]
×
F
n
κ
(
s
)
=
0
,
where
(32)
A
=
κ
(
κ
+
1
)
δ
2
+
(
2
κ
+
1
)
U
0
H
δ
+
U
0
H
2
+
γ
V
1
,
B
=
-
U
0
H
δ
-
γ
V
0
.
Comparing (31) and (25), we can easily obtain the coefficients c
i (i
=
1,2
,
3) and analytical expressions ξ
j (j
=
0,1
,
2). By using Table 1, we are obtained the coefficients c
i (i
=
4,5
,
.
.
.
,
13). The values of all coefficients are shown in Table 2.Table 2
The specific values for the parametric constants necessary for calculating the energy eigenvalues and eigenfunctions of the spin symmetry limit.
Constant
Analytic value
Constant
Analytic value
c
1
c
1
1
c
2
1
c
3
1
c
4
0
c
5
-
1
2
c
6
1
4
+
A
+
B
+
β
2
δ
2
c
7
-
B
+
2
β
2
δ
2
c
8
β
2
δ
2
c
9
1
4
+
A
δ
2
c
10
-
2
β
2
δ
2
c
11
2
η
κ
c
12
-
β
2
δ
2
c
13
η
κ
+
1
2
ξ
0
β
2
δ
2
ξ
1
B
+
2
β
2
δ
2
ξ
2
A
+
B
+
β
2
δ
2In Table2, η
κ
=
(
(
κ
+
1
/
2
)
2
+
(
U
0
H
/
δ
)
(
(
U
0
H
/
δ
)
+
2
κ
+
1
)
+
γ
V
1
/
δ
2
)
1
/
2, and by using (28), we can obtain the energy eigenvalues of the Manning-Rosen potential with spin symmetric limit as
(33)
(
β
2
δ
2
n
+
1
2
+
(
κ
+
1
2
)
2
+
U
0
H
δ
(
U
0
H
δ
+
2
κ
+
1
)
+
γ
V
1
δ
2
-
β
2
δ
2
)
2
=
γ
(
V
1
-
V
0
)
+
β
2
δ
2
+
(
κ
+
1
2
)
2
+
U
0
H
δ
(
U
0
H
δ
+
2
κ
)
-
1
4
.
To find the corresponding wave functions, referring to Table 2 and (29), we find the functions
(34)
ρ
(
s
)
=
s
-
2
β
2
/
δ
2
(
1
-
s
)
2
η
κ
,
ϕ
(
s
)
=
s
-
β
2
/
δ
2
(
1
-
s
)
η
κ
+
1
/
2
,
y
n
(
s
)
=
P
n
(
-
2
β
2
/
δ
2
,
2
η
κ
)
(
1
-
2
s
)
.
By using F
n
κ
(
s
)
=
ϕ
(
s
)
y
n
(
s
), we get the upper component of the Dirac spinor as
(35)
F
n
κ
(
r
)
=
N
n
κ
e
-
β
2
r
(
1
-
e
-
δ
r
)
η
κ
+
1
/
2
×
P
n
(
-
2
β
2
/
δ
2
,
2
η
κ
)
(
1
-
2
e
-
δ
r
)
,
where N
n
κ is normalization constant. The lower component of the Dirac spinor from (10) can be calculated as
(36)
G
n
κ
(
r
)
=
1
M
+
E
n
κ
-
C
s
(
d
d
r
+
κ
r
-
U
(
r
)
)
F
n
κ
(
r
)
,
where E
≠
-
M
+
C
s. Some numerical results of (33) are given in Table 3. In Table 3, we obtain numerical results in the spin symmetry in the presence and absence of tensor Hulthén form potential. We took a set of parameter values, M
=
5.0 fm−1, C
s
=
0, δ
=
0.025 fm−1, and α
=
0.01.Table 3
The bound state energy eigenvalues in unit off
m
-
1 of the spin symmetry Hulthén potential for several values of n and κ. Here the first and second states column according to (
l
,
j
).
|
n
,
l
,
κ
〉
State
α
=
0.01
U
0
H
=
0
α
=
0.01
U
0
H
=
0.2
|
n
,
l
,
κ
〉
State
α
=
0.01
U
0
H
=
0
α
=
0.01
U
0
H
=
0.1
E
n
,
κ (κ
<
0)
E
n
,
κ (κ
<
0)
E
n
,
κ (κ
>
0
)
E
n
,
κ (κ
>
0)
|
0
,
1
,
-
2
〉
0
p
3
/
2
-
4.9999843769
-
4.9999282575
|
0
,
1
,
1
〉
0
p
1
/
2
-
4.9999843769
-
4.9999566000
|
0
,
2
,
-
3
〉
0
d
5
/
2
-
4.9999843762
-
4.9999149381
|
0
,
2
,
2
〉
0
d
3
/
2
-
4.9999843762
-
4.999961417
8
|
0
,
3
,
-
4
〉
0
f
7
/
2
-
4.9999843759
-
4.9998943876
|
0
,
3
,
3
〉
0
f
5
/
2
-
4.9999843759
-
4.9999648452
|
0
,
4
,
-
5
〉
0
g
9
/
2
-
4.9999843757
-
4.9998593993
|
0
,
4
,
4
〉
0
g
7
/
2
-
4.9999843757
-
4.999967400
8
|
1
,
1
,
-
2
〉
1
p
3
/
2
-
4.9999149463
-
4.9997802943
|
1
,
1
,
1
〉
1
p
1
/
2
-
4.9999149463
-
4.999831330
7
|
1
,
2
,
-
3
〉
1
d
5
/
2
-
4.9999023580
-
4.9997675648
|
1
,
2
,
2
〉
1
d
3
/
2
-
4.9999023580
-
4.999834974
9
|
1
,
3
,
-
4
〉
1
f
7
/
2
-
4.9998943876
-
4.9997500354
|
1
,
3
,
3
〉
1
f
5
/
2
-
4.9998943876
-
4.999837781
9
|
1
,
4
,
-
5
〉
1
g
9
/
2
-
4.9998888999
-
4.9997244252
|
1
,
4
,
4
〉
1
g
7
/
2
-
4.9998888999
-
4.999840010
2We can observe that every pair of orbitals(
n
p
1
/
2
,
n
p
3
/
2
), (
n
d
3
/
2
,
n
d
5
/
2
), and (
n
f
5
/
2
,
n
f
7
/
2
) has the same energy in the absence of the tensor potential (
U
0
H
=
0
). Thus, they can be viewed as the spin doublets; that is, the state 1
p
1
/
2 with n
=
1 and κ
=
1 forms a spin doublet with the 1
p
3
/
2 state with n
=
1 and κ
=
-
2. On the other hand, in the presence of the tensor potential (
U
0
H
≠
0
), one can notice that degeneracy between every pair of spin doublets is removed. In Figure 2, we have investigated the effect of the tensor potential on the spin doublet splitting by considering the following pairs of orbital: (
1
p
3
/
2
,
1
p
1
/
2
), (
1
f
7
/
2
,
1
f
5
/
2
). From Figure 2, we observe that, in the case of (
U
0
H
=
0
) (no tensor interaction), members of spin doublets have the same energy. However, in the presence of the tensor potential (
U
0
H
≠
0
), these degeneracies are removed.Figure 2
Contribution of the tensor Hulthén potential parameter to the energy levels in the case of spin symmetry.


### 3.2. p-Spin Symmetry Solution

In this subsection we will obtain the energy eigenvalues and the corresponding wave functions for the p-spin symmetric limit of the Manning-Rosen potential, that is, solutions of (21). Substituting (22) and (23) into (21) gives
(37)
[
d
2
d
r
2
-
{
κ
(
κ
-
1
)
δ
2
+
(
2
κ
-
1
)
U
0
H
δ
+
U
0
H
2
+
γ
~
V
1
}
00
×
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
+
{
U
0
H
δ
-
γ
~
V
0
}
e
-
δ
r
1
-
e
-
δ
b
00
-
β
~
2
d
2
d
r
2
]
G
n
κ
(
r
)
=
0
,
and further making the change of variables s
=
e
-
δ
r, we obtain
(38)
d
2
G
n
κ
(
s
)
d
s
2
+
1
-
s
s
(
1
-
s
)
d
G
n
κ
(
s
)
d
s
+
1
δ
2
s
2
(
1
-
s
)
2
×
[
-
(
A
~
+
B
~
+
β
~
2
)
s
2
+
(
B
~
+
2
β
~
2
)
s
-
β
~
2
]
×
G
n
κ
(
s
)
=
0
,
where
(39)
A
~
=
κ
(
κ
-
1
)
δ
2
+
(
2
κ
-
1
)
U
0
H
δ
+
U
0
H
2
+
γ
~
V
1
B
~
=
U
0
H
δ
-
γ
~
V
0
.
On the other hand, to avoid repetition in the solution of (38), in the p-spin symmetric case,
(40)
(
β
~
2
δ
2
n
+
1
2
+
(
κ
-
1
2
)
2
+
U
0
H
δ
(
U
0
H
δ
+
2
κ
-
1
)
+
γ
~
V
1
δ
2
-
β
~
2
δ
2
)
2
=
γ
~
(
V
1
-
V
0
)
+
β
~
2
δ
2
+
(
κ
-
1
2
)
2
+
U
0
H
δ
(
U
0
H
δ
+
2
κ
)
-
1
4
.
Also, we get the lower component of the Dirac spinor as
(41)
G
n
κ
(
r
)
=
N
~
n
κ
e
-
β
~
2
r
(
1
-
e
-
δ
r
)
η
~
κ
+
1
/
2
P
n
(
-
2
β
~
2
/
δ
2
,
2
η
~
κ
)
×
(
1
-
2
e
-
δ
r
)
,
where N
~
n
κ is normalization constant and η
~
κ
=
(
κ
-
1
/
2
)
2
+
(
U
0
H
/
δ
)
(
U
0
H
/
δ
+
2
κ
-
1
)
+
(
γ
~
V
1
)
/
δ
2.The upper component of the Dirac spinor from (11) can be calculated as
(42)
F
n
κ
(
r
)
=
1
M
-
E
n
κ
+
C
ps
(
d
d
r
-
κ
r
+
U
(
r
)
)
G
n
κ
(
r
)
,
where E
≠
M
+
C
ps. In Table 4, we use the same parameters as in the previous subsection (with C
ps
=
0) and we obtain numerical results in the p-spin symmetry in the presence and absence of tensor Hulthén form potential. The existence of the degeneracy was observed in the following doublets: (
1
s
1
/
2
,
0
d
3
/
2
), (
1
p
3
/
2
,
0
f
5
/
2
), (
1
d
5
/
2
,
0
g
7
/
2
), (
1
f
7
/
2
,
0
h
9
/
2
), and so on, in the absence of the tensor potential (
U
0
H
=
0
), and they are considered as p-spin doublets. In Figure 3, we have investigated the effect of the tensor potential on the p-spin doublet splitting by considering the following pairs of orbital: (
1
d
5
/
2
,
0
g
7
/
2
), (
2
f
7
/
2
,
1
h
9
/
2
).Table 4
The bound state energy eigenvalues in unit off
m
-
1 of the p-spin symmetry Hulthén potential. Here the first states column according to (
l
,
j
) and second states column (
l
+
2
,
j
+
1
).
|
n
,
l
~
,
κ
〉
State
α
=
0.01
U
0
H
=
0
α
=
0.01
U
0
H
=
0.2
|
n
,
l
~
,
κ
〉
State
α
=
0.01
U
0
H
=
0
α
=
0.01
U
0
H
=
0.1
E
n
,
κ (κ
<
0)
E
n
,
κ (κ
<
0)
E
n
,
κ (κ
>
0
)
E
n
,
κ (κ
>
0)
|
1
,
1
,
-
1
〉
1
s
1
/
2
-
4.9868996028
-
4.9978434057
|
0
,
1
,
2
〉
0
d
3
/
2
-
4.9868996028
-
4.9971232725
|
1
,
2
,
-
2
〉
1
p
3
/
2
-
4.9921014723
-
4.9973604254
|
0
,
2
,
3
〉
0
f
5
/
2
-
4.9921014723
-
4.9976559580
|
1
,
3
,
-
3
〉
1
d
5
/
2
-
4.9945990389
-
4.9966407919
|
0
,
3
,
4
〉
0
g
7
/
2
-
4.9945990389
-
4.9980326519
|
1
,
4
,
-
4
〉
1
f
7
/
2
-
4.9960022563
-
4.9954882731
|
0
,
4
,
5
〉
0
h
9
/
2
-
4.9960022563
-
4.9983102937
|
2
,
1
,
-
1
〉
2
s
1
/
2
-
4.9913874762
-
4.9913874762
|
1
,
1
,
2
〉
1
d
3
/
2
-
4.9913874762
-
4.9970467176
|
2
,
2
,
-
2
〉
2
p
3
/
2
-
4.9938803034
-
4.9938803034
|
1
,
2
,
3
〉
1
f
5
/
2
-
4.9938803034
-
4.9974493983
|
2
,
3
,
-
3
〉
2
d
5
/
2
-
4.9953083139
-
4.9953083139
|
1
,
3
,
4
〉
1
g
7
/
2
-
4.9953083139
-
4.9977503418
|
2
,
4
,
-
4
〉
2
f
7
/
2
-
4.9962128946
-
4.9962128946
|
1
,
4
,
5
〉
1
h
9
/
2
-
4.9962128946
-
4.9979825424Figure 3
Contribution of the tensor Hulthén potential parameter to the energy levels in the case of p-spin symmetry.The energy spectrum of spin and p-spin symmetry given by Tables3 and 4 changes with the relevant quantum numbers as well as the screening parameter δ. The variation of the energy spectrum (
E
n
,
κ
) according to the screening parameter δ for some states is shown in Figures 4 and 5.Figure 4
The variation of the energy spectrumE
n
,
k for tensor Hulthén potential in the case of spin symmetry according to the screening parameter δ. Here we selected U
0
H
=
0.1, α
=
0.01.Figure 5
The variation of the energy spectrumE
n
,
k for tensor Hulthén potential in the case of p-spin symmetry according to the screening parameter δ. Here we selected U
0
H
=
0.1, α
=
0.01.


## 3.1. Spin Symmetry Solution

Substituting (22) and (23) into (18) gives
(30)
[
d
2
d
r
2
-
{
κ
(
κ
+
1
)
δ
2
+
(
2
κ
+
1
)
U
0
H
δ
+
U
0
H
2
+
γ
V
1
δ
2
}
00
×
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
+
{
-
U
0
H
δ
-
γ
V
0
}
00
×
e
-
δ
r
1
-
e
-
δ
b
-
β
2
]
F
n
κ
(
r
)
=
0
.
Defining a new variable of the form s
=
e
-
δ
r and using (30), we obtain
(31)
d
2
F
n
κ
(
s
)
d
s
2
+
1
-
s
s
(
1
-
s
)
d
F
n
κ
(
s
)
d
s
+
1
δ
2
s
2
(
1
-
s
)
2
×
[
-
(
A
+
B
+
β
2
)
s
2
+
(
B
+
2
β
2
)
s
-
β
2
]
×
F
n
κ
(
s
)
=
0
,
where
(32)
A
=
κ
(
κ
+
1
)
δ
2
+
(
2
κ
+
1
)
U
0
H
δ
+
U
0
H
2
+
γ
V
1
,
B
=
-
U
0
H
δ
-
γ
V
0
.
Comparing (31) and (25), we can easily obtain the coefficients c
i (i
=
1,2
,
3) and analytical expressions ξ
j (j
=
0,1
,
2). By using Table 1, we are obtained the coefficients c
i (i
=
4,5
,
.
.
.
,
13). The values of all coefficients are shown in Table 2.Table 2
The specific values for the parametric constants necessary for calculating the energy eigenvalues and eigenfunctions of the spin symmetry limit.
Constant
Analytic value
Constant
Analytic value
c
1
c
1
1
c
2
1
c
3
1
c
4
0
c
5
-
1
2
c
6
1
4
+
A
+
B
+
β
2
δ
2
c
7
-
B
+
2
β
2
δ
2
c
8
β
2
δ
2
c
9
1
4
+
A
δ
2
c
10
-
2
β
2
δ
2
c
11
2
η
κ
c
12
-
β
2
δ
2
c
13
η
κ
+
1
2
ξ
0
β
2
δ
2
ξ
1
B
+
2
β
2
δ
2
ξ
2
A
+
B
+
β
2
δ
2In Table2, η
κ
=
(
(
κ
+
1
/
2
)
2
+
(
U
0
H
/
δ
)
(
(
U
0
H
/
δ
)
+
2
κ
+
1
)
+
γ
V
1
/
δ
2
)
1
/
2, and by using (28), we can obtain the energy eigenvalues of the Manning-Rosen potential with spin symmetric limit as
(33)
(
β
2
δ
2
n
+
1
2
+
(
κ
+
1
2
)
2
+
U
0
H
δ
(
U
0
H
δ
+
2
κ
+
1
)
+
γ
V
1
δ
2
-
β
2
δ
2
)
2
=
γ
(
V
1
-
V
0
)
+
β
2
δ
2
+
(
κ
+
1
2
)
2
+
U
0
H
δ
(
U
0
H
δ
+
2
κ
)
-
1
4
.
To find the corresponding wave functions, referring to Table 2 and (29), we find the functions
(34)
ρ
(
s
)
=
s
-
2
β
2
/
δ
2
(
1
-
s
)
2
η
κ
,
ϕ
(
s
)
=
s
-
β
2
/
δ
2
(
1
-
s
)
η
κ
+
1
/
2
,
y
n
(
s
)
=
P
n
(
-
2
β
2
/
δ
2
,
2
η
κ
)
(
1
-
2
s
)
.
By using F
n
κ
(
s
)
=
ϕ
(
s
)
y
n
(
s
), we get the upper component of the Dirac spinor as
(35)
F
n
κ
(
r
)
=
N
n
κ
e
-
β
2
r
(
1
-
e
-
δ
r
)
η
κ
+
1
/
2
×
P
n
(
-
2
β
2
/
δ
2
,
2
η
κ
)
(
1
-
2
e
-
δ
r
)
,
where N
n
κ is normalization constant. The lower component of the Dirac spinor from (10) can be calculated as
(36)
G
n
κ
(
r
)
=
1
M
+
E
n
κ
-
C
s
(
d
d
r
+
κ
r
-
U
(
r
)
)
F
n
κ
(
r
)
,
where E
≠
-
M
+
C
s. Some numerical results of (33) are given in Table 3. In Table 3, we obtain numerical results in the spin symmetry in the presence and absence of tensor Hulthén form potential. We took a set of parameter values, M
=
5.0 fm−1, C
s
=
0, δ
=
0.025 fm−1, and α
=
0.01.Table 3
The bound state energy eigenvalues in unit off
m
-
1 of the spin symmetry Hulthén potential for several values of n and κ. Here the first and second states column according to (
l
,
j
).
|
n
,
l
,
κ
〉
State
α
=
0.01
U
0
H
=
0
α
=
0.01
U
0
H
=
0.2
|
n
,
l
,
κ
〉
State
α
=
0.01
U
0
H
=
0
α
=
0.01
U
0
H
=
0.1
E
n
,
κ (κ
<
0)
E
n
,
κ (κ
<
0)
E
n
,
κ (κ
>
0
)
E
n
,
κ (κ
>
0)
|
0
,
1
,
-
2
〉
0
p
3
/
2
-
4.9999843769
-
4.9999282575
|
0
,
1
,
1
〉
0
p
1
/
2
-
4.9999843769
-
4.9999566000
|
0
,
2
,
-
3
〉
0
d
5
/
2
-
4.9999843762
-
4.9999149381
|
0
,
2
,
2
〉
0
d
3
/
2
-
4.9999843762
-
4.999961417
8
|
0
,
3
,
-
4
〉
0
f
7
/
2
-
4.9999843759
-
4.9998943876
|
0
,
3
,
3
〉
0
f
5
/
2
-
4.9999843759
-
4.9999648452
|
0
,
4
,
-
5
〉
0
g
9
/
2
-
4.9999843757
-
4.9998593993
|
0
,
4
,
4
〉
0
g
7
/
2
-
4.9999843757
-
4.999967400
8
|
1
,
1
,
-
2
〉
1
p
3
/
2
-
4.9999149463
-
4.9997802943
|
1
,
1
,
1
〉
1
p
1
/
2
-
4.9999149463
-
4.999831330
7
|
1
,
2
,
-
3
〉
1
d
5
/
2
-
4.9999023580
-
4.9997675648
|
1
,
2
,
2
〉
1
d
3
/
2
-
4.9999023580
-
4.999834974
9
|
1
,
3
,
-
4
〉
1
f
7
/
2
-
4.9998943876
-
4.9997500354
|
1
,
3
,
3
〉
1
f
5
/
2
-
4.9998943876
-
4.999837781
9
|
1
,
4
,
-
5
〉
1
g
9
/
2
-
4.9998888999
-
4.9997244252
|
1
,
4
,
4
〉
1
g
7
/
2
-
4.9998888999
-
4.999840010
2We can observe that every pair of orbitals(
n
p
1
/
2
,
n
p
3
/
2
), (
n
d
3
/
2
,
n
d
5
/
2
), and (
n
f
5
/
2
,
n
f
7
/
2
) has the same energy in the absence of the tensor potential (
U
0
H
=
0
). Thus, they can be viewed as the spin doublets; that is, the state 1
p
1
/
2 with n
=
1 and κ
=
1 forms a spin doublet with the 1
p
3
/
2 state with n
=
1 and κ
=
-
2. On the other hand, in the presence of the tensor potential (
U
0
H
≠
0
), one can notice that degeneracy between every pair of spin doublets is removed. In Figure 2, we have investigated the effect of the tensor potential on the spin doublet splitting by considering the following pairs of orbital: (
1
p
3
/
2
,
1
p
1
/
2
), (
1
f
7
/
2
,
1
f
5
/
2
). From Figure 2, we observe that, in the case of (
U
0
H
=
0
) (no tensor interaction), members of spin doublets have the same energy. However, in the presence of the tensor potential (
U
0
H
≠
0
), these degeneracies are removed.Figure 2
Contribution of the tensor Hulthén potential parameter to the energy levels in the case of spin symmetry.


## 3.2. p-Spin Symmetry Solution

In this subsection we will obtain the energy eigenvalues and the corresponding wave functions for the p-spin symmetric limit of the Manning-Rosen potential, that is, solutions of (21). Substituting (22) and (23) into (21) gives
(37)
[
d
2
d
r
2
-
{
κ
(
κ
-
1
)
δ
2
+
(
2
κ
-
1
)
U
0
H
δ
+
U
0
H
2
+
γ
~
V
1
}
00
×
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
+
{
U
0
H
δ
-
γ
~
V
0
}
e
-
δ
r
1
-
e
-
δ
b
00
-
β
~
2
d
2
d
r
2
]
G
n
κ
(
r
)
=
0
,
and further making the change of variables s
=
e
-
δ
r, we obtain
(38)
d
2
G
n
κ
(
s
)
d
s
2
+
1
-
s
s
(
1
-
s
)
d
G
n
κ
(
s
)
d
s
+
1
δ
2
s
2
(
1
-
s
)
2
×
[
-
(
A
~
+
B
~
+
β
~
2
)
s
2
+
(
B
~
+
2
β
~
2
)
s
-
β
~
2
]
×
G
n
κ
(
s
)
=
0
,
where
(39)
A
~
=
κ
(
κ
-
1
)
δ
2
+
(
2
κ
-
1
)
U
0
H
δ
+
U
0
H
2
+
γ
~
V
1
B
~
=
U
0
H
δ
-
γ
~
V
0
.
On the other hand, to avoid repetition in the solution of (38), in the p-spin symmetric case,
(40)
(
β
~
2
δ
2
n
+
1
2
+
(
κ
-
1
2
)
2
+
U
0
H
δ
(
U
0
H
δ
+
2
κ
-
1
)
+
γ
~
V
1
δ
2
-
β
~
2
δ
2
)
2
=
γ
~
(
V
1
-
V
0
)
+
β
~
2
δ
2
+
(
κ
-
1
2
)
2
+
U
0
H
δ
(
U
0
H
δ
+
2
κ
)
-
1
4
.
Also, we get the lower component of the Dirac spinor as
(41)
G
n
κ
(
r
)
=
N
~
n
κ
e
-
β
~
2
r
(
1
-
e
-
δ
r
)
η
~
κ
+
1
/
2
P
n
(
-
2
β
~
2
/
δ
2
,
2
η
~
κ
)
×
(
1
-
2
e
-
δ
r
)
,
where N
~
n
κ is normalization constant and η
~
κ
=
(
κ
-
1
/
2
)
2
+
(
U
0
H
/
δ
)
(
U
0
H
/
δ
+
2
κ
-
1
)
+
(
γ
~
V
1
)
/
δ
2.The upper component of the Dirac spinor from (11) can be calculated as
(42)
F
n
κ
(
r
)
=
1
M
-
E
n
κ
+
C
ps
(
d
d
r
-
κ
r
+
U
(
r
)
)
G
n
κ
(
r
)
,
where E
≠
M
+
C
ps. In Table 4, we use the same parameters as in the previous subsection (with C
ps
=
0) and we obtain numerical results in the p-spin symmetry in the presence and absence of tensor Hulthén form potential. The existence of the degeneracy was observed in the following doublets: (
1
s
1
/
2
,
0
d
3
/
2
), (
1
p
3
/
2
,
0
f
5
/
2
), (
1
d
5
/
2
,
0
g
7
/
2
), (
1
f
7
/
2
,
0
h
9
/
2
), and so on, in the absence of the tensor potential (
U
0
H
=
0
), and they are considered as p-spin doublets. In Figure 3, we have investigated the effect of the tensor potential on the p-spin doublet splitting by considering the following pairs of orbital: (
1
d
5
/
2
,
0
g
7
/
2
), (
2
f
7
/
2
,
1
h
9
/
2
).Table 4
The bound state energy eigenvalues in unit off
m
-
1 of the p-spin symmetry Hulthén potential. Here the first states column according to (
l
,
j
) and second states column (
l
+
2
,
j
+
1
).
|
n
,
l
~
,
κ
〉
State
α
=
0.01
U
0
H
=
0
α
=
0.01
U
0
H
=
0.2
|
n
,
l
~
,
κ
〉
State
α
=
0.01
U
0
H
=
0
α
=
0.01
U
0
H
=
0.1
E
n
,
κ (κ
<
0)
E
n
,
κ (κ
<
0)
E
n
,
κ (κ
>
0
)
E
n
,
κ (κ
>
0)
|
1
,
1
,
-
1
〉
1
s
1
/
2
-
4.9868996028
-
4.9978434057
|
0
,
1
,
2
〉
0
d
3
/
2
-
4.9868996028
-
4.9971232725
|
1
,
2
,
-
2
〉
1
p
3
/
2
-
4.9921014723
-
4.9973604254
|
0
,
2
,
3
〉
0
f
5
/
2
-
4.9921014723
-
4.9976559580
|
1
,
3
,
-
3
〉
1
d
5
/
2
-
4.9945990389
-
4.9966407919
|
0
,
3
,
4
〉
0
g
7
/
2
-
4.9945990389
-
4.9980326519
|
1
,
4
,
-
4
〉
1
f
7
/
2
-
4.9960022563
-
4.9954882731
|
0
,
4
,
5
〉
0
h
9
/
2
-
4.9960022563
-
4.9983102937
|
2
,
1
,
-
1
〉
2
s
1
/
2
-
4.9913874762
-
4.9913874762
|
1
,
1
,
2
〉
1
d
3
/
2
-
4.9913874762
-
4.9970467176
|
2
,
2
,
-
2
〉
2
p
3
/
2
-
4.9938803034
-
4.9938803034
|
1
,
2
,
3
〉
1
f
5
/
2
-
4.9938803034
-
4.9974493983
|
2
,
3
,
-
3
〉
2
d
5
/
2
-
4.9953083139
-
4.9953083139
|
1
,
3
,
4
〉
1
g
7
/
2
-
4.9953083139
-
4.9977503418
|
2
,
4
,
-
4
〉
2
f
7
/
2
-
4.9962128946
-
4.9962128946
|
1
,
4
,
5
〉
1
h
9
/
2
-
4.9962128946
-
4.9979825424Figure 3
Contribution of the tensor Hulthén potential parameter to the energy levels in the case of p-spin symmetry.The energy spectrum of spin and p-spin symmetry given by Tables3 and 4 changes with the relevant quantum numbers as well as the screening parameter δ. The variation of the energy spectrum (
E
n
,
κ
) according to the screening parameter δ for some states is shown in Figures 4 and 5.Figure 4
The variation of the energy spectrumE
n
,
k for tensor Hulthén potential in the case of spin symmetry according to the screening parameter δ. Here we selected U
0
H
=
0.1, α
=
0.01.Figure 5
The variation of the energy spectrumE
n
,
k for tensor Hulthén potential in the case of p-spin symmetry according to the screening parameter δ. Here we selected U
0
H
=
0.1, α
=
0.01.


## 4. Bound States of the Manning-Rosen Potential with Coulomb Tensor Interaction

In this stage we will solve Dirac (13) and (14) that are obtained in Section 2, for Manning-Rosen potential including a Coulomb-like tensor interaction under consideration of spin and pseudospin symmetry case.Under spin symmetry, from (13), for the upper radial spinor component, we obtain
(43)
[
d
2
d
r
2
-
λ
κ
(
λ
κ
+
1
)
r
2
-
γ
(
V
0
e
-
δ
r
1
-
e
-
δ
b
+
V
1
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
)
-
β
2
d
2
d
r
2
]
F
n
κ
(
r
)
=
0
,
where λ
κ
=
κ
+
U
0
C, γ
=
M
+
E
n
κ
-
C
s, and β
2
=
(
M
-
E
n
κ
)
(
M
+
E
n
κ
-
C
s
). Also, κ
=
l and κ
=
-
l
-
1 for κ
<
0 and κ
>
0, respectively.Also under p-spin symmetry, from (14), for the lower radial spinor component, we have
(44)
[
d
2
d
r
2
-
λ
κ
(
λ
κ
-
1
)
r
2
-
γ
~
(
V
0
e
-
δ
r
1
-
e
-
δ
b
+
V
1
e
-
2
δ
r
(
1
-
e
-
δ
r
)
2
)
-
β
~
2
d
2
d
r
2
]
G
n
κ
(
r
)
=
0
,
where λ
κ
=
κ
+
U
0
C, γ
~
=
E
n
κ
-
M
-
C
ps, and β
~
2
=
(
M
+
E
n
κ
)
(
M
-
E
n
κ
+
C
ps
). Also, κ
=
-
l
~ and κ
=
l
~
+
1 for κ
<
0 and κ
>
0, respectively.


### 4.1. Spin Symmetry Solution

Substituting (22) and (23) into (43) and defining a new variable of the form s
=
e
-
δ
r, we get
(45)
d
2
F
n
κ
(
s
)
d
s
2
+
1
-
s
s
(
1
-
s
)
d
F
n
κ
(
s
)
d
s
+
1
δ
2
s
2
(
1
-
s
)
2
×
[
-
(
A
+
B
+
β
2
)
s
2
+
(
B
+
2
β
2
)
s
-
β
2
]
×
F
n
κ
(
s
)
=
0
.
Here
(46)
A
=
λ
κ
(
λ
κ
+
1
)
δ
2
+
γ
V
1
,
B
=
-
γ
V
0
.
Comparing (45) with relations that are obtained in Section 3.1, we can easily obtain that
(47)
η
κ
=
1
4
+
A
δ
2
=
(
κ
+
U
0
C
+
1
2
)
2
+
γ
V
1
δ
2
.
We can obtain the energy eigenvalues of the Manning-Rosen potential with spin symmetric limit as
(48)
(
n
+
1
2
+
(
κ
+
U
0
C
+
1
2
)
2
+
γ
V
1
δ
2
-
β
2
δ
2
)
2
=
γ
(
V
1
-
V
0
)
+
β
2
δ
2
+
(
κ
+
U
0
C
+
1
2
)
2
-
1
4
.
To find the corresponding wave functions, we can use (34), (35), and (36). Some numerical results of (48) are given in Table 5. In Table 5, we obtain numerical results in the spin symmetry in the presence and absence of tensor Coulomb form potential. We took a set of parameter values: M
=
5
.
0 fm−1, C
s
=
0, δ
=
0.025 fm−1, and α
=
0.01. In Figure 6, we have investigated the effect of the tensor potential on the spin doublet splitting by considering the following pairs of orbital: (
1
p
3
/
2
,
1
p
1
/
2
), (
1
f
7
/
2
,
1
f
5
/
2
).Table 5
The bound state energy eigenvalues in unit off
m
-
1 of the spin symmetry Coulomb potential for several values of n and κ. Here the first and second states column according to (
l
,
j
).
|
n
,
l
,
κ
〉
State
α
=
0.01
U
0
C
=
0
α
=
0.01
U
0
C
=
1
|
n
,
l
,
κ
〉
State
α
=
0.01
U
0
C
=
0
α
=
0.01
U
0
C
=
1
E
n
,
κ (κ
<
0
)
E
n
,
κ (κ
<
0)
E
n
,
κ (κ
>
0
)
E
n
,
κ (κ
>
0)
|
0
,
1
,
-
2
〉
0
p
3
/
2
-
4.9999843769
-
4.9999843788
|
0
,
1
,
1
〉
0
p
1
/
2
-
4.9999843769
-
4.9999843762
|
0
,
2
,
-
3
〉
0
d
5
/
2
-
4.9999843762
-
4.9999843769
|
0
,
2
,
2
〉
0
d
3
/
2
-
4.9999843762
-
4.9999843759
|
0
,
3
,
-
4
〉
0
f
7
/
2
-
4.9999843759
-
4.9999843762
|
0
,
3
,
3
〉
0
f
5
/
2
-
4.9999843759
-
4.9999843757
|
0
,
4
,
-
5
〉
0
g
9
/
2
-
4.9999843757
-
4.9999843759
|
0
,
4
,
4
〉
0
g
7
/
2
-
4.9999843757
-
4.9999843756
|
1
,
1
,
-
2
〉
1
p
3
/
2
-
4.9999149463
-
4.9999375152
|
1
,
1
,
1
〉
1
p
1
/
2
-
4.9999149463
-
4.999902358
0
|
1
,
2
,
-
3
〉
1
d
5
/
2
-
4.9999023580
-
4.9999149463
|
1
,
2
,
2
〉
1
d
3
/
2
-
4.9999023580
-
4.9998943876
|
1
,
3
,
-
4
〉
1
f
7
/
2
-
4.9998943876
-
4.999902358
0
|
1
,
3
,
3
〉
1
f
5
/
2
-
4.9998943876
-
4.9998888999
|
1
,
4
,
-
5
〉
1
g
9
/
2
-
4.9998888999
-
4.9998943876
|
1
,
4
,
4
〉
1
g
7
/
2
-
4.9998888999
-
4.9998848950Figure 6
Contribution of the tensor Coulomb potential parameter to the energy levels in the case of spin symmetry.


### 4.2. p-Spin Symmetry Solution

Substituting (22) and (23) into (44) and defining a new variable of the form s
=
e
-
δ
r, we get
(49)
d
2
G
n
κ
(
s
)
d
s
2
+
1
-
s
s
(
1
-
s
)
d
G
n
κ
(
s
)
d
s
+
1
δ
2
s
2
(
1
-
s
)
2
×
[
-
(
A
~
+
B
~
+
β
~
2
)
s
2
+
(
B
~
+
2
β
~
2
)
s
-
β
~
2
]
×
G
n
κ
(
s
)
=
0
.
Here
(50)
A
~
=
λ
κ
(
λ
κ
-
1
)
δ
2
+
γ
~
V
1
,
B
~
=
-
γ
~
V
0
.
Comparing (49) with relations that are obtained in Section 3.2, we can easily obtain that
(51)
η
~
κ
=
1
4
+
A
~
δ
2
=
(
κ
+
U
0
C
-
1
2
)
2
+
γ
~
V
1
δ
2
.
We can obtain the energy eigenvalues of the Manning-Rosen potential with p-spin symmetric limit as
(52)
(
n
+
1
2
+
(
κ
+
U
0
C
-
1
2
)
2
+
γ
~
V
1
δ
2
-
β
~
2
δ
2
)
2
=
γ
~
(
V
1
-
V
0
)
+
β
~
2
δ
2
+
(
κ
+
U
0
C
-
1
2
)
2
-
1
4
.
Also, we get the lower and upper component of the Dirac spinor as (41) and (42). In Table 6, we use the same parameters as in the previous subsection (with C
ps
=
0) and we obtain numerical results in the p-spin symmetry in the presence and absence of tensor Coulomb form potential. In Figure 7, we have investigated the effect of the tensor potential on the p-spin doublet splitting by considering some pairs of orbitals. Also the variation of the energy spectrum (
E
n
,
κ
) of spin and p-spin symmetry given by Tables 5 and 6, according to the screening parameter δ for some states, is shown in Figures 8 and 9.Table 6
The bound state energy eigenvalues in unit off
m
-
1 of the p-spin symmetry Coulomb potential. Here the first states column according to (
l
,
j
) and second states column (
l
+
2
,
j
+
1
).
|
n
,
l
~
,
κ
〉
State
α
=
0.01
U
0
C
=
0
α
=
0.01
U
0
C
=
1
|
n
,
l
~
,
κ
〉
State
α
=
0.01
U
0
C
=
0
α
=
0.01
U
0
C
=
1
E
n
,
κ (κ
<
0
)
E
n
,
κ (κ
<
0)
E
n
,
κ (κ
>
0
)
E
n
,
κ (κ
>
0)
|
1
,
1
,
-
1
〉
1
s
1
/
2
-
4.9868996028
-
4.9727532453
|
0
,
1
,
2
〉
0
d
3
/
2
-
4.9868996028
-
4.9921014723
|
1
,
2
,
-
2
〉
1
p
3
/
2
-
4.9921014723
-
4.9868996028
|
0
,
2
,
3
〉
0
f
5
/
2
-
4.9921014723
-
4.9945990389
|
1
,
3
,
-
3
〉
1
d
5
/
2
-
4.9945990389
-
4.9921014723
|
0
,
3
,
4
〉
0
g
7
/
2
-
4.9945990389
-
4.9960022563
|
1
,
4
,
-
4
〉
1
f
7
/
2
-
4.9960022563
-
4.9945990389
|
0
,
4
,
5
〉
0
h
9
/
2
-
4.9960022563
-
4.9968759153
|
2
,
1
,
-
1
〉
2
s
1
/
2
-
4.9913874762
-
4.9863374401
|
1
,
1
,
2
〉
1
d
3
/
2
-
4.9913874762
-
4.9938803034
|
2
,
2
,
-
2
〉
2
p
3
/
2
-
4.9938803034
-
4.9913874762
|
1
,
2
,
3
〉
1
f
5
/
2
-
4.9938803034
-
4.9953083139
|
2
,
3
,
-
3
〉
2
d
5
/
2
-
4.9953083139
-
4.9938803034
|
1
,
3
,
4
〉
1
g
7
/
2
-
4.9953083139
-
4.9962128946
|
2
,
4
,
-
4
〉
2
f
7
/
2
-
4.9962128946
-
4.9953083139
|
1
,
4
,
5
〉
1
h
9
/
2
-
4.9962128946
-
4.9968280175Figure 7
Contribution of the tensor Coulomb potential parameter to the energy levels in the case of p-spin symmetry.Figure 8
The variation of the energy spectrumE
n
,
k for tensor Coulomb potential in the case of spin symmetry according to the screening parameter δ. Here we selected U
0
C
=
0.1, α
=
0.01.Figure 9
The variation of the energy spectrumE
n
,
k for tensor Coulomb potential in the case of p-spin symmetry according to the screening parameter δ. Here we selected U
0
C
=
0.1, α
=
0.01.


## 4.1. Spin Symmetry Solution

Substituting (22) and (23) into (43) and defining a new variable of the form s
=
e
-
δ
r, we get
(45)
d
2
F
n
κ
(
s
)
d
s
2
+
1
-
s
s
(
1
-
s
)
d
F
n
κ
(
s
)
d
s
+
1
δ
2
s
2
(
1
-
s
)
2
×
[
-
(
A
+
B
+
β
2
)
s
2
+
(
B
+
2
β
2
)
s
-
β
2
]
×
F
n
κ
(
s
)
=
0
.
Here
(46)
A
=
λ
κ
(
λ
κ
+
1
)
δ
2
+
γ
V
1
,
B
=
-
γ
V
0
.
Comparing (45) with relations that are obtained in Section 3.1, we can easily obtain that
(47)
η
κ
=
1
4
+
A
δ
2
=
(
κ
+
U
0
C
+
1
2
)
2
+
γ
V
1
δ
2
.
We can obtain the energy eigenvalues of the Manning-Rosen potential with spin symmetric limit as
(48)
(
n
+
1
2
+
(
κ
+
U
0
C
+
1
2
)
2
+
γ
V
1
δ
2
-
β
2
δ
2
)
2
=
γ
(
V
1
-
V
0
)
+
β
2
δ
2
+
(
κ
+
U
0
C
+
1
2
)
2
-
1
4
.
To find the corresponding wave functions, we can use (34), (35), and (36). Some numerical results of (48) are given in Table 5. In Table 5, we obtain numerical results in the spin symmetry in the presence and absence of tensor Coulomb form potential. We took a set of parameter values: M
=
5
.
0 fm−1, C
s
=
0, δ
=
0.025 fm−1, and α
=
0.01. In Figure 6, we have investigated the effect of the tensor potential on the spin doublet splitting by considering the following pairs of orbital: (
1
p
3
/
2
,
1
p
1
/
2
), (
1
f
7
/
2
,
1
f
5
/
2
).Table 5
The bound state energy eigenvalues in unit off
m
-
1 of the spin symmetry Coulomb potential for several values of n and κ. Here the first and second states column according to (
l
,
j
).
|
n
,
l
,
κ
〉
State
α
=
0.01
U
0
C
=
0
α
=
0.01
U
0
C
=
1
|
n
,
l
,
κ
〉
State
α
=
0.01
U
0
C
=
0
α
=
0.01
U
0
C
=
1
E
n
,
κ (κ
<
0
)
E
n
,
κ (κ
<
0)
E
n
,
κ (κ
>
0
)
E
n
,
κ (κ
>
0)
|
0
,
1
,
-
2
〉
0
p
3
/
2
-
4.9999843769
-
4.9999843788
|
0
,
1
,
1
〉
0
p
1
/
2
-
4.9999843769
-
4.9999843762
|
0
,
2
,
-
3
〉
0
d
5
/
2
-
4.9999843762
-
4.9999843769
|
0
,
2
,
2
〉
0
d
3
/
2
-
4.9999843762
-
4.9999843759
|
0
,
3
,
-
4
〉
0
f
7
/
2
-
4.9999843759
-
4.9999843762
|
0
,
3
,
3
〉
0
f
5
/
2
-
4.9999843759
-
4.9999843757
|
0
,
4
,
-
5
〉
0
g
9
/
2
-
4.9999843757
-
4.9999843759
|
0
,
4
,
4
〉
0
g
7
/
2
-
4.9999843757
-
4.9999843756
|
1
,
1
,
-
2
〉
1
p
3
/
2
-
4.9999149463
-
4.9999375152
|
1
,
1
,
1
〉
1
p
1
/
2
-
4.9999149463
-
4.999902358
0
|
1
,
2
,
-
3
〉
1
d
5
/
2
-
4.9999023580
-
4.9999149463
|
1
,
2
,
2
〉
1
d
3
/
2
-
4.9999023580
-
4.9998943876
|
1
,
3
,
-
4
〉
1
f
7
/
2
-
4.9998943876
-
4.999902358
0
|
1
,
3
,
3
〉
1
f
5
/
2
-
4.9998943876
-
4.9998888999
|
1
,
4
,
-
5
〉
1
g
9
/
2
-
4.9998888999
-
4.9998943876
|
1
,
4
,
4
〉
1
g
7
/
2
-
4.9998888999
-
4.9998848950Figure 6
Contribution of the tensor Coulomb potential parameter to the energy levels in the case of spin symmetry.


## 4.2. p-Spin Symmetry Solution

Substituting (22) and (23) into (44) and defining a new variable of the form s
=
e
-
δ
r, we get
(49)
d
2
G
n
κ
(
s
)
d
s
2
+
1
-
s
s
(
1
-
s
)
d
G
n
κ
(
s
)
d
s
+
1
δ
2
s
2
(
1
-
s
)
2
×
[
-
(
A
~
+
B
~
+
β
~
2
)
s
2
+
(
B
~
+
2
β
~
2
)
s
-
β
~
2
]
×
G
n
κ
(
s
)
=
0
.
Here
(50)
A
~
=
λ
κ
(
λ
κ
-
1
)
δ
2
+
γ
~
V
1
,
B
~
=
-
γ
~
V
0
.
Comparing (49) with relations that are obtained in Section 3.2, we can easily obtain that
(51)
η
~
κ
=
1
4
+
A
~
δ
2
=
(
κ
+
U
0
C
-
1
2
)
2
+
γ
~
V
1
δ
2
.
We can obtain the energy eigenvalues of the Manning-Rosen potential with p-spin symmetric limit as
(52)
(
n
+
1
2
+
(
κ
+
U
0
C
-
1
2
)
2
+
γ
~
V
1
δ
2
-
β
~
2
δ
2
)
2
=
γ
~
(
V
1
-
V
0
)
+
β
~
2
δ
2
+
(
κ
+
U
0
C
-
1
2
)
2
-
1
4
.
Also, we get the lower and upper component of the Dirac spinor as (41) and (42). In Table 6, we use the same parameters as in the previous subsection (with C
ps
=
0) and we obtain numerical results in the p-spin symmetry in the presence and absence of tensor Coulomb form potential. In Figure 7, we have investigated the effect of the tensor potential on the p-spin doublet splitting by considering some pairs of orbitals. Also the variation of the energy spectrum (
E
n
,
κ
) of spin and p-spin symmetry given by Tables 5 and 6, according to the screening parameter δ for some states, is shown in Figures 8 and 9.Table 6
The bound state energy eigenvalues in unit off
m
-
1 of the p-spin symmetry Coulomb potential. Here the first states column according to (
l
,
j
) and second states column (
l
+
2
,
j
+
1
).
|
n
,
l
~
,
κ
〉
State
α
=
0.01
U
0
C
=
0
α
=
0.01
U
0
C
=
1
|
n
,
l
~
,
κ
〉
State
α
=
0.01
U
0
C
=
0
α
=
0.01
U
0
C
=
1
E
n
,
κ (κ
<
0
)
E
n
,
κ (κ
<
0)
E
n
,
κ (κ
>
0
)
E
n
,
κ (κ
>
0)
|
1
,
1
,
-
1
〉
1
s
1
/
2
-
4.9868996028
-
4.9727532453
|
0
,
1
,
2
〉
0
d
3
/
2
-
4.9868996028
-
4.9921014723
|
1
,
2
,
-
2
〉
1
p
3
/
2
-
4.9921014723
-
4.9868996028
|
0
,
2
,
3
〉
0
f
5
/
2
-
4.9921014723
-
4.9945990389
|
1
,
3
,
-
3
〉
1
d
5
/
2
-
4.9945990389
-
4.9921014723
|
0
,
3
,
4
〉
0
g
7
/
2
-
4.9945990389
-
4.9960022563
|
1
,
4
,
-
4
〉
1
f
7
/
2
-
4.9960022563
-
4.9945990389
|
0
,
4
,
5
〉
0
h
9
/
2
-
4.9960022563
-
4.9968759153
|
2
,
1
,
-
1
〉
2
s
1
/
2
-
4.9913874762
-
4.9863374401
|
1
,
1
,
2
〉
1
d
3
/
2
-
4.9913874762
-
4.9938803034
|
2
,
2
,
-
2
〉
2
p
3
/
2
-
4.9938803034
-
4.9913874762
|
1
,
2
,
3
〉
1
f
5
/
2
-
4.9938803034
-
4.9953083139
|
2
,
3
,
-
3
〉
2
d
5
/
2
-
4.9953083139
-
4.9938803034
|
1
,
3
,
4
〉
1
g
7
/
2
-
4.9953083139
-
4.9962128946
|
2
,
4
,
-
4
〉
2
f
7
/
2
-
4.9962128946
-
4.9953083139
|
1
,
4
,
5
〉
1
h
9
/
2
-
4.9962128946
-
4.9968280175Figure 7
Contribution of the tensor Coulomb potential parameter to the energy levels in the case of p-spin symmetry.Figure 8
The variation of the energy spectrumE
n
,
k for tensor Coulomb potential in the case of spin symmetry according to the screening parameter δ. Here we selected U
0
C
=
0.1, α
=
0.01.Figure 9
The variation of the energy spectrumE
n
,
k for tensor Coulomb potential in the case of p-spin symmetry according to the screening parameter δ. Here we selected U
0
C
=
0.1, α
=
0.01.


## 5. Conclusion

In this paper, we obtained the approximate analytical solutions of the Dirac equation for the Manning-Rosen potential with the consideration of Hulthén and Coulomb-like potentials as tensor interaction within the framework of spin and pseudospin symmetry limits using the generalized parametric NU technique. We have obtained the energy levels in a closed form and the corresponding wave functions in terms of the Jacobi polynomials. Also we obtained the energy eigenvalues in the absence(
U
0
H
=
0
) and the presence (
U
0
H
=
0.1
,
0.2
) of the Hulthén tensor interaction and the energy eigenvalues in the absence (
U
0
C
=
0
) and the presence (
U
0
C
=
1
) of the Coulomb-like tensor interaction for various values of the quantum numbers n and κ. The results are reported in Tables 3–6 under the condition of the spin and pseudospin symmetries, respectively, where we can see the way the tensor interaction affects the degeneracy of the system. Figures 2, 3, 6, and 7 present the effect of the Hulthén and Coulomb-like tensor interaction on the bound states, respectively. Obviously, the degeneracy between the members of doublet states in spin and pseudospin symmetries is removed by these tensor interaction potentials. We represent the effects of the screening parameter δ on the bound states for Hulthén and Coulomb-like tensor interaction in Figures 4, 5, 8, and 9, respectively. It is seen that if the δ-parameter increases, for both tensor interactions, the bound states become more bounded in both for symmetry limits. Therefore, the results of our work find many applications in both nuclear and Hadron physics and they obviously indicate the degeneracy-removing role of tensor interaction [104], which plays a significant role in the spectroscopes.



---

*Source: 870523-2014-03-20.xml*