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

International Tables for Crystallography (2010). Vol. B, ch. 4.6, pp. 598-602   | 1 | 2 |

Section 4.6.3.1. Incommensurately modulated structures (IMSs)

W. Steurera* and T. Haibacha

aLaboratory of Crystallography, Department of Materials, ETH Hönggerberg, HCI G 511, Wolfgang-Pauli-Strasse 10, CH-8093 Zurich, Switzerland
Correspondence e-mail:  walter.steurer@mat.ethz.ch

4.6.3.1. Incommensurately modulated structures (IMSs)

| top | pdf |

One-dimensionally modulated structures are the simplest representatives of IMSs. The vast majority of the one hundred or so IMSs known so far belong to this class (Cummins, 1990[link]). However, there is also an increasing number of IMSs with 2D or 3D modulation. The dimension d of the modulation is defined by the number of rationally independent modulation wavevectors (satellite vectors) [{\bf q}_{i}] (Fig. 4.6.3.1)[link]. The electron-density function of a dD modulated 3D crystal can be represented by the Fourier series [\rho ({\bf r}) = (1/V) {\textstyle\sum\limits_{\bf H}} F({\bf H}) \exp (- 2\pi i {\bf H} \cdot {\bf r}).]The Fourier coefficients (structure factors) [F({\bf H})] differ from zero only for reciprocal-space vectors [{\bf H} = {\textstyle\sum_{i = 1}^{3}} h_{i}{\bf a}_{i}^{*} + {\textstyle\sum_{j = 1}^{d}} m_{j}{\bf q}_{j} =] [{\textstyle\sum_{i = 1}^{3 + d}} h_{i}{\bf a}_{i}^{*}] with [h_{i}, m_{j} \in {\bb Z}]. The d satellite vectors are given by [{\bf q}_{j} = {\bf a}_{3 + j}^{*} = {\textstyle\sum_{i = 1}^{3}} \alpha_{ij}{\bf a}_{i}^{*}], with [\alpha_{ij}] a [3 \times d] matrix σ. In the case of an IMS, at least one entry to σ has to be irrational. The wavelength of the modulation function is [\lambda_{j} = 1/q_{j}]. The set of vectors H forms a Fourier module [M^{*} = \{{\bf H} = {\textstyle\sum_{i = 1}^{3 + d}} h_{i}{\bf a}_{i}^{*} |h_{i} \in {\bb Z}\}] of rank [n = 3 + d], which can be decomposed into a rank 3 and a rank d submodule [M^{*} = M_{1}^{*} \oplus M_{2}^{*}.\; M_{1}^{*} = \{h_{1}{\bf a}_{1}^{*} + h_{2}{\bf a}_{2}^{*} + h_{3}{\bf a}_{3}^{*}\}] corresponds to a [{\bb Z}] module of rank 3 in a 3D subspace (the physical space), [M_{2}^{*} = \{h_{4}{\bf a}_{4}^{*} + \ldots + h_{3 + d}{\bf a}_{3 + d}^{*}\}] corresponds to a [{\bb Z}] module of rank d in a dD subspace (perpendicular space). The submodule [M_{1}] is identical to the 3D reciprocal lattice [\Lambda^{*}] of the average structure. [M_{2}] results from the projection of the perpendicular-space component of the [(3 + d)\hbox{D}] reciprocal lattice [\Sigma^{*}] upon the physical space. Owing to the coincidence of one subspace with the physical space, the dimension of the embedding space is given as [(3 + d)\hbox{D}] and not as nD. This terminology points out the special role of the physical space.

[Figure 4.6.3.1]

Figure 4.6.3.1 | top | pdf |

