#### Abstract

We obtain exact scattering solutions of the Dirac equation in 1 + 1 dimensions for a double square barrier vector potential. The potential bottom between the two barriers is chosen to be higher than 2*mc*^{2}, whereas the top of the barriers is at least 2*mc*^{2} above the bottom. The relativistic version of the conventional double barrier transmission resonances is obtained for energies within ±*mc*^{2} from the height of the barriers. However, due to our judicious choice of potential configuration we also find two more (subbarrier) transmission resonance regions below the conventional one. Both are located within the two Klein energy zones and characterized by resonances that are broader than the conventional ones. The design of our double barrier so as to enable us to establish these two new subbarrier transmission resonance regions is our main finding.

#### 1. Introduction

The basic equation of relativistic quantum mechanics was formulated more than 80 years ago by Paul Dirac [1–3]. It describes the state of electrons in a way consistent with special relativity, requiring that electrons have spin and predicting the existence of an antiparticle partner to the electron (the positron). The physics and mathematics of the Dirac equation is very rich, illuminating, and provides a theoretical framework for different physical phenomena that are not present in the nonrelativistic regime such as the Klein paradox, supercriticality (supercritical transmission of a relativistic particle through a potential barrier) [1–7], and the anomalous quantum Hall effect in graphene [8, 9]. It is well known that the Dirac equation has positive as well as negative energy solutions [1–3]. The positive and negative energy subspaces are completely disconnected. This is a general feature of the solution space of the Dirac equation, which is sometimes overlooked. Since the equation is linear, then the complete solution must be a linear combination of the two. Incorporating the missing part of the negative energy solution, which is not taken into account in the traditional solution, one of the authors (A. D. Alhaidari) gave a new approach to the resolution of the famous Klein paradox within relativistic quantum mechanics [10]. This approach leads to the correct physical and mathematical interpretations of this phenomenon without resort to quantum field theory.

On the other hand, tunneling phenomena played an important role in nonrelativistic quantum mechanics due to its important application in electronic devices [11, 12]. It was Esaki who discovered a characteristic called negative differential resistance (NDR) whereby, for PN junction diodes, the current voltage characteristics has a sharp peak at a certain voltage associated with resonant tunneling. This constituted the first important confirmation that this phenomenon is due to the quantum mechanical tunneling effect of electrons [13–15]. Tunneling is a purely quantum phenomenon that happens in the classically forbidden region; its experimental observation constituted a very important support to the quantum theory. On the other hand, the study of tunneling of relativistic particles through one-dimensional potentials has been restricted to some simple configurations such as * δ*-potentials and square barriers, mainly, in the study of the possible relativistic corrections to mesoscopic conduction [16] and the analysis of resonant tunneling through multibarrier systems [17]. Relativistic studies of the quark and Dirac particles in 1D periodic potential were also reported separately in [18, 19].

Very recently, electron transport through electrostatic barriers in single and bilayer graphene has been studied using the Dirac equation, and barrier penetration effects analogous to the Klein paradox were noted [20]. The study of transmission resonances in relativistic wave equations in external potentials has been discussed extensively in the literature [21, 22]. In this case, for given values of the energy and shape of the barrier, the probability of transmission reaches unity even if the potential strength is larger than the energy of the particle, a phenomenon that is not present in the nonrelativistic case. The relation between low momentum resonances and supercritical transmission has been established by Dombey and Calogeracos [7] and Kennedy [23]. Some results on the scattering of Dirac particles by a one-dimensional potential exhibiting resonant behavior have also been reported [22–24].

However, recent studies have shown that inducing a finite bandgap in graphene by epitaxially growing it on a substrate is possible [25] and, therefore, its energy dispersion relation is no longer linear in momentum. This process generated gaps in graphene energy spectrum and resulted in a finite effective mass for its charge carriers and opened up nanoelectronic opportunities for graphene. This 2D massless system can then be mapped into an effectively massive 1D one [26] and, consequently, the problem of Klein tunneling of Dirac fermions across a potential can be put on the test, because the potential barriers can be seen as n-p-n junctions of graphene if they are high enough. Under these circumstances, it is interesting to characterize the system behavior by determining the full expressions of the corresponding reflection and transmission coefficients.

