arxiv: v1 [cond-mat.mes-hall] 13 Dec 2018

Transcription

1 Transmission in Graphene through Time-oscillating inear Barrier l Bouâzzaoui Choubabi a, Ahmed Jellal a,b, and Miloud Mekkaoui a a aboratory of Theoretical Physics, Faculty of Sciences, Chouaïb Doukkali University, arxiv:8.55v [cond-mat.mes-hall] 3 Dec 8 PO Box, l Jadida, Morocco b Saudi Center for Theoretical Physics, Dhahran, Saudi Arabia Abstract Transmission probabilities of Dirac fermions in graphene under linear barrier potential oscillating in time is investigated. Solving Dirac equation we end up with the solutions of the energy spectrum depending on several modes coming from the oscillations. These will be used to obtain a transfer matrix that allows to determine transmission amplitudes of all modes. Due to numerical difficulties in truncating the resulting coupled channel equations, we limit ourselves to low quantum channels, i.e. l =, ±, and study the three corresponding transmission probabilities. PACS numbers: b; b; 7.8.Rj Keywords: Graphene, time-oscillating linear potential, Dirac equation, transmission modes. a.jellal@ucd.ac.ma

2 Introduction Graphene is a stable planar monolayer of densely crystallized carbon atoms in a two-dimensional honeycomb lattice []. Since its realization in, graphene has generated an immense interest in study of their mechanical, electronic, optical, thermal and chemical properties [, 3, ]. The electronic properties are justified by his Hamiltonian, relativistic approximated in Tight-Binding model framework, at low energy to the neighborhoods of the six Dirac points at the Brillouin zone corners in the reciprocal lattice. These are represented by two inequivalent points K and K corresponding to the two atoms A and B constituting the pattern in the direct lattice. In addition, the Dirac fermions in graphene behave as chiral massless particles with a "light speed" equal to the crystal velocity Fermi v F c/3 with a gapless linear dispersion relation near Dirac points. This allows graphene to be a candidate for manufacturing the carbon-based nanoelectronic devices. Quantum transport in periodically driven quantum systems is an important subject not only of academic value but also for device and optical applications. In particular quantum interference within an oscillating time-periodic electromagnetic field gives rise to additional sidebands at energies ɛ + l ω l =, ±, in the transmission probability originating from the fact that electrons exchange energy quanta ω carried by photons of the oscillating field, ω being the frequency of the oscillating field. The standard model in this context is that of a time-modulated scalar potential in a finite region of space. It was studied earlier by Dayem and Martin [5] who provided the experimental evidence of photon assisted tunneling in experiments on superconducting films under microwave fields. ater on Tien and Gordon [6] provided the first theoretical explanation of these experimental observations. Further theoretical studies were performed later by many research groups, in particular Buttiker investigated the barrier traversal time of particles interacting with a time-oscillating barrier [7]. Wagner [8] gave a detailed treatment on photon-assisted tunneling through a strongly driven double barrier tunneling diode. He also studied the transmission probability of electrons traversing a quantum well subject to a harmonic driving force [9] where transmission side-bands have been predicted. Grossmann [], on the other hand, investigated the tunneling through a double-well perturbed by a monochromatic driving force which gave rise to unexpected modifications in the tunneling phenomenon. In our previous work [], we have solved the D Dirac equation describing graphene in the presence of a linear vector potential. The discretization of the transverse momentum due to the infinite mass boundary condition reduced our D Dirac equation to an effective massive D Dirac equation with an effective mass equal to the quantized transverse momentum. We have used both a numerical Poincaré map approach, based on space discretization of the original Dirac equation, and a direct analytical method. These two approaches have been used to study tunneling phenomena through a biased graphene strip. We have showed that the numerical results generated by the Poincaré map are in complete agreement with the analytical results. In the second one [], the energy spectrum of a graphene sheet subject to a single barrier potential having a time periodic oscillating height and subject to a magnetic field is analyzed. The corresponding transmission is studied as function of the incident energy and potential parameters. Quantum interference within the oscillating barrier has an important effect on quasiparticles tunneling. In particular the time-periodic electrostatic potential generates additional sidebands at energies in the transmission probability originating from the photon absorption or emission within the oscillating barrier.

