#### Abstract

We solve approximately Dirac equation for Eckart plus Hulthen potentials with Coulomb-like and Yukawa-like tensor interaction in the presence of spin and pseudospin symmetry for . The formula method is used to obtain the energy eigenvalues and wave functions. We also discuss the energy eigenvalues and the Dirac spinors for Eckart plus Hulthen potentials with formula method. To show the accuracy of the present model, some numerical results are shown in both pseudospin and spin symmetry limits.

#### 1. Introduction

One of the interesting problems in nuclear and high energy physics is to obtain analytical solution of the Klein-Gordon, Duffin-Kemmer-Petiau, and Dirac equations for mixed vector and scalar potentials [1]. The study of relativistic effect is always useful in some quantum mechanical systems [2, 3]. Therefore, the Dirac equation has become the most appealing relativistic wave equation for spin-1/2 particles. For example, in the relativistic treatment of nuclear phenomena, the Dirac equation is used to describe the behavior of the nuclei in nucleus and also to solve many problems of high energy physics and chemistry. For this reason, it has been used extensively to study the relativistic heavy ion collisions and heavy ion spectroscopy and more recently in laser-matter interaction (for a review, see [4] and references therein) and condensed matter physics [5, 6]. The idea about spin symmetry and pseudospin symmetry with the nuclear shell model has been introduced in 1969 by Arima et al. (1969) and Hecht and Adler (1969) [7, 8]. Spin and pseudospin symmetries are symmetries of a Dirac Hamiltonian with vector and scalar potentials. They are realized when the difference, , or the sum, , is constant. The near realization of these symmetries may explain degeneracy in some heavy meson spectra (spin symmetry) or in single-particle energy levels in nuclei (pseudospin symmetry), when these physical systems are described by relativistic mean-field theories (RMF) with scalar and vector potentials [9]. Recently, some authors have studied various types of potential with a tensor potential, under the conditions of pseudospin and spin symmetry. They have found out that the tensor interaction removes the degeneracy between two states in the pseudospin and spin doublet [10–13]. The pseudospin and spin symmetry appearing in nuclear physics refers to a quasidegeneracy of the single-nucleon doublets and can be characterized with the nonrelativistic quantum numbers () and (), where , , and are the single-nucleon radial, orbital, and total angular momentum quantum numbers for a single particle, respectively [7, 8]. These kinds of various methods have been used for the analytical solutions of the Klein-Gordon equation and Dirac equation such as the super symmetric quantum mechanics [14–16], asymptotic iteration method (AIM) [17, 18], factorization method [19, 20], Laplace transform approach [21], GPS method [22, 23] and the path integral method [24–26], and Nikiforov-Uvarov method [27–29]. The Klein-Gordon and Dirac wave equations are frequently used to describe the particle dynamics in relativistic quantum mechanics with some typical kinds of potential by using different methods [30]. For example, Kratzer potential [31, 32], Woods-Saxon potential [33, 34], Scarf potential [35, 36], Hartmann potential [37, 38], Rosen Morse potential, [39, 40], Hulthen potential [41], and Eckart potential [42, 43].

In this paper, we attempt to solve approximately Dirac wave equation for for Eckart plus Hulthen potentials for the spin and pseudospin symmetry with a tensor potential by using the formula method. The organization of this paper is as follows: in Section 2, the formula method is reviewed [44]. In Section 3 we review basic Dirac equations briefly. In Sections 3.1 and 3.2, solutions of Dirac wave equation for the spin and pseudospin symmetry of these potentials in the presence of Coulomb-like tensor interaction are presented, respectively. In Sections 3.3 and 3.4, solutions of Dirac wave equation for the spin and pseudospin symmetry of these potentials in the presence of Coulomb-like plus Yukawa-like tensor interaction are presented, respectively. In Section 4, we provide results and discussion. The conclusion is given in Section 5.

#### 2. Review of Formula Method

The formula method has been used to solve the Schrodinger, Dirac, and Klein-Gordon wave equations for a certain kind of potential. In this method, the differential equations can be written as follows [44]:For a given Schrödinger-like equation in the presence of any potential model which can be written in the form of (1), the energy eigenvalues and the corresponding wave function can be obtained by using the following formulas, respectively [44]:where And is the normalization constant. In special case where , the energy eigenvalues and the corresponding wave function can be obtained as [44]The solution provides a valuable means for checking and improving models and numerical methods introduced for solving complicated quantum systems.

