International
Tables for
Crystallography
Volume B
Reciprocal space
Edited by U. Shmueli

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

Section 4.2.4.  Quantitative interpretation

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.4. Quantitative interpretation

| top | pdf |

4.2.4.1. Introduction

| top | pdf |

In these sections quantitative interpretations of the elastic part of diffuse scattering (X-rays and neutrons) are outlined. Although similar relations are valid, magnetic scattering of neutrons is excluded. Obviously, all disorder phemomena are strongly temperature dependent if thermal equilibrium is reached. Consequently, the interpretation of diffuse scattering should include a statistical thermodynamical treatment. Unfortunately, no quantitative theory for the interpretation of structural phenomena is so far available: all quantitative solutions introduce formal order parameters such as correlation functions or distributions of defects. At low temperatures (low concentration of defects) the distribution function plays the dominant role in diffuse scattering. With increasing temperature the number of defects increases with corresponding strong interactions between them. Therefore, correlations become increasingly important, and phase transformations of first or higher order may occur which need a separate theoretical treatment. In many cases large fluctuations of structural properties occur which are closely related to the dynamical properties of the crystal. Theoretical approximations are possible but their presentation is far beyond the scope of this article. Hence we restrict ourselves to formal parameters in the following.

Point defects or limited structural units, such as molecules, clusters of finite size etc., may only be observed in diffraction for a sufficiently large number of defects. This statement is no longer true in high-resolution electron diffraction where single defects may be observed either by diffraction or by optical imaging if their contrast is high enough. Hence, electron microscopy and diffraction provide valuable methods for the interpretation of disorder phenomena.

The arrangement of a finite assembly of structural defects is described by its structure and its three-dimensional (3D) distribution function. Structures with a strict 1D periodicity (chain-like structures) need a 2D distribution function, while for structures with a 2D periodicity (layers) a 1D distribution function is sufficient. Since the distribution function is the dominant factor in statistics with correlations between defects, we define the dimensionality of disorder as that of the corresponding distribution function. This definition is more effective in diffraction problems because the dimension of the disorder problem determines the dimension of the diffuse scattering: 1D diffuse streaks, 2D diffuse layers, or a general 3D diffuse scattering.

Strictly speaking, completely random distributions cannot be realized as shown in Section 4.2.3.[link] They occur approximately if the following conditions are satisfied.

  • (1) The average volume of all defects including their surrounding strain fields [NcV_{d}] (N = number of unit cells, c = concentration of defects, Vd = volume of the defect with [V_{d}> V_{c}, V_{c}] = volume of the unit cell) is small in comparison with the total volume [NV_{c}] of the crystal, or [V_{c} \gg cV_{d}].

  • (2) Interactions between the defects are negligible. These conditions, however, are valid in very rare cases only, i.e. where small concentrations and vanishing strain fields are present. Remarkable exceptions from this rule are real point defects without interactions, such as isotope distribution (neutron diffraction!), or the system AuAg at high temperature.

As already mentioned, disorder phenomena may be observed in thermal equilibrium. Two completely different cases have to be considered.

  • (1) The concentration of defects is given by the chemical composition, i.e. impurities in a closed system.

  • (2) The number of defects increases with temperature and also depends on pressure or other parameters, i.e. interstitials, voids, static displacements of atoms, stacking faults, dislocations etc.

In many cases the defects do not occur in thermal equilibrium. Nevertheless, their diffuse scattering is temperature dependent because of the anomalous thermal movements at the boundary of the defect. Hence, the observation of a temperature-dependent behaviour of diffuse scattering cannot be taken as a definite criterion of thermal equilibrium without further careful interpretation.

Ordering of defects may take place in a very anisotropic manner. This is demonstrated by the huge number of examples of 1D disorder. As shown by Jagodzinski (1963)[link] this type of disorder cannot occur in thermal equilibrium for the infinite crystal. This type of disorder is generally formed during crystal growth or mechanical deformation. Similar arguments may be applied to 2D disorder. This is a further reason why the so-called Ising model can hardly be used in order to obtain interaction energies of structural defects. From these remarks it becomes clear that order parameters are more or less formal parameters without strict thermodynamical meaning.

The following section is organized as follows: first we discuss the simple case of 1D disorder where reliable solutions of the diffraction problem are available. Intensity calculations of diffuse scattering of 2D disorder by chain-like structures follow. Finally, the 3D case is treated, where formal solutions of the diffraction problem have been tried and applied successfully to metallic systems to some extent. A short concluding section concerns the special phenomenon of orientational disorder.

4.2.4.2. One-dimensional disorder of ordered layers

| top | pdf |

As has been pointed out above, it is often useful to start the interpretation of diffuse scattering by checking the diffraction pattern with respect to the dimensionality of the disorder concerned. Since each disordered direction in the crystal demands a violation of the corresponding Laue condition, this question may easily be answered by looking at the diffuse scattering. Diffuse streaks in reciprocal space are due to a one-dimensional violation of the Laue conditions, and will be called one-dimensional disorder. This kind of order is typical for layer structures, but it is frequently observed in cases where several sequences of layers do not differ in the interactions of next-nearest neighbours. Typical examples are structures which may be described in terms of close packing, e.g. hexagonal and cubic close packing.

For a quantitative interpretation of diffuse streaks we need one-dimensional correlation functions, which may uniquely be determined if a single independent correlation function is active. According to equation (4.2.3.28)[link] Fourier transformation yields the information required. In all other cases a specific model has to be suggested for a full interpretation of diffuse streaks. Another comment seems to be necessary: disorder parameters can be defined uniquely only if the diffraction pattern allows for a differentiation between long-range and short-range order. This question can at least partly be answered by studying the line width of sharp reflections with a very good resolution. Since integrated intensities of sharp reflections have to be separated from the diffuse scattering, this question is of outstanding importance in most cases. Inclusion of diffuse parts in the diffraction pattern during intensity integration of sharp reflections may lead to serious errors in the interpretation of the average structure.

The existence of diffuse streaks in more than one direction of reciprocal space means that the diffraction problem is no longer one-dimensional. Sometimes the problem may be treated independently, if the streaks are sharp, and no interference effects may be observed in the diffraction pattern in areas where the diffuse streaks do overlap. In all other cases there are correlations between the various directions of one-dimensional disorder which may be determined with the aid of a model covering more than one of the pertinent directions of disorder.

Before starting the discussion of the quantitative solution of the one-dimensional problem, some remarks should be made on the usefulness of quantitative disorder parameters. It is well known from statistical thermodynamics that a one-dimensional system cannot show long-range order above [T = 0\hbox{ K}]. Obviously, this statement is in contradiction with many experimental observations where long-range order is realized even in layer structures. The reason for this behaviour is given by the following arguments which are valid for any structure. Let us assume a structure with strong interactions at least in two directions. From the theoretical treatment of the two-dimensional Ising model it is known that such a system shows long-range order below a critical temperature [T_{c}]. This statement is true even if the layer is finite, although the strict thermodynamic behaviour is not really critical in the thermodynamical sense. A three-dimensional crystal can be constructed by adding layer after layer. Since each layer has a typical two-dimensional free energy, the full statistics of the three-dimensional crystal may be calculated by introducing a specific free energy for the various stackings of layers. Obviously, this additional energy has to include terms describing potential and entropic energies as well. They may be formally developed into contributions of next, overnext etc. nearest neighbours. Apparently, the contribution to entropy must include configurational and vibrational parts which are strongly coupled. As long as the layers are finite, there is a finite probability of a fault in the stacking sequence of layers which approaches zero with increasing extension of the layers. Consequently, the free energy of a change in the favourite stacking sequence becomes infinite quadratically with the size of the layer. Therefore, the crystal should be either completely ordered or disordered; the latter case can only be realized if the free energies of one or more stacking sequences are exactly equal (very rare, but possible over a small temperature range of phase transformations). An additional positive entropy associated with a deviation from the periodic stacking sequence may lead to a kind of competition between entropy and potential energy, in such a way that periodic sequences of faults result. Obviously, this situation occurs in the transition range of two structures differing only in their stacking sequence. On the other hand, one must assume that defects in the stacking sequence may be realized if the size of the layers is small. This situation occurs during crystal growth, but one should remember that the number of stacking defects should decrease with increasing size of the growing crystal. Apparently, this rearrangement of layers may be suppressed as a consequence of relaxation effects. The growth process itself may influence the propagation of stacking defects and, consequently, the determination of stacking-fault probabilities, aiming at an interpretation of the chemical bonding seems to be irrelevant in most cases.

The quantitative solution of the diffraction problem of one-dimensional disorder follows a method similar to the Ising model. As long as next-nearest neighbours alone are considered, the solution is very simple only if two possibilities of structure factors are to be taken into account. Introducing the probability of equal pairs 1 and 2, α, one arrives at the known solution for the a priori probability [p_{\mu}] and a posteriori probabilities [p_{\mu \mu'}({\bf m})], respectively. In the one-dimensional Ising model with two spins and the interaction energies [(U \mp \Delta U) / k_{B}T], defining the pair probability [[\alpha = p_{11} ({\bf 1})]] [\alpha = {\exp \{\pm \Delta U / k_{B}T\}\over [\exp \{+ \Delta U / k_{B}T\} + \exp \{- \Delta U / k_{B}T\}]}] the full symmetry is [p_{1} = p_{2} = {1\over 2}], and [p_{11} ({\bf m}) = p_{22} ({\bf m})].

Consequently: [p_{12} ({\bf m}) = p_{22} ({\bf m}) = 1 - p_{11} ({\bf m}).] The scattered intensity is given by [{I} ({\bf H}) = L (h, k) \textstyle\sum\limits_{{\bf m}} \langle FF_{m}^{+}\rangle (N - |m|) \exp \{- 2\pi iml\}, \eqno(4.2.4.1)] where [{\bf m} = m{\bf c}], [N = \hbox{number of unit cells}] in the c direction and [\langle FF^{+}\rangle] depends on [\lambda_{1}], [\lambda_{2}] which are the eigenvalues of the matrix [\pmatrix{\alpha &1 - \alpha\cr 1 - \alpha &\alpha\cr}.] From the characteristic equation [\lambda^{2} - 2\alpha \lambda - 1 + 2\alpha = 0 \eqno(4.2.4.2)] one has [\lambda_{1} = 1\hbox{;} \quad \lambda_{2} = 2\alpha - 1. \eqno(4.2.4.2a)] [\lambda_{1}] describes a sharp Bragg reflection (average structure) which need not be calculated. Its intensity is simply proportional to [\langle F({\bf H})\rangle]. The second characteristic value yields a diffuse reflection in the same position if the sign is positive [(\alpha\gt 0.5)], and in a position displaced by [{1\over 2}] in reciprocal space if the sign is negative [(\alpha\lt 0.5)]. Because of the symmetry conditions [p_{11} ({\bf m})] only is needed; it may be determined with the aid of the boundary conditions [p_{11} ({\bf 0}) = 1, \quad p_{11} ({\bf 1}) = \alpha,] and the general relation that [p_{\mu \mu'} ({\bf m})] is given by [p_{\mu \mu'} ({\bf m}) = c'_{\mu \mu'} \lambda_{1}^{m} + c''_{\mu \mu'} \lambda_{2}^{m}.] The final solution of our problem yields simply: [\eqalign{p_{11} ({\bf m}) &= {\textstyle{1\over 2}} + {\textstyle{1\over 2}} \lambda_{2}^{m} = p_{22} ({\bf m}),\cr p_{12} ({\bf m}) &= {\textstyle{1\over 2}} - {\textstyle{1\over 2}} \lambda_{2}^{m} = p_{21} ({\bf m}).}] The calculation of the scattered intensity is now performed with the general formula [\eqalignno{I({\bf H}) &= L(h, k) \textstyle\sum\limits_{m} \textstyle\sum\limits_{\mu, \, \mu'} p_{\mu} p_{\mu \mu'} ({\bf m}) F_{\mu} F_{\mu'}^{+} (N - |m|) &\cr &\quad \times \exp \{- 2 \pi iml\}. &(4.2.4.3)}] Evaluation of this expression yields [\eqalignno{I({\bf H}) &= L(h, k) \textstyle\sum\limits_{m} (N - |m|) \exp \{- 2 \pi iml\} &\cr &\quad \times [|{\textstyle{1\over 2}} (F_{1} + F_{2})|^{2} \lambda_{1}^{m} + |{\textstyle{1\over 2}} (F_{1} - F_{2})|^{2} \lambda_{2}^{m}]. &(4.2.4.4)}] Since the characteristic solutions of the problem are real: [\eqalignno{I({\bf H}) &= L({\bf h}) |(F_{1} + F_{2})/2|^{2} &\cr &\quad + L(h, k) |(F_{1} - F_{2})/2|^{2} &\cr &\quad \times {1 - |\lambda_{2}|^{2}\over 1 - 2 \lambda_{2} \cos 2 \pi l + |\lambda_{2}|^{2}}. &(4.2.4.5)}] This particle size effect has been neglected in (4.2.4.5)[link]. This result confirms the fact mentioned above that the sharp Bragg peaks are determined by the averaged structure factor and the diffuse one by its mean-square deviation.