Schematic diffraction patterns for IMSs with (a) 1D, (b) 2D and (c) 3D modulation. The satellite vectors correspond to [{\bf q} = \alpha_{1}{\bf a}_{1}^{*}] in (a), [{\bf q}_{1} = \alpha_{11}{\bf a}_{1}^{*} + (1/2){\bf a}_{2}^{*}] and [{\bf q}_{2} = -\alpha_{12}{\bf a}_{1}^{*} + (1/2){\bf a}_{2}^{*}], where [\alpha_{11} = \alpha_{12}], in (b), and [{\bf q}_{1} = \alpha_{11}{\bf a}_{1}^{*} + \alpha_{31}{\bf a}_{3}^{*}], [{\bf q}_{2} = \alpha_{12}(-{\bf a}_{1}^{*} + {\bf a}_{2}^{*})\ +] [ \alpha_{32}{\bf a}_{3}^{*}], [{\bf q}_{3} = -\alpha_{13}{\bf a}_{2}^{*} + \alpha_{33}{\bf a}_{3}^{*}], where [\alpha_{11} = \alpha_{12} = \alpha_{13}] and [\alpha_{31} =] [ \alpha_{32} = \alpha_{33}], in (c). The areas of the circles are proportional to the reflection intensities. Main (filled circles) and satellite (open circles) reflections are indexed (after Janner et al., 1983b[link]).

Hence the reciprocal-basis vectors [{\bf a}_{i}^{*}, i = 1, \ldots, 3 + d], can be considered to be physical-space projections of reciprocal-basis vectors [{\bf d}_{i}^{*}, i = 1, \ldots, 3 + d], spanning a [(3 + d)\hbox{D}] reciprocal lattice [\Sigma^{*}]: [\displaylines{\Sigma^{*} = \left\{{\bf H} = {\textstyle\sum\limits_{i = 1}^{3 + d}} h_{i}{\bf d}_{i}^{*}\Big|h_{i} \in {\bb Z}\right\},\cr {\bf d}_{i}^{*} = ({\bf a}_{i}^{*}, {\bf 0}),\, i = 1, \ldots, 3 \hbox{ and } {\bf d}_{3 + j}^{*} = ({\bf a}_{3 + j}^{*}, c{\bf e}_{j}^{*}), \,j = 1, \ldots, d.}]The first vector component of [{\bf d}_{i}^{*}] refers to the physical space, the second to the perpendicular space spanned by the mutually orthogonal unit vectors [{\bf e}_{j}]. c is an arbitrary constant which can be set to 1 without loss of generality.

A direct lattice Σ with basis [{\bf d}_{i}], [i = 1, \ldots, 3 + d] and [{\bf d}_{i} \cdot {\bf d}_{j}^{*} = \delta_{ij}], can be constructed according to [\displaylines{\Sigma = \left\{{\bf r} = {\textstyle\sum\limits_{i = 1}^{3 + d}} m_{i}{\bf d}_{i}\Big| m_{i} \in {\bb Z}\right\},\cr {\bf d}_{i} = \left({\bf a}_{i}, - {\textstyle\sum\limits_{j = 1}^{d}} \alpha_{ij}(1/c){\bf e}_{j}\right),\, i = 1,\ldots, 3\cr \hbox{ and } {\bf d}_{3 + j} = \left({\bf 0},(1/c){\bf e}_{j}^{*}\right),\, j = 1, \ldots, d.}]Consequently, the aperiodic structure in physical space [{\bf V}^{\parallel}] is equivalent to a 3D section of the [(3 + d)\hbox{D}] hypercrystal.

4.6.3.1.1. Indexing

| top | pdf |

The 3D reciprocal space [M^{*}] of a [(3 + d)\hbox{D}] IMS consists of two separable contributions, [M^{*} = \left\{{\bf H} = {\textstyle\sum\limits_{i = 1}^{3}} h_{i} {\bf a}_{i}^{*} + {\textstyle\sum\limits_{j = 1}^{d}} m_{j} {\bf q}_{j}\right\},]the set of main reflections [(m_{j} = 0)] and the set of satellite reflections [(m_{j} \neq 0)] (Fig. 4.6.3.1)[link]. In most cases, the modulation is only a weak perturbation of the crystal structure. The main reflections are related to the average structure, the satellites to the difference between the average and actual structure. Consequently, the satellite reflections are generally much weaker than the main reflections and can be easily identified. Once the set of main reflections has been separated, a conventional basis [{\bf a}_{i}^{*}, i = 1, \ldots, 3], for [\Lambda^{*}] is chosen.

