International
Tables for Crystallography Volume B Reciprocal space Edited by U. Shmueli © International Union of Crystallography 2006 |
International Tables for Crystallography (2006). Vol. B. ch. 4.1, pp. 401-402
Section 4.1.2.1. Equations of motion
aChemical Crystallography Laboratory, University of Oxford, 9 Parks Road, Oxford OX1 3PD, England |
As a result of thermal fluctuations, the atoms vibrate about their equilibrium positions, so that the actual position of the th atom in the lth primitive cell is given by
with r representing the equilibrium position and u the thermal displacement. (In lattice-dynamical theory it is advantageous to deal with the primitive cell, as it possesses the fewest degrees of freedom.) The kinetic energy of the vibrating crystal is
where
is the mass of atom
and the index
(
= 1, 2, 3) refers to the Cartesian components of the displacement. (The dot denotes the time derivative.)
If the adiabatic approximation is invoked, the potential energy V of the crystal can be expressed as a function of the instantaneous atomic positions. Expanding V in powers of , using the three-dimensional form of Taylor's series, we have
where
is the static (equilibrium) potential and
,
are given by
The subscript zero indicates that the derivatives are to be evaluated at the equilibrium configuration. In the harmonic approximation,
and all higher terms in the expansion are neglected.
At equilibrium the forces on an atom must vanish, so that Ignoring the static potential
, the quadratic term
only remains and the Hamiltonian for the crystal (the sum of the kinetic and potential energies) is then
where
is an element of the 3 × 3 `atomic force-constant matrix' and is defined (for distinct atoms
,
) by
It is the negative of the force in the
direction imposed on the atom
when atom
is displaced unit distance along
with all the remaining atoms fixed at their equilibrium sites.
is defined differently for the self-term with
and
:
Thus the self-matrix describes the force on
when the atom itself is displaced with all the remaining atoms kept stationary.
There are restrictions on the number of distinct force constants : these are imposed by symmetry and by the requirement that the potential energy is invariant under infinitesimal translations and rotations of the rigid crystal. Such constraints are discussed in the book by Venkataraman et al. (1975
).
Applying Hamilton's equations of motion to equation (4.1.2.1) now gives
These represent 3nN coupled differential equations, where n is the number of atoms per primitive cell
and N is the number of cells per crystal
.
By applying the periodic boundary conditions, the solutions of equation (4.1.2.2) can be expressed as running, or travelling, plane waves extending throughout the entire crystal. The number of independent waves (or normal modes) is 3nN. Effectively, we have transferred to a new coordinate system: instead of specifying the motion of the individual atoms, we describe the thermal motion in terms of normal modes, each of which contributes to the displacement of each atom. The general solution for the
component of the displacement of
is then given by the superposition of the displacements from all modes:
Here q is the wavevector of a mode (specifying both its wavelength and direction of propagation in the crystal) and
its frequency. There are N distinct wavevectors, occupying a uniformly distributed mesh of N points in the Brillouin zone (reciprocal cell); each wavevector is shared by 3n modes which possess, in general, different frequencies and polarization properties. Thus an individual mode is conveniently labelled (j q), where j is an index (
) indicating the branch . The scalar quantity
in equation (4.1.2.3
) is the amplitude of excitation of (j q) and
is the element of the eigenvector
referring to the displacement in the
direction of the atom
. The eigenvector itself, with dimensions n × 1, determines the pattern of atomic displacements in the mode (j q) and its magnitude is fixed by the orthonormality and closure conditions
and
with
indicating complex conjugate and
the Kronecker delta.
The pre-exponential, or amplitude, terms in (4.1.2.3) are independent of the cell number. This follows from Bloch's (1928
) theorem which states that, for corresponding atoms in different cells, the motions are identical as regards their amplitude and direction and differ only in phase. The theorem introduces an enormous simplification as it allows us to restrict attention to the 3n equations of motion of the n atoms in just one cell, rather than the 3nN equations of motion for all the atoms in the crystal.
Substitution of (4.1.2.3) into (4.1.2.2
) gives the equations of motion in the form
in which
is an element of the dynamical matrix D(q).
is defined by
where
is the position of atom
with respect to the cell origin, L is
and
is the separation between cells l and
. The element
is obtained by writing down the
component of the force constant between atoms
which are L cells apart and multiplying by the phase factor
; this term is then summed over those values of L covering the range of interaction of
and
.
The dynamical matrix is Hermitian and has dimensions 3n × 3n. Its eigenvalues are the squared frequencies of the normal modes and its eigenvectors
determine the corresponding pattern of atomic displacements. The frequencies of the modes in three of the branches, j, go to zero as q approaches zero: these are the acoustic modes. The remaining
branches contain the optic modes. There are N distinct q vectors, and so, in all, there are 3N acoustic modes and
optic modes. Thus copper has acoustic modes but no optic modes, silicon and rock salt have an equal number of both, and lysozyme possesses predominantly optic modes.
References