#### 3. Basic Dirac Equations

In the relativistic description, the Dirac equation of a single nucleon with the mass moving in an attractive scalar potential and a repulsive vector potential can be written as [45]where is the relativistic energy, is the mass of a single particle, and and are the 4 × 4 Dirac matrices. For a particle in a central field, the total angular momentum and commute with the Dirac Hamiltonian, where is the orbital angular momentum. For a given total angular momentum , the eigenvalues of the are , where negative sign is for aligned spin and positive sign is for unaligned spin. The wave functions can be classified according to their angular momentum and spin-orbit quantum number as follows:where and are upper and lower components and and are the spherical harmonic functions. is the radial quantum number and is the projection of the angular momentum on the -axis. The orbital angular momentum quantum numbers and represent the spin and pseudospin quantum numbers. Substituting (7) into (6), we obtain couple equations for the radial part of the Dirac equation as follows by :where and are the difference and the sum of the potentials and , respectively, and is a tensor potential. We obtain the second-order Schrodinger-like equation asWe consider bound state solutions that demand the radial components satisfying , and [45].

##### 3.1. Solution Spin Symmetric with Coulomb-Like Tensor Interaction

Under the condition of the spin symmetry—that is, or —the upper component Dirac equation can be written asThe potential is taken as the Eckart [42, 43] plus Hulthen potentials [41]: where the parameters , , , and are real parameters; these parameters describe the depth of the potential well, and the parameter is related to the range of the potential.

For the tensor term, we consider the Coulomb-like potential [46],where is the Coulomb radius and and stand for the charges of the projectile particle and the target nucleus , respectively.

By substituting (11) and (12) into (10), we obtain the upper radial equation of Dirac equation aswhere , and .

Equation (13) is exactly solvable only for the case of . In order to obtain the analytical solutions of (13), we employ the improved approximation scheme suggested by Greene and Aldrich [47, 48] and replace the spin-orbit coupling term with the expression that is valid for [49]:The behavior of the improved approximation is plotted in Figure 1. We can see the good agreement for small values.

By applying the transformation , (13) is brought into the formwhere the parameters , , and are considered as follows:Now by comparing (15) with (1), we can easily obtain the coefficients () as follows:The values of the coefficients () are also found from (4) as below:Thus, by the use of energy equation (2) for energy eigenvalues, we findIn Tables 1–3, we give the numerical results for the spin symmetric energy eigenvalues (in units of fm^{−1}).

Let us find the corresponding wave functions. In reference to (3) and (18), we can obtain the upper wave function aswhere is the normalization constant; on the other hand, the lower component of the Dirac spinor can be calculated fromWe have obtained the energy eigenvalues and the wave function of the radial Dirac equation for Eckart plus Hulthen potentials with Coulomb-like tensor interaction in the presence of the spin symmetry for .

##### 3.2. Solution Pseudospin Symmetry with Coulomb-Like Tensor Interaction

For pseudospin symmetry—that is, or —the lower component Dirac equation can be written asWe consider the scalar, vector, and tensor potentials as the following [41]: Substituting (23) into (22), we obtain the lower radial equation of Dirac equation aswhere , and .

By using the transformation and employing the improved approximation, (24) is brought into the formwhere the parameters , , and are considered as follows:We can easily obtain the coefficients () by comparing (25) with (1) as follows:The values of the coefficients () are also found from (4) as below:We have (29) for energy eigenvalues by the use of (2):In Tables 4–6, we give the numerical results for the pseudospin symmetric energy eigenvalues (in units of fm^{−1}).

By using (3) and (28), we can finally obtain the lower component of the Dirac spinor as below:where is the normalization constant; on the other hand, the upper component of the Dirac spinor can be calculated from (30) asWe have obtained the energy eigenvalues and the spinors of the radial Dirac equation for Eckart plus Hulthen potentials with Coulomb-like tensor interaction in the presence of the pseudospin symmetry for . We show in Figure 2 behavior energy for various in the spin and pseudospin symmetry.

**(a)**

**(b)**

##### 3.3. Solution Spin Symmetry with Coulomb Plus Yukawa-Like Tensor Interaction