The only ambiguity is in the assignment of rationally independent satellite vectors [{\bf q}_{i}]. They should be chosen inside the reciprocal-space unit cell (Brillouin zone) of [\Lambda^{*}] in such a way as to give a minimal number d of additional dimensions. If satellite vectors reach the Brillouin-zone boundary, centred [(3 + d)\hbox{D}] Bravais lattices are obtained. The star of satellite vectors has to be invariant under the point-symmetry group of the diffraction pattern. There should be no contradiction to a reasonable physical modulation model concerning period or propagation direction of the modulation wave. More detailed information on how to find the optimum basis and the correct setting is given by Janssen et al. (2004)[link] and Janner et al. (1983a[link],b[link]).

4.6.3.1.2. Diffraction symmetry

| top | pdf |

The Laue symmetry group [K^{\rm L} = \{R\}] of the Fourier module [M^{*}], [M^{*} = \left\{{\bf H} = {\textstyle\sum\limits_{i = 1}^{3}} h_{i} {\bf a}_{i}^{*} + {\textstyle\sum\limits_{j = 1}^{d}} m_{j} {\bf q}_{j} = {\textstyle\sum\limits_{i = 1}^{3 + d}} h_{i} {\bf a}_{i}^{*}\right\}, \Lambda^{*} = \left\{{\bf H} = {\textstyle\sum\limits_{i = 1}^{3}} h_{i} {\bf a}_{i}^{*}\right\},]is isomorphous to or a subgroup of one of the 11 3D crystallographic Laue groups leaving [\Lambda^{*}] invariant. The action of the point-group symmetry operators R on the reciprocal basis [{\bf a}_{i}^{*}, i = 1, \ldots, 3 + d], can be written as [R{\bf a}_{i}^{*} = {\textstyle\sum\limits_{j = 1}^{3 + d}} \Gamma_{ij}^{T} (R) {\bf a}_{j}^{*}, i = 1, \ldots, 3 + d.]

The [(3 + d) \times (3 + d)] matrices [\Gamma^{T}(R)] form a finite group of integral matrices which are reducible, since R is already an orthogonal transformation in 3D physical space. Consequently, R can be expressed as pair of orthogonal transformations [(R^{\parallel}, R^{\perp})] in 3D physical and dD perpendicular space, respectively. Owing to their mutual orthogonality, no symmetry relationship exists between the set of main reflections and the set of satellite reflections. [\Gamma^{T}(R)] is the transpose of [\Gamma (R)] which acts on vector components in direct space.

For the [(3 + d)\hbox{D}] direct-space (superspace) symmetry operator [(R_{\rm s}, {\bf t}_{\rm s})] and its matrix representation [\Gamma (R_{\rm s}, {\bf t}_{\rm s})] on Σ, the following decomposition can be performed: [\Gamma (R_{\rm s}) = \pmatrix{\Gamma^{\parallel}(R) &0\cr \Gamma^{\rm M}(R) &\Gamma^{\perp}(R)\cr} \hbox{ and } {\bf t}_{\rm s} = ({\bf t}_{3}, {\bf t}_{d}).][\Gamma^{\parallel}(R)] is a [3 \times 3] matrix, [\Gamma^{\perp}(R)] is a [d \times d] matrix and [\Gamma^{\rm M}(R)] is a [d \times 3] matrix. The translation operator [{\bf t}_{\rm s}] consists of a 3D vector [{\bf t}_{3}] and a dD vector [{\bf t}_{d}]. According to Janner & Janssen (1979)[link], [\Gamma^{\rm M}(R)] can be derived from [\Gamma^{\rm M}(R) = \sigma \Gamma^{\parallel}(R) - \Gamma^{\perp}(R) \sigma]. [\Gamma^{\rm M}(R)] has integer elements only as it contains components of primitive-lattice vectors of [\Lambda^{*}], whereas σ in general consists of a rational and an irrational part: [\sigma = \sigma^{\rm i} + \sigma^{\rm r}]. Thus, only the rational part gives rise to nonzero entries in [\Gamma^{\rm M}(R)]. With the order of the Laue group denoted by N, one obtains [\sigma^{\rm i} \equiv (1/N) {\textstyle\sum_{R}} \Gamma^{\perp}(R) \sigma \Gamma^{\parallel}(R)^{-1}], where [\Gamma^{\perp}(R) \sigma^{\rm i} \Gamma^{\parallel}(R)^{-1} = \sigma^{\rm i}], implying that [\Gamma^{\rm M}(R) =\sigma^{\rm r} \Gamma^{\parallel}(R)] [ - \ \Gamma^{\perp}(R) \sigma^{\rm r}] and [0 = \sigma^{\rm i} \Gamma^{\parallel}(R) - \Gamma^{\perp}(R) \sigma^{\rm i}].