Motivated by the above progress, we study in this work the resonant transmission of a beam of relativistic particles through two separated square barriers with elevated potential bottom in between and investigate transmission resonance in this structure [1–3]. Under special conditions, dictated by our judicious choice of potential bottom elevations and barrier heights, we demonstrate the occurrence of three (subbarrier) regions of transmission resonances. One of them is the relativistic extension of the conventional nonrelativistic double barrier transmission resonances for energies within from the height of the barriers. The other two are located within the two Klein energy zones where only positive and negative energy oscillatory solutions coexist at the same energy. The latter resonances are broader than the conventional ones.

#### 2. Scattering Solution of the Dirac Equation

The physical configuration associated with the double barrier problem in our study is shown in Figure 1. In the relativistic units , the one-dimensional stationary Dirac equation with vector potential coupling can be written as [1–3]
where *V*(*x*) is the time component of the vector potential whose space component vanishes (i.e., gauged away due to gauge invariance). The potential *V*(*x*) is defined by
where and are positive potential parameters (see Figure 1) such that and . We divide configuration space according to the piece-wise constant potential sections into three regions numbered 0 and ± corresponding to and , respectively. In regions 0, where the potential vanishes, the equation becomes the free Dirac equation that relates the two spinor components as follows:
This relationship is valid for . Since the problem is linear and because belongs to the tive energy spectrum, then (2.3) with the top/bottom sign is valid only for positive/negative energy, respectively. After choosing a sign in (2.3) then the other spinor component obeys the following Schrödinger-like second order differential equation:
We should emphasize that (2.4) does not give the two components of the spinor that belong to the same energy subspace. One has to choose one sign in (2.4) to obtain only one of the two components then substitute that into (2.3) with the corresponding sign to obtain the other component. Now, within the double barrier (the regions) the same analysis follows but with the substitution giving
where stands for either of the two potentials . Generally, in any region of constant potential *V*, positive/negative energy solutions occur for relativistic energies larger/smaller than . Of these, the oscillatory solutions of the form hold for , where . On the other hand, the exponential solutions of the form hold for .

The scattering solution, which is the subject of this work, pertains to energies . It is straightforward to write down the positive and negative energy solutions of (2.3)–(2.5). First, we write the wave vector associated with regions of space in which the potential equals to zero, , and , as
where and . This results in oscillatory solutions if is pure imaginary which happens when or (i.e., in the two grey regions of Figure 1). Otherwise, these solutions are exponentials (i.e., in the white areas of the figure). The oscillatory positive/negative energy solutions are located in the light/dark grey areas of Figure 1, respectively. We divide configuration space from left to right into five regions indexed by . The general positive energy solution in these regions (both oscillatory and exponentials) can be written as
whereas the negative energy solutions are of the form
and are constants (the complex amplitudes) associated with right and left “traveling” solutions in the * ν*th region, respectively. The energy parameters and are defined byNote that for real/imaginary values, respectively. The complex constant amplitudes will be determined by the boundary conditions. We should note that the oscillatory solutions, , in (2.7a) and (2.7b) represent a wave traveling in the ±

*x*direction for positive energy solutions and in the

*x*direction for negative energy solutions. The solution of the Dirac equation to the right of the double barrier consists of positive energy plane-wave solutions traveling in the ±

*x*directions. However, the physical boundary conditions of the problem allow only transmitted waves traveling to the right after passing through the double barrier (i.e., ). Moreover, and without loss of generality, we can normalize the incident beam to unit amplitude (i.e., ).

Matching the spinor wavefunctions at the four boundaries defined by and gives relations between in * ν*th region and those in the neighboring region. We prefer to express these relationships in terms of 2 × 2 transfer matrices between different regions, , with . Finally, we obtain the full transfer matrix over the whole double barrier which can be written, in an obvious notation, as follows:
where and we have set and ;

*R*and being the reflection and transmission amplitudes, respectively. We have assumed an incident wave from left normalized to unit amplitude (i.e., and ). The explicit form of the transfer matrices depends on the specific energy range. There are three such ranges for all . These ranges are (i) , (ii) , and (iii) . Therefore, we end up with the full set of twelve transfer matrices given in the Appendix. Equation (2.9) leads to Time reversal invariance and the relevant conservation laws dictate that the transfer matrix

*M*(

*E*) has a unit determinant along with the following symmetry properties (

*E*) = (

*E*)* and (

*E*) = (

*E*)*, where * stands for complex conjugation. These can easily be checked using the explicit forms given in the Appendix. Thus, symmetry considerations impose strong conditions on the structure of the transfer matrix. Using these properties in (2.10) gives the expected flux conservation . Moreover, from (2.10) we see that full transmission or resonance transmission occurs at energies where the condition is satisfied (equivalently, ).

#### 3. Results and Discussion

