Debye Theory: Discrepancies of Einstein Theory
Debye Theory: Discrepancies of Einstein Theory
Debye Theory: Discrepancies of Einstein Theory
Debye Theory
2
Debye Theory
Considers crystal lattice as harmonic oscillators coupled together
by interatomic forces.
Though not exact, Debye theory agrees with experiment
satisfactorily over a wide range of temperatures.
Works with the allowed modes of oscillations rather than the
individual atomic oscillators.
Internal Energy
Though the atomic oscillators are coupled together, the normal
modes are independent of one another.
T
Total
t l vibrational
ib ti l energy can be
b ffound
d bby fi
finding
di the
th internal
i t
l
energy associated with the normal modes having frequencies in a
range d about any given frequency , and integrating over
frequency.
The differential energy:
dU ( ) g ( )d
( ) : Mean energy
g ( )d : Number of modes in d about
4
Mean Energy
We derived the mean energy while discussing Einstein Theory of
specific heat:
1
2
( )
1
e
/ kbT
Density of States
In practice, it is easier to calculate g(k)dk instead of g()d.
dU ( ) g ( )d
( ) g (k )
dk
d
d
Density of States
Assume a rectangular crystal of dimensions Lx, Ly, Lz, with
volume V in periodic boundary conditions.
2 n y
2 nx
2 nz
Allowed propagation vectors: k x
, ky
, kz
Lx
Ly
Lz
2
Lz
2
Ly
2
Lx
8 3
8 3
Lx Ly Lz
V
7
Density of States
Spherical surfaces corresponding to constant values k and k + dk
kz
dk
4 k 2 dk k 2Vdk
g (k )dk
8 3 / V
2 2
ky
k
Considering one longitudinal
and two transverse modes
kx
3k 2Vdk
g (k )dk
2 2
8
dk/d
We have to know the dispersion relation.
For three-dimensional cases, the dispersion relation is not easily
expressed in closed form.
Debye assumptions:
Linear dispersion relation
k km
v0 k ,
v0 : Constant velocity
Cutoff frequency
4 3
8 3
km N
, or
3
V
1/3
km 6 2
V
Internal Energy
m
U dU
U
( ) g (k )
dk
d
d
3k 2V 1
/ k T
d
2
2 e b 1 2 v0
3V 1
2
2 v0
2
/ k T
2 d
2 e b 1 v0
3
3
0 2 e / kbT 1 d
3V m4
3V m 3
2 3
d
2 v0 8 2 2 v03 0 e / kbT 1
3V
2 3
2 v0
10
Internal Energy
3V m4 3V
U 2 3
2 v0 8 2 2 v03
e / kbT 1
4 3
8 3
6 2 N 6 2 Nv03
km N
V
3
V
km3
m3
Replacing
p
gV
9
9 N m 3 d
U N m 3
8
m 0 e / kbT 1
11
Specific Heat
U
cv
9
9 N m
3
N m 3
8
m 0 e / kb T 1
U
T
9N
3
m
9 N 2
3
m kbT 2
/ kbT
e / kb T
4 e / k T
b
e / kbT 1
d
kbT 2
12
Specific Heat
cv
Let
9 N 2
m3 kbT 2
4 e / k T
b
/ kb T
kT
x d b dx
kbT
Boundary values:
0x0
m x
T
cv 9 Nkb
m
, : Debye Temperature
kbT T
/T
x 4e x
e x 1
dx
13
T
9
9 N m 3
9
T
U N m 3
N
9
Nk
T
m
b
/
k
T
8
m 0 e b 1 8
Simplification:
Si lifi i
/T
x3
dx
ex 1
x3
x3
x2
x
e 1 1 x 1
9
T
U N m 9 NkbT
8
/T
/T
3
9
T x
x dx N m 9 NkbT
8
3 0
3
9
T 9
N m 3NkbT N m 3NkbT
8
T 8
3
cv 3 Nkb
14
T0
9
9 N m 3
9
T
N m 9 NkbT
U N m 3
/
k
T
b
1 8
8
m 0 e
/T
x3
dx
ex 1
Simplification:
Si lifi i
1
1 1
x
x
e 1 e 1 e x
/T
1
x
e
e x e x e nx
n0
n 1
n 1
3
x
x3
dx x
dx x3e nx dx
x
0 e 1
0
e 1
n 1
n 1
3 nx
xe
dx I n
n 1
15
Integration
In x e
0
3 nx
nx
e nx
3
2 e
dx x 3
3
x
dx 0 x 2 e nx dx
0
n 0
n
n 0
e nx
3
3 e nx
6
dx 0 2
2x
x2
n
n 0 n 0
n
n
6
2
n
e nx
6
x n n2
6
3
n
e nx
6
n n4
e nx
6
dx 0 3
n
n
xe nx dx
xe nx dx
16
T0
x3
1
dx
6
4
x
e 1
n 1 n
/T
Riemann
Ri
zeta function:
f
i
1
1 1 1
4
(4) 4 1 4 4 4
2 3 4
90
n 1 n
9
T
U N m 9 NkbT
8
/T
9
3 4
T
N m
NkbT
8
5
12 4
T
cv
Nkb
5
4
x3
9
T 6
dx N m 9 NkbT
ex 1
8
90
3
17