Example

In the case of a 3D IMS with 1D modulation [(d = 1)] the [3 \times d] matrix [\sigma = \pmatrix{\alpha_{1}\cr \alpha_{2}\cr \alpha_{3}\cr}]has the components of the wavevector [{\bf q} = {\textstyle\sum_{i = 1}^{3}} \alpha_{i} {\bf a}_{i}^{*} = {\bf q}^{\rm i} + {\bf q}^{\rm r}]. [\Gamma^{\perp}(R) = \varepsilon = \pm 1] because for [d = 1], q can only be transformed into [\pm {\bf q}]. Corresponding to [{\bf q}^{\rm i} \equiv (1/N) {\textstyle\sum_{R}} \varepsilon R {\bf q}], one obtains [R^{T} {\bf q}^{\rm i} \equiv \varepsilon {\bf q}^{\rm i}] (modulo [\Lambda^{*}]). The [3 \times 1] row matrix [\Gamma^{\rm M}(R)] is equivalent to the difference vector between [R^{T} {\bf q}] and [\varepsilon {\bf q}] (Janssen et al., 2004[link]).

For a monoclinic modulated structure with point group [2/m] for [M^{*}] (unique axis [{\bf a}_{3}]) and satellite vector [{\bf q} = (1/2) {{\bf a}_{1}}^{*} + \alpha_{3} {{\bf a}_{3}}^{*}], with [\alpha_{3}] an irrational number, one obtains [\eqalign{{\bf q}^{\rm i} &\equiv (1/N) {\textstyle\sum\limits_{R}} \varepsilon R {\bf q}\cr &= {\textstyle{1 \over 4}} \left(+1 \cdot \pmatrix{1 &0 &0\cr 0 &1 &0\cr 0 &0 &1\cr} \pmatrix{1/2\cr 0\cr \alpha_{3}\cr}\right.\cr &\quad+ 1 \cdot \pmatrix{\bar{1} &0 &0\cr 0 &\bar{1} &0\cr 0 &0 &1\cr} \pmatrix{1/2\cr 0\cr \alpha_{3}\cr}- 1 \cdot \pmatrix{\bar{1} &0 &0\cr 0 &\bar{1} &0\cr 0 &0 &\bar{1}\cr} \pmatrix{1/2\cr 0\cr \alpha_{3}\cr}\cr &\quad \left.- 1 \cdot \pmatrix{1 &0 &0\cr 0 &1 &0\cr 0 &0 &\bar{1}\cr} \pmatrix{1/2\cr 0\cr \alpha_{3}\cr}\right)\cr &= \pmatrix{0\cr 0\cr \alpha_{3}\cr}.}]From the relations [R^{T} {\bf q}^{\rm i} \equiv \varepsilon {\bf q}^{\rm i}\; (\hbox{modulo } \Lambda^{*})], it can be shown that the symmetry operations 1 and 2 are associated with the perpendicular-space transformations [\varepsilon = 1], and m and [\bar{1}] with [\varepsilon = -1]. The matrix [\Gamma^{\rm M}(R)] is given by [\eqalign{\Gamma^{\rm M}(2) &= \sigma^{\rm r} \Gamma^{\parallel}(2) - \Gamma^{\perp}(2) \sigma^{\rm r}\cr &= \pmatrix{1/2\cr 0\cr 0\cr} \pmatrix{\bar{1} &0 &0\cr 0 &\bar{1} &0\cr 0 &0 &1\cr} - (+1) \pmatrix{1/2\cr 0\cr 0\cr} = \pmatrix{\bar{1}\cr 0\cr 0\cr}}]for the operation 2, for instance.

