MTLE-6120: Advanced Electronic Properties of Materials Electron - - PowerPoint PPT Presentation

mtle 6120 advanced electronic properties of materials
SMART_READER_LITE
LIVE PREVIEW

MTLE-6120: Advanced Electronic Properties of Materials Electron - - PowerPoint PPT Presentation

1 MTLE-6120: Advanced Electronic Properties of Materials Electron transport: phonons and electron-phonon scattering Contents: Phonon density of states and heat capacity Thermal conductivity: electronic and lattice Electron-phonon


slide-1
SLIDE 1

MTLE-6120: Advanced Electronic Properties of Materials Electron transport: phonons and electron-phonon scattering

Contents:

◮ Phonon density of states and heat capacity ◮ Thermal conductivity: electronic and lattice ◮ Electron-phonon momentum relaxation time

Reading:

◮ Kasap: 4.10

1

slide-2
SLIDE 2

Phonons in 1D

◮ 1D chain of atoms of mass M connected by springs K ◮ Frequencies ω = 2

  • K

M

  • sin ka

2

  • ◮ Group velocity at k → 0 is vL = a
  • K

M , the sound velocity

2(K/M)1/2

  • π/a

+π/a Frequency ω Wavevector k

2

slide-3
SLIDE 3

Phonons in 3D (eg. Au)

10 20 Γ X W L Γ K Energy [meV]

◮ Phonons have a band structure, just like electrons ◮ Key difference: linear near k = 0(Γ) rather than quadratic ◮ In 3D isotropic mateirals: two sound speeds: longitudinal (vL) and

transverse (vT )

◮ Two tranverse modes degenerate near k = 0 only (three polarizations) ◮ Compare energy scale against electrons: factor of 100 smaller

3

slide-4
SLIDE 4

Phonon density of states

◮ Number of states per unit volume in

k per unit volume =

1 (2π)d ◮ Phonon dispersion relations ω = ωn(

k) (many bands)

◮ Therefore density of states g(ω) = n

  • d

k (2π)d δ(ω − ωn(

k))

◮ Need a simple model to evaluate analytically (analogous to free electrons)

4

slide-5
SLIDE 5

Debye model

◮ How many total phonon states in one band per unit cell? One! ◮ Debye model: keep linear dispersion with correct total number of states ◮ If ω = v|

k| for one band, density of states g(ω) =

  • d

k (2π)d δ(ω − v| k|) = xdkd−1dk (2π)d δ(ω − vk) = xdwd−1dw (2πv)d δ(ω − w) (w = vk) = xdωd−1 (2πv)d = 4πω2 (2πv)3 (in 3D)

5

slide-6
SLIDE 6

Debye frequency

◮ For linear phonon dispersion in 3D

g(ω) = 4πω2 (2πv)3

◮ If this is true for all ω, then total number of modes is

  • dωg(ω) = ∞

◮ Debye model: cutoff model at ωD to keep correct number of states ◮ Number of states per unit volume (one band) = nion (# ions / volume) ◮ Therefore impose condition:

nion = ωD dωg(ω) = ωD dω 4πω2 (2πv)3 = 4πω3

D

3(2πv)3 = ω3

D

6π2v3 ⇒ ωD = v(6π2nion)1/3

◮ Corresponding cutoff in wave vector, kD = ωD/v = (6π2nion)1/3 ◮ Note similarity to kF = (3π2n)1/3 for electrons

6

slide-7
SLIDE 7

Phonon density of states: real metals

60 120 180 240 300 DOS [1029 eV-1 m-3]

a) Al

  • Eq. 4 (this work)

Debye 60 120 180 240 300

b) Ag

60 120 180 240 0.02 0.04 0.06 DOS [1029 eV-1 m-3]

ε [eV]

c) Au

60 120 180 240 0.02 0.04 0.06

ε [eV]

d) Cu

◮ Using two Debye models for velocities vL and vT ◮ Note energy scale ED = kBTD = ωD

  • Phys. Rev. B 91, 075120 (2016)

7

slide-8
SLIDE 8

Bose statistics

◮ At temperature T each phononic state of energy E has average occupation

nph(E) = 1 exp

E kBT − 1 ◮ For massless bosons like photons and phonons, µ = 0 because number not

constrained

1 µ-2kBT µ µ+2kBT E f(E) Bose Classical Fermi

8

slide-9
SLIDE 9

Lattice heat capacity

◮ Density of phonon modes g(ω) = 4πω2 (2πv)3 Θ(ωD − ω) (Debye model) ◮ Phonon modes occupied by Bose function nph(ω) at temperature T ◮ Internal energy per unit volume:

U = ∞ dω(ω)g(ω)nph(ω) = ∞ dω(ω) 4πω2 (2πv)3 Θ(ωD − ω) 1 exp

ω kBT − 1

=

  • 2π2v3

ωD ω3dω exp

ω kBT − 1

9

slide-10
SLIDE 10

Lattice heat capacity: high T (T ≫ TD)

◮ At high temperatures ω ≪ kBT ⇒ exp ω kBT − 1 ≈ ω kBT ◮ Correspondingly, internal energy is approximately

U ≈

  • 2π2v3

ωD ω3dω

ω kBT

= 1 2π2v3 ωD ω2dωkBT = 1 2π2v3 ω3

D

3 kBT = k3

D

6π2 kBT = nionkBT

◮ Accounting for the three polarizations, U = 3nionkBT and CV = 3nionkB ◮ This is (of course) the equipartition result: high T → classical limit

10

slide-11
SLIDE 11

Lattice heat capacity: low T (T ≪ TD)

◮ At low temperatures,

U ≈

  • 2π2v3

ωD ω3dω exp

ω kBT − 1

=

  • 2π2v3

kBT

  • 4

ωD kBT

x3dx ex − 1 x ≡ ω kBT ≈ k4

BT 4

2π23v3 ∞ x3dx ex − 1 ωD kBT = TD T ≫ 1 = k4

BT 4

2π23v3 · π4 15 = π2k4

BT 4

303v3

◮ Accounting for one longitudinal and two transverse velocities:

U = T 4 π2k4

B

303 1 v3

L

+ 1 v3

T 1

+ 1 v3

T 2

  • ∝ T 4

◮ And corresponding heat capacity

CV ≡ dU dT = T 3 2π2k4

B

153 1 v3

L

+ 1 v3

T 1

+ 1 v3

T 2

  • ∝ T 3

11

slide-12
SLIDE 12

Lattice heat capacity: real metals

Debye model does a good job! Proportional to T 3 for T ≪ TD, constant 3nionkB for T ≫ TD (classical equipartition result)

10 20 30 40

Cl [105 J/m3K]

a) Al

  • Eq. 5 (this work)

Debye 10 20 30 40

b) Ag

10 20 30 0.5 1 1.5 2

Cl [105 J/m3K] Tl [103 K]

c) Au

10 20 30 0.5 1 1.5 2

Tl [103 K]

d) Cu

  • Phys. Rev. B 91, 075120 (2016)

12

slide-13
SLIDE 13

Heat capacity: electrons vs lattice

Mostly lattice! Except at very low T (< 10 K) when T wins over T 3

102 103 104 105 106 C [J/m3K] a) Al Total Lattice Electronic 102 103 104 105 106 b) Ag 102 103 104 105 106 100 101 102 103 C [ J/m3K] T [K] c) Au 102 103 104 105 106 100 101 102 103 T [K] d) Cu

13

slide-14
SLIDE 14

Electronic thermal conductivity setup

◮ Energy current due to one electron E

v (contrast with −e v for charge current)

◮ At energy E, energy flux g(E)f(E)(Ev(E)) in random directions,

averages to zero

◮ Non-uniform temperature ⇒ average energy flow in one direction ◮ Assume constant temperature gradient T(x) = T0 + x dT dx ◮ What is the net flux of energy across x = 0?

qx =

  • dx
  • dEg(E)f(E, T(x))(Ev(E))

· 1

2πd cos θ 4π

· e−(−x/ cos θ)/λ

λ

, x < 0 −

−1 2πd cos θ 4π

· e−(−x/ cos θ)/λ

λ

, x > 0

14

slide-15
SLIDE 15

Electronic thermal conductivity derivation

qx =

  • dx
  • dEg(E)f(E, T(x))(Ev(E))

1

d cos θ 2

· ex/(λ cos θ)

λ

, x < 0 −

−1 d cos θ 2

· ex/(λ cos θ)

λ

, x > 0 =

  • dx
  • dEg(E)
  • f(E, T0)
  • →0 (equilibrium)

