24.12.2014 Views

Set of supplementary notes.

Set of supplementary notes.

Set of supplementary notes.

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

52 CHAPTER 3. FROM ATOMS TO SOLIDS<br />

Figure 3.15: Comparison <strong>of</strong> Debye density <strong>of</strong> states (a) with that <strong>of</strong> a real material (b).<br />

There are 3 acoustic branches, and 3(m − 1) optical branches.<br />

It is convenient to start with a simple description <strong>of</strong> the optical branch(es), the Einstein<br />

model, which approximates the branch as having a completely flat dispersion ω(k) = ω 0 . In<br />

that case, the density <strong>of</strong> states in frequency is simply<br />

D E (ω) = Nδ(ω − ω 0 ) . (3.34)<br />

We have a different results for the acoustic modes, which disperse linearly with momentum<br />

as ω → 0. Using a dispersion ω = vk, and following the earlier argument used for electrons, we<br />

get the Debye model<br />

D D (ω) =<br />

4πk2 dk<br />

(2π/L) 3 dω = V ω2<br />

2π 2 v . (3.35)<br />

3<br />

Of course this result cannot apply once the dispersion curves towards, the zone boundary, and<br />

there must be an upper limit to the spectrum. In the Debye model, we cut <strong>of</strong>f the spectrum at a<br />

frequency ω D , which is determined so that the total number <strong>of</strong> states (N) is correctly counted,<br />

i.e. by choosing<br />

which yields<br />

∫ ωD<br />

0<br />

dωD D (ω) = N (3.36)<br />

ωD 3 = 6π2 v 3 N<br />

. (3.37)<br />

V<br />

Notice that this corresponds to replacing the correct cut<strong>of</strong>f in momentum space (determined<br />

by intersecting Brillouin zone planes) with a sphere <strong>of</strong> radius<br />

k D = ω D /v . (3.38)<br />

3.7 Lattice specific heat<br />

Phonons obey Bose-Einstein statistics, but their number is not conserved and so the chemical<br />

potential is zero, leading to the Planck distribution<br />

n(ω) =<br />

1<br />

exp(h¯ω/k B T ) − 1 . (3.39)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!