arxiv: v1 [cond-mat.mes-hall] 5 Nov 2012

Transcription

1 Spin-orital Texture in Topological Insulators Haijun Zhang 1, Chao-Xing Liu & Shou-Cheng Zhang 1 1 Department of Physics, McCullough Building, Stanford University, Stanford, CA Department of Physics, The Pennsylvania State University, University Park, Pennsylvania (Dated: Novemer 6, 01) arxiv: v1 [cond-mat.mes-hall] 5 Nov 01 Relativistic spin-orit coupling plays an essential role in the field of topological insulators and quantum spintronics. It gives rise to the topological non-trivial and structure and enales electric manipulation of the spin degree of freedom. Because of the spin-orit coupling, rich spin-orital coupled textures can exist oth in momentum and in real space. For three dimensional topological insulators in the Bi Se 3 family, topological surface states with p z oritals have a left-handed spin texture for the upper Dirac cone and a right-handed spin texture for the lower Dirac cone. In this work, we predict a new form of the spin-orital texture associated with the p x and p y oritals. For the upper Dirac cone, a left-handed (right-handed) spin texture is coupled to the radial ( tangential ) orital texture, whereas for the lower Dirac cone, the coupling of spin and orital textures is the exact opposite. The tangential ( radial ) orital texture is dominant for the upper (lower) Dirac cone, leading to the right-handed spin texture for the in-plane oritals of oth the upper and lower Dirac cones. A spinresovled and photon polarized angle-resolved photoemission spectroscopy experiment is proposed to oserve this novel spin-orital texture. PACS numers: , f,73.0.-r I. INTRODUCTION Three-dimensional topological insulators(tis) are new states of quantum matter with helical gapless surface states consisting of odd numer of Dirac cones inside the ulk and gap protected y time-reversal symmetry (TRS). 1 4 The underlying physical origin of the topological property of TIs is the strong spin-orit coupling (SOC), which plays a similar role as the Lorentz force in the Quantum Hall state. Due to the SOC interaction, the spin and momentum are locked to each other, forming a spin texture in the momentum space for the surface states of TIs 5 7. The spin texture has een directly oserved in the spin-resolved angle-resolved photon emission spectroscopy (spin-resolved ARPES) 8 1. The spin texture gives rise to a non-trivial Berry phase for the topological surface states and suppresses the ackscatterings under TRS, leading to possile device applications in spintronics. Besides the spin texture, it has also een shown recently that the atomic p oritals of the Bi Se 3 family of topological insulators form a pattern in the momentum space, dued as the orital texture, for the topological surface states. 13,14 In this work, we predict a coupled spin-orital texture for the topological surface states. Based on oth the effective k p theory and a-initio calculations, we find, esides the usual locking etween the electron spin and the crystal momentum, the spin texture is also locked to the atomic orital texture, which is dued as spin-orital texture. We show that p z oritals have left-handed spin texture for the upper Dirac cone and right-handed spin texture for the lower Dirac cone, sharing the same feature as the total spin texture of the surface states. In contrast, the in-plane oritals a k y k x FIG. 1. (color online) a,, The tangential orital texture with the right-handed helical spin texture (a) and the radial orital texture with the left-handed helical spin texture () for the upper Dirac cone. (p x and p y oritals) reveal more intriguing features: for the upper Dirac cone of surface states, a radial orital texture is coupled to a left-handed spin texture and a tangential orital texture is coupled to a right-handed spin texture. For the lower Dirac cone, the coupling etween spin and orital textures is exactly opposite. An electron spin-resolved and photon polarized ARPES experiment is proposed to oserve this novel spin-orital texture of the surface states of TIs. II. EFFECTIVE THEORY OF THE SPIN-ORBITAL TEXTURE The surface states of TIs are descried y the Dirac type of effective Hamiltonian 5,15 H surf (k x,k y ) = v f (σ x k y σ y k x ), (1) with the Fermi velocity v f and Pauli matrix σ. The salient feature of this effective Hamiltonian is the spin- k y k x

