Homework 5

Homework # 5
Chapter 5 Kittel
Phys 175A
Dr. Ray Kwok
SJSU
Prob. 1 Singularity in density of states
Victor Chikhani
(a) From the dispersion relation derived in Chapter 4 for a monatomic

linear lattice of N atoms with nearest-neighbor interactions, show that
the density of states is D( ) = 2N 1
1
( )2
m
2 2
where m is the maximum frequency.

Dispersion relation: (K ) = sin Ka 2
K ( ) = sin 1 Na dK
D( ) =
m
2 a m d
Therefore, dK 2 1
=
d a 2 2 1 2
(m )
2N 1
And, D( ) =
2 2 12
(m )
(b) Suppose that an optical phonon branch has the form
(K)= 0-AK2, near K=0 in three dimensions. Show that
L 2
3
for <0 and D()=0
1
D( ) = ( )
2 2 2
2 A 2
3 o
for >0 . Here the density of modes is discontinuous.
From the dispersion relation: (K ) = o AK 2

we get K ( ) = 11 ( 0 )12
2
A

and dK
=
1
1
1
1

d 2A 2 ( o ) 2
L 3 dK
For 3-D we know D( ) = (4 K 2 )
2 d
L 3 4 ( 0 ) 1 1
3
L 2 1
Therefore, D( ) = 1 1
= 3 o( ) 2
2 A A 2 2( o ) 2 2 A 2
Prob. 2 RMS thermal dilation of crystal cell
Daniel Wolpert
(a) Estimate for 300 K the RMS thermal dilation V/V for a primitive cell of sodium. Take the bulk modulus
as 7x1010 erg cm-3. Note that the Debye temperature 158 K is less than 300 K, so that the thermal energy is
of the order kBT.
kBT = B a3 (V/V)** T = 300K
V/V = kBT/Ba3 B = 7x1010 erg cm-3
V/V = .125 a = 4.225 x10-10 m

kB = 1.38062x10-16 erg/K
(b) Use this result to estimate the RMS thermal fluctuation a/a of the lattice parameter
V = a3
V/V = 3a/a ** B = P/(V/V)
U = PdV = W (Q=0, adiabatic, no external heat)
a/a = .125/3 = .04 or U = PV
average U = kT (longitudinal only bulk)
V = a3 (bcc, 2 atoms per cubic cell)
Prob. 3 Zero point lattice
displacement and strain
Adam Gray
A) In the Debye approximation, show that the mean square

displacement of an atom at absolute zero is <R2> =
3d2 / 8v3 , where v is the velocity of sound. Start
from the result (4.29) summed over the independent
lattice modes: <R2> = ( / 2V) -1 . We have included
a factor of to go from mean square amplitude to mean
square displacement.
<R2> = ( / 2V) -1
-1 can be expressed as:
d
0
dD() 1
Where d is the cutoff frequency, and D() is the density of states.
Using the following two equations in Kittel:
(5.20) D() = (VK2/22)(dK/d)
(5.21) = vK
We end up with D() = (VK2/22v)
Plugging this into the integral above results in:
-1 = Vd2 / 42v3
Plugging this into the above expression for <R2> :

<R2> = ( / 2V) (Vd2 / 42v3)
<R2> = ( d2 / 82v3)
The final missing factor of 3 comes from the phonon velocity being
independent of the polarization according to Kittel.
B) Show that -1 and <R2> diverge for a one-dimensional
lattice, but the mean square strain is finite. Consider
<(dR/dx)2)> = (1/4) K2u02 as the mean square strain,
and show that it is equal to d2L/4MNv3 for a line of N
atoms each of mass M, counting longitudinal modes
only. The divergence of R2 is not significant for any
physical measurement.
I think the equation for mean strain should be instead of

Consider R ~ uocos(t-kx), dR/dx would give you kR.
Averaging would give for each time average and spatial average.
To show the divergence, Equation (5.15) gives that
D() = L/v
If we use this to evaluate the same integral for -1 from
earlier, the result is
-1 = (L/v)(Ln(d)-Ln(0))
Which diverges at low values, since Ln(0) goes to .
Starting with <(dR/dx)2)> = (1/4) K2u02
Similar to part a) we can replace the sum by an integral over
the frequencies with K2u02 multiplied by D():
d

0 dD()K2u02(1/4)
From Kittel we have
(4.29) u02 = 4(n+1/2)/V
(5.15) D() = L/v
(5.21) K = /v
Plugging these in to the above integral yields
d
d(L/v)(/v)2(2/V)(1/4)
d
(L/Vv ) 0 d
0
3
Which results in: (Ld2/4Vv3)