3 Based on our previous work [, ], we consider a graphene sheet subjected to oscillating barrier potential along the x-direction while the carriers are free in the y-direction. In fact, the barrier height V of linear barrier potential oscillates sinusoidally with oscillation amplitude U and frequency ω. After getting the energy spectrum, we match all wavefunctions to end with a transfer matrix that allow to determine the transmission probability. Since we have many energy modes due to the barrier oscillations in time, we focus only on three channels and describe numerically their behaviors in terms of different physical parameters of the present system. The manuscript is organized as follows. In section, we present the theoretical model describing a graphene sheet in the presence of time-oscillating linear barrier potential. Subsequently, we explicitly determine the eigenvalues and eigenspinors for each region composing our system. In section 3, the scattering problem for Dirac fermions will be solved using the boundary conditions, which will be used to calculate the corresponding transmission probabilities for all energy modes. In section, we numerically discuss our results by making appropriate choices of the physical parameters characterizing our system. In the final section, we conclude our results. Theoretical formulation Our system is composed of three regions, successively, see Figure, where in regions I and III there is only pristine graphene, while the intermediate region II subjected to in the same time to time-harmonic potential of amplitude U and pulsation ω as well as linear potential generated by the two flat armatures sited at interfaces x = and x =. The armatures are identical, with the same surface spaced by a distance, and are perpendicular to the graphene sheet. The condenser two flat armatures is biased by a potential difference V V. We assume that the dimensional requirement allows the condenser to be considered as infinite in order to neglect the edge effects. In the condenser, the electric field prevailing between the two armatures comes from superposition of the electric fields created by each armature, the resultant field remains perpendicular to the armatures. Since the field spatially derives of a potential so the equipotential surfaces are parallel to the armatures and vary gradually linear between the two potentials reported in the armatures. We assume that the graphene sheet passes through grooves at the centers of armatures without disturbing the system Figure in such a way that the intermediate region of the sheet is subjected to a linear inter-armatures potential given by Figure V V V x = x + V = Sx + V, V > V. Our system is now made of massless Dirac fermions move along the x-direction passing through region II subjected to time-oscillating linear barrier V x, t = V x + U cosωt Figure, but they are free in the y-direction. We assume that the graphene sheet is characterized by very large length scale x-direction and width W y-direction. In input region I, at the interface x =, the incident fermions of energy and angle φ with respect to the x-direction are reflected with energies + m ω m =, ±, ±, and incident angles π φ m. In output region III, at the interface x =, after transmission the fermions have the same energies with transmission angle φ m.

4 Figure Color online Schematic of a graphene sheet whose region II subjected to a time-oscillating linear barrier together potential. Figure Color online Configuration of time-oscillating linear barrier potential on graphene sheet composed of of three regions. The Hamiltonian describing our system can be written as H = i v F σ + V x, ti where υ F is the Fermi velocity, σ = σ x, σ y are the usual Pauli matrices and I the unit matrix hereafter we use v F = =. Our system has finite width W with infinite mass boundary condition on the wavefunction at boundaries y = and y = W along the y-direction [3, ]. These result in a quantization of the transverse momentum = π W n +, n =,,. 3 On other hand, the separability requires to write eigenspinors as ψx, y, t = ψx, yψt and then 3

5 from the eigenvalue equation Hψ j x, y, t = i t ψ j x, y, t we obtain [ + U cos ωt] ψ j x, y, t = i t ψ j x, y, t which can be integrated to end up with the solution ψ j x, y, t = ψ j x, y, e it e iu sinωt/ω. 5 By introducing the expansion in terms of the m-th order Bessel function of the first kind J m, we can write the eigenspinors as ψ j x, y, t = ψ j x, y, e i t m= J m αe imωt 6 which correspond to the eigenvalues = + mω 7 where we have set α = U ω. Taking into account of the energy conservation, the eigenspinors that describe our fermions in the j-th region can be expressed as a linear combination of those at energies + lω l =, ±, ±,. Thus we have m,l=+ ψ j x, y, t = e ikyy m,l= ψ l jx, yj m l α e i+mωt 8 where the spinor ψj l x, y will be determined according to each region j. Indeed, in region I, by solving the time-dependent Dirac equation we can easily get the solution at energy for the incoming fermions ψ inc x, y, t = k + i α = s k + α e ik x e ikyy e it 9 = s e iφ where s = sgn, φ is the angle that the incident fermions make with the x-direction, k and are the x and y-components of the fermion wave vector, respectively. Reflected and transmitted waves have components at all energies + lω l =, ±, ±,. Then the eigenspinors ψ r x, y, t for reflected fermions are ψ r x, y, t = m,l= k l + i α l = s l kl + r l α l e ik lx+iy J m l α e i+mωt = s l e iφ l and correspond to the eigenvalues + lω = s l k l + 3

