Project supported by the Higher Education Project (Grant No. 698/UN27.11/PN/2015).

The Dirac equation for Eckart potential and trigonometric Manning Rosen potential with exact spin symmetry is obtained using an asymptotic iteration method. The combination of the two potentials is substituted into the Dirac equation, then the variables are separated into radial and angular parts. The Dirac equation is solved by using an asymptotic iteration method that can reduce the second order differential equation into a differential equation with substitution variables of hypergeometry type. The relativistic energy is calculated using Matlab 2011. This study is limited to the case of spin symmetry. With the asymptotic iteration method, the energy spectra of the relativistic equations and equations of orbital quantum number *l *can be obtained, where both are interrelated between quantum numbers. The energy spectrum is also numerically solved using the Matlab software, where the increase in the radial quantum number *n *_{r }causes the energy to decrease. The radial part and the angular part of the wave function are defined as hypergeometry functions and visualized with Matlab 2011. The results show that the disturbance of a combination of the Eckart potential and trigonometric Manning Rosen potential can change the radial part and the angular part of the wave function.

In 1928 Dirac investigated the relativistic wave equation covariance of the Schrodinger equation and proposed a matrix *α **β **I *^{[ 1 ] }The Dirac equation describes the motion of a relativistic particle with spin 1/2, which is widely used to solve the problems of nuclear physics and high energy physics. ^{[ 2 ] }One of the important tasks of relativistic quantum mechanics is to find an accurate and exact solution of the Dirac equation for some noncentral physical potentials, which are a mixture of attractive scalar potential *S *( *r *) and repulsive vector potential *V *( * r *). Spin symmetry occurs when

*S*(

*r*) =

