Mod4.2 - Electrons, Phonons, Photons
Mod4.2 - Electrons, Phonons, Photons
Mod4.2 - Electrons, Phonons, Photons
Venkat Viswanathan
May 5, 2015
Key Concepts:
Bose and Fermi statistics, bosons and fermions, electron gas, Fermi
energy, Fermi momentum, Fermi temperature, electron pressure, electron heat capacity, crystal lattice, lattice vibrations, phonons, Einstein
model, Debye model, Black Body radiation, statistical mechanics of photons, Bose occupation, Plancks law.
(1)
Q(T , V , N ) exp(N )
N =0
!
X
N =0 n1 ,n2 ,...
!
exp
n +
n N
!
n
Y X
n1 ,n2 ,...
en n
pV = kB T log = kB T
log
n n
(2)
e +
n
n = 0
1
1 e +
(3)
If the particles are Fermions, any given state can only have either
n = 0 or n = 1 particles, thus:
X
n =0
e +
n
= 1 + e +
(4)
(5)
(6)
e + 1, we have:
X
hN i
e + = e q and pV kB T e q = kB T hN i
X
X
1 N N
= epV /(kB T ) = exp e q =
Q(T , V , N )eN
e
q =
N!
N =0
thus Q =
1 N
N! q
N =0
h2
n2x + n2y + n2z
2
8mL
translational modes
which are the energy levels for a particle in a box for a particles with
mass m = 9.10938 1031 kg.
The average number of electrons in the ~n state is:
1
= F (E~n )
(7)
e (E~n ) + 1
h
i1
where we have defined the Fermi function F () = e () + 1
.
The total number of electrons is given by:
hn~n i =
hN i = 2
X
X
X
nx =1 ny =1 nz =1
1
e (E~n ) + 1
=2
X
X
F (E~n )
(8)
nx =1 ny =1 nz =1
0
Z
Z
Z
Z
3
2V
L 1
dkx
dky
dkz F [(~k )] =
d~kF [(~k )]
=2 3 3
2
(2 )3
momentum pF according to F =
found to be:
hN i =
p2F
2m
2
}2 k F
2m
. Thus, at T = 0, hN i is
2V 4 3
8 V (2m)3/2 3/2
kF =
F
F =
3
(2 ) 3
3
h3
3
8
2/3
h2
2m
e () + 1
1
(10)
2 2
k
where we have used = }2m
.
In the limit T 0(or ), we have 1 + e () e (F ) for
< F and 1 + e () 1 for > F .
Therefore, the pressure is written as:
Z
16V (2m)2/3 5/2
4V (2m)3/2 F
1/2
d
(
)
=
F
(11)
pV =
F
h3
15h3
0
e +
1 + e +
(12)
}2 k2 .
2m
V (2m)2/3 3/2
F
where
dF
d
e ()
= [e () +1]2 .
Figure 2:
the Fermi
The derivative
dF ()
function d
e () [e () + 1]
of
=
(15)
Therefore, the heat capacity for a metal in the limit of small temperature is given by:
CV =
3
5 2 T
2
T
hN iF
2 2 =
hN ikB
5
12 F
2
F
(16)
Crystalline Solid
The atoms in a crystal are arranged in a regular array of points in space
with a variety of possible crystalline lattices (Fig. 4). At zero temperature, the atomic coordinates are uniquely locked into spatial positions
that minimize the energy. At finite temperature, the atoms fluctuate
about the energy-minimum positions, leading to lattice vibrations that
govern the thermodynamic behavior (Fig. 5).
(0)
(0)
(0)
(0)
~ri ~rj V |0 : ~ri ~ri
~rj ~rj
+ ...
i=1 j =1
The linear term is zero (by definition), thus the energy is:
V V0 +
N N
1 XX X
si sj Kij
2
, =x,y,z
(17)
i=1 j =1
(0)
where si = e ~ri ~ri
and Kij = ~ri ~rj V |0 .
The Matrix Kij can be diagonalized into normal modes (eigenvectors) with effective elastic constants (eigenvalues). Since there are
3N 6 3N degrees of freedom, there are 3N normal modes.
The potential energy is written as:
3N
V = V0 +
1X 2
Kl l
2
(18)
l =1
(19)
l =1
where pl and m
l are the effective momentum and mass of the lth normal
mode, respectively.
The total energy E is decomposed into normal modes with an individualmode energy in the form of a harmonic oscillator. These normal modes
are called phonons.
Phonons act as quasi-particles which means they are distinguishable
and independent, i.e. they dont interact between each other.
The Hamiltonian (sum of kinetic and potential energy) of the lth
p2
(20)
...
j1 =0 j2 =0
= eV0
V0
P3N
l=1
l
K
m
l.
(jl + 12 )}l
j3N =0
3N X
e (jl + 2 )}l
l=1 jl =0
= eV0
1
3N
Y
e 2 }l
1 e }l
(21)
l =1
n=0 z
1
1z .
X
log Q
hEi =
= V0 +
l =1
1
}l e }l
+ }l
2
1 e }l
The next step is to analyze two competing models for the phonon
frequencies l .
Einstein model
In the Einstein model, we assume there is a single characteristic frequency of the crystal, defined as the Einstein frequency E . The average
energy for this model is given by:
hEi = V0 +
3N }E
3N }E e }l
+
2
1 e }E
hEi
T
V ,N
2
2 E
}
k E
E
E
BT
e
e T
3N
k
3N kB k}B
B
T
T
=
=
2
2
}E
E
kB T
T
1
e
1e
Debye model
The Debye model treats the solid as an elastic material. Vibrational
modes in an elastic solid correspond to sound waves, thus the frequencies
satisfy = ck, where c is the speed of sound in the solid and k = m/L,
where m = 1, 2, ...
m1 m2 m3
3 Z kc
Z kc
Z kc
L
dk1
dk2
dk3 (...)
=
0
0
0
3 Z kc
3
Z
L
L
1 D
= 4
dkk 2 (...) = 4
d 2 (...)
3
c
0
0
where kc is a cutoff wavemode (to be determined), and D = Ckc is
called the Debye frequency.
A complete conversion will include one longitudinal mode with cl
and two transverse modes with ct . This gives:
Z D
3
X
L
1
2
+ 3
d 2 (...)
(22)
(...) = 4
c3l
ct
0
l
P
To find D we must enforce that l 1 = 3N , thus we have:
3
Z D
X
L
2
1
1 = 4
+
d 2 (...)
c3l
c3t
0
l
1/3
3
3
1
2 D
9N
L
+ 3
= 3N D =
=
3
L 1
c3l
ct
+ c23
c3
l
(23)
Define the Debye temperature D = }kBD , which defines the temperature scale for vibrational fluctuations. To test this model, we find the
heat capacity:
CV = kB
hEi
T3
= 9N kB 3
D
dx
0
x4 ex
D /T
(1 ex )2
10
Statistical Thermodynamics
To derive Plancks radiation law directly from our statistical mechanics,
we note that number of photons in the gas is not fixed. The photons
are absorbed and re-emitted by the enclosure walls. Since the photons
are non-interacting it is by this absorption and re-emission that equilibrium is maintained in the gas. Since, also the free energy F (T , V , N ) is
constant in equilibrium (at constant T and V ) while N varies it follows
that F /N = 0, that is
F
=
=0
N V ,N
(24)
The photon gas is then a Bose gas with = 0 so that the canonical
partition function is given directly as
Q=
r
Y
(1 exp(s ))1
(25)
s=1
(26)
(27)
8
h
2
3
c (exp(h ) 1)
(28)
u( ) =
or
u( ) =
hc/
' kT
(exp(hc/kT ) 1)
(29)
=
hc
exp(hc/kT )
(30)
(31)
1
8hc
()g () '
exp(hc/kT )
V
5
(32)
11
8k 4
(hc)3
x3
T 4 = aT 4
exp(x) 1
(34)
8k 4 4
(hc)3 15
(35)
dx
0
Where
a=
(36)
2ck 4 4
erg
= 5.67 105 2
3
(hc) 15
cm Sec K 4
(37)
E=
As
=
Thermodynamic Properties
We may calculate all the thermodynamic properties of black body radiation using statistical mechanics through the partition function Q, where
F = kT log(Q)
(38)
And Q was derived as shown in Eq. 25. This is the basic method
of statistical thermodynamics. The aim is to reproduce all the thermodynamic properties with all factors and constants evaluated. This
gives
12
r
Y
F = kT log
(1 exp(s ))
s=1
r
X
F = kT
s=1
d
log(1 exp(s ))1
h3
F = 2kT
(39)
Where the s are the single photon states and the factor of 2 arises
from the two polarizations available to each photon. This can be integrated in a variety of ways. Perhaps the most direct is to integrate
over phase space (d = dV 4p2 dp) and write = pc. Introducing the
dimensionless variable x = = pc, the Helmholtz free energy is:
F =
1
3
8k 4
(hc)3
V T4
(40)
The dimensionless integral here can be transformed into that appearing in the constant of a of Eq. 35, by an integration by parts, i.e.,
Z
I=
x3 d[log(1 exp(x))]
(41)
I=
dx
0
x3
=
exp(x) 1
15
(42)
F
T
= 4 aV T 3
3
V
(44)
(45)
13
F
V
= 1 aT 4
3
T
(46)
(47)
References
J. C. Ho, H. H. Hamdeh, M. W. Barsoum, and T. El-Raghy. Low
temperature heat capacity of T i3 SiC2 . Journal of Applied Physics, 85
(11):7970, June 1999.
14