6 where r l is the reflection amplitude and J m l α = δ m,l because the modulation amplitude is V j = in this case. φ l = tan /k l, sign s l = sgn + lω refers to the conduction and valence bands of region I. The number k l can be obtained from 3 k l = s l + lω ky. Finally in region I x <, we have the eigenspinors ψ I x, y, t = e ikyy m,l= [δ l, αl e ik lx + r l α l ] e ik lx δ m,l e i+mωt. 5 For region II < x < a linear combination of eigenspinors at energies +lω l =, ±, ±, has to be taken and therefore we have ψ II x, y, t = l= ψ l x, y m= J m α e imωt i+lωt. 6 These can be written in terms of the parabolic cylinder function [] such that the first component is given by φ + II = c D ν Q l + c D ν Q l 7 where ɛ l = + lω V, ν = i S, Q lx = S eiπ/ Sx + ɛ l, c and c are constants. The second component takes the form φ II = c [ ɛ l + SxD ν Q l + ] Se iπ/ D ν+ Q l c Se iπ/ D ν Q l. 8 The components of the spinor solution of the Dirac equation in region II can be obtained from 7 and 8, such as ψ II x, y, t = φ + II + iφ II φ + II iφ II = e ikyy l= c,l η + l x η l x where the function η ± l x and ξ± l x are given by η ± l x = D ν Q l + c,l ξ + l x ξ l x J m l α e i+mωt m= 9 S e iπ/ D ν Q l ξ ± l x = ± Se iπ/ D ν+ Q l ± iɛ l ± isxd ν Q l. For region III x >, the eigenspinors ψ t x, y, t for transmitted fermions read as ψ t x, y, t = m,l= t l α l e ik lx+iy J m l α e i+mωt 5

7 where t l is the transmission amplitude, which can mapped in terms of null vector {b l } as ψ III x, y, t = e ikyy m,l= [t l αl e ik lx + b l α l e ik lx ] δ m,l e i+mωt. 3 The obtained solutions of the energy spectrum will be used to deal with different issues. More precisely, these results will be employed to determine the transmission probabilities for all energy modes. 3 Transmission through oscillating linear barrier Using the orthogonality of {e imωt } and the continuity of eigenspinors at boundaries x =, x =, to end up with a set of equations for r m, t m, b m, c,m, c,m δ m, + r m = l= δ m, α m r m α m = c,l η l + + c,lξ l + J m l α l= t m e ikm + b m e ikm = c,l ηl + c,lξl J m l α 5 l= t m α m e ikm b m α m e ikm = c,l η l + + c,lξ l + J m l α 6 l= c,l ηl c,lξl J m l α. 7 Since Dirac fermions pass through a region subjected to time-oscillating linear potential, transitions from the central band to sidebands channels at energies + mω m =, ±, ±, occur as fermions ex-change energy quanta with the oscillating field. expressed in the transfer matrix formalism, such as Ξ Ξ = M M M M Ξ Ξ = M This procedure is most conveniently Ξ where the total transfer matrix M = M, M, and transfer matrices Mj, j +, that couple the wave function in the j-th region to the wave function in the j + -th region, are Ξ 8 I I C + G + M, = N + N C + M, = G + C G C G I I N + N K + O O K 9 3 by setting the quantities N ± = ± α m ± δ m,l, m,l G ± τ m,l = ξ± l τj m l α, C ± τ m,l = η± l τj m l α 3 K ± m,l = ±e±ikm δ m,l 3 6

