Research Article | Open Access
Dirac Particle for the Position Dependent Mass in the Generalized Asymmetric Woods-Saxon Potential
The one-dimensional Dirac equation with position dependent mass in the generalized asymmetric Woods-Saxon potential is solved in terms of the hypergeometric functions. The transmission and reflection coefficients are obtained by considering the one-dimensional electric current density for the Dirac particle and the equation describing the bound states is found by utilizing the continuity conditions of the obtained wave function. Also, by using the generalized asymmetric Woods-Saxon potential solutions, the scattering states are found out without making calculation for the Woods-Saxon, Hulthen, cusp potentials, and so forth, which are derived from the generalized asymmetric Woods-Saxon potential and the conditions describing transmission resonances and supercriticality are achieved. At the same time, the data obtained in this work are compared with the results achieved in earlier studies and are observed to be consistent.
The scattering and bound states in nonrelativistic and relativistic quantum mechanics within an external potential are studied in order to describe the behavior of particles, atoms, and molecules in physics. Thus, the nonrelativistic and relativistic particle equations in an external potential have been widely handled for the scattering and bound states [1–16]. In the recent years, these external potentials have been generalized to asymmetric potentials in theoretical and experimental physics [13–22]. Also, last year, an international team at the ISOLDE radioactive-beam facility at CERN reported that atomic nuclei could suppose asymmetric “pear” shapes .
Another important issue is the concept of position dependent mass. This concept is studied in different applications of physics, for example, electronic properties of semiconductor materials , quantum liquids , quantum dots , and the nuclear forces of relativistic two nucleons [26, 27].
In this paper, the generalized asymmetric Woods-Saxon potential (GAWSp), that includes the Woods-Saxon potential, is considered to obtain the solutions of the Dirac equation with position dependent mass. The Woods-Saxon potential  has been widely used in the nonlinear theory of scalar mesons [29, 30], the nuclear shell model, and the distribution of nuclear densities [31, 32]. The GAWSp is introduced in  where is the Heaviside step function and is potential intensity and real. Parameters , , , , , , , and are real. The GAWSp is shown in Figure 1 and in special cases, the GAWSp turns into the potentials which are used in the nonrelativistic and relativistic quantum physics. Some special cases of the GAWSp are shown in Table 1.
The paper is organized as follows. In Section 2, the transmission and reflection coefficients are found for the effective mass Dirac particle. In Section 3, an equation is obtained for energy eigenvalues. In Section 4, the results and discussions and the plots constructed by using Mathematica Software are commented and the energy values are given numerically; moreover, the transmission resonance condition for the scattering states and supercriticality are derived and the data obtained in this work are compared with the results achieved in earlier studies. Finally, the present paper is summarized in Section 5.
2. Scattering Solutions of the GAWSp for the Variable Mass Dirac Particle
The Dirac equation for the relativistic free-particle is given as the following form  (in natural units ): where , the Dirac particle mass, depends on one spatially coordinate . In order to obtain the one-dimensional Dirac equation under an external potential , the gamma matrices and are reduced to the Pauli matrices and , respectively, where and are decomposed into upper and lower components of the two-component wave function and (3) turns into the two coupled differential equations as follows: Then the following two expressions are written as they were described by Flügge as  Putting these expressions into (4), the following equations are achieved: By using the two equations given above, for and , the following two independently second-order differential equations are obtained: The Dirac particle mass is supposed in  where the function is given by . The parameter is the rest mass of Dirac particle and is a real positive and small parameter; therefore, the derivative terms including the mass function are ignored in (8). So, (8) becomes It is clearly seen from (6) and (7) that it will be sufficient to solve one of (10) and (11) to obtain the solutions describing the scattering states. Therefore, we start calculation by solving (10). Let us now define a new variable in (10) for the region , in this instance, (10) becomes and getting a trial wave function , (12) reduces to the Gaussian differential equation  where The solution of (13) is found in the form of hypergeometric functions, as follows: so, the whole left-hand solution becomes Now, we can make the process to find the solution for the region . Choosing a new variable , (10) is written as and assuming a test wave function , the solution of (17) becomes the following form: In this case, the whole right-hand solution is where In order to investigate the transmission () and reflection () coefficients, the solutions obtained above have to be used asymptotic behaviors as and . If is chosen, the left-hand solution can be written as and if is chosen, the right-hand solution is The one-dimensional electric current density for the Dirac particle is given as the following form: so as to find electric current density, (21) and (22) are inserted into (6). In this case, the asymptotic behaviors for and become Using (21), (22), and (24) for electric current density given in (23), , , and are calculated and the reflection and transmission coefficients are found as follows, respectively: where , , and represent the incident, transmitted, and reflected currents, respectively. Finally, to investigate numerically the reflection and transmission coefficients, the continuity condition and the first derivative of the wave function are used. So, and , taking part in the above expressions, are found as where the constants are given in Table 2.
3. The Equation of Energy Eigenvalues
The purpose of this section is to derive an equation that is known as the bound state condition for the energy eigenvalues. For this, is put instead of in the potential given in (1) and the new variable is defined as variable in the region . As a result, (10) is written in the following form: Taking a wave function as , the solution of (27) is obtained as where Similarly, choosing the variable and wave function in the region , the positive region solution is found as where Because of the boundary conditions (which mean the wave functions go to zero at infinity), the wave functions given in (28) and (30) are written as In order to find an equation for energy eigenvalues, the continuity conditions of the wave function are used, and . After performing the necessary calculations, an equation for energy eigenvalues is obtained as follows: where the constants are given Table 3.
4. Results and Discussions
4.1. Interpretation of the Drawings and Comparison of the Results with Previous Ones
In Figure 2, while the parameter goes to zero and the GAWSp is a symmetric potential, one can clearly see that the unitary condition, , is provided in the GAWSp for the position dependent mass Dirac particle. On the other hand, if the parameter does not go to zero and the GAWSp is an asymmetric one, the unitary condition is not provided. Figure 3 shows the effect of position dependent mass for the Dirac particle within the GAWSp on the transmission coefficient. In this figure, it is seen that the presence of variable mass causes the formation of a wider line according to the previous graph that is one of the constant masses. One of the advantages of this work is that, by using the GAWSp solutions for the Dirac particle with variable mass, the transmission and reflection coefficients are directly obtained without doing calculation for the potentials which are given in Table 1. For example, in the numerical results of the GAWSp, if one chooses , , , and , one can acquire the scattering states of the Woods-Saxon potential, by putting , , , , and , the scattering states of the Hulthen potential are obtained, substituting , , , and , the scattering states of the cusp potential are found. Figures 4–6 are drawn using the GAWSp solutions for the Dirac particle. Figure 4 displays a plot of the transmission coefficient describing the behavior of the Dirac particle in the Woods-Saxon potential versus the energy and the potential strength for both the position dependent mass and the constant mass. From Figure 4, it is seen, in the presence of position dependent mass, that the zero region of the transmission coefficient expands. At the same time, the results of Figure 4 are compatible with [7, 20]. Figure 5 shows the effect of the position dependent Dirac particle for the Hulthen potential on the transmission coefficient. The data drawn from Figure 5 prove that the values of the transmission coefficient in the case of position dependent mass decrease. Furthermore, these results are similar to [12, 22]. Figure 6 displays the effect of the Dirac particle for the cusp potential on the transmission coefficient. Figure 7 shows a plot of the transmission coefficients obtained for the Dirac particle with the position dependent mass against the energy and the potential strength for the GAWS, the Woods-Saxon, the Hulthen, and the cusp potentials. In the left panel which is located in this figure, while the transmission resonance appears for the GAWS and Woods-Saxon potentials at lower values of energy, it does not exist for the Hulthen and the cusp potentials at same values of energy. In the right panel taking part in the same figure, while the transmission resonance peaks are more for the Woods-Saxon potential at higher values of potential strength (), the transmission resonance does not occur for the cusp potential at same values of . Moreover, if one wants to achieve the scattering solutions of the asymmetric Hulthen and the asymmetric cusp potentials, one can use the solutions describing the scattering states, of the GAWSp, found in this study.