By using the Pekeris approximation, we present solutions of the Dirac equation with the generalized Woods-Saxon potential with arbitrary spin-orbit coupling number under spin symmetry limit. We obtain energy eigenvalues and corresponding eigenfunctions in closed forms. Some numerical results are given too.

1. Introduction

The Dirac equation, which describes the motion of a spin 1/2 particle, has been used in solving many problems of nuclear and high-energy physics. One of the most important concepts for understanding the traditional magic numbers in stable nuclei is the spin symmetry breaking [14]. On the other hand, the -spin symmetry observed originally almost 40 years ago as a mechanism to explain different aspects of the nuclear structure is one of the most interesting phenomena in the relativistic quantum mechanics. Within the framework of Dirac equation, -spin symmetry used to feature deformed nuclei, superdeformation, and to establish an effective shell model [58], whereas spin symmetry is relevant for mesons [9, 10]. Spin symmetry occurs when and pseudospin symmetry occurs when , [11, 12], where is scalar potential and is vector potential. Pseudospin symmetry is exact under the condition of , and spin symmetry is exact under the condition of [13]. For the first time, the spin symmetry tested in the realistic nuclei [14], and then this symmetry is investigated by examining the radial wave functions [1517]. The pseudospin symmetry refers to a quasi-degeneracy of single nucleon doublets with nonrelativistic quantum numbers and , where , , and are single nucleon radial, orbital, and total angular quantum numbers, respectively. The total angular momentum is , where is pseudo-angular momentum and is pseudospin angular momentum [12, 1820].

On the other hand, some typical physical models have been studied like harmonic oscillator [19, 20], Woods-Saxon potential [21, 22], Morse potential [2325], Eckart potential [2628], Pöschl-Teller potential [29], Manning-Rosen potential [30], and so forth [3136].

The interactions between nuclei are commonly described by using a potential that consists of the Coulomb and the nuclear potentials. These potentials are usually taken to be of the Woods-Saxon form. The Coulomb plus Woods-Saxon potentials are well known as modified Woods-Saxon potential that plays a great role in nuclear physics. Fusion barriers for a large number of fusion reactions from light to heavy systems can be described well with this potential [3740]. But the modified Woods-Saxon potential has no exact or approximate solutions. The generalized Woods-Saxon is near the modified Woods-Saxon potential, and therefore we can solve the generalized Woods-Saxon instead of the modified Woods-Saxon potential. The form of the generalized Woods-Saxon potential is as follows [41, 42]: where and determine the potential depth, is the width of the potential, is the surface thickness, and  MeV [42]. In Figure 1, we illustrated the shape of the Woods-Saxon potential () and its generalized form in (1) ().

The purpose of this work is devoted to studying the spin symmetry solutions of the Dirac equation for arbitrary spin-orbit coupling quantum number with the generalized Woods-Saxon potential, which was not considered before.

This paper is organized as follows. In Section 2, we briefly introduced the Dirac equation with scalar and vector potential with arbitrary spin-orbit coupling number under spin symmetry limit. The generalized parametric Nikiforov-Uvarov method is presented in Section 3. The energy eigenvalue equations and corresponding eigenfunctions are obtained in Section 4. In this section, some remarks and numerical results are given too. Finally, conclusion is given in Section 5.

2. Dirac Equation under Spin and Pseudospin Symmetry Limit

The Dirac equation with scalar potential and vector potential is where is the relativistic energy of the system and is the three-dimensional momentum operator. and are the usual Dirac matrices given as where is Pauli matrix and is unitary matrix. The total angular momentum operator and spin-orbit operator , where is orbital angular momentum operator and commutes with Dirac Hamiltonian. The eigenvalues of spin-orbit coupling operator are and for unaligned spin and the aligned spin , respectively. can be taken as a complete set of the conservative quantities. Thus, the Dirac spinors can be written according to radial quantum number and spin-orbit coupling number as follows: where is the upper (large) component and is the lower (small) component of the Dirac spinors. and are the spherical harmonic functions coupled to the total angular momentum and its projection on the axis. Substituting (4) into (2) with the usual Dirac matrices, one obtains two coupled differential equations for the upper and the lower radial wave functions and as


Eliminating and from (5a) and (5b), we obtain the following second-order Schrödinger-like differential equations for the upper and lower components of the Dirac wave functions, respectively:

where and .

2.1. Spin Symmetry Limit

In the spin symmetry limit or   constant [14], then (7b) becomes where and for and , respectively. In (8), can be taken as generalized Woods-Saxon potential, and it is reduced to where

For the solution of (9), we will use the Nikiforov-Uvarov method which is briefly introduced in the following section.

2.2. Pekeris-Type Approximation to the Spin-Orbit Coupling Term

Because of the spin-orbit coupling term, that is, , (9) cannot be solved analytically, except for . Therefore, we shall use the Pekeris approximation [43] in order to deal with the spin-orbit coupling terms, and we may express the spin-orbit term as follows: In addition, we may also approximately express it in the following way: where , , and is constant () [4448]. If we expand the expression of (12) around up to the second-order term and next compare it with (11), we can obtain expansion coefficients , , and as follows: Now, we can take the potential (12) instead of the spin-orbit coupling potential (11).

3. The Parametric Generalization Nikiforov-Uvarov Method

This powerful mathematical tool solves second-order differential equations. Let us consider the following differential equation [4951]: According to the Nikiforov-Uvarov method, the eigenfunctions and eigenenergies, respectively, are where In the rather more special case of [50, 51],

and, from (15), we find for the wave function

4. Bound States of the Generalized Woods-Saxon Potential with Arbitrary

Substituting (12) into (9) and using transformation , we find the following second order deferential equation for the upper component of the Dirac spinor as Comparing the previous equation with (14), one can find the following parameters: and also From (16), (21), and (22), we obtained the closed form of the energy eigenvalues for the generalized Woods-Saxon potential in spin symmetry as Recalling and from (10a) and (10b) the previous equation becomes a quadratic algebraic equation in . Thus, the solution of this algebraic equation with respect to can be obtained in terms of particular values of and . In Table 1, we have calculated some energy levels of the generalized Woods-Saxon potential under the spin symmetry limit for different values of . The empirical values that can be found in [52], as  fm,  fm,  u, and (MeV), are used. Here, is the mass number of target nucleus, and . From Table 1, we can see that pairs ,, , , and so forth are degenerate. Thus, each pair is considered as spin doublet and has negative energy, and also when increases, the energy decreases. The reason is that when increases, the depth of the potential well becomes deeper and then the energy levels decrease. In Figure 2, the results are shown as a function of . As we saw in Table 1, when the depth of the potential well increases, then the energy levels decrease. In Figure 3, the results are presented as a function of width of the potential . Here, when the width of the potential increases, the energy levels decrease. Finally we showed the numerical results as a function of surface thickness in Figure 4, and one can observe that when the surface thickness increases, the energy levels decrease too.

To find corresponding wave functions, referring to (15), (21), and (22), we find the upper component of the Dirac spinor as or equivalently Finally, the lower component of the Dirac spinor can be calculated as where [20].

5. Conclusion

In this paper we have studied the spin symmetry of a Dirac nucleon subjected to scalar and vector generalized Woods-Saxon potentials. The quadratic energy equation and spinor wave functions for bound states have been obtained by parametric form of the Nikiforov-Uvarov method. It is shown that there exist negative-energy bound states in the case of exact spin symmetry (). We gave some numerical results of the energy eigenvalues too.


The authors thank the kind referee for positive and invaluable suggestions, which improved this paper greatly.