8 where the null matrix is denoted by O and I is the unit matrix. We are considering fermions propagating from left to right with energy, then, τ =,, Ξ = {δ,l } and Ξ = {b m} is the null vector, whereas Ξ = {t l } and Ξ = {r l} are the vectors of transmitted and reflected waves respectively. From the above considerations, one can easily obtain the relation Ξ = M Ξ. 33 The minimum number N of sidebands that need to be considered is determined by the strength of the oscillation, N > U ω. Then the infinite series for T can be truncated to consider a finite number of terms starting from N up to N. Furthermore, analytical results are obtained if we consider small values of α = U ω and include only the first two sidebands at energies ± ω along with the central band at energy. This gives the result where k =,,,, N and M t N+k = M [k +, N + ] 3 denotes the inverse matrix M. Then for α =, it remains only the transmission t = M [, ] for central bands, which gives exactly the result obtained in our previous work [] t = e ik [ + z ] [ ξ + η ξ η + ] [ξ + + z ξ ] [η z η + ] [η + + z η ] [ξ z ξ + ]. 35 On the other hand, we can use the reflected J ref and transmitted J tra currents to write the transmission T l and reflection R l probabilities corresponding to our system. These are T l = J tra,l J inc,, R l = J ref,l J inc, 36 where T l is describing the scattering of fermions with incident energy in the region I into the sideband with energy + lω in region III. Thus, the rank of the transfer matrix M increases with the amplitude of the time-oscillating potential. From the Hamiltonian system, we derive the current density J = v F ψ σ x ψ 37 which gives the incident, reflected and transmitted currents J inc, = v F α + α 38 J ref,l = v F rl r l α l + α l 39 J tra,l = v F t l t l α l + α l Then from the above results, we derive the transmission and reflection probabilities T l = k l k t l, R l = k l k r l Due to numerical difficulties, we truncate 8 retaining only terms corresponding to the central and 7

9 first sidebands, namely l =,,, which allow to write t = M [, ], t = M [, ], t = M [3, ]. To explore the above results and go deeply in order to underline our system behavior, we will introduce the numerical analysis. For this, we will focus only on few channels and choose different configurations of the physical parameters. Numerical results In Figure 3, we present the transmission probability versus incident energy. Figure 3a shows the general behavior of the transmission probability of a linear barrier α = limited by two values of potential V =, V = for = 3 and =, respectively. We observe that the transmission is forbidden below the energy = ky =, which behaves like an effective mass ky = m [5, 6]. For energies ky V ky the transmission shows oscillations and gives rise to the Klein paradox, i.e. T α = = Klein zone. In the range of energy V V +, the transmission decreases from = V until a minimum value where it oscillates and after increases by passing from = V to reach a maximal value at = V +. However for energy > V + our system shows a classical behavior. To give comparison and show the relevance of our finding, we plots two cases according to choice of the barrier heights V, V. Indeed, Figure 3b illustrates a particular case of a linear barrier V =, V = studied in our previous work [] where transmission corresponding to the Klein zone is omitted and transmission oscillates around a minimum, then it behaves in the same way as shown in 3a. Figure 3c presents the case of a simple square barrier V V = 3 where the Klein zone is conserved and transmission corresponding to energies V V + is replaced by another in range V V + without oscillations [7]. Finally, we observe that Figures 3 highlight the effect of the barrier width on transmission, as long as increases the minimal transmission of the intermediate zones decreases and the number of oscillations increases. a b c.75 V.75 V.75 V V Tα=.5.5 V Tα=.5.5 Tα=.5.5 V 3 V 5 6 V 3 5 V V 5 6 Figure 3 Color onlinetransmission probability of three barrier potentials as function of incident energy, with =, α =, = 3 blue color and = red color for a: V = and V =, b: V = and V = and c: V V = 3. In Figures we present the three channel transmissions together with total one summation over three channels versus incident energy for α =.3, α =.6 and α =.9. These concern the central transmission band T blue color, two first sidebands T l= red color, T l= green color, and total 8

