Itinerant Spin Dynamics in Structures of ... - Jacobs University

Itinerant Spin Dynamics in Structures of ... - Jacobs University Itinerant Spin Dynamics in Structures of ... - Jacobs University

jacobs.university.de
from jacobs.university.de More from this publisher
11.03.2014 Views

84 Chapter 5: Spin Hall Effect DOS, their average distance and the disorder strength [SBK + 10]. Another way for efficiently injecting spin currents into semiconductors is to make use of the spin Hall effect (SHE). This effect was first proposed by D’yakonov and Perel [DP71a] and describes in today’s terminology the extrinsic SHE which requires spin dependent impurity scattering. It was experimentally confirmed by the angle-resolved optical detection of spin polarization at the edges of a two-dimensional layer [WKSJ05, KMGA04]. The theory of the SHE has been developed in the last 10 years, as reviewed in Refs. [Sch06, ERH07]. The spin transport not only occurs due to the spin precession in the bulk, but is affected by the scattering from nonmagnetic impurities, which can depend on the direction of the spin itself due to the so called skew scattering and the side jump mechanism[Sch06]. In the first part of the following chapter, we will outline the formalism to calculate analytically the SHE which arises even in the absence of impurities, the so called intrinsic SHE, which is due to the bulk SOC. For the clean case we will include both, Rashba and Dresselhaus SOC. In the second part we focus on numerical methods: The first attempt is application of exact diagonalization which has of course strong limitations concerning the system size. To overcome this limitations, we apply the Kernel Polynomial Method (KPM) in the last part of this thesis. It is first applied to treat the metal-insulator transition (MIT) in a symplectic system, finally we calculate the SHE using the KPM. 5.1.1 About the Definition of Spin Current We have learned in Sec.2.3.4 that in presence of the spin-orbit interaction, the spin current components are not conserved[CSS + 04, SZXN06] even if we assume no spin relaxation: In the continuity equation ∂s z ∂t +D e∇J spin = T s − 1 (ˆτ s ) ij s j (5.1) = τ〈∇v F (B SO (k)×S) i 〉− 1 (ˆτ s ) ij s j . (5.2) anadditionaltorquetermT s appearsincontrasttoEq.2.12besidesthetermwhichdescribes the spin relaxation rate 1/τ s . It follows that Noether’s theorem is not applicable to define the spin current. This is a reason why the comparison between results done with Kubo formalism and calculations using Landauer-Büttiker approach is not trivial, where the spin

Chapter 5: Spin Hall Effect 85 current in the latter is defined by[MM05] where I s p,µ = 1 4π ∑ p≠q,ν Tr[Γ µ pG R Γ ν qG A ](V p −V q ), (5.3) Γ µ p = i(Σ µ p −(Σ µ p) † ), (5.4) with the retarded self energy Σ µ p due to coupling of lead p, which has voltage V p , and sample for spin channel µ and advanced and retarded Green’s function G A/R . However, in an experiment the SHE can be measured. This measurement has to be connected to the theoretical description using linear response to a transverse electrical field with the frequency ω which yields, Appendix B, Eq.(B.16), σ µν (ω) = i V ∑ m,n (f(E m )−f(E n ))〈m|j ν |n〉〈n|j µ |m〉 E n −E m E n −E m +ω +iη . (5.5) To calculate the spin Hall conductivity (SHC) the correlation function consists, in contrast to charge conductivity, of the charge current and a current which contains a spin operator. In the simplified model we use the spin current which is given by the anticommutator of the spin and the group velocity v = i[H,r], J z = 4 {σ z,v}. (5.6) This spin current does not differ from the one defined in Sec.2.3.[ESL05, BNnM04] Note that this quantity is dissipationless, because the spin current is even under the time-reversal operation. 5.2 SHE without Impurities: Exact Calculation We consider theHamiltonian for alattice which provides linear Rashba, Eq.(2.15), and linear and cubic Dresselhaus SOC, Eq.(2.14), as introduced in Sec.2.3.3. The confinement to generate the 2D electron gas is in [001] direction. In order to do calculations numerically, one needs to define a tight binding model on a discrete lattice of finite lattice spacing a. It has the following characterization: • square lattice,