For the following reason there are no examples for quantitative applications: two different structures generally have different lattice constants; hence the original assumption of an undisturbed lattice geometry is no longer valid. The only case known to the authors is the typical lamellar structure of plagioclases, reported by Jagodzinski & Korekawa (1965)[link]. The authors interpret the well known `Schiller effect' as a consequence of optical diffraction. Hence, the size of the lamellae is of the order of 2000 Å. This long-period superstructure cannot be explained in terms of next-nearest-neighbour interactions. In principle, however, the diffraction effects are similar: instead of the diffuse peak as described by the second term in equation (4.2.4.5)[link], satellites of first and second order, etc. accompanying the Bragg peaks are observed. The study of this phenomenon (Korekawa & Jagodzinski, 1967[link]) has not so far resulted in a quantitative interpretation.

Obviously, the symmetry relation used in the formulae discussed above is only valid if the structures described by the [F_{\mu}] are related by symmetries such as translations, rotations or combinations of both. The type of symmetry has an important influence on the diffraction pattern.

  • (1) Translation parallel to the ordered layers

    If the translation vector between the two layers in question is such that [2 \Delta {\bf r}] is a translation vector parallel to the layer, there are two relevant structure factors [F_{1}, F_{2} = F_{1} \exp (2 \pi i{\bf H} \cdot \Delta {\bf r}).] [{\bf H} \cdot \Delta {\bf r}] may be either an integer or an integer [+\ {\textstyle{1\over 2}}]. Since any integer may be neglected because of the translation symmetry parallel to the layer, we have [F_{1} = F_{2}] in the former case, and [F_{1} = - F_{2}] in the latter. As a consequence either the sharp reflections given in equation (4.2.4.4)[link] vanish, or the same is true for the diffuse ones. Hence, the reciprocal lattice may be described in terms of two kinds of lattice rows, sharp and diffuse, parallel to the reciprocal coordinate l.

    Disorder of this type is observed very frequently. One of the first examples was wollastonite, CaSiO3, published by Jefferey (1953)[link]. Here the reflections with [k = 2n] are sharp Bragg peaks without any diffuse scattering. Diffuse streaks parallel to (h00), however, are detected for [k = 2n + 1]. In the light of the preceding discussion, the translation vector is [{1\over 2}{\bf b}], and the plane of ordered direction (plane of intergrowth of the two domains) is (100). Hence the displacement is parallel to the said plane. Since the intensity of the diffuse lines does not vary according to the structure factor involved, the disorder cannot be random. The maxima observed are approximately in the position of a superstructure, generated by large domains without faults in the stacking sequence, mutually displaced by [{1\over 2}{\bf b}] (antiphase domains). This complicated ordering behaviour is typical for 1D order and may easily be explained by the above-mentioned fact that an infinitely extended interface between two domains causes an infinite unfavourable energy (Jagodzinski, 1964b[link], p. 188). Hence, a growing crystal should become increasingly ordered. This consideration explains why the agreement between a 1D disorder theory and experiment is often so poor.

    Examples where more than one single displacement vector occur are common. If these are symmetrically equivalent all symmetries have to be considered. The most important cases of displacements differing only by translation are the well known close-packed structures (see below). A very instructive example is the mineral maucherite (approximately Ni4As3). According to Jagodzinski & Laves (1947)[link] the structure has the following disorder parameters: interface (001), displacement vectors [000], [[{1\over 2}00]], [[0{1\over 2}0]], [[{1\over 2}{1\over 2}0]]. From equation (4.2.4.5)[link] we obtain: [\langle F({\bf H})\rangle = [1 + \exp \{\pi ih\} + \exp \{\pi ik\} + \exp \{\pi i(h + k)\}]/4.] Hence there are sharp reflections for [h, k = \hbox{ even}], and diffuse ones otherwise. Further conclusions may be drawn from the average structure.

  • (2) Translation perpendicular to the ordered layers

    If the translation is c/2 the structure factors are: [\eqalign{F_{2} &= F_{1} \exp \{2 \pi il\}\cr F_{1} &= F_{2} \hbox{ for } l = \hbox{even}\cr F_{1} &= - F_{2} \hbox{ for } l = \hbox{odd}.}] There are sharp [(l = 2n)] and diffuse [(l = 2n + 1)] reflections on all reciprocal-lattice rows discussed above.

    Since the sharp and diffuse reflections occur on the same reciprocal line there is a completely different behaviour compared with the preceding case. In general, a component of any displacement vector perpendicular to the interface gives rise to a change in chemical composition as shown in the next example: in a binary system consisting of A and B atoms with a tendency towards an alternating arrangement of A and B layers, any fault in the sequence BABAB|BABAB|B increases the number of B atoms (or A atoms). Generally such kinds of defects will lead to an interface with a different lattice constant, at least in the direction perpendicular to the interface. Consequently the exact displacement vectors of [{{1} \over {2}}], [{{1} \over {3}}], [{{1} \over {4}}] are rare. Since ordered structures should be realized in the 1D case, incommensurate superstructures will occur which are very abundant during ordering processes. An interesting example has been reported and interpreted by Cowley (1976a[link],b[link]) where the displacement vector has a translational period of [{{1} \over {4}}] perpendicular to the plane of intergrowth. Reflections 00l and 22l with [l = 4n] are sharp, all remaining reflections more or less diffuse. Since the maxima [(\bar{1}11), (\bar{1}33)] show a systematically different behaviour, there is also a displacement component parallel to the plane of intergrowth in question. The semi-quantitative interpretation has been given in his paper.

  • (3) Rotations

    The discussion concerning layers related by a twofold rotation parallel to c may easily be made by simply considering their structure factors. Since the layers do not obey the twofold symmetry their structure factors are generally different; unless they equalize accidentally there are sharp and diffuse reflections according to the values of [\langle F\rangle] and [\Delta F], respectively. Obviously, [F_{1} = F_{2}] is valid only if [h = k = 0]; consequently there is just one reciprocal-lattice row free from diffuse scattering.

  • (4) Asymmetric case

    In the asymmetric case the symmetry conditions used above are no longer valid: [p_{1} \neq p_{2}, p_{12} \neq p_{21}, p_{11} \neq p_{22}.] But there is one condition which may be derived from the invariance of the numbers of pairs in the relevant and its opposite direction: [p_{\mu} p_{\mu \mu'} ({\bf m}) = p_{\mu'} p_{\mu' \mu} (- {\bf m}).] This equation requires that [p_{\mu \mu'}({\bf m})] is not necessarily symmetric in m. The calculation of characteristic values yields [\lambda_{1} = 1, \qquad \lambda_{2} = (\alpha_{1} + \alpha_{2}) - 1. \eqno(4.2.4.6)] The a priori probabilities are now different from [{{1} \over {2}}], and may be calculated by considering [p_{\mu \mu'}({\bf m}) \rightarrow p_{\mu'} ({\bf m} \rightarrow \infty)]: [p_{1} = \alpha_{1}/(\alpha_{1} + \alpha_{2})\hbox{;} \qquad p_{2} = \alpha_{2}/(\alpha_{1} + \alpha_{2}).] The intensity is given by an expression very similar to (4.2.4.5)[link]: [\eqalignno{{I}({\bf H}) &= L({\bf h})|[\alpha_{1}/(\alpha_{1} + \alpha_{2})] F_{1} + [\alpha_{2}/(\alpha_{1} + \alpha_{2})] F_{2}|^{2} &\cr &\quad+ L(h, k)|[\alpha_{1}/(\alpha_{1} + \alpha_{2})] F_{1} - [\alpha_{2}/(\alpha_{1} + \alpha_{2})] F_{2}|^{2} &\cr &\quad\times (1 - |\lambda_{2}|^{2})/(1 - 2\lambda_{2} \cos 2\pi l + |\lambda_{2}|^{2}). &(4.2.4.7)}] Again there are sharp Bragg reflections and diffuse ones in the same positions, or in a displaced position depending on the sign of [\lambda_{2}].

    From a discussion of the next-nearest-neighbour Ising model one may conclude that the detailed study of the qualitative behaviour of sharp and diffuse reflections may give additional information on the symmetry of the layers involved.

    In the case of translations between neighbouring layers not fulfilling the condition [{\bf h} \cdot {\bf r}] = integer, where r is parallel to the layer, more than two structure factors have to be taken into account. If [n{\bf h} \cdot {\bf r}] = integer, where n is the smallest integer fulfilling the said condition, n different structure factors have to be considered. The characteristic equation has formally to be derived with the aid of an [n \times n] matrix containing internal symmetries which may be avoided by adding the phase factors [\varepsilon = \exp \{2 \pi i{\bf H} \cdot {\bf r}/n\}], [\varepsilon^{+} = \exp \{-2 \pi i{\bf H} \cdot {\bf r}/n\}] to the probability of pairs. The procedure is allowed if the displacements r and −r are admitted only for neighbouring layers. The matrix yielding the characteristic values may then be reduced to [\pmatrix{\alpha_{1}\varepsilon &(1 - \alpha_{1})\varepsilon^{+}\cr (1 - \alpha_{2})\varepsilon &\alpha_{2}\varepsilon^{+}\cr}] and yields the characteristic equation [\lambda^{2} - \lambda (\alpha_{1}\varepsilon + \alpha_{2}\varepsilon^{+}) - 1 + \alpha_{1} + \alpha_{2} = 0. \eqno(4.2.4.8)] Equation (4.2.4.8)[link] gives sharp Bragg reflections for [{\bf H} \cdot {\bf r}/n] = integer; the remaining diffuse reflections are displaced corresponding to the phase of the complex characteristic value. Equation (4.2.4.8)[link] has been used in many cases. Qualitative examples are the mixed-layer structures published by Hendricks & Teller (1942)[link]. An example of a four-layer-type structure is treated by Dubernat & Pezerat (1974)[link]. A first quantitative treatment with good agreement between theory and experimental data (powder diffraction) has been given by Dorner & Jagodzinski (1972)[link] for the binary system TiO2–SnO2. In the range of the so-called spinodal decomposition the chemical composition of the two domains and the average lengths of the two types of domains could be determined. Another quantitative application was reported by Jagodzinski & Hellner (1956)[link] for the transformation of [\hbox{RhSn}_{2}] into a very complicated mixed-layer type. A good agreement of measured and calculated diffuse scattering (asymmetric line profiles, displacement of maxima) could be found over a wide angular range of single-crystal diffraction.

4.2.4.2.1. Stacking disorder in close-packed structures

| top | pdf |

From an historical point of view stacking disorder in close-packed systems is most important. The three relevant positions of ordered layers are represented by the atomic coordinates [|0, 0|, |{{1} \over {3}}, {{2} \over {3}}|, |{{2} \over {3}}, {{1} \over {3}}|] in the hexagonal setting of the unit cell, or simply by the figures 1, 2, 3 in the same sequence. Structure factors [F_{1}, F_{2}, F_{3}] refer to the corresponding positions of the same layer: [\eqalign{F_{2} &= F_{1} \exp \{2 \pi i(h - k)/3\},\cr F_{3} &= F_{1} \exp \{- 2 \pi i(h - k)/3\},}] hence [F_{1} = F_{2} = F_{3}\quad \hbox{if } h - k \equiv 0\ (\hbox{mod } 3).] According to the above discussion the said indices define the reciprocal-lattice rows exhibiting sharp reflections only, as long as the distances between the layers are exactly equal. The symmetry conditions caused by the translation are normally: [\eqalign{&p_{1} = p_{2} = p_{3},\qquad p_{11} = p_{22} = p_{33}, \cr &p_{12} = p_{23} = p_{31},\quad p_{13} = p_{21} = p_{32}.}] For the case of close packing of spheres and some other problems any configuration of m layers determining the a posteriori probability [p_{\mu \mu'} ({\bf m}), \mu = \mu'], has a symmetrical counterpart where μ is replaced by [\mu' + 1] (if [\mu' = 3, \mu' + 1 = 1]).

In this particular case [p_{12}({\bf m}) = p_{13}({\bf m})], and equivalent relations generated by translation.

Nearest-neighbour interactions do not lead to an ordered structure if the principle of close packing is obeyed (no pairs in equal positions) (Hendricks & Teller, 1942[link]; Wilson, 1942[link]). Extension of the interactions to next-but-one or more neighbours may be carried out by introducing the method of matrix multiplication developed by Kakinoki & Komura (1954[link], 1965[link]), or the method of overlapping clusters (Jagodzinski, 1954[link]). The latter procedure is outlined in the case of interactions between four layers. A given set of three layers may occur in the following 12 combinations: [\eqalignno{&123, 231, 312\hbox{;}\quad 132, 213, 321\hbox{;}\cr &121, 232, 313\hbox{;}\quad 131, 212, 323.}] Since three of them are equivalent by translation, only four representatives have to be introduced: [123\hbox{;} \quad 132\hbox{;} \quad 121\hbox{;} \quad 131.] In the following the new indices 1, 2, 3, 4 are used for these four representatives for the sake of simplicity.

In order to construct the statistics layer by layer the next layer must belong to a triplet starting with the same two symbols with which the preceding one ended, e.g. 123 can only be followed by 231, or 232. In a similar way 132 can only be followed by 321 or 323. Since both cases are symmetrically equivalent the probabilities [\alpha_{1}] and [1 - \alpha_{1}] are introduced. In a similar way 121 may be followed by 212 or 213 etc. For these two groups the probabilities [\alpha_{2}] and [1 - \alpha_{2}] are defined. The different translations of groups are considered by introducing the phase factors as described above. Hence, the matrix for the characteristic equation may be set up as follows. As representative cluster of each group is chosen that one having the number 1 at the centre, e.g. 312 is representative for the group 123, 231, 312; in a similar way 213, 212 and 313 are the remaining representatives. Since this arrangement of three layers is equivalent by translation, it may be assumed that the structure of the central layer is not influenced by the statistics to a first approximation. The same arguments hold for the remaining three groups. On the other hand, the groups 312 and 213 are equivalent by rotation only. Consequently their structure factors may differ if the influence of the two neighbours has to be taken into account. A different situation exists for the groups 212 and 313 which are correlated by a centre of symmetry, which causes different corresponding structure factors. It should be pointed out, however, that the structure factor is invariant as long as there is no influence of neighbouring layers on the structure of the central layer. The latter is often observed in close-packed metal structures, or in compounds like ZnS, SiC and others. For the calculation of intensities [p_{\mu} p_{\mu \mu'}, F_{\mu} F_{\mu'}^{+}] is needed.

According to the following scheme of sequences any sequence of pairs is correlated with the same phase factor for [FF^{+}] due to translation, if both members of the pair belong to the same group. Consequently the phase factor may be attached to the sequence probability such that [FF^{+}] remains unchanged, and the group may be treated as a single element in the statistics. In this way the reduced matrix for the solution of the characteristic equation is given by

[F_{\mu}][F_{\mu'}^{+}]
(1)(2)(3)(4)
[312, 123 (\varepsilon^{+}), 231 (\varepsilon)][212, 323 (\varepsilon^{+}), 131 (\varepsilon)][213, 321 (\varepsilon^{+}), 132 (\varepsilon)][313, 121 (\varepsilon^{+}), 232 (\varepsilon)]
(1) [312, 123 (\varepsilon), 231 (\varepsilon^{+})] [\alpha_{1} \varepsilon^{+}] 0 0 [(1 - \alpha_{1}) \varepsilon^{+}]
(2) [212, 323 (\varepsilon), 131 (\varepsilon^{+}) ] [(1 - \alpha_{2}) \varepsilon^{+} ] 0 0 [\alpha_{2} \varepsilon^{+}]
(3) [213, 321 (\varepsilon), 132 (\varepsilon^{+})] 0 [(1 - \alpha_{2}) \varepsilon] [\alpha_{1} \varepsilon] 0
(4) [313, 121 (\varepsilon), 232 (\varepsilon^{+})] 0 [\alpha_{2} \varepsilon] [(1 - \alpha_{2})\varepsilon] 0

There are three solutions of the diffraction problem:

  • (1) If [h - k = 0] (mod 3), [\varepsilon = +1], there are two quadratic equations: [\eqalign{\lambda^{2} - (\alpha_{1} + \alpha_{2})\lambda - 1 + \alpha_{1} + \alpha_{2} &= 0\cr \lambda^{2} - (\alpha_{1} - \alpha_{2})\lambda + 1 - \alpha_{1} - \alpha_{2} &= 0\cr} \eqno(4.2.4.9)] with solutions [\displaylines{\lambda_{1} = 1, \qquad \lambda_{2} = \alpha_{1} + \alpha_{2} - 1\cr \hfill \lambda_{3/4} = {{\alpha_{1} - \alpha_{2}} \over {2}} \pm \left[{{(\alpha_{1} - \alpha_{2})^{2}} \over {4}} - 1 + \alpha_{1} - \alpha_{2}\right]^{1/2}. \hfill (4.2.4.10)}]

    [\lambda_{1}] and [\lambda_{2}] are identical with the solution of the asymmetric case of two kinds of layers [cf. equation (4.2.4.6)].[link] They yield sharp reflections for l = integer, and diffuse ones in a position determined by the sequence probabilities [\alpha_{1}] and [\alpha_{2}] (position either l = integer, l = ½ +integer, respectively). The remaining two characteristic values may be given in the form [\lambda = |\lambda| \exp \{2 \pi i\varphi\}], where ϕ determines the position of the reflection. If the structure factors of the layers are independent of the cluster, [\lambda_{2}, \lambda_{3}, \lambda_{4}] become irrelevant because of the new identity of the F's (no diffuse scattering). Weak diffuse intensities on the lattice rows [k - k = 0] (mod 3) may be explained in terms of this influence.

  • (2) The remaining two solutions for [\varepsilon = \exp \{\pm 2 \pi i(h - k)/3\}] are equivalent, and result in the same characteristic values. They have been discussed explicitly in the literature; the reader is referred to the papers of Jagodzinski (1949a[link],b[link],c[link], 1954[link]).

In order to calculate the intensities one has to reconsider the symmetry of the clusters, which is different from the symmetry of the layers. Fortunately, a threefold rotation axis is invariant against the translations, but this is not true for the remaining symmetry operations in the layer if there are any more. Since we have two pairs of inequivalent clusters, namely 312, 213 and 212, 313, there are only two different a priori probabilities [p_{1} = p_{3}] and [p_{2} = p_{4} = {{1} \over {2}} (1 - 2p_{1})].

The symmetry conditions of the new clusters may be determined by means of the so-called `probability trees' described by Wilson (1942)[link] and Jagodzinski (1949b[link], pp. 208–214). For example: [p_{11} = p_{33}], [p_{22} = p_{44}], [p_{13} = p_{31}], [p_{24} = p_{42}] etc.

It should be pointed out that clusters 1 and 3 describe a cubic arrangement of three layers in the case of simple close packing, while clusters 2 and 4 represent the hexagonal close packing. There may be a small change in the lattice constant c perpendicular to the layers. Additional phase factors then have to be introduced in the matrix for the characteristic equation, and a recalculation of the constants is necessary. As a consequence, the reciprocal-lattice rows [(h - k) \equiv 0\;(\hbox{mod } 3)] become diffuse if [l \neq 0], and the diffuseness increases with l. A similar behaviour results for the remaining reciprocal-lattice rows.

The final solution of the diffraction problem results in the following general intensity formula: [\eqalignno{I({\bf H}) &= L(h, k) N \textstyle\sum\limits_{\nu} \{A_{\nu}({\bf H}) (1 - |\lambda_{\nu}|^{2}) &\cr &\quad \times [1 - 2|\lambda_{\nu}| \cos 2 \pi (l - \varphi_{\nu}) + |\lambda_{\nu}|^{2}]^{-1} &\cr &\quad - 2 B_{\nu}({\bf H})|\lambda_{\nu}| \sin 2 \pi (l - \varphi_{\nu}) &\cr &\quad \times [1 - 2|\lambda_{\nu}| \cos 2 \pi (l - \varphi_{\nu}) + |\lambda_{\nu}|^{2}]^{-1}\}. &(4.2.4.11)}] Here [A_{\nu}] and [B_{\nu}] represent the real and imaginary part of the constants to be calculated with the aid of the boundary conditions of the problem. The first term in equation (4.2.4.11)[link] determines the symmetrical part of a diffuse reflection with respect to the maximum, and is completely responsible for the integrated intensity. The second term causes an antisymmetrical contribution to intensity profiles but does not influence the integrated intensities. These general relations enable a semi-quantitative interpretation of the sharp and diffuse scattering in any case, without performing the time-consuming calculations of the constants which may only be done in more complicated disorder problems with the aid of a computer program evaluating the boundary conditions of the problem.

This can be carried out with the aid of the characteristic values and a linear system of equations (Jagodzinski, 1949a[link],b[link],c[link]), or with the aid of matrix formalism (Kakinoki & Komura, 1954[link]; Takaki & Sakurai, 1976[link]). As long as only the line profiles and positions of the reflections are required, these quantities may be determined experimentally and fitted to characteristic values of a matrix. The size of this matrix is given by the number of sharp and diffuse maxima observed, while [|\lambda_{\nu}|] and [\exp\{2 \pi i \varphi_{\nu}\}] may be found by evaluating the line width and the position of diffuse reflections. Once this matrix has been found, a semi-quantitative model of the disorder problem can be given. If a system of sharp reflections is available, the averaged structure can be solved as described in Section 4.2.3.2[link]. The determination of the constants of the diffraction problem is greatly facilitated by considering the intensity modulation of diffuse scattering, which enables a phase determination of structure factors to be made under certain conditions.

The theory of closed-packed structures with three equivalent translation vectors has been applied very frequently, even to systems which do not obey the principle of close-packing. The first quantitative explanation was published by Halla et al. (1953)[link]. It was shown there that single crystals of C18H24 from the same synthesis may have a completely different degree of order. This was true even within the same crystal. Similar results were found for C, Si, CdI2, CdS2, mica and many other compounds. Quantitative treatments are less abundant [e.g. CdI2: Martorana et al. (1986)[link]; MX3 structures: Conradi & Müller (1986)[link]]. Special attention has been paid to the quantitative study of polytypic phase transformations in order to gain information about the thermodynamical stability or the mechanism of layer displacements, e.g. Co (Edwards & Lipson, 1942[link]; Frey & Boysen, 1981[link]), SiC (Jagodzinski, 1972[link]; Pandey et al., 1980a[link],b[link],c[link]), ZnS (Müller, 1952[link]; Mardix & Steinberger, 1970[link]; Frey et al., 1986[link]) and others.

Certain laws may be derived for the reduced integrated intensities of diffuse reflections. `Reduction' in this context means a division of the diffuse scattering along l by the structure factor, or the difference structure factor if [\langle F \rangle \neq 0]. This procedure is valuable if the number of stacking faults rather than the complete solution of the diffraction problem is required.

The discussion given above has been made under the assumption that the full symmetry of the layers is maintained in the statistics. Obviously, this is not necessarily true if external lower symmetries influence the disorder. An important example is the generation of stacking faults during plastic deformation. Problems of this kind need a complete reconsideration of symmetries. Furthermore, it should be pointed out that a treatment with the aid of an extended Ising model as described above is irrelevant in most cases. Simplified procedures describing the diffuse scattering of intrinsic, extrinsic, twin stacking faults and others have been described in the literature. Since their influence on structure determination can generally be neglected, the reader is referred to the literature for additional information.

4.2.4.3. Two-dimensional disorder of chains

| top | pdf |

In this section disorder phenomena are considered which are related to chain-like structural elements in crystals. This topic includes the so-called `1D crystals' where translational symmetry (in direct space) exists in one direction only – crystals in which highly anisotropic binding forces are responsible for chain-like atomic groups, e.g. compounds which exhibit a well ordered 3D framework structure with tunnels in a unique direction in which atoms, ions or molecules are embedded. Examples are compounds with platinum, iodine or mercury chains, urea inclusion compounds with columnar structures (organic or inorganic), 1D ionic conductors, polymers etc. Diffuse-scattering studies of 1D conductors have been carried out in connection with investigations of stability/instability problems, incommensurate structures, phase transitions, dynamic precursor effects etc. These questions are not treated here. For general reading of diffuse scattering in connection with these topics see, e.g., Comes & Shirane (1979[link], and references therein). Also excluded are specific problems related to polymers or liquid crystals (mesophases) (see Chapter 4.4[link] ) and magnetic structures with chain-like spin arrangements.

Trivial diffuse scattering occurs as 1D Bragg scattering (diffuse layers) by internally ordered chains. Diffuse phenomena in reciprocal space are due to `longitudinal' disordering within the chains (along the unique direction) as well as to `transverse' correlations between different chains over a restricted volume. Only static aspects are considered; diffuse scattering resulting from collective excitations or diffusion-like phenomena which are of inelastic or quasielastic origin are not treated here.

4.2.4.3.1. Scattering by randomly distributed collinear chains

| top | pdf |

As found in any elementary textbook of diffraction the simplest result of scattering by a chain with period c [l({\bf r}) = l(z) = \textstyle\sum\limits_{n_{3}} \delta (z - n_{3}c) \eqno(4.2.4.12)] is described by one of the Laue equations: [G(L) = |L(L)|^{2} = \sin^{2} \pi NL/\sin^{2} \pi L \eqno(4.2.4.13)] which gives broadened profiles for small N. In the context of phase transitions the Ornstein–Zernike correlation function is frequently used, i.e. (4.2.4.13[link]) is replaced by a Lorentzian: [1 / \{\xi^{2} + 4 \pi^{2}(L - l)^{2}\}, \eqno(4.2.4.14)] where ξ denotes the correlation length.

In the limiting case [N \rightarrow \infty], (4.2.4.13[link]) becomes [\textstyle\sum\limits_{l} \delta (L - l). \eqno(4.2.4.15)] The scattering by a real chain a(r) consisting of molecules with structure factor [F_{M}] is therefore determined by [F_{M}({\bf H}) = \textstyle\sum\limits_{j} f_{j} \exp \{2 \pi i(Hx_{j} + Ky_{j} + Lz_{j})\}. \eqno(4.2.4.16)] The Patterson function is: [\eqalignno{P({\bf r}) &= (1 / c) \textstyle\int \textstyle\int |F_{0}(H, K)|^{2} \cos 2 \pi (Hx + Ky) \;\hbox{d}H \;\hbox{d}K &\cr &\quad + (2 / c) \textstyle\sum\limits_{l} \textstyle\int \textstyle\int |F_{l}|^{2} \exp \{2 \pi i(Hx + Ky)\} &\cr &\quad \times \exp \{- 2 \pi ilz\} \;\hbox{d}H \;\hbox{d}K, &(4.2.4.17)}] where the index l denotes the only relevant position [L = l] (the subscript M is omitted).

The intensity is concentrated in diffuse layers perpendicular to [{\bf c}^{*}] from which the structural information may be extracted. Projections are: [\textstyle\int a({\bf r}) \;\hbox{d}z = \textstyle\int \textstyle\int F_{0}(H, K) \exp \{2 \pi i(Hx + Ky)\} \;\hbox{d}H \;\hbox{d}K \eqno(4.2.4.18)] [\textstyle\int \textstyle\int a({\bf r}) \;\hbox{d}x \;\hbox{d}y = (2 / c) \textstyle\sum\limits_{l} F_{l}(00l) \exp (- 2 \pi ilz). \eqno(4.2.4.19)] Obviously the z parameters can be determined by scanning along a meridian (00L) through the diffuse sheets (diffractometer recording). Owing to intersection of the Ewald sphere with the set of planes the meridian cannot be recorded on one photograph; successive equi-inclination photographs are necessary. Only in the case of large c spacings is the meridian well approximated in one photograph.

There are many examples where a tendency to cylindrical symmetry exists: chains with p-fold rotational or screw symmetry around the preferred direction or assemblies of chains (or domains) with statistical orientational distribution around the texture axis. In this context it should be mentioned that symmetry operations with rotational parts belonging to the 1D rod groups actually occur, i.e. not only p = 2, 3, 4, 6.

In all these cases a treatment in the frame of cylindrical coordinates is advantageous (see, e.g., Vainshtein, 1966[link]): [\let\normalbaselines\relax\openup3pt\matrix{\hbox{Direct space}\hfill &\hbox{Reciprocal space}\hfill\cr x = r \cos \psi \hfill&  H = H_{r} \cos \Psi\hfill\cr y = r \sin \psi \hfill &  K = H_{r} \sin \Psi\hfill\cr z = z\hfill& L = L\hfill}] [\eqalignno{a(r, \psi, z) &= \textstyle\int \textstyle\int \textstyle\int F({\bf H}) \exp \{- 2 \pi i[H_{r}r \cos (\psi - \Psi) + Lz]\}\cr &\quad \times H_{r} \;\hbox{d}H_{r} \;\hbox{d}\Psi \;\hbox{d}L &(4.2.4.20)}] [\eqalignno{F({\bf H}) &= \textstyle\int \textstyle\int \textstyle\int a(r, \psi, z) \exp \{2 \pi i[H_{r}r \cos (\psi - \Psi) + Lz]\}\cr &\quad \times r \;\hbox{d}r \;\hbox{d}\psi \;\hbox{d}z. &(4.2.4.21)}] The integrals may be evaluated by the use of Bessel functions: [J_{n}(u) = {\textstyle{1\over 2}} \pi i^{n} \textstyle\int \exp \{i(u \cos \varphi + n\varphi)\} \;\hbox{d}\varphi] [(u = 2 \pi rH_{r}\hbox{; } \varphi = \psi - \Psi)].

