Abstract
We investigate quantum tunneling through a single electric and/or magnetic barrier on the surface of a threedimensional topological insulator. We found that (1) the propagating behavior of electrons in such system exhibits a strong dependence on the direction of the incident electron wavevector and incident energy, giving the possibility to construct a wave vector and/or energy filter; (2) the spin orientation can be tuned by changing the magnetic barrier structure as well as the incident angles and energies.
PACS numbers: 72.25.Dc; 73.20.r; 73.23.b; 75.70.i.
1. Introduction
The recent discovery of a new quantum state of matter, topological insulator, has generated a lot of interest due to its great scientific and technological importance [15]. In a topological insulator, spinorbit coupling opens an energy gap in the bulk, and results in helical surface states residing in the bulk gap in the absence of magnetic fields. Such surface states are spindependent and are topologically protected by timereversal symmetry [47] and distinct from conventional surface states, which are fragile and depend sensitively on the details of the surface geometry and bonding. This discovery sparked intensive experimental and theoretical interests, both for its fundamental novel electronics properties as well as possible applications in a new generation of electric devices.
Very recently, the surface states in Bibased alloys, Bi_{1x}Sb_{x}, Bi_{2}Se_{3}, Bi_{2}Te_{3}, were theoretically predicted [6,8] and experimentally observed by using angleresolved photoemission spectroscopy (ARPES) [912]. These 3D topological insulators have robust and simple surface states consisting of a single Dirac cone at the Γ point [8]. Note that this Hamiltonian appears similar to graphene [13], but topological insulators have an odd number of massless Dirac cones on the surface, ensured by the Z_{2 }topological invariant of the bulk, while graphene has twofold massless Dirac cones at the K and K' valleys. Another essential difference is connected with spinrelated properties. In the surface, Hamiltonian of the 3D TI σ acts on the real spin of the charge carriers, while for graphene it stands for the pseudo spin, i.e., the A and B sublattices of graphene. Hence, it is natural to manipulate spin transport on the surface of a 3D topological insulator by controlling the electron orbital motion. Based on the topological surface Hamiltonian, it is clear that σ · k is a quantum conserved quantity which implies that spin and momentum of the electron are locked. For instance in a tunneling process, the reflected electron will reverse its spin due to the helical property of the surface states, i.e., the spinmomentum locking [14]. This feature will lead to some interesting phenomena, such as the spindependent conductance [15] and the twisted RKKY interaction [16].
In this study, we investigate electron tunneling through single electric and magnetic potential barriers which can be created by depositing a ferromagnetic metallic strip on the surface of a 3D topological insulator. We find that the inplane spin orientation of the transmitted and the reflected electrons can be rotated over certain angles that are determined by the incident angle and energy. Our results demonstrate that the magnetic field of the magnetic barrier bends the trajectory of the electrons, and therefore rotate the spin.
2. Theory
First we focus on the electron transmission through a single electric and magnetic potential barrier on the surface of a 3D topological insulator, as shown in Figure 1. We can create magnetic and electric potentials underneath a superconducting plate and a ferromagnetic strip [1720]. The relevant magnetic and electric fields are directed perpendicular to the surface of the 3D TI, i.e., V(x) = V[Θ(x)  Θ(x  d)]/2 and B_{z}(x) = B[Θ(x)  Θ(x  d)]/2, where Θ is the Heaviside step function. The magnetic field is perpendicular to the surface and the Landau gauge A = (0, Bx, 0) is adopted in our calculation. The lowenergy electrons near the Γ point of the Dirac cones can be well described by the effective Hamiltonian [21],
Figure 1. Schematic of an electric and magnetic potential barrier on the surface of a 3D TI.
where v_{F }is the Fermi velocity, σ^{i}(i = x, y, z) are the Pauli matrices, V is the gate voltage applied on the magnetic metal strips, and the last term H_{Z }≡ gμ_{B}σ · B is induced by Zeeman spin spitting. Note that for g = 23 in Bi_{2}Se_{3}, the Zeeman term affects the transmission slightly at low magnetic field. The momentum is π = p + eA, where the vector potential of the inhomogeneous magnetic field generated by the magnetic metal stripe, , is the unit vector normal to the surface. Note that , which implies that is a quantum conserved quantity during the tunneling processes, i.e., spinmomentum locking. If the incident electrons are spin polarized along the direction of the vector , the magnetic field bends the trajectory of the electrons, resulting in a rotation of the spin of the transmitted electrons. The reflected electrons suffer similar spin rotations accompanied by the reversal of the momentum p_{x}. Interestingly, the gate voltage can be used to control the reflection and transmission, and therefore tune the spin polarization of the reflected and transmitted electrons.
For simplicity, we introduce the dimensionless units: l_{B }= [ħ/eB_{0}]^{1/2}, B(x) → B_{0}B(x), A_{y}(x) → B_{0}l_{B}A_{y}(x), E → E_{0}E, r → l_{B}r, k → k/l_{B}, and rewrite the Hamiltonian as
In the presence of the magneticelectric potentials, the wavevector satisfies . Note that the translational invariance along the y direction gives rise to conservation of k_{y}, and thus the solutions can be written as . In the free region, the eigenvalues are E_{± }= ±v_{F}π, and the corresponding eigenvectors are
In the barrier region the vector potential is A = (0, Bx, 0) and the solution can be expressed in terms of parabolic cylinder functions D_{ν},
with complex coefficients c_{± }and υ ≡ (E_{F } V)^{2 } (gμ_{B}B)^{2}. After some lengthy algebra, all of the above coefficients of the wave functions can be obtained from the boundary conditions.
3. Spin and momentum filtering
First we consider that the incident electrons are spin polarized along the direction of the vector , i.e., perpendicular to the direction of the electron motion. We examine the transmission characteristics for a pure electric barrier (B = 0, V = 4), which are shown in Figure 2a,b. In Figure 2a we plot the transmission probability as a function of the incident angle and energy. The symmetry of the lowenergy effective Hamiltonian when A = 0 ensures the invariance of the transmission with respect to k_{y }→ k_{y}. Note that, perfect transmission always exists in the vicinity of normal incidence, i.e., this is the socalled Klein tunneling, induced by the helical property of the Dirac fermion. This perfect tunneling process can even occur for a low incident energy, and a high and wide barrier. Figure 2b shows the spin orientation as a function of the incident angles and the incident energies. From the spin orientation one can see that the transmitted electron spins are polarized along the same direction as the incident electron spin, indicating that a pure electric barrier without magnetic field will not affect the spin orientation, the reason is that an electric barrier cannot bend the trajectories of the electrons and the outgoing electrons will propagate in the same direction as the incident ones. The spin of the reflected electron, with reversed longitudinal wavevector and conserved transverse wavevector, is rotated due to the spinmomentum locking.
Figure 2. Transmission spectrum through the single magnetic/electric barrier and spin rotation of transmitted/reflected electrons. (a) The contour plot of the transmission probability as a function of the incident angle and the incident energy, for a fixed barrier width d = 1 and height V = 4, for magnetic field B = 0. (b) The spin orientation of the incident (black arrows), transmitted (red arrows), and reflected (green arrows) electrons as a function of the incident angle and energy. The lengths of the arrows indicate the transmitted or reflected probabilities. The magnetic field unit is B_{0 }= 1 T, the energy unit is E_{0 }= 16 meV, and the length unit is l_{B }= 26 nm.
Next we consider the tunneling process through a pure magnetic single barrier. The transmission probability is shown in Figure 3a for a magnetic barrier (B = 1) with width d = 1. Compared to the pure electric barrier case (see Figure 2a), the transmission becomes asymmetric with respect to the inplane momentum k_{y }parallel along the interface, since the magnetic field breaks the time reversal symmetry. Note that Figure 3a implies that for a certain incident angle ϕ, tunneling is forbidden, the boundary of the total reflection region (T = 0) can be approximately given by the relation k_{y }+ Bd = E, which implies that the transmitted wave is an evanescent mode which decays exponentially along the propagating direction. This total reflection is always present regardless of the incidence angle ϕ when E < Bd/2, and results in an asymmetric behavior of the transmission as a function of the incident angle. This total reflection can be understood semiclassically: when the width d of the barrier is large as compared to the Fermi wavelength λ_{F}, the electron inside the barrier will move on a cyclotron orbit, and if the cyclotron orbit radius R_{c }< d, the incident electron will exit the barrier region backwards eventually. This feature illustrates that the Dirac quasiparticles can be confined by the inhomogeneous magnetic field, and, more interestingly, the incident electron spinoriented parallel along the interface will be flipped in the total reflection region. Figure 3b shows the spin orientation as a function of the incident angles and energies, which is also asymmetric with respect to the incident angles, as a consequence of the perpendicular inhomogeneous magnetic field. Note that the transmitted electron spin no longer points along the same direction of the incident electron spin, but is rotated over an angle, which is determined by the incident angle, the electron energy and the barrier width. In addition, changing the length of the barrier and/or the magnetic field can tune the total reflection and the perfect transmission regions.
Figure 3. The same as Figure 2, but now for a pure magnetic barrier with B = 1, and V = 0.
We also examined the electron transmission through a combined electric and magnetic barrier. The transmission and spin orientation (see Figure 4a,b) become very different from that of the pure magnetic barrier case. The interplay between electric barrier and magnetic field strongly reduces the perfect transmission region, and provides us with an additional way to control the transmission and the spin orientation of the transmitted and reflected electrons.
Figure 4. The same as Figure 2, but now for a combined electric and magnetic barrier with B = 1, and V = 4.
The transmission feature shown above can be examined by the measurable quantities, the conductance G, and Fano factor F [22,23]. The ballistic conductance and Fano factor for a given Fermi energy at zero temperature are given by summing over the modes,
where is taken as the conductance unit. The dependence of the conductance G and the Fano factor F on the incident energy E_{F }are plotted in Figure 5. In lowenergy region, the magnetic barrier can suppress the electron transmission probability for any incident angles, and thus lead to vanished conductances. The conductance increases as the Fermi energy since the magnetic field is not strong enough to confine the electrons. For a pure electric barrier, the conductance is very large in low energy region, since the electrons can effectively tunnel through the barrier via the down branch of the surface state, i.e., the Klein tunneling. When the incident energy approaches the top of an electric barrier, the conductance decrease drastically. This is because the wave vectors in the electric barrier become imaginary if the electrons are not normal incident. The Fano factor shows good correspondence with the changes of conductance. We can see clearly a significant change when the barrier changes from transparent to opaque. These quantities are what experimentalists can actually measure.
Figure 5. Conductance and Fano factor of the single magnetic/electric barrier. (a) The conductance as a function of the incident energy with B = 1, V = 0 (black solid line); B = 1, V = 4 (red solid line); and B = 0, V = 4 (blue solid line). (b) The Fano factor as a function of the incident energy for the same magnetic/electric barrier.
4. Conclusion
In summary, we investigated theoretically quantum tunneling processes through a single electric and magnetic barrier on the surface of a 3D topological insulator. Our theoretical results show that the propagating behavior of electrons in such a 2D system can be controlled by electric and magnetic barriers. We have also observed the rotation of the electron spin, depending on the parameters of the magnetic barrier structure as well as the incident angle and energy. This investigation could be helpful to offer a promising functional unit for future wavevector, and/or spin filtering quantum devices.
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
ZW conceived of the study and carried out the numerical calculation. JL participated in establishing the physical model and developing the numerical code. All authors have participated in the interpretation of the numerical results. All authors read and approved the final manuscript.
Acknowledgements
The study was supported by NSFC Grant No. 11104232.
References