The physical content of particle scattering through the double barrier depends on the energy of the incoming particle, which can assume any value larger than . In order to allow for supercritical transmission through the potential, we need to impose certain conditions on the heights of the potential barriers. Our study concentrates on Klein energy zones where full transmission can take place. The situation of interest to our study concerns two Klein energy zones, which arises when and . As an example, we calculate the transmission coefficient as a function of energy for a given set of potential parameters. The result is shown in Figures 2 and 3. We combine Figures 1 and 2 together in such a way that Figure 2 is turned 90 degrees such that the energy axis of Figures 1 and 2 coincide. In this way, one can clearly see how the different transmission regions correspond to the different potential regions. In addition to the expected above-barrier full transmission for some values of energies larger than , one can clearly identify three subbarrier regions where transmission resonances occur. These are(i) the lower Klein energy zone (): a region of seven resonances,(ii) the higher Klein energy zone (): a region of four resonances,(iii) the relativistic version of the conventional nonrelativistic double barrier transmission resonance (): a region of four resonances.It is also clear that resonances in the conventional energy zone are very sharp (or very narrow) whereas those in the two Klein energy zones are broad (or wide). This means that resonance states corresponding to the former decay are much slower and have longer tunneling time than the latter. In Table 1, we list the resonance energies for this potential configuration to an accuracy of 10 decimal places. These results were obtained as solution of the equation . As further insight into the dynamics of this relativistic model, Figure 4 shows an animation of Figure 2 as the distance between the two barriers, , varies from to . The animation shows the following.(i) The density of resonances in each of the three subbarrier regions increases with , that is, the energy separation between resonances decreases with .(ii) As increases, resonance energies drop down (fall or dive) from the above-barrier region into the conventional resonance region then into the higher Klein energy zone.(iii) Additionally, as increases, resonance energies are created at the bottom of the spectrum (at ) then move up into the lower Klein energy zone.

Figure 4 gives another animation of Figure 2 as the width of the barriers, , varies from to . The animation shows that all resonances get sharper with an increase in but the number of resonances in the conventional region does not change (i.e., the population density of resonances in this region is independent of ). We would like to mention a related recent work by Villalba and Gonzalez-Arraga [27] who considered the resonant tunneling through a double square barrier and double cusp potentials. Our problem differs from that in [27] by the choice of an elevated bottom of the potential well, which gives rise to two Klein energy zones of resonance. This potential design gave rise to a peculiar energy dependence of the transmission with three resonance regions; one is due to the conventional quantum tunneling and two others are due to Klein tunneling. Another animation of Figure 2 is given in Figure 5.

For completeness, it is also interesting to consider the potential barrier with the configuration shown in Figure 6, where and . The associated transmission coefficient as a function of energy is shown in Figure 7 for the given potential parameters where the values of were interchanged. In Figure 8, we combine Figures 6 and 7 together in such a way that Figure 7 is turned 90 degrees such that the energy axis of Figures 6 and 7 coincide. In this way, one can clearly see how the different transmission regions correspond to the different potential regions. It looks as if the energy region of Figure 2 was flipped and placed in the energy region of Figure 7. Consequently, the sharp resonance region became sandwiched between the two Klein energy zones while the transmission structure of the higher Klein energy zone was reversed.

Finally, we note that the present work will not remain at this stage but will be pursued further. We plan to use the results obtained so far to deal with different issues related to transport properties in graphene. One of the main characteristics of Dirac fermions in graphene is the accuracy with which we can model their behavior by having extremely small mass (in fact, even massless). This implies that at any finite energy the model should be treated relativistically. This endows fermions in graphene with the ability to tunnel through a single potential barrier with probability one [20, 28–31]. It is then natural to extend that analysis to our two-barrier problem case and investigate the basic features of such a system. However, we would like to mention that graphene is a two-dimensional (2D) system, and 2D carrier tunneling through 1D barriers can be very complicated and direction dependent. Our present results correspond to the transmission of 2D carriers only when the carriers move perpendicular to the potential barriers.

#### Appendix

#### Transfer Matrices

If we define , and , then in the first energy interval, , the four transfer matrices at the boundaries and are given by

Then, in the second energy range, , we have

Note that and have the same form as in (A.1). Finally, in the energy range, , we obtain

#### Acknowledgments

This work is sponsored by the Saudi Center for Theoretical Physics (SCTP). The authors also acknowledge the support provided by King Fahd University of Petroleum & Minerals under Group Projects RG1108-1 and RG1108-2.