But, V=NM for the given problem.
So: <(dR/dx)2)> = (1/4) K2u02 = (Ld2/4v3)
Prob. 4 Heat capacity of layer lattice
Gregory Kaminsky
Heat capacity of layer lattice. (a) Consider a dielectric

crystal made up of layers of atoms, with rigid coupling
between layers, so that the motion of the atoms is
restricted to the plane of the layer. Show that the phonon
heat capacity in the Debye approximation in the low
temperature limit is proportional to T2.
(b) Adjacent layers are now weakly bound, what form
would the phonon heat capacity approach at extremely
low temperatures.
Energy is:
hw K , p
U= exp(hw
K p / ) 1
K,p
hw
U= p
dwDp ( w) exp(hw / ) 1
D(w) is the number of modes per unit frequency
range , also called the Density of states.
Apply periodic boundary conditions over N2 primitive cells
within a square of side L
K is determined by
exp[i ( K x x + K y y )] = exp[i ( K x ( x + L) + K y ( y + L)]
Whence 2 4 N
K x , K y = 0; ; ; ;
L L L
Therefore there is one allowed value

of K per area (2/L)^2 in K space
or area per K is (L/ 2)2 = A/42
The total number of modes with
wavevector less then K is (L/2)2
times the area of the sphere of
radius K.
N = (L/ 2)2K2
Density of states = number of

modes per unit frequency range
The density of states for each polarization is D(w) = dN/dw =

2*(L/ 2)2K(dK/dw) = (KA/2)*dK/dw
w = vK
D(w) = wA/(2v2)
If there are N primitive cells in the specimen, the total
number of acoustic phonons is N. A cutoff frequency wd is
determined by using equation for N.
N = (L/ 2)2K2 = (L/ 2v)2w2
wd = (2v/L)*(N/ )1/2
For each polarization type, energy is given by:
hw
U = dwD( w)
exp( hw / ) 1
wD wA hw
U= 0
dw
2v 2 exp( hw / ) 1
Assuming that phonon velocity is independent of the
polarization type, we multiply by the factor of 2.
Ah wD w2
U=
v 2 0
dw
exp(hw / ) 1 Clat =
U
T
hw hw D
With x= And xD = = /T
k BT kBT
3 3 2
Ak T xD x

B
U= 2 2 dw
v h 0 exp( x ) 1
3 2 2
U 3 Ak T xD x

B
Clat = = 2 2 dw
T v h 0 exp( x ) 1
hw D
At low temperature limit xD = = /T
kBT
xD
2
x
0dw exp( x) 1 This is some finite constant.
Thus phonon heat capacity is proportional to T2.

Adjacent layers are now weakly bound, what form would the phonon heat
capacity approach at extremely low temperatures.
If adjacent layers are the only thing that can

vibrate (at low T), then the problem now is
a 1D and phonon heat capacity will
approach T at low temperature,
extrapolating from 2D and 3D models.
[in-plane coupling is relatively too stiff so the only motion

would be between planes.]
Prob. 5 Gruneisen Constant
Michael Tuffley

Homework 5

Uploaded by

Copyright:

Available Formats

Homework 5

Uploaded by

Document Information

Original Description:

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Homework 5

Uploaded by

Copyright:

Available Formats

What is the dispersion relation derived in Chapter 4 for a monatomic linear lattice?

What is the dispersion relation derived in Chapter 4 for a monatomic linear lattice?

What is the density of states expression derived from the dispersion relation?

What is the density of states expression derived from the dispersion relation?

Homework # 5

(a) From the dispersion relation derived in Chapter 4 for a monatomic

where m is the maximum frequency.

for >0 . Here the density of modes is discontinuous.

From the dispersion relation: (K ) = o AK 2

kBT = B a3 (V/V)** T = 300K

V/V = kBT/Ba3 B = 7x1010 erg cm-3

V/V = .125 a = 4.225 x10-10 m

A) In the Debye approximation, show that the mean square

Plugging this into the above expression for <R2> :

I think the equation for mean strain should be instead of

Which results in: (Ld2/4Vv3)

Heat capacity of layer lattice. (a) Consider a dielectric

Therefore there is one allowed value

Density of states = number of

The density of states for each polarization is D(w) = dN/dw =

Thus phonon heat capacity is proportional to T2.

If adjacent layers are the only thing that can

[in-plane coupling is relatively too stiff so the only motion

You might also like