The 2D problem [a = a(r, \psi)] is treated first; an extension to the general case [a(r, \psi, z)] is easily made afterwards.

Along the theory of Fourier series one has: [\eqalign{a(r, \psi) &= \textstyle\sum\limits_{n} a_{n}(r) \exp \{in \psi\}\cr a_{n}(r) &= {1\over 2 \pi} \int a(r, \psi) \exp \{- in \psi\} \;\hbox{d}\psi\cr} \eqno(4.2.4.22)] or with: [\eqalign{\alpha_{n} &= {1\over 2 \pi} \int a(r, \psi) \cos (n\psi) \;\hbox{d}\psi\cr \beta_{n} &= {1\over 2 \pi} \int a(r, \psi) \sin (n\psi) \;\hbox{d}\psi\cr a_{n} (r) &= |a_{n} (r)| \exp \{- i\psi_{n} (r)\} \cr |a_{n} (r)| &= \sqrt{a_{n}^{2} + \beta_{n}^{2}} \cr \psi_{n} (r) &= \hbox{arctan } \beta_{n}/\alpha_{n}.\cr}] If contributions to anomalous scattering are neglected a(r, ψ) is a real function: [a (r, \psi) = \textstyle\sum\limits_{n} |a_{n} (r)| \cos [n\psi - \psi_{n} (r)]. \eqno(4.2.4.23)] Analogously, one has [F (H_{r}, \Psi) = \textstyle\sum\limits_{n} |F_{n} (H_{r})| \exp (in \Psi). \eqno(4.2.4.24)] [F (H_{r}, \Psi)] is a complex function; [F_{n} (H_{r})] are the Fourier coefficients which are to be evaluated from the [a_{n} (r)]: [\eqalignno{F_{n} (H_{r}) &= {1\over 2\pi} \int F (H_{r}, \Psi) \exp \{- in\Psi\} \;\hbox{d} \Psi &\cr \noalign{\vskip5pt} &= \exp \{in \pi/2\} \int a_{n} (r) J_{n} (2\pi r H_{r}) 2\pi r \;\hbox{d}r &\cr F (H_{r}, \Psi) &= \textstyle\sum\limits_{n} \exp \{in [\Psi + (\pi/2)]\} \textstyle\int a_{n} (r) &\cr &\quad \times J_{n} (2\pi r H_{r}) 2\pi r \;\hbox{d}r &(4.2.4.25)}] [\eqalignno{a (r, \psi) &= \textstyle\sum\limits_{n} \exp \{in [\Psi - (\pi/2)]\} \textstyle\int F_{n} (H_{r}) &\cr &\quad \times J_{n} (2\pi r H_{r}) 2\pi H_{r} \;\hbox{d}H_{r}. &(4.2.4.26)}] The formulae may be used for calculation of diffuse intensity distribution within a diffuse sheet, in particular when the chain molecule is projected along the unique axis [cf. equation (4.2.4.18)[link]].

Special cases are:

  • (a) Complete cylinder symmetry [\eqalignno{F (H_{r}) &= 2\pi \textstyle\int a (r) J_{0} (2\pi r H_{r}) r \;\hbox{d}r &(4.2.4.27)\cr a (r) &= 2\pi \textstyle\int F (H_{r}) J_{n} (2\pi r H_{r}) H_{r} \;\hbox{d}H_{r}. &(4.2.4.28)}%(4.2.4.28)]

  • (b) p-fold symmetry of the projected molecule [a (r, \psi) = {a [r, \psi + (2\pi/p)]}] [\eqalignno{F_{p} (H_{r}, \Psi) &= \textstyle\sum\limits_{n} \exp \{inp[\Psi + (\pi/2)]\} &\cr \noalign{\vskip5pt} &\quad \times \textstyle\int a_{np} (r) J_{np} (2\pi r H_{r}) 2\pi r \;\hbox{d}r &(4.2.4.29)\cr \noalign{\vskip5pt} a_{p} (r, \psi) &= \textstyle\sum\limits_{n} |a_{np} (r)| \cos [np\psi - \psi_{np} (r)]. &(4.2.4.30)}%(4.2.4.30)] Only Bessel functions [J_{0}, J_{p}, J_{2p}, \ldots] occur. In most cases [J_{2p}] and higher orders may be neglected.

  • (c) Vertical mirror planes

    Only cosine terms occur, i.e. all [\beta_{n} = 0] or [\psi_{n} (r) = 0].

The general 3D expressions valid for extended chains with period c [equation (4.2.4.12)[link]] are found in an analogous way: [a (r, \psi, z) = a_{M} (r, \psi, z) * l(z)] [\eqalignno{F ({\bf H}) &= F_{l} (H_{r}, \Psi, L) = F_{M} ({\bf H}) L (L) &\cr &= \textstyle\int \textstyle\int \textstyle\int a_{M} (r, \psi, z) \exp \{2\pi i [H_{r} r \cos (\psi - \Psi) + Lz]\} &\cr &\quad \times 2\pi r \;\hbox{d}r\;\hbox{d}\psi\;\hbox{d}z &(4.2.4.31)}] using a series expansion analogous to (4.2.4.23[link]) and (4.2.4.24[link]): [a_{nl} (r) = {1\over 2\pi} \int \int a_{M} \exp \{- i (n\psi - 2\pi lz)\} \;\hbox{d}\psi \;\hbox{d}z \eqno(4.2.4.32)] [F_{nl} (H_{r}) = \exp \{in \pi/2\} \textstyle\int a_{nl} (r) J_{n} (2\pi H_{r} r) 2\pi r \;\hbox{d}r \eqno(4.2.4.33)] one has: [{F_{l} ({\bf H}) = \textstyle\sum\limits_{n} \exp \{in [\Psi + (\pi/2)]\} \textstyle\int a_{nl} (r) J_{n} (2\pi H_{r} r) 2\pi r \;\hbox{d}r.} \eqno(4.2.4.34)] In practice the integrals are often replaced by discrete summation of j atoms at positions: [r = r_{j}], [\psi = \psi_{j}], [z = z_{j}] [(0 \leq z_{j} \lt  c)]: [\eqalignno{F_{l} ({\bf H}) &= \textstyle\sum\limits_{j} \textstyle\sum\limits_{n} f_{j} J_{n} (2\pi H_{r} r_{j}) \exp \{-in \psi_{j}\} &\cr &\quad \times \exp (2\pi ilz_{j}) \exp \{in [\Psi + (\pi/2)]\} &(4.2.4.35)}] or [F_{l} ({\bf H}) = \textstyle\sum\limits_{n} (\alpha_{n} + i\beta_{n}) \exp \{in \Psi\}] [\eqalign{\alpha_{n} &= \textstyle\sum\limits_{j} f_{j} J_{n} (2\pi H_{r} r_{j}) \cos \{n [(\pi/2) - \psi_{j}] + 2\pi lz_{j}\} \cr \beta_{n} &= \textstyle\sum\limits_{j} f_{j} J_{n} (2\pi H_{r} r_{j}) \sin \{n [(\pi/2) - \psi_{j}] + 2\pi lz_{j}\}.}] Intensity in the lth diffuse layer is given by [\eqalignno{I_{l} &= \textstyle\sum\limits_{n} \textstyle\sum\limits_{n'} [(\alpha_{n} \alpha_{n'} + \beta_{n} \beta_{n'}) + i (\alpha_{n'} \beta_{n} - \alpha_{n} \beta_{n'})] &\cr &\quad \times \exp \{i (n - n') \Psi\}. &(4.2.4.36)}]

  • (a) Cylinder symmetry (free rotating molecules around the chain axis or statistical averaging with respect to ψ over an assembly of chains). Only component [F_{0l}] occurs: [F_{0l} (H_{r}, L) = 2\pi \textstyle\int \textstyle\int \langle a_{M} \rangle J_{0} (2\pi H_{r} r) \exp \{2\pi ilz\} r \;\hbox{d}r \;\hbox{d}z] or [F_{0l} (H_{r}, L) = \textstyle\sum\limits_{j} f_{j} J_{0} (2\pi H_{r} r_{j}) \exp \{2\pi ilz_{j}\}.] In particular, [F_{00} (H_{r})] determines the radial component of the molecule projected along z: [F_{00} (H_{r}) = \textstyle\sum\limits_{j} f_{j} J_{0} (2\pi H_{r} r_{j}).]

  • (b) p-fold symmetry of a plane molecule (or projected molecule) as outlined previously: only components np instead of n occur. Bessel functions [J_{0}] and [J_{p}] are sufficient in most cases.

  • (c) Vertical mirror plane: see above.

  • (d) Horizontal mirror plane (perpendicular to the chain): Exponentials [\exp \{2\pi ilz\}] in equation (4.2.4.32)[link] may be replaced by [\cos 2\pi lz].

  • (e) Twofold symmetry axis perpendicular to the chain axis (at positions [\psi = 0], [2\pi/p, \ldots]). Exponentials in equation (4.2.4.32),[link] [\exp \{- i (np\psi - 2\pi lz)\}], are replaced by the corresponding cosine term [\cos (np\psi + 2\pi lz)].

Formulae concerning the reverse method (Fourier synthesis) are not given here (see, e.g., Vainshtein, 1966[link]). Usually there is no practical use in diffuse-scattering work because it is very difficult to separate out a single component [F_{nl}]. Every diffuse layer is affected by all components [F_{nl}]. There is a chance if one diffuse layer corresponds predominantly to one Bessel function.

4.2.4.3.2. Disorder within randomly distributed collinear chains

| top | pdf |

Deviations from strict periodicities in the z direction within one chain may be due to loss of translational symmetry of the centres of the molecules along z and/or due to varying orientations of the molecules with respect to different axes, such as azimuthal misorientation, tilting with respect to the z axis or combinations of both types. As in 3D crystals, there may or may not exist 1D structures in an averaged sense.

4.2.4.3.2.1. General treatment

| top | pdf |

All formulae given in this section are only special cases of a 3D treatment (see, e.g., Guinier, 1963[link]). The 1D lattice (4.2.4.12)[link] is replaced by a distribution: [\eqalignno{d (z) &= \textstyle\sum\limits_{\nu} \delta (z - z_{\nu}) \cr \noalign{\vskip5pt}D (L) &= \textstyle\sum\limits_{\nu} \exp \{2\pi iLz_{\nu}\} &(4.2.4.37)\cr \noalign{\vskip5pt}F ({\bf H}) &= F_{M} ({\bf H}) D (L).}] The Patterson function is given by [P ({\bf r}) = [a_{M} ({\bf r}) * a_{M} - ({\bf r})] * [d (z) * d - (z)]. \eqno(4.2.4.38)] Because the autocorrelation function [w = d * d] is centrosymmetric [{w (z) = N\delta (z) + \textstyle\sum\limits_{\nu} \textstyle\sum\limits_{\mu} \delta [z - (z_{\nu} - z_{\mu})]}\ {+ \textstyle\sum\limits_{\nu} \textstyle\sum\limits_{\mu} \delta [z + (z_{\nu} - z_{\mu})],} \eqno(4.2.4.39)] the interference function [W (L)] [(= |D (L)|^{2})] is given by [\eqalignno{W (L) &= N + 2 \textstyle\sum\limits_{\nu} \textstyle\sum\limits_{\mu} \cos 2\pi [L (z_{\nu} - z_{\mu})] &(4.2.4.40)\cr \noalign{\vskip5pt}I ({\bf H}) &= |F_{M} ({\bf H})|^{2} W (L). &(4.2.4.41)}%(4.2.4.41)] Sometimes, e.g. in the following example of orientational disorder, there is an order only within domains. As shown in Section 4.2.3[link], this may be treated by a box or shape function [b(z) = 1] for [z \leq z_{N}] and 0 elsewhere. [\eqalignno{d (z) &= d_{\infty} b (z) \cr \noalign{\vskip5pt}a ({\bf r}) &= a_{M} ({\bf r}) * [d_{\infty} b (z)] &(4.2.4.42)\cr \noalign{\vskip5pt}F ({\bf H}) &= F_{M} ({\bf H}) [D_{\infty} * B (L)]}] with [\eqalign{&\ \ b (z) * b - (z) \leftrightarrow |B (L)|^{2}\cr &I = |F_{M} ({\bf H})|^{2}| D_{\infty} * B (L)|^{2}.} \eqno(4.2.4.43)] If the order is perfect within one domain one has [D_{\infty} (L) \simeq \sum \delta (L - l)]; [(D_{\infty} * B) = \sum D (L - l)]; i.e. each reflection is affected by the shape function.

4.2.4.3.2.2. Orientational disorder

| top | pdf |

A misorientation of the chain molecules with respect to one another is taken into account by different structure factors [F_{M}]. [I ({\bf H}) = \textstyle\sum\limits_{\nu} \textstyle\sum\limits_{\mu} F_{\nu} ({\bf H}) F_{\mu} ({\bf H})^{+} \exp \{2\pi iL (z_{\nu} - z_{\mu})\}. \eqno(4.2.4.44)] A further discussion follows the same arguments outlined in Section 4.2.3[link]. For example, a very simple result is found in the case of uncorrelated orientations. Averaging over all pairs [F_{\nu} F_{\mu}^{+}] yields [I ({\bf H}) = N (\langle |F|^{2}\rangle - \langle |F|\rangle^{2}) + \langle |F|\rangle^{2} L (L), \eqno(4.2.4.44a)] where [\eqalign{\langle |F|\rangle^{2} &= 1/N^{2} \langle F_{\nu} F_{\mu}^{+}\rangle \cr \noalign{\vskip5pt}&= \textstyle\sum\limits_{\nu} \alpha_{\nu} F_{\nu} ({\bf H}) \textstyle\sum\limits_{\mu} \alpha_{\mu} F_{\mu}^{+} ({\bf H}) \quad (\nu \neq \mu) \cr \noalign{\vskip5pt}\langle |F|^{2}\rangle &= 1/N \langle F_{\nu} F_{\nu}^{+}\rangle = \textstyle\sum\limits_{\nu} \alpha_{\nu} |F_{\nu} ({\bf H})|^{2}.}] Besides the diffuse layer system there is a diffuse background modulated by the H dependence of [[\langle |F ({\bf H})|^{2}\rangle - \langle |F ({\bf H})|\rangle^{2}]].

4.2.4.3.2.3. Longitudinal disorder

| top | pdf |

