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
Bloch, F. (1928). Über die Quantenmechanik der Elektronen in Kristallgittern. Z. Phys. 52, 555–600.Google ScholarVenkataraman, G., Feldkamp, L. A. & Sahni, V. C. (1975). Dynamics of perfect crystals. Cambridge Mass.: MIT Press.Google Scholar