Trieste MaX Conference, 31. Jan. 2018
Peter Grünberg Institute and Institute for Advanced Simulation
Stefan Blügel
Spin-Orbit Interaction – A Path to Topological Matter in Real and Momentum Space
PAGE 1
Spin-Orbit Interaction A Path to Topological Matter in Real and - - PowerPoint PPT Presentation
Spin-Orbit Interaction A Path to Topological Matter in Real and Momentum Space Peter Grnberg Institute and Institute for Advanced Simulation Stefan Blgel PAGE 1 Trieste MaX Conference, 31. Jan. 2018 Topology of electrons in an
Trieste MaX Conference, 31. Jan. 2018
Peter Grünberg Institute and Institute for Advanced Simulation
Stefan Blügel
Spin-Orbit Interaction – A Path to Topological Matter in Real and Momentum Space
PAGE 1
Trieste MaX Conference, 31. Jan. 2018
Fibre bundle theory
integer Chern number – topological invariant
Nash and Sen, Topology and Geometry for Physicists
H(k)
PAGE 2
Trieste MaX Conference, 31. Jan. 2018
Topological matter
momentum energy edge states
valence band conduction band
2D
valence band conduction band
Topology of Bloch wavefunction
PAGE 3
Topological classification
Z2 = 1 Z2 = 0
Trieste MaX Conference, 31. Jan. 2018
PAGE 4
Dissipationless edge states Topological matter
momentum energy edge states
valence band conduction band
2D
Quantum Spin Hall Effect defect
Trieste MaX Conference, 31. Jan. 2018
, j E, j ˆ m ˆ m z z T odd T
even
(b)
Spin Orbit Torque
3D Topological Insulators Response Properties
Aguilera et al., PRB 88, 045206 (2013)
Energy Relativistic GW Bi2Se3
electronic structure causes non- trivial topological invariants Berry curvature
Z ∼ ZZ Ω ˆ
mkxdkxd ˆ
m
Bi2Te3 Bi2Se3 Sb2Te3
QAnomalous HE
M
Generation of Spin-Currents
Goal: § Exploration of topological phase space
PAGE 5
Trieste MaX Conference, 31. Jan. 2018
Im Ma Juba Da Sa
from Bertrand Dupé
PAGE 6
Trieste MaX Conference, 31. Jan. 2018
Im Ma Juba Da Sa
from Bertrand Dupé from Karin Everschor-Sitte
PAGE 7
Trieste MaX Conference, 31. Jan. 2018
PAGE 8
Skyrmion= non-trivial, smooth mapping from Sd to order parameter space (“trivial winding at infinity”) magnetization direction Smooth mapping Here d=2, S2 → S2 hedgehog vector field of magnetization direction m(x,y)= M/M
Q = 1 4π Z
R2 m ·
✓∂m ∂x × ∂m ∂y ◆ dxdy
m(x,y)
Trieste MaX Conference, 31. Jan. 2018
Ultrathin films with induced chirality (Fe/Ir, Mn/W, Pd/Fe/Ir) Layers of materials with intrinsic chirality (cubic helimagnets FeGe, MnSi, Fe1-xCoxSi)
Lorentz Transmission Electron Microscopy
X.Z. Yu et al. Nature 465, 90 (2010)
Magnetic Force Microscopy
Spin-Polarized Scanning Tunneling Microscopy
Bapp ≠ 0
PAGE 9
Trieste MaX Conference, 31. Jan. 2018
Multiscale modeling v Spin-Lattice Model: v Micromagnetic-model:
H = 1 2 X
ij
Jij mimj+ X
ij
Dij
c
z }| { mi × mj + X
i
mi K mi+ X
ij
1 r 3
ij
⇥ mimj − (mi ˆ ei)(mj ˆ ei) ⇤
v DFT-model: E DFT
tot (q, ˆ
erot) = E DFT
noSOC(q) + ∆E DFT SOC (q, ˆ
erot)
E(m) = Z
R2
⇥ A |rm|2 + D : (rm ⇥ m) + m · K · m B m · ˆ ez ⇤ dr
From total energy calculation to
PAGE 10
Trieste MaX Conference, 31. Jan. 2018
Multiscale modeling v Spin-Lattice Model: v Micromagnetic-model:
H = 1 2 X
ij
Jij mimj+ X
ij
Dij
c
z }| { mi × mj + X
i
mi K mi+ X
ij
1 r 3
ij
⇥ mimj − (mi ˆ ei)(mj ˆ ei) ⇤
v DFT-model: E DFT
tot (q, ˆ
erot) = E DFT
noSOC(q) + ∆E DFT SOC (q, ˆ
erot)
q ˆ erot ˆ erot
E(m) = Z
R2
⇥ A |rm|2 + D : (rm ⇥ m) + m · K · m B m · ˆ ez ⇤ dr
PAGE 11
Trieste MaX Conference, 31. Jan. 2018
Ab-initio A, D, K v Spin-Lattice Model: v Micromagnetic-model:
H = 1 2 X
ij
Jij mimj+ X
ij
Dij
c
z }| { mi × mj + X
i
mi K mi+ X
ij
1 r 3
ij
⇥ mimj − (mi ˆ ei)(mj ˆ ei) ⇤
v DFT-model: E DFT
tot (q, ˆ
erot) = E DFT
noSOC(q) + ∆E DFT SOC (q, ˆ
erot)
§ Spin Stiffness: § Spiralization (micromagnetic D)
A = ∂2 ∂q2 EDFT
tot (q) ∝
X
j>0
J0jR2
0j
D = ∂ ∂qEDFT
tot (q) ∝
X
j>0
D0j ⊗ R0j
E(m) = Z
R2
⇥ A |rm|2 + D : (rm ⇥ m) + m · K · m B m · ˆ ez ⇤ dr
PAGE 12
Trieste MaX Conference, 31. Jan. 2018
Schematic representation of workflow in KKRnano
Calculate charge density Calculate new potential Mixing potentials i f Setting up reference system Preconditioned iterative solution of sparse linear equation (Dyson-equation)
PAGE 13
Trieste MaX Conference, 31. Jan. 2018
(GaAs, InSb, interfaces, surfaces, ...) Time reversal + space inversion symmetry: Time reversal only , Effective spin-orbit (“magnetic”) field Ω: Time reversal symmetry:
PAGE 14
Trieste MaX Conference, 31. Jan. 2018
v spin-orbit coupling has fascinating realizations and ramifications in solids Examples:
PAGE 15
Trieste MaX Conference, 31. Jan. 2018
Energy Storage Memory
permanent magnets magneto-caloric materials hard disk drive MRAM TMR
IoT
magnetic sensors
PAGE 16
Trieste MaX Conference, 31. Jan. 2018
PAGE 17
Trieste MaX Conference, 31. Jan. 2018
18
MOST GW WORKS PUBLISHED a posteriori SOC:
LDA (without SOC) + GW (without SOC) + SOC(LDA)
GW+SOC
Trieste MaX Conference, 31. Jan. 2018
19
OUR WORK full SOC:
LDA (with SOC) + GW (with SOC)
(more accurate but ~10 times more time-consuming)
MOST GW WORKS PUBLISHED a posteriori SOC:
LDA (without SOC) + GW (without SOC) + SOC(LDA)
Sakuma et al., PRB 84 085144 (2011)
GSOCWSOC GW+SOC
Trieste MaX Conference, 31. Jan. 2018
20
OUR WORK full SOC:
LDA (with SOC) + GW (with SOC)
(more accurate but ~10 times more time-consuming)
Sakuma et al., PRB 84 085144 (2011)
GW+SOC GSOCWSOC
Aguilera, Friedrich, Blügel, PRB 88, 165136 (2013)
LDA+SOC
GSOCWSOC
GW+SOC
MOST GW WORKS PUBLISHED a posteriori SOC:
LDA (without SOC) + GW (without SOC) + SOC(LDA)
Trieste MaX Conference, 31. Jan. 2018
min
"GW "
max
contribution
"LDA" (~100 nm)
21
Trieste MaX Conference, 31. Jan. 2018
max
contribution
min
Dispersion of the lower Dirac cone?
"LDA" "GW " (~100 nm)
22
Comparison with ARPES: Bi2Se3
23
Comparison with ARPES: Bi2Se3
ARPES 0.3 Å-1 0.3 Å-1 1 eV "GW " "LDA" 0.3 Å-1 1 eV
24
Trieste MaX Conference, 31. Jan. 2018
PAGE 25
Trieste MaX Conference, 31. Jan. 2018
manipulation at reasonable fields and currents
Albert Fert, Vincent Cross and João Sampaio, Nature Nanotechnology 8, 152 (2013)
The Fert criteria
PAGE 26
Trieste MaX Conference, 31. Jan. 2018
Multiscale modeling v Spin-Lattice Model: v Micromagnetic-model:
H = 1 2 X
ij
Jij mimj+ X
ij
Dij
c
z }| { mi × mj + X
i
mi K mi+ X
ij
1 r 3
ij
⇥ mimj − (mi ˆ ei)(mj ˆ ei) ⇤
E(m) = Z
R2
⇥ A |rm|2 + D : (rm ⇥ m) + m · K · m B m · ˆ ez ⇤ dr
PAGE 27
Trieste MaX Conference, 31. Jan. 2018
Spin-polarized STM image Theory result
20 15 10 5
0.0 0.4 0.8
SR SOC
2.5 (nm)
4 2
E
T O T
) m
a / V e m (
0.0 0.2
1.3
K110 K100
Ferriani et al., PRL 101 027201 (2008)
Mn/W(100)
PAGE 28
H = 1 2 X
ij
Jij mimj+ X
ij
Dij
c
z }| { mi × mj + X
i
mi K mi+ X
ij
1 r 3
ij
⇥ mimj − (mi ˆ ei)(mj ˆ ei) ⇤
Trieste MaX Conference, 31. Jan. 2018
SS phase SkL phase iSk phase
<010> <100> Spontaneous nucleation of individual skyrmion with finite life-time
Nandy, Kiselev, Blügel, PRL.116, 177202 (2016)
PAGE 29
Trieste MaX Conference, 31. Jan. 2018
Interlayer exchange coupling (IEC) between reference and free magnetic layer may compensate the required magnetic field.
Nandy, Kiselev Blügel PRL.116, 177202 (2016)
PAGE 30
Trieste MaX Conference, 31. Jan. 2018
Mitglied der Helmholtz-Gemeinschaft
SS phase SkL phase iSk phase
KCo = 2.0 meV/Co, Beff = 20 T Mn/W7/Co4/Pt/W(001) KCo = 2.7 meV/Co,
Beff = 15 T
PAGE 31
Trieste MaX Conference, 31. Jan. 2018
Unrelaxed Relaxed
Steep slope at the Fermi energy System is extremely sensitive on lattice relaxations Energy shifts due to Hybridization effects
J1,5(E) = −1 3Trα,β 1 π Im Z EF dE Tr ⇥ G1,5 tα
5 G5,1 tβ 1
⇤
j(E) = dJ dE
Unrelaxed Relaxed
Bauer, Mavropoulos, Zeller, Blügel to be published
Fe/Ir(111)
PAGE 32
Trieste MaX Conference, 31. Jan. 2018
Ir(111) Fe Pd
Nature Comm. 7, 11779 (2016)
PAGE 34
Trieste MaX Conference, 31. Jan. 2018
Pd Ir(111) Fe Pd Ir Ir Fe Fe Pd Pd Pd
Nature Comm. 7, 11779 (2016)
PAGE 35
Trieste MaX Conference, 31. Jan. 2018
Ir(111) Fe Pd Rh Ir Ir Fe Fe RhxPd1-x Rh
Nature Comm. 7, 11779 (2016)
PAGE 36
Tuning of the exchange Tuning of the DM
Trieste MaX Conference, 31. Jan. 2018
Dupe, Bihlmayer, Blügel, Heinze , Nature Comm. 7, 11779 (2016)
ESS(B) E
S k X
(B) EFM ( B )
PAGE 37
Trieste MaX Conference, 31. Jan. 2018
PAGE 38
Trieste MaX Conference, 31. Jan. 2018
GSky = GFM + GFM ∆ V GSky
Ir(111) Fe Pd
Spin-Polarized Scanning Tunneling Microscopy
Pd/Fe/Ir(111)
PAGE 39
Trieste MaX Conference, 31. Jan. 2018
* * * * * *
PAGE 40
Nature Comm. 6, 8541 (2015)
Trieste MaX Conference, 31. Jan. 2018
* * *
*
TXMR(r) = LDOSvac
FM − LDOSvac {S}
LDOSvac
FM
× 100%
PAGE 41
Trieste MaX Conference, 31. Jan. 2018
* *
Nature Comm. 6, 8541 (2015)
PAGE 42
Trieste MaX Conference, 31. Jan. 2018
Quantum Phenomena for the New Information Age § Spinorbitronics
§ Spintextures for neuro-inspired computing § Ultrafast and antiferromagnetic spintronics § 3D nanoscale magnetic textures & dynamics
§ Quantum materials
§ Emergent complex phase space topology § Topological superconductors for QC
§ Materials discovery lab – Computer
§ Cognitive Materials and Functionality Discovery
From Nicola Marzari PAGE 43
Trieste MaX Conference, 31. Jan. 2018
Daniel Wortmann Gustav Bihlmayer Gregor Michalicek Uliana Alekseeva Im Ma Juba Da Sa
KKRnano
Rudolf Zeller Roman Kovacik Marcel Bornemann Paul Baumeister Dirk Pleiter Jens Bröder Daniel Wortmann
PAGE 44