10 T magenta color T = l=,, k l k t l. 3 Since the time-oscillating linear barrier vibrates sinusoidally and longitudinally along the propagating x-direction of Dirac fermions, with amplitude U and frequency ω, we notice that this causes a change in the effective mass from to ± ω. For low values of α the central transmission band T is dominant compared to those of lateral a. As long as α increases the central transmission band l = decreases while the two first sidebands l = ± increase. In Figure b we have co-dominance of the three bands but in Figure c the two first sidebands l = ± dominate. We observe that the total transmission decreases if α increases. Since α = U ω of T might be caused by the missing modes. is barrier height dependent, the decreasing a b c T, T-,, α= V U cosω t V T, T-,, α= T, T-,, α= V 3 V 5 6 V 3 V 5 6 V 3 V 5 6 Figure Color online Transmission probabilities of time-oscillating linear barrier potential as a function of incident energy, with =, = 3, V =, V = and ω = for a: α =.3, b: α =.6 and c: α =.9. T blue color, T l= red color, T l= green color, T magenta color. Figure 5 shows transmission probabilities T,, α =.8, T α = versus barrier width under the conditions indicated in its caption. All transmission probabilities have a periodicity of = π, except T α =. The transmissions in Figure 5a, corresponding to energy < < V, maintain sinusoidal oscillations, while the two sidebands transmission T, α =.8 are almost equal throughout the periodicity of. For [kπ, kπ + a[ ]k + π a, k + π], k N and a [, π] see 5a, T α =.8 who dominates, for ]kπ+a, k+π a[ of T, α =.8 who dominates, however for = kπ+a or = k+π a the three transmissions T,, α =.8 are co-dominant. Throughout the variation of, T α = alternately takes values close to unity. Figure 5b presents T,, α =.8, T α = for = V = where the transmissions are waning in oscillation as long as increases. For [kπ, kπ + a[ ]k + π a, k + π] we have the dominance of T α =.8 until = 8π. Beyond this value, the transmission T α =.8 dominates T α =.8 and the two transmissions are maintained with the increase of, but T α =.8 vanishes at = π. The transmission probability T α = decays sinusoidally and vanishes at = π. In Figure 5c V < < V, all transmission probabilities decay and vanish before = 5π. The transmission probabilities in Figure 5d = V = are similar to those of Figure 5b = V =, except that the roles of T α =.8 and T α =.8 are inverted. Above = V =, from 5d = we see that all transmissions begin to behave in similar way as in 5a. The transmission probabilities T,, α =.8, T α = for = 5 > V in Figure 5e are similar to those of the Figure 5a whose energy is = 5 < V and the difference is that those of Figure 5e are smooth and thin 9

11 compared to those in the Figure 5a. a b c T-,,α=.8, T α= T-,,α=.8, T α= T-,,α=.8, T α= a π π 6π 8π π π π π 6π 8π π π π π 6π 8π π π d e f T-,,α=.8, T α= T-,,α=.8, T α= T-,,α=.8, T α= π π 6π 8π π π π π 6π 8π π π π π 6π 8π π π Figure 5 Color online ffect of barrier width on time-periodic linear barrier transmission probabilities in terms of, with =, V =, V =, α =,.8, T α = magenta color, T α =.5 blew color, T α =.5 red color, and T α =.5 green color, for a: = 5, b: =, c: = 6, d: =, e: =, f: = 5. Now we fix the barrier width as = π, π and vary to show its effect on transmissions probabilities for the three modes. For this, we divide the interval { of the wave vectors into two regions where 6a, 6b, 6d, 6e < V correspond to regions R =], [ and R =], [ with } { : dashed red line, = : dashed black line and 6c, 6f correspond to regions R 3 =], 3 [ } and R =]3, [ with 3 < : dashed orange line, = : dashed magenta line. Figure 6 shows that for any value of energy, we have always the constraint <. It is clearly seen that in R and R 3 there is manifestation of all transmission modes, but in R the T α =.8 mode is missing, and in R all modes T,, α =.8, T α = are forbidden. In R 6a, T α = is the first dominant, followed by T, α =.8, which are co-dominant at the beginning and finish by dominance of T α =.8 with respect to T α =.8. Finally the central band T α = dominates T α =.8 only just before. In R 6a, T α =, T α =.8, T α =.8 show peaks one peak, two peaks, two peaks respectively, the first peak of T α = dominates the first peak of T α =.8 that also dominates the first peak of T α =.8 and the second peak of T α =.8 dominates the second peak of T α =.8. At the beginning of R 6d, T α =, T α =.8 are co-dominant between them and dominate together the two transmissions T, α =.8 which are also co-dominant between them. After =, each of the four transmissions generates a first peak, where that of T α = is greater than that of T α =.8, which also greater than that of T, α =.8. The transmissions in R 6d behave in the same way as those in R 6a, except towards the end of region one has a generation of a supplementary peak for each transmission with T α = > T α =.8 > T α =.8. We notice that in R 6b, R 3 6c, R 6e, and R 3 6f, each transmission probability has a single peak. We observe that there are T α = > T α =.8 > T α =.8 > T α =.8

12 in R 6b, T α = > T α =.8 > T α =.8 > T α =.8 in R 3 6c, T α = = T α =.8 > T α =.8 T α =.8 in R 6e and R 3 6f. a b c T-,,α=.8, T α= T-,,α=.8, T α= T-,,α=.8, T α= d e f T-,,α=.8, T α= T-,,α=.8, T α= T-,,α=.8, T α= Figure 6 Color online Transmission probabilities versus wave vector with α =,.8, V =, V =, T α = magenta color, T α =.5 blew color, T α =.5 red color, and T α =.5 green color, for a: = 5 and = π, b: = and = π, c: = and = π, d: = 5 and = π, e: = and = π, f: = and = π. 5 Conclusion We have considered Dirac fermions in graphene sheet subjected to a linear potential barrier oscillating harmonically in time along x-direction but free along y-direction. In fact, we have applied between two of sheet a potential barrier type V x = V V x + V + U cosωt, which is oscillating sinusoidally with an amplitude U and frequency ω. By solving the eigenvalue equation, we have obtained the solutions of the energy spectrum in terms of different physical parameters involved in the Hamiltonian of system. Moreover, we have showed that the barrier in time generated an infinite modes giving rise to several energy modes. Subsequently, we have studied the transport properties of our system through transmission probabilities. Indeed, using the transfer matrix method, we have calculated the corresponding transitions and reflections coefficients. These were used to show that how the transmission probabilities for the central band l = and two first sidebands l = ± are affected by various physical parameters such that incident energy, barrier width and transverse wave vector. In summary, our numerical results support the assertion that quantum interference has an important effect on fermions through graphene based linear barrier driven by time-oscillations. Since most optical applications in electronic devices are based on interference phenomena then we expect that the results of our computations might be of interest to designers of graphene-based electronic devices.

13 Acknowledgments The generous support provided by the Saudi Center for Theoretical Physics SCTP is highly appreciated by all authors. References [] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, Science 36, 666. [] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and A. A. Firsov, Nature 38, [3] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, Rev. Mod. Phys. 8, 9 9. [] C. W. Beenakker, Rev. Mod. Phys. 8, [5] A. H. Dayem and R. J. Martin, Phys. Rev. ett. 8, [6] P. K. Tien and J. P. Gordon, Phys. Rev. 9, [7] M. Moskalets and M. Buttiker, Phys. Rev. B 66, [8] M. Wagner, Phys. Rev. A 5, [9] M. Wagner, Phys. Rev. B 9, [] F. Grossmann, T. Dittrich, P. Jung, and P. Hanggi, Phys. Rev. ett. 67, [] H. Bahlouli,. B. Choubabi, A. l Mouhafid, and A. Jellal, Solid State Communications 5, 39. [] A. Jellal, M. Mekkaoui,. B. Choubabi, and H. Bahlouli, ur. Phys. J. B 87, 3. [3] J. Tworzydlo, B. Trauzettel, M. Titov, A. Rycerz, and C. W. J. Beenakker, Phys. Rev. ett. 96, [] M. V. Berry and R. J. Modragon, Proc. R. Soc. ondon Ser. A, [5] A. Jellal,. B. Choubabi, H. Bahlouli, and A. Ajaafari, J ow Temp. Phys. 68,. [6] A. D. Alhaidari, A. Jellal,. B. Choubabi, and H. Bahlouli, Quntum matter, 3. [7] H. Bahlouli,. B. Choubabi, A. Jellal, and M. Mekkaoui, J ow Temp. Phys. 59, 5.