Lec 4
Lec 4
Lec 4
H
N
m
({R
I
}) = [
I
-
h
2
2M
I
2
I
+E
0
({R
I
})]
m
({R
I
}) = E
m
m
({R
I
}) (1)
Dene displacements of the atoms from their equilibrium positions, R
I
= R
0
I
+u
I
, and
expand the potential energy in powers of u
I
,
E
0
({R
I
}) = E
0
({R
0
I
}) +
1
2
IJ
u
I
u
J
+. . . (2)
where the linear terms vanish, the Ds are second derivatives, etc. The exact eigenstates
can be written as
m
({u
I
}).
Harmonic Approximation
In the harmonic approximation, where we neglect all higher order derivatives, the equations
are exactly soluble as independent oscillators, even though it is a strongly coupled system
of nuclei. This is a very good approximation at low temperatures for all solids except H
and He, where zero point motion is large. (Note that the Lindeman criterion for melting is
_
u
2
0.2a, where a is a typical interatomic distance, so the harmonic approximation is
often good for many purposes up to the melting temperature.
The quantum system can be treated by solving the classical equations of motion for
the classical oscillator eigenmodes and quantizing the oscillators. From the Bloch theo-
rem the eigenmodes can be classied by the k, e.g., for one atom per cell, (considering a
normalization to a volume of N unit cells)
u
I
=
1
k
u
k
e
ikR
0
I
; p
I
= M
du
I
dt
=
1
k
p
k
e
ikR
0
I
, (3)
The classical hamiltonian becomes
H =
k
_
1
2M
p
k
p
k
+
M
2
2
k
u
k
u
k
_
. (4)
where
M
2
k
J
D
IJ
e
ik(R
0
I
R
0
J
)
. (5)
Dening new operators
Q
k
=
_
M
k
2
-
h
_
1/2
u
k
, P
k
=
_
1
2
-
hM
k
_
1/2
u
k
, (6)
2 561 F 2005 Lecture 4
the hamiltonian becomes
H =
k
-
h
k
[P
k
P
k
+Q
k
Q
k
] . (7)
Quantizing the independent oscillators,the creation and annihilation operators are
a
k
= (Q
k
iP
k
), a
k
= (Q
k
+iP
k
), (8)
and the hamiltonian becomes
H =
1
2
k
-
h
k
_
a
k
a
k
+a
k
a
k
_
=
k
-
h
k
_
a
k
a
k
+
1
2
_
. (9)
Note that the order of the operators is crucial!
The state of the nuclei in the harmonic approximation is fully specied by the number
of quanta in each of the eigenmodes modes of vibration. Since the quantized modes are
bosons, the many body state is the symmetrized combination, which is simply a product
wavefunction.
It is straightforward to show that in d dimensions there are d modes that have frequency
that vanishes as k 0.
This is the rst example of elementary excitations. This gives the exact many-body
description in the harmonic approximation, and it denes a complete set of eigenstates
which are a convenient basis for describing the system even when anharmonic interactions
are included.
Collective Excitation
A phonon at a particular k, is a collective excitation in which all the nuclei move
in unison. The harmonic case is the special case where the equations for this collective
excitation are easily solved to nd non-interacting excitations (phonons). In general this is
not the case!
Qualitative examples of phonons
See Phillips Ch. 10 for the example of a chain of atoms.
Since the nuclei are charged, long ranged Coulomb interactions lead to longitudinal
plasma modes with that is nite as k 0 unless the electrons screen the long range
interactions. This happens in a metal, and a homework problem is to derive the qualitative
result that the phonons properly have accoustic modes with k, and qualitatively relate
the sound velocity to the Fermi velocity.
Anharmonic Interactions
In the presence of anharmonic interactions, the problem is no longer exactly soluble.
The energy can always be written in an expansion
E
0
({R
I
}) = E
0
({R
0
I
}) +
1
2!
IJ
D
IJ
u
I
u
J
+
1
3!
IJK
D
IJK
u
I
u
J
u
K
+. . . , (10)
561 F 2005 Lecture 4 3
where we have dropped the vector index for simplicity. It is most advantageous to work in
the basis of the elementary excitations, (the phonons) since the harmonic terms are diagonal
in this basis. Since the u
H =
p
p
c
+
p
c
p
+
1
2
k
V
k
pp
c
+
p+k
c
+
p
c
p
c
p
(11)
=
p
p
c
+
p
c
p
+
1
2
k
V
k
(
k
k
N), (12)
where V
k
=
1
4e
2
k
2
and
k
=
p
c
+
p+k
c
p
. In general one always has sums of the form
1
k
1
(2)
D
_
d
D
k, where D is the dimension. (We will sometimes set the normalization
volume = 1.)
This hamiltonian is valid only for homogeneous Jellium because we have assumed the
single particle terms are diagonal in momentum index p. However, a general form can be
made for for atoms, crystals, and molecules if one changes only the single particle part,
since the electron-electron interactions always have the Coulombic form shown.
For example, the exchange energy can be derived as
E
x
=
1
2
k
V
k
pp
0|c
+
p+k
c
+
p
c
p
c
p
|0 (13)
=
1
N
(
1
2
k,k
,
V
k
)
k+k
k,
(14)
=
1
N
(
1
2
k,k
<k
F
V
kk
) =
0.916
r
s
Ryd. (15)
Hartree-Fock Theory for Jellium
We have previously derived the form of exchange in the rst-quantized notation. The
key ingredients were the pair distribution function dened by
g(r, r
) =
1
(r)
1
(r
)
p(r, r
) (16)
where
p(r, r
) = |
i=j
(r r
i
)(r r
j
)| (17)
4 561 F 2005 Lecture 4
is the joint probability of simultaneously nding an electron at points r and r
.
In second quantized form in terms of the elementary excitations (plane waves) it is most
convenient to work in k-space. For a homogeneous system one can dene the static structure
factor S(k) (see Pines and Fulde)
S(k) =
_
d
3
rp(r)exp(ikr) (18)
=
1
N
|
i,j
exp(ik(r
i
r
j
)| (19)
=
1
N
|
k
k
| (20)
which includes the self-term. In terms of creation and annihilation operators, this may be
written
S(k) =
1
N
0|
p
c
pk
c
p
+k
c
p
|0, (21)
which after some algebra can be shown to be
S(k) = N
k=0
+
1
N
p
(1 n
pk
)n
p
, (22)
where is the occupation n
p
= 1 for k < k
F
and 0 for k > k
F
. (See Pines p 75 and class
notes.)
Thus the Hartree-Fock structure factor is very simple in k space (See Pines p 75 and
class notes.) For k = 0, dene x =
k
2k
F
, in which case
S(k) =
3
2
x
1
2
x
3
, x < 1 (23)
S(k) = 1, x > 1 (24)
A main point is that S(k) = 1 for k > 2k
F
, so that the the particles are uncorrelated except
for the Pauli exclusion between states with k < k
F
.
For a homogeneous system, one can Fourier transform this expression
g(r) =
1
N
k
(S(k) 1)exp(ikr) (25)
to recover the same results as derived before for the exchange hole in real space. (Note: For
an uncorrelated system S(k) = 1, k = 0, and S(0) = N, (self term). Some texts dene S
without the self term, so that S(0) = 0, e.g., Mahan, p. 74.)
Beyond Hartree-Fock: Correlations
In the class notes we give the main reasons why it is dicult to go beyond H-F. In nite
systems, one nds the conguration-interaction expansion which grows exponentially with
the number of particles. In condensed matter expansions formally diverge (Gellmann-
Bruckner). This is easily seen in Jellium just by examining the integrals. What to do?
This is the subject of the many-body theory of the course.
Conclusions
561 F 2005 Lecture 4 5
Elementary excitations: Phonons and electrons
Second quantization convenient for specifying hamiltonian and interactions
Creation and annihilation operators dened to create and annihilate eigenstates of a
single particle hamiltonian: harmonic phonons; non-interacting electrons
Nuclear motion is a strongly interacting problem - the motion of nuclei is highly
correlated - but in the harmonic approximation one can transform to non-interacting
elementary excitations - phonons - a phonon is a collective excitation of the nuclei
Interacting electron hamiltonian
Structure factor and Exchange hole in Hartree-Fock approximation revisited in second
quantization
How to go beyond Hartree-Fock? Perturbation expansion in the interaction ( e
2
)
formally diverges - can be summed in a convergent way - leads to non-analytic behavior
(as it must!)
Gell-Mann Brueckner results; Wigner interpolation formula; improved results