In this context the structure factor of a chain molecule is neglected. Irregular distances between the molecules within a chain occur owing to the shape of the molecules, intrachain interactions and/or interaction forces via a surrounding matrix. A general discussion is given by Guinier (1963)[link]. It is convenient to reformulate the discrete Patterson function, i.e. the correlation function (4.2.4.39[link]). [\eqalignno{w (z) &= N\delta (z) + \textstyle\sum\limits_{\nu} \delta [z \pm (z_{\nu} - z_{\nu + 1})] &\cr &\quad + \textstyle\sum\limits_{\nu} \delta [z \pm (z_{\nu} - z_{\nu + 2})] + \ldots &(4.2.4.39a)}] in terms of continuous functions [a_{\mu} (z)] which describe the probability of finding the µth neighbour within an arbitrary distance [\eqalignno{w' (z) &= w (z)/N = \delta (z) + a_{1} (z) + a_{-1} (z) + \ldots \cr &\quad + a_{\mu} (z) a_{-\mu} (z) + \ldots &(4.2.4.45)}] [[\textstyle\int a_{\mu} (z) \;\hbox{d}z = 1, a_{\mu} (z) = a_{-\mu} (-z)].]

There are two principal ways to define [a_{\mu} (z)]. The first is the case of a well defined one-dimensional lattice with positional fluctuations of the molecules around the lattice points, i.e. long-range order is retained: [a_{\mu} (z) = \mu c_{0} + z_{\mu}], where [z_{\mu}] denotes the displacement of the µth molecule in the chain. Frequently used are Gaussian distributions: [c' \exp \{- (z - \mu c_{0})^{2}/2\Delta^{2}\}] ([c'] = normalizing constant; [\Delta ] = standard deviation). Fourier transformation [equation (4.2.4.45)[link]] gives the well known result [I_{d} \sim (1 - \exp \{- L^{2} \Delta^{2}\}),] i.e. a monotonically increasing intensity with L (modulation due to a molecular structure factor neglected). This result is quite analogous to the treatment of the scattering of independently vibrating atoms. If (short-range) correlations exist between the molecules the Gaussian distribution is replaced by a multivariate normal distribution where correlation coefficients [\kappa^{\mu}] [(0  \lt  \kappa  \lt  1)] between a molecule and its µth neighbour are incorporated. [\kappa^{\mu}] is defined by the second moment: [\langle z_{0} z_{\mu}\rangle/\Delta^{2}]. [a_{\mu} (z) = c'' \exp \{- (z - \mu c_{0})^{2}/2\Delta^{2} (1 - \kappa^{\mu})\}.] Obviously the variance increases if the correlation diminishes and reaches an upper bound of twice the single site variance. Fourier transformation gives an expression for diffuse intensity (Welberry, 1985)[link]: [\eqalignno{I_{d} (L) &\sim \exp \{- L^{2} \Delta^{2}\} \textstyle\sum\limits_{j} (- L^{2} \Delta^{2})\hskip 2pt^{j}/j! &\cr &\quad \times (1 - \kappa^{2j})/(1 + \kappa^{2j} - 2\kappa\hskip 2pt^{j} \cos 2\pi Lc_{0}). &(4.2.4.46)}] For small Δ, terms with [j \gt 1] are mostly neglected. The terms become increasingly important with higher values of L. On the other hand, [\kappa\hskip 2pt^{j}] becomes smaller with increasing j, each additional term in equation (4.2.4.46)[link] becomes broader and, as a consequence, the diffuse planes in reciprocal space become broader with higher L.

In a different way – in the paracrystal method – the position of the second and subsequent molecules with respect to some reference zero point depends on the actual position of the predecessor. The variance of the position of the µth molecule relative to the first becomes unlimited. There is a continuous transition to a fluid-like behaviour of the chain molecules. This 1D paracrystal (sometimes called distortions of second kind) is only a special case of the 3D paracrystal concept (see Hosemann & Bagchi, 1962[link]; Wilke, 1983[link]). Despite some difficulties with this concept (Brämer, 1975[link]; Brämer & Ruland, 1976[link]) it is widely used as a theoretical model for describing diffraction of highly distorted lattices. One essential development is to limit the size of a paracrystalline grain so that fluctuations never become too large (Hosemann, 1975[link]).

If this concept is used for the 1D case, [a_{\mu}(z)] is defined by convolution products of [a_{1}(z)]. For example, the probability of finding the next-nearest molecule is given by [a_{2}(z) = \textstyle\int a_{1}(z')a_{1}(z - z') \;\hbox{d}z' = a_{1}(z) * a_{1}(z)] and, generally: [a_{\mu}(z) = a_{1}(z) * a_{1}(z) * a_{1}(z) * \ldots * a_{1}(z)] (µ-fold convolution).

The mean distance between next-nearest neighbours is [\langle c \rangle = \textstyle\int z' a_{1}(z') \;\hbox{d}z'] and between neighbours of the µth order: [\mu \langle c \rangle]. The average value of [a_{\mu} = 1/\langle c \rangle], which is also the value of w(z) for [z > z_{k}], where the distribution function is completely smeared out. The general expression for the interference function G(L) is [G(L) = 1 + \textstyle\sum\limits_{\mu} \{F^{\mu} + F^{+\mu}\} = \hbox{Re} \{(1 + F)/(1 - F)\} \eqno(4.2.4.47)] with [F(L) \leftrightarrow a_{1}(z), F^{\mu}(L) \leftrightarrow a_{\mu}(z)].

With [F = |F|e^{i\chi}\ (\chi = L\langle c \rangle)], equation (4.2.4.47)[link] is written: [G(L) = [1 - |F(L)|^{2}]/[1 - 2|F(L)| \cos \chi + |F(L)|^{2}]. \eqno(4.2.4.47a)] [Note the close similarity to the diffuse part of equation (4.2.4.5),[link] which is valid for 1D disorder problems.]

This function has maxima of height [(1 + |F|)/(1 - |F|)] and minima of height [(1 - |F|)/(1 + |F|)] at positions [lc^{*}] and [\left(l + {1\over 2}\right)c^{*}], respectively. With decreasing [|F|] the oscillations vanish; a critical L value (corresponding to [z_{k}]) may be defined by Gmax/Gmin ≲ 1.2. Actual values depend strongly on [F({\bf H})].

The paracrystal method is substantiated by the choice [a_{1}(z)], i.e. the disorder model. Again, frequently used is a Gaussian distribution: [\eqalign{a_{1}(z) &= 1/\sqrt{2\pi} \Delta \cdot \exp \{- (z - \langle c \rangle)^{2}/2\Delta^{2}\}\cr a_{\mu}(z) &= 1/\sqrt{\mu} \cdot 1/\sqrt{2\pi} \Delta \cdot \exp \{- (z - \mu \langle c \rangle)^{2}/2\mu \Delta^{2}\}} \eqno(4.2.4.48)] with the two parameters [\langle c \rangle, \Delta].

There are peaks of height [1/[\pi^{2} L^{2} (\Delta/\langle c \rangle)^{2}]] which obviously decrease with [L^{2}] and [(\Delta/\langle c \rangle)^{2}]. The oscillations vanish for [|F| \simeq 0.1], i.e. [1/\langle c \rangle \simeq 0.25/\Delta]. The width of the mth peak is [\Delta_{m} = \sqrt{m} \Delta]. The integral reflectivity is approximately [1/\langle c \rangle [1 - \pi^{2} L^{2} (\Delta/\langle c \rangle)^{2}]] and the integral width (defined by integral reflectivity divided by peak reflectivity) (background subtracted!) [1/\langle c \rangle \pi^{2} L^{2} (\Delta/\langle c \rangle)^{2}] which, therefore, increases with [L^{2}]. In principle the same results are given by Zernike & Prins (1927)[link]. In practice a single Gaussian distribution is not fully adequate and modified functions must be used (Rosshirt et al., 1985[link]).

A final remark concerns the normalization [equation (4.2.4.39)].[link] Going from (4.2.4.39)[link] to (4.2.4.45)[link] it is assumed that N is a large number so that the correct normalization factors [(N - |\mu|)] for each [a_{\mu}(z)] may be approximated by a uniform N. If this is not true then [\eqalignno{G(L) &= N + \textstyle\sum\limits_{\mu} (N - |\mu|)(F^{\mu} + F^{+\mu}) \cr &= N \;\hbox{Re}\; \{(1 + |F|)/(1 - |F|)\} &(4.2.4.49)\cr &\quad -2 \;\hbox{Re}\; \{|F|(1 - |F|^{N})/(1 - |F|^{2})\}.}] The correction term may be important in the case of relatively small (1D) domains. As mentioned above, the structure factor of a chain molecule was neglected. The H dependence of [F_{m}], of course, obscures the intensity variation of the diffuse layers as described by (4.2.4.47a)[link].

The matrix method developed for the case of planar disorder was adapted to 2D disorder by Scaringe & Ibers (1979)[link]. Other models and corresponding expressions for diffuse scattering are developed from specific microscopic models (potentials), e.g. in the case of Hg3−δAsF6 (Emery & Axe, 1978[link]; Radons et al., 1983[link]), hollandites (Beyeler et al., 1980[link]; Ishii, 1983[link]), iodine chain compounds (Endres et al., 1982[link]) or urea inclusion compounds (Forst et al., 1987[link]).

4.2.4.3.3. Correlations between different almost collinear chains

| top | pdf |

In real cases there are more or less strong correlations between different chains at least within small domains. Deviations from a strict (3D) order of chain-like structural elements are due to several reasons: shape and structure of the chains, varying binding forces, thermodynamical or kinetic considerations.

Many types of disorder occur. (1) Relative shifts parallel to the common axis while projections along this axis give a perfect 2D ordered net (`axial disorder'). (2) Relative fluctuations of the distances between the chains (perpendicular to the unique axis) with short-range order along the transverse a and/or b directions. The net of projected chains down to the ab plane is distorted (`net distortions'). Disorder of types (1) and (2) is sometimes correlated owing to non-uniform cross sections of the chains. (3) Turns, twists and torsions of chains or parts of chains. This azimuthal type of disorder may be treated similarly to the case of azimuthal disorder of single-chain molecules. Correlations between axial shifts and torsions produce `screw shifts' (helical structures). Torsion of chain parts may be of dynamic origin (rotational vibrations). (4) Tilting or bending of the chains in a uniform or non-uniform way (`conforming/non-conforming'). Many of these types and a variety of combinations between them are found in polymer and liquid crystals and are treated therefore separately. Only some simple basic ideas are discussed here in brief.

For the sake of simplicity the paracrystal concept in combination with Gaussians is used again. Distribution functions are given by convolution products of next-nearest-neighbour distribution functions. As long as averaged lattice directions and lattice constants in a plane perpendicular to the chain axis exist, only two functions [a_{100} = a_{1}(xyz)] and [a_{010} = a_{2}(xyz)] are needed to describe the arrangement of next-nearest chains. Longitudinal disorder is treated as before by a third distribution function [a_{001} = a_{3}(xyz)]. The phenomena of chain bending or tilting may be incorporated by an x and y dependence of [a_{3}]. Any general fluctuation in the spatial arrangement of chains is given by [a_{mpq} = a_{1} * \ldots * a_{1} * a_{2} * \ldots * a_{2} * a_{3} * \ldots * a_{3}. \eqno(4.2.4.50)] (m-fold, p-fold, q-fold self-convolution of [a_{1}, a_{2}, a_{3}], respectively.) [w({\bf r}) = \delta ({\bf r}) \textstyle\sum\limits_{m} \textstyle\sum\limits_{p} \textstyle\sum\limits_{q} [a_{mpq} ({\bf r}) + a_{-mpq} ({\bf r})]. \eqno(4.2.4.51)] [a_{\nu}\ (\nu = 1, 2, 3)] are called fundamental functions. If an averaged lattice cannot be defined, more fundamental functions [a_{\nu}] are needed to account for correlations between them.

By Fourier transformations the interference function is given by [\eqalign{G({\bf H}) &= \textstyle\sum\limits_{m} \textstyle\sum\limits_{p} \textstyle\sum\limits_{q} F_{1}^{m} F_{2}^{p} F_{3}^{q} = G_{1} G_{2} G_{3}\hbox{;}\cr G_{\nu} &= \hbox{Re} \{(1 + |F_{\nu}|)/(1 - |F_{\nu}|)\}.} \eqno(4.2.4.52)] If Gaussian functions are assumed, simple pictures are derived. For example: [\eqalignno{a_{1}({\bf r} + \langle {\bf a} \rangle) &= 1/(2\pi)^{3/2} \cdot 1/(\Delta_{11} \Delta_{12} \Delta_{13}) &\cr &\quad \times \exp \left\{- {\textstyle{1\over 2}}[(x^{2}/\Delta_{11}^{2}) + (y^{2}/\Delta_{12}^{2}) + (z^{2}/\Delta_{13}^{2})]\right\} &&\cr&&(4.2.4.53)}] describes the distribution of neighbours in the x direction (mean distance [\langle a \rangle]). Parameter [\Delta_{13}] concerns axial, [\Delta_{11}] and [\Delta_{12}] radial and tangential fluctuations, respectively. Pure axial distribution along c is given by projection of [a_{1}] on the z axis, pure net distortions by projection on the [x - y] plane. If the chain-like structure is neglected the interference function [G_{1}({\bf H}) = \exp \{-2 \pi^{2}(\Delta_{11}^{2} H^{2} + \Delta_{12}^{2} K^{2} + \Delta_{13}^{2} L^{2})\} \eqno(4.2.4.54)] describes a set of diffuse planes perpendicular to [{\bf a}^{*}] with mean distance [1/\langle a \rangle]. These diffuse layers broaden along H with [m\Delta_{11}] and decrease in intensity along K and L monotonically. There is an ellipsoidal-shaped region in reciprocal space defined by main axes of length [1/\Delta_{11}, 1/\Delta_{12}, 1/\Delta_{13}] with a limiting surface given by [|F| \simeq 0.1], beyond which the diffuse intensity is completely smeared out. The influence of [a_{2}] may be discussed in an analogous way.

If the chain-like arrangement parallel to c [equation (4.2.4.12)][link] is taken into consideration, [l(z) = \textstyle\sum\limits_{n_{3}} \delta (z - n_{3} c)\hbox{;}] the set of planes perpendicular to [{\bf a}^{*}] (and/or [{\bf b}^{*}]) is subdivided in the L direction by a set of planes located at [l \cdot 1/c] [equation (4.2.4.15)].[link]

Longitudinal disorder is given by [a_{3}(z)] [equation (4.2.4.48),[link] [\Delta_{33} = \Delta]] and leads to two intersecting sets of broadened diffuse layer systems.

Particular cases like pure axial distributions [(\Delta_{11}, \Delta_{12} \sim 0)], pure tangential distributions (net distortions: [\Delta_{11}, \Delta_{13} \sim 0]), uniform bending of chains or combinations of these effects are discussed in the monograph of Vainshtein (1966)[link].

4.2.4.4. Disorder with three-dimensional correlations (defects, local ordering and clustering)

| top | pdf |

4.2.4.4.1. General formulation (elastic diffuse scattering)

| top | pdf |

In this section general formulae for diffuse scattering will be derived which may best be applied to crystals with a well ordered average structure, characterized by (almost) sharp Bragg peaks. Textbooks and review articles concerning defects and local ordering are by Krivoglaz (1969)[link], Dederichs (1973)[link], Peisl (1975)[link], Schwartz & Cohen (1977)[link], Schmatz (1973[link], 1983[link]), Bauer (1979)[link], and Kitaigorodsky (1984)[link]. A series of interesting papers on local order is given by Young (1975)[link] and also by Cowley et al. (1979)[link]. Expressions for polycrystalline sample material are given by Warren (1969)[link] and Fender (1973)[link].

Two general methods may be applied:

  • (a) the average difference cluster method, where a representative cluster of scattering differences between the average structure and the cluster is used; and

  • (b) the method of short-range-order correlation functions where formal parameters are introduced.

Both methods are equivalent in principle. The cluster method is generally more convenient in cases where a single average cluster is a good approximation. This holds for small concentrations of clusters with sufficient space in between. The method of short-range-order parameters is optimal in cases where isolated clusters are not realized and the correlations do not extend to long distances. Otherwise periodic solutions are more convenient in most cases.

In any case, the first step towards the solution of the diffraction problem is the accurate determination of the average structure. As described in Section 4.2.3.2[link] important information on fractional occupations, interstitials and displacements (unusual thermal parameters) of atoms may be derived. Unfortunately all defects contribute to diffuse scattering; hence one has to start with the assumption that the disorder to be interpreted is predominant. Fractional occupancy of certain lattice sites by two or more kinds of atoms plays an important role in the literature, especially in metallic or ionic structures. Since vacancies may be treated as atoms with zero scattering amplitude, structures containing vacancies may be formally treated as multi-component systems.

Since the solution of the diffraction problem should not be restricted to metallic systems with a simple (primitive) structure, we have to consider the structure of the unit cell – as given by the average structure – and the propagation of order according to the translation group separately. In simple metallic systems this difference is immaterial. It is well known that the thermodynamic problem of propagation of order in a three-dimensional crystal can hardly be solved analytically in a general way. Some solutions have been published with the aid of the so-called Ising model using next-nearest-neighbour interactions. They are excellent for an understanding of the principles of order–disorder phenomena, but they can scarcely be applied quantitatively in practical problems. Hence, methods have been developed to derive the propagation of order from the diffraction pattern by means of Fourier transformation. This method has been described qualitatively in Section 4.2.3.1[link], and will be used here for a quantitative application. In a first approximation the assumption of a small number of different configurations of the unit cell is made, represented by the corresponding number of structure factors. Displacements of atoms caused by the configurations of the neighbouring cells are excluded. This problem will be treated subsequently.

The finite number of structures of the unit cell in the disordered crystal is given by [F_{\nu} ({\bf r}) = \textstyle\sum\limits_{j} \textstyle\sum\limits_{\mu} \pi_{j\mu}^{\nu}\; f_{\mu} ({\bf r - r}_{j}). \eqno(4.2.4.55)] Note that [F_{\nu}({\bf r})] is defined in real space, and [{\bf r}_{j}] gives the position vector of site j; [\pi_{j\mu}^{\nu} = 1] if in the νth structure factor the site j is occupied by an atom of kind µ, and 0 elsewhere.

In order to apply the laws of Fourier transformation adequately, it is useful to introduce the distribution function of [F_{\nu}] [\pi_{\nu} ({\bf r}) = \textstyle\sum\limits_{{\bf n}} \pi_{{\bf n}\nu} \delta ({\bf r} - {\bf n}) \eqno(4.2.4.56)] with [\pi_{{\bf n}\nu} = 1], if the cell [{\bf n} = n_{1}{\bf a} + n_{2}{\bf b} + n_{3}{\bf c}] has the [F_{\nu}] structure, and [\pi_{{\bf n}\nu} = 0] elsewhere.

In the definitions given above [\pi_{{\bf n}\nu}] are numbers (scalars) assigned to the cell. Since all these are occupied we have [\textstyle\sum\limits_{\nu} \pi_{\nu} ({\bf r}) = l({\bf r})] with [l({\bf r})] = lattice in real space.

The structure of the disordered crystal is given by [\textstyle\sum\limits_{\nu} \pi_{\nu} ({\bf r}) * F_{\nu} ({\bf r}). \eqno(4.2.4.57)] [\pi_{\nu} ({\bf r})] consists of [\alpha_{\nu} N] points, where [N = N_{1} N_{2} N_{3}] is the total (large) number of unit cells and [\alpha_{\nu}] denotes the a priori probability (concentration) of the νth cell occupation.

It is now useful to introduce [\Delta \pi_{\nu} ({\bf r}) = \pi_{\nu} ({\bf r}) - \alpha_{\nu} l({\bf r}) \eqno(4.2.4.58)] with [\textstyle\sum\limits_{\nu} \Delta \pi_{\nu} ({\bf r}) = \textstyle\sum\limits_{\nu} \pi_{\nu} ({\bf r}) - l({\bf r}) \textstyle\sum\limits_{\nu} \alpha_{\nu} = l({\bf r}) - l({\bf r}) = 0.] Introducing (4.2.4.58)[link] into (4.2.4.57)[link]: [\eqalignno{\textstyle\sum\limits_{\nu} \pi_{\nu} ({\bf r}) * F_{\nu} ({\bf r}) &= \textstyle\sum\limits_{\nu} \Delta \pi_{\nu} ({\bf r}) * F_{\nu} ({\bf r}) + l({\bf r}) \textstyle\sum\limits_{\nu} \alpha_{\nu} F_{\nu} ({\bf r}) &\cr \noalign{\vskip5pt}&= \textstyle\sum\limits_{\nu} \Delta \pi_{\nu} ({\bf r}) * F_{\nu} ({\bf r}) + l({\bf r}) \langle F({\bf r})\rangle. &(4.2.4.59)}] Similarly: [\eqalign{&\phantom{\sum_{\nu} \alpha_{\nu} \Delta F_{\nu}} \Delta F_{\nu} ({\bf r}) = F_{\nu} ({\bf r}) - \langle F({\bf r})\rangle\cr &\textstyle\sum\limits_{\nu} \alpha_{\nu} \Delta F_{\nu} ({\bf r}) = \textstyle\sum\limits_{\nu} \alpha_{\nu} F_{\nu} ({\bf r}) - \textstyle\sum\limits_{\nu} \alpha_{\nu} \langle F({\bf r})\rangle = 0.} \eqno(4.2.4.60)] Using (4.2.4.60)[link] it follows from (4.2.4.58)[link] that [\eqalignno{\textstyle\sum\limits_{\nu} \pi_{\nu} ({\bf r}) * F_{\nu} ({\bf r}) &= \textstyle\sum\limits_{\nu} \Delta \pi_{\nu} ({\bf r}) * \Delta F_{\nu} ({\bf r}) &\cr \noalign{\vskip5pt}&\quad + \textstyle\sum\limits_{\nu} \Delta \pi_{\nu} ({\bf r}) * \langle F({\bf r})\rangle + l({\bf r}) * \langle F({\bf r})\rangle &\cr \noalign{\vskip5pt}&= \textstyle\sum\limits_{\nu} \Delta \pi_{\nu} ({\bf r}) * \Delta F_{\nu} ({\bf r}) + l({\bf r}) * \langle F({\bf r})\rangle. &\cr&&(4.2.4.61)}] Comparison with (4.2.4.59)[link] yields [\textstyle\sum\limits_{\nu} \Delta \pi_{\nu} ({\bf r}) * F_{\nu} ({\bf r}) = \textstyle\sum\limits_{\nu} \Delta \pi_{\nu} ({\bf r}) * \Delta F_{\nu} ({\bf r}).] Fourier transformation of (4.2.4.61)[link] gives [\textstyle\sum\limits_{\nu} \Pi_{\nu} ({\bf H}) F_{\nu} ({\bf H}) = \textstyle\sum\limits_{\nu} \Delta \Pi_{\nu} ({\bf H}) \Delta F_{\nu} ({\bf H}) + L({\bf H}) \langle F({\bf H})\rangle] with [\textstyle\sum\limits_{\nu} \Delta \Pi_{\nu} ({\bf H}) = 0\hbox{;} \quad \textstyle\sum\limits_{\nu} \Delta F_{\nu} ({\bf H}) = 0.] The expression for the scattered intensity is therefore [\eqalignno{I({\bf H}) &= \left|\textstyle\sum\limits_{\nu} \Delta \Pi_{\nu} ({\bf H}) \Delta F_{\nu} ({\bf H})\right|^{2} + |L({\bf H}) \langle F({\bf H})\rangle |^{2} &\cr \noalign{\vskip5pt}&\quad + L({\bf H}) \left\{\langle F^{+} ({\bf H})\rangle \textstyle\sum\limits_{\nu} \Delta \Pi_{\nu} ({\bf H}) F_{\nu} ({\bf H})\right. &\cr \noalign{\vskip5pt}&\quad + \left.\langle F^{+} ({\bf H})\rangle \textstyle\sum\limits_{\nu} \Delta \Pi_{\nu}^{+} ({\bf H}) F_{\nu}^{+} ({\bf H})\right\}. &(4.2.4.62)}] Because of the multiplication by L(H) the third term in (4.2.4.62)[link] contributes to sharp reflections only. Since they are correctly given by the second term in (4.2.4.62)[link], the third term vanishes, Hence, the diffuse part is given by [I_{d} ({\bf H}) = \left|\textstyle\sum\limits_{\nu} \Delta \Pi_{\nu} ({\bf H}) \Delta F_{\nu} ({\bf H})\right|^{2}. \eqno(4.2.4.63)] For a better understanding of the behaviour of diffuse scattering it is useful to return to real space: [\eqalignno{i_{d} ({\bf r}) &= \textstyle\sum\limits_{\nu} \Delta \pi_{\nu} ({\bf r}) * \Delta F_{\nu} ({\bf r}) * \textstyle\sum\limits_{\nu'} \Delta \pi_{\nu'} (-{\bf r}) * \Delta F_{\nu'} (-{\bf r}) &\cr \noalign{\vskip5pt}&= \textstyle\sum\limits_{\nu} \textstyle\sum\limits_{\nu'} \Delta \pi_{\nu}   ({\bf r}) * \Delta \pi_{\nu'} (-{\bf r}) * \Delta F_{\nu} ({\bf r}) * \Delta F_{\nu'} (-{\bf r}) &\cr&&(4.2.4.64)}] and with (4.2.4.58)[link]: [\eqalignno{i_{d} ({\bf r}) &= \textstyle\sum\limits_{\nu} \textstyle\sum\limits_{\nu'} [\pi_{\nu} ({\bf r}) - \alpha_{\nu} l({\bf r})] * [\pi_{\nu'} (-{\bf r}) - \alpha_{\nu'} l({\bf r})] &\cr \noalign{\vskip5pt}&\quad * \Delta F_{\nu} ({\bf r}) * \Delta F_{\nu'} (-{\bf r}). &(4.2.4.65)}] Evaluation of this equation for a single term yields [\eqalignno{&[\pi_{\nu} ({\bf r}) * \pi_{\nu'} (-{\bf r}) - \alpha_{\nu} l({\bf r}) * \pi_{\nu'} (-{\bf r}) &\cr \noalign{\vskip5pt}&\quad - \alpha_{\nu'} l({\bf r}) * \pi_{\nu} ({\bf r}) + \alpha_{\nu} \alpha_{\nu'} l({\bf r}) * l({\bf r})] &\cr \noalign{\vskip5pt}&\quad * \Delta F_{\nu} ({\bf r}) * \Delta F_{\nu'} ({\bf r}). &(4.2.4.66)}] Since l(r) is a periodic function of points, all convolution products with l(r) are also periodic. For the final evaluation the decrease of a number of overlapping points (maximum N) in the convolution products with increasing displacements of the functions is neglected (no particle-size effect). Then (4.2.4.66)[link] becomes [\eqalignno{&[\pi_{\nu} ({\bf r}) * \pi_{\nu'} (-{\bf r}) - N\alpha_{\nu} \alpha_{\nu'} l({\bf r}) - N\alpha_{\nu} \alpha_{\nu'} l({\bf r}) + N\alpha_{\nu} \alpha_{\nu'} l({\bf r})] &\cr \noalign{\vskip5pt}&\quad * \Delta F_{\nu} ({\bf r}) * \Delta F_{\nu'} (-{\bf r}) &\cr \noalign{\vskip5pt}&\quad\quad = [\pi_{\nu} ({\bf r}) * \pi_{\nu'} (-{\bf r}) - N\alpha_{\nu} \alpha_{\nu'} l({\bf r})] &\cr \noalign{\vskip5pt}&\quad\quad\quad * \Delta F_{\nu} ({\bf r}) * \Delta F_{\nu'} (-{\bf r}). &(4.2.4.67)}] If the first term in (4.2.4.67)[link] is considered, the convolution of the two functions for a given distance n counts the number of coincidences of the function [\pi_{\nu} ({\bf r})] with [\pi_{\nu'} (-{\bf r})]. This quantity is given by [Nl({\bf r}) \alpha_{\nu} p_{\nu \nu'} (-{\bf r})], where [\alpha_{\nu} p_{\nu \nu'} ({\bf r})] is the probability of a pair occupation in the r direction.

Equation (4.2.4.67)[link] then reads: [\eqalignno{&[Nl({\bf r}) \alpha_{\nu} p_{\nu \nu'} (-{\bf r}) - N\alpha_{\nu} \alpha_{\nu'} l({\bf r})] &\cr \noalign{\vskip5pt}&\quad * \Delta F_{\nu} ({\bf r}) * \Delta F_{\nu'} (-{\bf r}) &\cr \noalign{\vskip5pt}&\quad\quad = N\alpha_{\nu} [p'_{\nu \nu'} (-{\bf r}) l({\bf r})] &\cr \noalign{\vskip5pt}&\quad\quad\quad * [\Delta F_{\nu} ({\bf r}) * \Delta F_{\nu'} (-{\bf r})] &(4.2.4.68)}] with [p'_{\nu \nu'} ({\bf r}) = p_{\nu \nu'} ({\bf r}) - \alpha_{\nu'}]. The function [\alpha_{\nu} p'_{\nu \nu'} ({\bf r})] is usually called the pair-correlation function [g\leftrightarrow \alpha_{\nu \nu'|{\bf n n}'}] in the physical literature.

The following relations hold: [\displaylines{\hfill \textstyle\sum\limits_{\nu} \alpha_{\nu} = 1 \hfill (4.2.4.69a)\cr \noalign{\vskip5pt}\hfill \textstyle\sum\limits_{\nu'} \alpha_{\nu} p_{\nu \nu'} ({\bf r}) = \alpha_{\nu'}\hbox{;}\quad \textstyle\sum\limits_{\nu'} \alpha_{\nu} p'_{\nu \nu'} ({\bf r}) = 0 \hfill (4.2.4.69b)\cr \noalign{\vskip5pt}\hfill \textstyle\sum\limits_{\nu} \alpha_{\nu} p_{\nu \nu'} ({\bf r}) = \alpha_{\nu}\hbox{;}\quad \textstyle\sum\limits_{\nu} \alpha_{\nu} p'_{\nu \nu'} ({\bf r}) = 0 \hfill (4.2.4.69c)\cr \noalign{\vskip5pt}\hfill\alpha_{\nu} p_{\nu \nu'} ({\bf r}) = \alpha_{\nu'} p_{\nu' \nu} (-{\bf r}). \hfill (4.2.4.69d)}] Also, functions normalized to unity are in use. Obviously the following relation is valid: [p'_{\nu \nu'} (0) = \delta_{\nu \nu'} - \alpha_{\nu'}].

Hence: [\alpha_{\nu \nu'|{\bf n n}'} = \alpha_{\nu'} p'_{\nu \nu'} ({\bf r}) / (\delta_{\nu \nu'} - \alpha_{\nu'})] is unity for [{\bf r} = 0\; ({\bf n} = {\bf n}')]. This property is especially convenient in binary systems.

With (4.2.4.68)[link], equation (4.2.4.64)[link] becomes [{i_{d} ({\bf r}) = N \textstyle\sum\limits_{\nu} \textstyle\sum\limits_{\nu'} [\alpha_{\nu} p'_{\nu \nu'} (-{\bf r}) l({\bf r})] * [\Delta F_{\nu} ({\bf r}) * \Delta F_{\nu'} (-{\bf r})]} \eqno(4.2.4.70)] and Fourier transformation yields [{I_{d} ({\bf H}) = N \textstyle\sum\limits_{\nu} \textstyle\sum\limits_{\nu'} [\alpha_{\nu} P'_{\nu \nu'} (-{\bf H}) * L({\bf H})] \Delta F_{\nu} ({\bf H}) \Delta F_{\nu'}^{+} ({\bf H}).} \eqno(4.2.4.71)] It may be concluded from equations (4.2.4.69)[link] that all functions [p'_{\nu \nu'} ({\bf r})] may be expressed by [p'_{11} ({\bf r})] in the case of two structure factors [F_{1}], [F_{2}]. Then all [p'_{\nu \nu'} ({\bf r})] are symmetric in r; the same is true for the [P'_{\nu \nu'} ({\bf H})]. Consequently, the diffuse reflections described by (4.2.4.71)[link] are all symmetric. The position of the diffuse peak depends strongly on the behaviour of [p'_{\nu \nu'} ({\bf r})]; in the case of cluster formation Bragg peaks and diffuse peaks coincide. Diffuse superstructure reflections are observed if the [p'_{\nu \nu'} ({\bf r})] show some damped periodicities.

It should be emphasized that the condition [p'_{\nu \nu'} ({\bf r}) = p'_{\nu \nu'} (-{\bf r})] may be violated for [\nu \neq \nu'] if more than two cell occupations are involved. As shown below, the possibly asymmetric functions may be split into symmetric and antisymmetric parts. From equation (4.2.3.8)[link] it follows that the Fourier transform of the antisymmetric part of [p'_{\nu \nu'} ({\bf r})] is also antisymmetric. Hence, the convolution in the two terms in square brackets in (4.2.4.71)[link] yields an antisymmetric contribution to each diffuse peak, generated by the convolution with the reciprocal lattice L(h).

Obviously, equation (4.2.4.71)[link] may also be applied to primitive lattices, occupied by two or more kinds of atoms. Then the structure factors [F_{\nu}] are merely replaced by the atomic scattering factors [f_{\nu}] and the [\alpha_{\nu}] are equivalent to the concentrations of atoms [c_{\nu}]. In terms of the [\alpha_{\nu \nu'|{\bf n n}'}] (Warren short-range-order parameters) equation (4.2.4.71)[link] reads [{I_{d} ({\bf H}) = N(\;\overline{f^{2}} - \bar{f}^{2}) \textstyle\sum\limits_{\bf n} \textstyle\sum\limits_{{\bf n}'} \alpha_{\nu \nu'|{\bf n n}'} \exp \{2 \pi i {\bf H} ({\bf n} - {{\bf n}'})\}.} \eqno(4.2.4.71a)] In the simplest case of a binary system A, B [\displaylines{\alpha_{{\bf nn}'} = (1 - p_{AB|{\bf n n}'}) / c_{B} = (1 - p_{BA|{\bf n n'}}) / c_{A}\hbox{;} \cr \noalign{\vskip5pt} c_{A} p_{AB} = c_{B} p_{BA}\hbox{;} \quad p_{AA} = 1 - p_{AB}\hbox{;}}] [\eqalignno{I_{d} ({\bf H}) &= Nc_{A} c_{B} (\;f_{A} - f_{B})^{2} \textstyle\sum\limits_{\bf n} \textstyle\sum\limits_{{\bf n}'} \alpha_{{\bf nn}'} \cr \noalign{\vskip5pt} &\quad \times \exp \{2\pi i {\bf H} ({\bf n - n}')\}. &(4.2.4.71b)}] [The exponential in (4.2.4.71b)[link] may even be replaced by a cosine term owing to the centrosymmetry of this particular case.]

It should be mentioned that the formulations of the problem in terms of pair probabilities, pair correlation functions, short-range-order parameters or concentration waves (Krivoglaz, 1969[link]) are equivalent. Using continuous electron (or nuclear) density functions where site occupancies are implied, the Patterson function may be used, too (Cowley, 1981[link]).

4.2.4.4.2. Random distribution

| top | pdf |

As shown above in the case of random distributions all [p'_{\nu \nu'} ({\bf r})] are zero, except for [{\bf r} = 0]. Consequently, [p'_{\nu \nu'} ({\bf r}) l({\bf r})] may be replaced by [\alpha_{\nu} p'_{\nu \nu'} ({\bf r}) = \alpha_{\nu} \delta_{\nu \nu'} - \alpha_{\nu} \alpha_{\nu'}. \eqno(4.2.4.72)] According to (4.2.4.59)[link] and (4.2.4.61)[link] the diffuse scattering can be given by the Fourier transformation of [\eqalign{&\textstyle\sum\limits_{\nu} \textstyle\sum\limits_{\nu'} \Delta \pi_{\nu} ({\bf r}) * \Delta \pi_{\nu'} (-{\bf r}) * F_{\nu} ({\bf r}) * F_{\nu} (-{\bf r})\cr \noalign{\vskip5pt}&\quad = \textstyle\sum\limits_{\nu} \textstyle\sum\limits_{\nu'} p'_{\nu \nu'} ({\bf r}) * F_{\nu} ({\bf r}) * F_{\nu'} (-{\bf r})}] or with (4.2.4.72)[link]: [i_{d} ({\bf r}) = N \textstyle\sum\limits_{\nu} \textstyle\sum\limits_{\nu'} [\alpha_{\nu} \delta_{\nu \nu'} - \alpha_{\nu} \alpha_{\nu'}] * F_{\nu} ({\bf r}) * F_{\nu'} (-{\bf r}).] Fourier transformation gives [\eqalignno{I_{d} ({\bf H}) &= N \left\{\textstyle\sum\limits_{\nu} \alpha_{\nu} |F_{\nu} ({\bf H})|^{2} - \textstyle\sum\limits_{\nu} \alpha_{\nu} F_{\nu} ({\bf H}) \textstyle\sum\limits_{\nu'} \alpha_{\nu'} F_{\nu'}^{+} ({\bf H})\right\} &\cr \noalign{\vskip5pt}&= N \{\langle |F({\bf H})|^{2}\rangle - |\langle F({\bf H})\rangle |^{2}\}. &(4.2.4.73)}] This is the most general form of any diffuse scattering of systems ordered randomly (`Laue scattering'). Occasionally it is called `incoherent scattering' (see Section 4.2.2[link]).

4.2.4.4.3. Short-range order in multi-component systems

| top | pdf |

The diffuse scattering of a disordered binary system without displacements of the atoms has already been discussed in Section 4.2.4.4.1.[link] It could be shown that all distribution functions [p'_{\nu \nu'} ({\bf r})] are mutually dependent and may be replaced by a single function [cf. (4.2.4.69)[link]]. In that case [p'_{\nu \nu'} ({\bf r}) = p'_{\nu \nu'} (-{\bf r})] was valid for all. This condition, however, may be violated in multi-component systems. If a tendency towards an [F_{1} F_{2} F_{3}] order in a ternary system is assumed, for example, [p_{12} ({\bf r})] is apparently different from [p_{12} (-{\bf r})]. In this particular case it is useful to introduce [\eqalign{\langle p'_{\nu \nu'} ({\bf r})\rangle &= {\textstyle{1\over 2}} [p'_{\nu \nu'} ({\bf r}) + p'_{\nu \nu'} (-{\bf r})]\hbox{;}\cr \Delta p'_{\nu \nu'} ({\bf r}) &= {\textstyle{1\over 2}} [p'_{\nu \nu'} ({\bf r}) - p'_{\nu \nu'} (-{\bf r})]}] and their Fourier transforms [\langle P'_{\nu \nu'} ({\bf H})\rangle, \Delta P'_{\nu \nu'} ({\bf H})], respectively.

