International
Tables for Crystallography Volume D Physical properties of crystals Edited by A. Authier © International Union of Crystallography 2013 |
International Tables for Crystallography (2013). Vol. D, ch. 1.9, pp. 235-245
Section 1.9.4. Graphical representation^{a}GZG Abt. Kristallographie, Goldschmidtstrasse 1, 37077 Göttingen, Germany |
Atomic displacement tensors (ADTs) described by their tensor coefficients may be represented graphically to clarify their physical meaning. Different graphical representations exist and will be discussed separately for second- and higher-order tensors in the following.
Numerous examples of graphical representations of thermal-motion tensors (or, more generally speaking, atomic displacement tensors) have appeared in the literature since the early days of the computer program ORTEP written by C. K. Johnson (1965), yet the equal-probability surface usually displayed is only one of the possible representations of a second-order atomic displacement tensor. Representation surfaces are usually calculated in a Cartesian coordinate system. Accordingly, one has to transform the second-order ADT b into U_{C} described in a Cartesian frame:The transformation matrix depends on the choice of Cartesian axes e_{i} with respect to the reciprocal-cell axes a^{i} (or equally well with respect to the direct axes a_{i}). Choosing e_{1} along a^{1}, e_{2} in the a^{1}a^{2} plane and e_{3} completing the right-handed set, one obtains for the transformation matrix F (see also Willis & Pryor, 1975)
Clearly, there is an infinite number of possible choices for relating a Cartesian frame to a crystallographic coordinate system with correspondingly different transformation matrices F (see e.g. Chapter 1.1 of IT B). The most useful representation surface of the second-order atomic displacement tensor U_{C} is the representation quadric defined by the tensor invariantwhere u is a displacement vector; U^{−1} is often called the variance–covariance matrix and has (in a general axes frame) covariant components. Under the conditions of positive definiteness,the surface of the representation quadric is an ellipsoid whose semi-major axes (for ) are of lengths equal to the root-mean-square displacements (r.m.s.d.'s) along the axes directions. The thermal vibration ellipsoids calculated in ORTEP are related to this surface; considering the discussion in Section 1.9.1, they should more appropriately be called atomic displacement ellipsoids or simply ORTEP ellipsoids. One notes that the Fourier transform of the atomic DWF, the atomic probability density function P(u), is given in the case of a second-order tensor as a trivariate Gaussian distribution,
On comparing (1.9.4.3) and (1.9.4.5), it is evident that (1.9.4.3) defines a surface of constant probability of finding a (displaced) atom. The integral of (1.9.4.5) over the volume inside the ellipsoid is a constant. For with the integration limit (2.5003), the integral is equal to one half (nine tenths), and the ellipsoid is then called a 50 (90) per cent probability ellipsoid.
Other representation surfaces can be defined and are useful for special considerations. The quantities of interest are either the r.m.s.d.'s or the mean-square displacements (m.s.d.'s) defined in direct space. Here a distinction has to be made between the averaged squared displacement along a certain direction and the average for all squared displacements of an atom projected onto a given direction. Representation surfaces may also be calculated in reciprocal space, related to surfaces in direct space by Fourier transformation. For further details, see Nelmes (1969) and Hummel et al. (1990).
Representation surfaces of higher-order tensors may be calculated from their invariants. While for second-order tensors surfaces can be found that fully describe the directional aspects of the tensor involved, higher-order tensors need several different surfaces for a full description (see e.g. Wondratschek, 1958; Sirotin, 1961). This makes the graphical representation of the displacements somewhat cumbersome and it is therefore rarely used. Instead, the probability density functions [given in equations (6.1.1.46 ), (6.1.1.48 ) or (6.1.1.49 ) of IT C] are calculated from the tensor coefficients and displayed in sections or as three-dimensional surfaces. If the higher-order terms are small, it is more appropriate to display only the difference between the total p.d.f. and the related Gaussian p.d.f., which may be calculated from the second-order displacement tensor using equation (1.9.4.5). Here, the second-order terms that were refined together with the higher-order terms are usually used (not the best-fitting second-order terms of a fit in the harmonic approximation):The resulting anharmonic deformation densities (or disorder deformation densities in the case of static disorder) P_{deformation}(u) may be displayed in a similar way to the total p.d.f.'s P_{general}(u). The graphical representations appropriate for displaying those densities are similar to those used for electronic deformation densities (see e.g. Smith et al., 1977). A number of examples of displacement deformation densities of high symmetry are shown in Fig. 1.9.4.1 as three-dimensional contour maps.
References
Hummel, W., Raselli, A. & Bürgi, H.-B. (1990). Analysis of atomic displacement parameters and molecular motion in crystals. Acta Cryst. B46, 683–692.Johnson, C. K. (1965). ORTEP – a FORTRAN thermal ellipsoid plot program. Report ORNL-3794. Oak Ridge National Laboratory, Tennessee, USA.
Nelmes, R. J. (1969). Representational surfaces for thermal motion. Acta Cryst. A25, 523–526.
Sirotin, Yu. I. (1961). Plotting tensors of a given symmetry. Sov. Phys. Crystallogr. 6, 263–271.
Smith, V. H. Jr, Price, P. F. & Absar, I. (1977). Representations of electron density and its topographical features. Isr. J. Chem. 16, 187–197.
Willis, B. T. M. & Pryor, A. W. (1975). Thermal vibrations in crystallography. Cambridge University Press.
Wondratschek, H. (1958). Über die Möglichkeit der Beschreibung kristallphysikalischer Eigenschaften durch Flächen. Z. Kristallogr. 110, 127–135.