11.03.2014 Views

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

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

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

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

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

0.10<br />

V⩵2.8<br />

Σ SH e 8Π<br />

0.05<br />

0.00<br />

0.05<br />

0.10<br />

0.15<br />

0.20<br />

Ρ<br />

0.15<br />

0.10<br />

0.05<br />

2⋆Σ SH e 8Π<br />

4 2 0 2 4<br />

(a)<br />

E<br />

0.00<br />

10 5 0 5 10<br />

(b)<br />

E<br />

Figure 5.4: (a) SHC, as a function <strong>of</strong> Fermi energy <strong>in</strong> units <strong>of</strong> t, <strong>in</strong> presence <strong>of</strong> impurities<br />

<strong>of</strong> b<strong>in</strong>ary type calculated us<strong>in</strong>g exact diagonalization. The impurity strength is V = −2.8t<br />

with a concentration <strong>of</strong> 10%. The system size is L 2 = 32 2 , and the SOC is Rashba type<br />

with α 2 = 1.2t with cut<strong>of</strong>f η = 0.06. (b) Comparison <strong>of</strong> a) with DOS (blue curve).<br />

KPM <strong>in</strong> a Nutshell<br />

The Kernel Polynomial Method (KPM) was first proposed by Silver et al.[SR94]<br />

to calculate DOS <strong>of</strong> large systems. It is a method to expand <strong>in</strong>tegrable functions def<strong>in</strong>ed<br />

on a f<strong>in</strong>ite <strong>in</strong>terval f : [a,b] −→ R <strong>in</strong> terms <strong>of</strong> Chebyshev polynomials <strong>of</strong> the first,<br />

T n (x) = cos(narccos(x)), (5.31)<br />

or second k<strong>in</strong>d<br />

i.e., we can write for <strong>in</strong>stance<br />

f(x) =<br />

U n (x) = s<strong>in</strong>((n+1)arccos(x)) ,<br />

s<strong>in</strong>(arccos(x))<br />

[<br />

1<br />

π √ µ 0 +2<br />

1−x 2<br />

]<br />

∞∑<br />

µ n T n (x) , (5.32)<br />

if we assume that the function f has been rescaled to ˜f : [−1,1] −→ R and can be expanded<br />

us<strong>in</strong>g the polynomials <strong>of</strong> the first k<strong>in</strong>d, which are (as the one <strong>of</strong> the second k<strong>in</strong>d) def<strong>in</strong>ed<br />

n=1<br />

on the <strong>in</strong>terval [−1,1]. If so, the coefficients are given by<br />

µ n =<br />

∫ 1<br />

−1<br />

dx ˜f(x)T n (x). (5.33)

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

Saved successfully!

Ooh no, something went wrong!