The asymmetric correlation functions are therefore expressed by [\eqalign{p'_{\nu \nu'} ({\bf r}) &= \langle p'_{\nu \nu'} ({\bf r})\rangle + \Delta p'_{\nu \nu'} ({\bf r})\hbox{;}\cr \noalign{\vskip5pt}p'_{\nu \nu'} (-{\bf r}) &= \langle p'_{\nu \nu'} ({\bf r})\rangle - \Delta p'_{\nu \nu'} ({\bf r})\hbox{;}\cr \noalign{\vskip5pt}\Delta p'_{\nu \nu} ({\bf r}) &= 0.}] Consequently, [i_{d} ({\bf r})] (4.2.4.70)[link] and [I_{d} ({\bf H})] (4.2.4.71)[link] may be separated according to the symmetric and antisymmetric contributions. The final result is: [\eqalignno{I_{d} ({\bf H}) &= N \Bigl\{\textstyle\sum\limits_{\nu} \alpha_{\nu} [P'_{\nu \nu} ({\bf H}) * L({\bf H})] |F_{\nu} ({\bf H})|^{2} &\cr \noalign{\vskip5pt}&\quad + \textstyle\sum\limits_{\nu > \nu'} \alpha_{\nu} [\langle P'_{\nu \nu'} ({\bf H})\rangle * L({\bf H})] &\cr \noalign{\vskip5pt}&\quad \times [\Delta F_{\nu} ({\bf H}) \Delta F_{\nu'}^{+} ({\bf H}) + \Delta F_{\nu}^{+} ({\bf H}) \Delta F_{\nu'} ({\bf H})] &\cr \noalign{\vskip5pt}&\quad + \textstyle\sum\limits_{\nu > \nu'} \alpha_{\nu} [\Delta P'_{\nu \nu'} (-{\bf H}) * L({\bf H})] &\cr \noalign{\vskip5pt}&\quad \times [\Delta F_{\nu} ({\bf H}) \Delta F_{\nu'}^{+} ({\bf H}) - \Delta F_{\nu}^{+} ({\bf H}) \Delta F_{\nu'} ({\bf H})]\Bigr\}. &(4.2.4.74)}] Obviously, antisymmetric contributions to line profiles will only occur if structure factors of acentric cell occupations are involved. This important property may be used to draw conclusions with respect to structure factors involved in the statistics. It should be mentioned here that the Fourier transform of the antisymmetric function [\Delta p'_{\nu \nu'} ({\bf r})] is imaginary and antisymmetric. Since the last term in (4.2.4.74)[link] is also imaginary, the product of the two factors in brackets is real, as it should be.

4.2.4.4.4. Displacements: general remarks

| top | pdf |