+∂f(E, T ∂T xdT dx

  • (Ev(E))

· 1

d cos θ 2

· ex/(λ cos θ)

λ

, x < 0 −1

d cos θ 2

· ex/(λ cos θ)

λ

, x > 0 = dT dx

  • dEg(E)v(E)E ∂f(E, T)

∂T

  • 1

d cos θ 2 −∞ dxx ex/(λ cos θ) λ

− −1

d cos θ 2

∞ dxx ex/(λ cos θ)

λ

  • = dT

dx

  • dEg(E)v(E)E ∂f(E, T)

∂T

1

d cos θ 2

λ cos2 θ − −1

d cos θ 2

λ cos2 θ

  • = −dT

dx

  • dEg(E)v(E)E ∂f

∂T λ/3 κ = τ 3

  • dEg(E)v2(E)E ∂f(E, T)

∂T (λ = vτ)

15

slide-16
SLIDE 16

Electronic thermal conductivity: result

Substitute Fermi function in thermal conductivity expression: κ = τ 3

  • dEg(E)v2(E)E ∂f(E, T)

∂T ≈ v2

F τ

3

  • dEg(E)E ∂f(E, T)

∂T ∂f ∂T sharply peaked

  • = v2

F τ

3 CV = v2

F τ

3 · g(EF )π2k2

BT

3 = v2

F g(EF )τ

3 · π2k2

BT

3 Note similarity to electrical conductivity σ = v2

F g(EF )τ

3

· e2. Therefore expect Lorenz number L ≡ κ σT = π2k2

B

3e2 ≈ 2.44 × 10−8(V/K)2 to be same temperature-independent constant across metals (not just free-electron-like ones)

16

slide-17
SLIDE 17

Wiedemann-Franz law

L ≡ κ σT = π2k2

B

3e2 ≈ 2.44 × 10−8(V/K)2 Metal L at T = 100 K (V/K)2 L at T = 273 K (V/K)2 Copper 1.9 × 10−8 2.3 × 10−8 Gold 2.0 × 10−8 2.4 × 10−8 Aluminum 1.5 × 10−8 2.2 × 10−8 Zinc 1.8 × 10−8 2.3 × 10−8 Lead 2.0 × 10−8 2.5 × 10−8

17

slide-18
SLIDE 18

Lattice thermal conductivity

◮ Phonons cannot transport charge, but they can transport heat. ◮ Our derivation for electrons didn’t assume electrons till this point:

κe = τ 3

  • dEg(E)v2(E)E ∂f(E, T)

∂T

◮ Lattice contribution: phonon DOS and Bose occupations instead

κL = τph 3

  • dωg(ω)v2(ω)ω ∂nph(ω, T)

∂T

◮ Which is larger? ◮ For semiconductors / insulators: few electrons (and holes) ⇒ κL dominates ◮ For metals:

◮ n ≈ nion ◮ vF ∼ 106 m/s ≫ vL, vT ∼ 103 − 104 m/s (electrons win) ◮ τph ∼ ps, while τe ∼ 10 fs (phonons win)

◮ Net result: κe dominates at room temperature, κL important at high T ◮ Small κL ⇒ Lorenz number (based only on κe) close to ideal

18

slide-19
SLIDE 19

Electron-phonon scattering rate

◮ All transport coefficients depend on scattering time τ ◮ For electrons: τ determined by electron-phonon scattering ◮ Crude argument during Drude discussion:

◮ Electrons only scatter against displaced ions (now we know why: band

theory)

◮ Cross-section ∝ mean-squared displacement ∝ kBT (equipartition) ◮ Therefore τ −1 ∝ T, and ρ =

m ne2τ ∝ T

◮ Why is this not true at low T? ◮ Equipartion no longer valid for T ≪ TD

100 200 300

T emperature [K]

1 2 3 4

Resistivity [10-8 m]

P u r e c

  • p

p e r C

  • l

d w

  • r

k e d + 1 % N i + 2 % N i

19

slide-20
SLIDE 20

Electron-phonon scattering rate estimation

◮ Process: electron in one state absorbs (or emits) a phonon to go into

another state

◮ Phonon absorption rate of electron in state i (Fermi’s Golden rule):

τ −1 = 2π

  • j

|j|H′|i|2δ(Ej − Ei − ω) (using j for final state to prevent mix up with Fermi energy EF )

◮ Ingredient 1: sum over states j ◮ Count up electronic states using

  • dEjg(Ej)

◮ State must be unoccupied (Fermi statistics): (1 − f(Ej))

(was not relevant before because only one electron)

◮ Also need to sum over phonons available to absorb:

  • dωg(ω)nph(ω)

◮ Net result:

  • j

  • dEjg(Ej)(1 − f(Ej))
  • dωg(ω)nph(ω)

20

slide-21
SLIDE 21

Electron-phonon scattering: selection rule

◮ So far assumed that final states counted as:

  • j

  • dEjg(Ej)(1 − f(Ej))
  • dωg(ω)nph(ω)

◮ But many pairs of final states will not satisfy

ki + kph = kj

◮ Instead of all pairs, only count those with correct momentum ◮ RHS on sphere of radius kj ≈ ki (phonon energies negligible) ◮ LHS on sphere of radius kph centered on

ki

◮ Total pairs of

kj, kph: (4πk2

j) · (4πk2 ph) ◮ Momentum conserving set: 2πkph ◮ Fraction conserving momentum:

2πkph (4πk2

j) · (4πk2 ph) =

1 8πk2

jkph

≈ vph 8πk2

i ω ◮ Net final state counting:

  • j

  • dEjg(Ej)(1 − f(Ej))
  • dωg(ω)nph(ω) vph

8πk2

i ω

21

slide-22
SLIDE 22

Electron-phonon scattering: momentum relaxation

◮ Does electron-phonon scattering randomize momentum?

(We assumed this about τ in Drude theory and its quantum extension.)

◮ Change in angle θ between

ki and kj given by cos θ = k2

i + k2 j − k2 ph

2kikj ≈ 1 − k2

ph

2k2

i ◮ Extent of randomization proportional to:

1 − cos θ ≈ ω2 2k2

i v2 ph ◮ To calculate momentum-relaxtion time instead of

mean-free time, include this (1 − cos θ) factor:

  • j

  • dEjg(Ej)(1 − f(Ej))
  • dωg(ω)nph(ω) vph

8πk2

i ω ·

ω2 2k2

i v2 ph

22

slide-23
SLIDE 23

Electron-phonon scattering: matrix element

◮ Matrix element:

j|H′|i =

  • d

rψ∗

j (

r)∆Vph( r)ψi( r) where ∆V ( r) is change in electron potential due to phonon

◮ For estimate, make similar to argument made in classical case ◮ Instead of equipartition, apply to single phonon mode ◮ Mean-squared atom displacement ∝ ω ◮ ∆V (

r) ∝ displacement of atoms (not squared)

◮ Therefore assume

j|H′|i ∝ √ ω ≈

  • ce-phω(say)

(constant of proportionality ce-ph has dimensions of energy)

23

slide-24
SLIDE 24

Electron-phonon scattering: Fermi golden rule

◮ Collecting ingredients together:

τ −1 ≈ 2π

  • dEjg(Ej)(1 − f(Ej))
  • dωg(ω)nph(ω)

· vph 8πk2

i ω ·

ω2 2k2

i v2 ph

·

  • ce-phω

2 δ(Ej − Ei − ω)

◮ Neglect phonon energies ω ≪ Ej, Ei (two orders smaller)

τ −1 ≈ 2π

  • dEjg(Ej)(1 − f(Ej))
  • dωg(ω)nph(ω)

· ce-phω2 16πk4

i vph

δ(Ej − Ei) = 2π g(Ei)(1 − f(Ei))

  • dω g(ω)
  • ∝ω2

nph(ω) ce-ph 16πk4

i vph

ω2 ∝ g(Ei) ωD dω ω4 exp

ω kBT − 1

24

slide-25
SLIDE 25

Electron-phonon scattering: temperature dependence

◮ For high temperatures T ≫ TD, exp ω kBT − 1 ≈ ω kBT :

τ −1 ∝ g(EF ) ωD dω ω4

ω kBT

∝ g(EF )T

◮ For low temperature T ≪ TD,

τ −1 ∝ g(EF ) ωD dω ω4 exp

ω kBT − 1

≈ g(EF ) kBT

  • 5 ∞

dx x4 ex − 1 ∝ g(EF )T 5

25

slide-26
SLIDE 26

Electron-phonon scattering: conclusions

◮ Temperature dependence:

ρ = 3 g(EF )e2v2

F

× τ −1 ∝ v−2

F

  • T,

T ≫ TD T 5, T ≪ TD explains low temperature deviations

◮ Reduction from classical result because few phonon quanta

  • ccupied at low T (similar to blackbody spectrum reason)

◮ Transition metals ⇒ low d-band velocities at EF ⇒ high resistivity

100 200 300

T emperature [K]

1 2 3 4

Resistivity [10-8 m]

P u r e c

  • p

p e r C

  • l

d w

  • r

k e d + 1 % N i + 2 % N i

26

slide-27
SLIDE 27

DFT + Fermi Golden rule

Quantitative prediction of resistivities

Metal DFT ¯ τ [fs] DFT ρ [Ωm] Expt ρ [Ωm] Al 11.5 2.79 × 10−8 2.71 × 10−8 Cu 35.6 1.58 × 10−8 1.71 × 10−8 Ag 36.4 1.58 × 10−8 1.62 × 10−8 Au 26.3 2.23 × 10−8 2.26 × 10−8

Far from isotropic even for nominally free electron metals

ACS Nano 10, 957 (2016)

27