In this section, for the spin symmetry, we consider , , and as the following:where and are the real parameters; substitution of (32) into (10) yieldswhere the parameters , , and are considered as follows:By comparing (33) with (1), we can easily obtain the coefficients () as follows:The values of the coefficients () are also found from (4) as below:By the use of energy equation (2) for energy eigenvalues, we haveIn reference to (3) and using (36), we can obtain the upper wave function:where is the normalization constant; on the other hand, the upper component of the Dirac spinor can be calculated from (38) asThe effects of the Yukawa-like tensor interactions on the upper and lower components for the spin symmetry are shown in Figure 3.

**(a)**

**(b)**

We have obtained the energy eigenvalues and the spinors of the radial Dirac equation for Eckart plus Hulthen potentials with the spin symmetry for in the presence and absence of Yukawa tensor.

##### 3.4. Solution Pseudospin Symmetry with Coulomb Plus Yukawa-Like Tensor Interaction

In this section, for the pseudospin symmetry, we consider , , and as the following: where and are the real parameters. Substitution of (40) into (22) yieldswhere the parameters , , and are considered as follows:By comparing (41) with (1), we can easily obtain the coefficients () as follows:The values of the coefficients () are also found from (4) as below:By using energy equation (2) for energy eigenvalues, we haveBy the use of (3) and (44), we can finally obtain the lower component of the Dirac spinor as below: where is the normalization constant; on the other hand, the upper component of the Dirac spinor can be calculated from (46) asThe effects of the Yukawa-like tensor interactions on the upper and lower components for the pseudospin symmetry are shown in Figure 4.

**(a)**

**(b)**

We have obtained the energy eigenvalues and the spinors of the radial Dirac equation for Eckart plus Hulthen potentials with the pseudospin symmetry for in the presence and absence of Yukawa tensor interaction. In Figures 5 and 6, we show the effect of Coulomb-like tensor interaction and Yukawa plus Coulomb-like tensor interaction in the remove of degeneracy.

**(a)**

**(b)**

**(a)**

**(b)**

In Figure 7, we show the comparison between Coulomb and Yukawa-like tensor interaction for spin and pseudospin symmetries.

**(a)**

**(b)**

In Figures 5, 6, and 7, we showed that degeneracy is removed by tensor interaction. Furthermore, the amount of the energy difference between the two states in the doublets increases with increasing and .

#### 4. Results and Discussion

We obtained the energy eigenvalues in the absence and the presence of the Coulomb-like tensor potential for various values of the quantum numbers and . In Tables 1 and 4, in the absence of the tensor interaction (), the degeneracy between spin doublets and pseudospin doublets is observed. For example, we observe the degeneracy in (, ), (, ), and so forth in the spin symmetry, and we observe the degeneracy in (, ), (, ), and so forth in the pseudospin symmetry. When we consider the tensor interaction, for example, by parameter , the degeneracy is removed. In Tables 2 and 3 for the spin symmetry and also in Tables 5 and 6 for the pseudospin symmetry, we show that degeneracy exists between spin doublets for several of the parameters and , and we show that degeneracy is removed in the presence of tensor interaction. In Figure 2, the degeneracy is removed by tensor interaction effect in spin symmetry and pseudospin symmetry; also the amount of the energy difference between the two states in the doublets increases with increasing parameter . The effects of the Yukawa-like tensor interactions on the upper and lower components of radial Dirac equation for the symmetries are shown in Figures 3 and 4. The sensitiveness of the pseudospin doublets (, ) and (, ) and spin doublets (, ) and (, ) for the effects of the Yukawa-like tensor interactions and for Yukawa plus Coulomb-like tensor interaction is given in Figures 5 and 6. In Figure 7 we show that the Coulomb-tensor interaction is stronger than Yukawa-like tensor interaction to remove degeneracy.

#### 5. Conclusions

In this paper, we have discussed approximately the solutions of the Dirac equation for Eckart plus Hulthen potentials with spin symmetry and pseudospin symmetry for . We obtained the energy eigenvalues and the wave function in terms of the generalized polynomials functions via the formula method. To show the accuracy of the present model, some numerical values of the energy levels are shown in Figures 3, 4, 5, and 6. We have showed that the energy degeneracy in pseudospin and spin doublets is removed by the tensor interaction effect.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.