Even small displacements may have an important influence on the problem of propagation of order. Therefore, no structural treatments other than the introduction of formal parameters (e.g. Landau's theory) have been published in the literature. Most of the examples with really reliable results refer to binary systems, and even these represent very crude approximations, as will be shown below. For this reason we shall restrict ourselves here to binary systems, although general formulae where displacements are included may be developed in a formal way.

Two kinds of atoms, [f_{1}({\bf r})] and [f_{2}({\bf r})], are considered. Obviously, the position of any given atom is determined by its surroundings. Their extension depends on the forces acting on the atom under consideration. These may be very weak in the case of metals (repulsive forces, so-called `size effect'), but long-range effects have to be expected in ionic crystals. For the development of formulae authors have assumed that small displacements [\Delta_{\nu \nu'}({\bf r})] may be assigned to the pair correlation functions [p'_{\nu \nu'}({\bf r})] by adding a phase factor [\exp \{2 \pi i {\bf H} \cdot \Delta_{\nu \nu'} {\bf r}\}] which is then expanded in the usual way: [{\exp \{2 \pi i {\bf H} \cdot \Delta_{\nu \nu'} {\bf r}\} \simeq 1 + 2 \pi i {\bf H} \cdot \Delta_{\nu \nu'} {\bf r} - 2 [\pi {\bf H} \cdot \Delta_{\nu \nu'} {\bf r}]^{2}.} \eqno(4.2.4.75)] The displacements and correlation probabilities are separable if the change of atomic scattering factors in the angular range considered may be neglected. The formulae in use are given in the next section. As shown below, this method represents nothing other than a kind of average over certain sets of displacements. For this purpose the correct solution of the problem has to be discussed. In the simplest model the displacements are due to next-nearest neighbours only. It is assumed further that the configurations rather than the displacements determine the position of the central atom and a general displacement of the centre of the first shell does not occur (no influence of a strain field). Obviously, the formal correlation function of pairs is not independent of displacements. This difficulty may be avoided either by assuming that the pair correlation function has already been separated from the diffraction data, or by theoretical calculations of the correlation function (mean-field method) (Moss, 1966[link]; de Fontaine, 1972[link], 1973[link]). The validity of this procedure is subject to the condition that the displacements have no influence on the correlation functions themselves.

The observation of a periodic average structure justifies the definition of a periodic array of origins which normally depends on the degree of order. Local deviations of origins may be due to fluctuations in the degree of order and due to the surrounding atoms of a given site occupation. For example, a b.c.c. lattice with eight nearest neighbours is considered. It is assumed that only these have an influence on the position of the central atom owing to different forces of the various configurations. With two kinds of atoms, there are [2^{9} = 512] possible configurations of the cluster (central atom plus 8 neighbours). Symmetry considerations reduce this number to 28. Each is characterized by a displacement vector. Hence, their a priori probabilities and the propagation of 28 different configurations have to be determined. Since each atom has to be considered as the centre once, this problem may be treated by introducing 28 different atomic scattering factors as determined from the displacements: [f_{\nu} ({\bf r}) \exp \{2 \pi i {\bf H} \cdot \Delta {\bf r}_{\nu \nu'}\}]. The diffraction problem has to be solved with the aid of the propagation of order of overlapping clusters. This is demonstrated by a two-dimensional model with four nearest neighbours (Fig. 4.2.4.1)[link]. Here the central and the neighbouring cluster (full and broken lines) overlap with two sites in the, e.g., x direction. Hence, only neighbouring clusters with the same overlapping pairs are admitted. These restrictions introduce severe difficulties into the problem of propagation of cluster ordering which determines the displacement field. Since it was assumed that the problem of pair correlation had been solved, the cluster probabilities may be derived by calculating [\alpha_{\nu} l ({\bf r}) \textstyle\prod\limits_{{\bf n} \neq 0} p'_{\nu \nu'} ({\bf r - n}). \eqno(4.2.4.76)] In the product only next-nearest neighbours have to be included. This must be performed for the central cluster [({\bf r} = 0)] and for the reference cluster at [{\bf r} = n_{1} {\bf a} + n_{2} {\bf b} + n_{3} {\bf c}], because all are characterized by different displacements. So far, possible displacement of the centre has not been considered; this may also be influenced by the problem of propagation of cluster ordering. These displacement factors should best be attached to the function describing the propagation of order which determines, in principle, the local fluctuations of the lattice constants (strain field etc.). This may be understood by considering a binary system with a high degree of order but with atoms of different size. Large fluctuations of lattice constants are involved in the case of exsolution of the two components because of their different lattice parameters, but they become small in the case of superstructure formation where a description in terms of antiphase domains is reasonable (equal lattice constants). This example demonstrates the mutual dependence of ordering and displacements which is mostly neglected in the literature.

[Figure 4.2.4.1]

Figure 4.2.4.1 | top | pdf |

Construction of the correlation function in the method of overlapping clusters.

The method of assigning phase factors to the pair correlation function is now discussed. Pair correlation functions average over all pairs of clusters having the same central atom. An analogous argument holds for displacements: using pair correlations for the determination of displacements means nothing other than averaging over all displacements caused by various clusters around the same central atom. There remains the general strain field due to the propagation of order, whereas actual displacements of atoms are realized by fluctuations of configurations. Since large fluctuations of this type occur in highly disordered crystals, the displacements become increasingly irrelevant. Hence, the formal addition of displacement factors to the pair correlation function does not yield too much information about the structural basis of the displacements. This situation corresponds exactly to the relationship between a Patterson function and a real structure: the structure has to be found which explains the more or less complicated function completely, and its unique solution is rather difficult. These statements seem to be necessary because in most publications related to this subject these considerations are not taken into account adequately. Displacements usually give rise to antisymmetric contributions to diffuse reflections. As pointed out above, the influence of displacements has to be considered as phase factors which may be attached either to the structure factors or to the Fourier transforms [P'_{\nu \nu'} ({\bf H})] of the correlation functions in equation (4.2.4.71)[link]. As has been mentioned in the context of equation (4.2.4.74)[link] antisymmetric contributions will occur if acentric structure factors are involved. Apparently, this condition is met by the phase factors of displacements. In consequence, antisymmetric contributions to diffuse reflections may also originate from the displacements. This fact can also be demonstrated if the assignment of phase factors to the Fourier transforms of the correlation functions is advantageous. In this case equations (4.2.4.69a,b)[link] are no longer valid because the functions [p'_{\nu \nu'} ({\bf r})] become complex. The most important change is the relation corresponding to (4.2.4.69)[link]: [\alpha_{\nu} p'_{\nu \nu'} ({\bf r}) = \alpha_{\nu'} p_{\nu \nu'}^{'+} (-{\bf r}) \leftrightarrow \alpha_{\nu} P'_{\nu \nu'} ({\bf H}) = \alpha_{\nu'} P_{\nu'}^{'+} ({\bf H}). \eqno(4.2.4.77)] Strictly speaking we have to replace the a priori probabilities [\alpha_{\nu}] by complex numbers [\alpha_{\nu} \exp (2 \pi i \Delta {\bf r}_{\nu} \cdot {\bf H})] which are determined by the position of the central atom. In this way all correlations between displacements may be included with the aid of the clusters mentioned above. To a rough approximation it may be assumed that no correlations of this kind exist. In this case the complex factors may be assigned to the structure factors involved. Averaging over all displacements results in diffraction effects which are very similar to a static Debye–Waller factor for all structure factors. On the other hand, the thermal motion of atoms is treated similarly. Obviously both factors affect the sharp Bragg peaks. Hence, this factor can easily be determined by the average structure which contains a Debye–Waller factor including static and thermal displacements. It should be pointed out, however, that these static displacements cause elastic diffuse scattering which cannot be separated by inelastic neutron scattering techniques.

A careful study of the real and imaginary parts of [\langle p'_{\nu \nu'}({\bf r}) \rangle = \langle p'_{\nu \nu'}({\bf r}) \rangle_{R} + \langle p'_{\nu \nu'}({\bf r}) \rangle_{I}] and [\Delta p'_{\nu \nu'}({\bf r}) = \Delta p'_{\nu \nu'}({\bf r})_{R} + \Delta p'_{\nu \nu'}({\bf r})_{I}] and their Fourier transforms results, after some calculations, in the following relation for diffuse scattering: [\eqalignno{I_{d} &\simeq N \textstyle\sum\limits_{\nu} \alpha_{\nu} |\Delta F_{\nu}({\bf H})|^{2} \{[\langle P'_{\nu \nu}({\bf H}) \rangle - \Delta P'_{\nu \nu}({\bf H})] * L({\bf H})\} &\cr &\quad+ 2 N \textstyle\sum\limits_{\nu > \nu'} \alpha_{\nu} (\Delta F_{\nu} \Delta F_{\nu'}^{+})_{R} &\cr &\quad\times \{[\langle P'_{\nu \nu'}({\bf H}) \rangle_{R} - \Delta P'_{\nu \nu'}({\bf H})_{I}] * L({\bf H})\} &\cr &\quad+ 2 N \textstyle\sum\limits_{\nu > \nu'} \alpha_{\nu} (\Delta F_{\nu} \Delta F_{\nu'}^{+})_{I} &\cr &\quad\times \{[\langle P'_{\nu \nu'}({\bf H}) \rangle_{I} - \Delta P'_{\nu \nu'}({\bf H})_{R}] * L({\bf H})\}. &(4.2.4.78)}] It should be noted that all contributions are real. This follows from the properties of Fourier transforms of symmetric and antisymmetric functions. All [\Delta P_{\nu \nu'}({\bf H})] are antisymmetric; hence they generate antisymmetric contributions to the line profiles. In contrast to equation (4.2.4.75)[link], the real and the imaginary parts of the structure factors contribute to the asymmetry of the line profiles.

4.2.4.4.5. Distortions in binary systems

| top | pdf |

In substitutional binary systems (primitive cell with only one sublattice) the Borie–Sparks method is widely used (Sparks & Borie, 1966[link]; Borie & Sparks, 1971[link]). The method is formulated in the short-range-order-parameter formalism. The diffuse scattering may be separated into two parts (a) owing to short-range order and (b) owing to static displacements.

Corresponding to the expansion (4.2.4.75)[link], [I_{d} = I_{\rm sro} + I_{2} + I_{3}], where [I_{\rm sro}] is given by equation (4.2.4.71b)[link] and the correction terms [I_{2}] and [I_{3}] relate to the linear and the quadratic term in (4.2.4.75)[link]. The intensity expression will be split into terms of A–A, A–B, … pairs. More explicitly [\Delta_{\nu \nu'}{\bf r} = {\bf u}_{{\bf n} \nu} - {\bf u}_{{\bf n}' \nu'}] and with the following abbreviations: [\eqalignno{\delta_{{\bf nn}'|AA} &= {\bf u}_{{\bf n}|A} - {\bf u}_{{\bf n}'|A} = x_{{\bf nn}'|AA} {\bf a} + y_{{\bf nn}'|AA} {\bf b} + z_{{\bf nn}'|AA} {\bf c}\cr \delta_{{\bf nn}'|AB} &= {\bf u}_{{\bf n}|A} - {\bf u}_{{\bf n}'|B} = \ldots \cr F_{{\bf nn}'|AA} &= f_{A}^{2} / (\;f_{A} - f_{B})^{2} \cdot [(c_{A} / c_{B}) + \alpha_{{\bf nn}'}]\cr F_{{\bf nn}'|BB} &= f_{B}^{2} / (\;f_{A} - f_{B})^{2} \cdot [(c_{B} / c_{A}) + \alpha_{{\bf nn}'}]\cr F_{{\bf nn}'|AB} &= 2f_{A} f_{B} / (\;f_{A} - f_{B})^{2} \cdot (1 - \alpha_{{\bf nn}'}) = F_{{\bf nn}'|BA}}] one finds (where the short-hand notation is self-explanatory): [\eqalignno{I_{2} &= 2 \pi ic_{A} c_{B} (\;f_{A} - f_{B})^{2} \textstyle\sum\limits_{\bf n} \textstyle\sum\limits_{{\bf n}'} \{H \cdot [F_{{\bf nn}'|AA} \langle x_{{\bf nn}'|AA} \rangle &\cr &\quad + F_{{\bf nn}'|BB} \langle x_{{\bf nn}'|BB} \rangle + F_{{\bf nn}'|AB} \langle x_{{\bf nn}'|AB} \rangle] + K \cdot [\hbox{`}y\hbox{'}] &\cr &\quad + L \cdot [\hbox{`}z\hbox{'}]\} \exp \{2 \pi i{\bf H} \cdot ({\bf n} - {\bf n}')\} &(4.2.4.79)}] [\eqalignno{I_{3} &= c_{A} c_{B} (\;f_{A} - f_{B})^{2} (-2 \pi)^{2} \textstyle\sum\limits_{\bf n} \textstyle\sum\limits_{{\bf n}'} \{H^{2} [F_{{\bf nn}'} \langle x_{{\bf nn}'|AA}^{2} \rangle &\cr &\quad + F_{{\bf nn}'|BB} \langle x_{{\bf nn}'|BB}^{2} \rangle + F_{{\bf nn}'|AB} \langle x_{{\bf nn}'|AB}^{2} \rangle] &\cr &\quad + K^{2} \cdot [\hbox{`}y^{2}\hbox{'}] + L^{2} \cdot [\hbox{`}z^{2}\hbox{'}] &\cr &\quad + HK [F_{{\bf nn}'|AA} \langle (xy)_{{\bf nn}'|AA} \rangle + F_{{\bf nn}'|BB} \langle (xy)_{{\bf nn}'|BB} \rangle &\cr &\quad + F_{{\bf nn}'|AB} \langle (xy)_{{\bf nn}'|AB} \rangle] &\cr &\quad + KL [\hbox{`}(yz)\hbox{'}] + LH [\hbox{`}(zx)\hbox{'}]\} &\cr &\quad \times \exp \{2 \pi i{\bf H} \cdot ({\bf n} - {\bf n}')\}. &(4.2.4.80)}] With further abbreviations [\eqalign{\gamma_{{\bf nn}'|x} &= 2 \pi (F_{{\bf nn}'|AA} \langle x_{{\bf nn}'|AA} \rangle + F_{{\bf nn}'|BB} \langle x_{{\bf nn}'|BB} \rangle \cr &\quad + F_{{\bf nn}'|AB} \langle x_{{\bf nn}'|AB} \rangle) \cr \gamma_{{\bf nn}'|y} &= \ldots \cr \gamma_{{\bf nn}'|z} &= \ldots \cr \delta_{{\bf nn}'|x} &= (-2 \pi^{2}) (F_{{\bf nn}'|AA} \langle x_{{\bf nn}'|AA}^{2} \rangle + F_{{\bf nn}'|BB} \langle x_{{\bf nn}'|BB}^{2} \rangle \cr &\quad + F_{{\bf nn}'|AB} \langle x_{{\bf nn}'|AB}^{2} \rangle) \cr \delta_{{\bf nn}'|y} &= \ldots \cr \delta_{{\bf nn}'|z} &= \ldots \cr \varepsilon_{{\bf nn}'|xy} &= (-4 \pi^{2}) (F_{{\bf nn}'|AA} \langle (xy)_{{\bf nn}'|AA} \rangle + F_{{\bf nn}'|BB} \langle (xy)_{{\bf nn}'|BB} \rangle \cr &\quad + F_{{\bf nn}'|AB} \langle (xy)_{{\bf nn}'|AB} \rangle) \cr \varepsilon_{{\bf nn}'|yz} &= \ldots \cr \varepsilon_{{\bf nn}'|zx} &= \ldots \cr {I}_{2} &= c_{A} c_{B} (\;f_{A} - f_{B})^{2} \textstyle\sum\limits_{\bf n} \textstyle\sum\limits_{{\bf n}'} i (\gamma_{{\bf nn}'|x} + \gamma_{{\bf nn}'|y} + \gamma_{{\bf nn}'|z}) \cr &\quad \times \exp \{2 \pi i{\bf H} \cdot ({\bf n} - {\bf n}')\}\cr I_{3} &= c_{A} c_{B} (\;f_{A} - f_{B})^{2} \textstyle\sum\limits_{\bf n} \textstyle\sum\limits_{{\bf n}'} (\delta_{{\bf nn}'|x} H^{2} + \delta_{{\bf nn}'|y} K^{2}\cr &\quad + \delta_{{\bf nn}'|z} L^{2} + \varepsilon_{{\bf nn}'|xy} HK + \varepsilon_{{\bf nn}'|yz} KL\cr &\quad + \varepsilon_{{\bf nn}'|zx} LH \exp \{2 \pi i {\bf H} \cdot ({\bf n} - {\bf   n}')\}.}] If the [F_{{\bf nn}'|AA}, \ldots] are independent of [|{\bf H}|] in the range of measurement which is better fulfilled with neutrons than with X-rays (see below), γ, δ, [epsilon] are the coefficients of the Fourier series: [\eqalign{Q_{x} &= \textstyle\sum\limits_{\bf n} \textstyle\sum\limits_{{\bf n}'} i\gamma_{{\bf nn}'|x} \exp \{2 \pi i {\bf H} \cdot ({\bf n} - {\bf n}')\}\hbox{;}\cr Q_{y} &= \ldots\hbox{;} \qquad Q_{z} = \ldots\hbox{;}\cr R_{x} &= \textstyle\sum\limits_{\bf n} \textstyle\sum\limits_{{\bf n}'} \delta_{{\bf nn}'|x} \exp \{2 \pi i {\bf H} \cdot ({\bf n} - {\bf n}')\}\hbox{;}\cr R_{y} &= \ldots\hbox{;} \qquad R_{z} = \ldots\hbox{;}\cr S_{xy} &= \textstyle\sum\limits_{\bf n} \textstyle\sum\limits_{{\bf n}'} \varepsilon_{{\bf nn}'|xy} \exp \{2 \pi i {\bf H} \cdot ({\bf n} - {\bf n}')\}\hbox{;}\cr S_{yz} &= \ldots\hbox{;} \qquad S_{zx} = \ldots.}] The functions Q, R, S are then periodic in reciprocal space.

The double sums over n, n′ may be replaced by [N \sum_{m, \, n, \, p}] where m, n, p are the coordinates of the interatomic vectors [({\bf n} - {\bf n}')] and [I_{2}] becomes [\eqalignno{I_{2} &= - Nc_{A} c_{B} (\;f_{A} - f_{B})^{2} \textstyle\sum\limits_{m} \textstyle\sum\limits_{n} \textstyle\sum\limits_{p} (H\gamma_{lmn|x} + \ldots + \ldots) &\cr &\quad \times \sin 2 \pi (Hm + Kn + Lp). &(4.2.4.81)}] The intensity is therefore modulated sinusoidally and increases with scattering angle. The modulation gives rise to an asymmetry in the intensity around a Bragg peak. Similar considerations for [I_{3}] reveal an intensity contribution [h_{i}^{2}] times a sum over cosine terms which is symmetric around the Bragg peaks. This term shows quite an analogous influence of local static displacements and thermal movements: an increase of diffuse intensity around the Bragg peaks and a reduction of Bragg intensities, which is not discussed here. The second contribution [I_{2}] has no analogue owing to the non-vanishing average displacement. The various diffuse intensity contributions may be separated by symmetry considerations. Once they are separated, the single coefficients may be determined by Fourier inversion. Owing to the symmetry constraints there are relations between the displacements [\langle x \ldots \rangle] and, in turn, between the γ and Q components. The same is true for the δ, [epsilon], and R, S components. Consequently, there are symmetry conditions for the individual contributions of the diffuse intensity which may be used to distinguish them. Generally the total diffuse intensity may be split into only a few independent terms. The single components of Q, R, S may be expressed separately by combinations of diffuse intensities which are measured in definite selected volumes in reciprocal space. Only a minimum volume must be explored in order to reveal the behaviour over the whole reciprocal space. This minimum repeat volume is different for the single components: [I_{\rm sro}], Q, R, S or combinations of them.

The Borie–Sparks method has been applied very frequently to binary and even ternary systems; some improvements have been communicated by Bardhan & Cohen (1976)[link]. The diffuse scattering of the historically important metallic compound Cu3Au has been studied by Cowley (1950a[link],b[link]), and the pair correlation parameters could be determined. The typical fourfold splitting was found by Moss (1966)[link] and explained in terms of atomic displacements. The same splitting has been found for many similar compounds such as Cu3Pd (Ohshima et al., 1976[link]), Au3Cu (Bessière et al., 1983[link]), and Ag1−xMgx (x = 0.15–0.20) (Ohshima & Harada, 1986[link]). Similar pair correlation functions have been determined. In order to demonstrate the disorder parameters in terms of structural models, computer programs were used (e.g. Gehlen & Cohen, 1965[link]). A similar microdomain model was proposed by Hashimoto (1974[link], 1981[link], 1983[link], 1987[link]). According to approximations made in the theoretical derivation the evaluation of diffuse scattering is generally restricted to an area in reciprocal space where the influence of displacements is of the same order of magnitude as that of the pair correlation function. The agreement between calculation and measurement is fairly good but it should be remembered that the amount and quality of the experimental information used is low. No residual factors are so far available; these would give an idea of the reliability of the results.

The more general case of a multi-component system with several atoms per lattice point was treated similarly by Hayakawa & Cohen (1975)[link]. Sources of error in the determination of the short-range-order coefficients are discussed by Gragg et al. (1973)[link]. In general the assumption of constant [F_{{\bf nn}'|AA},\ldots] produces an incomplete separation of the order- and displacement-dependent components of diffuse scattering. By an alternative method, by separation of the form factors from the Q, R, S functions and solving a large array of linear relationships by least-squares methods, the accuracy of the separation of the various contributions is improved (Tibbals, 1975[link]; Georgopoulos & Cohen, 1977[link]; Wu et al., 1983[link]). The method does not work for neutron diffraction. Also, the case of planar short-range order with corresponding diffuse intensity along rods in reciprocal space may be treated along the Borie & Sparks method (Ohshima & Moss, 1983[link]).

Multi-wavelength methods taking advantage of the variation of the structure factor near an absorption edge (anomalous dispersion) are discussed by Cenedese et al. (1984)[link]. The same authors show that in some cases the neutron method allows for a contrast variation by using samples with different isotope substitution.

4.2.4.4.6. Powder diffraction

| top | pdf |

Evaluation of diffuse-scattering data from powder diffraction follows the same theoretical formulae developed for the determination of the radial distribution function for glasses and liquids (Debye & Menke, 1931[link]; Warren & Gingrich, 1934[link]). The final formula for random distributions may be given as (Fender, 1973[link]) [{I_{d}^{p} = \{\langle |F ({\bf H})|^{2}\rangle - |\langle F ({\bf H})\rangle |^{2}\} \textstyle\sum\limits_{i} s_{i} \sin (2 \pi Hr_{i}) / (2 \pi Hr_{i}).} \eqno(4.2.4.82)] [s_{i}] represents the number of atoms at distance [r_{i}] from the origin. An equivalent expression for a substitutional binary alloy is [{I_{d}^{p} = \alpha (1 - \alpha) \{|\;f_{2} ({\bf H}) - f_{2} ({\bf H})|^{2}\} \textstyle\sum\limits_{i} s_{i} \sin (2 \pi Hr_{i}) / (2 \pi Hr_{i}).} \eqno(4.2.4.83)]

4.2.4.4.7. Small concentrations of defects

| top | pdf |

In the literature small concentrations are treated in terms of fluctuations of the functions [\pi_{{\bf n}\nu}] as defined in equation (4.2.4.56)[link]. Generally we prefer the introduction of the distribution function of the defects or clusters. Since this problem has already been treated in Section 4.2.4.4.3[link] only some very brief remarks are given here. The most convenient way to derive the distribution function correctly from experimental data is the use of low-angle scattering which generally shows one or more clear maxima caused by partly periodic properties of the distribution function. For the deconvolution of the distribution function, received by Fourier transformation of the corrected diffused low-angle scattering, the reader is referred to the relevant literature. Since deconvolutions are not unique some reasonable assumptions are necessary for the final solution. Anomalous scattering may be very helpful if applicable.

4.2.4.4.8. Cluster method

| top | pdf |

As mentioned above, the cluster method may be useful for the interpretation of disorder problems. In the general formula of diffuse scattering of random distributions equation (4.2.2.13)[link] may be used. Here [|\langle F ({\bf H})\rangle |^{2}] describes the sharp Bragg maxima, while [|\Delta F ({\bf H})|^{2} = \langle |F ({\bf H})|^{2}\rangle - |\langle F ({\bf H})\rangle |^{2}] represents the contribution to diffuse scattering. Correlation effects can also be taken into account by using clusters of sufficient size if their distribution may be considered as random in good approximation. The diffuse intensity is then given by [I_{d} ({\bf H}) = \textstyle\sum\limits_{\nu} p_{\nu} |F_{\nu} ({\bf H})|^{2} - |\textstyle\sum\limits_{\nu} p_{\nu} F_{\nu} ({\bf H})|^{2}, \eqno(4.2.4.84)] where [F_{\nu} ({\bf H})] represents the difference structure factor of the νth cluster and [p_{\nu}] is its a priori probability. Obviously equation (4.2.4.84)[link] is of some use in two cases only. (1) The number of clusters is sufficiently small and meets the condition of nearly random distribution. In principle, its structure may then be determined with the aid of refinement methods according to equation (4.2.4.84)[link]. Since the second term is assumed to be known from the average structure, the first term may be evaluated by introducing as many parameters as there are clusters involved. A special computer program for incoherent refinement has to be used if more than one representative cluster has to be introduced. In the case of more clusters, constraints are necessary. (2) The number of clusters with similar structures is not limited. It may be assumed that their size distribution may be expressed by well known analytical expressions, e.g. Gaussians or Lorentzians. The distribution is still assumed to be random.

An early application of the cluster method was the calculation of the diffuse intensity of Guinier–Preston zones, where a single cluster is sufficient (see, e.g., Gerold, 1954[link]; Bubeck & Gerold, 1984[link]). Unfortunately no refinements of cluster structures have so far been published. The full theory of the cluster method was outlined by Jagodzinski & Haefner (1967)[link].

Some remarks on the use of residual factors should be added here. Obviously the diffuse scattering may be used for refinements in a similar way as in conventional structure determination. For this purpose a sufficiently small reciprocal lattice has to be defined. The size of the reciprocal cell has to be chosen with respect to the maximum gradient of diffuse scattering. Then the diffuse intensity may be described by a product of the real intensity distribution and the small reciprocal lattice. Fourier transformation yields the convolution of the real disordered structure and a large unit cell. In other words, the disordered structure is subdivided into large units and subsequently superimposed (`projected') in a single cell. In cases where a clear model of the disorder could be determined, a refinement procedure for atomic and other relevant parameters can be started. In this way a residual factor may be determined. A first approach has been elaborated by Epstein & Welberry (1983)[link] in the case of substitutional disorder of two molecules.

The outstanding limiting factor is the collection of weak intensity data. The amount increases rapidly with the complexity of the structure and could even exceed by far the amount which is needed in the case of protein structure refinement. Hence, it seems to be reasonable to restrict the measurement to distinct areas in reciprocal space. Most of these publications, however, use too little information when compared with the minimum of data which would be necessary for the confirmation of the proposed model. Hence, physical and chemical considerations should be used as an additional source of information.

4.2.4.4.9. Comparison between X-ray and neutron methods

| top | pdf |

Apart from experimental arguments in favour of either method, there are some specific points which should be mentioned in this context. The diffuse scattering in question must be separated from Bragg scattering and from other diffuse-scattering contributions. Generally both methods are complementary: neutrons are preferable in cases where X-rays show only a small scattering contrast: (heavy) metal hydrides, oxides, carbides, Al–Mg distribution etc. In favourable cases it is possible to suppress (nuclear) Bragg scattering of neutrons when isotopes are used so that [\sum_{\nu} c_{\nu} f_{\nu} = 0] for all equivalent positions. Another way to separate Bragg peaks is to record the diffuse intensity, if possible, at low [|{\bf H}|] values. This can be achieved either by measurement at low θ angles or by using long wavelengths. For reasons of absorption the latter point is the domain of neutron scattering. Exceeding the Bragg cut-off, Bragg scattering is ruled out. In this way `diffuse' background owing to multiple Bragg scattering is avoided. Other diffuse-scattering contributions which increase with the [|{\bf H}|] value are thus also minimized: thermal diffuse scattering (TDS) and scattering due to long-range static displacements. On the other hand, lattice distortions, Huang scattering, … should be measured at large values of [|{\bf H}|]. TDS can be separated by purely elastic neutron methods within the limits given by the energy resolution of an instrument. This technique is of particular importance at higher temperatures where TDS becomes remarkably strong. Neutron scattering is a good tool only in cases where (isotope/spin-)incoherent scattering is not too strong. In the case of magnetic materials confusion with paramagnetic diffuse scattering could occur. This is also important when electrons are trapped by defects which themselves act as paramagnetic centres.

As mentioned in Section 4.2.4.4.4[link] the evaluations of the γ, δ, [epsilon] depend on the assumption that the f's do not depend on [|{\bf H}|] strongly within the range of measurement. Owing to the atomic form factor, this is not always well approximated in the X-ray case and is one of the main sources of error in the determination of the short-range-order parameters.

4.2.4.4.10. Dynamic properties of defects

| top | pdf |

Some brief remarks concerning the dynamic properties of defects as discussed in the previous sections now follow. Mass defects (impurity atoms), force-constant defects etc. influence the dynamic properties of the undistorted lattice and one could think of a modified TDS as discussed in Chapter 4.1[link] . In the case of low defect concentrations special vibrational modes characterized by large amplitudes at the defect with frequency shifts and reduced lifetimes (resonant modes) or vibrational modes localized in space may occur. Other modes with frequencies near these particular modes may also be affected. Owing to the very low intensity of these phenomena their influence on the normal TDS is negligible and may be neglected in diffuse-scattering work. Theoretical treatments of crystals with higher defect concentrations are extremely difficult and not developed so far. For further reading see Böttger (1983)[link].

4.2.4.5. Orientational disorder

| top | pdf |

Molecular crystals show in principle disorder phenomena similar to those discussed in previous sections (substitutional or displacement disorder). Here we have to replace the structure factors [F_{\nu}({\bf H})], used in the previous sections, by the molecular structure factors in their various orientations. Usually these are rapidly varying functions in reciprocal space which may obscure the disorder diffuse scattering. Disorder in molecular crystals is treated by Guinier (1963)[link], Amorós & Amorós (1968)[link], Flack (1970)[link], Epstein et al. (1982)[link], Welberry & Siripitayananon (1986[link], 1987[link]), and others.

A particular type of disorder is very common in molecular and also in ionic crystals: the centres of masses of molecules or ionic complexes form a perfect 3D lattice but their orientations are disordered. Sometimes these solids are called plastic crystals. For comparison, the liquid-crystalline state is characterized by an orientational order in the absence of long-range positional order of the centres of the molecules. A clear-cut separation is not possible in cases where translational symmetry occurs in low dimension, e.g. in sheets or parallel to a few directions in crystal space. For discussion of these mesophases see Chapter 4.4[link] .

An orientationally disordered crystal may be imagined in a static picture by freezing molecules in different sites in one of several orientations. Local correlations between neighbouring molecules and correlations between position and orientation may be responsible for orientational short-range order. Often thermal reorientations of the molecules are related to an orientationally disordered crystal. Thermal vibrations of the centres of masses of the molecules, librational or rotational excitations around one or more axes of the molecules, jumps between different equilibrium positions or diffusion-like phenomena are responsible for diffuse scattering of dynamic origin. As mentioned above the complexity of molecular structures and the associated large number of thermal modes complicate a separation from static disorder effects.

Generally high Debye–Waller factors are typical for scattering of orientationally disordered crystals. Consequently only a few Bragg reflections are observable. A large amount of structural information is stored in the diffuse background. It has to be analysed with respect to an incoherent and coherent part, elastic, quasielastic or inelastic nature, short-range correlations within one and the same molecule and between orientations of different molecules, and cross correlations between positional and orientational disorder scattering. Combined X-ray and neutron methods are therefore highly recommended.

4.2.4.5.1. General expressions

| top | pdf |

On the assumption of a well ordered 3D lattice, a general expression for the scattering by an orientationally disordered crystal with one molecule per unit cell may be given. This is a very common situation. Moreover, orientational disorder is frequently related to molecules with an overall `globular' shape and consequently to crystals of high (in particular, averaged) spherical symmetry. In the following the relevant equations are given for this situation; these are discussed in some detail in a review article by Fouret (1979)[link]. The orientation of a molecule is characterized by a parameter [\omega_{l}], e.g. the set of Eulerian angles of three molecular axes with respect to the crystal axes: [\omega_{l} = 1, \ldots, D] (D possible different orientations). The equilibrium position of the centre of mass of a molecule in orientation [\omega_{l}] is given by [{\bf r}_{l}], the equilibrium position of atom k within a molecule l in orientation [\omega_{l}] by [{\bf r}_{lk}] and a displacement from this equilibrium position by [{\bf u}_{lk}]. Averaging over a long time, i.e. supposing that the lifetime of a discrete configuration is long compared with the period of atomic vibrations, the observed intensity may be deduced from the intensity expression corresponding to a given configuration at time t: [\eqalignno{I({\bf H}, t) &= \textstyle\sum\limits_{l} \textstyle\sum\limits_{l'} F_{l}({\bf H}, t) F_{l'}^{+}({\bf H}, t) &\cr &\quad \times \exp \{2 \pi i{\bf H} \cdot ({\bf r}_{l} - {\bf r}_{l'})\} &(4.2.4.85)\cr F_{l}({\bf H}, t) &= \textstyle\sum\limits_{k} f_{k} \exp \{2 \pi i{\bf H} \cdot ({\bf r}_{lk} + {\bf u}_{lk})\}. &(4.2.4.86)}%(4.2.4.86)] Averaging procedures must be carried out with respect to the thermal vibrations (denoted by an overbar) and over all configurations (symbol [\langle\; \rangle]). The centre-of-mass translational vibrations and librations of the molecules are most important in this context. (Internal vibrations of the molecules are assumed to be decoupled and remain unconsidered.) [\eqalignno{I({\bf H}, t) &= \textstyle\sum\limits_{l} \textstyle\sum\limits_{l'} {\langle \overline{F_{l}({\bf H}, t) F_{l'}^{+}({\bf H}, t)}\rangle} &\cr &\quad \times \exp \{2 \pi i{\bf H} \cdot ({\bf r}_{l} - {\bf r}_{l'})\}. &(4.2.4.85a)}] Thermal averaging gives (cf. Chapter 4.1[link] ) [\eqalignno{I &= \textstyle\sum\limits_{l} \textstyle\sum\limits_{l'} \overline{F_{l}F_{l'}^{+}} \exp \{2 \pi i{\bf H} \cdot ({\bf r}_{l} - {\bf r}_{l'})\} &\cr \overline{F_{l}F_{l'}^{+}} &= \textstyle\sum\limits_{k} \textstyle\sum\limits_{k'} f_{k}\;f_{k'} \exp \{2 \pi i{\bf H} \cdot ({\bf r}_{lk} - {\bf r}_{l'k'})\} &\cr &\quad \times \overline{\exp \{2 \pi i{\bf H} \cdot ({\bf u}_{lk} - {\bf u}_{l'k'})\}}. &(4.2.4.87)}] In the harmonic approximation [\exp\{\overline{2 \pi i{\bf H} \cdot \Delta {\bf u}}\}] is replaced by [\exp\{{1\over 2}\overline{|2 \pi {\bf H} \cdot \Delta {\bf u}|^{2}}\}]. This is, however, a more or less crude approximation because strongly anharmonic vibrations are quite common in an orientationally disordered crystal. In this approximation [\overline{F_{l}F_{l'}^{+}}] becomes [\eqalignno{\overline{F_{l}F_{l'}^{+}} &= \textstyle\sum\limits_{k} \textstyle\sum\limits_{k'} f_{k}\;f_{k'} \exp \{- B_{k}(\omega_{l})\} &\cr &\quad \times \exp \{- B_{k'}(\omega_{l'})\} \exp \{D_{lk;\ l'k'}\}. &(4.2.4.88)}] [B_{k}] is equal to [{1\over 2}\overline{(2 \pi {\bf H} \cdot {\bf u}_{lk})^{2}}] (Debye–Waller factor) and depends on the specific configuration [\omega_{l}]. [D_{lk;\ l'k'} = \overline{(2 \pi {\bf H} \cdot {\bf u}_{lk})(2 \pi {\bf H} \cdot {\bf u}_{l'k'})}] includes all the correlations between positions, orientations and vibrations of the molecules.

Averaging over different configurations demands a knowledge of the orientational probabilities. The probability of finding molecule l in orientation [\omega_{l}] is given by [p(\omega_{l})]. The double probability [p(\omega_{l}, \omega_{l'})] gives the probability of finding two molecules [l, l'] in different orientations [\omega_{l}] and [\omega_{l'}], respectively. In the absence of correlations between the orientations we have: [p(\omega_{l}, \omega_{l'}) = p(\omega_{l})p(\omega_{l'})]. If correlations exist: [p(\omega_{l}, \omega_{l'}) = p(\omega_{l})p'(\omega_{l} | \omega_{l'})] where [p'(\omega_{l} | \omega_{l'})] defines the conditional probability that molecule [l'] has the orientation [\omega_{l}] if molecule l has the orientation [\omega_{l'}]. For long distances between l and [l'\ p'(\omega_{l} | \omega_{l'})] tends to [p(\omega_{l'})].

The difference [\Delta (\omega_{l} | \omega_{l'}) = p'(\omega_{l} | \omega_{l'}) - p(\omega_{l'})] characterizes, therefore, the degree of short-range orientational correlation. Note that this formalism corresponds fully to the [p_{\mu}], [p_{\mu \mu'}] used in the context of translational disorder.

The average structure factor, sometimes called averaged form factor, of the molecule is given by [\langle F_{l}\rangle = \textstyle\sum\limits_{\omega_{l}} p(\omega_{l}) F_{l}(\omega_{l}). \eqno(4.2.4.89)]

  • (a) Negligible correlations between vibrations of different molecules (Einstein model): [D_{lk;\ l'k'} = 0 \hbox{ for } l \neq l'.] From (4.2.4.88)[link] it follows (the prime symbol takes the Debye–Waller factor into account): [\eqalignno{\overline{\langle I\rangle} &= N^{2}|\langle F'\rangle|^{2} L({\bf H}) &\cr &\quad + N \Bigl\{\textstyle\sum\limits_{k} \textstyle\sum\limits_{k'} \textstyle\sum\limits_{\omega_{l}} p(\omega_{l}) f_{k}f_{k'} \exp \bigl\{2 \pi i{\bf H} \cdot ({\bf r}_{lk} - {\bf r}_{lk'})\bigr\} &\cr &\quad \times \exp \{D_{lk;\ lk'}\} - |\langle F'\rangle|^{2}\Bigr\} &\cr &\quad + N \textstyle\sum\limits_{\Delta l\neq 0} \textstyle\sum\limits_{\omega_{l}} \textstyle\sum\limits_{\omega_{l'}} p(\omega_{l}) \Delta (\omega_{l} | \omega_{l'}) &\cr &\quad \times F'_{l}(\omega_{l}) F_{l'}^{+}(\omega_{l'}) &\cr &\quad \times \exp \{2 \pi i{\bf H} \cdot ({\bf r}_{l} - {\bf r}_{l'})\}. &(4.2.4.90)}] [L({\bf H})] is the reciprocal lattice of the well defined ordered lattice. The first term describes Bragg scattering from an averaged structure. The second term governs the diffuse scattering in the absence of short-range orientational correlations. The last term takes the correlation between the orientations into account.

    If rigid molecules with centre-of-mass translational displacements and negligible librations are assumed, which is a first approximation only, [|\langle F\rangle|^{2}] is no longer affected by a Debye–Waller factor.

    In this approximation the diffuse scattering may therefore be separated into two parts: [{N(\langle F^{2}\rangle - |\langle F'\rangle|^{2}) = N(F^{2} - |\langle F\rangle|^{2}) + N(|\langle F\rangle|^{2} - |\langle F'\rangle|^{2})} \eqno(4.2.4.91)] with [\eqalignno{\langle F^{2}\rangle &= \textstyle\sum\limits_{\omega_{l}} \textstyle\sum\limits_{\omega_{l'}} \textstyle\sum\limits_{k} \textstyle\sum\limits_{k'} f_{k}(\omega_{l}) f_{k'}(\omega_{l'}) p(\omega_{l}) &\cr &\quad \times \exp \{2 \pi i{\bf H} \cdot ({\bf r}_{lk} - {\bf r}_{l'k'})\}. &(4.2.4.92)}] The first term in (4.2.4.91)[link] gives the scattering from equilibrium fluctuations in the scattering from individual molecules (diffuse scattering without correlations), the second gives the contribution from the centre-of-mass thermal vibrations of the molecules.

  • (b) If intermolecular correlations between the molecules cannot be neglected, the final intensity expression for diffuse scattering is very complicated. In many cases these correlations are caused by dynamical processes (see Chapter 4.1[link] ). A simplified treatment assumes the molecule to be a rigid body with a centre-of-mass displacement [{\bf u}_{l}] and neglects vibrational–librational and librational–librational correlations: [D_{l;\ l'} = \overline{(2 \pi {\bf H} \cdot {\bf u}_{l})(2 \pi {\bf H} \cdot {\bf u}_{l'})}] [(l \neq l')]. The following expression approximately holds: [\eqalignno{\overline{\langle I\rangle} &= N^{2}|\langle F'\rangle|^{2} L({\bf H}) &\cr &\quad + \langle \textstyle\sum\limits_{l} \textstyle\sum\limits_{l'} F'_{l}(\omega_{l}) F_{l'}^{'+}(\omega_{l'}) \exp \{D_{l;\ l'}\}\rangle \{2 \pi i{\bf H} \cdot ({\bf r}_{l} - {\bf r}_{l'})\} &\cr &\quad + N \Big\{\textstyle\sum\limits_{\omega_{l}} \textstyle\sum\limits_{k, \, k'} p(\omega_{l}) f_{k}\;f_{k'} \exp \{2 \pi i{\bf H} \cdot ({\bf r}_{lk} - {\bf r}_{l'k'})\} &\cr &\quad \times \exp \{D_{lk;\ l'k'}\} - \textstyle\sum\limits_{\omega_{l}} \textstyle\sum\limits_{\omega_{l'}} \textstyle\sum\limits_{k} \textstyle\sum\limits_{k'} p(\omega_{l}) p(\omega_{l'}) f_{k}\;f_{k'} &\cr &\quad \times \exp \{2 \pi i{\bf H} \cdot ({\bf r}_{lk} - {\bf r}_{l'k'})\} &\cr &\quad \vphantom{\textstyle\sum\limits_{l\neq l'}} \times \exp \{D_{lk;\ l'k'}\}\Big\} + \textstyle\sum\limits_{l\neq l'} \textstyle\sum\limits_{\omega_{l}} \sum\limits_{\omega_{l}} p(\omega_{l}) \Delta (\omega_{l} | \omega_{l'}) &\cr &\quad \times F'_{l}(\omega_{l}) F_{l'}^{'+}(\omega_{l'}) &\cr &\quad \times \exp \{2 \pi i{\bf H} \cdot ({\bf r}_{l} - {\bf r}_{l'})\} \exp \{D_{l;\ l'}\}. &(4.2.4.93)}] Again the first term describes Bragg scattering and the second corresponds to the average thermal diffuse scattering in the disordered crystal. Because just one molecule belongs to one unit cell only acoustic waves contribute to this part. To an approximation, the result for an ordered crystal may be used by replacing F by [\langle F'\rangle] [Chapter 4.1, equation (4.1.3.4)[link] ]. The third term corresponds to random-disorder diffuse scattering. If librations are neglected this term may be replaced by [N(\langle F^{2}\rangle - \langle F\rangle^{2})]. The last term in (4.2.4.93)[link] describes space correlations. Omission of [\exp\{D_{l;\ l'}\}] or expansion to [\sim (1 + D_{l;\ l'})] are further simplifying approximations.

    In either (4.2.4.90)[link] or (4.2.4.93)[link] the diffuse-scattering part depends on a knowledge of the conditional probability [\Delta (\omega_{l} | \omega_{l'})] and the orientational probability [p(\omega_{l})]. The latter may be found, at least in principle, from the average structure factor.

4.2.4.5.2. Rotational structure (form) factor

| top | pdf |

In certain cases and with simplifying assumptions, [\langle F\rangle] [equation (4.2.4.89)[link]] and [\langle \Delta F^{2}\rangle] [equation (4.2.4.92)[link]] may be calculated. Assuming only one molecule per unit cell and treating the molecule as a rigid body, one derives from the structure factor of an ordered crystal [F_{l}] [\langle F\rangle = \textstyle\sum\limits_{k} f_{k} \langle \exp \{2 \pi i{\bf H} \cdot {\bf r}_{lk}\}\rangle \eqno(4.2.4.94)] and [\eqalignno{\langle \Delta F^{2}\rangle &= \textstyle\sum\limits_{k} \textstyle\sum\limits_{k'} f_{k}\;f_{k'} [\langle \exp \{2 \pi i{\bf H} \cdot ({\bf r}_{lk} - {\bf r}_{l'k'})\}\rangle &\cr &\quad - \langle \exp \{2 \pi i{\bf H} \cdot {\bf r}_{lk}\}\rangle \langle \exp \{2 \pi i{\bf H} \cdot {\bf r}_{l'k'}\}\rangle]. &(4.2.4.95)}] If the molecules have random orientation in space the following expressions hold [see, e.g., Dolling et al. (1979)[link]]: [\langle F\rangle = \textstyle\sum\limits_{k} f_{k}\;j_{0}({\bf H} \cdot {\bf r}_{k}) \eqno(4.2.4.96)] [\eqalignno{\langle |\Delta F|^{2}\rangle &= \textstyle\sum\limits_{k} \textstyle\sum\limits_{k'} f_{k}\;f_{k'} \{\;j_{0} [{\bf H} \cdot ({\bf r}_{k} - {\bf r}_{k'})] &\cr &\quad - j_{0}({\bf H} \cdot {\bf r}_{k})j_{0}({\bf H} \cdot {\bf r}_{k'})\}. &(4.2.4.97)}] [j_{0}(z)] is the zeroth order of the spherical Bessel functions and describes an atom k uniformly distributed over a shell of radius [r_{k}].

In practice the molecules perform more or less finite librations about the main orientation. The structure factor may then be found by the method of symmetry-adapted functions [see, e.g., Press (1973)[link], Press & Hüller (1973)[link], Dolling et al. (1979)[link], Prandl (1981[link], and references therein)]. [\langle F\rangle = \textstyle\sum\limits_{k} f_{k} 4 \pi \textstyle\sum\limits_{\nu} \textstyle\sum\limits_{\mu =- \nu}^{+ \nu} i^{\nu}j_{\nu}({\bf H} \cdot {\bf r}_{k}) C_{\nu \mu}^{(k)} Y_{\nu \mu}(\theta, \varphi). \eqno(4.2.4.98)] [j_{\nu}(z)] is the νth order of spherical Bessel functions, the coefficients [C_{\nu \mu}^{(k)}] characterize the angular distribution of [{\bf r}_{k}], [Y(\theta, \varphi)] are the spherical harmonics where [|{\bf H}|, \theta, \varphi] denote polar coordinates of H.

The general case of an arbitrary crystal, site and molecular symmetry and the case of several symmetrically equivalent orientationally disordered molecules per unit cell are treated by Prandl (1981)[link]; an example is given by Hohlwein et al. (1986)[link]. As mentioned above, cubic plastic crystals are common and therefore mostly studied up to now. The expression for [\langle F\rangle] may then be formulated as an expansion in cubic harmonics, [K_{\nu \mu}(\theta, \varphi)]: [\langle F\rangle = \textstyle\sum\limits_{k} f_{k} 4 \pi \textstyle\sum\limits_{\nu} \textstyle\sum\limits_{\mu} i^{\nu} j_{\nu}({\bf H} \cdot {\bf r}_{k}) C_{\nu \mu}^{'(k)} K_{\nu \mu}(\theta, \varphi). \eqno(4.2.4.99)] ([C'_{\nu \mu}] are modified expansion coefficients.)

Taking into account isotropic centre-of-mass translational displacements, which are not correlated with the librations, we obtain: [\langle F'\rangle = \langle F\rangle \exp \{- {\textstyle{1 \over 6}} H^{2} \langle U^{2}\rangle \}. \eqno(4.2.4.100)] U is the mean-square translational displacement of the molecule. Correlations between translational and vibrational displacements are treated by Press et al. (1979)[link].

Equivalent expressions for crystals with symmetry other than cubic may be found from the same concept of symmetry-adapted functions [tables are given by Bradley & Cracknell (1972)[link]].

4.2.4.5.3. Short-range correlations

| top | pdf |

The final terms in equations (4.2.4.90)[link] and (4.2.4.93)[link] concern correlations between the orientations of different molecules. Detailed evaluations need a knowledge of a particular model. Examples are compounds with nitrate groups (Wong et al., 1984[link]; Lefebvre et al., 1984[link]), CBr4 (More et al., 1980[link], 1984[link]), and many others (see Sherwood, 1979[link]). The situation is even more complicated when a modulation wave with respect to the occupation of different molecular orientations is superimposed. A limiting case would be a box-like function describing a pattern of domains. Within one domain all molecules have the same orientation. This situation is common in ferroelectrics where molecules exhibit a permanent dipole moment. The modulation may occur in one or more directions in space. The observed intensity in this type of orientationally disordered crystal is characterized by a system of more or less diffuse satellite reflections. The general scattering theory of a crystal with occupational modulation waves follows the same lines as outlined in Section 4.2.3.1.[link]

References

Amorós, J. L. & Amorós, M. (1968). Molecular crystals: their transforms and diffuse scattering. New York: John Wiley.
Bardhan, P. & Cohen, J. B. (1976). X-ray diffraction study of short-range-order structure in a disordered Au3Cu alloy. Acta Cryst. A32, 597–614.
Bauer, G. S. (1979). Diffuse elastic neutron scattering from nonmagnetic materials. In Treatise on materials science and technology, Vol. 15, edited by G. Kostorz, pp. 291–336. New York: Academic Press.
Bessière, M., Lefebvre, S. & Calvayrac, Y. (1983). X-ray diffraction study of short-range order in a disordered Au3Cu alloy. Acta Cryst. B39, 145–153.
Beyeler, H. U., Pietronero, L. & Strässler, S. (1980). Configurational model for a one-dimensional ionic conductor. Phys. Rev. B, 22, 2988–3000.
Borie, B. & Sparks, C. J. (1971). The interpretation of intensity distributions from disordered binary alloys. Acta Cryst. A27, 198–201.
Bradley, C. J. & Cracknell, A. P. (1972). The mathematical theory of symmetry in solids, pp. 51–76. Oxford: Clarendon Press.
Brämer, R. (1975). Statistische Probleme der Theorie des Parakristalls. Acta Cryst. A31, 551–560.
Brämer, R. & Ruland, W. (1976). The limitations of the paracrystalline model of disorder. Macromol. Chem. 177, 3601–3617.
Bubeck, E. & Gerold, V. (1984). An X-ray investigation in the small and wide angle range from G.P.I. zones in Al–Cu. In Microstructural characterization of materials by non-microscopical techniques, edited by N. H. Andersen, M. Eldrup, N. Hansen, R. J. Jensen, T. Leffers, H. Lilholt, O. B. Pedersen & B. N. Singh. Roskilde: Risø National Laboratory.
Böttger, H. (1983). Principles of the theory of lattice dynamics. Weinheim: Physik Verlag.
Cenedese, P., Bley, F. & Lefebvre, S. (1984). Diffuse scattering in disordered ternary alloys: neutron measurements of local order in a stainless steel Fe0.56Cr0.21Ni0.23. Acta Cryst. A40, 228–240.
Comes, R. & Shirane, G. (1979). X-ray and neutron scattering from one-dimensional conductors. In Highly conducting one dimensional solids, ch. 2, edited by J. T. Devreese, R. P. Evrard & V. E. Van Doren. New York: Plenum.
Conradi, E. & Müller, U. (1986). Fehlordnung bei Verbindungen mit Schichtstrukturen. II. Analyse der Fehlordnung in Wismuttriiodid. Z. Kristallogr. 176, 263–269.
Cowley, J. M. (1950a). X-ray measurement of order in single crystals of Cu3Au. J. Appl. Phys. 21, 24–36.
Cowley, J. M. (1950b). An approximate theory of order in alloys. Phys. Rev. 77, 664–675.
Cowley, J. M. (1976a). Diffraction by crystals with planar faults. I. General theory. Acta Cryst. A32, 83–87.
Cowley, J. M. (1976b). Diffraction by crystals with planar faults. II. Magnesium fluorogermanate. Acta Cryst. A32, 88–91.
Cowley, J. M. (1981). Diffraction physics, 2nd ed. Amsterdam: North-Holland.
Cowley, J. M., Cohen, J. B., Salamon, M. B. & Wuensch, B. J. (1979). Modulated structures. AIP Conf. Proc. No. 53.
Debye, B. & Menke, H. (1931). Untersuchung der molekularen Ordnung in Flüssigkeiten mit Röntgenstrahlung. Ergeb. Tech. Roentgenkd. II, 1–22.
Dederichs, P. H. (1973). The theory of diffuse X-ray scattering and its application to the study of point defects and their clusters. J. Phys. F, 3, 471–496.
Dolling, G., Powell, B. M. & Sears, V. F. (1979). Neutron diffraction study of the plastic phases of polycrystalline SF6 and CBr4. Mol. Phys. 37, 1859–1883.
Dorner, C. & Jagodzinski, H. (1972). Entmischung im System SnO2–TiO2. Krist. Tech. 7, 427–444.
Dubernat, J. & Pezerat, H. (1974). Fautes d'empilement dans les oxalates dihydratés des métaux divalents de la série magnésienne (Mg, Fe, Co, Ni, Zn, Mn). J. Appl. Cryst. 7, 387–393.
Edwards, O. S. & Lipson, H. (1942). Imperfections in the structure of cobalt. I. Experimental work and proposed structure. Proc. R. Soc. London Ser. A, 180, 268–277.
Emery, V. J. & Axe, J. D. (1978). One-dimensional fluctuations and the chain-ordering transformation in Hg3 − δAsF6. Phys. Rev. Lett. 40, 1507–1511.
Endres, H., Pouget, J. P. & Comes, R. (1982). Diffuse X-ray scattering and order–disorder effects in the iodine chain compounds N,N′-diethyl-N,N′-dihydrophenazinium iodide, E2PI1.6 and N,N′-dibenzyl-N,N′-dihydrophenazinium iodide, B2PI1.6. J. Phys. Chem. Solids, 43, 739–748.
Epstein, J. & Welberry, T. R. (1983). Least-squares analysis of diffuse scattering from substitutionally disordered crystals: application to 2,3-dichloro-6,7-dimethylanthracene. Acta Cryst. A39, 882–892.
Epstein, J., Welberry, T. R. & Jones, R. (1982). Analysis of the diffuse X-ray scattering from substitutionally disordered molecular crystals of monoclinic 9-bromo-10-methylanthracene. Acta Cryst. A38, 611–618.
Fender, B. E. F. (1973). Diffuse scattering and the study of defect solids. In Chemical applications of thermal neutron scattering, ch. 11, edited by B. T. M. Willis. Oxford University Press.
Flack, H. D. (1970). Short-range order in crystals of anthrone and in mixed crystals of anthrone–anthraquinone. Philos. Trans. R. Soc. London Ser. A, 266, 583–591.
Fontaine, D. de (1972). An analysis of clustering and ordering in multicomponent solid solutions. I. Stability criteria. J. Phys. Chem. Solids, 33, 297–310.
Fontaine, D. de (1973). An analysis of clustering and ordering in multicomponent solid solutions. II. Fluctuations and kinetics. J. Phys. Chem. Solids, 34, 1285–1304.
Forst, R., Jagodzinski, H., Boysen, H. & Frey, F. (1987). Diffuse scattering and disorder in urea inclusion compounds OC(NH2)2 + CnH2n + 2. Acta Cryst. B43, 187–197.
Fouret, P. (1979). Diffuse X-ray scattering by orientationally disordered solids. In The plastically crystalline state, ch. 3, edited by J. N. Sherwood. New York: John Wiley.
Frey, F. & Boysen, H. (1981). Disorder in cobalt single crystals. Acta Cryst. A37, 819–826.
Frey, F., Jagodzinski, H. & Steger, W. (1986). On the phase transformation zinc blende to wurtzite. Bull. Minéral. Crystallogr. 109, 117–129.
Gehlen, P. & Cohen, J. B. (1965). Computer simulation of the structure associated with local order in alloys. Phys. Rev. A, 139, 844–855.
Georgopoulos, P. & Cohen, J. B. (1977). The determination of short range order and local atomic displacements in disordered binary solid solutions. J. Phys. (Paris) Colloq. 38(C7), 191–196.
Gerold, V. (1954). Röntgenographische Untersuchungen über die Aushärtung einer Aluminium–Kupfer-Legierung mit Kleinwinkel-Schwenkaufnahmen. Z. Metallkd. 45, 593–607.
Gragg, J. E., Hayakawa, M. & Cohen, J. B. (1973). Errors in qualitative analysis of diffuse scattering from alloys. J. Appl. Cryst. 6, 59–66.
Guinier, A. (1963). X-ray diffraction in crystals, imperfect solids and amorphous bodies. San Francisco: Freeman.
Halla, F., Jagodzinski, H. & Ruston, W. R. (1953). One-dimensional disorder in dodecahydrotriphenylene, C18H24. Acta Cryst. 6, 478–488.
Hashimoto, S. (1974). Correlative microdomain model for short range ordered alloy structures. I. Diffraction theory. Acta Cryst. A30, 792–798.
Hashimoto, S. (1981). Correlative microdomain model for short range ordered alloy structures. II. Application to special cases. Acta Cryst. A37, 511–516.
Hashimoto, S. (1983). Correlative microdomain model for short range ordered alloy structures. III. Analysis for diffuse scattering from quenched CuAu alloy. Acta Cryst. A39, 524–530.
Hashimoto, S. (1987). Intensity expression for short-range order diffuse scattering with ordering energies in a ternary alloy system. J. Appl. Cryst. 20, 182–186.
Hayakawa, M. & Cohen, J. B. (1975). Experimental considerations in measurements of diffuse scattering. Acta Cryst. A31, 635–645.
Hendricks, S. B. & Teller, E. (1942). X-ray interference in partially ordered layer lattices. J. Chem. Phys. 10, 147–167.
Hohlwein, D., Hoser, A. & Prandl, W. (1986). Orientational disorder in cubic CsNO2 by neutron powder diffraction. Z. Kristallogr. 177, 93–102.
Hosemann, R. (1975). Micro paracrystallites and paracrystalline superstructures. Macromol. Chem. Suppl. 1, pp. 559–577.
Hosemann, R. & Bagchi, S. N. (1962). Direct analysis of diffraction by matter. Amsterdam: North-Holland.
Ishii, T. (1983). Static structure factor of Frenkel–Kontorova-systems at high temperatures. Application to K-hollandite. J. Phys. Soc. Jpn, 52, 4066–4073.
Jagodzinski, H. (1949a). Eindimensionale Fehlordnung und ihr Einfluß auf die Röntgeninterferenzen. I. Berechnung des Fehlordnungsgrades aus den Röntgeninterferenzen. Acta Cryst. 2, 201–208.
Jagodzinski, H. (1949b). Eindimensionale Fehlordnung und ihr Einfluß auf die Röntgeninterferenzen. II. Berechnung der fehlgeordneten dichtesten Kugelpackungen mit Wechsel wirkungen der Reichweite 3. Acta Cryst. 2, 208–214.
Jagodzinski, H. (1949c). Eindimensionale Fehlordnung und ihr Einfluß auf die Röntgeninterferenzen. III. Vergleich der Berechungen mit experimentellen Ergebnissen. Acta Cryst. 2, 298–304.
Jagodzinski, H. (1954). Der Symmetrieeinfluß auf den allgemeinen Lösungsansatz eindimensionaler Fehlordnungsprobleme. Acta Cryst. 7, 17–25.
Jagodzinski, H. (1963). On disorder phenomena in crystals. In Crystallography and crystal perfection, edited by G. N. Ramachandran, pp. 177–188. London: Academic Press.
Jagodzinski, H. (1964b). Diffuse disorder scattering by crystals. In Advanced methods of crystallography, edited by G. N. Ramachandran, pp. 181–219. London: Academic Press.
Jagodzinski, H. (1972). Transformation from cubic to hexagonal silicon carbide as a solid state reaction. Kristallografiya, 16, 1235–1246. [Translated into English in Sov. Phys. Crystallogr. 16, 1081–1090.]
Jagodzinski, H. & Haefner, K. (1967). On order–disorder in ionic non-stoichiometric crystals. Z. Kristallogr. 125, 188–200.
Jagodzinski, H. & Hellner, E. (1956). Die eindimensionale Phasenumwandlung des RhSn2. Z. Kristallogr. 107, 124–149.
Jagodzinski, H. & Korekawa, M. (1965). Supersatelliten im Beugungsbild des Labradorits (Ca2Na)(Si2Al)2O8. Naturwissenschaften, 52, 640–641.
Jagodzinski, H. & Laves, R. (1947). Über die Deutung der Entmischungsvorgänge in Mischkristallen unter besonderer Berücksichtigung der Systeme Aluminium–Kupfer und Aluminium–Silber. Z. Metallkd. 40, 296–305.
Jefferey, J. W. (1953). Unusual X-ray diffraction effects from a crystal of wollastonite. Acta Cryst. 6, 821–826.
Kakinoki, J. & Komura, Y. (1954). Intensity of X-ray diffraction by a one-dimensionally disordered crystal. I. General derivation in the cases of `Reichweite' s = 0 and s = 1. J. Phys. Soc. Jpn, 9, 169–183.
Kakinoki, J. & Komura, Y. (1965). Diffraction by a one-dimensionally disordered crystal. I. The intensity equation. Acta Cryst. 19, 137–147.
Kitaigorodsky, A. I. (1984). Mixed crystals. Springer series in solid state science, Vol. 33, chs. 6.4, 6.5, 8.4. Berlin: Springer.
Korekawa, M. & Jagodzinski, H. (1967). Die Satelliten-reflexe des Labradorits. Schweiz. Mineral. Petrogr. Mitt. 47, 269–278.
Krivoglaz, M. A. (1969). Theory of X-ray and thermal neutron scattering by real crystals. Part I. New York: Plenum.
Lefebvre, J., Fouret, R. & Zeyen, C. (1984). Structure determination of sodium nitrate near the order–disorder phase transition. J. Phys. (Paris), 45, 1317–1327.
Mardix, S. & Steinberger, I. T. (1970). Tilt and structure transformation in ZnS. J. Appl. Phys. 41, 5339–5341.
Martorana, A., Marigo, A., Toniolo, L. & Zenetti, R. (1986). Stacking faults in the β-form of magnesium dichloride. Z. Kristallogr. 176, 1–12.
More, M., Lefebvre, J. & Hennion, B. (1984). Quasi-elastic coherent neutron scattering in the disordered phase of CBr4. Experimental evidence of local order and rotational dynamics of molecules. J. Phys. (Paris), 45, 303–307.
More, M., Lefebvre, J., Hennion, B., Powell, B. M. & Zeyen, C. (1980). Neutron diffuse scattering in the disordered phase of CBr4. I. Experimental. Elastic and quasi-elastic coherent scattering in single crystals. J. Phys. C, 13, 2833–2846.
Moss, S. C. (1966). Local order in solid alloys. In Local atomic arrangements studied by X-ray diffraction, ch. 3, edited by J. B. Cohen & J. E. Hilliard, pp. 95–122. New York: Gordon and Breach.
Müller, H. (1952). Die eindimensionale Umwandlung Zinkblende-Wurtzit und die dabei auftretenden Anomalien. Neues Jahrb. Mineral. Abh. 84, 43–76.
Ohshima, K. & Harada, J. (1986). X-ray diffraction study of short-range ordered structure in a disordered Ag–15 at.% Mg alloy. Acta Cryst. B42, 436–442.
Ohshima, K. & Moss, S. C. (1983). X-ray diffraction study of basal-(ab)-plane structure and diffuse scattering from silver atoms in disordered stage-2 Ag0.18TiS2. Acta Cryst. A39, 298–305.
Ohshima, K., Watanabe, D. & Harada, J. (1976). X-ray diffraction study of short-range order diffuse scattering from disordered Cu–29.8% Pd alloy. Acta Cryst. A32, 883–892.
Pandey, D., Lele, S. & Krishna, P. (1980a). X-ray diffraction from one dimensionally disordered 2H-crystals undergoing solid state transformation to the 6H structure. I. The layer displacement mechanism. Proc. R. Soc. London Ser. A, 369, 435–439.
Pandey, D., Lele, S. & Krishna, P. (1980b). X-ray diffraction from one dimensionally disordered 2H-crystals undergoing solid state transformation to the 6H structure. II. The deformation mechanism. Proc. R. Soc. London Ser. A, 369, 451–461.
Pandey, D., Lele, S. & Krishna, P. (1980c). X-ray diffraction from one dimensionally disordered 2H-crystals undergoing solid state transformation to the 6H structure. III. Comparison with experimental observations on SiC. Proc. R. Soc. London Ser. A, 369, 463–477.
Peisl, J. (1975). Diffuse X-ray scattering from the displacement field of point defects and defect clusters. J. Appl. Cryst. 8, 143–149.
Prandl, W. (1981). The structure factor of orientational disordered crystals: the case of arbitrary space, site, and molecular point-group. Acta Cryst. A37, 811–818.
Press, W. (1973). Analysis of orientational disordered structures. I. Method. Acta Cryst. A29, 252–256.
Press, W., Grimm, H. & Hüller, A. (1979). Analysis of orientational disordered structures. IV. Correlations between orientation and position of a molecule. Acta Cryst. A35, 881–885.
Press, W. & Hüller, A. (1973). Analysis of orientational disordered structures. II. Examples: Solid CD4, p-D2 and NDBr4. Acta Cryst. A29, 257–263.
Radons, W., Keller, J. & Geisel, T. (1983). Dynamical structure factor of a 1-d harmonic liquid: comparison of different approximation methods. Z. Phys. B, 50, 289–296.
Rosshirt, E., Frey, F., Boysen, H. & Jagodzinski, H. (1985). Chain ordering in E2PI1.6 (5,10-diethylphenazinium iodide). Acta Cryst. B41, 66–76.
Scaringe, P. R. & Ibers, J. A. (1979). Application of the matrix method to the calculation of diffuse scattering in linearly disordered crystals. Acta Cryst. A35, 803–810.
Schmatz, W. (1973). X-ray and neutron scattering on disordered crystals. In Treatise on materials science and technology, Vol. 2, ch. 3.1, edited by H. Hermans. New York: Academic Press.
Schmatz, W. (1983). Neutron scattering studies of lattice defects: static properties of defects. In Methods of experimental physics, solid state: nuclear physics, Vol. 21, ch. 3.1, edited by J. N. Mundy, S. J. Rothman, M. J. Fluss & L. C. Smedskajew. New York: Academic Press.
Schwartz, L. H. & Cohen, J. B. (1977). Diffraction from materials. New York: Academic Press.
Sherwood, J. N. (1979). The plastically crystalline state. New York: John Wiley.
Sparks, C. J. & Borie, B. (1966). Methods of analysis for diffuse X-ray scattering modulated by local order and atomic displacements. In Atomic arrangements studied by X-ray diffraction, ch. 1, edited by J. B. Cohen & J. E. Hilliard, pp. 5–50. New York: Gordon and Breach.
Takaki, Y. & Sakurai, K. (1976). Intensity of X-ray scattering from one-dimensionally disordered crystal having the multilayer average structure. Acta Cryst. A32, 657–663.
Tibbals, J. E. (1975). The separation of displacement and substitutional disorder scattering: a correction for structure factor ratio variation. J. Appl. Cryst. 8, 111–114.
Vainshtein, B. K. (1966). Diffraction of X-rays by chain molecules. Amsterdam: Elsevier.
Warren, B. E. (1969). X-ray diffraction. Reading: Addison-Wesley.
Warren, B. E. & Gingrich, N. S. (1934). Fourier integral analysis of X-ray powder patterns. Phys. Rev. 46, 368–372.
Welberry, T. R. (1985). Diffuse X-ray scattering and models of disorder. Rep. Prog. Phys. 48, 1543–1593.
Welberry, T. R. & Siripitayananon, J. (1986). Analysis of the diffuse scattering from disordered molecular crystals: application to 1,4-dibromo-2,5-diethyl-3,6-dimethylbenzene at 295 K. Acta Cryst. B42, 262–272.
Welberry, T. R. & Siripitayananon, J. (1987). Analysis of the diffuse scattering from disordered molecular crystals: application to 1,3-dibromo-2,5-diethyl-4,6-dimethylbenzene at 295 K. Acta Cryst. B43, 97–106.
Wilke, W. (1983). General lattice factor of the ideal paracrystal. Acta Cryst. A39, 864–867.
Wilson, A. J. C. (1942). Imperfections in the structure of cobalt. II. Mathematical treatment of proposed structure. Proc. R. Soc. London Ser. A, 180, 277–285.
Wong, S. F., Gillan, B. E. & Lucas, B. W. (1984). Single crystal disorder diffuse X-ray scattering from phase II ammonium nitrate, NH4NO3. Acta Cryst. B40, 342–346.
Wu, T. B., Matsubara, E. & Cohen, J. B. (1983). New procedures for qualitative studies of diffuse X-ray scattering. J. Appl. Cryst. 16, 407–414.
Young, R. A. (1975). Editor. International discussion meeting on studies of lattice distortions and local atomic arrangements by X-ray, neutron and electron diffraction. J. Appl. Cryst. 8, 79–191.
Zernike, F. & Prins, J. A. (1927). Die Beugung von Röntgenstrahlen in Flüssigkeiten als Effekt der Molekülanordnung. Z. Phys. 41, 184–194.








































to end of page
to top of page