The matrix representations [\Gamma^{T}(R_{\rm s})] of the symmetry operators R in reciprocal [(3 + d)\hbox{D}] superspace decompose according to [\Gamma^{T}(R_{\rm s}) = \pmatrix{\Gamma^{\|T}(R) &\Gamma^{{\rm M}T}(R)\cr 0 &\Gamma^{\perp T}(R)\cr}.]Phase relationships between modulation functions of symmetry-equivalent atoms can give rise to systematic extinctions of different classes of satellite reflections. The extinction rules may include indices of both main and satellite reflections. A full list of systematic absences is given in the table of [(3 + 1)\hbox{D}] superspace groups (Janssen et al., 2004[link]). Thus, once point symmetry and systematic absences are found, the superspace group can be obtained from the tables in a way analogous to that used for regular 3D crystals. A different approach for the symmetry description of IMSs from the 3D Fourier-space perspective has been given by Dräger & Mermin (1996)[link].

4.6.3.1.3. Structure factor

| top | pdf |

The structure factor of a periodic structure is defined as the Fourier transform of the density distribution [\rho ({\bf r})] of its unit cell (UC): [{F} ({\bf H}) = {\textstyle\int\limits_{\rm UC}} \rho ({\bf r}) \exp (2 \pi i {\bf H} \cdot {\bf r})\ \hbox{d}{\bf r}.]The same is valid in the case of the (3 + d)D description of IMSs. The parallel- and perpendicular-space components are orthogonal to each other and can be separated. The Fourier transform of the parallel-space component of the electron-density distribution of a single atom gives the usual atomic scattering factors [f_{k} ({\bf H}^{\parallel})]. For the structure-factor calculation, one does not need to use [\rho ({\bf r})] explicitly. The hyperatoms correspond to the convolution of the electron-density distribution in 3D physical space with the modulation function in dD perpendicular space. Therefore, the Fourier transform of the (3 + d)D hyperatoms is simply the product of the Fourier transform [f_{k} ({\bf H}^{\parallel})] of the physical-space component with the Fourier transform of the perpendicular-space component, the modulation function.

For a general displacive modulation one obtains for the ith coordinate [x_{ik}] of the kth atom in 3D physical space [x_{ik} = \bar{x}_{ik} + u_{ik} (\bar{x}_{4}, \ldots, \bar{x}_{3 + d}),\ i = 1, \ldots, 3,]where [\bar{x}_{ik}] are the basic-structure coordinates and [u_{ik} (\bar{x}_{4}, \ldots, \bar{x}_{3 + d})] are the modulation functions with unit periods in their arguments (Fig. 4.6.3.2)[link]. The arguments are [\bar{x}_{3 + j} = \alpha_{ij} \bar{x}_{ik}^{0} + t_{j},\ j = 1, \ldots, d], where [\bar{x}_{ik}^{0}] are the coordinates of the kth atom referred to the origin of its unit cell and [t_{j}] are the phases of the modulation functions. The modulation functions [u_{ik} (\bar{x}_{4}, \ldots, \bar{x}_{3 + d})] themselves can be expressed in terms of a Fourier series as [\eqalign{&u_{ik} \left(\bar{x}_{4}, \ldots, \bar{x}_{3 + d}\right)\cr &\quad = {\textstyle\sum\limits_{n_{1} = 1}^{\infty}} \ldots {\textstyle\sum\limits_{n_{d} = 1}^{\infty}} \left\{^{u} C_{ik}^{\,n_{1} \ldots n_{d}}\cos \left[2 \pi \left(n_{1} \bar{x}_{4} + \ldots + n_{d} \bar{x}_{3 + d}\right)\right]\right.\cr &\quad\quad \left. +\, ^{u}S_{ik}^{\,n_{1} \ldots n_{d}}\sin \left[2 \pi \left(n_{1} \bar{x}_{4} + \ldots + n_{d} \bar{x}_{3 + d}\right)\right]\right\},}]where [n_{j}] are the orders of harmonics for the jth modulation wave of the ith component of the kth atom and their amplitudes are [^{u}C_{ik}^{\,n_{1} \ldots n_{d}}] and [^{u}S_{ik}^{\,n_{1} \ldots n_{d}}].