Kane CL, Mele EJ: Z_{2 }Topological order and the quantum spin Hall effect.
Phys Rev Lett 2005, 95:146802.
Kane CL, Mele EJ: Quantum spin Hall effect in graphene. Phys Rev Lett 2005, 95:226801
PubMed Abstract  Publisher Full Text 
Bernevig BA, Zhang SC: Quantum spin Hall effect.
Phys Rev Lett 2006, 96:106802. PubMed Abstract  Publisher Full Text

Bernevig BA, Hughes T, Zhang SC: Quantum spin Hall effect and topological phase transition in HgTe quantum wells.
Science 2006, 314:1757. PubMed Abstract  Publisher Full Text

Moore JE, Balents L: Topological invariants of timereversalinvariant band structures.

Fu L, Kane CL, Mele EJ: Topological insulators in three dimensions.
Phys Rev Lett 2007, 98:106803. PubMed Abstract  Publisher Full Text

Fu L, Kane CL: Topological insulators with inversion symmetry.

Qi XL, Hughes TL, Zhang SC: Topological field theory of timereversal invariant insulators.

Zhang HJ, Liu CX, Qi XL, Dai X, Fang Z, Zhang SC: Topological insulators in Bi_{2}Se_{3}, Bi_{2}Te_{3 }and Sb_{2}Te_{3 }with a single Dirac cone on the surface.
Nat Phys 2009, 5:438. Publisher Full Text