*V*(

*r*

*S*(

*r*) = −

*V*(

*r*

^{[ 3 ] }

The Dirac equation for the case of exact spin symmetric occurs when the difference in the magnitude betwen the repulsive vector potential and the attractive scalar potential is zero and the sum of the magnitudes of the repulsive vector potential and attractive scalar potential is equal to a given potential. The exact pseudo-spin symmetry occurs when the sum of the magnitude of the repulsive vector potential and the attractive scalar potential is zero and the difference between the vector potential and scalar potential is equal to a given potential, which is central or non-central. ^{[ 4 ] }

Solutions of the Dirac equation for some typical potentials under special cases of spin symmetry and pseudospin symmetry have been investigated. For spin symmetry *κ *( *κ *+ 1) = *l *( *l *+ 1) that gives *κ *= *l *= *j *+ 1/2 for *κ *> 0 and *κ *= −( *l *+ 1) = −( *j *+ 1/2) for *κ *< 0, where *κ *is the eigenvalue of the spin orbit coupling operator, *l *is orbital quantum number, and *j *is the total angular momentum quantum number. For the pseudospin symmetry case, *κ *( *κ *− 1) = *l̃ *( *l̃ *+ 1) that gives *κ *= *l̃ *+ 1 = *j *+ 1/2 for *κ *> 0 and *κ *= − *l̃ *= −( *j *+ 1/2) for *κ *< 0, where *l̃ *is the quantum number of the pseudospin orbital. In nuclear physics, spin symmetry and pseudospin symmetry concepts have been used to study the aspect of deformed and super deformation nuclei. The concept of spin symmetry has been applied to the spectra of meson and antinucleon. ^{[ 5 ] }

Some researchers have investigated the Dirac equation by using a variety of potentials and different methods, such as the spin symmetry in the antinucleon spectrum and tensortype Coulomb potential with spin–orbit number *κ *in a state of spin symmetry and p-spin symmetry, ^{[ 6 ] }bound states of the Dirac equation with position-dependent mass for the Eckart potential, ^{[ 7 ] }the exact solution of Klein–Gordon with the Pöschl–Teller double-ring-shaped Coulomb potential, ^{[ 8 ] }the exact solution of the Dirac equation for the Coulomb potential plus NAD potential by using the Nikorov–Uvarov method, ^{[ 9 ] }the potential Deng–Fan and the Coulomb potential tensor using the asymptotic iteration method (AIM), ^{[ 10 ] }the potential Poschl–Teller plus the Manning Rosen radial section with the hypergeometry method, ^{[ 11 ] }the solution of Klein–Gordon equation for Hulthen non-central potential in radial part with Romanovski polynomial ^{[ 12 ] }and the solution of the Schrodinger equation with the Hulthen plus Manning– Rosen potential, ^{[ 13 ] }the Scarf potential with the new tensor coupling potential for spin and pseudospin symmetries using Romanovski polynomials, ^{[ 14 ] }for the q-deformed hyperbolic Poschl–Teller potential and the trigonometric Scarf II noncentral potential by using AIM, ^{[ 15 ] }eigensolutions of the deformed Woods–Saxon potential via AIM, ^{[ 16 ] }approximate solutions of the Klein Gordon equation with an improved Manning Rosen potential in *D *-dimensions using SUSYQM, ^{[ 17 ] }and eigen spectra of the Dirac equation for a deformed Woods–Saxon potential via the similarity transformation. ^{[ 18 ] }

In this study, we will solve the Dirac equation for the Eckart potential plus the trigonometric Manning Rosen potential in a state of spin symmetry using an asymptotic iteration method. The asymptotic iteration method is a method of solving the second-order differential equations. ^{[ 20 ] }

The rest of this paper is organized as follows. The asymptotic iteration method will be briefly reviewed in Section 2. In Section 3, we review the Eckart potential and trigonometric Manning Rosen potential, give a brief introduction to the Dirac equation with equal scalar and vector potentials, and apply the separation of variables in spherical coordinates. In Section 4, we solve the radial and angular parts of the Dirac equation by using the modified Eckart potential combined with trigonometric Manning Rosen potential and obtain the relativistic energy spectrum and wavefunction via the asymptotic iteration method. In Section 5, we present graphically some wavefunctions of the Dirac equation, present several relativistic energy spectra, and discuss some consequences of the results obtained. Finally, some conclusions are drawn from the present study in Section 6.

This method is used to solve differential equation in the following form:

*λ*

_{k }and

*s*

_{k },

^{[ 19 ] }

^{[ 19 ] }

*α*can be solved by applying the asymptotic iteration

*k*

*y*

_{n }(

*x*) is given by

*C′*is the normalized constant, and

_{2 }F

_{1 }is the hipergeometry function.

^{[ 20 ] }

The Dirac equation with scalar potential *S *( *r *) and vector potential *V *( *r **c *= 1) is given by

*E*is the relativistic energy,

*p*

*p*

*α*

*β*

*I*

*σ*

^{[ 14 ] }whose components are expressed as

*f*

_{nκ }(

*r*

*g*

_{nκ }(

*r*

*F*

_{n r }(

*r*) is the Dirac spin component,

*G*

_{n r }(

*r*) is the Dirac pseudospin component,

*l*is the orbital quantum number,

*l̄*is the quantum number of the pseudo orbital, and

*m*is the projection of the angular momentum on the

*z*axis.

By substituting Eq. (

*F*

_{n r }(

*r*) and

*G*

_{n r }(

*r*) in Eqs. (

*F*

_{n r }(

*r*) and

*G*

_{n r }(

*r*) components of the Dirac wavefunction as follows:

*κ*is the quantum number of the spin orbital,

*l*is the orbital quantum number, and

*l̃*the quantum number of the pseudospin orbital.

In the case of exact spin symmetry d *Δ *( *r *)/d *r *= * Δ *(

*r*

*C*= const, equation (

*Σ*(

*r*) is the Eckart potensial, which is defined as

*αr*) = 1 − 2

*z*into Eq. (

*B*

_{s }+

*E*

_{s }= 4

*δ*

^{2 }and

*B*

_{s }−

*E*

_{s }= 4

*γ*

^{2 }.

By comparing this equation with Eq. ( *λ *_{0 }( *z *) and *s *_{0 }( *z *) values. By means of Eq. ( *λ *_{k }( *z *) and *s *_{k }( *z *)

*ɛ*

_{n r }=

*A*

_{s }/4.

When the above expressions are generalized into

*E*as follows:

*n*

_{r }is the radial quantum number (

*n*

_{r }= 0,1,2…),

*l*is orbital quantum number which is obtained from the angular part solution. The radial wavefunction can be obtained by using Eqs. (

*z*= [1-coth (

*αr*)]/2, so

*C*

_{n r }is the radial normalization constant,

_{2 }F

_{1 }is the hypergeometry function, and (2

*δ*+ 1)

_{n r }is a Pochamer symbol.

For the angular part in Eq. ( *P *( *θ *)/sin ( *θ *), we have

*l*being the orbital quantum number.

Equation ( *z *(1 − *z *), then we have

*λ*

_{k }and

*s*

_{k }, where

*k*is the stated iteration, by using Eq. (

When the above expressions are generalized

*Δ*

_{k }, which yields

*l*is the orbital quantum number and

*n*

_{l }is the angular quantum number.

Angular wavefunction can be obtained by using Eqs. (

*z*= [1 + i cot (

*θ*)]/2, so the angular part wavefunction can be obtained as follows:

*C*

_{n l }being the angular normalization constant.

In this section, we discuss several results obtained in the previous section. From the energy eigen value in Eq. ( *v *= 0.25, *α *= 0.5, *V *_{0 }= 0.4, *V *_{1 }= 0.6, and *M *= 5 fm ^{−1 }. The positive value is taken due to the spin symmetric limit. Table *q *and *n *_{r }in the same quantum state cause the energy eigenvalue to decrease.

Table *V *_{0 }and *V *_{1 }, which are plotted in Fig. *V *_{0 }and *V *_{1 }. Figure *V *_{0 }and *V *_{1 }cause the energy to decrease.

By varying parameters corresponding to values of *δ *and *γ *, some of the radial wavefunctions are listed in Table *V *_{0 }and *V *_{1 }, and by parameter *q *as well. The parameter *q *has a dimension of inversed distance in space, which describes the effect range of the trigonometric Manning Rosen potential in space. If *q *is enlarged, it physically means that the potential range is smaller in a space. Table *n *_{r }for the radial wavefunction induces the distribution of particles to move away from or approach to the core of the atom. From Fig. *n *_{r }= 1 < *n *_{r }= 2, which means that at *n *_{r }= 2, the distribution of particles moves away from the nucleus compared with the case of smaller *n *_{r }. This indicates that the higher the value of *n *_{r }is, the greater the distance of the particle from the core is.

For the angular solution of the wavefunction in Eq. ( *n *_{l }= 4 but different values of *α *and *β *which corespond to parameters *V *_{0 }and *V *_{1 }in Fig. *α *< *β *; in Figs. *α *> *β *. For the same value of *n *_{l }, the sum of waves does not change.

In this paper, we study the Dirac equation for particle spin-1/2 in the Eckart potential combined with the trigonometric Manning Rosen potential in condition that the scalar potential equals the vector potential. The radial part of the spinor wavefunction is obtained approximately from Eq. ( *l *as expressed in Eq. ( *n *_{r }causes the energy to decrease.

**Reference**

1 | |

2 | Hindawi Publishing Corporation Advances in High Energy Physics 2013 987632 |

3 | |

4 | ITB Journal Publisher 46 205 |

5 | AIP Conf. Proc. 1615 111 |

6 | |

7 | Chin. Phys. B 22 010301 |

8 | Acta Phys. Polon. A 126 647 |

9 | |

10 | Indian J. Phys. 88 405 |

11 | |

12 | |

13 | |

14 | |

15 | |

16 | |

17 | |

18 | |

19 | |

20 |