[Figure 4.6.3.2]

Figure 4.6.3.2 | top | pdf |

The relationships between the coordinates [x_{1k}, x_{4k}, \bar{x}_{1}, \bar{x}_{4}] and the modulation function [u_{1k}] in a special section of the [(3 + d)\hbox{D}] space.

Analogous expressions can be derived for a density modulation, i.e., the modulation of the occupation probability [p_{k} \left(\bar{x}_{4}, \ldots, \bar{x}_{3 + d}\right)]: [\eqalign{&p_{k} \left(\bar{x}_{4}, \ldots, \bar{x}_{3 + d}\right)\cr &\quad= {\textstyle\sum\limits_{n_{1} = 1}^{\infty}} \ldots {\textstyle\sum\limits_{n_{d} = 1}^{\infty}} \left\{^{p} C_{k}^{\,n_{1} \ldots n_{d}} \cos \left[2 \pi \left(n_{1} \bar{x}_{4} + \ldots + n_{d} \bar{x}_{3 + d}\right)\right] \right.\cr &\quad \quad\left. + \,\,^{p}S_{k}^{\,n_{1} \ldots n_{d}} \sin \left[2 \pi \left(n_{1} \bar{x}_{4} + \ldots + n_{d} \bar{x}_{3 + d}\right)\right]\right\},}]and for the modulation of the tensor of thermal parameters [B_{ijk} \left(\bar{x}_{4}, \ldots, \bar{x}_{3 + d}\right)]: [\eqalign{&B_{ijk} \left(\bar{x}_{4}, \ldots, \bar{x}_{3 + d}\right)\cr &= {\textstyle\sum\limits_{n_{1} = 1}^{\infty}} \ldots {\textstyle\sum\limits_{n_{d} = 1}^{\infty}} \left\{^{B} C_{ijk}^{\,n_{1} \ldots n_{d}} \cos \left[2 \pi \left(n_{1} \bar{x}_{4} + \ldots + n_{d} \bar{x}_{3 + d}\right)\right]\right.\cr &\quad \left. + \,\,^{B}S_{ijk}^{\,n_{1} \ldots n_{d}} \sin \left[2 \pi \left(n_{1} \bar{x}_{4} + \ldots + n_{d} \bar{x}_{3 + d}\right)\right]\right\}.}]The resulting structure-factor formula is [\eqalign{F({\bf H}) &= \textstyle\sum\limits_{k = 1}^{N'} \textstyle\sum\limits_{(R, \, t)} \textstyle\int\limits_{0}^{1} {\rm d} \bar{x}_{4, \, k} \ldots \textstyle\int\limits_{0}^{1} {\rm d} \bar{x}_{3 + d, \, k}f_{k}({\bf H}^{\parallel})p_{k}\cr &\quad \times \exp \left(- \textstyle\sum\limits_{i, \, j = 1}^{3 + d} h_{i} \left[RB_{ijk} R^{T}\right]h_{j} + 2 \pi i \textstyle\sum\limits_{j = 1}^{3 + d} h_{j} Rx_{jk} + h_{j}t_{j}\right)}]for summing over the set (R, t) of superspace symmetry operations and the set of N′ atoms in the asymmetric unit of the [(3 + d)\hbox{D}] unit cell (Yamamoto, 1982[link]). Different approaches without numerical integration based on analytical expressions including Bessel functions have also been developed. For more information see Paciorek & Chapuis (1994)[link], Petricek, Maly & Cisarova (1991)[link], and references therein.

For illustration, some fundamental IMSs will be discussed briefly (see Korekawa, 1967[link]; Böhm, 1977[link]).

Harmonic density modulation. A harmonic density modulation can result on average from an ordered distribution of vacancies on atomic positions. For an IMS with N atoms per unit cell one obtains for a harmonic modulation of the occupancy factor [p_{k} = (p_{k}^{0}/2) \left\{1 + \cos \left[2 \pi \left(\bar{x}_{4, \, k} + \varphi_{k}\right)\right]\right\},\quad 0 \leq p_{k}^{0} \leq 1,]the structure-factor formula for the mth order satellite [(0 \leq m \leq 1)][\eqalign{F_{0}({\bf H}) &= (1/2) {\textstyle\sum\limits_{k = 1}^{N}}\, f_{k}({\bf H}^{\parallel})T_{k}({\bf H}^{\parallel}) \exp (2 \pi i{\bf H} \cdot {\bf r}_{k}),\cr F_{m}({\bf H}) &= (1/2) {\textstyle\sum\limits_{k = 1}^{N}}\, f_{k}({\bf H}^{\parallel})T_{k}({\bf H}^{\parallel})(p_{k}^{0}/2)^{|m|} \exp \left[2 \pi i \left({\textstyle\sum\limits_{i = 1}^{3}} h_{i}x_{ik} + m \varphi_{k}\right)\right].}]Thus, a linear correspondence exists between the structure-factor magnitudes of the satellite reflections and the amplitude of the density modulation. Furthermore, only first-order satellites exist, since the modulation wave consists only of one term. An important criterion for the existence of a density modulation is that a pair of satellites around the origin of the reciprocal lattice exists (Fig. 4.6.3.3)[link].

[Figure 4.6.3.3]

Figure 4.6.3.3 | top | pdf |

Schematic diffraction patterns for 3D IMSs with (a) 1D harmonic and (b) rectangular density modulation. The modulation direction is parallel to [{\bf a}_{2}]. In (a) only first-order satellites exist; in (b), all odd-order satellites can be present. In (c), the diffraction pattern of a harmonic displacive modulation along [{\bf a}_{1}] with amplitudes parallel to [{\bf a}_{2}^{*}] is depicted. Several reflections are indexed. The areas of the circles are proportional to the reflection intensities.

Symmetric rectangular density modulation. The box-function-like modulated occupancy factor can be expanded into a Fourier series, [\displaylines{ p_{k} = p_{k}^{0} (4/\pi) \left\{{\textstyle\sum\limits_{n = 1}^{\infty}} [(-1)^{n + 1}/(2n - 1)] \cos \left[2\pi (2n - 1) (\bar{x}_{4, \, k} + \varphi_{k})\right]\right\},\cr \hfill 0 \leq p_{k}^{0} \leq 1,\cr}]and the resulting structure factor of the mth order satellite is [\eqalign{F_{0} ({\bf H}) &= (1/2) {\textstyle\sum\limits_{k = 1}^{N}}\, f_{k} ({\bf H}^{\parallel}) T_{k} ({\bf H}^{\parallel}) \exp \left(2\pi i {\textstyle\sum\limits_{i = 1}^{3}} h_{i} x_{ik}\right),\cr F_{m} ({\bf H}) &= (1/\pi m) \sin (m\pi / 2) {\textstyle\sum\limits_{k = 1}^{N}}\, f_{k} ({\bf H}^{\parallel}) T_{k} ({\bf H}^{\parallel}) p_{k}^{0} \cr&\quad\times{\rm exp}\left[2\pi i \left({\textstyle\sum\limits_{i = 1}^{3}} h_{i} x_{ik} + m \varphi_{k}\right)\right].}]According to this formula, only odd-order satellites occur in the diffraction pattern. Their structure-factor magnitudes decrease linearly with the order [|m|] (Fig. 4.6.3.3[link]b).

Harmonic displacive modulation. The displacement of the atomic coordinates is given by the function [x_{ik} = x_{ik}^{0} + A_{ik} \cos \left[2\pi (\bar{x}_{4, \, k} + \varphi_{k})\right],\quad i = 1, \ldots, 3,]and the structure factor by [\eqalign{F_{0} ({\bf H}) &= {\textstyle\sum\limits_{k = 1}^{N}} \,f_{k} ({\bf H}^{\parallel}) T_{k} ({\bf H}^{\parallel}) J_{0} (2\pi {\bf H}^{\parallel} \cdot {\bf A}_{k}) \exp \left(2\pi i {\textstyle\sum\limits_{i = 1}^{3}} h_{i} x_{ik}\right),\cr F_{m} ({\bf H}) &= {\textstyle\sum\limits_{k = 1}^{N}}\, f_{k} ({\bf H}^{\parallel}) T_{k} ({\bf H}^{\parallel}) J_{m} (2\pi {\bf H}^{\parallel} \cdot {\bf A}_{k}) \cr&\quad\times {\rm exp}\left[2\pi i \left({\textstyle\sum\limits_{i = 1}^{3}} h_{i} x_{ik} + m \varphi_{k}\right)\right].}]The structure-factor magnitudes of the mth-order satellite reflections are a function of the mth-order Bessel functions. The arguments of the Bessel functions are proportional to the scalar products of the amplitude and the diffraction vector. Consequently, the intensity of the satellites will vary characteristically as a function of the length of the diffraction vector. Each main reflection is accompanied by an infinite number of satellite reflections (Figs. 4.6.3.3[link]c and 4.6.3.4[link]).

[Figure 4.6.3.4]

Figure 4.6.3.4 | top | pdf |

The relative structure-factor magnitudes of mth-order satellite reflections for a harmonic displacive modulation are proportional to the values of the mth-order Bessel function [J_{m}(x)].

References

Böhm, H. (1977). Eine erweiterte Theorie der Satellitenreflexe und die Bestimmung der modulierten Struktur des Natriumnitrits. Habilitation thesis, University of Munster.
Cummins, H. Z. (1990). Experimental studies of structurally incommensurate crystal phases. Phys. Rep. 185, 211–409.
Dräger, J. & Mermin, N. D. (1996). Superspace groups without the embedding: the link between superspace and Fourier-space crystallography. Phys. Rev. Lett. 76, 1489–1492.
Janner, A. & Janssen, T. (1979). Superspace groups. Physica A, 99, 47–76.
Janner, A., Janssen, T. & de Wolff, P. M. (1983a). Bravais classes for incommensurate crystal phases. Acta Cryst. A39, 658–666.
Janner, A., Janssen, T. & de Wolff, P. M. (1983b). Determination of the Bravais class for a number of incommensurate crystals. Acta Cryst. A39, 671–678.
Janssen, T., Janner, A., Looijenga-Vos, A. & de Wolff, P. M. (2004). Incommensurate and commensurate modulated crystal structures. In International Tables for Crystallography, Vol. C, edited by E. Prince, ch. 9.8. Dordrecht: Kluwer Academic Publishers.
Korekawa, M. (1967). Theorie der Satellitenreflexe. Habilitation thesis, University of Munich.
Paciorek, W. A. & Chapuis, G. (1994). Generalized Bessel functions in incommensurate structure analysis. Acta Cryst. A50, 194–203.
Petricek, V., Maly, K. & Cisarova, I. (1991). The computing system `JANA'. In Methods of Structural Analysis of Modulated Structures and Quasicrystals, edited by J. M. Pérez-Mato, F. J. Zuniga & G. Madriaga, pp. 262–267. Singapore: World Scientific.
Yamamoto, A. (1982). Structure factor of modulated crystal structures. Acta Cryst. A38, 87–92.








































to end of page
to top of page