2 a c FIG.. (color online) a, The Dirac cone of Bi Se 3 on the surface with the normal direction [0001] with the spin texture marked y lue arrows., c, The projection of p z orital and the related in-plane spin texture for upper () and lower (c) Dirac cones. More red means more p z character. The red arrows represent the in-plane spin texture related to the p z oritals. The insets are the schematics of the spin texture marked y green arrows. momentum locking, which means for a fixed momentum k, the spin, denoted y the Pauli matrix σ, has a fixed direction for the eigenstate of the Hamiltonian. Since the spin is always perpendicular to the momentum, we can introduce a helicity operator, defined as ĥ = kẑ ( 1 k σ) which commutates with the Hamiltonian, to determine the handness of the spin texture. For the upper Dirac cone of surface states, the helical operator ĥ = 1, leading to a left-handed spin texture in the momentum space while for the lower Dirac cone, ĥ = 1 yields a a right-handed spin texture. However, one should note that here the spin is not the real spin, ut the total angular momentum J = S + L, which is a comination of the real spin S and the orital angular momentum L due to SOC. Consequently, the asis of the surface effective Hamiltonian (1) are denoted as Ψ Jz=± 1 with the lower indices ± 1 representing the total angular momentum along z direction. In order to understand what is the texture for the real spin S, it is necessary to write down the explicit form of the asis wavefunction Ψ ± 1. The form of the asis Ψ ± 1 can e constructed y symmetry considerations. Generally the asis Ψ ± 1 depends on the momentum k and we can expand it up to the first order in k as Ψ ± 1 = Ψ(0) + Ψ (1) Here we ± 1 ±. 1 areonly interested in the p oritalsof Bi and Se atoms in the topological insulator Bi Se 3 and can decompose the zeroth-order wavefunction as Ψ (0) = ± 1 α [u 0,α α,p z, ( ) +v 0,α α,p ±, ( ) ] () and the first-order wavefunction as Ψ (1) = ± 1 α [±k ± (iu 1,α α,p, ( ) +iv 1,α α,p z, ( ) ) (3) iw 1,α k α,p ±, ( ) ] where k ± = k x ±ik y, and denote the spin, p z and p ± = 1 ( p x ±i p y ) denote different p oritals, and α denotes indices other than the spin and orital, such as atom indices. Here u 0(1),α, v 0(1),α and w 1,α are material-dependent parameters. By comparing with the a-initio calculations, we find that we can take them to e real. Ψ 1 and Ψ 1 are related to each other y TRS. The expressions of the asis () and (3) can e

