Tables for
Volume B
Reciprocal space
Edited by U. Shmueli

International Tables for Crystallography (2006). Vol. B, ch. 1.2, pp. 10-24   | 1 | 2 |

Chapter 1.2. The structure factor

P. Coppensa*

aDepartment of Chemistry, Natural Sciences & Mathematics Complex, State University of New York at Buffalo, Buffalo, New York 14260-3000, USA
Correspondence e-mail:

This chapter summarizes the mathematical development of the structure-factor formalism. It starts with the definition of the structure factor appropriate for X-ray and neutron scattering, and includes the derivation of the appropriate expressions. Beyond the isolated-atom case, the atom-centred spherical harmonic (multipole) model is treated in detail. Expressions for thermal motion in the harmonic approximation and for treatments including anharmonicity are given and their relative merits are discussed.

1.2.1. Introduction

| top | pdf |

The structure factor is the central concept in structure analysis by diffraction methods. Its modulus is called the structure amplitude. The structure amplitude is a function of the indices of the set of scattering planes h, k and l, and is defined as the amplitude of scattering by the contents of the crystallographic unit cell, expressed in units of scattering. For X-ray scattering, that unit is the scattering by a single electron [(2.82 \times 10^{-15}\;\hbox{m})], while for neutron scattering by atomic nuclei, the unit of scattering length of [10^{-14}\;\hbox{m}] is commonly used. The complex form of the structure factor means that the phase of the scattered wave is not simply related to that of the incident wave. However, the observable, which is the scattered intensity, must be real. It is proportional to the square of the scattering amplitude (see, e.g., Lipson & Cochran, 1966[link]).

The structure factor is directly related to the distribution of scattering matter in the unit cell which, in the X-ray case, is the electron distribution, time-averaged over the vibrational modes of the solid.

In this chapter we will discuss structure-factor expressions for X-ray and neutron scattering, and, in particular, the modelling that is required to obtain an analytical description in terms of the features of the electron distribution and the vibrational displacement parameters of individual atoms. We concentrate on the most basic developments; for further details the reader is referred to the cited literature.

1.2.2. General scattering expression for X-rays

| top | pdf |