Hsieh D, Qian D, Wray L, Xia Y, Hor YS, Cava RJ, Hasan MZ: A topological Dirac insulator in a quantum spin hall phase.
Nature 2008, 452:970. PubMed Abstract  Publisher Full Text

Hsieh D, Xia Y, Wray L, Qian D, Pal A, Dil JH, Osterwalder J, Meier F, Bihlmayer G, Kane CL, Hor YS, Cava RJ, Hasan MZ: Observation of unconventional quantum spin textures in topological insulators.
Science 2009, 323:919. PubMed Abstract  Publisher Full Text

Xia Y, Qian D, Hsieh D, Wray L, Pal A, Lin H, Bansil A, Grauer D, Hor YS, Cava RJ, Hasan MZ: Observation of a largegap topologicalinsulator class with a single Dirac cone on the surface.
Nat Phys 2009, 5:398. Publisher Full Text

Chen YL, Analytis JG, Chu JH, Liu ZK, Mo SK, Qi XL, Zhang HJ, Lu DH, Dai X, Fang Z, Zhang SC, Fisher IR, Hussain Z, Shen ZX: Experimental realization of a threedimensional topological insulator, Bi_{2}Te_{3}.
Science 2009, 325:178. PubMed Abstract  Publisher Full Text

Castro Neto AH, Guinea F, Peres NMR, Novoselov KS, Geim AK: The electronic properties of graphene.
Rev Mod Phys 2009, 81:109. Publisher Full Text