84 Chapter 5: <strong>Sp<strong>in</strong></strong> Hall Effect<br />

DOS, their average distance and the disorder strength [SBK + 10].<br />

Another way for efficiently <strong>in</strong>ject<strong>in</strong>g sp<strong>in</strong> currents <strong>in</strong>to semiconductors is to make<br />

use <strong>of</strong> the sp<strong>in</strong> Hall effect (SHE). This effect was first proposed by D’yakonov and Perel<br />

[DP71a] and describes <strong>in</strong> today’s term<strong>in</strong>ology the extr<strong>in</strong>sic SHE which requires sp<strong>in</strong> dependent<br />

impurity scatter<strong>in</strong>g. It was experimentally confirmed by the angle-resolved optical<br />

detection <strong>of</strong> sp<strong>in</strong> polarization at the edges <strong>of</strong> a two-dimensional layer [WKSJ05, KMGA04].<br />

The theory <strong>of</strong> the SHE has been developed <strong>in</strong> the last 10 years, as reviewed <strong>in</strong> Refs.<br />

[Sch06, ERH07].<br />

The sp<strong>in</strong> transport not only occurs due to the sp<strong>in</strong> precession <strong>in</strong> the bulk, but is affected<br />

by the scatter<strong>in</strong>g from nonmagnetic impurities, which can depend on the direction <strong>of</strong> the<br />

sp<strong>in</strong> itself due to the so called skew scatter<strong>in</strong>g and the side jump mechanism[Sch06]. In the<br />

first part <strong>of</strong> the follow<strong>in</strong>g chapter, we will outl<strong>in</strong>e the formalism to calculate analytically<br />

the SHE which arises even <strong>in</strong> the absence <strong>of</strong> impurities, the so called <strong>in</strong>tr<strong>in</strong>sic SHE, which<br />

is due to the bulk SOC. For the clean case we will <strong>in</strong>clude both, Rashba and Dresselhaus<br />

SOC. In the second part we focus on numerical methods: The first attempt is application <strong>of</strong><br />

exact diagonalization which has <strong>of</strong> course strong limitations concern<strong>in</strong>g the system size. To<br />

overcome this limitations, we apply the Kernel Polynomial Method (KPM) <strong>in</strong> the last part<br />

<strong>of</strong> this thesis. It is first applied to treat the metal-<strong>in</strong>sulator transition (MIT) <strong>in</strong> a symplectic<br />

system, f<strong>in</strong>ally we calculate the SHE us<strong>in</strong>g the KPM.<br />

5.1.1 About the Def<strong>in</strong>ition <strong>of</strong> <strong>Sp<strong>in</strong></strong> Current<br />

We have learned <strong>in</strong> Sec.2.3.4 that <strong>in</strong> presence <strong>of</strong> the sp<strong>in</strong>-orbit <strong>in</strong>teraction, the<br />

sp<strong>in</strong> current components are not conserved[CSS + 04, SZXN06] even if we assume no sp<strong>in</strong><br />

relaxation: In the cont<strong>in</strong>uity equation<br />

∂s z<br />

∂t +D e∇J sp<strong>in</strong> = T s − 1<br />

(ˆτ s ) ij<br />

s j (5.1)<br />

= τ〈∇v F (B SO (k)×S) i<br />

〉− 1<br />

(ˆτ s ) ij<br />

s j . (5.2)<br />

anadditionaltorquetermT s appears<strong>in</strong>contrasttoEq.2.12besidesthetermwhichdescribes<br />

the sp<strong>in</strong> relaxation rate 1/τ s . It follows that Noether’s theorem is not applicable to def<strong>in</strong>e<br />

the sp<strong>in</strong> current. This is a reason why the comparison between results done with Kubo<br />

formalism and calculations us<strong>in</strong>g Landauer-Büttiker approach is not trivial, where the sp<strong>in</strong>

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!