Tables for
Volume B
Reciprocal space
Edited by U. Shmueli

International Tables for Crystallography (2006). Vol. B, ch. 4.2, pp. 410-420   | 1 | 2 |

Section 4.2.3. General treatment

H. Jagodzinskia and F. Freyb

aInstitut für Kristallographie und Mineralogie, Universität, Theresienstrasse 41, D-8000 München 2, Germany, and  bInstitut für Kristallographie und Mineralogie, Universität, Theresienstrasse 41, D-8000 München 2, Germany

4.2.3. General treatment

| top | pdf | Qualitative interpretation of diffuse scattering

| top | pdf |

Any structure analysis of disordered structures should start with a qualitative interpretation of diffuse scattering. This problem may be facilitated with the aid of Fourier transforms and their algebraic operations (see, e.g., Patterson, 1959[link]). For simplicity the following modified notation is used in this section: functions in real space are represented by small letters, e.g. a(r), b(r), … except for F(r) and P(r) which are used as general symbols for a structure and the Patterson function, respectively; functions in reciprocal space are given by capital letters A(H), B(H); r and H are general vectors in real and reciprocal space, respectively, [Hz + Ky + Lz] is the scalar product H · r; dr and dH indicate integrations in three dimensions in real and reciprocal space, respectively. Even for X-rays the electron density [\rho ({\bf r})] will generally be replaced by the scattering potential a(r). Consequently, anomalous contributions to scattering may be included if complex functions a(r) are admitted. In the neutron case a(r) refers to a quasi-potential. Using this notation we obtain the scattered amplitude and phase [A({\bf H}) \exp i\varphi] [\eqalignno{A({\bf H}) &= \textstyle\int\limits_{{\bf r}} a({\bf r}) \exp \{2 \pi i {\bf H} \cdot {\bf r}\} \;\hbox{d}{\bf r} &(\cr a({\bf r}) &= \textstyle\int\limits_{{\bf H}} A({\bf H}) \exp \{-2 \pi i {\bf H} \cdot {\bf r}\} \;\hbox{d}{\bf H} &(}] (constant factors are omitted).

a(r) and A(H) are reversibly and uniquely determined by Fourier transformation. Consequently equations ([link] may simply be replaced by [a({\bf r}) \leftrightarrow A({\bf H})], where the double-headed arrow represents the two integrations given by ([link] and means: A(H) is the Fourier transform of a(r), and vice versa. The following relations may easily be derived from ([link]: [a({\bf r}) + b({\bf r}) \leftrightarrow A({\bf H}) + B({\bf H})\quad (\hbox{law of addition}) \eqno(] [\beta a({\bf r}) \leftrightarrow \beta A({\bf H})\quad (\hbox{law of scalar multiplication}) \eqno(] ([\beta ] = scalar quantity).

On the other hand, the multiplication of two functions does not yield a relation of similar symmetrical simplicity: [\eqalignno{a({\bf r}) b({\bf r}) &\leftrightarrow \textstyle\int A({\bf H}') B({\bf H} - {\bf H}') \;\hbox{d}{\bf H}' &\cr &= A({\bf H}) * B({\bf H}) &(\cr a({\bf r}) * b({\bf r}) &= \textstyle\int a({\bf r}') b({\bf r} - {\bf r}') \;\hbox{d}{\bf r}' &\cr &\leftrightarrow A({\bf H}) B({\bf H}) &(}] (laws of convolution and multiplication).

Since [a({\bf r}) b({\bf r}) = b({\bf r}) a({\bf r})]: [\textstyle\int A({\bf H}') B({\bf H} - {\bf H}') \;\hbox{d}{\bf H}' = \textstyle\int B({\bf H}') A({\bf H} - {\bf H}') \;\hbox{d}{\bf H}'] and vice versa. The convolution operation is commutative in either space.

For simplicity the symbol [a({\bf r}) * b({\bf r})] instead of the complete convolution integral is used. The distribution law [a(b + c) = ab + ac] is valid for the convolution as well: [a({\bf r}) * [b({\bf r}) + c({\bf r})] = a({\bf r}) * b({\bf r}) + a({\bf r}) * c({\bf r}). \eqno(] The associative law of multiplication does not hold if mixed products (convolution and multiplication) are used: [a({\bf r}) * [b({\bf r}) c({\bf r})] \neq [a({\bf r}) * b({\bf r})] c({\bf r}). \eqno(] From equations ([link] one has: [\eqalign{a({\bf r} - {\bf r}_{0}) &\leftrightarrow A({\bf H}) \exp (2 \pi i {\bf H} \cdot {\bf r}_{0})\cr A({\bf H} - {\bf H}_{0}) &\leftrightarrow a({\bf r}) \exp (-2 \pi i {\bf H}_{0} \cdot {\bf r}) \cr} \eqno(] (law of displacements).

Since symmetry operations are well known to crystallographers in reciprocal space also, the law of inversion is mentioned here only: [a(-{\bf r}) \leftrightarrow A(-{\bf H}). \eqno(] Consequently, if [a({\bf r}) = a(-{\bf r})], then [A({\bf H}) = A(-{\bf H})]. In order to calculate the intensity the complex conjugate [A^{+}({\bf H})] is needed: [\eqalignno{a^{+}({\bf r}) &\leftrightarrow A^{+}(-{\bf H}) &(\cr a^{+}(-{\bf r}) &\leftrightarrow A^{+}({\bf H}). &(}] Equations ([link] yield the relationship [A^{+}(-{\bf H}) = A({\bf H})] (`Friedel's law') if a(r) is a real function. The multiplication of a function with its conjugate is given by: [a({\bf r}) * a^{+}(-{\bf r}) \leftrightarrow |A({\bf H})|^{2},] with [a ({\bf r}) * a^{+} (-{\bf r}) = \textstyle\int a ({\bf r'}) a ({\bf r' - r}) \;\hbox{d}{\bf r'} = P({\bf r}). \eqno(] Note that [P({\bf r}) = P(-{\bf r})] is not valid if a(r) is complex. Consequently [|A(-{\bf H})|^{2} \neq |A({\bf H})|^{2}]. This is shown by evaluating [A(-{\bf H})A^{+}(-{\bf H})] [A(-{\bf H})A^{+} (-{\bf H}) \leftrightarrow a(-{\bf r}) * a^{+} ({\bf r}) = P(-{\bf r}). \eqno(] Equation ([link] is very useful for the determination of the contribution of anomalous scattering to diffuse reflections.

Most of the diffuse diffraction phenomena observed may be interpreted qualitatively or even semi-quantitatively in a very simple manner using a limited number of important Fourier transforms, given below. Fourier transforms

| top | pdf |

  • (1) Normalized Gaussian function

    [(\pi^{3/2} \alpha \beta \gamma)^{-1} \exp \{-(x / \alpha)^{2} - (y / \beta)^{2} - (z / \gamma)^{2}\}. \eqno(] This plays an important role in statistics. Its Fourier transform is again a Gaussian: [\exp \{-\pi^{2} (\alpha^{2} H^{2} + \beta^{2} K^{2} + \gamma^{2} L^{2})\}. \eqno(] The three parameters α, β, γ determine the width of the curve. Small values of α, β, γ represent a broad maximum in reciprocal space but a narrow one in real space, and vice versa. The constant has been chosen such that the integral of the Gaussian is unity in real space. The product of two Gaussians in reciprocal space [\eqalignno{&\exp \{- \pi^{2} (\alpha_{1}^{2} H^{2} + \beta_{1}^{2} K^{2} + \gamma_{1}^{2} L^{2})\} &\cr &\quad \times \exp \{- \pi^{2} (\alpha_{2}^{2} H^{2} + \beta_{2}^{2} K^{2} + \gamma_{2}^{2} L^{2})\} &\cr &\quad \quad = \exp \{- \pi^{2} [(\alpha_{1}^{2} + \alpha_{2}^{2})H^{2} + (\beta_{1}^{2} + \beta_{2}^{2}) K^{2} &\cr &\qquad \quad + (\gamma_{1}^{2} + \gamma_{2}^{2}) L^{2}]\} &(}] again represents a Gaussian of the same type, but with a sharper profile. Consequently its Fourier transform given by the convolution of the transforms of the two Gaussians is itself a Gaussian with a broader maximum. It may be concluded from this discussion that the Gaussian with [\alpha, \beta, \gamma \rightarrow 0] is a δ function in real space, and its Fourier transform is unity in reciprocal space. The convolution of two δ functions is again a δ function.

  • (2) Lattices

    Lattices in real and reciprocal space may be described by δ functions [l({\bf r}) = \textstyle\sum\limits_{{\bf n}} \delta ({\bf r - n})] and [L({\bf H}) = \textstyle\sum\limits_{{\bf h}} \delta ({\bf H - h}),] where n, h represent the components of the displacement vectors in real and reciprocal space, respectively. The Fourier transforms of lattices with orthogonal basis vectors of unit length and an infinite number of points in all three dimensions correspond to each other. In the following the relation [l({\bf r}) \leftrightarrow L({\bf H})] is used in this generalized sense.

    The Fourier transforms of finite lattices are given by [{\sin \pi N_{1} H \over \sin \pi H} {\sin \pi N_{2} K \over \sin \pi K} {\sin \pi N_{3} L \over \sin \pi L}, \eqno(] which is a periodic function in reciprocal space, but, strictly speaking, non-periodic in real space. It should be pointed out that the correspondence of lattices in either space is valid only if the origin coincides with a δ function. This fact may easily be understood by applying the law of displacement given in equation ([link].

  • (3) Box functions

    The Fourier transform of a box function b(r) with unit height is: [b({\bf r}) \leftrightarrow {\sin \pi \alpha H \over \pi H} {\sin \pi \beta K \over \pi K} {\sin \pi \gamma L \over \pi L}. \eqno(] α, β, γ describe its extension in the three dimensions. This function is real as long as the centre of symmetry is placed at the origin, otherwise the law of displacement has to be used. The convolution of the box function is needed for the calculation of intensities: [\eqalignno{t({\bf r}) &= b({\bf r}) * b(-{\bf r}) &\cr &\leftrightarrow \left({\sin \pi \alpha H \over \pi H}\right)^{2} \left({\sin \pi \beta K \over \pi K}\right)^{2} \left({\sin \pi \gamma L \over \pi L}\right)^{2}. &(}] t(r) is a generalized three-dimensional `pyramid' of doubled basal length when compared with the corresponding length of the box function. The top of the pyramid has a height given by the number of unit cells covered by the box function. Obviously, the box function generates a particle size in real space by multiplying the infinite lattice l(r) by b(r). Fourier transformation yields a particle-size effect well known in diffraction. Correspondingly, the termination effect of a Fourier synthesis is caused by multiplication by a box function in reciprocal space, which causes a broadening of maxima in real space.

  • (4) Convolutions

    It is often very useful to elucidate the convolution given in equations ([link] by introducing the corresponding pictures in real or reciprocal space. Since [1\cdot f({\bf r}) = f({\bf r}), \delta ({\bf H}) * F({\bf H}) = F({\bf H})] the convolution with a δ function must result in an identical picture of the second function, although the function is used as [f(-{\bf r})] in the integrals of equations ([link], [f({\bf r} - {\bf r}')] with r′ as variable in the integral of convolution. The convolution with [f(-{\bf r})] brings the integral into the form [\textstyle\int f({\bf r}')\;f({\bf r}' - {\bf r}) \;\hbox{d}{\bf r}', \eqno(] which is known as the Patterson function (or self- or auto-convolution) and represents the generalized Patterson function including anomalous scattering [cf. equation (][link].

    The change of the variable in the convolution integral may sometimes lead to confusion if certain operations are applied to the arguments of the functions entering the integral. Hence, it seems to be useful to mention the invariance of the convolution integral with respect to a change of sign, or a displacement, respectively, if applied to r′ in both functions. Consequently, the convolution with the inverted function [a({\bf r}) * b(-{\bf r})] may be determined as follows: [\eqalignno{b'({\bf r}) &= b(-{\bf r}) &\cr a({\bf r}) * b(-{\bf r}) &= a({\bf r}) * b'({\bf r}) = \textstyle\int a({\bf r}') b' ({\bf r} -{\bf r}') \;\hbox{d}{\bf r}' &\cr &= \textstyle\int a({\bf r}') b({\bf r}' -{\bf r}) \;\hbox{d}{\bf r}' = P'({\bf r}). &(}] This equation means that the second function is displaced into the positive direction by r, then multiplied by the first function and integrated. In the original meaning of the convolution the operation represents a displacement of the second function into the positive direction and an inversion at the displaced origin before multiplication and subsequent integration. On comparing both operations it may be concluded that [P'({\bf r}) \neq P'(-{\bf r})] if the second function is acentric. For real functions both have to be acentric. In a similar way it may be shown that the convolution of [\eqalignno{a({\bf r - m}) * b({\bf r} - {\bf m}') &= \textstyle\int\limits_{{\bf r}'} a({\bf r}' - {\bf m}) b({\bf r} - {\bf m}' - {\bf r}') \;\hbox{d}{\bf r}' &\cr &= \textstyle\int\limits_{{\bf r}''} a({\bf r}'') b({\bf r} - {\bf m}' - {\bf m} - {\bf r}'') \;\hbox{d}{\bf r}''. &\cr&&(}] Equation ([link] indicates a displacement by [{\bf m}' + {\bf m}] with respect to the convolution of the undisplaced functions. Consequently [\delta ({\bf r - m}) * \delta ({\bf r} - {\bf m}') = \delta ({\bf r} - {\bf m} - {\bf m}'). \eqno(] Obviously, the commutative law of convolution is obeyed; on the other hand, the convolution with the inverted function yields [\delta ({\bf r - m' + m}),] indicating that the commutative law (interchange of m and m′) is violated because of the different signs of m and m′.

    The effectiveness of the method outlined above may be greatly improved by introducing further Fourier transforms of useful functions in real and reciprocal space (Patterson, 1959[link]). Applications

| top | pdf |

  • (1) Clusters in a periodic lattice (low concentrations)

    The exsolution of clusters of equal sizes is considered. The lattice of the host is undistorted and the clusters have the same lattice but a different structure. A schematic drawing is shown in Fig.[link] Two different structures are introduced by [\eqalign{F_{1} ({\bf r}) &= \textstyle\sum\limits_{\nu} \delta ({\bf r - r}_{\nu}) * F_{\nu} ({\bf r})\cr F_{2} ({\bf r}) &= \textstyle\sum\limits_{\mu} \delta ({\bf r - r}_{\mu}) * F_{\mu} ({\bf r}).}] Their Fourier transforms are the structure factors [F_{1}({\bf H})], [F_{2}({\bf H})]. The bold lines in Fig.[link] indicate the clusters, which may be represented by box functions [b({\bf r})] in the simplest case. It should be pointed out, however, that a more complicated shape means nothing other than a replacement of [b({\bf r})] by another shape function [b'({\bf r})] and its Fourier transform [B'({\bf H})]. The distribution of clusters is represented by [d({\bf r}) = \textstyle\sum\limits_{{\bf m}} \delta ({\bf r - m}),] where m refers to the centres (crosses in Fig.[link] of the box functions. The problem is therefore defined by: [l({\bf r}) * F_{1} ({\bf r}) + [l({\bf r}) b({\bf r})] * [F_{2}({\bf r}) - F_{1} ({\bf r})] * d({\bf r}). \eqno(] The incorrect addition of [F_{1}({\bf r})] to the areas of clusters [F_{2}({\bf r})] is compensated by subtracting the same contribution from the second term in equation ([link]. In order to determine the diffuse scattering the Fourier transformation of ([link] is performed: [L({\bf H}) F_{1} ({\bf H}) + [L({\bf H}) * B({\bf H})][F_{2} ({\bf H}) - F_{1} ({\bf H})] D({\bf H}). \eqno(] The intensity is given by [{|L({\bf H}) F_{1} ({\bf H}) + [L({\bf H}) * B({\bf H})] [F_{2} ({\bf H}) - F_{1} ({\bf H})] D({\bf H})|^{2}.} \eqno(] Evaluation of equation ([link] yields three terms (c.c. means complex conjugate): [\displaylines{(\hbox{i}) \qquad\ |L ({\bf H}) F_{1} ({\bf H})|^{2} \hfill\cr (\hbox{ii}) \qquad \{[L ({\bf H}) F_{1} ({\bf H})][L ({\bf H}) * B ({\bf H})]\hfill\cr \qquad\ \phantom{[L ({\bf H}) F_{1}]} \times [F_{2} ({\bf H}) - F_{1} ({\bf H})] D ({\bf H}) + \hbox{c.c.}\}\hfill\cr (\hbox{iii}) \qquad |[L ({\bf H}) * B ({\bf H})][F_{2} ({\bf H}) - F_{1} ({\bf H})] D ({\bf H})|^{2}.\hfill}] The first two terms represent modulated lattices [multiplication of [L({\bf H})] by [F_{1}({\bf H})]]. Consequently, they cannot contribute to diffuse scattering which is completely determined by the third term. Fourier transformation of this term gives [[l({\bf r}) = l(-{\bf r})]; [b({\bf r}) = b(-{\bf r})]; [\Delta F = F_{2} - F_{1}]]: [\eqalignno{&[l({\bf r}) b({\bf r})] * \Delta F({\bf r}) * d({\bf r}) * [l({\bf r}) b({\bf r})] * \Delta F^{+} (-{\bf r}) * d(-{\bf r}) &\cr &\quad = [l({\bf r}) b({\bf r})] * [l({\bf r}) b({\bf r})] * \Delta F({\bf r}) * \Delta F^{+} (-{\bf r}) * d({\bf r}) * d(-{\bf r}) &\cr &\quad = [l({\bf r}) t({\bf r})] * \Delta F({\bf r}) * \Delta F^{+} (-{\bf r}) * d({\bf r}) * d(-{\bf r}). &(}] According to equation ([link] and its subsequent discussion the convolution of the two expressions in square brackets was replaced by l(r)t(r), where t(r) represents the `pyramid' of n-fold height discussed above and n is the number of unit cells within b(r). [d({\bf r}) * d(-{\bf r})] is the Patterson function of the distribution function d(r). Its usefulness may be recognized by considering the two possible extreme solutions, namely the random and the strictly periodic distribution.


    Figure | top | pdf |

    Model of the two-dimensional distribution of point defects, causing changes in the surroundings.

    If no fluctuations of domain sizes are admitted the minimum distance between two neighbouring domains is equal to the length of the domain in the corresponding direction. This means that the distribution function cannot be completely random. In one dimension the solution of a random distribution of particles of a given size on a finite length shows that the distribution functions exhibit periodicities depending on the average free volume of one particle (Zernike & Prins, 1927[link]). Although the problem is more complicated in three dimensions, there should be no fundamental difference in the exact solutions.

    On the other hand, it may be shown that the convolution of a pseudo-random distribution may be obtained if the average free volume is large. This is shown in Fig.[link](a) for the particular case of a cluster smaller than one unit cell. A strictly periodic distribution function (superstructure) may result, however, if the volume of the domain and the average free volume are equal. Obviously, the practical solution for the self-convolution of the distribution function (= Patterson function) lies somewhere in between, as shown in Fig.[link](b). If a harmonic periodicity damped by a Gaussian is assumed, this self-convolution of the distribution in real space may be considered to consist of two parts, as shown in Figs.[link](c), (d). Note that the two different solutions result in completely different diffraction patterns:


    Figure | top | pdf |

    One-dimensional Patterson functions of various point-defect distributions: (a) random distribution; (b) influence of finite volume of defects on the distribution function; (c), (d) decomposition of (b) into a periodic (c) and a convergent (d) part; (e) Fourier transform of (c) + (d); (f) changes of (e), if the centres of the defects show major deviations from the origins of the lattice.

    • (i) The geometrically perfect lattice extends to distances which are large when compared with the correlation length of the distribution function. Then the Patterson function of the distribution function concentrates at the positions of the basic lattice, which is given by multiplication by the lattice l(r). The corresponding convolution in reciprocal space gives the same contribution to all Bragg reflections (Fig.[link]e).

    • (ii) There is no perfect lattice geometry. In this case a continuous Patterson function results. Fourier transformation yields an influence which is now restricted primarily to the reflection 000, i.e. to the low-angle diffraction range.

    Figs.[link](e), (f) show the different diffraction patterns of the diffuse scattering which is concentrated around the Bragg maxima. Although the discussion of the diffuse scattering was restricted to the case of identical domains, the introduction of a distribution of domain sizes does not influence the diffraction pattern essentially, as long as the fluctuation of sizes is small compared with the average volume of domain sizes and no strong correlation exists between domains of any size (size-independent random distribution).

    The complete qualitative discussion of the diffraction pattern may be made by investigating the Fourier transform of ([link]: [[L({\bf H}) * T({\bf H})]|\Delta F ({\bf H})|^{2}|D ({\bf H})|^{2}. \eqno(] The first factor in ([link] describes the particle-size effect of a domain containing the influence of a surrounding strain field and the new structure of the domains precipitated from the bulk. D(H) has its characteristic variation near the Bragg peaks (Figs.,f), and is less important in between. For structure determination of domains, intensities near the Bragg peaks should be avoided. Note that equation ( may be used for measurements applying anomalous scattering in both the centric and the acentric case.

    Solution of the diffraction problem. In equation ([link] [\Delta F ({\bf H})] is replaced by its average [\langle \Delta F ({\bf H}) \rangle = \textstyle\sum\limits_{\mu} p_{\mu} \Delta F_{\mu} ({\bf H}),] where [p_{\mu}] represents the a priori probability of a domain of type µ. This replacement becomes increasingly important if small clusters (domains) have to be considered. Applications of the formulae to Guinier–Preston zones are given by Guinier (1942)[link] and Gerold (1954)[link]; a similar application to clusters of vacancies in spinels with an excess in Al2O3 was outlined by Jagodzinski & Haefner (1967)[link].

    Although refinement procedures are possible in principle, the number of parameters entering the diffraction problem becomes increasingly large if more clusters or domains (different sizes) have to be introduced. Another difficulty results from the large number of diffraction data which must be collected to perform a reliable structure determination. There is no need to calculate the first two terms in equation ([link] which do contribute to the sharp Bragg peaks only, because their intensity is simply described by the averaged structure factor [|\langle F ({\bf H})\rangle |^{2}]. These terms may therefore be replaced by [|L({\bf H})|^{2}| \langle F({\bf H})\rangle |^{2}] with [|\langle F({\bf H})\rangle|^{2} = \left|\textstyle\sum\limits_{\mu} p_{\mu} F_{\mu} ({\bf H})\right|^{2} \eqno(] where [p_{\mu}] is the a priori probability of the structure factor [F_{\mu}({\bf H})]. It should be emphasized here that ([link] is independent of the distribution function d(r), or its Fourier transform D(H).

  • (2) Periodic lamellar domains

    Here d(r) is one-dimensional, and can easily be calculated: a periodic array of two types of lamellae having the same basic lattice l(r), but a different structure, is shown in Fig.[link] The size of the two types of lamellae may be different. The structure of the first domain type is given by a convolution with [F_{1} ({\bf r})] (Fig.[link]a) and that of the second domain type by [F_{2} ({\bf r})] (Fig.[link]b). Introducing [\langle F ({\bf r})\rangle] and [\Delta F ({\bf r})], the structure in real space is described by: [\eqalignno{&[l ({\bf r}) b_{1} ({\bf r})] * d ({\bf r}) * F_{1} ({\bf r}) + [l ({\bf r}) b_{2} ({\bf r})] * d ({\bf r}) * F_{2} ({\bf r}) &\cr &\quad = \{[l ({\bf r}) b_{1} ({\bf r}) + l ({\bf r}) b_{2} ({\bf r})] * d ({\bf r})\} * \langle F ({\bf r})\rangle &\cr &\qquad + [l ({\bf r}) b_{1} ({\bf r}) - l ({\bf r}) b_{2} ({\bf r})] * d ({\bf r}) * \Delta F ({\bf r}). &(}] Obviously the first term in curly brackets in equation ([link] is no more than l(r) itself and d(r) is strictly periodic. [b_{1} ({\bf r})] and [b_{2} ({\bf r})] are box functions, mutually displaced by [\pm (n_{1} + n_{2}) / 2] unit cells in the c direction [[n_{1}, n_{2}] are the numbers of cells covered by [b_{1} ({\bf r})] and [b_{2} ({\bf r})], respectively].


    Figure | top | pdf |

    Periodic array of domains consisting of two different atoms, represented by different heights. (a) Distribution of domain type 1, (b) distribution of domain type 2.

    Fourier transformation of equation ([link] yields [{L ({\bf H}) \langle F ({\bf H})\rangle + \{L ({\bf H}) * [B_{1} ({\bf H}) - B_{2} ({\bf H})]\} D ({\bf H}) \Delta F ({\bf H}).} \eqno(] The first term in equation ([link] gives the normal sharp reflections of the average structure, while the second describes superlattice reflections [sublattice [L_{s} ({\bf H}) = D ({\bf H})] in reciprocal space], multiplied by [\Delta F ({\bf H})] and another `structure factor' generated by the convolution of the reciprocal lattice L(h) with [[B_{1} ({\bf H}) - B_{2} ({\bf H})]] (cf. Fig.[link] b). Since the centres of [b_{1} ({\bf r})] and [b_{2} ({\bf r})] are mutually displaced, the expression in square brackets includes extinctions if [b_{1} ({\bf r})] and [b_{2} ({\bf r})] represent boxes equal in size. These extinctions are discussed below. It should be pointed out that [L_{s} ({\bf H})] and its Fourier transform [l_{s} ({\bf r})] are commensurate with the basic lattice, as long as no change of the translation vector at the interface of the lamellae occurs. Obviously, [L_{s} ({\bf H})] becomes incommensurate in the general case of a slightly distorted interface. Considerations of this kind play an important role in the discussion of modulated structures.

    No assumption has been made so far for the position of the interface. This point is meaningless only in the case of a strictly periodic array of domains (no diffuse scattering). Therefore it seems to be convenient to introduce two basic vectors parallel to the interface in real space which demand a new reciprocal vector perpendicular to them defined by [({\bf a}' \times {\bf b}') / V'], where a′, b′ are the new basic vectors and V′ is the volume of the supercell. As long as the new basic vectors are commensurate with the original lattice, the direction of the new reciprocal vector [{\bf c}'^{*}], perpendicular to a′, b′, passes through the Bragg points of the original reciprocal lattice and the reciprocal lattice of the superlattice remains commensurate as long as V′ is a multiple of V [(V' = mV,\ m = \hbox{integer})]. Since the direction of c is arbitrary to some extent, there is no clear rule about the assignment of superlattice reflections to the original Bragg peaks. This problem becomes very important if extinction rules of the basic lattice and the superlattice have to be described together.


    We consider a b.c.c. structure with two kinds of atoms (1, 2) with a strong tendency towards superstructure formation (CsCl-type ordering). According to equations ([link],c[link]) and ([link]) we obtain, in the case of negligible short-range order, the following expressions for sharp and diffuse scattering (c = concentration): [\eqalign{I_{s} &= |c F_{1} ({\bf H}) + (1 - c) F_{2} ({\bf H})|^{2}\quad \hbox{for } h + k + l = 2n \cr I_{d} &= c (1 - c)| F_{1} ({\bf H}) - F_{2} ({\bf H})|^{2}\quad \hbox{elsewhere}.}]

    With increasing short-range order the sharp reflections remain essentially unaffected, while the diffuse ones concentrate into diffuse maxima at h with [h + k + l = 2n + 1]. This process is treated more extensively below. As long as the domains exhibit no clear interface, it is useful to describe the ordering process with the two possible cell occupations of a pair of different atoms; then contributions of equal pairs may be neglected with increasing short-range order. Now the two configurations 1, 2 and 2, 1 may be given with the aid of the translation [{1\over 2} ({\bf a} + {\bf b} + {\bf c})]. Hence the two structure factors are [F_{1} \hbox{ and } F_{2} = F_{1} \exp \{\pi i (h + k + l)\}.] Since both structure factors occur with the same probability, the equations for sharp and diffuse reflections become [\eqalign{I_{s} &= {\textstyle{1\over 4}} |F_{1} ({\bf H})|^{2} [1 + \exp \{\pi i (h + k + l)\}]^{2} \cr I_{d} &= {\textstyle{1\over 4}} |F_{1} ({\bf H})|^{2} [1 - \exp \{\pi i (h + k + l)\}]^{2}.}] It is well recognized that no sharp reflections may occur for [h + k + l = 2n + 1], and the same holds for the diffuse scattering if [h + k + l = 2n]. This extinction rule for diffuse scattering is due to suppression of contributions of equal pairs. The situation becomes different for lamellar structures. Let us first consider the case of lamellae parallel to (100). The ordered structure is formed by an alternating sequence of monoatomic layers, consisting of atoms of types 1 and 2, respectively. Hence, the interface between two neighbouring domains is a pair of equal layers 1,1 or 2,2 which are not equivalent. Each interface of type 1 (2) may be described by an inserted layer of type 1 (2), and the chemical composition differs from 1:1 if one type of interface is preferred. Since the contribution of equal pairs has been neglected in deriving the extinction rule of diffuse scattering (see above) this rule is no longer valid. Because of the lamellar structure the diffuse intensity is concentrated into streaks parallel to (00l). Starting from the diffuse maximum (010), the diffuse streak passes over the sharp reflection 011 to the next diffuse one 012 etc., and the extinction rule is violated as long as one of the two interfaces is predominant. Hence, the position of the interface determines the extinction rule in this orientation.

    A completely different behaviour is observed for lamellae parallel to (110). This structure is described by a sequence of equal layers containing 1 and 2 atoms. The interface between two domains (exchange of the two different atoms) is now nothing other than the displacement parallel to the layer of the original one in the ordered sequence. Calculation of the two structure factors would involve displacements [\pm {1\over 2} ({\bf a} - {\bf b})]. Starting from the diffuse reflection 001, the diffuse streak parallel to (HH0) passes through (111), (221), (331), …; i.e. through diffuse reflections only. On the other hand, rows (HH2) going through (002), (112), (222), … do not show any diffuse scattering. Hence, we have a new extinction rule for diffuse scattering originating from the orientation of interfaces. This fact is rather important in structure determination. For various reasons, lamellar interfaces show a strong tendency towards a periodic arrangement. In diffraction the diffuse streak then concentrates into more or less sharp superstructure reflections. These are not observed on those rows of the reciprocal lattice which are free from diffuse scattering. The same extinction law is not valid in the case of the (100) orientation of the interfaces. Summarizing, we may state that three types of extinction rules have to be considered:

    • (a) Normal extinctions for the average structure.

    • (b) Extinction of the difference structure factors for diffuse scattering.

    • (c) Extinctions caused by the ordering process itself.

  • (3) Lamellar system with two different structures, where [\langle F ({\bf H})\rangle] and [\Delta F ({\bf H})] do not obey any systematic extinction law

    The convolution of the second term in equation ([link] (cf. Fig.[link] may be represented by a convolution of the Fourier transform of a box function [B_{1} ({\bf H})] with the reciprocal superlattice. Since [B_{1} ({\bf H})] is given by [\sin(\pi m_{s} H) / (\pi H)], where [m_{s} = \hbox{number of cells}] of the supercell, the reader might believe that the result of the convolution may easily be determined quantitatively: this assumption is not correct because of the slow convergence of [B_{1} ({\bf H})]. The systematic concidences of the maxima, or minima, of [B_{1} ({\bf H})] with the points of the superlattice in the commensurate case cause considerable changes in intensities especially in the case of a small thickness of the domains. For this reason an accurate calculation of the amplitudes of satellites is necessary (Jagodzinski & Penzkofer, 1981[link]):

    • (a) Bragg peaks of the basic lattice [I \sim |\langle F ({\bf H})\rangle|^{2}\semi \eqno(]

    • (b) satellites: [\nu = 2n] (n = integer except 0) [I \sim |2 \sin \pi \nu C / [\sin \pi \nu / (N_{1} + N_{2})] \Delta F ({\bf H})|^{2}\hbox{;} \eqno(]

    • (c) satellites: [\nu = 2n + 1] [I \sim |2 \cos \pi \nu C / [\sin \pi \nu / (N_{1} + N_{2})] \Delta F ({\bf H})|^{2}. \eqno(] [\nu] = order of satellites, [C = {1\over 2} (N_{1} - N_{2}) / (N_{1} + N_{2})], [N_{1}, N_{2}=] number of cells within [b_{1} ({\bf r})] and [b_{2} ({\bf r})], respectively.

    Obviously, there is again a systematic extinction rule for even satellites if [N_{1} = N_{2}].

    Equation ([link] indicates an increasing intensity of first even-order satellites with increasing C. Intensities of first even and odd orders become nearly equal if [N_{2} \simeq {1\over 2} N_{1}]. Smaller values of [N_{2}] result in a decrease of intensities of both even and odd orders (no satellites if [N_{2} = 0]). The denominators in equations ([link],c[link]) indicate a decrease in intensity with increasing order of satellites. The quantitative behaviour of the intensities needs a more detailed discussion of the numerator in equations ([link],b[link]) with increasing order of satellites. Obviously, there are two kinds of extinction rules to be taken into account: systematic absences for the various orders of satellites, and the usual extinctions for [\langle F ({\bf H})\rangle] and [\Delta F ({\bf H})]. Both have to be considered separately in order to arrive at reliable conclusions. This different behaviour of the superlattice reflections (satellites) and that of the basic lattice may well be represented by a multi-dimensional group-theoretical representation as has been shown by de Wolff (1974[link]), Janner & Janssen (1980a[link],b[link]), de Wolff et al. (1981[link]), and others.

  • (4) Non-periodic system (qualitative discussion)

    Following the discussion of equations ([link] [link] [link] one may conclude that the fluctuations of domain sizes cause a broadening of satellites, if the periodic distribution function has to be replaced by a statistical one. In this case the broadening effect increases with the order of satellites. The intensities, however, are completely determined by the distribution function and can be estimated by calculating the intensities of the perfectly ordered array, as approximated by the distribution function.

    A careful check of [\langle F ({\bf H})\rangle] and [\Delta F ({\bf H})] in equations ([link] [link] [link] shows that the position of the interface plays an important role for the intensities of satellites. Since this position determines the origin of the unit cells in the sublattice, we have to choose this origin for the calculation of F(H) and [\Delta F({\bf H})]. This involves phase factors which are meaningless for integral values of H, (i) if the average [\langle F ({\bf H})\rangle] refers to different structures with arbitrary origin, or (ii), which is important for practical cases, where no change occurs in the origin of related structures for neighbouring domains which are bound to an origin by general convention (e.g. a centre of symmetry). This statement is no longer true for non-integral values of H which are needed for the calculation of intensities of satellites. The intensities of satellites become different for different positions of the interface, even in the absence of a relative displacement between neighbouring domains with respect to an origin by convention. This statement may be extended to non-periodic distribution functions. Consequently, one may conclude that the study of diffuse scattering yields information on the interfacial scattering. For slightly different structures at the interface two cases are important:

    • (i) the two structures are related by symmetry (e.g. by a twin law); and

    • (ii) the difference between the two structures cannot be described by a symmetry operation.

    In structures based on the same sublattice, the first case seems to be more important, because two different structures with the same sublattice are improbable. In the first case there is an identical sublattice if the symmetry operation in question does not influence the plane of intergrowth, e.g. a mirror plane should coincide with the plane of intergrowth. Since we have two inequivalent mirror planes in any sublattice, there are two such planes. It is assumed that no more than one unit cell of both domains at the interface has a slightly different structure without any change of geometry of the unit cell, and the number of unit cells is equal because of the equivalence of both domain structures (twins). Fig. [link] shows a picture of this model; Figs., (c) [link] explain that this structure may be described by two contributions:


    Figure | top | pdf |

    Influence of distortions at the boundary of domains, and separation into two parts; for discussion see text.

    • (i) The first term is already given by equation ([link] [link] [link] for [N_{1} = N_{2}], consequently odd orders of satellites only are observed.

    • (ii) The second term may be described by a superlattice containing [2N_{1}] cells with an alternating arrangement of interfaces, correlated by the relevant plane of symmetry.

    In real space the second term may be constructed by convolution of the one-dimensional superlattice with two difference structures displaced by [\mp N_{1} / 2] units of the sublattice; its Fourier transformation yields [\eqalignno{&L_{s} ({\bf H}) \{\Delta F_{i} ({\bf H}) \exp \{2 \pi i N_{1} H / 2\} &\cr &\quad \qquad + \Delta F'_{i} ({\bf H}) \exp \{-2 \pi i N_{1} H / 2\}\}, &(}] where [\Delta F_{i}, \Delta F'_{i}] correspond to the Fourier transforms of the contributions shown in Fig.[link](c). Since [H = \nu / N_{1}] there are alternating contributions to the νth satellite, which may be calculated more accurately by taking into account the symmetry operations. The important difference between equations ([link] [link] [link] and ([link] is the missing decrease in intensity with increasing order of satellites. Consequently one may conclude that the interface contributes to low- and high-order satellites as well, but its influence prevails for high-order satellites. Similar considerations may be made for two- and three-dimensional distributions of domains. A great variety of extinction rules may be found depending on the type of order approximated by the distribution under investigation.

  • (5) Two kinds of lamellar domains with variable size distribution

    Obviously the preceding discussion of the diffuse scattering from domains is restricted to more or less small fluctuations of domain sizes. This is specifically valid if the most probable domain size does not differ markedly from the average size. The condition is violated in the case of order–disorder phenomena. It may happen that the smallest ordered area is the most probable one, although the average is considerably larger. This may be shown for a lamellar structure of two types of layers correlated by a (conditional) pair probability [p_{\mu \mu'}({\bf 1})]. As shown below, a pair at distance m occurs with the probability [p_{\mu} p_{\mu \mu'} ({\bf m})] which may be derived from the pair-probability of nearest neighbours [p_{\mu} p_{\mu \mu'}({\bf 1})]. (In fact only one component of vector m is relevant in this context.) The problem will be restricted to two kinds of layers [(\mu, \mu' = 1, 2)]. Furthermore, it will be symmetric in the sense that the pair probabilities obey the following rules [p_{11} ({\bf m}) = p_{22} ({\bf m}), \quad p_{12} ({\bf m}) = p_{21} ({\bf m}). \eqno(]

    It may be derived from equation ([link] that the a priori probabilities [p_{\mu}] of a single layer are [{1\over 2}] and [p_{11} ({\bf 0}) = p_{22} ({\bf 0}) = 1, \quad p_{12} ({\bf 0}) = p_{21} ({\bf 0}) = 0.] With these definitions and the general relation [p_{11} ({\bf m}) + p_{12} ({\bf m}) = p_{22} ({\bf m}) + p_{21} ({\bf m}) = 1] the a priori probability of a domain containing m layers of type 1 may be calculated with the aid of [p_{11} ({\bf 1})] [[0 \leq p_{11} ({\bf 1}) \leq 1]]: [p_{\mu} = {\textstyle{1\over 2}}\; p_{11} ({\bf 1})^{m - 1} [1 - p_{11} ({\bf 1})]. \eqno(] Hence the most probable size of domains is a single layer because a similar relation holds for layers of type 2. Since the average thickness of domains is strongly dependent on [p_{11}({\bf 1})] [infinite for [p_{11} ({\bf 1}) = 1], and one layer for [p_{11} ({\bf 1}) = 0]] it may become very large in the latter case. Consequently there are extremely large fluctuations if [p_{11} ({\bf 1})] is small, but different from zero.

    It may be concluded from equation ([link] that the function [p_{11} ({\bf m})] decreases monotonically with increasing m, approaching [{1\over 2}] with [{\bf m} \rightarrow \infty]. Apparently this cannot be true for a finite crystal if [p_{11} ({\bf m})] is unity (structure of two types of domains) or zero (superstructure of alternating layers). In either case the crystal should consist of a single domain of type 1 or 2, or one of the possible superstructures 1212 …, 2121 …, respectively. Hence one has to differentiate between long-range order, where two equivalent solutions have to be considered, and short-range order, where [p_{11} ({\bf m})] approaches the a priori probability [{1\over 2}] for large m. This behaviour of [p_{11} ({\bf m})] and [p_{12} ({\bf m})], which may also be expressed by equivalent correlation functions, is shown in Figs.[link](a) (short-range order) and[link](b) (long-range order). [p_{11} ({\bf m})] approaches [{1\over 2} + s] for large m with [s = 0] in the case of short-range order, while [p_{12} ({\bf m})] becomes [{1\over 2} - s]. Obviously a strict correlation between [p_{11} ({\bf 1})] and s exists which has to be calculated. For a qualitative interpretation of the diffraction pictures this correlation may be derived from the diffraction pattern itself. The [p_{\mu \mu'} ({\bf m})] are separable into a strictly periodic and a monotonically decreasing term approaching zero in both cases. This behaviour is shown in Figs., (b)[link]. The periodic term contributes to sharp Bragg scattering. In the case of short-range order the symmetry relations given in equation ([link] are valid. The convolution in real space yields with factors [t({\bf r})] (equations[link] [link] [link]): [\eqalign{&\textstyle{1\over 2} t ({\bf r}) [\textstyle\sum\limits_{{\bf m}} \delta ({\bf r} + {\bf m}) p'_{11} ({\bf m})] * F_{1} ({\bf r}) * F_{1} (-{\bf r})\cr &\quad + \textstyle{1\over 2} t ({\bf r}) [\textstyle\sum\limits_{{\bf m}} \delta ({\bf r} + {\bf m}) p'_{12} ({\bf m})] * F_{1} ({\bf r}) * F_{2} (-{\bf r})\cr &\quad + \textstyle{1\over 2} t ({\bf r}) [\textstyle\sum\limits_{{\bf m}} \delta ({\bf r} + {\bf m}) p'_{21} ({\bf m})] * F_{2} ({\bf r}) * F_{1} (-{\bf r})\cr &\quad + \textstyle{1\over 2} t ({\bf r}) [\textstyle\sum\limits_{{\bf m}} \delta ({\bf r} + {\bf m}) p'_{22} ({\bf m})] * F_{2} ({\bf r}) * F_{2} (-{\bf r}),}] where [p'_{\mu \mu'} ({\bf m})] are factors attached to the δ functions: [\eqalign{p'_{11} ({\bf m}) &= p_{11} ({\bf m}) - {\textstyle{1\over 2}} = p'_{22} ({\bf m})\cr p'_{12} ({\bf m}) &= p'_{21} ({\bf m}) = -p'_{11} ({\bf m}).}] The positive sign of n in the δ functions results from the convolution with the inverted lattice [cf. Patterson (1959[link], equation 32)]. Fourier transformation of the four terms given above yields the four corresponding expressions [(\mu , \mu' !)]: [{\textstyle{1\over 2} [T({\bf H}) * \textstyle\sum\limits_{{\bf m}} p'_{\mu \mu'} ({\bf m}) \exp \{-2 \pi i {\bf H} \cdot {\bf m}\}] F_{\mu} ({\bf H}) F_{\mu'}^{+} ({\bf H}).} \eqno(] Now the summation over m may be replaced by an integral if the factor [l({\bf m})] is added to [p'_{\mu \mu'} ({\bf m})], which may then be considered as the smoothest continuous curve passing through the relevant integer values of m: [\textstyle\sum \rightarrow \textstyle\int l({\bf m}) p'_{\mu \mu'} ({\bf m}) \exp \{-2 \pi i {\bf H} \cdot {\bf m}\} \;\hbox{d}{\bf m}] since both [l({\bf m})] and [p'_{\mu \mu'} ({\bf m})] are symmetric in our special case we obtain [\textstyle\sum = L ({\bf H}) * P'_{\mu \mu'} ({\bf H}).] Insertion of the sum in equation ([link] results in [{\textstyle{1\over 2}} [L ({\bf H}) * T({\bf H}) * P'_{\mu \mu'} ({\bf H})] F_{\mu} ({\bf H}) F_{\mu'}^{+} ({\bf H}). \eqno(] Using all symmetry relations for [p'_{\mu \mu'} ({\bf m})] and [P'_{\mu \mu'} ({\bf H})], respectively, we obtain for the diffuse scattering after summing over μ, μ′ [I_{d} \approx [L ({\bf H}) * T({\bf H}) * P'_{11} ({\bf H})]|\Delta F ({\bf H})|^{2} \eqno(] with [\Delta F ({\bf H}) = {\textstyle{1\over 2}} [F_{1} ({\bf H}) - F_{2} ({\bf H})]].


    Figure | top | pdf |

    Typical distributions of mixed crystals (unmixing): (a) upper curve: short-range order only; (b) lower curve: long-range order.


    Figure | top | pdf |

    Decomposition of Fig.[link](a) into a periodic and a rapidly convergent part.

    It should be borne in mind that [P'_{11} ({\bf H})] decreases rapidly if [p'_{11} ({\bf r})] decreases slowly and vice versa. It is interesting to compare the different results from equations ([link] and ([link]. Equation ([link] indicates diffuse maxima at the positions of the sharp Bragg peaks, while the multiplication by [D({\bf H})] causes satellite reflections in the neighbourhood of Bragg maxima. Both equations contain the factor [|\Delta F ({\bf H})|^{2}] indicating the same influence of the two structures. More complicated formulae may be derived for several cell occupations. In principle, a result similar to equation ([link] will be obtained, but more interdependent correlation functions [p_{\mu \mu'} ({\bf r})] have to be introduced. Consequently, the behaviour of diffuse intensities becomes more differentiated in so far as all [p'_{\mu \mu'}({\bf r})] are now correlated with the corresponding [\Delta F_{\mu} ({\bf r}), \Delta F_{\mu'} (-{\bf r})]. Hence the method of correlation functions becomes increasingly ineffective with increasing number of correlation functions. Here the cluster method seems to be more convenient and is discussed below.

  • (6) Lamellar domains with long-range order: tendency to exsolution

    The Patterson function of a disordered crystal exhibiting long-range order is shown in Fig.[link](b). Now [p_{11} (\infty)] converges against [{1\over 2} + s], the a priori probability changes correspondingly. Since [p_{12} (\infty)] becomes [{1\over 2} - s], the symmetry relation given in equation ([link] is violated: [p_{11} ({\bf r}) \neq p_{22} ({\bf r})] for a finite crystal; it is evident that another crystal shows long-range order with the inverted correlation function, [p_{22} (\underline{\infty}) = {1\over 2} + s], [p_{21} (\underline{\infty}) = {1\over 2} - s], respectively, such that the symmetry [p_{11} ({\bf r}) = p_{22} ({\bf r})] is now valid for an assembly of finite crystals only. According to Fig.[link](b) there is a change in the intensities of the Bragg peaks. [\eqalign{&I_{1} \sim |({\textstyle{1\over 2}} + s) F_{1} ({\bf H}) + ({\textstyle{1\over 2}} - s) F_{2} ({\bf H})|^{2} \cr &I_{2} \sim |({\textstyle{1\over 2}} + s) F_{2} ({\bf H}) + ({\textstyle{1\over 2}} - s) F_{1} ({\bf H})|^{2},\cr} \eqno(] where [I_{1}, I_{2}] represent the two solutions discussed for the assembly of crystals which have to be added with the probability [{1\over 2}]; the intensities of sharp reflections become [I = (I_{1} + I_{2})/2. \eqno(] Introducing equation ([link] into ([link] we obtain [I \sim |{\textstyle{1\over 2}} [F_{1} ({\bf H}) + F_{2} ({\bf H})]|^{2} + s^{2} |F_{1} ({\bf H}) - F_{2} ({\bf H})^{2}|^{2}. \eqno(] [s = 0] corresponds to the well known behaviour of sharp reflections, [s = {1\over 2}] (maximum long-range order) gives [I \sim {\textstyle{1\over 2}} [|F_{1} ({\bf H})|^{2} + |F_{2} ({\bf H})|^{2}]. \eqno(] This result reveals some difficulties for structure determination of the averaged structure as long as s is different from zero or [{1\over 2}], since in the former case the use of integrated sharp Bragg intensities yields a correct average structure. If [s = {1\over 2}], a correct structure determination can only be performed with a refinement allowing for an incoherent superposition of two different structures. Having subtracted all periodic contributions to [p_{\mu \mu'} ({\bf r})], new functions which describe the remaining non-periodic parts have to be introduced (Fig.[link] b). In order to obtain a clear overview of intensities, [p'_{\mu \mu'} ({\bf r})] is again defined: [p'_{\mu \mu'} ({\bf r}) = cp_{\mu \mu'} ({\bf r}) - p_{\mu \mu'} (\infty),] where c should be chosen such that [p_{\mu \mu'} ({\bf 0}) = 1]. By this definition a very simple behaviour of the diffuse scattering is obtained: [\eqalign{&p'_{11} ({\bf r}): {\textstyle{1\over 2}} - s\hbox{;} \quad p'_{12} ({\bf r}): - ({\textstyle{1\over 2}} - s)\hbox{;}\cr &p'_{22} ({\bf r}): {\textstyle{1\over 2}} + s\hbox{;} \quad p'_{21} ({\bf r}): - ({\textstyle{1\over 2}} + s).}] With the definitions introduced above it is found that: [p'_{11} ({\bf r}) = p'_{22} ({\bf r}).] The diffuse scattering is given by: [I_d({\bf H}) = ({\textstyle{1 \over 4}}-s^2)|F_1({\bf H})-F_2({\bf H})|^2[P_{11}^{\prime}({\bf H})*L({\bf H})].\eqno(] Since equation ([link] is symmetrical with respect to an interchange of [F_{1}] and [F_{2}], the same result is obtained for [I_{2}]. Diffuse reflections occur in the positions of the sharp ones; the integrated intensities of sharp and diffuse reflections are independent of the special shape of [ P'_{11}({\bf H})]: [p_{11}({\bf 0}) = 1]; hence [1 = \textstyle\int P'_{11}({\bf H}) \exp \{2 \pi i{\bf 0\cdot H}\} \;\hbox{d}{\bf H} = \textstyle\int P'_{11}({\bf H}) \;\hbox{d}{\bf H}.]

  • (7) Lamellar domains with long-range order: tendency to superstructure

    So far it has been tacitly assumed that the crystal shows a preference for equal neighbours. If there is a reversed tendency (pairs of unequal neighbours are more probable) the whole procedure outlined above may be repeated as shown in Fig.[link] for the one-dimensional example. With the same probability of an unlike pair as used for the equal pair in the preceding example, the order process approaches an alternating structure such that the even-order neighbours have the same pair probabilities, while the odd ones are complementary for equal pairs (Fig.[link]. In order to calculate intensities, it is necessary to introduce a new lattice with the doubled lattice constant and the corresponding reciprocal lattice with [{b^{*}}'= b^{*}/2]. In order to describe the probability [p_{\mu \mu'} ({\bf r})], one has to introduce two lattices in real space – the normal lattice with the undisplaced origin and the displaced one. Fourier transformation of the new functions yields the following very similar results:


    Figure | top | pdf |

    The same distribution (cf. Fig.[link]) in the case of superstructure formation.

    • Sharp Bragg reflections

      (a) [k' = \hbox{even}] [I = |{\textstyle{1\over 2}}[F_{1} ({\bf H}) + F_{2} ({\bf H})]|^{2} \eqno(]

    • (b) [k' = \hbox{odd}] [I = s^{2}|{\textstyle{1\over 2}}[F_{1} ({\bf H}) - F_{2} ({\bf H})]|^{2} \eqno(]

    • Diffuse reflections

      (c) [k' = \hbox{odd}] [I = ({\textstyle{1\over 4}} - s^{2}) |{\textstyle{1\over 2}}[F_{1} ({\bf H}) - F_{2} ({\bf H})]|^{2}. \eqno(]

    Obviously, there is a better situation for determination of the averaged structure which may be performed without any difficulty, regardless of whether s is different from zero or not. For this purpose even reflections (or reflections in the old setting) may be used. The inclusion of odd reflections in the structure determination of the superstructure is also possible if convenient H-independent scaling factors are introduced in order to compensate for the loss in intensity which is unavoidable for the integration of the diffuse scattering.

    A few comments should be made on the physical meaning of the formulae derived above. All formulae may be applied to the general three-dimensional case, where long-range and short-range order is a function of the relevant thermodynamical parameters. In practice, long-range order will never be realized in a real crystal consisting of mosaic blocks which may behave as small subunits in order–disorder transitions. Another reason to assume partly incoherent areas in single crystals is the presence of possible strains or other distortions at the interfaces between domains which should cause a decrease of the averaged areas of coherent scattering. All these effects may lead to diffuse scattering in the neighbourhood of Bragg peaks, similar to the diffuse scattering caused by domain structures. For this reason an incoherent treatment of domains is probably more efficient, although considerable errors in intensity measurements may occur. A very careful study of line profiles is generally useful in order to decide between the various possibilities.

  • (8) Order–disorder in three dimensions

    Correlation functions in three dimensions may have very complicated periodicities; hence a careful study is necessary as to whether or not they may be interpreted in terms of a superlattice. If so, extinction rules have to be determined in order to obtain information on the superspace group. In the literature these are often called modulated structures because a sublattice, as determined by the basic lattice, and a superlattice may well be defined in reciprocal space: reflections of a sublattice including (000) are formally described by a multiplication by a lattice having larger lattice constants (superlattice) in reciprocal space; in real space this means a convolution with the Fourier transform of this lattice (sublattice). In this way the averaged structure is generated in each of the subcells (superposition or `projection' of all subcells into a single one). Obviously, the Patterson function of the averaged structure contains little information in the case of small subcells. Hence it is advisable to include the diffuse scattering of the superlattice reflections at the beginning of any structure determination.

    N subcells in real space are assumed, each of them representing a kind of a complicated `atom' which may be equal by translation or other symmetry operation. Once a superspace group has been determined, the usual extinction rules of space groups may be applied, remembering that the `atoms' themselves may have systematic extinctions. Major difficulties arise from the existence of different symmetries of the subgroup and the supergroup. Since the symmetry of the supergroup is lower in general, all missing symmetry elements may cause domains, corresponding to the missing symmetry element: translations cause antiphase domains in their generalized sense, other symmetry elements cause twins generated by rotations, mirror planes or the centre of symmetry. If all these domains are small enough to be detected by a careful study of line profiles, using diffraction methods with a high resolution power, the structural study may be facilitated by a reasonable explanation of scaling factors to be introduced for groups of reflections affected by the possible domain structures.

  • (9) Density modulations

    A density modulation of a structure in real space leads to pairs of satellites in reciprocal space. Each main reflection is accompanied by a pair of satellites in the directions [\mp \Delta {\bf H}] with phases [\mp 2 \pi \varphi]. The reciprocal lattice may then be written in the following form [(0 \leq \alpha \leq 1)]: [\eqalign{L({\bf H}) &+ {\alpha\over 2} L({\bf H} + \Delta {\bf H}) \exp \{2 \pi i \Phi\}\cr &+ {\alpha\over 2} L({\bf H} - \Delta {\bf H}) \exp \{2 \pi i (- \Phi)\}.}] Fourier transformation yields [\eqalignno{&l({\bf r}) \left[1 + {\alpha\over 2} \exp \{2 \pi i (\Delta {\bf H} \cdot {\bf r} + \Phi)\}\right.&\cr &\quad \left. + {\alpha\over 2} \exp \{- 2 \pi i (\Delta {\bf H} \cdot {\bf r} + \Phi)\}\right]&\cr &\quad\quad= l({\bf r}) [1 + \alpha \cos (2 \pi \Delta {\bf H} \cdot {\bf r} + \Phi)]. &(}] Equation ([link] describes the lattice modulated by a harmonic density wave. Since phases cannot be determined by intensity measurements, there is no possibility of obtaining any information on the phase relative to the sublattice. From ([link] it is obvious that the use of higher orders of harmonics does not change the situation. If [\Delta {\bf H}] is not rational such that [n \Delta {\bf H}] (n = integer) does not coincide with a main reflection in reciprocal space, the modulated structure is incommensurate with the basic lattice, and the phase of the density wave becomes meaningless. The same is true for the relative phases of the various orders of harmonic modulations of the density. This uncertainty even remains valid for commensurate density modulations of the sublattice, because coinciding higher-order harmonics in reciprocal space cause the same difficulty; higher-order coefficients cannot uniquely be separated from lower ones, consequently structure determination becomes impossible unless phase-determination methods are applied. Fortunately, density modulations of pure harmonic character are impossible for chemical reasons; they may be approximated by disorder phenomena for the averaged structure only. If diffuse scattering is taken into account the situation is changed considerably: A careful study of the diffuse scattering alone, although difficult in principle, will yield the necessary information about the relative phases of density waves (Korekawa, 1967[link]).

  • (10) Displacement modulations

    Displacement modulations are more complicated, even in a primitive structure. The Fourier transform of a longitudinal or a transverse displacement wave has to be calculated and this procedure does not result in a function of similar simplicity. A set of satellites is generated whose amplitudes are described by Bessel functions of νth order, where ν represents the order of the satellites. With α as amplitude of the displacement wave the intensity of the satellites increases with the magnitude of the product [H\alpha]. This means that a single harmonic displacement causes an infinite number of satellites. They may be unobservable at low diffraction angles as long as the amplitudes are small. If the displacement modulation is incommensurate there are no coincidences with reflections of the sublattice. Consequently, the reciprocal space is completely covered with an infinite number of satellites, or, in other words, with diffuse scattering. This is a clear indication that incommensurate displacement modulations belong to the category of disordered structures. Statistical fluctuations of amplitudes of the displacement waves cause additional diffuse scattering, regardless of whether the period is commensurate or incommensurate (Overhauser, 1971[link]; Axe, 1980[link]). Fluctuations of `phases' (= periods) cause a broadening of satellites in reciprocal space, but no change of their integrated intensities as long as the changes are not correlated with fluctuation periods. The broadening of satellite reflections increases with the order of satellites and [({\bf H} \cdot {\bf a})]. Obviously, there is no fundamental difference in the calculation of diffuse scattering with an ordered supercell of sufficient size.

    The use of optical transforms has been revived recently, although its efficiency is strongly dependent on the availability of a useful computer program capable of producing masks for optical diffraction. An atlas of optical transforms is available (Wooster, 1962[link]; Harburn et al., 1975[link]), but the possibility cannot be excluded that the diffuse scattering observed does not fit well into one of the diffraction pictures shown. Yet one of the major advantages of this optical method is the simple experimental setup and the high brilliance owing to the use of lasers. This method is specifically useful in disordered molecular structures where only a few orientations of the molecules have to be considered. It should be borne in mind, however, that all optical masks must correspond to projections of the disorder model along one specific direction which generates the two-dimensional diffraction picture under consideration. An important disadvantage is caused by the difficulty in simulating the picture of an atom. This situation may be improved by using computer programs with a high-resolution matrix printer representing electron densities by point densities of the printer. This latter method seems to be very powerful because of the possibility of avoiding `ghosts' in the diffraction picture. Guideline to solve a disorder problem

| top | pdf |

Generally, structure determination of a disordered crystal should start in the usual way by solving the average structure. The effectiveness of this procedure strongly depends on the distribution of integrated intensities of sharp and diffuse reflections. In cases where the integrated intensities of Bragg peaks is predominant, the maximum information can be drawn from the averaged structure. The observations of fractional occupations of lattice sites, split positions and anomalous temperature factors are indications of the disorder involved. Since these aspects of disorder phenomena in the averaged structure may be interpreted very easily, a detailed discussion of this matter is not given here (see any modern textbook of X-ray crystallography).

Difficulties may arise from the intensity integration which should be carried out using a high-resolution diffraction method. The importance of this may be understood from the following argument. The averaged structure is determined by the coherent superposition of different structure factors. This interpretation is true if there is a strictly periodic subcell with long-range order which allows for a clear separation of sharp and diffuse scattering. There are important cases, however, where this procedure cannot be applied without loss of information.

  • (a) The diffuse scattering is concentrated near the Bragg peaks for a large number of reflections. Because of the limited resolution power of conventional single-crystal methods the separation of sharp and diffuse scattering is impossible. Hence, the conventional study of integrated intensities does not really lead to an averaged structure. In this case a refinement should be tried using an incoherent superposition of different structure factors. Application of this procedure is subject to conditions which have to be checked very carefully before starting the refinement: first, it is necessary to estimate the amount of diffuse scattering not covered by intensity integration of the `sharp' reflections. Since loss in intensity, hidden in the background scattering, is underestimated very frequently, it should be checked whether nearly coinciding sharp and diffuse maxima are modulated by the same structure factor. It may be difficult to meet this condition in some cases; e.g. this condition is fulfilled for antiphase domains, but the same is not true for twin domains.

  • (b) The concentration of diffuse maxima near Bragg peaks is normally restricted to domain structures with a strictly periodic sublattice. Cases deviating from this rule are possible. Since they are rare, they are omitted here. Even structures with small deviations from the average structure do not lead to structure factors for diffuse scattering which are proportional to that of the average structure. This has been shown in the case of a twin structure correlated by a mirror plane where the reflections of a zone only have equal structure factors (Cowley & Au, 1978[link]). This effect causes even more difficulties for orthogonal lattices where the two twins have reflections in exactly the same positions, although differing in their structure factors. In this particular case, the incoherent or coherent treatment in refinements may be seriously hampered by strains originating from the boundary. Unsatisfactory refinements may be explained in this way but this does not improve their reliability.

The integrated intensity of any structure is independent of atomic positions if the atomic form factors remain unchanged by structural fluctuations. Small deviations of atomic form factors owing to electron-density changes of valence electrons are neglected. Consequently, the integrated diffuse intensities remain unchanged if the average structure is not altered by the degree of order. The latter condition is obeyed in cases where a geometrical long-range order of the lattice is independent of the degree of order, and no long-range order in the structure exists. This law is extremely useful for the interpretation of diffuse scattering. Unfortunately, intensity integration over coinciding sharp and diffuse maxima does not necessarily lead to a structure determination of the corresponding undistorted structure. This integration may be useful for antiphase domains without major structural changes at the boundaries. In all other cases the deviations of domains (or clusters) from the averaged structure determine the intensities of maxima which are no longer correlated with those of the average structure.

If the integrated intensity of diffuse scattering is comparable with, or even larger than, those of the Bragg peaks it is useful to begin the interpretation with a careful statistical study of the diffuse intensities. Intensity statistics can be applied in a way similar to the intensity statistics in classical structure determination. The following rules are briefly discussed in order to enable a semi-quantitative interpretation of the essential features of disorder to be realized.

  • (1) First, it is recommended that the integrated intensities be studied in certain areas of reciprocal space.

  • (2) Since low-angle scattering is very sensitive to fluctuations of densities, the most important information can be drawn from its intensity behaviour. If there is at least a one-dimensional sublattice in reciprocal space without diffuse scattering, it may often be concluded that there is no important low-angle scattering either. This law is subject to the condition of a sufficient number of reflections obeying this extinction rule without any exception.

  • (3) If the diffuse scattering shows maxima and minima, it should be checked whether the maxima observed may be approximately assigned to a lattice in reciprocal space. Obviously, this condition can hardly be met exactly if these maxima are modulated by a kind of structure factor, which causes displacements of maxima proportional to the gradient of this structure factor. Hence this influence may well be estimated from a careful study of the complete diffuse diffraction pattern.

    It should then be checked whether the corresponding lattice represents a sub- or a superlattice of the structure. An increase of the width of reflections as a function of growing [|{\bf H}|] indicates strained clusters of sub- or superlattice.

  • (4) The next step is the search for extinction rules of diffuse scattering. The simplest is the lack of low-angle scattering which has already been mentioned above. Since diffuse scattering is generally given by equation ([link] [\eqalign{I_{d} ({\bf H}) &= \langle |F ({\bf H})|^{2}\rangle - |\langle F ({\bf H})\rangle |^{2} \cr &= \textstyle\sum\limits_{\mu} p_{\mu} |F_{\mu} ({\bf H})|^{2} - |\textstyle\sum\limits_{\mu} p_{\mu} F_{\mu} ({\bf H})|^{2},}] it may be concluded that this condition is fulfilled in cases where all structural elements participating in disorder differ by translations only (stacking faults, antiphase domains etc.). They add phase factors to the various structure factors, which may become [n2\pi] (n = integer) for specific values of the reciprocal vector H. If all [p_{\mu}] are equivalent by symmetry: [p \textstyle\sum\limits_{\mu} |F_{\mu} ({\bf H})|^{2} - [p \textstyle\sum\limits_{\mu} F_{\mu} ({\bf H})] [p \textstyle\sum\limits_{\mu} F_{\mu}^{+} ({\bf H})] = 0.]

    Other possibilities of vanishing diffuse scattering may be derived in a similar manner for special reflections if glide operations are responsible for disorder. Since we are concerned with disordered structures, these glide operations need not necessarily be a symmetry operation of the lattice. It should be pointed out, however, that all these extinction rules of diffuse scattering are a kind of `anti'-extinction rule, because they are valid for reflections having maximum intensity for the sharp reflections unless the structure factor itself vanishes.

  • (5) Furthermore, it is important to plot the integrated intensities of sharp and diffuse scattering as a function of the reciprocal coordinates at least in a semi-quantitative way. If the ratio of integrated intensities remains constant in the statistical sense, we are predominantly concerned with a density phenomenon. It should be pointed out, however, that a particle-size effect of domains behaves like a density phenomenon (density changes at the boundary!).

    If the ratio of `diffuse' to `sharp' intensities increases with diffraction angle, we have to take into account atomic displacements. A careful study of this ratio yields very important information on the number of displaced atoms. The result has to be discussed separately for domain structures if the displacements are equal in the subcells of a single domain, but different for the various domains. In the case of two domains with displacements of all atoms the integrated intensities of sharp and diffuse reflections become statistically equal for large [|{\bf H}|]. Other rules may be derived from statistical considerations.

  • (6) The next step of a semi-quantitative interpretation is the check of the intensity distribution of diffuse reflections in reciprocal space. Generally this modulation is simpler than that of the sharp reflections. Hence it is frequently possible to start a structure determination with diffuse scattering. This method is extremely helpful for one- and two-dimensional disorder where partial structure determinations yield valuable information, even for the evaluation of the average structure.

  • (7) In cases where no sub- or superlattice belonging to diffuse scattering can be determined, a careful check of integrated intensities in the surroundings of Bragg peaks should again be performed. If systematic absences are found, the disorder is most probably restricted to specific lattice sites which may also be found in the average structure. The accuracy, however, is much lower here. If no such effects correlated with the average structure are observed, the disorder problem is related to a distribution of molecules or clusters with a structure differing from the average structure. As pointed out in Section[link] the problem of the representative structure(s) of the molecule(s) or the cluster(s) should be solved. Furthermore their distribution function(s) is (are) needed. In this particular case it is very useful to start with a study of diffuse intensity at low diffraction angles in order to acquire the information about density effects. Despite the contribution to sharp reflections, one should remember that the information derived from the average structure may be very low (e.g. small displacements, low concentrations etc.).

  • (8) As pointed out above, a Patterson picture – or strictly speaking a difference Patterson ([|\Delta F|^{2}]-Fourier synthesis) – may be very useful in this case. This method is promising in the case of disorder in molecular structures where the molecules concerned are at least partly known. Hence the interpretation of the difference Patterson may start with some internal molecular distances. Non-molecular structures show some distances of the average structure. Consequently a study of the important distances will yield information on displacements or replacements in the average structure. For a detailed study of this matter the reader is referred to the literature (Schwartz & Cohen, 1977[link]).

Although it is highly improbable that exactly the same diffraction picture will really be found, the use of an atlas of optical transforms (Wooster, 1962[link]; Harburn et al., 1975[link]; Welberry & Withers, 1987[link]) may be very helpful at the beginning of any study of diffuse scattering. The most important step is the separation of the distribution function from the molecular scattering. Since this information may be derived from a careful comparison of low-angle diffraction with the remaining sharp reflections, this task is not too difficult. If the influence of the distribution function is unknown, the reader is strongly advised to disregard the immediate neighbourhood of Bragg peaks in the first step of the interpretation. Obviously information may be lost in this way but, as has been shown in the past, much confusion caused by the attempt to interpret the scattering near the Bragg peaks with specific structural properties of a cluster or molecular model is avoided. The inclusion of this part of diffuse scattering can be made after the complete interpretation of the change of the influence of the distribution function on diffraction in the far-angle region.


Axe, J. D. (1980). Debye–Waller factors for incommensurate structures. Phys. Rev. B, 21, 4181–4190.Google Scholar
Cowley, J. M. & Au, A. Y. (1978). Diffraction by crystals with planar faults. III. Structure analysis using microtwins. Acta Cryst. A34, 738–743.Google Scholar
Gerold, V. (1954). Röntgenographische Untersuchungen über die Aushärtung einer Aluminium–Kupfer-Legierung mit Kleinwinkel-Schwenkaufnahmen. Z. Metallkd. 45, 593–607.Google Scholar
Guinier, A. (1942). Le mécanisme de la précipitation dans un cristal de solution solide métallique. Cas des systémes aluminum–cuivre et aluminum–argent. J. Phys. Radium, 8, 124–136.Google Scholar
Harburn, G., Taylor, C. A. & Welberry, T. R. (1975). An atlas of optical transforms. London: Bell.Google Scholar
Jagodzinski, H. & Haefner, K. (1967). On order–disorder in ionic non-stoichiometric crystals. Z. Kristallogr. 125, 188–200.Google Scholar
Jagodzinski, H. & Penzkofer, B. (1981). A new type of satellite in plagioclases. Acta Cryst. A37, 754–762.Google Scholar
Janner, A. & Janssen, T. (1980a). Symmetry of incommensurate crystal phases. I. Commensurate basic structures. Acta Cryst. A36, 399–408.Google Scholar
Janner, A. & Janssen, T. (1980b). Symmetry of incommensurate crystal phases. II. Incommensurate basic structures. Acta Cryst. A36, 408–415.Google Scholar
Korekawa, M. (1967). Theorie der Satellitenreflexe. Habilitationsschrift der Naturwissenschaftlichen Fakultät der Universität München, Germany.Google Scholar
Overhauser, A. W. (1971). Observability of charge-density waves by neutron diffraction. Phys. Rev. B, 3, 3173–3182.Google Scholar
Patterson, A. L. (1959). Fourier theory. In International tables for X-ray crystallography, Vol. II, ch. 2.5. Birmingham: Kynoch Press.Google Scholar
Schwartz, L. H. & Cohen, J. B. (1977). Diffraction from materials. New York: Academic Press.Google Scholar
Welberry, T. R. & Withers, R. L. (1987). Optical transforms of disordered systems displaying diffuse intensity loci. J. Appl. Cryst. 20, 280–288.Google Scholar
Wolff, P. M. de (1974). The pseudo-symmetry of modulated crystal structures. Acta Cryst. A30, 777–785.Google Scholar
Wolff, P. M. de, Janssen, T. & Janner, A. (1981). The superspace groups for incommensurate crystal structures with a one-dimensional modulation. Acta Cryst. A37, 625–636.Google Scholar
Wooster, W. A. (1962). Diffuse X-ray reflections from crystals, chs. IV, V. Oxford: Clarendon Press.Google Scholar
Zernike, F. & Prins, J. A. (1927). Die Beugung von Röntgenstrahlen in Flüssigkeiten als Effekt der Molekülanordnung. Z. Phys. 41, 184–194.Google Scholar

to end of page
to top of page