Hsieh D, Xia Y, Qian D, Wray L, Dil JH, Meier F, Osterwalder J, Patthey L, Checkelsky JG, Ong NP, Fedorov AV, Lin H, Bansil A, Grauer D, Hor YS, Cava RJ, Hasan MZ: A tunable topological insulator in the spin helical Dirac transport regime.
Nature 2009, 460:1101. PubMed Abstract  Publisher Full Text

Zhang LB, Cheng F, Zhai F, Chang K: Electrical switching of the edge channel transport in HgTe quantum wells with an inverted band structure.

Zhu JJ, Yao D, Zhang SC, Chang K: Electrically Controllable Surface Magnetism on the Surface of Topological Insulators.
Phys Rev Lett 2011, 106:097201. PubMed Abstract  Publisher Full Text

Ye PD, Weiss D, Gerhardts RR, Seeger M, von Klitzing K, Eberl K, Nickel H: Electrons in a periodic magnetic field induced by a regular array of micromagnets.
Phys Rev Lett 1995, 74:3013. PubMed Abstract  Publisher Full Text

Matulis A, Peeters FM, Vasilopoulos P: Wavevectordependent tunneling through magnetic barriers.
Phys Rev Lett 1994, 72:1518. PubMed Abstract  Publisher Full Text

Hong J, Joo S, Kim TS, Rhie K, Kim KH, Kim SU, Lee BC, Shinc KH: Local Hall effect in hybrid ferromagnetic/semiconductor devices.
Appl Phys Lett 2007, 90:023510. Publisher Full Text

Cerchez M, Hugger S, Heinzel T, Schulz N: Effect of edge transmission and elastic scattering on the resistance of magnetic barriers: Experiment and theory.
Phys Rev B 2007., 75
035341 and references therein

Ion G, Franz M: Inverse spingalvanic effect in the interface between a topological insulator and a ferromagnet.
Phys Rev Lett 2010, 104:146802. PubMed Abstract  Publisher Full Text

Berry MV, Mondragon RJ: Neutrino billiards: timereversal symmetrybreaking without magnetic fields.
Proc R Soc London, Ser A 1987, 412:53. Publisher Full Text

Tworzyd J, Trauzettel B, Titov M, Rycerz A, Beenakker CWJ: Subpoissonian shot noise in graphene.
Phys Rev Lett 2006, 96:246802. PubMed Abstract  Publisher Full Text