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. 400-402
Section 4.1.2. Dynamics of three-dimensional crystals
aChemical Crystallography Laboratory, University of Oxford, 9 Parks Road, Oxford OX1 3PD, England |
For modes of vibration of very long wavelength, the crystal can be treated as a homogeneous elastic continuum without referring to its crystal or molecular structure. The theory of the propagation of these elastic waves is based on Hooke's law of force and on Newton's equations of motion. As the wavelength of the vibrations becomes shorter and shorter and approaches the separation of adjacent atoms, the calculation of the vibrational properties requires a knowledge of the crystal structure and of the nature of the forces between adjacent atoms. The three-dimensional treatment is based on the formulation of Born and von Kármán, which is discussed in detail in the book by Born & Huang (1954) and in more elementary terms in the books by Cochran (1973
) and by Willis & Pryor (1975
).
Before setting up the equations of motion, it is necessary to introduce three approximations:
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.
Quantum concepts are not required in solving the equations of motion (4.1.2.4) to determine the frequencies and displacement patterns of the normal modes. The only place where quantum mechanics is necessary is in calculating the energy of the mode, and from this the amplitude of vibration
.
It is possible to discuss the theory of lattice dynamics from the beginning in the language of quantum mechanics (Donovan & Angress, 1971). Instead of treating the modes as running waves, they are conceived as an assemblage of indistinguishable quasi-particles called phonons. Phonons obey Bose–Einstein statistics and are not limited in number. The number of phonons, each with energy
in the vibrational state specified by q and j, is given by
and the mode energy
by
Thus the quantum number
describes the degree of excitation of the mode (j q). The relation between
and the amplitude
is
Equations (4.1.2.6
) to (4.1.2.8
) together determine the value of
to be substituted into equation (4.1.2.3
) to give the atomic displacement in terms of the absolute temperature and the properties of the normal modes.
In solving the lattice-dynamical problem using the Born–von Kármán analysis, the first step is to set up a force-constant matrix describing the interactions between all pairs of atoms. This is followed by the assembly of the dynamical matrix D, whose eigenvalues give the frequencies of the normal modes and whose eigenvectors determine the patterns of atomic displacement for each mode.
Before considering the extension of this treatment to molecular crystals, we shall comment briefly on the less rigorous treatments of Einstein and Debye.
In the Einstein model it is assumed that each atom vibrates in its private potential well, entirely unaffected by the motion of its neighbours. There is no correlation between the motion of different atoms, whereas correlated motion – in the form of collective modes propagating throughout the crystal – is a central feature in explaining the characteristics of the TDS. Nevertheless, the Einstein model is occasionally used to represent modes belonging to flat optic branches of the dispersion relations, with the frequency written symbolically as (constant).
In the Debye model the optic branches are ignored. The dispersion relations for the remaining three acoustic branches are assumed to be the same and represented by where
is a mean sound velocity. The Brillouin zone is replaced by a sphere with radius
chosen to ensure the correct number of modes. The linear relationship (4.1.2.9
) holds right up to the boundary of the spherical zone. In an improved version of the Debye model, (4.1.2.9
) is replaced by the expression
which is the same as (4.1.2.9
) at q = 0 but gives a sinusoidal dispersion relation with zero slope at the zone boundary.
The full Born–von Kármán treatment becomes excessively cumbersome when applied to most molecular crystals. For example, for naphthalene with two molecules or 36 atoms in the primitive cell, the dynamical matrix has dimensions 108 × 108. Moreover, the physical picture of molecules or of groups of atoms, vibrating in certain modes as quasi-rigid units, is lost in the full treatment.
To simplify the setting up of the dynamical matrix, it is assumed that the molecules vibrate as rigid units in the crystal with each molecule possessing three translational and three rotational (librational) degrees of freedom. The motion of these rigid groups as a whole is described by the external modes of motion, whereas the internal modes arise from distortions within an individual group. The frequencies of these internal modes, which are largely determined by the strong intramolecular forces, are unaffected by the phase of the oscillation between neighbouring cells: the modes are taken, therefore, to be equivalent to those of the free molecule. The remaining external modes are calculated by applying the Born–von Kármán procedure to the crystal treated as an assembly of rigid molecules.
The dynamical matrix D(q) now has dimensions , where
is the number of molecules in the primitive cell: for naphthalene, D is reduced to 12 × 12. The elements of D can be expressed in the same form as equation (4.1.2.5
) for an atomic system.
,
refer to molecules which are L cells apart and the indices
,
(
) label the six components of translation and rotation.
in equation (4.1.2.5
) is replaced by
where
represents the 3 × 3 molecular-mass matrix for
and the 3 × 3 moment-of-inertia matrix referred to the principal axes of inertia for
. The 6 × 6 force-torque constant matrices
are derived by taking the second derivative of the potential energy of the crystal with respect to the coordinates of translation and rotation.
References