3 a 3 c d FIG. 3. (color online) a,, c, d, The p x projection on the states of upper (a) and lower (c) Dirac cones, and the p y projection of upper () and lower (d) Dirac cones. More red means more p x character in (a) and (c), and more red means more p y character in () and (d). The red arrows indicate the in-plane spin texture related to the oritals. The insets are the schematics of the spin texture. sustituted into the eigen wavefunctions of the ] Hamiltonian (1), Φ ± = 1 [±ie iθ k Ψ 1 + Ψ 1, yielding the following forms of the wavefunctions Φ + = α Φ = α [(u 0,α v 1,α k) α,p z, θ i (v 0,α u 1,α k w 1,α k) α,p r, θ (4) + 1 (v 0,α u 1,α k +w 1,α k) α,p t, θ ] [(u 0,α +v 1,α k) α,p z, θ + i (v 0,α +u 1,α k +w 1,α k) α,p r, θ 1 (v 0,α +u 1,α k w 1,α k) α,p t, θ ]. (5) Here θ ( θ ) = 1 (+( )ie iθ k + ) stands for the left-handed(right-handed)helicalspintextureand p r = cosθ k p x +sinθ k p y, p t = sinθ k p x +cosθ k p y are the radial and tangential orital textures, as shown in Fig. 1a,, respectively. From the expressions (4) and (5), we can clearly see that p z orital is coupled to the lefthanded spin texture θ for the upper Dirac cone and the right-handed spin texture θ for the lower Dirac cone. Furthermore, for the upper Dirac cone, the radial orital texture p r is always coupled to the left-handed spin texture θ and the tangential orital texture p t is always coupled to the right-handed spin texture θ. The situation is exactly opposite for the lower Dirac cone. The expressions (4) and (5) are the main analytical results of this paper, which show explicitly the spinorital texture. To confirm our analytical results, ainitio method is adopted to calculate the projection of surface states on the spin and orital asis, defined y the quantity D ± i,η = Φ ± ( p i p i s η ) Φ ±, (6) where p i = p x,p y,p z for the three p oritals, s 0 = 1 denotes the charge part and s x,y,z denote the three Pauli matrices for the spin. In the following, we will compare the analytical calculation of the intensity D ± i,η with ainitio calculations.

4 a 4 c Character Density p x upper Dirac cone p y upper Dirac cone p x lower Dirac cone p y lower Dirac cone 0 1/π π 3/π π θ d Polarization Energy (ev) FIG. 4. (color online) a,, The tangential orital texture with the related in-plane spin texture for the upper Dirac cone (a) and the radial orital texture for the lower () from a-initio calculations. More red means more p x character, and more lue means more p y character. The red arrows represent the total in-plane spin texture related to p x and p y oritals. c, The p x and p y projections onto upper and lower Dirac cones. The solid curves are the p x and p y projections for the upper Dirac cone at energy level 0.10eV, and the dashed for the lower Dirac cone at the energy 0.07eV. The asic feature is the π period which exactly agrees with the prediction y the effective model. d, Orital polarization P px. The energy of Dirac point is shifted to e zero. The positive value of the orital polarization represents the radial orital texture, and the negative value represents the tangential orital texture. In order to plot more visually, the P px for the upper Dirac cone is inversed, marked y the red color. P px is exact zero at Dirac point, which indicates the transition point etween the tangential and radial orital textures. III. SPIN-ORBITAL TEXTURE FROM AB-INITIO CALCULATIONS The Vienna A-initio Simulation Package (VASP) 16,17 is employed to carry out a-initio calculations with the framework of the Perdew-Burke-Ernzerhof-type(PBE) 18 generalized gradient approximation (GGA) of density functional theory 19. Projector augmented wave (PAW) pseudo-potentials are used for all of the calculations in this work and are used for k-grid of ulk and free-standing calculations, respectively. The kinetic energy cutoff is fixed to 450eV. 6 quintuple layers (QLs) are fixed in the supercell for free-standing calculations, and the thickness of vacuum is taken to e 50Å. The lattice constant and the atomic parameters are directly otained from experiments. SOC is included with the non self-consistent calculation. In order to compare with the result of ARPES experiments, the projections of all the oritals are only for the first Se and Bi atoms on the top surface of the free-standing model. The surface states of Bi Se 3 consist of a single Dirac cone at Γ point on one surface inside the large ulk and gap ( 0.3eV) 5,1, which provides an ideal material to study the coupling of spin and orital textures of surface states. As the starting point, we compare the ulk and structure of Bi Se 3 with the previous calculation 5 and find good agreements. The surface states are otained from the calculation of a free-standing structure with the normal direction [0001], as shown in Fig. a. The lue arrows represent the spin texture, where the spin is mainly lying in plane near the Dirac point. The spin texture is left-handed for the upper Dirac cone and right-handed for the lower one, the same as the total angular momentum texture. To understand the underlying physics, we calculate the spin texture for different atomic oritals. For p z oritals, a left-handed helical

5 5 spin texture is found for the upper Dirac cone and a right-handed texture for the lower Dirac cone, as shown in Fig.,c. The schematic of the spin texture are shown in the inset. Here the ackground color indicates the projection of p z oritals, which is isotropic, and the red arrows represent the corresponding in-plane spin texture. The spin texture of p z oritals can e reproduced with the expressions [D ± p z,σ x,d ± p z,σ y,d ± p z,σ z ] = ± α (u 0,α v 1,α k) [sinθ k, cosθ k,0] with ± for the upper and lower Dirac cone and k = (k,θ k ) in the polar coordinate. The spin textures for in-plane oritals are shown in Fig. 3a, for the upper Dirac cone and in Fig. 3c,d for the lower Dirac cone, respectively. We find that the associated spins for p x and p y oritals don t rotate clockwise or anti-clockwise around the Dirac point as in the case of p z oritals, ut instead, they take the form [D ± p x,x,d± p x,y,d± p x,z ] = α [D ± p y,x,d± p y,y,d± p y,z ] = ± α v 0,α [sinθ k,cosθ k,0] (7) v 0,α [sinθ k,cosθ k,0] (8) for small k around the Γ point. The corresponding spin textures are shown schematically in the inset of Fig. 3a,c for p x oritals in upper and lower Dirac cones and in the inset of Fig. 3,d for p y oritals. Unlike p z oritals, the amplitude of p x and p y oritals for the surface states is not isotropic, ut has θ k angular dependence around the Fermi surface, as shown in Fig. 4c. We may take the difference of the amplitudes etween p x and p y oritals, as shown y colors in Fig. 4a,. Here more red means more p x character, and more lue means more p y character. The angular dependence indicates a tangential orital texture for the upper Dirac cone and a radial orital texture for the lower Dirac cone, as schematically shown y the inset of Fig. 4a,, respectively. This orital texture was experimentally oserved recently 14. Furthermore, we also plot the total spin textures for in-plane oritals on the same figure, which show a right-handed texture for oth upper and lower Dirac cones. All these salient feature can e understood y the wavefunctions (4) and (5). For the upper Dirac cone, although oth θ p r and θ p t terms exist in the wavefunction(4), their associated coefficients are unequal. When α (v 0,α u 1,α k+w 1,α k) > α (v 0,α u 1,α k w 1,α k), θ p t term dominates over θ p r term, dominantlygivingatangentialoritaltexture coupled to a right-handed spin texture. Similarly, for the lower Dirac cone, when α (v 0,α+u 1,α k+w 1,α k) > α (v 0,α + u 1,α k w 1,α k), θ p r term in the wavefunction (5) is dominant, yielding a radial orital texture coupled to a right-handed spin texture. The difference etween p x and p y oritals can e calculated directly as D p ± x,0 D± p = cosθ y,0 k α [(v 0,α ku 1,α )kw 1,α ], which indeed shows a θ k angular dependence, and the total spin textures for in-plane oritals can e otained as [D x ±,D± y ] = [D± p + x,x D± p y,x,d± p + x,y D± p ] = y,y 4[ sinθ k,cosθ k ] α (v 0,α ku 1,α )kw 1,α, which shows a right-handedspintexture when α (v 0,α ku 1,α )kw 1,α > 0. Especially, if k gets close zero, oth the total spin texture [D x ±,D± y ] of the in-plane oritals and the difference etweenp x andp y oritalsd p ± x,0 D± p approacheszero, y,0 also as shown in Fig. 4a,. Therefore, there is a transition from a tangential orital texture in the upper Dirac cone to a radial orital texture in the lower Dirac cone, switching exactly at the Dirac point. To quantitatively descrie this transition, we introduce a polarization quantity P px (±) = D p x,0(±,θ = 0) D px,0(±,θ = 90) D px,0(±,θ = 0)+D px,0(±,θ = 90) (9) with ± for upper and lower Dirac cones. The plot of P px (±) is shown in Fig. 4d where the energy level of the Dirac point is shifted to zero. In order to show the plot more visually, we reverse the value of the P px (+) for the upper Dirac cone plotted with the red. The feature of P px (±) undoutedly indicates that the state of the lower Dirac cone forms a radialorital texture, and the state of the upper Dirac cone forms a tangential orital texture. The Dirac point is shown to e the exact transition point from the tangential to radial orital texture. This is exactly the ehavior oserved in a recent experiment and explained within the first principle calculations 14. The numerical results fit well to the analytical calculation, with the expression P px (±) = α α (v 0,α Eu 1,α / v f )Ew 1,α / v [ f (v 0,α Eu 1,α / v f ) +E w1,α / vf (10) with the energy E. For small E around Dirac point, P px (±) α v0,αw1,α E α (v0,α) v f shows the linear dependence on energy, as found in Fig. 4d. Although in-plane oritals show different spin textures compared to p z oritals, we stress the p z oritals (50%) dominate the states near the Dirac point with p x and p y only around 30%. Therefore, the spin texture for the whole states show left-handed for the upper Dirac cone and right-handed for the lower Dirac cone, the same as thatofp z oritals,aswellasthetotalangularmomentum texture. IV. DISCUSSION In order to detect the spin texture of electrons, the spin-resolved ARPES technology has een developed y taking advantage of spin-dependent scatting processes and precisely measuring the magnitude of the asymmetry in the spin-dependent intensity with perfect spin-polarimeters. 1 The non-trivial spin texture of surface states of TIs has een clearly oserved y experiments. 8 1 In addition, the orital character can e detected through the photon polarization selection ]

6 6 rules ased on the symmetry analysis. With this technology, the orital texture of surface states of Bi Se 3 was reported recently y a polarized ARPES experiment. 14 Therefore, it is possile to comine these two technologies together in an electron spin-resolved and photon polarized ARPES experiment, with oth the spin and orital textures extracted in the same measurement. The predicted spin-orital texture can e directly confirmed in this type of experiment, which can explicitly reveal how SOC plays a role in the real material at the atomic level. We would like to thank Dr. Dan Dessau for sharinghis experimental data and for useful discussions, which partially motivated this work. This work is supported y the Defense Advanced Research Projects Agency Microsystems Technology Office, MesoDynamic Architecture Program (MESO) through the contract numer N and y the DARPA Program on Topological Insulators Solid State Chemistry, New Materials and Properties, under the award numer N V. ACKNOWLEDGMENTS 1 X.-L. Qi and S.-C. Zhang, Physics Today 63, 33 (010). J. E. Moore, Nature 464, 194 (010). 3 M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 8, 3045 (010). 4 X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 (011). 5 H. Zhang, C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. Zhang, Nature Physics 5, 438 (009), /nphys W. Zhang, R. Yu, H.-J. Zhang, X. Dai, and Z. Fang, New Journal of Physics 1, (010). 7 O. V. Yazyev, J. E. Moore, and S. G. Louie, Phys. Rev. Lett. 105, (010). 8 D. Hsieh, Y. Xia, L. Wray, D. Qian, A. Pal, J. H. Dil, J. Osterwalder, F. Meier, G. Bihlmayer, C. L. Kane, Y. S. Hor, R.J.Cava, andm.z.hasan,science33,919(009), 9 S. Souma, K. Kosaka, T. Sato, M. Komatsu, A. Takayama, T. Takahashi, M. Kriener, K. Segawa, and Y. Ando, Phys. Rev. Lett. 106, (011). 10 S.-Y. Xu, L. A. Wray, Y. Xia, F. v. Rohr, Y. S. Hor, J. H. Dil, F. Meier, B. Slomski, J. Osterwalder, M. Neupane, H. Lin, A. Bansil, A. Fedorov, R. J. Cava, and M. Z. Hasan, arxiv: cond-mat/ (011). 11 Z.-H. Pan, E. Vescovo, A. V. Fedorov, D. Gardner, Y. S. Lee, S. Chu, G. D. Gu, and T. Valla, Phys. Rev. Lett. 106, (011). 1 C. Jozwiak, Y. L. Chen, A. V. Fedorov, J. G. Analytis, C. R. Rotundu, A. K. Schmid, J. D. Denlinger, Y.-D. Chuang, D.-H. Lee, I. R. Fisher, R. J. Birgeneau, Z.-X. Shen, Z. Hussain, and A. Lanzara, Phys. Rev. B 84, (011). 13 S. R. Park, J. Han, C. Kim, Y. Y. Koh, C. Kim, H. Lee, H. J. Choi, J. H. Han, K. D. Lee, N. J. Hur, M. Arita, K. Shimada, H. Namatame, and M. Taniguchi, Phys. Rev. Lett. 108, (01). 14 Y. Cao, J. A. Waugh, X.-W. Zhang, J.-W. Luo, Q. Wang, T. J. Reer, S. K. Mo, Z. Xu, A. Yang, J. Schneeloch, G. Gu, M. Brahlek, N. Bansal, S. Oh, A. Zunger, and D.S. Dessau, ArXiv e-prints (01), arxiv: [condmat.mtrl-sci]. 15 C.-X. Liu, X.-L. Qi, H. Zhang, X. Dai, Z. Fang, and S.-C. Zhang, Phys. Rev. B 8, 0451 (010). 16 G. Kresse and J. Hafner, Phys. Rev. B 47, 558 (1993). 17 G. Kresse and D. Jouert, Phys. Rev. B 59, 1758 (1999). 18 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996). 19 P. Hohenerg and W. Kohn, Phys. Rev. 136, B864 (1964). 0 P. E. Blöchl, Phys. Rev. B 50, (1994). 1 Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H.Lin, A.Bansil, D. Grauer, Y. S. Hor, R. J. Cava, and M. Z. Hasan, Nature Physics 5, 398 (009), /nphys174. A. Damascelli, Z. Hussain, and Z.-X. Shen, Rev. Mod. Phys. 75, 473 (003).