Solid Math, Aalborg May 28, 2016 1 ELECTRONIC TRANSPORT in APERIODIC SOLIDS Jean BELLISSARD Westfälische Wilhelms-Universität, Münster Department of Mathematics Georgia Institute of Technology, Atlanta School of Mathematics & School of Physics e-mail: [email protected]Sponsoring
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Solid Math, Aalborg May 28, 2016 1
ELECTRONIC TRANSPORTin
APERIODIC SOLIDSJean BELLISSARD
Westfälische Wilhelms-Universität, MünsterDepartment of Mathematics
Georgia Institute of Technology, AtlantaSchool of Mathematics & School of Physics
Main ReferencesJ. Bellissard, H. Schulz-Baldes, A. van Elst,The Non Commutative Geometry of the Quantum Hall Effect, J. Math. Phys., 35, 5373-5471, (1994).
J. Bellissard, H. Schulz-Baldes, J. Stat. Phys., 91, 991-1026, (1998).
H. Schulz-Baldes, J. Bellissard, Rev. Math. Phys., 10, 1-46 (1998).
J. Bellissard, Coherent and dissipative transport in aperiodic solids,Lecture Notes in Physics, 597, Springer (2003), pp. 413-486.
G. Androulakis, J. Bellissard, C. Sadel, J. Stat. Phys., 147, (2012), 448-486.
Y. Xue, E. Prodan, Noncommutative Kubo formula: Applications to transport in disorderedtopological insulators with and without magnetic fields, Phys. Rev. B, 86, 155445, (2012).
E. Prodan, J. Bellissard, Mapping the Current-Current Correlation FunctionNear a Quantum Critical Point, Ann. of Phys., 368, 1-15, (2016).
Solid Math, Aalborg May 28, 2016 4
Content1. Dissipation, Kubo’s Formula
2. Anomalous Transport
3. Numerics
Solid Math, Aalborg May 28, 2016 5
A No-Go TheoremLet H = H∗ be bounded (one-electron Hamiltonian),
Let ~R = (R1, · · · ,Rd) be the position operator(selfadjoint, commuting coordinates)
Then the electronic current is
~J = −eı~
[H, ~R] ,
Adding a force ~F at time t = 0 leads to a new evolution withHamiltonian HF = H − ~F · ~R.
Solid Math, Aalborg May 28, 2016 6
A No-Go TheoremThe 0-frequency component of the current is
~j = limt→∞
∫ t
0
dst
eısHF/~ ~J e−ısHF/~ ,
Solid Math, Aalborg May 28, 2016 7
A No-Go TheoremThe 0-frequency component of the current is
~j = limt→∞
∫ t
0
dst
eısHF/~ ~J e−ısHF/~ ,
Simple algebra shows that (since ‖H‖ < ∞)
~F · ~j = const. limt→∞
H(t) −Ht
= 0 ,
Solid Math, Aalborg May 28, 2016 8
A No-Go TheoremThe 0-frequency component of the current is
~j = limt→∞
∫ t
0
dst
eısHF/~ ~J e−ısHF/~ ,
Simple algebra shows that (since ‖H‖ < ∞)
~F · ~j = const. limt→∞
H(t) −Ht
= 0 ,
WHY ?
Solid Math, Aalborg May 28, 2016 9
A No-Go TheoremThis is called Bloch’s Oscillations. It was observed in simulationsusing ultracold atoms in an artificial lattice produced by lasers.
Solid Math, Aalborg May 28, 2016 10
A No-Go TheoremThis is called Bloch’s Oscillations. It was observed in simulationsusing ultracold atoms in an artificial lattice produced by lasers.
To get a non trivial current we need
DISSIPATION !
Namely loss of information.
Solid Math, Aalborg May 28, 2016 11
The Drude Model (1900)Assumptions :
1. Electrons in a metal are free classical particles of mass m∗ andcharge q.
2. Let n denotes the electron density.
3. They experience collisions at random Poissonnian times · · · < tn <tn+1 < · · ·, with average relaxation time τrel.
4. If pn is the electron momentum between times tn and tn+1, thenthe pn+1’s is updated according to the Maxwell distribution attemperature T.
Solid Math, Aalborg May 28, 2016 12
The Drude Model (1900)
The Drude Kinetic Model
Solid Math, Aalborg May 28, 2016 13
The Drude Model (1900)An elementary calculation leads to the Drude formula
σ =q2nm∗
τrel
Solid Math, Aalborg May 28, 2016 14
The Drude Model (1900)An elementary calculation leads to the Drude formula
σ =q2nm∗
τrel
Heat conductivity can also be computed leading to
λ =3n
2m∗k2
BT τrel
Solid Math, Aalborg May 28, 2016 15
The Drude Model (1900)An elementary calculation leads to the Drude formula
σ =q2nm∗
τrel
Heat conductivity can also be computed leading to
λ =3n
2m∗k2
BT τrel
The ratio gives the Wiedemann-Franz Law (1853)
λσ
=32
(kB
q
)2T
Solid Math, Aalborg May 28, 2016 16
Aperiodicity
1. If the charges evolve in an aperiodic environment, their one-particle Hamiltonian is actually a family (Hω)ω∈Ω of self-adjointoperators depending on a parameter ω characterizing the de-gree of aperiodicity (disorder parameter).
2. The aperiodicity can be ordered like in quasicrystals (long rangeorder), or disordered like in semiconductors, glasses or liquids(short range correlations).
3. The space Ω of the disorder parameters is called the Hull. It isalways compact and metrizable.
4. The translation group G acts on Ω by homeomorphisms ta, a ∈ G.
Solid Math, Aalborg May 28, 2016 17
Aperiodicity1. Covariance: if G is the translation group, if U(a) represents the
translation by a ∈ G in the Hilbert space of quantum states,then
3. Trace per Unit Volume: ifP is a G-invariant ergodic probabilityon Ω then, for P-almost every ω
TP(
f (H))
=
∫Ω
dP(ω)〈x| f (Hω) x〉 = limΛ↑Rd
1|Λ|
Tr(
f (Hω) Λ)
Solid Math, Aalborg May 28, 2016 18
A Quantum Drude ModelAssumptions :
1. Replace the classical dynamics by the quantum one with one-particle Hamiltonian H = (Hω)ω∈Ω.
2. Collisions occur at random Poissonnian times · · · < tn < tn+1 <· · ·, with average relaxation time τrel.
3. At each collision, the density matrix is updated to the equilib-rium one. (Relaxation Time Approximation).
4. Electrons and Holes are Fermions: use the Fermi-Dirac distribu-tion to express the equilibrium density matrix.
Solid Math, Aalborg May 28, 2016 19
A Quantum Drude ModelA straightforward calculation leads to the Kubo formula(JB, Schulz-Baldes, Van Elst ’94)
σi, j =q2
~TP
(∂ j
(1
1 + eβ(H−µ)
)1
1/τrel − LH∂iH
)
Where
Solid Math, Aalborg May 28, 2016 20
A Quantum Drude Model
1. ∂iA = ı[Ri,A] is the quantum derivative w.r.t. the momentum.
2. LH(A) = ı/~ [H,A] is called the Liouvillian.
3. β = 1/kBT and µ is the chemical potential fixed by the electrondensity, namely
n = TP
(1
1 + eβ(H−µ)
)4. TP denotes the trace per unit volume, where P provides the way
the average over the volume is defined.
Solid Math, Aalborg May 28, 2016 21
The Work of E. Prodan: Numerical ResultsE. Prodan, “Quantum transport in disordered systems under magnetic fields:a study based on operator algebras”, arXiv:1204.6490. Appl. Math. Res. Express, (2012)
Numerical implementation of the previous Kubo Formula for dis-ordered systems was provided by E. Prodan. The formula gives anaccurate algorithm which is very stable against disorder.
He used this algorithm to investigate more thoroughly the plateauxof conductivity in the Quantum Hall Effect (QHE) with his collab-orators after 2012.
Solid Math, Aalborg May 28, 2016 22
The Work of E. Prodan: Numerical ResultsQuantum Hall Effect:
DoS (left) and colored mapof the Hall conductivity(right) for W = 3.The regions of quantizedHall conductivity, whichappear as well definedpatches of same color, areindicated at the right.
Solid Math, Aalborg May 28, 2016 23
The Work of E. Prodan: Numerical ResultsHall Plateaux:
First row (Second row): Thediagonal and the Hallresistivities as function ofFermi energy (density) at fixedmagnetic flux φ, temperatureT and disorder strength W
φ = 0.1 h/ekBT = 1/τrel = 0.025W = 1, 2, 3.
Each panel compares the dataobtained on the 100 × 100lattice (circles) and on the120 × 120 lattice (squares).
Solid Math, Aalborg May 28, 2016 24
The Work of E. Prodan: Numerical ResultsMetal-Insulator transition between Hall plateaux
Transition fromChern(PF) = 0 toChern(PF) = 1
The simulated (a) σxy and (b)σxx, as functions of EF atdifferent temperatures.(Song & Prodan ‘12)
It shows a fixed point atEF = Ec
F where
σxxT↓0→ σxy = e2/2h
Solid Math, Aalborg May 28, 2016 25
The Work of E. Prodan: Numerical ResultsMetal-Insulator transition between Hall plateaux: resistivity
−4 −3.5 −3 −2.50
0.5
1
1.5
2
2.5
3
3.5
4
EF
kT=0.01kT=0.02kT=0.03kT=0.04kT=0.06kT=0.08kT=0.10
ρxx (h/e2)
ρxy (h/e2)
Transition fromChern(PF) = 0 toChern(PF) = 1
ρxy as function of EF atdifferent temperatures. Thecurves at lower temperaturesdisplay quantized values wellbeyond the critical point,which is marked by the verticaldotted line. For conveniencewe also show the data for ρxx.(Song & Prodan ‘12)
Solid Math, Aalborg May 28, 2016 26
The Work of E. Prodan: Numerical ResultsMetal-Insulator transition between Hall plateaux: Scaling Law
Transition fromChern(PF) = 0 toChern(PF) = 1
The simulated ρxx as functionof EF (a) before and (b) afterthe horizontal axis wasrescaled as:
EF → EcF + (EF − Ec
F)( TT0
)−κwith Ec
F = −3.15 , kBT0 = .08and κ = .2 leading to p = 1(Song & Prodan ‘12)
Solid Math, Aalborg May 28, 2016 27
What is Coherent Transport ?Coherent transport corresponds to charge transport (electrons orholes) ignoring dissipation sources such as electron-phonon or electron-electron interactions.
1. The independent electrons approximation is justified.
2. The one-particle Hamiltonian (Hω)ω∈Ω is sufficient.
3. The wave packets diffuses through the medium at a rate depend-ing upon how much Bragg reflections are produced (quantuminterferences).
Solid Math, Aalborg May 28, 2016 28
What is Coherent Transport ?
WAVE DIFFUSION
but
NO CURRENT !
Solid Math, Aalborg May 28, 2016 29
Spectral Measures
1. The density of state (DOS)∫ +∞
−∞
dNP(E) f (E) = TP(
f (H))
2. The spectral measure relative to a given state ψ in the Hilbertspace, called local density of state (LDOS)∫ +∞
−∞
dµω,ψ(E) f (E) = 〈ψ| f (Hω)ψ〉
3. The current-current correlation (CCC) describes transport prop-erties.∫
R2dm(E,E′) f (E)g(E′) =
d∑i=1
TP f (H) ∂iH g(H) ∂iH
Solid Math, Aalborg May 28, 2016 30
The Current-Current Correlation Measure: NumericsE. Prodan, J. Bellissard, Ann. of Phys., 368, 1-15, (2016).
The QHE scaling laws come from a singularity in the Current-Current Measure. Here dm(E,E′) = f (E,E′)dEdE′
Left: intensity plot of thecurrent-current correlationdistribution f (E,E′).Right: level sets of f (E,E′)The calculation was made on a120 × 120 lattice and the adtawere averaged over 100random configurations.(Prodan & Bellissard ‘16)
Solid Math, Aalborg May 28, 2016 31
The Current-Current Correlation Measure: NumericsThe scaling law observed in the QHE resistivity at the metal-insulatortransition Ec can be explained by an expression of the form:
f (E,E′) = g(E + E′ − 2Ec
|E − E′|κ/p
), E,E′ ' Ec .
It leads to
σii =e2
h
∫∞
0
4π1 + y2 g
(EF
(Γ|E − E′|)κ/p
)dy , Γ =
1τrel.
Solid Math, Aalborg May 28, 2016 32
The Current-Current Correlation Measure: NumericsThe function g(t) can be computed also and fits well with a Gaus-sian curve.
4π2 g
(t)
t
Left: The trace of theasymptotic region where thescaling invariance of thecurrent-current correlationfunction occurs.Right: Plot of 10 values of thefunction g(t), together with aGaussian fit. (Prodan &Bellissard ‘15)
Solid Math, Aalborg May 28, 2016 33
Local ExponentsGiven a positive measure µ on R:
α±µ(E) = lim
supinf
ε ↓ 0
ln∫ E + ε
E − εdµ
ln ε
For ∆ a Borel subset of R:
α±µ(∆) = µ−ess
supinf
E∈∆
α±µ(E)
Solid Math, Aalborg May 28, 2016 34
Local ExponentsProperties:
1. For all E, α±µ(E) ≥ 0. In addition, α±µ(E) ≤ 1 for µ-almost all E.
2. If µ is ac on ∆ then α±µ(∆) = 1, if µ is pp on ∆ then α±µ(∆) = 0.
3. If µ and ν are equivalent measures on ∆, then α±µ(E) = α±ν (E)µ-almost surely.
4. α+µ coincides with the packing dimension.α−µ coincides with the Hausdorff dimension.
Solid Math, Aalborg May 28, 2016 35
Local Exponents
1. The LDOS exponent α±LDOS is defined as the maximum over thestate ψ of the local exponent associated with µψ.
2. The DOS exponents α±DOS is the local exponent associated withNP.
3. It follows that
α±LDOS(∆) ≤ α±DOS(∆)
Solid Math, Aalborg May 28, 2016 36
Transport Exponents
1. For ∆ ⊂ R a Borel subset, let P∆, ω be the corresponding spectralprojection of Hω. Set
~Rω(t) = eıtHω ~R e−ıtHω
2. The averaged spread of a typical wave packet with energy in ∆is measured by
L(p)∆ (t) =
∫ t
0
dst
∫Ω
dP 〈x|P∆, ω|~Rω(t) − ~R|pP∆, ω|x〉
1/p
3. Define β = β±p (∆) similarly so that L(p)∆ (t) ∼ tβ
Solid Math, Aalborg May 28, 2016 37
Transport ExponentsProperties:
• β−p (∆) ≤ β+p (∆) and β±p (∆) are non decreasing in p.
• The transport exponent is the spectral exponent of the Liouvil-lian LH localized around energies in ∆ near the eigenvalue 0(diagonal of the current-current correlation).
Solid Math, Aalborg May 28, 2016 38
Transport ExponentsHeuristic:
1. β = 0→ absence of diffusion: (ex: localization)
5. For 0 ≤ βF < 1/2, σ ↓ 0 as T ↓ 0: the system behaves as aninsulator.
Solid Math, Aalborg May 28, 2016 42
Transport in QuasicrystalsLectures on Quasicrystals,F. Hippert & D. Gratias Eds., Editions de Physique, Les Ulis, (1994),S. Roche, D. Mayou and G. Trambly de Laissardiere,Electronic transport properties of quasicrystals, J. Math. Phys., 38, 1794-1822 (1997).
Quasicrystalline alloys :
Metastable QC’s: AlMn(Shechtman D., Blech I., Gratias D. & Cahn J., PRL 53, 1951 (1984))