**Lattice dynamics- starting assumptions**

- Adiabatic approximation- assume electrons are attached rigidly to the nucleus.
- Assume the amplitude of the vibrations is small.
- Harmonic approximations- terms of order higher than
*u*in the interatomic potential are neglected.^{2}

**One-dimensional monatomic chain:** we’ll start by only including nearest neighbour forces.

positions *r _{n}*, equilibrium atomic separation

*r*

_{n}-r_{n-1}=adisplacements from equilibrium *u _{n}*

Interatomic pair potential

Potential energy of *nth* atom

Taylor expansion about *r _{n}-r_{n-1}=a*

*Φ’(a)=*0, so

Force

Also note that *Φ’’(a) *can be written as the force constant, *C*.

In fact, having done all that, we could have just started from this point by remembering Hooke’s Law *F=Cx* (where *x *is displacement from equilibrium).

Aside: just to check that we aren’t fiddling our force constants to make everything neat, we can use Newton III to show that *C _{n}=C_{-n}*

If we wanted to include more distant neighbour interactions (up to the Nth), our equation at this point would have the form

Going back to nearest neighbours, try a solution of the form

->

**Dispersion relation**

Nearest neighbours only

In the long wavelength limit (small *k*),

speed of sound

Considering more distant neighbours

Rule of thumb- number of modulations = number of further neighbour interactions considered.

**Brillouin Zones:** all the physics of a system is contained within a Brillouin Zone. The standard definition of the BZ is that it is “the Wigner-Seitz unit cell in reciprocal space”, but there’s no point just quoting this if you don’t understand (or can’t readily explain) what it means. For that reason, it also helps to remember this alternate definition-

The first Brillouin Zone is the set of points closer to the origin in reciprocal space than any other reciprocal lattice point.

The *n*th Brillouin Zone is the set of points reached from the origin in reciprocal space by crossing *n-1* Bragg planes and no fewer.

Each Brillouin Zone has the same volume. Any normal mode can be represented by a wavevector in the first Brillouin Zone.

**1D diatomic chain- intuitive method**

For a diatomic chain, the repeat distance is doubled in real space, so it is halved in *k-*space (i.e. the size of the Brillouin Zone is halved). At first glance, this seems to mean we have lost the information contained in the region. If we want to retain it, we will have to translate it back into the new 1^{st} Brillouin Zone.

If the two types of atom have different masses, a gap opens up between the two branches.

As you’ll see from the graphs, we’ve suddenly introduced the terms acoustic and optic branch, but don’t panic. You didn’t miss anything- these terms will be described below once we go through the maths of the 1D diatomic chain.

**1D diatomic chain- the maths:** use all the same assumptions as the monatomic chain.

Equation of motion for white atoms

Equation of motion for black atoms

Try solutions *m _{1}*

*m _{2}*

->

Eliminating *α* leads to

An alternate (and quite neat) method to get this answer without too much trouble is to use matrices as follows.

Start with trial solutions of the form

This gives

We can rewrite this in matrix form

For non-trivial solutions

This gives us a quadratic in *ω ^{2}* which we can solve easily enough to get the answer given above.

**Acoustic and optic modes:** first, we’ll consider what is going on at points A to D. Throughout this section, *m _{1}* >

*m*.

_{2}A. *ka<<1,* ,

*m _{1}* and

*m*oscillate in antiphase and with different amplitudes.

_{2}B. , ,

Heavier atoms are stationary, lighter atoms oscillate.

C. , ,

Lighter atoms are stationary, heavier atoms oscillate.

D. *ka<<1*, *ω~v _{s}k,* $latex \frac{B}{A}\simeq 1

Both masses oscillate in phases with the same amplitude.

On the acoustic branch, atoms are in phase near D -> sound waves.

On the optic branch, atoms are in antiphase near A. There is an overall dipole moment so optical transitions can be excited.

**In 3D:** for *p* atoms in the basis, we have 3*p* equations of motion in 3D -> 3*p* branches of the dispersion curve. Three of these branches are acoustic modes (1 longitudinal, 2 transverse), leaving 3*p*-3 optic branches.

Since it is harder to shear than to compress, *v _{s}* is smaller for the transverse modes.

**Lattice quantisation and normal modes:** if a system is in a normal mode, all components of the system are oscillating at the same frequency. If the Hamilton of the system contains only harmonic terms, these modes are uncoupled.

Take a system with a purely harmonic potential so that the normal modes look like harmonic oscillators.

Energy levels

If we want to go ahead and work out these summations, then proceed as follows.

(geometric progression)

as in thermodynamics.

**Lattice heat capacity:**

For phonons,

**Einstein model:** assume each atom is in an isolated potential at frequency *ω _{E}*. All atoms are at the same frequency, and each atom vibrates independently (no coupling).

(in 3D we have 3*N* oscillators)

Define the **Einstein temperature**

The Einstein model is not good for low temperature, although it does work fairly well at high temperature, where it tends to the classical result. Overall, though, this is not a good model- Einstein didn’t get it right all the time!

**Debye model: **approximate the dispersion curve as *ω=v _{s}k* all the way up to

*k=±π/a*.

3N oscillators (still working in 3D)

for *ω≤ ω _{D}*

otherwise

*ω _{D}* cut-off determined by (1)

Heat capacity

Put and define the **Debye temperature**

Now we can write

where

Substitute for *v _{s}* using (1) and tidy up to get

At low temperatures, *x* is large and upper states are unoccupied, so we can let the upper limit tend to infinity.

The integral is a standard one which we can look up to see that the answer is

At high temperatures *x* is small, so we can expand *e ^{x}* and ignore higher order terms.

(Dulong-Petit law)

The Debye model is a good approximation at low () and high temperatures, but it is not very good for intermediate temperatures.

**Actual density of states:** as we know, the real dispersion relation has a lower *ω _{max}* than the one given by the Debye model- this means that

*g(ω)*peaks earlier than in the Debye model.

Below this peak, as Debye says.

Considering transverse modes-

We would also see higher peaks if optic modes are present.

**Density of states summary**

Dimensions |
g(k) |
g(ω) |
Heat capacity |
||

electrons |
photons/phonons / |
low temp due to phonons |
due to electrons |
||

1 |
const |
~ε |
const |
~T |
~T |

2 |
~k |
const |
~ε |
~T |
~T |

3 |
~k |
~ε |
~ε |
~T |
~T |