The total scattering of X-rays contains both elastic and inelastic components. Within the first-order Born approximation (Born, 1926[link]) it has been treated by several authors (e.g. Waller & Hartree, 1929[link]; Feil, 1977[link]) and is given by the expression [I_{\rm total} ({\bf S}) = I_{\rm classical} {\textstyle\sum\limits_{n}} \left|{\textstyle\int} \psi_{n}^{*} \exp (2\pi i{\bf S}\cdot {\bf r}_{j}) \psi_{0}\; \hbox{d}{\bf r}\right|^{2}, \eqno(] where [I_{\rm classical}] is the classical Thomson scattering of an X-ray beam by a free electron, which is equal to [(e^{2}/mc^{2})^{2} (1 + \cos^{2} 2\theta)/2] for an unpolarized beam of unit intensity, ψ is the n-electron space-wavefunction expressed in the 3n coordinates of the electrons located at [{\bf r}_{j}] and the integration is over the coordinates of all electrons. S is the scattering vector of length [2\sin \theta/\lambda].

The coherent elastic component of the scattering, in units of the scattering of a free electron, is given by [I_{\rm coherent, \, elastic} ({\bf S}) = \left|{\textstyle\int} \psi_{0}^{*}\right| {\textstyle\sum\limits_{j}} \exp (2\pi i{\bf S}\cdot {\bf r}_{j}) |\psi_{0} \;\hbox{d}{\bf r}|^{2}. \eqno(]

If integration is performed over all coordinates but those of the jth electron, one obtains after summation over all electrons [I_{\rm coherent, \, elastic} ({\bf S}) = |{\textstyle\int} \rho ({\bf r}) \exp (2\pi i{\bf S}\cdot {\bf r}) \;\hbox{d}{\bf r}|^{2}, \eqno(] where [\rho({\bf r})] is the electron distribution. The scattering amplitude [A({\bf S})] is then given by [A({\bf S}) = {\textstyle\int} \rho ({\bf r}) \exp (2\pi i{\bf S}\cdot {\bf r}) \;\hbox{d}{\bf r} \eqno(] or [A({\bf S}) = \hat{F} \{\rho ({\bf r})\}, \eqno(] where [\hat{F}] is the Fourier transform operator.

1.2.3. Scattering by a crystal: definition of a structure factor

| top | pdf |

In a crystal of infinite size, [\rho ({\bf r})] is a three-dimensional periodic function, as expressed by the convolution [\rho_{\rm crystal} ({\bf r}) = {\textstyle\sum\limits_{n}} {\textstyle\sum\limits_{m}} {\textstyle\sum\limits_{p}} \rho_{\rm unit\; cell} ({\bf r})\ast \delta ({\bf r} - n{\bf a} - m{\bf b} - p{\bf c}), \eqno(] where n, m and p are integers, and δ is the Dirac delta function.

Thus, according to the Fourier convolution theorem, [A({\bf S}) = \hat{F} \{\rho ({\bf r})\} = {\textstyle\sum\limits_{n}} {\textstyle\sum\limits_{m}} {\textstyle\sum\limits_{p}} \hat{F} \{\rho_{\rm unit\; cell} ({\bf r})\} \hat{F} \{\delta ({\bf r} - n{\bf a} - m{\bf b} - p{\bf c}\}, \eqno(] which gives [A({\bf S}) = \hat{F} \{\rho_{\rm unit\; cell} ({\bf r})\} {\textstyle\sum\limits_{h}} {\textstyle\sum\limits_{k}} {\textstyle\sum\limits_{l}} \delta ({\bf S} - h{\bf a}^{*} - k{\bf b}^{*} - l{\bf c}^{*}). \eqno(]

Expression ([link] is valid for a crystal with a very large number of unit cells, in which particle-size broadening is negligible. Furthermore, it does not account for multiple scattering of the beam within the crystal. Because of the appearance of the delta function, ([link] implies that [{\bf S} = {\bf H}] with [{\bf H} = h{\bf a}^{*} + k{\bf b}^{*} + l{\bf c}^{*}].

The first factor in ([link], the scattering amplitude of one unit cell, is defined as the structure factor F: [F({\bf H}) = \hat{F} \{\rho_{\rm unit\; cell} ({\bf r})\} = {\textstyle\int_{\rm unit\; cell}} \rho ({\bf r}) \exp (2\pi i{\bf H}\cdot {\bf r}) \hbox{ d}{\bf r}. \eqno(]

1.2.4. The isolated-atom approximation in X-ray diffraction

| top | pdf |

To a reasonable approximation, the unit-cell density can be described as a superposition of isolated, spherical atoms located at [{\bf r}_{j}]. [\rho_{\rm unit\; cell} ({\bf r}) = {\textstyle\sum\limits_{j}} \rho_{{\rm atom}, \, j} ({\bf r})\ast\delta ({\bf r} - {\bf r}_{j}). \eqno(] Substitution in ([link] gives [{F({\bf H}) = {\textstyle\sum\limits_{j}} \hat{F} \{\rho_{{\rm atom}, \, j}\} \hat{F} \{\delta ({\bf r} - {\bf r}_{j})\} = {\textstyle\sum\limits_{j}}\; f_{j} \exp (2\pi i{\bf H}\cdot {\bf r}_{j})} \eqno(] or [\eqalignno{F(h, k, l) &= {\textstyle\sum\limits_{j}}\; f_{j} \exp 2\pi i(hx_{j} + ky_{j} + lz_{j})\cr &= {\textstyle\sum\limits_{j}}\; f_{j} \{\cos 2\pi (hx_{j} + ky_{j} + lz_{j})\cr &\quad + i \sin 2\pi (hx_{j} + ky_{j} + lz_{j})\}. &(}] [f_{j}(S)], the spherical atomic scattering factor, or form factor, is the Fourier transform of the spherically averaged atomic density [\rho_{j}(r)], in which the polar coordinate r is relative to the nuclear position. [f_{j}(S)] can be written as (James, 1982[link]) [\eqalignno{f_{j}(S) &= {\int\limits_{\rm atom}} \rho_{j} (r) \exp (2\pi i{\bf S}\cdot {\bf r}) \;\hbox{d}{\bf r}\cr &= {\int\limits_{\upsilon = 0}^{\pi}}\; {\int\limits_{\varphi = 0}^{2\pi}} \;{\int\limits_{r = 0}^{\infty}} \rho_{j} (r) \exp (2\pi i Sr \cos \vartheta) r^{2} \sin \vartheta\; \hbox{d}r \;\hbox{d}\vartheta \;\hbox{d}\varphi\cr &= {\int\limits_{0}^{r}} 4\pi r^{2} \rho_{j} (r) {\sin 2\pi Sr \over 2\pi Sr}\; \hbox{d}r \equiv {\int\limits_{0}^{r}} 4\pi r^{2} \rho_{ j} (r) j_{0} (2\pi Sr)\; \hbox{d}r\cr & \equiv \langle \;j_{0}\rangle, &(}] where [j_{0} (2\pi Sr)] is the zero-order spherical Bessel function.

[\rho_{j}(r)] represents either the static or the dynamic density of atom j. In the former case, the effect of thermal motion, treated in Section 1.2.9[link] and following, is not included in the expression.

When scattering is treated in the second-order Born approximation, additional terms occur which are in particular of importance for X-ray wavelengths with energies close to absorption edges of atoms, where the participation of free and bound excited states in the scattering process becomes very important, leading to resonance scattering. Inclusion of such contributions leads to two extra terms, which are both wavelength- and scattering-angle-dependent: [f_{j} (S, \lambda) = {f_{j}}^{0} (S) + f'_{j} (S, \lambda) + if''_{j} (S, \lambda). \eqno(]

The treatment of resonance effects is beyond the scope of this chapter. We note however (a) that to a reasonable approximation the S-dependence of j′ and j″ can be neglected, (b) that j′ and j″ are not independent, but related through the Kramers–Kronig transformation, and (c) that in an anisotropic environment the atomic scattering factor becomes anisotropic, and accordingly is described as a tensor property. Detailed descriptions and appropriate references can be found in Materlick et al. (1994)[link] and in Section 4.2.6[link] of IT C (2004)[link].

The structure-factor expressions ([link] [link] [link] [link] can be simplified when the crystal class contains non-trivial symmetry elements. For example, when the origin of the unit cell coincides with a centre of symmetry [(x, y, z \rightarrow -x, -y, -z)] the sine term in ([link] cancels when the contributions from the symmetry-related atoms are added, leading to the expression [F({\bf H}) = 2 {\textstyle\sum\limits_{j = 1}^{N/2}}\; f_{j} \cos 2\pi (hx_{j} + ky_{j} + lz_{j}), \eqno(] where the summation is over the unique half of the unit cell only.

Further simplifications occur when other symmetry elements are present. They are treated in Chapter 1.4[link] , which also contains a complete list of symmetry-specific structure-factor expressions valid in the spherical-atom isotropic-temperature-factor approximation.

1.2.5. Scattering of thermal neutrons

| top | pdf | Nuclear scattering

| top | pdf |

The scattering of neutrons by atomic nuclei is described by the atomic scattering length b, related to the total cross section [\sigma_{\rm total}] by the expression [\sigma_{\rm total} = 4\pi b^{2}]. At present, there is no theory of nuclear forces which allows calculation of the scattering length, so that experimental values are to be used. Two types of nuclei can be distinguished (Squires, 1978[link]). In the first type, the scattering is a resonance phenomenon and is associated with the formation of a compound nucleus (consisting of the original nucleus plus a neutron) with an energy close to that of an excited state. In the second type, the compound nucleus is not near an excited state and the scattering length is essentially real and independent of the energy of the incoming neutron. In either case, b is independent of the Bragg angle θ, unlike the X-ray form factor, since the nuclear dimensions are very small relative to the wavelength of thermal neutrons.

The scattering length is not the same for different isotopes of an element. A random distribution of isotopes over the sites occupied by that element leads to an incoherent contribution, such that effectively [\sigma_{\rm total} = \sigma_{\rm coherent} + \sigma_{\rm incoherent}]. Similarly for nuclei with non-zero spin, a spin incoherent scattering occurs as the spin states are, in general, randomly distributed over the sites of the nuclei.

For free or loosely bound nuclei, the scattering length is modified by [b_{\rm free} = [M/(m + M)]b], where M is the mass of the nucleus and m is the mass of the neutron. This effect is of consequence only for the lightest elements. It can, in particular, be of significance for hydrogen atoms. With this in mind, the structure-factor expression for elastic scattering can be written as [F({\bf H}) = {\textstyle\sum\limits_{j}} b_{j, \, {\rm coherent}} \exp 2\pi i(hx_{j} + ky_{j} + lz_{j}) \eqno(] by analogy to ([link]. Magnetic scattering

| top | pdf |

The interaction between the magnetic moments of the neutron and the unpaired electrons in solids leads to magnetic scattering. The total elastic scattering including both the nuclear and magnetic contributions is given by [|F({\bf H}){|^{2}_{\rm total}} = |F_{N} ({\bf H}) + {\bf Q}({\bf H})\cdot {\hat\boldlambda}|^{2}, \eqno(] where the unit vector [\hat{{\boldlambda}}] describes the polarization vector for the neutron spin, [F_{N}({\bf H})] is given by ([link] and Q is defined by [{\bf Q} = {mc \over eh} {\int} {\bf \widehat{H}} \times [{\bf M}({\bf r}) \times {\bf \widehat{H}}] \exp (2\pi i{\bf H}\cdot {\bf r})\ \hbox{d}{\bf r}. \eqno(] [{\bf M}({\bf r})] is the vector field describing the electron-magnetization distribution and [\widehat{{\bf H}}] is a unit vector parallel to H.

Q is thus proportional to the projection of M onto a direction orthogonal to H in the plane containing M and H. The magnitude of this projection depends on [\sin \alpha], where α is the angle between Q and H, which prevents magnetic scattering from being a truly three-dimensional probe. If all moments [{\bf M}({\bf r})] are collinear, as may be achieved in paramagnetic materials by applying an external field, and for the maximum signal (H orthogonal to M), ([link] becomes [{\bf Q} = {\bf M}({\bf H}) = {mc \over eh} {\int} {\bf M}({\bf r}) \exp (2\pi i{\bf H}\cdot {\bf r})\ \hbox{d}{\bf r} \eqno(] and ([link] gives [{|F|^{2}}_{\!\!\!\rm total} = |F_{N}({\bf H}) - M({\bf H})|^{2} \eqno(] and [|F{|^{2}_{\rm total}} = |F_{N} ({\bf H}) + M({\bf H})|^{2}] for neutrons parallel and antiparallel to [{\bf M} ({\bf H})], respectively.

1.2.6. Effect of bonding on the atomic electron density within the spherical-atom approximation: the kappa formalism

| top | pdf |

A first improvement beyond the isolated-atom formalism is to allow for changes in the radial dependence of the atomic electron distribution.

Such changes may be due to electronegativity differences which lead to the transfer of electrons between the valence shells of different atoms. The electron transfer introduces a change in the screening of the nuclear charge by the electrons and therefore affects the radial dependence of the atomic electron distribution (Coulson, 1961[link]). A change in radial dependence of the density may also occur in a purely covalent bond, as, for example, in the H2 molecule (Ruedenberg, 1962[link]). It can be expressed as [{\rho'_{\rm valence}} (r) = \kappa^{3} \rho_{\rm valence} (\kappa r) \eqno(] (Coppens et al., 1979[link]), where ρ′ is the modified density and κ is an expansion/contraction parameter, which is > 1 for valence-shell contraction and < 1 for expansion. The [\kappa^{3}] factor results from the normalization requirement.

The valence density is usually defined as the outer electron shell from which charge transfer occurs. The inner or core electrons are much less affected by the change in occupancy of the outer shell and, in a reasonable approximation, retain their radial dependence.

The corresponding structure-factor expression is [\eqalignno{F({\bf H}) &= {\textstyle\sum\limits_{j}} [\{P_{j, \, {\rm core}}\;f_{j, \, {\rm core}} (H) + P_{j, \, {\rm valence}}\;f_{j, \, {\rm valence}} (H/\kappa)\}\cr &\quad \times \exp (2\pi i{\bf H}\cdot {\bf r}_j)], &(}] where [P_{j, \, {\rm core}}] and [P_{j, \, {\rm valence}}] are the number of electrons (not necessarily integral) in the core and valence shell, respectively, and the atomic scattering factors [f_{j, \, {\rm core}}] and [f_{j, \, {\rm valence}}] are normalized to one electron. Here and in the following sections, the anomalous-scattering contributions are incorporated in the core scattering.

1.2.7. Beyond the spherical-atom description: the atom-centred spherical harmonic expansion

| top | pdf | Direct-space description of aspherical atoms

| top | pdf |

Even though the spherical-atom approximation is often adequate, atoms in a crystal are in a non-spherical environment; therefore, an accurate description of the atomic electron density requires non-spherical density functions. In general, such density functions can be written in terms of the three polar coordinates r, θ and ϕ. Under the assumption that the radial and angular parts can be separated, one obtains for the density function: [\Phi (r,\theta,\varphi) = R(r)\Theta (\theta,\varphi). \eqno(]

The angular functions Θ are based on the spherical harmonic functions [Y_{lm}] defined by [Y_{lm} (\theta, \varphi) = (-1)^{m} \left[\left({2l + 1 \over 4\pi}\right) {(l - |m|)! \over (l + |m|)!}\right]^{1/2} P_{l}^{m} (\cos \theta) \exp (im \varphi), \eqno(] with [-l \leq m \leq l], where [P_{l}^{m} (\cos \theta)] are the associated Legendre polynomials (see Arfken, 1970[link]). [\eqalign{P_{l}^{m} (x) &= (1 - x^{2})^{|m|/2} {\hbox{d}^{|m|}P_{l}(x) \over \hbox{d}x^{|m|}},\cr P_{l} (x) &= {1 \over l!2^{l}} {\hbox{d}^{l} \over \hbox{d}x^{l}} \left[(x^{2} - 1)^{l}\right].}]

The real spherical harmonic functions [y_{lmp}], [0 \leq m \leq l], [p = + \hbox{ or } -] are obtained as a linear combination of [Y_{lm}]: [\eqalignno{y_{lm+} (\theta, \psi) &= \left[{(2l + 1)(l - |m|)! \over 2\pi (1 + \delta_{m0}) (l + | m |)!}\right]^{1/2} P_{l}^{m} (\cos \theta) \cos m\varphi\cr &= N_{lm} P_{l}^{m} (\cos \theta)\cos m\varphi\cr & = (-1)^{m} (Y_{lm} + Y_{l, \, -m}) &(}] and [\eqalignno{y_{lm-} (\theta, \psi) &= N_{lm} P_{l}^{m} (\cos \theta)\sin m\varphi\cr &= (-1)^{m} (Y_{lm} - Y_{l, \, -m})/2i.& (\cr}] The normalization constants [N_{lm}] are defined by the conditions [{\textstyle\int} y_{lmp}^{2} \hbox{ d}\Omega = 1, \eqno(] which are appropriate for normalization of wavefunctions. An alternative definition is used for charge-density basis functions: [{\textstyle\int} |d_{lmp}| \hbox{ d}\Omega = 2 \hbox{ for } l \;\gt\; 0 \hbox{ and } \textstyle\int |d_{lmp}| \hbox{ d}\Omega = 1 \hbox{ for } l = 0. \eqno(] The functions [y_{lmp}] and [d_{lmp}] differ only in the normalization constants. For the spherically symmetric function [d_{00}], a population parameter equal to one corresponds to the function being populated by one electron. For the non-spherical functions with [l \;\gt\; 0], a population parameter equal to one implies that one electron has shifted from the negative to the positive lobes of the function.

The functions [y_{lmp}] and [d_{lmp}] can be expressed in Cartesian coordinates, such that [y_{lmp} = M_{lm} c_{lmp} \eqno(] and [d_{lmp} = L_{lm} c_{lmp}, \eqno(] where the [c_{lmp}] are Cartesian functions. The relations between the various definitions of the real spherical harmonic functions are summarized by [Scheme scheme1] in which the direction of the arrows and the corresponding conversion factors [X_{lm}] define expressions of the type ([link] [link]. The expressions for [c_{lmp}] with [l \leq 4] are listed in Table[link], together with the normalization factors [M_{lm}] and [L_{lm}].

Table| top | pdf |
Real spherical harmonic functions (x, y, z are direction cosines)

lSymbol CAngular function, [c_{lmp}]Normalization for wavefunctions, [M_{lmp}]§Normalization for density functions, [L_{lmp}]
ExpressionNumerical valueExpressionNumerical value
0 00 1 1 [(1/4\pi)^{1/2}] 0.28209 [1/4\pi] 0.07958
1 [\!\matrix{11+\cr 11-\cr 10\hfill\cr}] [\!\matrix{1\cr 1\cr 1\cr}] [\left.\!\matrix{x\cr y\cr z\cr}\right\}] [(3/4\pi)^{1/2}] 0.48860 [1/\pi] 0.31831
2 20 [1/2] [3z^{2} - 1] [(5/16\pi)^{1/2}] 0.31539 [\displaystyle{3\sqrt{3} \over 8\pi}] 0.20675
[\!\matrix{21+\cr 21-\cr 22+\cr 22-\cr}] [\!\matrix{3\cr 3\cr 6\cr 6\cr}] [\left.\!\matrix{xz\cr yz\cr (x^{2} - y^{2})/2\cr xy\cr}\right\}] [(15/4\pi)^{1/2}] 1.09255 [3/4] 0.75
3 30 [1/2] [5z^{3} - 3z] [(7/16\pi)^{1/2}] 0.37318 [\displaystyle{10 \over 13\pi}] 0.24485
[\!\matrix{31+\cr 31-\cr}] [\!\matrix{3/2\cr 3/2\cr}] [\left.\!\matrix{x[5z^{2} - 1]\cr y[5z^{2} - 1]\cr}\right\}] [(21/32\pi)^{1/2}] 0.45705 [\displaystyle\left(\hbox{ar} + {14 \over 5} - {\pi \over 4}\right)^{-1}] †† 0.32033
[\!\matrix{32+\cr 32-\cr}] [\!\matrix{15\cr 15\cr}] [\left.\!\matrix{(x^{2} - y^{2})z\cr 2xyz\cr}\right\}] [(105/16\pi)^{1/2}] 1.44531 1 1
[\!\matrix{33+\cr 33-\cr}] [\!\matrix{15\cr 15\cr}] [\left.\!\matrix{x^{3} - 3xy^{2}\cr -y^{3} + 3x^{2}y\cr}\right\}] [(35/32\pi)^{1/2}] 0.59004 [4/3\pi] 0.42441
4 40 [1/8] [35z^{4} - 30z^{2} + 3] [(9/256\pi)^{1/2}] 0.10579 ‡‡ 0.06942
[\!\matrix{41+\cr 41-\cr}] [\!\matrix{5/2\cr 5/2\cr}] [\left.\!\matrix{x[7z^{3} - 3z]\cr y[7z^{3} - 3z]\cr}\right\}] [(45/32\pi)^{1/2}] 0.66905 [\displaystyle{735 \over 512\sqrt{7} + 196}] 0.47400
[\!\matrix{42+\cr 42-\cr}] [\!\matrix{15/2\cr 15/2\cr}] [\left.\!\matrix{(x^{2} - y^{2})[7z^{2} - 1]\cr 2xy[7z^{2} - 1]\cr}\right\}] [(45/64\pi)^{1/2}] 0.47309 [\displaystyle{105\sqrt{7} \over 4(136 + 28\sqrt{7})}] 0.33059
[\!\matrix{43+\cr 43-\cr}] [\!\matrix{105\cr 105\cr}] [\left.\!\matrix{(x^{3} - 3xy^{2})z\cr (-y^{3} + 3x^{2}y)z\cr}\right\}] [(315/32\pi)^{1/2}] 1.77013 [5/4] 1.25
[\!\matrix{44+\cr 44-\cr}] [\!\matrix{105\cr 105\cr}] [\left.\!\matrix{x^{4} - 6x^{2}y^{2} + y^{4}\cr 4x^{3}y - 4xy^{3}\cr}\right\}] [(315/256\pi)^{1/2}] 0.62584 [15/32] 0.46875
5 50 [1/8] [63z^{5} - 70z^{3} - 15z] [(11/256\pi)^{1/2}] 0.11695 0.07674
[\!\matrix{51+\cr 51-\cr}] [15/8] [\left.\!\matrix{(21z^{4} - 14z^{2} + 1)x\cr (21z^{4} - 14z^{2} + 1)y\cr}\right\}] [(165/256\pi)^{1/2}] 0.45295 0.32298
[\!\matrix{52+\cr 52-\cr}] [105/2] [\left.\!\matrix{(3z^{3} - z) (x^{2} - y^{2})\cr 2xy(3z^{3} - z)\cr}\right\}] [(1155/64\pi)^{1/2}] 2.39677 1.68750
[\!\matrix{53+\cr 53-\cr}] [105/2] [\left.\!\matrix{(9z^{2} - 1) (x^{3} - 3xy^{2})\cr (9z^{2} - 1) (3x^{2}y - y^{3})\cr}\right\}] [(385/512\pi)^{1/2}] 0.48924 0.34515
[\!\matrix{54+\cr 54-\cr}] [945] [\left.\!\matrix{z(x^{4} - 6x^{2}y^{2} + y^{4})\cr z(4x^{3}y - 4xy^{3})\cr}\right\}] [(3465/256\pi)^{1/2}] 2.07566 1.50000
[\!\matrix{55+\cr 55-\cr}] [945] [\left.\!\matrix{x^{5} - 10x^{3}y^{2} + 5xy^{4}\cr 5x^{4} y - 10x^{2}y^{3} + y^{5}\cr}\right\}] [(693/512\pi)^{1/2}] 0.65638 0.50930
6 60 [1/16] [231z^{6} - 315z^{4} + 105z^{2} - 5] [(13/1024\pi)^{1/2}] 0.06357 0.04171
[\!\matrix{61+\cr 61-\cr}] [21/8] [\left.\!\matrix{(33z^{5} - 30z^{3} + 5z)x\cr (33z^{5} - 30z^{3} + 5z)y\cr}\right\}] [(273/256\pi)^{1/2}] 0.58262 0.41721
[\!\matrix{62+\cr 62-\cr}] [105/8] [\left.\!\matrix{(33z^{4} - 18z^{2} + 1) (x^{2} - y^{2})\cr 2xy (33z^{4} - 18z^{2} + 1)\cr}\right\}] [(1365/2048\pi)^{1/2}] 0.46060 0.32611
[\!\matrix{63+\cr 63-\cr}] [315/2] [\left.\!\matrix{(11z^{3} - 3z) (x^{3} - 3xy^{2})\cr (11z^{3} - 3z) (3x^{2}y - 3y)\cr}\right\}] [(1365/512\pi)^{1/2}] 0.92121 0.65132
[\!\matrix{64+\cr 64-\cr}] [945/2] [\left.\!\matrix{(11z^{2} - 1) (x^{4} - 6x^{2}y^{2} + y^{4})\cr (11z^{2} - 1) (4x^{3}y - 4xy^{3})\cr}\right\}] [(819/1024\pi)^{1/2}] 0.50457 0.36104
[\!\matrix{65+\cr 65-\cr}] 10395 [\left.\!\matrix{z(x^{5} - 10x^{3}y^{2} + 5xy^{4})\cr z(5x^{4}y - 10x^{2}y^{3} + y^{5})\cr}\right\}] [(9009/512\pi)^{1/2}] 2.36662 1.75000
[\!\matrix{66+\cr 66-\cr}] 10395 [\left.\!\matrix{x^{6} - 15x^{4}y^{2} + 15x^{2}y^{4} - y^{6}\cr 6x^{5}y - 20x^{3}y^{3} + 6xy^{5}\cr}\right\}] [(3003/2048\pi)^{1/2}] 0.68318 0.54687
7 70 [1/16] [429z^{7} - 693z^{5} + 315z^{3} - 35z] [(15/1024\pi)^{1/2}] 0.06828 0.04480
[\!\matrix{71+\cr 71-\cr}] [7/16] [\left.\!\matrix{(429z^{6} - 495z^{4} + 135z^{2} - 5)x\cr (429z^{6} - 495z^{4} + 135z^{2} - 5)y\cr}\right\}] [(105/4096\pi)^{1/2}] 0.09033 0.06488
[\!\matrix{72+\cr 72-\cr}] [63/8] [\left.\!\matrix{(143z^{5} - 110z^{3} + 15z) (x^{2} - y^{2})\cr 2xy(143z^{5} - 110z^{3} + 15z)\cr}\right\}] [(315/2048\pi)^{1/2}] 0.22127 0.15732
[\!\matrix{73+\cr 73-\cr}] [315/8] [\left.\!\matrix{(143z^{4} - 66z^{2} + 3) (x^{3} - 3xy^{2})\cr (143z^{4} - 66z^{2} + 3) (3x^{2}y - y^{3})\cr}\right\}] [(315/4096\pi)^{1/2}] 0.15646 0.11092
[\!\matrix{74+\cr 74-\cr}] [3465/2] [\left.\!\matrix{(13z^{3} - 3z) (x^{4} - 6x^{2}y^{2} + y^{4})\cr (13z^{3} - 3z) (4x^{3}y - 4xy^{3})\cr}\right\}] [(3465/1024\pi)^{1/2}] 1.03783 0.74044
[\!\matrix{75+\cr 75-\cr}] [10395/2] [\left.\!\matrix{(13z^{3} - 1) (x^{5} - 10x^{3}y^{2} + 5xy^{4})\cr (13z^{3} - 1) (5x^{4}y - 10x^{2}y^{3} + y^{5})\cr}\right\}] [(3465/4096\pi)^{1/2}] 0.51892 0.37723
[\!\matrix{76+\cr 76-\cr}] 135135 [\left.\!\matrix{z(x^{6} - 15x^{4}y^{2} + 15x^{2}y^{4} - y^{6})\cr z(6x^{5}y + 20x^{3}y^{3} - 6xy^{5})\cr}\right\}] [(45045/2048\pi)^{1/2}] 2.6460 2.00000
[\!\matrix{77+\cr 77-\cr}] 135135 [\left.\!\matrix{x^{7} - 21x^{5}y^{2} + 35x^{3}y^{4} - 7xy^{6}\cr 7x^{6}y - 35x^{4}y^{3} + 21x^{2}y^{5} - y^{7}\cr}\right\}] [(6435/4096\pi)^{1/2}] 0.70716 0.58205
Common factor such that [C_{lm}c_{lmp} = P_{l}^{m} (\cos \theta)_{\sin m\varphi}^{\cos m\varphi}.]
[x = \sin \theta \cos \varphi], [y = \sin \theta \sin \varphi], [z = \cos \theta].
§As defined by [y_{lmp} = M_{lmp}c_{lmp}] where [c_{lmp}] are Cartesian functions.
Paturle & Coppens (1988)[link], as defined by [d_{lmp} = L_{lmp}c_{lmp}] where [c_{lmp}] are Cartesian functions.
††ar = arctan (2).
‡‡[N_{\rm ang} = \{(14A_{-}^{5} - 14A_{+}^{5} + 20A_{+}^{3} - 20A_{-}^{3} + 6A_{-} - 6A_{+}) 2\pi\}^{-1}] where [A_{\pm} = [(30\pm \sqrt{480})/70]^{1/2}].

The spherical harmonic functions are mutually orthogonal and form a complete set, which, if taken to sufficiently high order, can be used to describe any arbitrary angular function.

The spherical harmonic functions are often referred to as multipoles since each represents the components of the charge distribution [\rho ({\bf r})], which gives non-zero contribution to the integral [\Theta_{lmp} = {\textstyle\int} \rho ({\bf r}) c_{lmp} r^{l}\; \hbox{d}{\bf r}], where [\Theta_{lmp}] is an electrostatic multipole moment. Terms with increasing l are referred to as monopolar [(l = 0)], dipolar [(l = 1)], quadrupolar [(l = 2)], octapolar [(l = 3)], hexadecapolar [(l = 4)], triacontadipolar [(l = 5)] and hexacontatetrapolar [(l = 6)].

Site-symmetry restrictions for the real spherical harmonics as given by Kara & Kurki-Suonio (1981)[link] are summarized in Table[link].

Table| top | pdf |
Index-picking rules of site-symmetric spherical harmonics (Kara & Kurki-Suonio, 1981[link])

λ, μ and j are integers.

SymmetryChoice of coordinate axesIndices of allowed [y_{lmp}], [d_{lmp}]
1 Any [{\rm All}\;(l, m, \pm)]
[\bar{1}] Any [(2\lambda, m, \pm)]
2 [2\!\!\parallel\!\!z] [(l, 2\mu, \pm)]
m [m\perp z] [(l, l-2j, \pm)]
[2/m] [2\!\!\parallel\!\! z, m\perp z] [(2\lambda, 2\mu, \pm)]
222 [2\!\!\parallel\!\! z, 2\!\!\parallel\!\! y] [(2\lambda, 2\mu, +)], [(2\lambda + 1, 2\mu, -)]
mm2 [2\!\!\parallel\!\! z, m\perp y] [(l, 2\mu, +)]
mmm [m\perp z, m\perp y, m\perp x] [(2\lambda, 2\mu, +)]
4 [4\!\!\parallel\!\! z] [(l, 4\mu, \pm)]
[\bar{4}] [\bar{4}\!\!\parallel\!\! z] [(2\lambda, 4\mu, \pm)], [(2\lambda + 1, 4\mu + 2, \pm)]
[4/m] [4\!\!\parallel\!\! z, m\perp z] [(2\lambda, 4\mu, \pm)]
422 [4\!\!\parallel\!\! z, 2\!\!\parallel\!\! y] [(2\lambda, 4\mu, +)], [(2\lambda + 1, 4\mu, -)]
4mm [4\!\!\parallel\!\! z, m\perp y] [(l, 4\mu, +)]
[\bar{4}]2m [\bar{4}\!\!\parallel\!\! z, 2\!\!\parallel\!\! x] [(2\lambda, 4\mu, +)], [(2\lambda + 1, 4\mu + 2, -)]
  [m\perp y] [(2\lambda, 4\mu, +)], [(2\lambda + 1, 4\mu + 2, +)]
[4/mmm] [4\!\!\parallel\!\! z, m\perp z, m\perp x] [(2\lambda, 4\mu, +)]
3 [3\!\!\parallel\!\! z] [(l, 3\mu, \pm)]
[\bar{3}] [\bar{3}\!\!\parallel\!\! z] [(2\lambda, 3\mu, \pm)]
32 [3\!\!\parallel\!\! z, 2\!\!\parallel\!\! y] [(2\lambda, 3\mu, +), (2\lambda + 1, 3\mu, -)]
  [2\!\!\parallel\!\! x] [(3\mu + 2j, 3\mu, +)],
    [(3\mu + 2j + 1, 3\mu, -)]
3m [3\!\!\parallel\!\! z, m \perp y] [(l, 3\mu, +)]
  [m \perp x] [(l, 6\mu, +), (l, 6\mu + 3, -)]
[\bar{3}m] [\bar{3}\!\!\parallel\!\! z, m \perp y] [(2\lambda, 3\mu, +)]
  [m \perp x] [(2\lambda, 6\mu, +), (2\lambda, 6\mu + 3, -)]
6 [6\!\!\parallel\!\! z] [(l, 6\mu, \pm)]
[\bar{6}] [\bar{6}\!\!\parallel\!\! z] [(2\lambda, 6\mu, \pm), (2\lambda + 1, 6\mu + 3, \pm)]
[6/m] [6\!\!\parallel\!\! z, m \perp z] [(2\lambda, 6\mu, \pm)]
622 [6\!\!\parallel\!\! z, 2\!\!\parallel\!\! y] [(2\lambda, 6\mu, +), (2\lambda + 1, 6\mu, -)]
6mm [6\!\!\parallel\!\! z, m\!\!\parallel\!\! y] [(l, 6\mu, +)]
[\bar{6}m2] [\bar{6}\!\!\parallel\!\! z, m \perp y] [(2\lambda, 6\mu, +), (2\lambda + 1, 6\mu + 3, +)]
  [m \perp x] [(2\lambda, 6\mu, +), (2\lambda + 1, 6\mu + 3, -)]
[6/mmm] [6\!\!\parallel\!\! z, m \perp z, m \perp y] [(2\lambda, 6\mu, +)]

In cubic space groups, the spherical harmonic functions as defined by equations ([link] [link] [link] are no longer linearly independent. The appropriate basis set for this symmetry consists of the `Kubic Harmonics' of Von der Lage & Bethe (1947)[link]. Some low-order terms are listed in Table[link] Both wavefunction and density-function normalization factors are specified in Table[link].

Table| top | pdf |
`Kubic Harmonic' functions

(a) Coefficients in the expression [K_{lj} = {\textstyle\sum\limits_{mp}} k_{mpj}^{l} y_{lmp}] with normalization [{\textstyle\int_{0}^{\pi}} {\textstyle\int_{0}^{2\pi}} |K_{lj}|^{2} \sin \theta\;\hbox{d}\theta \;\hbox{d}\varphi = 1] (Kara & Kurki-Suonio, 1981[link]).

Even l mp
l j0+2+4+6+8+10+
0 1 1          
4 1 [\textstyle{1 \over 2}\left({7 \over 3}\right)^{1/2}]   [\textstyle{1 \over 2}\left({5 \over 3}\right)^{1/2}]      
0.76376   0.64550      
6 1 [\textstyle{1 \over 2}\left({1 \over 2}\right)^{1/2}]   [\textstyle-{1 \over 2}\left({7 \over 2}\right)^{1/2}]      
0.35355   −0.93541      
6 2   [\textstyle{1 \over 4}11^{1/2}]   [\textstyle- {1 \over 4} 5^{1/2}]    
  0.82916   −0.55902    
8 1 [\textstyle{1 \over 8}33^{1/2}]   [\textstyle{1 \over 4}\left({7 \over 3}\right)^{1/2}]   [\textstyle{1 \over 8}\left({65 \over 3}\right)^{1/2}]  
0.71807   0.38188   0.58184  
10 1 [\textstyle{1 \over 8}\left({65 \over 6}\right)^{1/2}]   [\textstyle- {1 \over 4}\left({11 \over 2}\right)^{1/2}]   [\textstyle- {1 \over 8}\left({187 \over 6}\right)^{1/2}]  
0.41143   −0.58630   −0.69784  
10 2   [\textstyle{1 \over 8}\left({247 \over 6}\right)^{1/2}]   [\textstyle{1 \over 16}\left({19 \over 3}\right)^{1/2}]   [\textstyle{1 \over 16}85^{1/2}]
  0.80202   0.15729   0.57622
l j   2− 4− 6− 8−
3 1   1      
7 1   [\textstyle{1 \over 2}\left({13 \over 6}\right)^{1/2}]   [\textstyle{1 \over 2}\left({11 \over 16}\right)^{1/2}]  
  0.73598   0.41458  
9 1   [\textstyle{1 \over 4}3^{1/2}]   [\textstyle- {1 \over 4} 13^{1/2}]  
  0.43301   −0.90139  
9 2   [\textstyle{1 \over 2}\left({17 \over 6}\right)^{1/2}]   [\textstyle- {1 \over 2}\left({7 \over 6}\right)^{1/2}]  
  0.84163   −0.54006  

(b) Coefficients [k_{mpj}^{l}] and density normalization factors [N_{lj}] in the expression [K_{lj} = N_{lj} {\textstyle\sum\limits_{mp}} k_{mpj}^{l} u_{lmp}] where [u_{lm \pm} =P_{l}^{m} (\cos \theta)^{\cos m\varphi}_{\sin m\varphi}] (Su & Coppens, 1994[link]).

Even l [N_{lj}] mp
l j   0+ 2+ 4+ 6+ 8+ 10+
0 1 [1/4\pi = 0.079577] 1          
4 1 0.43454 1   [+1/168]      
6 1 0.25220 1   [-1/360]      
6 2 0.020833   1   [-1/792]    
8 1 0.56292 1   1/5940   [\textstyle{1 \over 672} \times {1 \over 5940}]  
10 1 0.36490 1   1/5460   [\textstyle{1 \over 4320} \times {1 \over 5460}]  
10 2 0.0095165 1     [1/43680]   [\textstyle- {1 \over 456} \times {1 \over 43680}]
l j     2− 4− 6− 8−
3 1 0.066667   1      
7 1 0.014612   1   [1/1560]  
9 1 0.0059569   1   [1/2520]  
9 2 0.00014800     1   [-1/4080]

(c) Density-normalized Kubic harmonics as linear combinations of density-normalized spherical harmonic functions. Coefficients in the expression [K_{lj} = {\textstyle\sum\limits_{mp}} k^{''l}_{mpj} d_{lmp}]. Density-type normalization is defined as [{\textstyle\int_{0}^{\pi}} {\textstyle\int_{0}^{2\pi}} |K_{lj}| \sin \theta\ \hbox{d} \theta\ \hbox{d} \varphi = 2 - \delta_{l0}].

Even l mp
l j0+2+4+6+8+10+
0 1 1          
4 1 0.78245   0.57939      
6 1 0.37790   −0.91682      
6 2   0.83848   −0.50000    
l j 2− 4− 6− 8−
3 1 1      
7 1 0.73145   0.63290  

(d) Index rules for cubic symmetries (Kurki-Suonio, 1977[link]; Kara & Kurki-Suonio, 1981[link]).

lj23 [m\bar{3}]432 [\bar{4}3m] [m\bar{3}m]
T [T_{h}] O [T_{d}] [O_{h}]
0 1 × × × × ×
3 1 ×     ×  
4 1 × × × × ×
6 1 × × × × ×
6 2 × ×      
7 1 ×     ×  
8 1 × × × × ×
9 1 ×     ×  
9 2 ×   ×    
10 1 × × × × ×
10 2 × ×      

A related basis set of angular functions has been proposed by Hirshfeld (1977)[link]. They are of the form [\cos^{n} \theta_{k}], where [\theta_{k}] is the angle with a specified set of [(n + 1)(n + 2)/2] polar axes. The Hirshfeld functions are identical to a sum of spherical harmonics with [l = n], [n - 2], [n - 4,\ldots (0, 1)] for [n \;\gt\; 1], as shown elsewhere (Hirshfeld, 1977[link]).

The radial functions [R(r)] can be selected in different manners. Several choices may be made, such as [{R_{l}(r) = {\zeta^{n_{l} + 3} \over (n_{l} + 2)!} r^{n(l)}\exp (-\zeta_{l}r) \qquad \hbox{(Slater type function)},} \eqno(] where the coefficient [n_{l}] may be selected by examination of products of hydrogenic orbitals which give rise to a particular multipole (Hansen & Coppens, 1978[link]). Values for the exponential coefficient [\zeta_{l}] may be taken from energy-optimized coefficients for isolated atoms available in the literature (Clementi & Raimondi, 1963[link]). A standard set has been proposed by Hehre et al. (1969)[link]. In the bonded atom, such values are affected by changes in nuclear screening due to migrations of charge, as described in part by equation ([link].

Other alternatives are: [R_{l} (r) = {\alpha^{n + 1} \over n!} r^{n} \exp (-\alpha r^{2})\qquad (\hbox{Gaussian function}) \eqno(] or [{R_{l} (r) = r^{l} L_{n}^{2l + 2} (\gamma r) \exp \left(-{\gamma r \over 2}\right)\quad (\hbox{Laguerre function}),} \eqno(] where L is a Laguerre polynomial of order n and degree [(2l + 2)].

In summary, in the multipole formalism the atomic density is described by [\eqalignno{\rho_{\rm atomic}({\bf r}) &= P_{c} \rho_{\rm core} + P_{\nu} \kappa^{3} \rho_{\rm valence} (\kappa r)\cr &\quad + {\textstyle\sum\limits_{l = 0}^{l_{\max}}} \kappa'^{3} R_{l}(\kappa' r) {\textstyle\sum\limits_{m = 0}^{l}} {\textstyle\sum\limits_{p}} P_{lmp} d_{lmp} ({\bf r}/r), &(}] in which the leading terms are those of the kappa formalism [expressions ([link], ([link]]; the subscript p is either + or −.

The expansion in ([link] is frequently truncated at the hexadecapolar [(l = 4)] level. For atoms at positions of high site symmetry the first allowed functions may occur at higher l values. For trigonally bonded atoms in organic molecules the [l = 3] terms are often found to be the most significantly populated deformation functions. Reciprocal-space description of aspherical atoms

| top | pdf |

The aspherical-atom form factor is obtained by substitution of ([link] in expression ([link]: [f_{j}({\bf S}) = {\textstyle\int} \rho_{j}({\bf r}) \exp (2\pi i{\bf S} \cdot {\bf r})\ \hbox{d}{\bf r}. \eqno(] In order to evaluate the integral, the scattering operator [{\exp (2\pi i{\bf S} \cdot {\bf r})}] must be written as an expansion of products of spherical harmonic functions. In terms of the complex spherical harmonic functions, the appropriate expression is (Weiss & Freeman, 1959[link]; Cohen-Tannoudji et al., 1977[link]) [{\exp (2\pi i{\bf S} \cdot {\bf r}) = 4\pi {\textstyle\sum\limits_{l = 0}^{\infty}}\; {\textstyle\sum\limits_{m = -l}^{l}} i^{l} j_{l} (2\pi Sr) Y_{lm} (\theta, \varphi) Y_{lm}^{*} (\beta, \gamma).} \eqno(]

The Fourier transform of the product of a complex spherical harmonic function with normalization [{\textstyle\int} |Y_{lm}|^{2}\ \hbox{d}\Omega = 1] and an arbitrary radial function [R_{l}(r)] follows from the orthonormality properties of the spherical harmonic functions, and is given by [{{\textstyle\int} Y_{lm} R_{l}(r) \exp (2\pi i {\bf S} \cdot {\bf r})\ \hbox{d}\tau = 4\pi i^{l}{\textstyle\int}j_{l} (2\pi {S}r) R_{l}(r) r^{2}\ \hbox{d}r Y_{lm} (\beta, \gamma),} \eqno(] where [j_{l}] is the lth-order spherical Bessel function (Arfken, 1970[link]), and θ and ϕ, β and γ are the angular coordinates of r and S, respectively.

For the Fourier transform of the real spherical harmonic functions, the scattering operator is expressed in terms of the real spherical harmonics: [\eqalignno{\exp (2\pi i{\bf S} \cdot {\bf r}) &= {\sum\limits_{l = 0}^{\infty}} i^{l} j_{l} (2\pi Sr) (2 - \delta_{m0}) (2l + 1) {\sum\limits_{m = 0}^{l}} {(l - m)! \over (l + m)!}\cr &\quad \times P_{l}^{m} (\cos \theta) P_{l}^{m} (\cos \beta) \cos [m(\phi - \gamma)], &(}] which leads to [{\textstyle\int} y_{lmp} (\theta, \varphi) R_{l}(r) \exp (2\pi i{\bf S} \cdot {\bf r})\ \hbox{d}\tau = 4\pi i^{l}\langle j_{l}\rangle y_{lmp} (\beta, \gamma). \eqno(] Since [y_{lmp}] occurs on both sides, the expression is independent of the normalization selected. Therefore, for the Fourier transform of the density functions [d_{lmp}] [{\textstyle\int} d_{lmp} (\theta, \varphi) R_{l}(r) \exp (2\pi i{\bf S} \cdot {\bf r})\ \hbox{d}\tau = 4\pi i^{l}\langle j_{l}\rangle d_{lmp} (\beta, \gamma). \eqno(]

In ([link] and ([link], [\langle j_{l}\rangle], the Fourier–Bessel transform, is the radial integral defined as [\langle j_{l}\rangle = {\textstyle\int} j_{l}(2\pi Sr) R_{l}(r) r^{2}\ \hbox{d}r \eqno(] of which [\langle j_{0}\rangle] in expression ([link] is a special case. The functions [\langle j_{l}\rangle] for Hartree–Fock valence shells of the atoms are tabulated in scattering-factor tables (IT IV, 1974[link]). Expressions for the evaluation of [\langle j_{l}\rangle] using the radial function ([link][link]c[link]) have been given by Stewart (1980)[link] and in closed form for ([link] by Avery & Watson (1977)[link] and Su & Coppens (1990)[link]. The closed-form expressions are listed in Table[link].

Table| top | pdf |
Closed-form expressions for Fourier transform of Slater-type functions (Avery & Watson, 1977[link]; Su & Coppens, 1990[link])

[\langle j_{k}\rangle \equiv {\textstyle\int_{0}^{\infty}} r^{N} \exp(-Zr)j_{k}(Kr)\;\hbox{d}r, K = 4\pi \sin \theta/\lambda.]

k 1 2 3 4 5 6 7 8
0 [\displaystyle{1 \over K^{2} + Z^{2}}] [\displaystyle{2Z \over (K^{2} + Z^{2})^{2}}] [\displaystyle{2(3Z^{2} - K^{2}) \over (K^{2} + Z^{2})^{3}}] [\displaystyle{24Z(Z^{2} - K^{2}) \over (K^{2} + Z^{2})^{4}}] [\displaystyle{24(5Z^{2} - 10K^{2}Z^{2} + K^{4}) \over (K^{2} + Z^{2})^{5}}] [\displaystyle{240Z(K^{2} - 3Z^{2}) (3K^{2} - Z^{2}) \over (K^{2} + Z^{2})^{6}}] [\displaystyle{720(7Z^{6} - 35K^{2}Z^{4} + 21K^{4}Z^{2} - K^{6}) \over (K^{2} + Z^{2})^{7}}] [\displaystyle{40320(Z^{7} - 7K^{2}Z^{5} + 7K^{4}Z^{3} - K^{6}Z) \over (K^{2} + Z^{2})^{8}}]
1   [\displaystyle{2K \over (K^{2} + Z^{2})^{2}}] [\displaystyle{8KZ \over (K^{2} + Z^{2})^{3}}] [\displaystyle{8K(5Z^{2} - K^{2}) \over (K^{2} + Z^{2})^{4}}] [\displaystyle{48KZ(5Z^{2} - 3K^{2}) \over (K^{2} + Z^{2})^{5}}] [\displaystyle{48K(35Z^{4} - 42K^{2}Z^{2} + 3K^{4}) \over (K^{2} + Z^{2})^{6}}] [\displaystyle{1920KZ(7Z^{4} - 14K^{2}Z^{2} + 3K^{4}) \over (K^{2} + Z^{2})^{7}}] [\displaystyle{5760K(21Z^{6} - 63K^{2}Z^{4} + 27K^{4}Z^{2} - K^{6}) \over (K^{2} + Z^{2})^{8}}]
2     [\displaystyle{8K^{2} \over (K^{2} + Z^{2})^{3}}] [\displaystyle{48K^{2}Z \over (K^{2} + Z^{2})^{4}}] [\displaystyle{48K^{2}(7Z^{2} - K^{2}) \over (K^{2} + Z^{2})^{5}}] [\displaystyle{384K^{2}Z(7Z^{2} - 3K^{2}) \over (K^{2} + Z^{2})^{6}}] [\displaystyle{1152K^{2}(21Z^{4} - 18K^{2}Z^{2} + K^{4}) \over (K^{2} + Z^{2})^{7}}] [\displaystyle{11520K^{2}Z(21Z^{4} - 30K^{2}Z^{2} + 5K^{4}) \over (K^{2} + Z^{2})^{8}}]
3       [\displaystyle{48K^{3} \over (K^{2} + Z^{2})^{4}}] [\displaystyle{384K^{3}Z \over (K^{2} + Z^{2})^{5}}] [\displaystyle{384K^{3}(9Z^{2} - K^{2}) \over (K^{2} + Z^{2})^{6}}] [\displaystyle{11520K^{3}Z(3Z^{2} - K^{2}) \over (K^{2} + Z^{2})^{7}}] [\displaystyle{11520K^{3}(33Z^{4} - 22K^{2}Z^{2} + K^{4}) \over (K^{2} + Z^{2})^{8}}]
4         [\displaystyle{384K^{4} \over (K^{2} + Z^{2})^{5}}] [\displaystyle{3840K^{4}Z \over (K^{2} + Z^{2})^{6}}] [\displaystyle{3840K^{4}(11Z^{2} - K^{2}) \over (K^{2} + Z^{2})^{7}}] [\displaystyle{46080K^{4}Z(11Z^{2} - 3K^{2}) \over (K^{2} + Z^{2})^{8}}]
5           [\displaystyle{3840K^{5} \over (K^{2} + Z^{2})^{6}}] [\displaystyle{46080K^{5}Z \over (K^{2} + Z^{2})^{7}}] [\displaystyle{40680K^{5}(13Z^{2} - K^{2}) \over (K^{2} + Z^{2})^{8}}]
6             [\displaystyle{46080K^{6} \over (K^{2} + Z^{2})^{7}}] [\displaystyle{645120K^{6}Z \over (K^{2} + Z^{2})^{8}}]
7               [\displaystyle{645120K^{7} \over (K^{2} + Z^{2})^{8}}]

Expressions ([link] [link] show that the Fourier transform of a direct-space spherical harmonic function is a reciprocal-space spherical harmonic function with the same l, m, or, in other words, the spherical harmonic functions are Fourier-transform invariant.

The scattering factors [f_{lmp}({\bf S})] of the aspherical density functions [R_{l}(r)d_{lmp}(\theta, \phi)] in the multipole expansion ([link] are thus given by [f_{lmp}({\bf S}) = 4\pi i^{l} \langle j_{l}\rangle d_{lmp} (\beta, \gamma). \eqno(]

The reciprocal-space spherical harmonic functions in this expression are identical to the functions given in Table[link], except for the replacement of the direction cosines x, y and z by the direction cosines of the scattering vector S.

1.2.8. Fourier transform of orbital products

| top | pdf |

If the wavefunction is written as a sum over normalized Slater determinants, each representing an antisymmetrized combination of occupied molecular orbitals [\chi_{i}] expressed as linear combinations of atomic orbitals [\varphi_{\nu}], i.e. [\chi_{i} = {\textstyle\sum\limits_{\nu}} \hbox{c}_{i\nu} \varphi_{\nu}], the electron density is given by (Stewart, 1969a[link]) [\rho ({\bf r}) = {\textstyle\sum\limits_{i}} n_{i} \chi_{i}^{2} = {\textstyle\sum\limits_{\mu}} {\textstyle\sum\limits_{\nu}} P_{\mu \nu} \varphi_{\mu} ({\bf r}) \varphi_{\nu} ({\bf r}), \eqno(] with [n_{i} = 1\hbox{ or }2]. The coefficients [P_{\mu \nu}] are the populations of the orbital product density functions [\phi_{\mu} ({\bf r}) \varphi_{\nu} ({\bf r})] and are given by [P_{\mu \nu} = {\textstyle\sum\limits_{i}} n_{i} c_{i\mu} c_{i\nu}. \eqno(]

For a multi-Slater determinant wavefunction the electron density is expressed in terms of the occupied natural spin orbitals, leading again to ([link] but with non-integer values for the coefficients [n_{i}].

The summation ([link] consists of one- and two-centre terms for which [\varphi_{\mu}] and [\varphi_{\nu}] are centred on the same or on different nuclei, respectively. The latter represent the overlap density, which is only significant if [\varphi_{\mu} ({\bf r})] and [\varphi_{\nu} ({\bf r})] have an appreciable value in the same region of space. One-centre orbital products

| top | pdf |

If the atomic basis consists of hydrogenic type s, p, d, f, … orbitals, the basis functions may be written as [\varphi (r, \theta, \varphi) = R_{l} (r) Y_{lm} (\theta, \varphi) \eqno(] or [\varphi (r, \theta, \varphi) = R_{l} (r) y_{lmp} (\theta, \varphi), \eqno(] which gives for corresponding values of the orbital products [\varphi_{\mu} ({\bf r}) \varphi_{\nu} ({\bf r}) = R_{l} (r) R_{l^{'}} (r) Y_{lm} (\theta, \varphi) Y_{l^{'}m^{'}} (\theta, \varphi) \eqno(] and [\varphi_{\mu} ({\bf r}) \varphi_{\nu} ({\bf r}) = R_{l} (r) R_{l^{'}} (r) y_{lmp} (\theta, \varphi) y_{l^{'}m^{'}p^{'}} (\theta, \varphi), \eqno(] respectively, where it has been assumed that the radial function depends only on l.

Because the spherical harmonic functions form a complete set, their products can be expressed as a linear combination of spherical harmonics. The coefficients in this expansion are the ClebschGordan coefficients (Condon & Shortley, 1957[link]), defined by [Y_{lm} (\theta, \varphi) Y_{l^{'}m^{'}} (\theta, \varphi) = {\textstyle\sum\limits_{L}} {\textstyle\sum\limits_{M}} C_{Lll^{'}}^{M mm^{'}} Y_{LM} (\theta, \varphi) \eqno(] or the equivalent definition [C_{Lll^{'}}^{M mm^{'}} = {\textstyle\int\limits_{0}^{\pi}} \sin \theta \;\hbox{d} \theta {\textstyle\int\limits_{0}^{2\pi}}\; \hbox{d}\varphi Y_{LM}^{*} (\theta, \varphi) Y_{lm} (\theta, \varphi) Y_{l^{'}m^{'}} (\theta, \varphi). \eqno(] The [C_{Lll^{'}}^{M mm^{'}}] vanish, unless [L + l + l^{'}] is even, [|l - l^{'} | \;\lt\; L \;\lt\; l + l^{'}] and [M = m + m^{'}].

The corresponding expression for [y_{lmp}] is [y_{lmp} (\theta, \varphi) y_{l^{'}m^{'}p^{'}} (\theta, \varphi) = {\textstyle\sum\limits_{L}} {\textstyle\sum\limits_{M}} C^{'} {\textstyle {\openup-4pt{\matrix{_{M mm'}\hfill\cr _{Lll'}\hfill\cr _{P}\hfill}}}} y_{LMP} (\theta, \varphi), \eqno(] with [M = |m + m^{'}|] and [|m - m^{'}|] for [p = p^{'}], and [M = - |m + m^{'}|] and [- |m - m^{'}|] for [p = - p^{'}] and [P = p \times p^{'}].

Values of C and [C^{'}] for [l \leq 2] are given in Tables[link] and[link] They are valid for the functions [Y_{lm}] and [y_{lmp}] with normalization [{\textstyle\int} |Y_{lm} |^{2}\; \hbox{d} \Omega = 1] and [{\textstyle\int} y_{lmp}^{2}\; \hbox{d} \Omega = 1].

Table| top | pdf |
Products of complex spherical harmonics as defined by equation ([link]

Y 00 Y 00 = 0.28209479Y 00
Y 10 Y 00 = 0.28209479Y 10
Y 10 Y 10 = 0.25231325Y 20 + 0.28209479Y 00
Y 11 Y 00 = 0.28209479Y 11
Y 11 Y 10 = 0.21850969Y 21
Y 11 Y 11 = 0.30901936Y 22
Y 11 Y 11− = −0.12615663Y 20 + 0.28209479Y 00
Y 20 Y 00 = 0.28209479Y 20
Y 20 Y 10 = 0.24776669Y 30 + 0.25231325Y 10
Y 20 Y 11 = 0.20230066Y 31 − 0.12615663Y 11
Y 20 Y 20 = 0.24179554Y 40 + 0.18022375Y 20 + 0.28209479Y 00
Y 21 Y 00 = 0.28209479Y 21
Y 21 Y 10 = 0.23359668Y 31 + 0.21850969Y 11
Y 21 Y 11 = 0.26116903Y 32
Y 21 Y 11− = −0.14304817Y 30 + 0.21850969Y 10
Y 21 Y 20 = 0.22072812Y 41 + 0.09011188Y 21
Y 21 Y 21 = 0.25489487Y 42 + 0.22072812Y 22
Y 21 Y 21− = −0.16119702Y 40 + 0.09011188Y 20 + 0.28209479Y 00
Y 22 Y 00 = 0.28209479Y 22
Y 22 Y 10 = 0.18467439Y 32
Y 22 Y 11 = 0.31986543Y 33
Y 22 Y 11− = −0.08258890Y 31 + 0.30901936Y 11
Y 22 Y 20 = 0.15607835Y 42 − 0.18022375Y 22
Y 22 Y 21 = 0.23841361Y 43
Y 22 Y 21− = −0.09011188Y 41 + 0.22072812Y 21
Y 22 Y 22 = 0.33716777Y 44
Y 22 Y 22− = 0.04029926Y 40 − 0.18022375Y 20 + 0.28209479Y 00

Table| top | pdf |
Products of real spherical harmonics as defined by equations ([link] and ([link]

y 00 y 00 = 0.28209479y 00
y 10 y 00 = 0.28209479y 10
y 10 y 10 = 0.25231325y 20 + 0.28209479y 00
y 11± y 00 = 0.28209479y 11±
y 11± y 10 = 0.21850969y 21±
y 11± y 11± = 0.21850969y 22+ − 0.12615663y 20 + 0.28209479y 00
y 11+ y 11− = 0.21850969y 22−
y 20 y 00 = 0.28209479y 20
y 20 y 10 = 0.24776669y 30 + 0.25231325y 10
y 20 y 11± = 0.20230066y 31± − 0.12615663y 11±
y 20 y 20 = 0.24179554y 40 + 0.18022375y 20 + 0.28209479y 00
y 21± y 00 = 0.28209479y 21±
y 21± y 10 = 0.23359668y 31± + 0.21850969y 11±
y 21± y 11± = ± 0.18467439y 32+ − 0.14304817y 30 + 0.21850969y 10
y 21± y 11∓ = 0.18467469y 32−
y 21± y 20 = 0.22072812y 41± + 0.09011188y 21±
y 21± y 21± = ± 0.18022375y 42+ ± 0.15607835y 22+ − 0.16119702y 40 + 0.09011188y 20 + 0.28209479y 00
y 21+ y 21− = −0.18022375y 42− + 0.15607835y 22−
y 22± y 00 = 0.28209479y 22±
y 22± y 10 = 0.18467439y 32±
y 22± y 11± = ± 0.22617901y 33+ − 0.05839917y 31+ + 0.21850969y 11+
y 22± y 11∓ = 0.22617901y 33− ± 0.05839917y 31− ∓ 0.21850969y 11−
y 22± y 20 = 0.15607835y 42± − 0.18022375y 22±
y 22± y 21± = ± 0.16858388y 43+ − 0.06371872y 41+ + 0.15607835y 21+
y 22± y 21∓ = 0.16858388y 43− ± 0.06371872y 41− ∓ 0.15607835y 21−
y 22± y 22± = ± 0.23841361y 44+ + 0.04029926y 40 − 0.18022375y 20 + 0.28209479y 00
y 22+ y 22− = 0.23841361y 44−

By using ([link] or ([link], the one-centre orbital products are expressed as a sum of spherical harmonic functions. It follows that the one-centre orbital product density basis set is formally equivalent to the multipole description, both in real and in reciprocal space. To obtain the relation between orbital products and the charge-density functions, the right-hand side of ([link] has to be multiplied by the ratio of the normalization constants, as the wavefunctions [y_{lmp}] and charge-density functions [d_{lmp}] are normalized in a different way as described by ([link] and ([link]. Thus [y_{lmp} (\theta, \varphi) y_{l^{'}m^{'}p^{'}} (\theta, \varphi) = {\textstyle\sum\limits_{L}} {\textstyle\sum\limits_{M}} R_{LMP} C^{'} {\let\normalbaselines\relax\openup-4pt{\matrix{_{M mm'}\hfill\cr_{Lll'}\hfill\cr_{P}\hfill}}} d_{LMP}(\theta, \varphi), \eqno(] where [R_{LMP} = M_{LMP} \hbox{ (wavefunction)}/L_{LMP}\hbox{ (density function)}]. The normalization constants [M_{lmp}] and [L_{lmp}] are given in Table[link], while the coefficients in the expressions ([link] are listed in Table[link].

Table| top | pdf |
Products of two real spherical harmonic functions [y_{lmp}] in terms of the density functions [d_{lmp}] defined by equation ([link]

y 00 y 00 = 1.0000d 00
y 10 y 00 = 0.43301d 10
y 10 y 10 = 0.38490d 20 + 1.0d 00
y 11± y 00 = 0.43302d 11±
y 11± y 10 = 0.31831d 21±
y 11± y 11± = 0.31831d 22+ − 0.19425d 20 + 1.0d 00
y 11+ y 11− = 0.31831d 22−
y 20 y 00 = 0.43033d 20
y 20 y 10 = 0.37762d 30 + 0.38730d 10
y 20 y 11± = 0.28864d 31± − 0.19365d 11±
y 20 y 20 = 0.36848d 40 + 0.27493d 20 + 1.0d 00
y 21± y 00 = 0.41094d 21±
y 21± y 10 = 0.33329d 31± + 0.33541d 11±
y 21± y 11± = ±0.26691d 32+ − 0.21802d 30 + 0.33541d 10
y 21± y 11∓ = −0.26691d 32−
y 21± y 20 = 0.31155d 41± + 0.13127d 21±
y 21± y 21± = ±0.25791d 42+ ± 0.22736d 22+ − 0.24565d 40 + 0.13747d 20 + 1.0d 00
y 21+ y 21− = 0.25790d 42− + 0.22736d 22−
y 22± y 00 = 0.41094d 22±
y 22± y 10 = 0.26691d 32±
y 22± y 11± = ± 0.31445d 33+ − 0.083323d 31+ + 0.33541d 11+
y 22± y 11∓ = 0.31445d 33− ± 0.083323d 31− ∓ 0.33541d 11−
y 22± y 20 = 0.22335d 42± − 0.26254d 22±
y 22± y 21± = ± 0.23873d 43+ − 0.089938d 41+ + 0.22736d 21+
y 22± y 21∓ = 0.23873d 43− ± 0.089938d 41− ∓ 0.22736d 21−
y 22± y 22± = ± 0.31831d 44+ + 0.061413d 40 − 0.27493d 20 + 1.0d 00
y 22+ y 22− = 0.31831d 44− Two-centre orbital products

| top | pdf |

Fourier transform of the electron density as described by ([link] requires explicit expressions for the two-centre orbital product scattering. Such expressions are described in the literature for both Gaussian (Stewart, 1969b[link]) and Slater-type (Bentley & Stewart, 1973[link]; Avery & Ørmen, 1979[link]) atomic orbitals. The expressions can also be used for Hartree–Fock atomic functions, as expansions in terms of Gaussian- (Stewart, 1969b[link], 1970[link]; Stewart & Hehre, 1970[link]; Hehre et al., 1970[link]) and Slater-type (Clementi & Roetti, 1974[link]) functions are available for many atoms.

1.2.9. The atomic temperature factor

| top | pdf |

Since the crystal is subject to vibrational oscillations, the observed elastic scattering intensity is an average over all normal modes of the crystal. Within the Born–Oppenheimer approximation, the theoretical electron density should be calculated for each set of nuclear coordinates. An average can be obtained by taking into account the statistical weight of each nuclear configuration, which may be expressed by the probability distribution function [P({\bf u}_{1}, \ldots , {\bf u}_{N})] for a set of displacement coordinates [{\bf u}_{1}, \ldots , {\bf u}_{N}].

In general, if [\rho ({\bf r},{\bf u}_{1}, \ldots , {\bf u}_{N})] is the electron density corresponding to the geometry defined by [{\bf u}_{1}, \ldots , {\bf u}_{N}], the time-averaged electron density is given by [\langle \rho ({\bf r})\rangle = {\textstyle\int} \rho ({\bf r}, {\bf u}_{1}, \ldots , {\bf u}_{N}) P ({\bf u}_{1}, \ldots , {\bf u}_{N})\; \hbox{d} {\bf u}_{1} \ldots \hbox{d} {\bf u}_{N}. \eqno(]

When the crystal can be considered as consisting of perfectly following rigid entities, which may be molecules or atoms, expression ([link] simplifies: [{\langle \rho_{\rm rigid\;group} ({\bf r})\rangle = {\textstyle\int} \rho_{\rm r.g.,\;static} ({\bf r} - {\bf u}) P ({\bf u})\ \hbox{d} {\bf u} = \rho_{\rm r.g.,\;static} * P ({\bf u}).} \eqno(]

In the approximation that the atomic electrons perfectly follow the nuclear motion, one obtains [\langle \rho_{\rm atom} ({\bf r})\rangle = \rho_{\rm atom, \, static} ({\bf r}) * P ({\bf u}). \eqno(] The Fourier transform of this convolution is the product of the Fourier transforms of the individual functions: [\langle f ({\bf H})\rangle = f({\bf H})T({\bf H}). \eqno(] Thus [T({\bf H})], the atomic temperature factor, is the Fourier transform of the probability distribution [P({\bf u})].

1.2.10. The vibrational probability distribution and its Fourier transform in the harmonic approximation

| top | pdf |

For a harmonic oscillator, the probability distribution averaged over all populated energy levels is a Gaussian, centred at the equilibrium position. For the three-dimensional isotropic harmonic oscillator, the distribution is [P(u) = (2\pi \langle u^{2}\rangle)^{-3/2} \exp \{-| u |^{2} / 2 \langle u^{2}\rangle\}, \eqno(] where [\langle u^{2}\rangle] is the mean-square displacement in any direction.

The corresponding trivariate normal distribution to be used for anisotropic harmonic motion is, in tensor notation, [P({\bf u}) = {|{\boldsigma}^{-1}|^{1/2} \over (2\pi)^{3/2}} \exp \{-\textstyle{1 \over 2} {\boldsigma}_{jk}^{-1} (u\hskip 2pt^{j}u^{k})\}. \eqno(] Here σ is the variance–covariance matrix, with covariant components, and [|{\boldsigma}^{-1}|] is the determinant of the inverse of σ. Summation over repeated indices has been assumed. The corresponding equation in matrix notation is [P({\bf u}) = {|{\boldsigma}^{-1}|^{1/2} \over (2\pi)^{3/2}} \exp \{\textstyle-{1 \over 2} ({\bf u})^{T} {\boldsigma}^{-1} ({\bf u})\}, \eqno(] where the superscript T indicates the transpose.

The characteristic function, or Fourier transform, of [P({\bf u})] is [T({\bf H}) = \exp \{-2\pi^{2} \sigma\hskip 2pt^{jk} h_{j}h_{k}\} \eqno(] or [T({\bf H}) = \exp \{-2\pi^{2} {\bf H}^{T} {\bf \boldsigma H}\}. \eqno(] With the change of variable [b\hskip 2pt^{jk} = 2\pi^{2} \sigma\hskip 2pt^{jk}], ([link] becomes [T({\bf H}) = \exp \{-b\hskip 2pt^{jk} h_{j}h_{k}\}.]

1.2.11. Rigid-body analysis

| top | pdf |

The treatment of rigid-body motion of molecules or molecular fragments was developed by Cruickshank (1956)[link] and expanded into a general theory by Schomaker & Trueblood (1968)[link]. The theory has been described by Johnson (1970b)[link] and by Dunitz (1979)[link]. The latter reference forms the basis for the following treatment.

The most general motions of a rigid body consist of rotations about three axes, coupled with translations parallel to each of the axes. Such motions correspond to screw rotations. A libration around a vector [{\boldlambda}\ (\lambda_{1},\ \lambda_{2},\ \lambda_{3})], with length corresponding to the magnitude of the rotation, results in a displacement [\delta {\bf r}], such that [\specialfonts \delta {\bf r} = ({\boldlambda} \times {\bf r}) = {{\bsf D}}{\bf r} \eqno(] with [\specialfonts{{\bsf D}} = \left[\matrix{0 &-\lambda_{3} &\lambda_{2}\cr \lambda_{3} &0 &-\lambda_{1}\cr -\lambda_{2} &\lambda_{1} &0\cr}\right], \eqno(] or in tensor notation, assuming summation over repeated indices, [\delta r_{i} = D_{ij}r_{j} = - \varepsilon_{ijk} \lambda_{k} r_{j} \eqno(] where the permutation operator [\varepsilon_{ijk}] equals +1 for i, j, k a cyclic permutation of the indices 1, 2, 3, or −1 for a non-cyclic permutation, and zero if two or more indices are equal. For [i = 1], for example, only the [\varepsilon_{123}] and [\varepsilon_{132}] terms occur. Addition of a translational displacement gives [\delta r_{i} = D_{ij}r_{j} + t_{i}. \eqno(]

When a rigid body undergoes vibrations the displacements vary with time, so suitable averages must be taken to derive the mean-square displacements. If the librational and translational motions are independent, the cross products between the two terms in ([link] average to zero and the elements of the mean-square displacement tensor of atom n, [{U^{n}_{ij}}], are given by [\eqalign{U_{11}^{n} &= + L_{22} r_{3}^{2} + L_{33} r_{2}^{2} - 2L_{23} r_{2} r_{3} + T_{11}\cr U_{22}^{n} &= + L_{33} r_{1}^{2} + L_{11} r_{3}^{2} - 2L_{13} r_{1} r_{3} + T_{22}\cr U_{33}^{n} &= + L_{11} r_{2}^{2} + L_{22} r_{1}^{2} - 2L_{12} r_{1} r_{2} + T_{33}\cr U_{12}^{n} &= - L_{33} r_{1} r_{2} - L_{12} r_{3}^{2} + L_{13} r_{2} r_{3} + L_{23} r_{1} r_{3} + T_{12}\cr U_{13}^{n} &= - L_{22} r_{1} r_{3} + L_{12} r_{2} r_{3} - L_{13} r_{2}^{2} + L_{23} r_{1} r_{2} + T_{13}\cr U_{23}^{n} &= - L_{11} r_{2} r_{3} + L_{12} r_{1} r_{3} - L_{13} r_{1} r_{2} - L_{23} r_{1}^{2} + T_{23},\cr} \eqno(] where the coefficients [L_{ij} = \langle \lambda_{i} \lambda_{j}\rangle] and [T_{ij} = \langle t_{i} t_{j}\rangle] are the elements of the [3 \times 3] libration tensor [ \specialfonts{{\bsf L}}] and the [3 \times 3] translation tensor [ \specialfonts{{\bsf T}}], respectively. Since pairs of terms such as [\langle t_{i} t_{j}\rangle] and [\langle t_{j} t_{i}\rangle] correspond to averages over the same two scalar quantities, the [\specialfonts {{\bsf T}}] and [\specialfonts {{\bsf L}}] tensors are symmetrical.

If a rotation axis is correctly oriented, but incorrectly positioned, an additional translation component perpendicular to the rotation axes is introduced. The rotation angle and the parallel component of the translation are invariant to the position of the axis, but the perpendicular component is not. This implies that the [\specialfonts {{\bsf L}}] tensor is unaffected by any assumptions about the position of the libration axes, whereas the [\specialfonts {{\bsf T}}] tensor depends on the assumptions made concerning the location of the axes.

The quadratic correlation between librational and translational motions can be allowed for by including in ([link] cross terms of the type [\langle D_{ik} t_{j}\rangle], or, with ([link], [\eqalignno{U_{ij} &= \langle D_{ik} D_{jl}\rangle r_{k} r_{l} + \langle D_{ik} t_{j} + D_{ji} t_{i}\rangle r_{k} + \langle t_{i} t_{j}\rangle \cr &= A_{ijkl} r_{k} r_{l} + B_{ijk} r_{k} + \langle t_{i} t_{j}\rangle,&(}] which leads to the explicit expressions such as [\eqalignno{U_{11} = \langle \delta r_{1}\rangle^{2} &= \langle \lambda_{3}^{2}\rangle r_{2}^{2} + \langle \lambda_{2}^{2}\rangle r_{3}^{2} - 2 \langle \lambda_{2} \lambda_{3}\rangle r_{2} r_{3}&\cr &\quad - 2 \langle \lambda_{3} t_{1}\rangle r_{2} - 2 \langle \lambda_{2} t_{1}\rangle r_{3} + \langle t_{1}^{2}\rangle,&\cr U_{12} = \langle \delta r_{1} \delta r_{2}\rangle &= -\langle \lambda_{3}^{2}\rangle r_{1} r_{2} + \langle \lambda_{1} \lambda_{3}\rangle r_{2} r_{3} + \langle \lambda_{2} \lambda_{3}\rangle r_{1} r_{3}&\cr &\quad - \langle \lambda_{1} \lambda_{2}\rangle r_{3}^{2} + \langle \lambda_{3} t_{1}\rangle r_{1} - \langle \lambda_{1} t_{1}\rangle r_{3} &\cr &\quad - \langle \lambda_{3} t_{2}\rangle + r_{2} \langle \lambda_{2} t_{2}\rangle r_{3} + \langle t_{1} t_{2}\rangle.&\cr&&(}]

The products of the type [\langle \lambda_{i}t_{j}\rangle] are the components of an additional tensor, [\specialfonts {{\bsf S}}], which unlike the tensors [\specialfonts {{\bsf T}}] and [ \specialfonts{{\bsf L}}] is unsymmetrical, since [\langle \lambda_{i}t_{j}\rangle] is different from [\langle \lambda_{j}t_{i}\rangle]. The terms involving elements of [\specialfonts {{\bsf S}}] may be grouped as [\langle \lambda_{3}t_{1}\rangle r_{1} - \langle \lambda_{3}t_{2}\rangle r_{2} + (\langle \lambda_{2}t_{2}\rangle - \langle \lambda_{1}t_{1}\rangle) r_{3} \eqno(] or [S_{31}r_{1} - S_{32}r_{2} + (S_{22} - S_{11}) r_{3}.] As the diagonal elements occur as differences in this expression, a constant may be added to each of the diagonal terms without changing the observational equations. In other words, the trace of [\specialfonts {{\bsf S}}] is indeterminate.

In terms of the [ \specialfonts{{\bsf L},{\bsf T}}] and [\specialfonts {{\bsf S}}] tensors, the observational equations are [U_{ij} = G_{ijkl}L_{kl} + H_{ijkl}S_{kl} + T_{ij}. \eqno(] The arrays [G_{ijkl}] and [H_{ijkl}] involve the atomic coordinates [(x,\ y,\ z) = (r_{1},\ r_{2},\ r_{3})], and are listed in Table[link] Equations ([link] for each of the atoms in the rigid body form the observational equations, from which the elements of [\specialfonts {{\bsf T},{\bsf L}}] and [\specialfonts {{\bsf S}}] can be derived by a linear least-squares procedure. One of the diagonal elements of [\specialfonts {{\bsf S}}] must be fixed in advance or some other suitable constraint applied because of the indeterminacy of [\specialfonts \hbox{Tr}({{\bsf S}})]. It is common practice to set [\specialfonts \hbox{Tr}({{\bsf S}})] equal to zero. There are thus eight elements of [\specialfonts {{\bsf S}}] to be determined, as well as the six each of [\specialfonts {{\bsf L}}] and [ \specialfonts{{\bsf T}}], for a total of 20 variables. A shift of origin leaves [\specialfonts {{\bsf L}}] invariant, but it intermixes [ \specialfonts{{\bsf T}}] and [\specialfonts {{\bsf S}}].

Table| top | pdf |
The arrays [G_{ijkl}] and [H_{ijkl}] to be used in the observational equations [U_{ij} = G_{ijkl} L_{kl} + H_{ijkl} S_{kl} + T_{ij}] [equation ([link]]


ij kl
11 0 [z^{2}] [y^{2}] [-2yz] 0 0
22 [z^{2}] 0 [x^{2}] 0 [-2xz] 0
33 [y^{2}] [x^{2}] 0 0 0 [-2xy]
23 yz 0 0 [-x^{2}] xy xz
31 0 xz 0 xy [-y^{2}] yz
12 0 0 xy xz yz [-z^{2}]


ij kl
11 0 0 0 0 [-2y] 0 0 0 2z
22 0 0 0 0 0 [-2z] 2x 0 0
33 0 0 0 [-2x] 0 0 0 2y 0
23 0 x x 0 0 y 0 z 0
31 y 0 y z 0 0 0 0 x
12 z z 0 0 x 0 y 0 0

If the origin is located at a centre of symmetry, for each atom at r with vibration tensor [\specialfonts {{\bsf U}}^{n}] there will be an equivalent atom at −r with the same vibration tensor. When the observational equations for these two atoms are added, the terms involving elements of [\specialfonts {{\bsf S}}] disappear since they are linear in the components of r. The other terms, involving elements of the [\specialfonts {{\bsf T}}] and [ \specialfonts{{\bsf L}}] tensors, are simply doubled, like the [\specialfonts {{\bsf U}}^{n}] components.

The physical meaning of the [\specialfonts {{\bsf T}}] and [\specialfonts {{\bsf L}}] tensor elements is as follows. [T_{ij}l_{i}l_{j}] is the mean-square amplitude of translational vibration in the direction of the unit vector l with components [l_{1},\ l_{2},\ l_{3}] along the Cartesian axes and [L_{ij} l_{i} l_{j}] is the mean-square amplitude of libration about an axis in this direction. The quantity [S_{ij} l_{i} l_{j}] represents the mean correlation between libration about the axis l and translation parallel to this axis. This quantity, like [T_{ij}l_{i}l_{j}], depends on the choice of origin, although the sum of the two quantities is independent of the origin.

The non-symmetrical tensor [\specialfonts {{\bsf S}}] can be written as the sum of a symmetric tensor with elements [{S^{S}_{ij}} = (S_{ij} + S_{ji})/2] and a skew-symmetric tensor with elements [{S^{A}_{ij}} = (S_{ij} - S_{ji})/2]. Expressed in terms of principal axes, [\specialfonts {{\bsf S}}^{S}] consists of three principal screw correlations [\langle \lambda_{I}t_{I}\rangle]. Positive and negative screw correlations correspond to opposite senses of helicity. Since an arbitrary constant may be added to all three correlation terms, only the differences between them can be determined from the data.

The skew-symmetric part [\specialfonts {{\bsf S}}^{A}] is equivalent to a vector [({\boldlambda} \times {\bf t})/2] with components [({\boldlambda} \times {\bf t})_{i}/2 = (\lambda_{j}t_{k} - \lambda_{k}t_{j})/2], involving correlations between a libration and a perpendicular translation. The components of [\specialfonts {{\bsf S}}^{A}] can be reduced to zero, and [\specialfonts {{\bsf S}}] made symmetric, by a change of origin. It can be shown that the origin shift that symmetrizes [\specialfonts {{\bsf S}}] also minimizes the trace of [\specialfonts {{\bsf T}}]. In terms of the coordinate system based on the principal axes of [\specialfonts {{\bsf L}}], the required origin shifts [\widehat{\rho}_{i}] are [\widehat{\rho}_{1} = {\widehat{S}_{23} - \widehat{S}_{32} \over \widehat{L}_{22} + \widehat{L}_{33}} \quad \widehat{\rho}_{2} = {\widehat{S}_{31} - \widehat{S}_{13} \over \widehat{L}_{11} + \widehat{L}_{33}} \quad \widehat{\rho}_{3} = {\widehat{S}_{12} - \widehat{S}_{21} \over \widehat{L}_{11} + \widehat{L}_{22}}, \eqno(] in which the carets indicate quantities referred to the principal axis system.

The description of the averaged motion can be simplified further by shifting to three generally non-intersecting libration axes, one each for each principal axis of [\specialfonts {{\bsf L}}]. Shifts of the [\specialfonts {{\bsf L}}_{1}] axis in the [\specialfonts {{\bsf L}}_{2}] and [\specialfonts {{\bsf L}}_{3}] directions by [^{1}\widehat{\rho}_{2} = - \widehat{S}_{13}/\widehat{L}_{11} \hbox{ and } ^{1}\widehat{\rho}_{3} = \widehat{S}_{12}/\widehat{L}_{11}, \eqno(] respectively, annihilate the [S_{12}] and [S_{13}] terms of the symmetrized [\specialfonts {{\bsf S}}] tensor and simultaneously effect a further reduction in [\specialfonts \hbox{Tr}({{\bsf T}})] (the presuperscript denotes the axis that is shifted, the subscript the direction of the shift component). Analogous equations for displacements of the [\specialfonts {{\bsf L}}_{2}] and [\specialfonts {{\bsf L}}_{3}] axes are obtained by permutation of the indices. If all three axes are appropriately displaced, only the diagonal terms of [\specialfonts {{\bsf S}}] remain. Referred to the principal axes of [\specialfonts {{\bsf L}}], they represent screw correlations along these axes and are independent of origin shifts.

The elements of the reduced [\specialfonts {{\bsf T}}] are [\eqalignno{^{r}T_{II} &= \widehat{T}_{II} - {\textstyle\sum\limits_{K \neq I}} (\widehat{S}_{KI})^{2}/\widehat{L}_{KK}\cr ^{r}T_{IJ} &= \widehat{T}_{IJ} - {\textstyle\sum\limits_{K}} \widehat{S}_{KI} \widehat{S}_{KJ}/\widehat{L}_{KK},\quad J \neq I. &(}]

The resulting description of the average rigid-body motion is in terms of six independently distributed instantaneous motions – three screw librations about non-intersecting axes (with screw pitches given by [\widehat{S}_{11}/\widehat{L}_{11}] etc.) and three translations. The parameter set consists of three libration and three translation amplitudes, six angles of orientation for the principal axes of [\specialfonts {{\bsf L}}] and [\specialfonts {{\bsf T}}], six coordinates of axis displacement, and three screw pitches, one of which has to be chosen arbitrarily, again for a total of 20 variables.

Since diagonal elements of [\specialfonts {{\bsf S}}] enter into the expression for [^{r}T_{IJ}], the indeterminacy of [\specialfonts \hbox{Tr}({{\bsf S}})] introduces a corresponding indeterminacy in [\specialfonts^{r} {{\bsf T}}]. The constraint [\specialfonts \hbox{Tr}({{\bsf S}}) = 0] is unaffected by the various rotations and translations of the coordinate systems used in the course of the analysis.

1.2.12. Treatment of anharmonicity

| top | pdf |

The probability distribution ([link] [link] is valid in the case of rectilinear harmonic motion. If the deviations from Gaussian shape are not too large, distributions may be used which are expansions with the Gaussian distribution as the leading term. Three such distributions are discussed in the following sections. The Gram–Charlier expansion

| top | pdf |

The three-dimensional Gram–Charlier expansion, introduced into thermal-motion treatment by Johnson & Levy (1974)[link], is an expansion of a function in terms of the zero and higher derivatives of a normal distribution (Kendall & Stuart, 1958[link]). If [D_{j}] is the operator [\hbox{d/d}u\hskip 2pt^{j}], [\eqalignno{P({\bf u}) &= [1 - c\hskip 2pt^{j}D_{j} + {1 \over 2!} c\hskip 2pt^{jk}D_{j}D_{k} - {1 \over 3!} c\hskip 2pt^{jkl}D_{j}D_{k}D_{l} + \ldots\cr &\quad + (-1)^{r} {c^{\alpha_{1}} \ldots c^{\alpha_{r}} \over r!} D_{\alpha_{1}} D_{\alpha_{r}}] P_{0}({\bf u}), &(}] where [P_{0}({\bf u})] is the harmonic distribution, [\alpha_{1} = 1, 2] or 3, and the operator [D_{\alpha_{1}} \ldots D_{\alpha_{r}}] is the rth partial derivative [\partial^{r}/(\partial u^{\alpha 1} \ldots \partial u^{\alpha r})]. Summation is again implied over repeated indices.

The differential operators D may be eliminated by the use of three-dimensional Hermite polynomials [H_{\alpha_{1}\ldots \alpha_{2}}] defined, by analogy with the one-dimensional Hermite polynomials, by the expression [{D_{\alpha_{1}}\ldots D_{\alpha_{r}} \exp (- {\textstyle {1 \over 2}} \sigma_{jk}^{-1} u\hskip 2pt^{j}u^{k}) = (-1)^{r} H_{\alpha_{1}\ldots \alpha_{r}} ({\bf u}) \exp (- {\textstyle{1 \over 2}} \sigma_{jk}^{-1} u\hskip 2pt^{j}u^{k}),} \eqno(] which gives [\eqalignno{P({\bf u}) &= \left[1 + {1 \over 3!} c\hskip 2pt^{jkl}H_{jkl}({\bf u}) + {1 \over 4!} c\hskip 2pt^{jklm}H_{jklm} ({\bf u}) + {1 \over 5!} c\hskip 2pt^{jklmn}H_{jklmn} ({\bf u})\right.\cr &\quad\left. + {1 \over 6!} c\hskip 2pt^{jklmnp}H_{jklmnp} ({\bf u}) + \ldots\right] P_{0}({\bf u}), &(}] where the first and second terms have been omitted since they are equivalent to a shift of the mean and a modification of the harmonic term only. The permutations of [j, k, l \ldots] here, and in the following sections, include all combinations which produce different terms.

The coefficients c, defined by ([link] and ([link], are known as the quasimoments of the frequency function [P(\bf u)] (Kutznetsov et al., 1960[link]). They are related in a simple manner to the moments of the function (Kendall & Stuart, 1958[link]) and are invariant to permutation of indices. There are 10, 15, 21 and 28 components of c for orders 3, 4, 5 and 6, respectively. The multivariate Hermite polynomials are functions of the elements of [\sigma_{jk}^{-1}] and of [u^{k}], and are given in Table[link] for orders [\leq 6] (IT IV, 1974[link]; Zucker & Schulz, 1982[link]).

Table| top | pdf |
Some Hermite polynomials (Johnson & Levy, 1974[link]; Zucker & Schulz, 1982[link])

H(u) = 1
Hj(u) = wj
Hjk(u) = wjwkpjk
Hjkl(u) = wjwkwl − (wjpkl + wkplj + wlpjk) = wjwkwl3w( jpkl)
Hjklm(u) = wjwkwlwm6w( jwkplm) + 3pj( kplm)
Hjklmn(u) = wjwkwlwmwn10w( lwmwnpjk) + 15w( npjkplm)
Hjklmnp(u) = wjwkwlwmwnwp − 15w( jwkwlwmpjk) + 45w( jwkplmpnp) − 15pj( kplmpnp)
where [w_{j}\equiv p{_{jk}}u^{k} \hbox{ and } p_{jk}] are the elements of [\sigma^{-1}], defined in expression ([link] [link]. Indices between brackets indicate that the term is to be averaged over all permutations which produce distinct terms, keeping in mind that [p_{jk} = p_{kj} \hbox{ and } w_{j}w_{k} = w_{k}w_{j}] as illustrated for [H_{jkl}].

The Fourier transform of ([link] is given by [\eqalignno{T({\bf H}) &= \left[1 - {4 \over 3} \pi^{3}ic\hskip 2pt^{jkl}h_{j}h_{k}h_{l} + {2 \over 3} \pi^{4}c\hskip 2pt^{jklm}h_{j}h_{k}h_{l}h_{m}\right.\cr &\quad + {4 \over 15} \pi^{5}ic\hskip 2pt^{jklmn}h_{j}h_{k}h_{l}h_{m}h_{n} &\cr &\quad\left. - {4 \over 45} \pi^{6}c\hskip 2pt^{jklmnp}h_{j}h_{k}h_{l}h_{m}h_{n}h_{p} + \ldots\right] T_{0}({\bf H}), &(}] where [T_{0}({\bf H})] is the harmonic temperature factor. [T({\bf H})] is a power-series expansion about the harmonic temperature factor, with even and odd terms, respectively, real and imaginary. The cumulant expansion

| top | pdf |

A second statistical expansion which has been used to describe the atomic probability distribution is that of Edgeworth (Kendall & Stuart, 1958[link]; Johnson, 1969[link]). It expresses the function [P({\bf u})] as [\eqalignno{P({\bf u}) &= \exp \left(\kappa\hskip 2pt^{j}D_{j} + {1 \over 2!} \kappa\hskip 2pt^{jk}D_{j}D_{k} - {1 \over 3!} \kappa\hskip 2pt^{jkl}D_{j}D_{k}D_{l}\right.\cr &\quad\left. + {1 \over 4!} \kappa\hskip 2pt^{jklm}D_{j}D_{k}D_{l}D_{m} - \ldots\right)P_{0}({\bf u}). &(}]

Like the moments μ of a distribution, the cumulants κ are descriptive constants. They are related to each other (in the one-dimensional case) by the identity [{\exp\left\{\kappa_{1}t + {\kappa_{2}t^{2} \over 2!} + \ldots {\kappa_{r}t^{r} \over r!} + \ldots\right\} = 1 + \mu_{1}t + {\mu_{2}t^{2} \over 2!} + \ldots + {\mu_{r}t^{r} \over r!}.} \eqno(] When it is substituted for t, ([link] is the characteristic function, or Fourier transform of [P(t)] (Kendall & Stuart, 1958[link]).

The first two terms in the exponent of ([link] can be omitted if the expansion is around the equilibrium position and the harmonic term is properly described by [P_{0}({\bf u})].

The Fourier transform of ([link] is, by analogy with the left-hand part of ([link] (with t replaced by [2\pi ih]), [\eqalignno{T({\bf H}) &= \exp \left[{(2\pi i)^{3} \over 3!} \kappa\hskip 2pt^{jkl}h_{j}h_{k}h_{l} + {(2\pi i)^{4} \over 4!} \kappa\hskip 2pt^{jklm}h_{j}h_{k}h_{l}h_{m} + \ldots \right] T_{0}({\bf H})\cr &= \exp \left[- {4 \over 3} \pi^{3}i \kappa\hskip 2pt^{jkl}h_{j}h_{k}h_{l} + {2 \over 3} \pi^{4} \kappa\hskip 2pt^{jklm}h_{j}h_{k}h_{l}h_{m} + \ldots \right] T_{0}({\bf H}),\cr &&(}] where the first two terms have been omitted. Expression ([link] is similar to ([link] except that the entire series is in the exponent. Following Schwarzenbach (1986)[link], ([link] can be developed in a Taylor series, which gives [\eqalignno{T({\bf H}) &= \left\{1 + {(2\pi i)^{3} \over 3!} \kappa\hskip 2pt^{jkl}h_{j}h_{k}h_{l} + {(2\pi i)^{4} \over 4!} \kappa\hskip 2pt^{jklm}h_{j}h_{k}h_{l}h_{m} + \ldots\right.\cr &\quad + {(2\pi i)^{6} \over 6!} \left[\vphantom{{\textstyle\sum\limits_{1}^{2}}}\kappa\hskip 2pt^{jklmp} + {6! \over 2!(3!)^{2}} \kappa\hskip 2pt^{jkl}\kappa^{mnp}\right] h_{j}h_{k}h_{l}h_{m}h_{n}h_{p}\cr &\quad + \left.\vphantom{{\textstyle\sum\limits_{1}^{2}}}\hbox{higher-order terms}\right\} T_{0}({\bf H}). &(}]

This formulation, which is sometimes called the Edgeworth approximation (Zucker & Schulz, 1982[link]), clearly shows the relation to the Gram–Charlier expansion ([link], and corresponds to the probability distribution [analogous to ([link]] [\eqalignno{P({\bf u}) &= P_{0}({\bf u}) \left\{1 + {1 \over 3!} \kappa\hskip 2pt^{jkl}H_{jkl}({\bf u}) + {1 \over 4!} \kappa\hskip 2pt^{jklm}H_{jklm}({\bf u}) + \ldots\right.\cr &{\hbox to 7.25pt{}} + {1 \over 6!} \left[\vphantom{{\textstyle\sum\limits_{1}^{2}}}\kappa\hskip 2pt^{jklmnp} + 10 \kappa\hskip 2pt^{jkl} \kappa^{mnp} \vphantom{{\textstyle\sum\limits_{1}^{2}}}\right] H_{jklmnp}\cr &\quad \left.+ \hbox{ higher-order terms}\vphantom{{\textstyle\sum\limits_{1}^{2}}}\right\}. &(}]

The relation between the cumulants [\kappa\hskip 2pt^{jkl}] and the quasimoments [c\hskip 2pt^{jkl}] are apparent from comparison of ([link] and ([link]: [\eqalignno{c\hskip 2pt^{jkl} &= \kappa\hskip 2pt^{jkl}\cr c\hskip 2pt^{jklm} &= \kappa\hskip 2pt^{jklm}\cr c\hskip 2pt^{jklmn} &= \kappa\hskip 2pt^{jklmn}\cr c\hskip 2pt^{jklmnp} &= \kappa\hskip 2pt^{jklmnp} + 10\kappa\hskip 2pt^{jkl} \kappa^{mnp}. &(}]

The sixth- and higher-order cumulants and quasimoments differ. Thus the third-order cumulant [\kappa\hskip 2pt^{jkl}] contributes not only to the coefficient of [H_{jkl}], but also to higher-order terms of the probability distribution function. This is also the case for cumulants of higher orders. It implies that for a finite truncation of ([link], the probability distribution cannot be represented by a finite number of terms. This is a serious difficulty when a probability distribution is to be derived from an experimental temperature factor of the cumulant type. The one-particle potential (OPP) model

| top | pdf |

When an atom is considered as an independent oscillator vibrating in a potential well [V({\bf u})], its distribution may be described by Boltzmann statistics. [P({\bf u}) = N \exp \{- V({\bf u})/kT\}, \eqno(] with N, the normalization constant, defined by [{\textstyle\int} P({\bf u})\ \hbox{d}{\bf u} = 1]. The classical expression ([link] is valid in the high-temperature limit for which [kT \gg {\it V}({\bf u})].

Following Dawson (1967)[link] and Willis (1969)[link], the potential function may be expanded in terms of increasing order of products of the contravariant displacement coordinates: [{V = V_{0} + \alpha_{j}u\hskip 2pt^{j} + \beta_{jk}u\hskip 2pt^{j}u^{k} + \gamma_{jkl}u\hskip 2pt^{j}u^{k}u^{l} + \delta_{jklm}u\hskip 2pt^{j}u^{k}u^{l}u^{m} + \ldots .} \eqno(] The equilibrium condition gives [\alpha_{j} = 0]. Substitution into ([link] leads to an expression which may be simplified by the assumption that the leading term is the harmonic component represented by [\beta_{jk}]: [\eqalignno{P({\bf u}) &= N \exp \{- \beta^{'}_{jk}u\hskip 2pt^{j}u^{k}\}\cr &\quad \times \{1 - \gamma^{'}_{jkl}u\hskip 2pt^{j}u^{k}u^{l} - \delta^{'}_{jklm}u\hskip 2pt^{j}u^{k}u^{l}u^{m} - \ldots \}, &(}] in which [\beta^{'} = \beta / kT] etc. and the normalization factor N depends on the level of truncation.

The probability distribution is related to the spherical harmonic expansion. The ten products of the displacement parameters [u\hskip 2pt^{j}u^{k}u^{l}], for example, are linear combinations of the seven octapoles [(l = 3)] and three dipoles [(l = 1)] (Coppens, 1980[link]). The thermal probability distribution and the aspherical atom description can be separated only because the latter is essentially confined to the valence shell, while the former applies to all electrons which follow the nuclear motion in the atomic scattering model.

The Fourier transform of the OPP distribution, in a general coordinate system, is (Johnson, 1970a[link]; Scheringer, 1985a[link]) [\eqalignno{T({\bf H}) &= T_{0} ({\bf H}) \left[1 - {4 \over 3} \pi^{3} i \gamma^{'}_{jkl} G\hskip 2pt^{jkl} ({\bf H}) + {2 \over 3} \pi^{4} \delta^{'}_{jklm} G\hskip 2pt^{jklm} ({\bf H})\right.\cr &\quad\left. + {4 \over 15} \pi^{5} i \varepsilon^{'}_{jklmn} G\hskip 2pt^{jklmn} ({\bf H}) - {4 \over 45} \pi^{6} i \varphi^{'}_{jklmnp} G\hskip 2pt^{jklmnp} ({\bf H}) \ldots\right],\cr& &(}] where [T_{0}] is the harmonic temperature factor and G represents the Hermite polynomials in reciprocal space.

If the OPP temperature factor is expanded in the coordinate system which diagonalizes [\beta_{jk}], simpler expressions are obtained in which the Hermite polynomials are replaced by products of the displacement coordinates [u\hskip 2pt^{j}] (Dawson et al., 1967[link]; Coppens, 1980[link]; Tanaka & Marumo, 1983[link]). Relative merits of the three expansions

| top | pdf |

The relative merits of the Gram–Charlier and Edgeworth expansions have been discussed by Zucker & Schulz (1982)[link], Kuhs (1983)[link], and by Scheringer (1985b)[link]. In general, the Gram–Charlier expression is found to be preferable because it gives a better fit in the cases tested, and because its truncation is equivalent in real and reciprocal space. The latter is also true for the one-particle potential model, which is mathematically related to the Gram–Charlier expansion by the interchange of the real- and reciprocal-space expressions. The terms of the OPP model have a specific physical meaning. The model allows prediction of the temperature dependence of the temperature factor (Willis, 1969[link]; Coppens, 1980[link]), provided the potential function itself can be assumed to be temperature independent.

It has recently been shown that the Edgeworth expansion ([link] always has negative regions (Scheringer, 1985b[link]). This implies that it is not a realistic description of a vibrating atom.

1.2.13. The generalized structure factor

| top | pdf |

In the generalized structure-factor formalism developed by Dawson (1975)[link], the complex nature of both the atomic scattering factor and the generalized temperature factor are taken into account. We write for the atomic scattering factor: [\eqalignno{f_{j} ({\bf H}) &= f_{j, \, c} ({\bf H}) + if_{j, \, a} ({\bf H}) + f^{'}_{j} + if^{''}_{j} &(\cr T_{j} ({\bf H}) &= T_{j, \, c} ({\bf H}) + iT_{j, \, a} ({\bf H}) &(}] and [F({\bf H}) = A({\bf H}) + iB({\bf H}), \eqno(] where the subscripts c and a refer to the centrosymmetric and acentric components, respectively. Substitution in ([link] [link] [link] [link] gives for the real and imaginary components A and B of [F({\bf H})] [\eqalignno{A ({\bf H}) &= {\textstyle\sum\limits_{j}} (f_{j, \, c} + f^{'}_{j}) [\cos (2\pi {\bf H} \cdot {\bf r}_{j}) T_{c} - \sin (2\pi {\bf H} \cdot {\bf r}_{j}) T_{a}]\cr &\quad - (f_{j, \, a} + f^{''}_{j}) [\cos (2\pi {\bf H} \cdot {\bf r}_{j}) T_{a} + \sin (2\pi {\bf H} \cdot {\bf r}_{j}) T_{c}]\cr&&(}] and [\eqalignno{B ({\bf H}) &= {\textstyle\sum\limits_{j}} (f_{j, \, c} + f^{'}_{j}) [\cos (2\pi {\bf H} \cdot {\bf r}_{j}) T_{a} + \sin (2\pi {\bf H} \cdot {\bf r}_{j}) T_{c}]\cr &\quad + (f_{j, \, a} + f^{''}_{j}) [\cos (2\pi {\bf H} \cdot {\bf r}_{j}) T_{c} - \sin (2\pi {\bf H} \cdot {\bf r}_{j}) T_{a}]\cr& &(}] (McIntyre et al., 1980[link]; Dawson, 1967[link]).

Expressions ([link] [link] illustrate the relation between valence-density anisotropy and anisotropy of thermal motion.

1.2.14. Conclusion

| top | pdf |

This chapter summarizes mathematical developments of the structure-factor formalism. The introduction of atomic asphericity into the formalism and the treatment of thermal motion are interlinked. It is important that the complexities of the thermal probability distribution function can often be reduced by very low temperature experimentation. Results obtained with the multipole formalism for atomic asphericity can be used to derive physical properties and d-orbital populations of transition-metal atoms (IT C, 2004[link]). In such applications, the deconvolution of the charge density and the thermal vibrations is essential. This deconvolution is dependent on the adequacy of the models summarized here.


The author would like to thank several of his colleagues who gave invaluable criticism of earlier versions of this manuscript. Corrections and additions were made following comments by P. J. Becker, D. Feil, N. K. Hansen, G. McIntyre, E. N. Maslen, S. Ohba, C. Scheringer and D. Schwarzenbach. Z. Su contributed to the revised version of the manuscript. Support of this work by the US National Science Foundation (CHE8711736 and CHE9317770) is gratefully acknowledged.


Arfken, G. (1970). Mathematical models for physicists, 2nd ed. New York, London: Academic Press.
Avery, J. & Ørmen, P.-J. (1979). Generalized scattering factors and generalized Fourier transforms. Acta Cryst. A35, 849–851.
Avery, J. & Watson, K. J. (1977). Generalized X-ray scattering factors. Simple closed-form expressions for the one-centre case with Slater-type orbitals. Acta Cryst. A33, 679–680.
Bentley, J. & Stewart, R. F. (1973). Two-centre calculations for X-ray scattering. J. Comput. Phys. 11, 127–145.
Born, M. (1926). Quantenmechanik der Stoszvorgänge. Z. Phys. 38, 803.
Clementi, E. & Raimondi, D. L. (1963). Atomic screening constants from SCF functions. J. Chem. Phys. 38, 2686–2689.
Clementi, E. & Roetti, C. (1974). Roothaan–Hartree–Fock atomic wavefunctions. At. Data Nucl. Data Tables, 14, 177–478.
Cohen-Tannoudji, C., Diu, B. & Laloe, F. (1977). Quantum mechanics. New York: John Wiley and Paris: Hermann.
Condon, E. V. & Shortley, G. H. (1957). The theory of atomic spectra. London, New York: Cambridge University Press.
Coppens, P. (1980). Thermal smearing and chemical bonding. In Electron and magnetization densities in molecules and solids, edited by P. J. Becker, pp. 521–544. New York: Plenum.
Coppens, P., Guru Row, T. N., Leung, P., Stevens, E. D., Becker, P. J. & Yang, Y. W. (1979). Net atomic charges and molecular dipole moments from spherical-atom X-ray refinements, and the relation between atomic charges and shape. Acta Cryst. A35, 63–72.
Coulson, C. A. (1961). Valence. Oxford University Press.
Cruickshank, D. W. J. (1956). The analysis of the anisotropic thermal motion of molecules in crystals. Acta Cryst. 9, 754–756.
Dawson, B. (1967). A general structure factor formalism for interpreting accurate X-ray and neutron diffraction data. Proc. R. Soc. London Ser. A, 248, 235–288.
Dawson, B. (1975). Studies of atomic charge density by X-ray and neutron diffraction – a perspective. In Advances in structure research by diffraction methods. Vol. 6, edited by W. Hoppe & R. Mason. Oxford: Pergamon Press.
Dawson, B., Hurley, A. C. & Maslen, V. W. (1967). Anharmonic vibration in fluorite-structures. Proc. R. Soc. London Ser. A, 298, 289–306.
Dunitz, J. D. (1979). X-ray analysis and the structure of organic molecules. Ithaca and London: Cornell University Press.
Feil, D. (1977). Diffraction physics. Isr. J. Chem. 16, 103–110.
Hansen, N. K. & Coppens, P. (1978). Testing aspherical atom refinements on small-molecule data sets. Acta Cryst. A34, 909–921.
Hehre, W. J., Ditchfield, R., Stewart, R. F. & Pople, J. A. (1970). Self-consistent molecular orbital methods. IV. Use of Gaussian expansions of Slater-type orbitals. Extension to second-row molecules. J. Chem. Phys. 52, 2769–2773.
Hehre, W. J., Stewart, R. F. & Pople, J. A. (1969). Self-consistent molecular orbital methods. I. Use of Gaussian expansions of Slater-type atomic orbitals. J. Chem. Phys. 51, 2657–2664.
Hirshfeld, F. L. (1977). A deformation density refinement program. Isr. J. Chem. 16, 226–229.
International Tables for Crystallography (2004). Vol. C. Mathematical, physical and chemical tables, edited by E. Prince. Dordrecht: Kluwer Academic Publishers.
International Tables for X-ray Crystallography (1974). Vol. IV. Birmingham: Kynoch Press. (Present distributor Kluwer Academic Publishers, Dordrecht.)
James, R. W. (1982). The optical principles of the diffraction of X-rays. Woodbridge, Connecticut: Oxbow Press.
Johnson, C. K. (1969). Addition of higher cumulants to the crystallographic structure-factor equation: a generalized treatment for thermal-motion effects. Acta Cryst. A25, 187–194.
Johnson, C. K. (1970a). Series expansion models for thermal motion. ACA Program and Abstracts, 1970 Winter Meeting, Tulane University, p. 60.
Johnson, C. K. (1970b). An introduction to thermal-motion analysis. In Crystallographic computing, edited by F. R. Ahmed, S. R. Hall & C. P. Huber, pp. 207–219. Copenhagen: Munksgaard.
Johnson, C. K. & Levy, H. A. (1974). Thermal motion analysis using Bragg diffraction data. In International tables for X-ray crystallography (1974), Vol. IV, pp. 311–336. Birmingham: Kynoch Press. (Present distributor Kluwer Academic Publishers, Dordrecht.)
Kara, M. & Kurki-Suonio, K. (1981). Symmetrized multipole analysis of orientational distributions. Acta Cryst. A37, 201–210.
Kendall, M. G. & Stuart, A. (1958). The advanced theory of statistics. London: Griffin.
Kuhs, W. F. (1983). Statistical description of multimodal atomic probability structures. Acta Cryst. A39, 148–158.
Kurki-Suonio, K. (1977). Symmetry and its implications. Isr. J. Chem. 16, 115–123.
Kutznetsov, P. I., Stratonovich, R. L. & Tikhonov, V. I. (1960). Theory Probab. Its Appl. (USSR), 5, 80–97.
Lipson, H. & Cochran, W. (1966). The determination of crystal structures. London: Bell.
McIntyre, G. J., Moss, G. & Barnea, Z. (1980). Anharmonic temperature factors of zinc selenide determined by X-ray diffraction from an extended-face crystal. Acta Cryst. A36, 482–490.
Materlik, G., Sparks, C. J. & Fischer, K. (1994). Resonant anomalous X-ray scattering. Theory and applications. Amsterdam: North-Holland.
Paturle, A. & Coppens, P. (1988). Normalization factors for spherical harmonic density functions. Acta Cryst. A44, 6–7.
Ruedenberg, K. (1962). The nature of the chemical bond. Phys. Rev. 34, 326–376.
Scheringer, C. (1985a). A general expression for the anharmonic temperature factor in the isolated-atom-potential approach. Acta Cryst. A41, 73–79.
Scheringer, C. (1985b). A deficiency of the cumulant expansion of the anharmonic temperature factor. Acta Cryst. A41, 79–81.
Schomaker, V. & Trueblood, K. N. (1968). On the rigid-body motion of molecules in crystals. Acta Cryst. B24, 63–76.
Schwarzenbach, D. (1986). Private communication.
Squires, G. L. (1978). Introduction to the theory of thermal neutron scattering. Cambridge University Press.
Stewart, R. F. (1969a). Generalized X-ray scattering factors. J. Chem. Phys. 51, 4569–4577.
Stewart, R. F. (1969b). Small Gaussian expansions of atomic orbitals. J. Chem. Phys. 50, 2485–2495.
Stewart, R. F. (1970). Small Gaussian expansions of Slater-type orbitals. J. Chem. Phys. 52, 431–438.
Stewart, R. F. (1980). Electron and magnetization densities in molecules and solids, edited by P. J. Becker, pp. 439–442. New York: Plenum.
Stewart, R. F. & Hehre, W. J. (1970). Small Gaussian expansions of atomic orbitals: second-row atoms. J. Chem. Phys. 52, 5243–5247.
Su, Z. & Coppens, P. (1990). Closed-form expressions for Fourier– Bessel transforms of Slater-type functions. J. Appl. Cryst. 23, 71–73.
Su, Z. & Coppens, P. (1994). Normalization factors for Kubic harmonic density functions. Acta Cryst. A50, 408–409.
Tanaka, K. & Marumo, F. (1983). Willis formalism of anharmonic temperature factors for a general potential and its application in the least-squares method. Acta Cryst. A39, 631–641.
Von der Lage, F. C. & Bethe, H. A. (1947). A method for obtaining electronic functions and eigenvalues in solids with an application to sodium. Phys. Rev. 71, 612–622.
Waller, I. & Hartree, D. R. (1929). Intensity of total scattering X-rays. Proc. R. Soc. London Ser. A, 124, 119–142.
Weiss, R. J. & Freeman, A. J. (1959). X-ray and neutron scattering for electrons in a crystalline field and the determination of outer electron configurations in iron and nickel. J. Phys. Chem. Solids, 10, 147–161.
Willis, B. T. M. (1969). Lattice vibrations and the accurate determination of structure factors for the elastic scattering of X-rays and neutrons. Acta Cryst. A25, 277–300.
Zucker, U. H. & Schulz, H. (1982). Statistical approaches for the treatment of anharmonic motion in crystals. I. A comparison of the most frequently used formalisms of anharmonic thermal vibrations. Acta Cryst. A38, 563–568.

to end of page
to top of page