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

International Tables for Crystallography (2010). Vol. B, ch. 4.6, pp. 590-624   | 1 | 2 |
doi: 10.1107/97809553602060000778

Chapter 4.6. Reciprocal-space images of aperiodic crystals

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

This chapter presents an aid to the characterization of aperiodic crystals based on information from diffraction experiments and gives a survey of aperiodic crystals from the viewpoint of the experimentally accessible reciprocal space. Characteristic features of the diffraction patterns of the different types of aperiodic crystals are shown. A standard way of determining the metrics and finding the optimum nD embedding is described. Structure-factor formulae for general and special cases are given.

4.6.1. Introduction

| top | pdf |

The discovery of materials with icosahedral diffraction symmetry (Shechtman et al., 1984[link]) was the main reason for the reassessment of the definition of crystallinity and for the introduction of the concept of aperiodic crystals. The first aperiodic crystal, i.e. a material with Bragg reflections not located only at reciprocal-lattice nodes, was identified long before (Dehlinger, 1927[link]). In the following decades a wealth of incommensurately modulated phases and composite crystals were discovered. Nevertheless, only a few attempts have been made to develop a crystallography of aperiodic crystals; the most powerful of these was the higher-dimensional approach (see de Wolff, 1974[link], 1977[link]; Janner & Janssen, 1979[link], 1980a[link],b[link]; de Wolff et al., 1981[link]). In fact, incommensurate structures can be easily described using the higher-dimensional approach and also, fully equivalently, in a dual way: as a three-dimensional (3D) combination of one or more periodic basic structures and one or several modulation waves (de Wolff, 1984[link]). However, with the discovery of quasicrystals and their noncrystalline symmetries, the latter approach failed and geometrical crystallography including the higher-dimensional approach received new attention. For more recent reviews of the crystallography of all three types of aperiodic crystals see van Smaalen (1995)[link], of incommensurately modulated structures see Cummins (1990)[link], of quasicrystals see Steurer (1990[link], 1996[link]), of quasicrystals and their crystalline approximants see Goldman & Kelton (1993)[link] and Kelton (1995)[link]. Textbooks on quasicrystals have been written by Janot (1994)[link] and Senechal (1995)[link].

According to the traditional crystallographic definition, an ideal crystal corresponds to an infinite 3D periodic arrangement of identical structure motifs. Its symmetry can be described by one of the 230 3D space groups. Mathematically, a periodic structure can be generated by the convolution of a function representing the structure motif with a lattice function. The structure motif can be given, for instance, by the electron-density distribution [\rho ({\bf r})] of one primitive unit cell of the structure. The lattice function [g({\bf r})] is represented by a set of δ functions at the nodes [{\bf r} = {\textstyle\sum_{i = 1}^{3}} k_{i} {\bf a}_{i}] of a 3D lattice Λ with basis [{\bf a}_{i}], [i = 1, \ldots, 3], and [k_{i} \in {\bb Z}] ([{\bb Z}] is the set of integer numbers). In reciprocal space, this convolution corresponds to the product of the Fourier transform [G({\bf H})] of the lattice function [g({\bf r})] and the Fourier transform [F({\bf H}) = {\textstyle\int_{v}} \rho ({\bf r}) \exp (2\pi i {\bf H} \cdot {\bf r}) \hbox{ d}{\bf r}] of the structure motif [\rho ({\bf r})]. [G({\bf H})] is represented by the reciprocal lattice [\Lambda^{*}] decorated with δ functions on the reciprocal-lattice nodes [{\bf H} = {\textstyle\sum_{i = 1}^{3}} h_{i} {\bf a}_{i}^{*}], with the reciprocal-basis vectors [{\bf a}_{i}^{*}], [i = 1, \ldots, 3], defined by [{\bf a}_{i} \cdot {\bf a}_{j}^{*} = \delta_{ij}] and [h_{i} \in {\bb Z}]. The product [G({\bf H}) \times F({\bf H})] is called the weighted reciprocal lattice; the weights are given by the structure factors [F({\bf H})]. Thus, the characteristic feature of an ideal crystal in direct and reciprocal space is the existence of a lattice. In direct space, this lattice is decorated with identical structure motifs preserving translational and point symmetry in the framework of space-group symmetry. In reciprocal space, only the point symmetry between structure factors is maintained. The Fourier spectrum (or Fourier image, i.e. the Fourier transform) of the electron-density distribution of an ideal crystal consists of a countably infinite set of discrete Bragg peaks with a strictly defined minimum distance.

This crystal definition can be generalized to [n \,\gt\, 3] dimensions. A d-dimensional (dD) ideal aperiodic crystal can be defined as a dD irrational section of an n-dimensional (nD, [n \,\gt\, d]) crystal with nD lattice symmetry. The intersection of the nD hypercrystal with the dD physical space is equivalent to a projection of the weighted nD reciprocal lattice [\Sigma^{*} = \left\{{\bf H} = {\textstyle\sum_{i = 1}^{n}} h_{i} {\bf d}_{i}^{*} |h_{i} \in {\bb Z}\right\}] onto the dD physical space. The resulting set (Fourier module) [M^{*} = \left\{{\bf H}^{\parallel} = {\textstyle\sum_{i = 1}^{n}} h_{i} {\bf a}_{i}^{*} |h_{i} \in {\bb Z}\right\}] is countably dense. Countably dense means that the dense set of Bragg peaks can be mapped one-to-one onto the set of natural numbers. Hence, the Bragg reflections can be indexed with integer indices on an appropriate basis. The Fourier module of the projected reciprocal-lattice vectors [{\bf H}^{\parallel}] has the structure of a [{\bb Z}] module of rank n. A [{\bb Z}] module is a free Abelian group, its rank n is given by the number of free generators (rationally independent vectors). The dimension of a [{\bb Z}] module is that of the vector space spanned by it. The vectors [{\bf a}_{i}^{*}] are the images of the vectors [{\bf d}_{i}^{*}] projected onto the physical space [{\bf V}^{\parallel}]. Thus, by definition, the 3D reciprocal space of an ideal aperiodic crystal consists of a countably dense set of Bragg reflections only. Contrary to an ideal crystal, a minimum distance between Bragg reflections does not exist in an aperiodic one. In summary, it may be stressed that the terms aperiodic and periodic refer to properties of crystal structures in dD space. In nD space, as considered here, lattice symmetry is always present and, therefore, the term crystal is used.

Besides the aperiodic crystals mentioned above, other classes of aperiodic structures with strictly defined construction rules exist (see Axel & Gratias, 1995[link]). Contrary to the kind of aperiodic crystals dealt with in this chapter, the Fourier spectra of aperiodic structures considered in the latter reference are continuous and contain only in a few cases additional sharp Bragg reflections (δ peaks).

Experimentally, the borderline between aperiodic crystals and their periodic approximations (crystalline approximants) is not sharply defined. Finite crystal size, static and dynamic disorder, chemical impurities and defects broaden Bragg peaks and cause diffuse diffraction phenomena. Furthermore, the resolution function of the diffraction equipment is limited.

However, the concept of describing an aperiodic structure as a dD physical-space section of an nD crystal (see Section 4.6.2[link]) is only useful if it significantly simplifies the description of its structural order. Thus, depending on the shape of the atomic surfaces, which gives information on the atomic ordering, incommensurately modulated structures (IMSs, Sections 4.6.2.2[link] and 4.6.3.1[link]), composite structures (CSs, Sections 4.6.2.3[link] and 4.6.3.2[link]), or quasiperiodic structures (QSs, Sections 4.6.2.4[link] and 4.6.3.3[link]) can be obtained from irrational cuts. The atomic surfaces are continuous [(n - d)]-dimensional objects for IMSs and CSs, and discrete [(n - d)]-dimensional objects for QSs. A class of aperiodic crystals with discrete fractal atomic surfaces also exists (Section 4.6.2.5[link]). In this case the Hausdorff dimension (Hausdorff, 1919[link]) of the atomic surface is not an integer number and smaller than [n - d]. The most outstanding characteristic feature of a fractal is its scale invariance: the object appears similar to itself `from near as from far, that is, whatever the scale' (Gouyet, 1996[link]).

To overcome the problems connected with experimental resolution, the translational symmetry of periodic crystals is used as a hard constraint in the course of the determination of their structures. Hence, space-group symmetry is taken for granted and only the local atomic configuration in a unit cell (actually, asymmetric unit) remains to be determined. In reciprocal space, this assumption corresponds to a condensation of Bragg reflections with finite full width at half maximum (FWHM) to δ peaks accurately located at the reciprocal-lattice nodes. Diffuse diffraction phenomena are mostly neglected. This extrapolation to the existence of an ideal crystal is generally out of the question even if samples of very poor quality (high mosaicity, microdomain structure, defects, …) are investigated.

The same practice is convenient for the determination of real aperiodic structures once the type of idealized aperiodic ordering is `known'. Again, the global ordering principle is taken as a hard constraint. For instance, the question of whether a structure is commensurately or incommensurately modulated can only be answered within a given experimental resolution. Experimentally, the ratio of the wavelength of a modulation to the period of the underlying lattice can always be determined as a rational number only. Saying that a structure is incommensurately modulated, with the above ratio being an irrational number, simply means that the experimental results can be better understood, modelled and interpreted assuming an incommensurate modulation. For example, an incommensurate charge-density wave can be moved through an ideal crystal without changing the energy of the crystal. This is not so for a commensurate modulation. In some cases, the modulation period changes with temperature in discrete steps (`devil's staircase'), generating a series of commensurate superstructures (`lock-in structures'); in other cases, a continuous variation can be observed within the experimental resolution. The latter case will be described best by an incommensurately modulated structure.

However, if only the local structure of an aperiodic crystal is of interest, a structure analysis does not take much more experimental effort than for a regular crystal. In contrast, for the analysis of the global structure, i.e. the characterization of the type of its `aperiodicity', diffraction experiments with the highest possible resolution are essential. Some problems connected with the structure analysis of aperiodic crystals are dealt with in Section 4.6.4[link].

To determine the long-range order – whether a real `quasicrystal' is perfectly quasiperiodic, on average quasiperiodic, a crystalline approximant or a nanodomain structure – requires information from experiments that are sensitive to changes of the global structure. Hence, one needs diffraction experiments that allow the accurate determination of the spatial intensity distribution. Consequently, the limiting factors for such experiments are the maximum spatial and intensity resolution of the diffraction and detection equipment, as well as the size and quality of the sample. Nevertheless, the resolution available on state-of-the-art standard synchrotron-beamline equipment is sufficient to test whether the ordering of atoms in an aperiodic crystal reaches the same degree of perfection as found in high-quality silicon. Of course, the higher the sample quality the more necessary it is to account for dynamical diffraction effects such as reflection broadening and displacement. Otherwise, a misinterpretation may bias the global structure modelling.

The following sections present an aid to the characterization of aperiodic crystals based on information from diffraction experiments and give a survey of aperiodic crystals from the viewpoint of the experimentally accessible reciprocal space. Characteristic features of the diffraction patterns of the different types of aperiodic crystals are shown. A standard way of determining the metrics and finding the optimum nD embedding is described. Structure-factor formulae for general and special cases are given.

4.6.2. The n-dimensional description of aperiodic crystals

| top | pdf |

4.6.2.1. Basic concepts

| top | pdf |

An incommensurate modulation of a lattice-periodic structure destroys its translational symmetry in direct and reciprocal space. In the early seventies, a method was suggested by de Wolff (1974)[link] for restoring the lost lattice symmetry by considering the diffraction pattern of an incommensurately modulated structure (IMS) as a projection of an nD reciprocal lattice upon the physical space. n, the dimension of the superspace, is always larger than or equal to d, the dimension of the physical space. This leads to a simple method for the description and characterization of IMSs as well as a variety of new possibilities in their structure analysis. The nD embedding method is well established today and can be applied to all aperiodic crystals with reciprocal-space structure equivalent to a [{\bb Z}] module with finite rank n (Janssen, 1988[link]). The dimension of the embedding space is determined by the rank of the [{\bb Z}] module, i.e. by the number of reciprocal-basis vectors necessary to allow for indexing all Bragg reflections with integer numbers. The point symmetry of the 3D reciprocal space (Fourier spectrum) constrains the point symmetry of the nD reciprocal lattice and restricts the number of possible nD symmetry groups.

In the following sections, the nD descriptions of the four main classes of aperiodic crystals are demonstrated on simple 1D examples of incommensurately modulated phases, composite crystals, quasicrystals and structures with fractally shaped atomic surfaces. The main emphasis is placed on quasicrystals that show scaling symmetry, a new and unusual property in structural crystallography. A detailed discussion of the different types of 3D aperiodic crystals follows in Section 4.6.3[link].

4.6.2.2. 1D incommensurately modulated structures

| top | pdf |

A periodic deviation of atomic parameters from a reference structure (basic structure, BS) is called a modulated structure (MS). In the case of mutual incommensurability of the basic structure and the modulation period, the structure is called incommensurately modulated. Otherwise, it is called commensurately modulated. The modulated atomic parameters may be one or several of the following:

  • (a) coordinates,

  • (b) occupancy factors,

  • (c) thermal displacement parameters,

  • (d) orientation of the magnetic moment.

An incommensurately modulated structure can be described in a dual way by its basic structure [s({\bf r})] and a modulation function [f(t)]. This allows the structure-factor formula to be calculated and a full symmetry characterization employing representation theory to be performed (de Wolff, 1984[link]). A more general method is the nD description: it relates the dD aperiodic incommensurately modulated structure to a periodic structure in nD space. This simplifies the symmetry analysis and structure-factor calculation, and allows more powerful structure-determination techniques.

The nD embedding method is demonstrated in the following 1D example of a displacively modulated structure. A basic structure [s({\bf r}) = s({\bf r} + n{\bf a})], with period a and [n \in {\bb Z}], is modulated by a function [f(t) = f({\bf q \cdot r}) = f(\alpha r) = f [\alpha r + (na/\alpha)]], with the satellite vector [{\bf q} = \alpha {\bf a}^{*}], period [\lambda = 1/q = a/\alpha], and α a rational or irrational number yielding a commensurately or incommensurately modulated structure [s_{m}({\bf r})] (Fig. 4.6.2.1)[link].

[Figure 4.6.2.1]

Figure 4.6.2.1 | top | pdf |

The combination of a basic structure [s({\bf r})], with period a, and a sinusoidal modulation function [f(t)], with amplitude A, period λ and [t = {\bf q}\cdot {\bf r}], gives a modulated structure (MS) [s_{m}({\bf r})]. The MS is aperiodic if a and λ are on incommensurate length scales. The filled circles represent atoms.

If the 1D IMS and its 1D modulation function are properly combined in a 2D parameter space [{\bf V} = ({\bf V}^{\parallel}, {\bf V}^{\perp})], a 2D lattice-periodic structure results (Fig. 4.6.2.2)[link]. The actual atoms are generated by the intersection of the 1D physical (external, parallel) space [{\bf V}^{\parallel}] with the continuous hyperatoms. The hyper­atoms have the shape of the modulation function along the perpendicular (internal, complementary) space [{\bf V}^{\perp}]. They result from a convolution of the physical-space atoms with their modulation functions.

[Figure 4.6.2.2]

Figure 4.6.2.2 | top | pdf |

2D embedding of the sinusoidally modulated structure illustrated in Fig. 4.6.2.1[link]. The correspondence between the actual displacement of an atom in the 1D structure and the modulation function defined in one additional dimension is illustrated in part (a). Adding to each atom its modulation function in this orthogonal dimension (perpendicular space [{\bf V}^{\perp}]) yields a periodic arrangement in 2D space V, part (b). The MS results as a special section of the 2D periodic structure along the parallel space [{\bf V}^{\parallel}]. It is obvious from a comparison of (b) and (c) that the actual MS is independent of the perpendicular-space scale.

A basis [{\bf d}_{1}, {\bf d}_{2}] (D basis) of the 2D hyperlattice [\Sigma = \{{\bf r} = {\textstyle\sum_{i = 1}^{2}} n_{i} {\bf d}_{i} | n_{i} \in {\bb Z}\}] is given by [{\bf d}_{1} = \pmatrix{a\cr -\alpha/c\cr}_{V}, {\bf d}_{2} = \pmatrix{0\cr 1/c\cr}_{V},]where a is the translation period of the BS and c is an arbitrary constant. The components of the basis vectors are given on a 2D orthogonal coordinate system (V basis). The components of the basis vector [{\bf d}_{1}] are simply the parallel-space period a of the BS and α times the perpendicular-space component of the basis vector [{\bf d}_{2}]. The vector [{\bf d}_{2}] is always parallel to the perpendicular space and its length is one period of the modulation function in arbitrary units (this is expressed by the arbitrary factor [1/c]). An atom at position r of the BS is displaced by an amount given by the modulation function [f(t)], with [f(t) = f({\bf q \cdot r})]. Hence, the perpendicular-space variable t has to adopt the value [{\bf q \cdot r} =] [\alpha {\bf a}^{*} \cdot r{\bf a} = \alpha r] for the physical-space variable r. This can be achieved by assigning the slope α to the basis vector [{\bf d}_{1}]. The choice of the parameter c has no influence on the actual MS, i.e. the way in which the 2D structure is cut by the parallel space (Fig. 4.6.2.2[link]c).

The basis of the lattice [\Sigma^{*} = \{{\bf H} = {\textstyle\sum_{i = 1}^{2}} h_{i} {\bf d}_{i}^{*} | h_{i} \in {\bb Z}\}], reciprocal to Σ, can be obtained from the condition [{\bf d}_{i} \cdot {\bf d}_{j}^{*} = \delta_{ij}]: [{\bf d}_{1}^{*} = \pmatrix{a^{*}\cr 0\cr}_{V}, {\bf d}_{2}^{*} = \pmatrix{\alpha a^{*}\cr c\cr}_{V},]with [a^{*} = 1/a]. The metric tensors for the reciprocal and direct 2D lattices for [c = 1] are [G^{*} = \pmatrix{a^{*2} &\alpha a^{*2}\cr \alpha a^{*2} &1 + \alpha^{2} a^{*2}\cr} \quad\hbox{and}\quad G = \pmatrix{a^{2} + \alpha^{2} &-\alpha\cr -\alpha &1\cr}.]The choice of an arbitrary number for c has no influence on the metrics of the physical-space components of the IMS in direct or reciprocal space.

The Fourier transform of the hypercrystal depicted in Fig. 4.6.2.2[link] gives the weighted reciprocal lattice shown in Fig. 4.6.2.3[link]. The 1D diffraction pattern [M^{*} = \{{\bf H}^{\parallel} = {\textstyle\sum_{i = 1}^{2}} h_{i} {\bf a}_{i}^{*} | h_{i} \in {\bb Z}\}] in physical space is obtained by a projection of the weighted 2D reciprocal lattice [\Sigma^{*}] along [{\bf V}^{\perp}] as the Fourier transform of a section in direct space corresponds to a projection in reciprocal space and vice versa: [M^{*} = \{{\bf H}^{\parallel}\}\quad {\lower2pt\hbox{${\buildrel {\hbox{projection } {\displaystyle\pi}^{\parallel} \hbox{ onto } \hbox{\bf V}^{\parallel}} \over {\longleftarrow\!\!\!\!-\!\!\!\!-\!\!\!\!-\!\!\!\!-\!\!\!\!-\!\!\!\!-\!\!\!\!-\!\!\!\!-\!\!\!\!-\!\!\!\!-\!\!\!\!-\!\!\!\!-\!\!\!\!-\!\!\!\!-\!\!\!\!-\!\!\!\!-\!\!\!\!-\!\!\!\!-\!\!\!\!-\!\!\!\!-}}$}} \quad\Sigma^{*} = \{{\bf H} = ({\bf H}^{\parallel}, {\bf H}^{\perp})\}.]Reciprocal-lattice points lying in physical space are referred to as main reflections, all others as satellite reflections. All Bragg reflections can be indexed with integer numbers [h_{1}, h_{2}] in the 2D description [{\bf H} = h_{1} {\bf d}_{1}^{*} + h_{2} {\bf d}_{2}^{*}]. In the physical-space description, the diffraction vector can be written as [{\bf H}^{\parallel} = h{\bf a}^{*} + m{\bf q} =] [{\bf a}^{*} (h_{1} + \alpha h_{2})], with [{\bf q} = \alpha {\bf a}^{*}] for the satellite vector and [m \in {\bb Z}] the order of the satellite reflection. For a detailed discussion of the embedding and symmetry description of IMSs see, for example, Janssen et al. (2004)[link].

[Figure 4.6.2.3]

Figure 4.6.2.3 | top | pdf |

Schematic representation of the 2D reciprocal-space embedding of the 1D sinusoidally modulated structure depicted in Figs. 4.6.2.1[link] and 4.6.2.2[link]. Main reflections are marked by filled circles and satellite reflections by open circles. The sizes of the circles are roughly related to the reflection intensities. The actual 1D diffraction pattern of the 1D MS results from a projection of the 2D reciprocal space onto the parallel space. The correspondence between 2D reciprocal-lattice positions and their projected images is indicated by dashed lines.

A commensurately modulated structure with [\alpha' = m/n] and [\lambda = (n/m) a], [m, n \in {\bb Z}], and with [c = 1], can be generated by shearing the 2D lattice Σ with a shear matrix [S_{m}]: [\displaylines{{\bf d}'_{i} = {\textstyle\sum\limits_{j = 1}^{2}} S_{mij} {\bf d}_{j}, \hbox{ with } S_{m} = \pmatrix{1 &-x\cr 0 &1\cr}_{D} \hbox{ and } x = \alpha' - \alpha,\cr {\bf d}'_{1} = {\bf d}_{1} - x{\bf d}_{2} = \pmatrix{a\cr -\alpha\cr}_{V} - (\alpha' - \alpha) \pmatrix{0\cr 1\cr}_{V} = \pmatrix{a\cr -\alpha'\cr}_{V}, \cr {\bf d}'_{2} = {\bf d}_{2} = \pmatrix{0\cr 1\cr}_{V}.\cr}]The subscript D (V) following the shear matrix indicates that it is acting on the D (V) basis. The shear matrix does not change the distances between the atoms in the basic structure. In reciprocal space, using the inverted and transposed shear matrix, one obtains [\displaylines{{{\bf d}_{i}^{*}}' = {\textstyle\sum\limits_{j = 1}^{2}} (S_{m}^{-1})_{ij}^{T} {\bf d}_{j}^{*}, \hbox{ with } (S_{m}^{-1})^{T} = \pmatrix{1 &0\cr x &1\cr}_{D} \hbox{ and } x = \alpha' - \alpha,\cr {{\bf d}_{1}^{*}}' = {\bf d}_{1}^{*} = \pmatrix{a^{*}\cr 0\cr}_{V},\cr {{\bf d}_{2}^{*}}' = x{\bf d}_{1}^{*} + {\bf d}_{2}^{*} = (\alpha' - \alpha) \pmatrix{a^{*}\cr 0\cr}_{V} + \pmatrix{\alpha a^{*}\cr 1\cr}_{V} = \pmatrix{\alpha' a^{*}\cr 1\cr}_{V}.}]

4.6.2.3. 1D composite structures

| top | pdf |

In the simplest case, a composite structure (CS) consists of two intergrown periodic structures with mutually incommensurate lattices. Owing to mutual interactions, each subsystem may be modulated with the period of the other. Consequently, CSs can be considered as coherent intergrowths of two or more incommensurately modulated substructures. The substructures have at least the origin of their reciprocal lattices in common. However, in all known cases, at least one common reciprocal-lattice plane exists. This means that at least one particular projection of the composite structure exhibits full lattice periodicity.

The unmodulated (basic) 1D subsystems of a 1D incommensurate intergrowth structure can be related to each other in a 2D parameter space [{\bf V} = ({\bf V}^{\parallel}, {\bf V}^{\perp})] (Fig. 4.6.2.4)[link]. The actual atoms result from the intersection of the physical space [{\bf V}^{\parallel}] with the hypercrystal. The hyperatoms correspond to a convolution of the real atoms with infinite lines parallel to the basis vectors [{\bf d}_{1}] and [{\bf d}_{2}] of the 2D hyperlattice [\Sigma = \{{\bf r} = {\textstyle\sum_{i = 1}^{2}} n_{i} {\bf d}_{i} | n_{i} \in {\bb Z}\}].

[Figure 4.6.2.4]

Figure 4.6.2.4 | top | pdf |

2D embedding of a 1D composite structure without mutual interaction of the subsystems. Filled and empty circles represent the atoms of the unmodulated substructures with periods [a_{1}] and [a_{2}], respectively. The atoms result from the parallel-space cut of the linear atomic surfaces parallel to [{\bf d}_{1}] and [{\bf d}_{2}].

An appropriate basis is given by [{\bf d}_{1} = \pmatrix{a_{1}\cr -c\cr}_{V}, {\bf d}_{2} = \pmatrix{0\cr c (a_{2}/a_{1})\cr}_{V},]where [a_{1}] and [a_{2}] are the lattice parameters of the two substructures and c is an arbitrary constant. Taking into account the interactions between the subsystems, each one becomes modulated with the period of the other. Consequently, in the 2D description, the shape of the hyperatoms is determined by their modulation functions (Fig. 4.6.2.5)[link].

[Figure 4.6.2.5]

Figure 4.6.2.5 | top | pdf |

2D embedding of a 1D composite structure with mutual interaction of the subsystems causing modulations. Filled and empty circles represent the modulated substructures with periods [a_{1}] and [a_{2}] of the basic substructures, respectively. The atoms result from the parallel-space cut of the sinusoidal atomic surfaces running parallel to [{\bf d}_{1}] and [{\bf d}_{2}].

A basis of the reciprocal lattice [{\Sigma}^{*} = \{{\bf H} = {\textstyle\sum_{i = 1}^{2}} h_{i} {\bf d}_{i}^{*} | h_{i} \in {\bb Z}\}] can be obtained from the condition [{\bf d}_{i} \cdot {\bf d}_{j}^{*} = \delta_{ij}]: [{\bf d}_{1}^{*} = \pmatrix{a_{1}^{*}\cr 0\cr}_{V}, {\bf d}_{2}^{*} = \pmatrix{a_{2}^{*}\cr (a_{2}^{*}/ca_{1}^{*})\cr}_{V}.]The metric tensors for the reciprocal and the direct 2D lattices for [c = 1] are [G^{*} = \let\normalbaselines\relax\openup2pt\pmatrix{a_{1}^{*2} &a_{1}^{*}a_{2}^{*}\cr a_{1}^{*}a_{2}^{*} &(1 + a_{1}^{*2})(a_{2}^{*}/a_{1}^{*})^{2}\cr} \hbox{ and } G = \pmatrix{1 + a_{1}^{2} &-a_{2}/a_{1}\cr -a_{2}/a_{1} &(a_{2}/a_{1})^{2}\cr}.]The Fourier transforms of the hypercrystals depicted in Figs. 4.6.2.4[link] and 4.6.2.5[link] correspond to the weighted reciprocal lattices illustrated in Figs. 4.6.2.6[link] and 4.6.2.7[link]. The 1D diffraction patterns [M^{*} = \{{\bf H}^{\parallel} = {\textstyle\sum_{i = 1}^{2}} h_{i}{\bf a}_{i}^{*} | h_{i} \in {\bb Z}\}] in physical space are obtained by a projection of the weighted 2D reciprocal lattices [\Sigma^{*}] upon [{\bf V}^{\parallel}]. All Bragg reflections can be indexed with integer numbers [h_{1}, h_{2}] in both the 2D description [{\bf H} = h_{1} {\bf d}_{1}^{*} + h_{2} {\bf d}_{2}^{*}] and in the 1D physical-space description with two parallel basis vectors [{\bf H}^{\parallel} = h_{1} {\bf a}_{1}^{*} + h_{2} {\bf a}_{2}^{*}].

[Figure 4.6.2.6]

Figure 4.6.2.6 | top | pdf |

Schematic representation of the reciprocal space of the embedded 1D composite structure depicted in Fig. 4.6.2.4[link]. Filled and empty circles represent the reflections generated by the substructures with periods [a_{1}] and [a_{2}], respectively. The actual 1D diffraction pattern of the 1D CS results from a projection of the 2D reciprocal space onto the parallel space. The correspondence between 2D reciprocal-lattice positions and their projected images is indicated by dashed lines.

[Figure 4.6.2.7]

Figure 4.6.2.7 | top | pdf |

Schematic representation of the reciprocal space of the embedded 1D composite structure depicted in Fig. 4.6.2.5[link]. Filled and empty circles represent the main reflections of the two subsystems. The satellite reflections generated by the modulated substructures are shown as grey circles. The diameters of the circles are roughly proportional to the intensities of the reflections. The actual 1D diffraction pattern of the 1D CS results from a projection of the 2D reciprocal space onto the parallel space. The correspondence between 2D reciprocal-lattice positions and their projected images is indicated by dashed lines.

The reciprocal-lattice points [{\bf H} = h_{1} {\bf d}_{1}^{*}] and [{\bf H} = h_{2} {\bf d}_{2}^{*}], [h_{1}], [h_{2} \in {\bb Z}], on the main axes [{\bf d}_{1}^{*}] and [{\bf d}_{2}^{*}] are the main reflections of the two substructures. All other reflections are referred to as satellite reflections. Their intensities differ from zero only in the case of modulated substructures. Each reflection of one subsystem coincides with exactly one reflection of the other subsystem.

4.6.2.4. 1D quasiperiodic structures

| top | pdf |

The Fibonacci sequence, the best investigated example of a 1D quasiperiodic structure, can be obtained from the substitution rule σ: [\hbox{S} \rightarrow \hbox{L}], [\hbox{L} \rightarrow \hbox{LS}], replacing the letter S by L and the letter L by the word LS (see e.g. Luck et al., 1993[link]). Applying the substitution matrix [S = \pmatrix{0 &1\cr 1 &1\cr}]associated with σ, this rule can be written in the form [\pmatrix{\hbox{S}\cr \hbox{L}\cr} \rightarrow \pmatrix{0 &1\cr 1 &1\cr} \pmatrix{\hbox{S}\cr \hbox{L}\cr} = \pmatrix{\hbox{L}\cr \hbox{L} + \hbox{S}\cr}.]S gives the sum of letters, [\hbox{L} + \hbox{S} = \hbox{S} + \hbox{L}], and not their order. Consequently, the same substitution matrix can also be applied, for instance, to the substitution [\sigma']: [\hbox{S} \rightarrow \hbox{L}], [\hbox{L} \rightarrow \hbox{SL}]. The repeated action of S on the alphabet [A = \{\hbox{S, L}\}] yields the words [A_{n} = \sigma^{n}(\hbox{S})] and [B_{n} = \sigma^{n}(\hbox{L}) = A_{n + 1}] as illustrated in Table 4.6.2.1[link]. The frequencies of letters contained in the words [A_{n}] and [B_{n}] can be calculated by applying the nth power of the transposed substitution matrix on the unit vector. From [\openup3pt\pmatrix{\nu_{n + 1}^{A}\cr \nu_{n + 1}^{B}\cr} = S^{T} \pmatrix{\nu_{n}^{A}\cr \nu_{n}^{B}\cr}]it follows that [\openup3pt\pmatrix{\nu_{n}^{A}\cr \nu_{n}^{B}\cr} = (S^{T})^{n} \pmatrix{1\cr 1\cr}.]In the case of the Fibonacci sequence, [v_{n}^{B}] gives the total number of letters S and L, and [v_{n}^{A}] the number of letters L.

Table 4.6.2.1| top | pdf |
Expansion of the Fibonacci sequence [B_{n} = \sigma^{n}(\hbox{L})] by repeated action of the substitution rule σ: [\hbox{S} \rightarrow \hbox{L}], [\hbox{L} \rightarrow \hbox{LS}]

[\nu_{\rm L}], [\nu_{\rm S}] are the frequencies of the letters L and S in word [B_{n}].

[B_{n}][\nu_{\rm L}][\nu_{\rm S}]n
L 1 0 0
LS 1 1 1
LSL 2 1 2
LSLLS 3 2 3
LSLLSLSL 5 3 4
LSLLSLSLLSLLS 8 5 5
LSLLSLSLLSLLSLSLLSLSL 13 8 6
  [\vdots] [\vdots] [\vdots]
  [F_{n + 1}] [F_{n}] n

An infinite Fibonacci sequence, i.e. a word [B_{n}] with [n \rightarrow \infty], remains invariant under inflation (deflation). Inflation (deflation) means that the number of letters L, S increases (decreases) under the action of the (inverted) substitution matrix S. Inflation and deflation represent self-similarity (scaling) symmetry operations on the infinite Fibonacci sequence. A more detailed discussion of the scaling properties of the Fibonacci chain in direct and reciprocal space will be given later.

The Fibonacci numbers [F_{n} = F_{n - 1} + F_{n - 2}] form a series with [\lim_{n \rightarrow \infty} (F_{n + 1}/F_{n}) = \tau] {the golden mean [\tau=[1 + (5)^{1/2}]/2] [= 2 \cos(\pi /5) = 1.618 \ldots\]}. The ratio of the frequencies of L and S in the Fibonacci sequence converges to τ if the sequence goes to infinity. The continued fraction expansion of the golden mean τ, [\tau = 1 + {1 \over 1 + {\displaystyle{1 \over 1 + {\displaystyle{1 \over 1 + \ldots}}}}},]contains only the number 1. This means that τ is the `most irrational' number, i.e. the irrational number with the worst truncated continued fraction approximation to it. This might be one of the reasons for the stability of quasiperiodic systems, where τ plays a role. The strong irrationality may impede the lock-in into commensurate systems (rational approximants).

By associating intervals (e.g. atomic distances) with length ratio τ to 1 to the letters L and S, a quasiperiodic structure [s({\bf r})] (Fibonacci chain) can be obtained. The invariance of the ratio of lengths [\hbox{L/S} = (\hbox{L} + \hbox{S})/\hbox{L} = \tau] is responsible for the invariance of the Fibonacci chain under scaling by a factor [\tau^{n}, n \in {\bb Z}]. Owing to a minimum atomic distance S in real crystal structures, the full set of inverse symmetry operators [\tau^{-n}] does not exist. Consequently, the set of scaling operators [s = \{\tau^{0} = 1, \tau^{1}, \ldots \}] forms only a semi-group, i.e. an associative groupoid. Groupoids are the most general algebraic sets satisfying only one of the group axioms: the associative law. The scaling properties of the Fibonacci sequence can be derived from the eigenvalues of the scaling matrix S. For this purpose the equation [\det |S - \lambda I| = 0]with eigenvalue λ and unit matrix I has to be solved. The evaluation of the determinant yields the characteristic polynomial [\lambda^{2} - \lambda - 1 = 0,]yielding in turn the eigenvalues [\lambda_{1} = [1 + (5)^{1/2}]/2 = \tau], [\lambda_{2} = [1 - (5)^{1/2}]/2 = - 1/\tau] and the eigenvectors[\hbox{\bf w}_{1} = \pmatrix{1\cr \tau\cr},\quad \hbox{\bf w}_{2} = \pmatrix{1\cr - 1/\tau\cr}.]Rewriting the eigenvalue equation gives for the first (i.e. the largest) eigenvalue [\pmatrix{0 &1\cr 1 &1\cr} \pmatrix{1\cr \tau\cr} = \pmatrix{\tau\cr 1 + \tau\cr} = \pmatrix{\tau\cr \tau^{2}\cr} = \tau \pmatrix{1\cr \tau\cr}.]Identifying the eigenvector[\pmatrix{1\cr \tau\cr}]with[\pmatrix{\hbox{S}\cr \hbox{L}\cr}]shows that an infinite Fibonacci sequence [s({\bf r})] remains invariant under scaling by a factor τ. This scaling operation maps each new lattice vector [\tau {\bf r}] upon a vector r of the original lattice: [s(\tau {\bf r}) = s({\bf r}).]Considering periodic lattices, these eigenvalues are integer numbers. For quasiperiodic `lattices' (quasilattices) they always correspond to algebraic numbers (Pisot numbers). A Pisot number is the solution of a polynomial equation with integer coefficients. It is larger than one, whereas the modulus of its conjugate is smaller than unity: [\lambda_{1} \,\gt \,1] and [|\lambda_{2}| \,\lt \,1] (Luck et al., 1993[link]). The total lengths [l_{n}^{A}] and [l_{n}^{B}] of the words [A_{n}, B_{n}] can be determined from the equations [l_{n}^{A} = \lambda_{1}^{n} l^{A}] and [l_{n}^{B} = \lambda_{1}^{n} l^{B}] with the eigenvalue [\lambda_{1}]. The left Perron–Frobenius eigenvector [{\bf w}_{1}] of S, i.e. the left eigenvector associated with [\lambda_{1}], determines the ratio S:L to 1:[\tau]. The right Perron–Frobenius eigenvector [{\bf w}_{1}] of S associated with [\lambda_{1}] gives the relative frequencies, 1 and τ, for the letters S and L (for a definition of the Perron–Frobenius theorem see Luck et al., 1993[link], and references therein).

The general case of an alphabet [A = \{\hbox{L}_{1} \ldots \hbox{L}_{k}\}] with k letters (intervals) [\hbox{L}_{i}], of which at least two are on incommensurate length scales and which transform with the substitution matrix S, [{\rm L}'_{i} \rightarrow {\textstyle\sum\limits_{j = 1}^{k}} S_{ij} {\rm L}_{j},]can be treated analogously. S is a [k \times k] matrix with non-negative integer coefficients. Its eigenvalues are solutions of a polynomial equation of rank k with integer coefficients (algebraic or Pisot numbers). The dimension n of the embedding space is generically equal to the number of letters (intervals) k involved by the substitution rule. From substitution rules, infinitely many different 1D quasiperiodic sequences can be generated. However, their atomic surfaces in the nD description are generically of fractal shape (see Section 4.6.2.5[link]).

The quasiperiodic 1D density distribution [\rho ({\bf r})] of the Fibonacci chain can be represented by the Fourier series [\rho ({\bf r}) = (1/V) \,{\textstyle\sum\limits_{{\bf H}^{\parallel}}}\, F ({\bf H}^{\parallel}) \exp (-2\pi i {\bf H}^{\parallel} \cdot {\bf r}),]with [{\bf H}^{\parallel} \in {\bb R}] (the set of real numbers). The Fourier coefficients [F({\bf H}^{\parallel})] form a Fourier module [M^{*} = \{{\bf H}^{\parallel} = {\textstyle\sum_{i = 1}^{2}} h_{i} {\bf a}_{i}^{*} | h_{i} \in {\bb Z}\}] equivalent to a [{\bb Z}] module of rank 2. Thus a periodic function in 2D space can be defined by [\rho ({\bf r}^{\parallel}, {\bf r}^{\perp}) = (1/V) \,{\textstyle\sum\limits_{{\bf H}}}\, F({\bf H}) \exp [-2\pi i ({\bf H}^{\parallel} \cdot {\bf r}^{\parallel} + {\bf H}^{\perp} \cdot {\bf r}^{\perp})],]where [{\bf r} = ({\bf r}^{\parallel}, {\bf r}^{\perp}) \in \Sigma] and [{\bf H} = ({\bf H}^{\parallel}, {\bf H}^{\perp}) \in \Sigma^{*}] are, by construction, direct- and reciprocal-lattice vectors (Figs. 4.6.2.8[link] and 4.6.2.9[link]): [\displaylines{{\bf r} = n_{1} {\bf d}_{1} + n_{2} {\bf d}_{2}, \hbox{ with } {\bf d}_{1} = {1 \over a^{*} (2 + \tau)} \pmatrix{1\cr -\tau\cr}_{V},\cr {\bf d}_{2} = {1 \over a^{*} (2 + \tau)} \pmatrix{\tau\cr 1\cr}_{V}{\semi} \cr {\bf H} = h_{1} {\bf d}_{1}^{*} + h_{2} {\bf d}_{2}^{*}, \hbox{ with } {\bf d}_{1}^{*} = a^{*} \pmatrix{1\cr - \tau\cr}_{V}, {\bf d}_{2}^{*} = a^{*} \pmatrix{\tau\cr 1\cr}_{V}. \cr}]The 1D Fibonacci chain results from the cut of the parallel (physical) space with the 2D lattice Σ decorated with line elements for the atomic surfaces (acceptance domains). In this description, the atomic surfaces correspond simply to the projection of one 2D unit cell upon the perpendicular-space coordinate. This satisfies the condition that each unit cell contributes exactly to one point of the Fibonacci chain (primitive unit cell). The physical space [{\bf V}^{\parallel}] is related to the eigenspace of the substitution matrix S associated with its eigenvalue [\lambda_{1} = \tau]. The perpendicular space [{\bf V}^{\perp}] corresponds to the eigenspace of the substitution matrix S associated with its eigenvalue [\lambda_{2} = - 1/\tau]. Thus, the physical space scales to powers of τ and the perpendicular space to powers of [- 1/\tau].

[Figure 4.6.2.8]

Figure 4.6.2.8 | top | pdf |

2D embedding of the Fibonacci chain. The short and long distances S and L, generated by the intersection of the atomic surfaces with the physical space [{\bf V}^{\parallel}], are indicated. The atomic surfaces are represented by bars parallel to [{\bf V}^{\perp}]. Their lengths correspond to the projection of one unit cell (shaded) upon [{\bf V}^{\perp}].

[Figure 4.6.2.9]

Figure 4.6.2.9 | top | pdf |

Schematic representation of the reciprocal space of the embedded Fibonacci chain depicted in Fig. 4.6.2.8[link]. The physical-space reciprocal basis [{\bf a}_{1}^{*}] and [{\bf a}_{2}^{*}] is marked. The diameters of the filled circles are roughly proportional to the reflection intensities. One 2D reciprocal-lattice unit cell is shadowed. The actual 1D diffraction pattern of the 1D Fibonacci chain results from a projection of the 2D reciprocal space onto the parallel space. The correspondence between 2D reciprocal-lattice positions and their projected images is indicated by dashed lines.

By block-diagonalization, the reducible substitution (scaling) matrix S can be decomposed into two non-equivalent irreducible representations. These can be assigned to the two 1D orthogonal subspaces [{\bf V}^{\parallel}] and [{\bf V}^{\perp}] forming the 2D embedding space [{\bf V} = {\bf V}^{\parallel} \oplus {\bf V}^{\perp}]. Thus, using [WSW^{-1} = S_{V} = S_{V}^{\parallel} \oplus S_{V}^{\perp}], where [W = \pmatrix{1 &\tau\cr -\tau &1\cr} = \pmatrix{{\bf d}_{1}^{*} &{\bf d}_{2}^{*}},]one obtains [\eqalign{\pmatrix{1 &\tau\cr -\tau &1\cr} \pmatrix{0 &1\cr 1 &1\cr}_{D} \pmatrix{1 &\tau\cr -\tau &1\cr}^{-1} &= \pmatrix{\tau{\,\,\vrule height 6pt depth 6pt} &{\hskip -13pt}0\cr\noalign{\vskip -5pt}\noalign{\vskip 5pt}\noalign{\hrule}\cr\noalign{\vskip -5pt}0{\,\,\vrule height 10pt depth 3pt} &{\hskip -8pt}-1/\tau\cr}_{V}\cr &= \pmatrix{S^{\parallel}{\hskip 3.6pt}{\vrule height 6pt depth 6pt} &{\hskip -5.8pt}0\cr\noalign{\hrule} 0{\hskip 7.9pt\vrule height 10pt depth 4pt} &S^{\perp}\cr}_{V},}]the scaling operations [S^{\parallel}] and [S^{\perp}] in parallel and in perpendicular space as indicated by the partition lines.

The metric tensors for the reciprocal and the direct 2D square lattices read [G^{*} = |a^{*}|^{2} (2 + \tau) \pmatrix{1 &0\cr 0 &1\cr} \hbox{ and } G = {1 \over |a^{*}|^{2} (2 + \tau)} \pmatrix{1 &0\cr 0 &1\cr}.]The short distance S of the Fibonacci sequence is related to [a^{*}] by [\eqalign{S &= 1/[a^{*} (2 + \tau)]\cr &= \min \left\{\big|\pi^{\parallel} ({\bf d}_{i} - {\bf d}_{j})\big| \bigg| \big|\pi^{\perp} ({\bf d}_{i} - {\bf d}_{j})\big| \,\lt\, \Omega_{\rm AS} \wedge i, j \in {\bb Z}\right\},\cr}]with the projectors [\pi^{\parallel}] and [\pi^{\perp}] onto [{\bf V}^{\parallel}] and [{\bf V}^{\perp}]. The point density [\rho_{\rm p}] of the Fibonacci chain, i.e. the number of vertices per unit length, can be calculated using the formula [\rho_{\rm p} = {\Omega_{\rm AS} \over \Omega_{\rm UC}} = {(1 + \tau)/[a^{*} (2 + \tau)] \over 1/[|a^{*}|^{2} (2 + \tau)]} = a^{*} \tau^{2},]where [\Omega_{\rm AS}] and [\Omega_{\rm UC}] are the areas of the atomic surface and of the 2D unit cell, respectively.

For an infinite Fibonacci sequence generated from the intervals S and L an average distance d can be calculated: [d = \lim\limits_{n \rightarrow \infty} {F_{n}\hbox{S} + F_{n + 1}\hbox{L} \over F_{n} + F_{n + 1}} = \lim\limits_{n \rightarrow \infty} {F_{n} (\hbox{S} + \tau \hbox{L}) \over F_{n} (1 + \tau)} = {\hbox{S}(1 + \tau^{2}) \over (1 + \tau)} = (3 - \tau)\hbox{S}.]Therefrom, the point density can also be calculated: [\rho_{\rm p} = 1/d = 1/[(3 - \tau)\hbox{S}] = [a^{*} (2 + \tau)]/(3 - \tau) = a^{*} \tau^{2}.]

An approximant structure of the Fibonacci sequence with a unit cell containing m intervals L and n intervals S can be generated by shearing the 2D lattice Σ by the shear matrix [S_{m}], [S_{m} = {1 \over 2 + \tau} \pmatrix{\tau^{2} + x\tau + 1 &-x\cr x\tau^{2} &\tau^{2} -x\tau + 1\cr}_{D},]where [x = (n\tau - m)/(m\tau + n)]: [\eqalign{{\bf d}'_{i} &= {\textstyle\sum\limits_{j = 1}^{2}} S_{mij} {\bf d}_{j}\semi\cr {\bf d}'_{1} &= {1 \over 2 + \tau} \left[\vphantom{1 \over 2 + \tau}(\tau^{2} + x\tau + 1) {\bf d}_{1} - x{\bf d}_{2}\right]\cr &= {1 \over (2 + \tau)a^{*}} \pmatrix{1\cr -\tau -x\cr}_{V}\cr &= {1 \over (2 + \tau)a^{*}} {\let\normalbaselines\relax\openup8pt\pmatrix{1\cr - {\displaystyle{2n\tau + m\tau \over m\tau + n}}\cr}_{V}},\cr {\bf d}'_{2} &= {1 \over 2 + \tau} \left[\vphantom{1 \over 2 + \tau}x\tau^{2} {\bf d}_{1} + (\tau^{2} - x\tau + 1){\bf d}_{2}\right]\cr &= {1 \over (2 + \tau)a^{*}} \pmatrix{\tau\cr -x\tau + 1\cr}_{V}\cr &= {1 \over (2 + \tau)a^{*}} {\let\normalbaselines\relax\openup8pt\pmatrix{\tau\cr {\displaystyle{2m\tau - n\tau \over m\tau + n}}\cr}}_{V}.}]This shear matrix does not change the magnitudes of the intervals L and S. In reciprocal space the inverted and transposed shear matrix is applied on the reciprocal basis, [(S_{m}^{-1})^{T} = {1 \over 2 + \tau} \pmatrix{\tau^{2} -x\tau + 1 &-x\tau^{2}\cr x &\tau^{2} + x\tau + 1\cr}_{D},]where [x = (n\tau - m)/(m\tau + n)]: [\eqalign{{\bf d}_{i}^{*'} &= {\textstyle\sum\limits_{j = 1}^{2}} (S_{m}^{-1})_{ij}^{T} {\bf d}_{j}^{*}\semi\cr {\bf d}_{1}^{*'} &= {1 \over 2 + \tau} \left[\vphantom{1 \over 2 + \tau}(\tau^{2} - x\tau + 1){\bf d}_{1}^{*} - x\tau^{2} {\bf d}_{2}^{*}\right]\cr &= a^{*} \pmatrix{1 - x\tau\cr -\tau\cr}_{V}\cr &= a^{*} {\let\normalbaselines\relax\openup8pt\pmatrix{{\displaystyle{2m\tau - n\tau \over m\tau + n}}\cr -\tau\cr}}_{V},\cr {\bf d}_{2}^{*'} &= {1 \over 2 + \tau} \left[\vphantom{1 \over 2 + \tau}x{\bf d}_{1}^{*} + (\tau^{2} + x\tau + 1){\bf d}_{2}^{*}\right]\cr &= a^{*} \pmatrix{\tau + x\cr 1\cr}_{V}\cr &= a^{*} {\let\normalbaselines\relax\openup8pt\pmatrix{{\displaystyle{2n\tau + m\tau \over m\tau + n}}\cr 1\cr}}_{V}.}]

The point [x_{n}(t)] of the nth interval L or S of an infinite Fibonacci sequence is given by [x_{n}(t) = \{x_{0} + n(3 - \tau) - (\tau - 1) [\hbox{frac} (n\tau + t) - (1/2)]\}\hbox{S},]where t is the phase of the modulation function [y(t) =] [{(\tau - 1)[\hbox{frac} (n\tau + t) - (1/2)]}] (Janssen, 1986[link]). Thus, the Fibonacci sequence can also be dealt with as an incommensurately modulated structure. This is a consequence of the fact that for 1D structures only the crystallographic point symmetries 1 and [\bar{1}] allow the existence of a periodic average structure.

The embedding of the Fibonacci chain as an incommensurately modulated structure can be performed as follows:

  • (1) select a subset [ \Lambda^{*} \subset M^{*}] of strong reflections for main reflections [{\bf H} = h{\bf a}^{*}, h \in {\bb Z}];

  • (2) define a satellite vector [{\bf q} = \alpha {\bf a}^{*}] pointing from each main reflection to the next satellite reflection.

One possible way of indexing based on the same [{\bf a}^{*}] as defined above is illustrated in Fig. 4.6.2.10[link]. The scattering vector is given by [{\bf H}^{\parallel} = h(\tau + 1){\bf a}^{*} + mq], where [{\bf q} = \tau {\bf a}^{*}], or, in the 2D representation, [{\bf H} = h_{1}{\bf d}_{1}^{*} + h_{2}{\bf d}_{2}^{*}], where[{\bf d}_{1}^{*} = a^{*} \pmatrix{1 + \tau\cr 0\cr}_{V}]and[{\bf d}_{2}^{*} = a^{*} \pmatrix{\tau\cr 1\cr}_{V},]with the direct basis [{\bf d}_{1} = {1 \over a^{*} (1 + \tau)} \pmatrix{1\cr - \tau\cr}_{V}, \quad{\bf d}_{2} = {1 \over a^{*}} \pmatrix{0\cr 1\cr}_{V}.]The modulation function is saw-tooth-like (Fig. 4.6.2.11)[link].

[Figure 4.6.2.10]

Figure 4.6.2.10 | top | pdf |

Reciprocal space of the embedded Fibonacci chain as a modulated structure. Several main and satellite reflections are indexed. The square reciprocal lattice of the quasicrystal description illustrated in Fig. 4.6.2.9[link] is indicated by grey lines. The reflections located on [{\bf V}^{\parallel}] can be considered to be projected either from the 2D square lattice of the embedding as for a QS or from the 2D oblique lattice of the embedding as for an IMS.

[Figure 4.6.2.11]

Figure 4.6.2.11 | top | pdf |

2D direct-space embedding of the Fibonacci chain as a modulated structure. The average period is [(3 - \tau)\hbox{S}]. The square lattice in the quasicrystal description shown in Fig. 4.6.2.8[link] is indicated by grey lines. The rod-like atomic surfaces are now inclined relative to [{\bf V}^{\parallel}] and arranged so as to give a saw-tooth modulation wave.

4.6.2.5. 1D structures with fractal atomic surfaces

| top | pdf |

A 1D structure with a fractal atomic surface (Hausdorff dimension 0.9157…) can be derived from the Fibonacci sequence by squaring its substitution matrix S: [\displaylines{ \pmatrix{\hbox{S}\cr \hbox{L}\cr} \rightarrow \pmatrix{1 &1\cr 1 &2\cr} \pmatrix{\hbox{S}\cr \hbox{L}\cr} = \pmatrix{\hbox{S} + \hbox{L}\cr \hbox{S} + 2\hbox{L}\cr}\cr \hbox{with } S^{2} = \pmatrix{1 &1\cr1 &2\cr},\cr}]corresponding to the substitution rule [\hbox{S } \rightarrow \hbox{ SL}], [\hbox{L } \rightarrow \hbox{ LLS}] as well as two other non-equivalent ones (see Janssen, 1995[link]). The eigenvalues [\lambda_{i}] are obtained by calculating [\det |S - \lambda I| = 0.]The evaluation of the determinant gives the characteristic polynomial [\lambda^{2} - 3\lambda + 1 = 0,]with the solutions [\lambda_{1, \, 2} = [3 \pm (5)^{1/2}]/2], with [\lambda_{1} = \tau^{2}] and [\lambda_{2} = 1/\tau^{2} = 2 - \tau], and the same eigenvectors[{\bf w}_{1} = \pmatrix{1\cr \tau\cr},\quad {\bf w}_{2} = \pmatrix{1\cr -1/\tau\cr}]as for the Fibonacci sequence. Rewriting the eigenvalue equation gives [\pmatrix{1 &1\cr 1 &2\cr} \pmatrix{1\cr \tau\cr} = \pmatrix{\tau + 1\cr 2\tau + 1\cr} = \pmatrix{\tau^{2}\cr \tau^{3}\cr} = \tau^{2} \pmatrix{1\cr \tau\cr}.]Identifying the eigenvector[\pmatrix{1\cr \tau\cr}]with[\pmatrix{\hbox{S}\cr \hbox{L}\cr}]shows that the infinite 1D sequence [s({\bf r})] multiplied by powers of its eigenvalue [\tau^{2}] (scaling operation) remains invariant (each new lattice point coincides with one of the original lattice): [s(\tau^{2} {\bf r}) = s({\bf r}).]The fractal sequence can be described on the same reciprocal and direct bases as the Fibonacci sequence. The only difference in the 2D direct-space description is the fractal character of the perpendicular-space component of the hyperatoms (Fig. 4.6.2.12)[link] (see Zobetz, 1993[link]).

[Figure 4.6.2.12]

Figure 4.6.2.12 | top | pdf |

(a) Three steps in the development of the fractal atomic surface of the squared Fibonacci sequence starting from an initiator and a generator. The action of the generator is to cut a piece from each side of the initiator and to add it where the initiator originally ended. This is repeated, cutting thinner and thinner pieces each time from the generated structures. (b) Magnification sequence of the fractal atomic surface illustrating its self-similarity. Each successive figure represents a magnification of a selected portion of the previous figure (from Zobetz, 1993[link]).

4.6.3. Reciprocal-space images

| top | pdf |

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)].

4.6.3.2. Composite structures (CSs)

| top | pdf |

Composite structures consist of N mutually incommensurate substructures with N basic sublattices [\Lambda_{\nu} = \{{\bf a}_{1\nu}, {\bf a}_{2\nu}, {\bf a}_{3\nu}\}], with [\nu = 1, \ldots, N]. The reciprocal sublattices [\Lambda_{\nu}^{*} = \{{\bf a}_{1\nu}^{*}, {\bf a}_{2\nu}^{*}, {\bf a}_{3\nu}^{*}\}], with [\nu = 1, \ldots, N], have either only the origin of the reciprocal lattice or one or two reciprocal-lattice directions in common. Thus, one needs [(3 + d)\, \lt\, 3 N] reciprocal-basis vectors for integer indexing of diffraction patterns that show Bragg reflections at positions given by the Fourier module [M^{*}]. The CSs discovered to date have at least one lattice direction in common and consist of a maximum number of [N = 3] substructures. They can be divided in three main classes: channel structures, columnar packings and layer packings (see van Smaalen, 1992[link], 1995[link]).

In the following, the approach of Janner & Janssen (1980b)[link] and van Smaalen (1992[link], 1995[link], and references therein) for the description of CSs is used. The set of diffraction vectors of a CS, i.e. its Fourier module [M^{*} = \{{\textstyle\sum_{i = 1}^{3 + d}} h_{i}{\bf a}_{i}^{*}\}], can be split into the contributions of the ν subsystems by employing [3 \times (3 + d)] matrices [Z_{ik\nu}] with integer coefficients [{\bf a}_{i\nu}^{*} = {\textstyle\sum_{k = 1}^{3 + d}} Z_{ik\nu}{\bf a}_{k}^{*}], [i = 1, \ldots, 3]. In the general case, each subsystem will be modulated with the periods of the others due to their mutual inter­actions. Thus, in general, CSs consist of several intergrown incommensurately modulated substructures. The satellite vectors [{\bf q}_{j\nu},\ j = 1, \ldots, d], referred to the νth subsystem can be obtained from [M^{*}] by applying the [d \times (3 + d)] integer matrices [V_{jk\nu}]: [{\bf q}_{j\nu} = {\textstyle\sum_{k = 1}^{3 + d}} V_{jk\nu} {\bf a}_{k}^{*}], [j = 1, \ldots , d]. The matrices consisting of the components [\sigma_{\nu}] of the satellite vectors [{\bf q}_{j\nu}] with regard to the reciprocal sublattices [\Lambda_{\nu}^{*}] can be calculated by [\sigma_{\nu} = (V_{3\nu} + V_{d\nu} \sigma)(Z_{3\nu} + Z_{d\nu} \sigma)^{-1}], where the subscript 3 refers to the [3 \times 3] submatrix of physical space and the subscript d to the [d \times d] matrix of the internal space. The juxtaposition of the [3 \times (3 + d)] matrix [Z_{\nu}] and the [d \times (3 + d)] matrix [V_{\nu}] defines the non-singular [(3 + d) \times (3 + d)] matrix [W_{\nu}], [W_{\nu} = \pmatrix{Z_{\nu}\cr V_{\nu}\cr}.]This matrix allows the reinterpretation of the Fourier module [M^{*}] as the Fourier module [M_{\nu}^{*} = M^{*} W_{\nu}] of a d-dimensionally modulated subsystem ν. It also describes the coordinate transformation between the superspace basis Σ and [\Sigma_{\nu}].

The superspace description is obtained analogously to that for IMSs (see Section 4.6.3.1[link]) by considering the 3D Fourier module [M^{*}] of rank [3 + d] as the projection of a [(3 + d)\hbox{D}] reciprocal lattice [\Sigma^{*}] upon the physical space. Thus, one obtains for the definition of the direct and reciprocal [(3 + d)] lattices (Janner & Janssen, 1980b[link]) [\eqalign{\Sigma^{*}: \,&\cases{{\bf a}_{i}^{*} & $= ({\bf a}_{i}^{*}, {\bf 0}){\hbox to 48pt{}} i = 1, \ldots, 3$\cr {\bf a}_{3 + j}^{*} & $= ({\bf a}_{3 + j}^{*}, {\bf e}_{j}^{*}){\hbox to 39pt{}} j = 1, \ldots, d$\cr}\cr \Sigma: \,&\cases{{\bf a}_{i}{\hbox to 3pt{}} &$= ({\bf a}_{i}, - {\sum\limits_{j = 1}^{d}} \sigma_{ji} {\bf e}_{j}){\hbox to 20pt{}} i = 1, \ldots, 3$\cr {\bf a}_{3 + j} &$= ({\bf 0}, {\bf e}_{j}) {\hbox to 51.5pt{}}j = 1, \ldots, d.$\cr}}]

4.6.3.2.1. Indexing

| top | pdf |

The indexing of diffraction patterns of composite structures can be performed in the following way:

  • (1) find the minimum number of reciprocal lattices [\Lambda_{\nu}^{*}] necessary to index the diffraction pattern;

  • (2) find a basis for [M^{*}], the union of sublattices [\Lambda_{\nu}^{*}];

  • (3) find the appropriate superspace embedding.

The [(3 + d)] vectors [{\bf a}_{i}^{*}] forming a basis for the 3D Fourier module [M^{*} = \{{\textstyle\sum_{i = 1}^{3 + d}} h_{i}{\bf a}_{i}^{*}\}] can be chosen such that [{\bf a}_{1}^{*}], [{\bf a}_{2}^{*}] and [{\bf a}_{3}^{*}] are linearly independent. Then the remaining d vectors can be described as a linear combination of the first three, defining the [d \times 3] matrix σ: [{\bf a}_{3 + j}^{*} = {\textstyle\sum_{i = 1}^{3 + d}} \sigma_{ji}{\bf a}_{i}^{*},\; j = 1, \ldots, d]. This is formally equivalent to the reciprocal basis obtained for an IMS (see Section 4.6.3.1[link]) and one can proceed in an analogous way to that for IMSs.

4.6.3.2.2. Diffraction symmetry

| top | pdf |

The symmetry of CSs can be described with basically the same formalism as used for IMSs. This is a consequence of the formally equivalent applicability of the higher-dimensional approach, in particular of the superspace-group theory developed for IMSs [see Janner & Janssen (1980a[link],b[link]); van Smaalen (1991[link], 1992[link]); Yamamoto (1992a[link])].

4.6.3.2.3. Structure factor

| top | pdf |

The structure factor [F({\bf H})] of a composite structure consists of the weighted contributions of the subsystem structure factors [F_{\nu}({\bf H}_{\nu})]:[\eqalign{F({\bf H}) &= \textstyle\sum\limits_{\nu}| J_{\nu} |F_{\nu}({\bf H}_{\nu})\semi\cr F_{\nu}({\bf H}) &= \textstyle\sum\limits_{(R^{\nu}, \, {\bf t}^{\nu})}^{\infty} \textstyle\sum\limits_{k = 1}^{N'} \textstyle\int\limits_{0}^{1} \hbox{d}\bar{x}_{4, \, k}^{\nu}\ldots \textstyle\int\limits_{0}^{1} \hbox{d}\bar{x}_{3 + d, \, k}^{\nu} f_{k}^{\nu} ({\bf H}^{\parallel})p_{k}^{\nu}\cr &\quad \times \exp {\left(- \textstyle\sum\limits_{i, \, j = 1}^{3 + d} h_{i}^{\nu} [R^{\nu} B_{ijk}^{\nu}R^{\nu T}]h_{j}^{\nu}+\,2\pi i \textstyle\sum\limits_{j = 1}^{3 + d} h_{j}^{\nu} R^{\nu} x_{jk}^{\nu} + h_{j}^{\nu} t_{j}^{\nu}\right),}}]with coefficients similar to those for IMSs.

The weights are the Jacobians of the transformations from [{\bf t}_{\nu}] to t, and [{\bf H}_{\nu}] are the reflection indices with respect to the subsystem Fourier modules [M_{\nu}^{*}] (van Smaalen, 1995[link], and references therein). The relative values of [|J_{\nu}|], where [J_{\nu} = \det \left[\left(V_{d\nu} - \sigma_{\nu} \cdot Z_{d\nu}\right)^{-1}\right],] are related to the volume ratios of the contributing subsystems. The subsystem structure factors correspond to those for IMSs (see Section 4.6.3.1[link]). Besides this formula, based on the publications of Yamamoto (1982)[link] and van Smaalen (1995)[link], different structure-factor equations have been discussed (Kato, 1990[link]; Petricek, Maly, Coppens et al., 1991[link]).

4.6.3.3. Quasiperiodic structures (QSs)

| top | pdf |

4.6.3.3.1. 3D structures with 1D quasiperiodic order

| top | pdf |

Structures with quasiperiodic order in one dimension and lattice symmetry in the other two dimensions are the simplest representatives of quasicrystals. A few phases of this structure type have been identified experimentally (see Steurer, 1990[link]). Since the Fibonacci chain represents the most important model of a 1D quasiperiodic structure, it will be used in this section to represent the quasiperiodic direction of 3D structures with 1D quasiperiodic order. As discussed in Section 4.6.2.4[link], 1D quasiperiodic structures are on the borderline between quasiperiodic and incommensurately modulated structures. They can be described using either of the two approaches. In the following, the quasiperiodic description will be preferred to take account of the scaling symmetry.

The electron-density-distribution function [\rho ({\bf r})] of a 1D quasiperiodically ordered 3D crystal can be represented by a 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}^{\parallel} {\bf a}_{i}^{*}] with [{h_{1}}^{\parallel} \in {\bb R}], [{h_{2}}^{\parallel}, {h_{3}}^{\parallel} \in {\bb Z}] or with integer indexing [{\bf H} = {\textstyle\sum_{i = 1}^{4}} h_{i} {\bf a}_{i}^{*}] with [h_{i} \in {\bb Z}]. The set of all vectors H forms a Fourier module [M^{*} = \{{\bf H}^{\parallel} = {\textstyle\sum_{i = 1}^{4}} h_{i}{\bf a}_{i}^{*} | h_{i} \in {\bb Z}\}] of rank 4 which can be decomposed into two rank 2 submodules [M^{*} = M_{1}^{*} \oplus M_{2}^{*}]. [M_{1}^{*} = \{h_{1}{\bf a}_{1}^{*} + h_{2}{\bf a}_{2}^{*}\}] corresponds to a [{\bb Z}] module of rank 2 in a 1D subspace, [M_{2}^{*} = \{h_{3}{\bf a}_{3}^{*} + h_{4}{\bf a}_{4}^{*}\}] corresponds to a [{\bb Z}] module of rank 2 in a 2D subspace. Consequently, the first submodule can be considered as a projection from a 2D reciprocal lattice, [M_{1}^{*} = \pi^{\parallel} (\Sigma^{*})], while the second submodule is of the form of a reciprocal lattice, [M_{2}^{*} = \Lambda^{*}].

Hence, the reciprocal-basis vectors [{\bf a}_{i}^{*}], [i = 1, \ldots, 4], can be considered to be projections of reciprocal-basis vectors [{\bf d}_{i}^{*}], [i = 1, \ldots, 4], spanning a 4D reciprocal lattice, onto the physical space [\Sigma^{*} = \{{\bf H} = {\textstyle\sum_{i = 1}^{4}} h_{i}{\bf d}_{i}^{*} | h_{i} \in {\bb Z}\}], with [{\bf d}_{1}^{*} = a_{1}^{*} \pmatrix{1\cr -\tau\cr 0\cr 0\cr}, \ {\bf d}_{2}^{*} = a_{1}^{*} \pmatrix{\tau\cr 1\cr 0\cr 0\cr}, \ {\bf d}_{3}^{*} = a_{3}^{*} \pmatrix{0\cr 0\cr 1\cr 0\cr}, \ {\bf d}_{4}^{*} = a_{4}^{*} \pmatrix{0\cr 0\cr 0\cr 1\cr}.]A direct lattice Σ with basis [{\bf d}_{i}, i = 1, \ldots, 4] and [{\bf d}_{i} \cdot {\bf d}_{j}^{*} = \delta_{ij}], can be constructed according to (compare Fig. 4.6.2.8)[link] [\Sigma =] [\{{\bf r} = {\textstyle\sum_{i = 1}^{4}} m_{i}{\bf d}_{i} | m_{i} \in {\bb Z}\}], with [\displaylines{{\bf d}_{1} = {1 \over a_{1}^{*} (2 + \tau)} \pmatrix{1\cr -\tau\cr 0\cr 0\cr}, {\bf d}_{2} = {1 \over a_{1}^{*} (2 + \tau)} \pmatrix{\tau\cr 1\cr 0\cr 0\cr}, \cr {\bf d}_{3}^{*} = {1 \over a_{3}^{*}} \pmatrix{0\cr 0\cr 1\cr 0\cr}, \ {\bf d}_{4}^{*} = {1 \over a_{4}^{*}} \pmatrix{0\cr 0\cr 0\cr 1\cr}.}]Consequently, the structure in physical space [{\bf V}^{\parallel}] is equivalent to a 3D section of the 4D hypercrystal.

4.6.3.3.1.1. Indexing

| top | pdf |

The reciprocal space of the Fibonacci chain is densely filled with Bragg reflections (Figs. 4.6.2.9[link] and 4.6.3.5[link]). According to the nD embedding method, the shorter the parallel-space distance [\Delta {\bf H}^{\parallel} = {\bf H}_{2}^{\parallel} - {\bf H}_{1}^{\parallel}] between two Bragg reflections, the larger the corresponding perpendicular-space distance [\Delta {\bf H}^{\perp} = {\bf H}_{2}^{\perp} - {\bf H}_{1}^{\perp}] becomes. Since the structure factor [F({\bf H})] decreases rapidly as a function of [{\bf H}^{\perp}] (Fig. 4.6.3.6)[link], `neighbouring' reflections of strong Bragg peaks are extremely weak and, consequently, the reciprocal space appears to be filled with discrete Bragg peaks even for low-resolution experiments.

[Figure 4.6.3.5]

Figure 4.6.3.5 | top | pdf |

The structure factors [F({\bf H})] (below) and their magnitudes [|F({\bf H})|] (above) of a Fibonacci chain decorated with equal point atoms are shown as a function of the parallel-space component [|{\bf H}^{\parallel}|] of the diffraction vector. The short distance in the Fibonacci chain is [\hbox{S} = 2.5\;\hbox{\AA}], all structure factors within [0 \leq |{\bf H}| \leq 2.5\;\hbox{\AA}^{-1}] have been calculated and normalized to [F(00) = 1].

[Figure 4.6.3.6]

Figure 4.6.3.6 | top | pdf |

The structure factors [F({\bf H})] (below) and their magnitudes [|F({\bf H})|] (above) of a Fibonacci chain decorated with equal point atoms are shown as a function of the perpendicular-space component [|{\bf H}^{\perp}|] of the diffraction vector. The short distance in the Fibonacci chain is [\hbox{S} = 2.5\;\hbox{\AA}], all structure factors within [0 \leq |{\bf H}| \leq 2.5\;\hbox{\AA}^{-1}] have been calculated and normalized to [F(00) = 1].

This property allows an unambiguous identification of a correct set of reciprocal-basis vectors. However, infinitely many sets allowing a correct indexing of the diffraction pattern with integer indices exist. Nevertheless, an optimum basis (low indices are assigned to strong reflections) can be derived: the intensity distribution, not the metrics, characterizes the best choice of indexing. Once the minimum distance S in the structure is identified from chemical considerations, the reciprocal basis should be chosen as described in Section 4.6.2.4[link]. It has to be kept in mind, however, that the identification of the metrics is not sufficient to distinguish in the 1D aperiodic case between an incommensurately modulated structure, a quasiperiodic structure or special kinds of structures with fractally shaped atomic surfaces.

A correct set of reciprocal-basis vectors can be identified in the following way:

  • (1) Find pairs of strong reflections whose physical-space diffraction vectors are related to each other by the factor τ.

  • (2) Index these reflections by assigning an appropriate value to [a^{*}]. This value should be derived from the shortest interatomic distance S expected in the structure.

  • (3) The reciprocal basis is correct if all observable Bragg reflections can be indexed with integer numbers.

4.6.3.3.1.2. Diffraction symmetry

| top | pdf |

The possible Laue symmetry group [K^{\rm 3D}] of the Fourier module [M^{*} = \{{\bf H}^{\parallel} = {\textstyle\sum_{i = 1}^{4}} h_{i}{\bf a}_{i}^{*} | h_{i} \in {\bb Z}\}] is any one of the direct product [K^{\rm 3D} = K^{\rm 2D} \otimes K^{\rm 1D} \otimes \bar{1}]. [K^{\rm 2D}] corresponds to one of the ten crystallographic 2D point groups, [K^{\rm 1D} = \{1\}] in the general case of a quasiperiodic stacking of periodic layers. Consequently, the nine Laue groups [\bar{1}], [2/m], mmm, 4/m, 4/mmm, [\bar{3}], [\bar{3}m], [6/m] and [6/mmm] are possible. These are all 3D crystallographic Laue groups except for the two cubic ones.

The (unweighted) Fourier module shows only 2D lattice symmetry. In the third dimension, the submodule [M_{1}^{*}] remains invariant under the scaling symmetry operation [S^{n} M_{1}^{*} = \tau^{n} M_{1}^{*}] with [n \in {\bb Z}]. The scaling symmetry operators [S^{n}] form an infinite group [s = \{\ldots, S^{-1}, S^{0}, S^{1}, \ldots \}] of reciprocal-basis transformations [S^{n}] in superspace, [\displaylines{S^{n} = \left(\matrix{0 &1 &0 &0\cr 1 &1 &0 &0\cr 0 &0 &1 &0\cr 0 &0 &0 &1\cr}\right)_{D}^{n}, \ S^{-1} = \left(\matrix{\bar{1} &1 &0 &0\cr 1 &0 &0 &0\cr 0 &0 &1 &0\cr 0 &0 &0 &1\cr}\right)_{D}, \cr S^{0} = \left(\matrix{1 &0 &0 &0\cr 0 &1 &0 &0\cr 0 &0 &1 &0\cr 0 &0 &0 &1\cr}\right)_{D},}]and act on the reciprocal basis [{\bf d}_{i}^{*}] in superspace.

4.6.3.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 nD description of a quasiperiodic structure. The parallel- and perpendicular-space components are orthogonal to each other and can be separated. In the case of the 1D Fibonacci sequence, the Fourier transform of the parallel-space component of the electron-density distribution of a single atom gives the usual atomic scattering factor [f({\bf H}^{\parallel})]. Parallel to [x^{\perp}], [\rho({\bf r})] adopts values [\neq 0] only within the interval [- (1 + \tau)/] [[2a^{*}(2 + \tau)] \leq {x}^{\perp} \leq (1 + \tau)/[2a^{*}(2 + \tau)]] and one obtains [\eqalign{ F({\bf H}) &= f({\bf H}^{\parallel})[a^{*} (2 + \tau)]/(1 + \tau) \cr&\quad\times\textstyle\int\limits_{-(1 + \tau)/[2a^{*} (2 + \tau)]}^{+(1 + \tau)/[2a^{*} (2 + \tau)]} \exp(2\pi i{\bf H}^{\perp} \cdot x^{\perp})\,\hbox{d}x^{\perp}.}]The factor [[a^{*}(2 + \tau)]/(1 + \tau)] results from the normalization of the structure factors to [F({\bf 0}) = f(0)]. With [\eqalign{{\bf H} &= h_{1}{\bf d}_{1}^{*} + h_{2}{\bf d}_{2}^{*} + h_{3}{\bf d}_{3}^{*} + h_{4}{\bf d}_{4}^{*}\cr &= h_{1}a_{1}^{*} \pmatrix{1\cr -\tau\cr 0\cr 0\cr} + h_{2}a_{1}^{*} \pmatrix{\tau\cr 1\cr 0\cr 0\cr} + h_{3}a_{3}^{*} \pmatrix{0\cr 0\cr 1\cr 0\cr} + h_{4}a_{4}^{*} \pmatrix{0\cr 0\cr 0\cr 1\cr}}]and [{\bf H}^{\perp} = a_{1}^{*} (-\tau h_{1} + h_{2})] the integrand can be rewritten as [\eqalign{ F({\bf H}) &= f({\bf H}^{\parallel}) [a^{*}(2 + \tau)]/(1 + \tau)\cr&\quad\times \textstyle{\int\limits_{- (1 + \tau)/[2a^{*}(2 + \tau)]}^{+ (1 + \tau)/[2a^{*}(2 + \tau)]}} \exp [2\pi i (- \tau h_{1} + h_{2}) x^{\perp}] \,\hbox{d}x^{\perp},}]yielding [\eqalign{F({\bf H}) &= f({\bf H}^{\parallel}) (2 + \tau)/[2\pi i (-\tau h_{1} + h_{2})(1 + \tau)]\cr&\quad \times\exp [2\pi i(- \tau h_{1} + h_{2}) x^{\perp}] \big|_{- (1 + \tau)/[2a^{*}(2 + \tau)]}^{+ (1 + \tau)/[2a^{*}(2 + \tau)]}.\cr}]Using [\sin x = (\hbox{e}^{ix} - \hbox{e}^{- ix})/2i] gives [\eqalign{F({\bf H}) &= f({\bf H}^{\parallel}) (2 + \tau)/[\pi(- \tau h_{1} + h_{2})(1 + \tau)] \cr&\quad\times\sin[\pi (1 + \tau)(- \tau h_{1} + h_{2})]/(2 + \tau).}]Thus, the structure factor has the form of the function [\sin (x)/x] with x a perpendicular reciprocal-space coordinate. The upper and lower limiting curves of this function are given by the hyperbolae [\pm 1/x] (Fig. 4.6.3.6)[link]. The continuous shape of [F({\bf H})] as a function of [{\bf H}^{\perp}] allows the estimation of an overall temperature factor and atomic scattering factor for reflection-data normalization (compare Figs. 4.6.3.6[link] and 4.6.3.7[link]).

[Figure 4.6.3.7]

Figure 4.6.3.7 | top | pdf |

The structure factors [F({\bf H})] of the Fibonacci chain decorated with aluminium atoms (Uoverall = 0.005 Å2) as a function of the parallel (above) and the perpendicular (below) component of the diffraction vector. The short distance is [\hbox{S} = 2.5\;\hbox{\AA}], all structure factors within [0 \leq |{\bf H}| \leq 2.5\;\hbox{\AA}^{-1}] have been calculated and normalized to [F(00) = 1].

In the case of a 3D crystal structure which is quasiperiodic in one direction, the structure factor can be written in the form [F({\bf H}) = {\textstyle\sum\limits_{k = 1}^{n}} \left[T_{k}({\bf H})f_{k}({\bf H}^{\parallel})g_{k}({\bf H}^{\perp}) \exp (2\pi i {\bf H} \cdot {\bf r}_{k})\right].]The sum runs over all n averaged hyperatoms in the 4D unit cell of the structure. The geometric form factor [g_{k}({\bf H}^{\perp})] corresponds to the Fourier transform of the kth atomic surface, [g_{k} ({\bf H}^{\perp}) = (1/A_{\rm UC}^{\perp}) {\textstyle\int\limits_{A_{k}}} \exp (2\pi i {\bf H}^{\perp} \cdot {\bf r}^{\perp})\ \hbox{d}{\bf r}^{\perp},]normalized to [A_{\rm UC}^{\perp}], the area of the 2D unit cell projected upon [{\bf V}^{\perp}], and [A_{k}], the area of the kth atomic surface.

The atomic temperature factor [T_{k}({\bf H})] can also have perpendicular-space components. Assuming only harmonic (static or dynamic) displacements in parallel and perpendicular space one obtains, in analogy to the usual expression (Willis & Pryor, 1975[link]), [\eqalign{T_{k}({\bf H}) &= T_{k}({\bf H}^{\parallel},{\bf H}^{\perp})\cr& = \exp (-2\pi^{2}{\bf H}^{\| T} \langle {\bf u}_{i}^{\parallel} {\bf u}_{j}^{\| T}\rangle {\bf H}^{\parallel}) \exp (-2\pi^{2}{\bf H}^{\perp T} \langle {\bf u}_{i}^{\perp} {\bf u}_{j}^{\perp T}\rangle {\bf H}^{\perp}),\cr}]with [\displaylines{ \langle {\bf u}_{i}^{\parallel} {\bf u}_{j}^{\| T}\rangle = \pmatrix{\langle {\bf u}_{1}^{\parallel 2}\rangle &\langle {\bf u}_{1}^{\parallel} \cdot {\bf u}_{2}^{\| T}\rangle &\langle {\bf u}_{1}^{\parallel} \cdot {\bf u}_{3}^{\| T}\rangle\cr \langle {\bf u}_{2}^{\parallel} \cdot {\bf u}_{1}^{\| T}\rangle &\langle {\bf u}_{2}^{\parallel 2}\rangle &\langle {\bf u}_{2}^{\parallel} \cdot {\bf u}_{3}^{\| T}\rangle\cr \langle {\bf u}_{3}^{\parallel} \cdot {\bf u}_{1}^{\| T}\rangle &\langle {\bf u}_{3}^{\parallel} \cdot {\bf u}_{2}^{\| T}\rangle &\langle {\bf u}_{3}^{\parallel 2}\rangle\cr}\cr\cr \hbox{ and } \langle {\bf u}_{i}^{\perp} {\bf u}_{j}^{\perp T}\rangle = \langle {\bf u}_{4}^{\perp}\rangle.}]The elements of the type [\langle {\bf u}_{i} \cdot {\bf u}_{j}^{T}\rangle] represent the average values of the atomic displacements along the ith axis times the displacement along the jth axis on the V basis.

4.6.3.3.1.4. Intensity statistics

| top | pdf |

In the following, only the properties of the quasiperiodic component of the 3D structure, namely the Fourier module [M_{1}^{*}], are discussed. The intensities [I({\bf H})] of the Fibonacci chain decorated with point atoms are only a function of the perpendicular-space component of the diffraction vector. [|F({\bf H})|] and [F({\bf H})] are illustrated in Figs. 4.6.3.5[link] and 4.6.3.6[link] as a function of [{\bf H}^{\parallel}] and of [{\bf H}^{\perp}]. The distribution of [|F({\bf H})|] as a function of their frequencies clearly resembles a centric distribution, as can be expected from the centrosymmetric 2D subunit cell. The shape of the distribution function depends on the radius [H_{\max}] of the limiting sphere in reciprocal space. The number of weak reflections increases with the square of [H_{\max}], that of strong reflections only linearly (strong reflections always have small [{\bf H}^{\perp}] components).

The weighted reciprocal space of the Fibonacci sequence contains an infinite number of Bragg reflections within a limited region of the physical space. Contrary to the diffraction pattern of a periodic structure consisting of point atoms on the lattice nodes, the Bragg reflections show intensities depending on the perpendicular-space components of their diffraction vectors.

The reciprocal space of a sequence generated from hyperatoms with fractally shaped atomic surfaces (squared Fibonacci sequence) is very similar to that of the Fibonacci sequence (Figs. 4.6.3.8[link] and 4.6.3.9[link]). However, there are significantly more weak reflections in the diffraction pattern of the `fractal' sequence, caused by the geometric form factor.

[Figure 4.6.3.8]

Figure 4.6.3.8 | top | pdf |

The structure factors [F({\bf H})] (below) and their magnitudes [|F({\bf H})|] (above) of the squared Fibonacci chain decorated with equal point atoms are shown as a function of the parallel-space component [|{\bf H}^{\parallel}|] of the diffraction vector. The short distance is [\hbox{S} = 2.5\;\hbox{\AA}], all structure factors within [0 \leq |{\bf H}| \leq 2.5\;\hbox{\AA}^{-1}] have been calculated and normalized to [F(00) = 1].

[Figure 4.6.3.9]

Figure 4.6.3.9 | top | pdf |

The structure factors [F({\bf H})] (below) and their magnitudes [|F({\bf H})|] (above) of the squared Fibonacci chain decorated with equal point atoms are shown as a function of the perpendicular-space component [|{\bf H}^{\perp}|] of the diffraction vector. The short distance is [\hbox{S} = 2.5\;\hbox{\AA}], all structure factors within [0 \leq |{\bf H}| \leq 2.5\;\hbox{\AA}^{-1}] have been calculated and normalized to [F(00) = 1].

4.6.3.3.1.5. Relationships between structure factors at symmetry-related points of the Fourier image

| top | pdf |

The two possible point-symmetry groups in the 1D quasi­periodic case, [K^{\rm 1D} = 1] and [K^{\rm 1D} = \bar{1}], relate the structure factors to [\eqalign{1:\qquad \quad F ({\bf H}) &= -F (\bar{\bf H}),\cr \bar{1}:\qquad \quad F ({\bf H}) &= F (\bar{\bf H}).}]A 3D structure with 1D quasiperiodicity results from the stacking of atomic layers with distances following a quasiperiodic sequence. The point groups [K^{\rm 3D}] describing the symmetry of such structures result from the direct product [K^{\rm 3D} = K^{\rm 2D} \otimes K^{\rm 1D}.] [K^{\rm 2D}] corresponds to one of the ten crystallographic 2D point groups, [K^{\rm 1D}] can be [\{1\}] or [\{1, m\}]. Consequently, 18 3D point groups are possible.

Since 1D quasiperiodic sequences can be described generically as incommensurately modulated structures, their possible point and space groups are equivalent to a subset of the [(3 + 1)\hbox{D}] superspace groups for IMSs with satellite vectors of the type [(00\gamma)], i.e. [{\bf q} = \gamma {\bf c}^{*}], for the quasiperiodic direction [001] (Janssen et al., 2004[link]).

From the scaling properties of the Fibonacci sequence, some relationships between structure factors can be derived. Scaling the physical-space structure by a factor [\tau^{n}], [n \in {\bb Z}], corresponds to a scaling of the perpendicular space by the inverse factor [(-\tau)^{-n}]. For the scaling of the corresponding reciprocal subspaces, the inverse factors compared to the direct spaces have to be applied.

The set of vectors r, defining the vertices of a Fibonacci sequence [s({\bf r})], multiplied by a factor τ coincides with a subset of the vectors defining the vertices of the original sequence (Fig. 4.6.3.10)[link]. The residual vertices correspond to a particular decoration of the scaled sequence, i.e. the sequence [\tau^{2}s({\bf r})]. The Fourier transform of the sequence [s({\bf r})] then can be written as the sum of the Fourier transforms of the sequences [\tau s({\bf r})] and [\tau^{2}s({\bf r})]; [{\textstyle\sum\limits_{k}} \exp (2\pi i{\bf H}\cdot {\bf r}_{k}) = {\textstyle\sum\limits_{k}} \exp (2\pi i {\bf H}\tau {\bf r}_{k}) + {\textstyle\sum\limits_{k}} \exp [2\pi i{\bf H} (\tau^{2}{\bf r}_{k} + \tau)].]In terms of structure factors, this can be reformulated as [F({\bf H}) = F(\tau {\bf H}) + \exp (2\pi i\tau {\bf H}) F(\tau^{2}{\bf H}).]

[Figure 4.6.3.10]

Figure 4.6.3.10 | top | pdf |

Part … LSLLSLSL … of a Fibonacci sequence [s({\bf r})] before and after scaling by the factor τ. L is mapped onto [\tau \hbox{L}], S onto [\tau \hbox{S} = \hbox{L}]. The vertices of the new sequence are a subset of those of the original sequence (the correspondence is indicated by dashed lines). The residual vertices [\tau^{2}s({\bf r})], which give when decorating [\tau s({\bf r})] the Fibonacci sequence [s({\bf r})], form a Fibonacci sequence scaled by a factor [\tau^{2}].

Hence, phases of structure factors that are related by scaling symmetry can be determined from each other.

Further scaling relationships in reciprocal space exist: scaling a diffraction vector[{\bf H} = h_{1}{\bf d}_{1}^{*} + h_{2}{\bf d}_{2}^{*} = h_{1}a^{*} \pmatrix{1\cr -\tau\cr}_{V} +\ h_{2}a^{*} \pmatrix{\tau\cr 1\cr}_{V}]with the matrix[S = \pmatrix{0 &1\cr 1 &1\cr}_{D},][\eqalign{\pmatrix{0 &1\cr 1 &1\cr}_{D} \pmatrix{h_{1}\cr h_{2}\cr}_{D} &= \pmatrix{F_{n} &F_{n + 1}\cr F_{n + 1} &F_{n + 2}\cr}_{D} \pmatrix{h_{1}\cr h_{2}\cr}_{D} \cr&= \pmatrix{F_{n}h_{1} + F_{n + 1} h_{2}\cr F_{n + 1}h_{1} + F_{n + 2}h_{2}\cr}_{D},}]increases the magnitudes of structure factors assigned to this particular diffraction vector H,[\big|F(S^{n}{\bf H})\big| \,\gt\, \big|F(S^{n - 1}{\bf H}) \big| \,\gt\, \ldots \,\gt\, \big| F(S{\bf H})\big| \,\gt\, \big| F({\bf H})\big|.]

This is due to the shrinking of the perpendicular-space component of the diffraction vector by powers of [(-\tau)^{-n}] while expanding the parallel-space component by [\tau^{n}] according to the eigenvalues τ and [-\tau^{-1}] of S acting in the two eigenspaces [{\bf V}^{\parallel}] and [{\bf V}^{\perp}]:[\displaylines{\eqalign{\pi^{\parallel}(S{\bf H}) &= \left(h_{2} + \tau\left(h_{1} + h_{2}\right)\right) a^{*} = \left(\tau h_{1} + h_{2}\left(\tau + 1 \right)\right) a^{*}\cr&= \tau \left(h_{1} + \tau h_{2}\right) a^{*},\cr \pi^{\perp}(S{\bf H}) &= \left(- \tau h_{2} + h_{1} + h_{2}\right) a^{*} = \left(h_{1} - h_{2}\left(\tau - 1\right)\right) a^{*} \cr&= - (1/\tau) \left(-\tau h_{1} + h_{2}\right) a^{*},\cr}\cr \big| F(\tau^{n}{\bf H}^{\parallel}) \big|\, \gt\, \big| F(\tau^{n - 1}{\bf H}^{\parallel}) \big| \,\gt\, \ldots\, \gt\, \big| F(\tau{\bf H}^{\parallel}) \big| \,\gt\, \big| F({\bf H}^{\parallel}) \big|.\cr}]Thus, for scaling n times we obtain[\eqalign{\pi^{\perp}(S^{n}{\bf H}) &= \left( - \tau \left(F_{n}h_{1} + F_{n + 1}h_{2}\right) + \left(F_{n + 1}h_{1} + F_{n + 2}h_{2}\right)\right) a^{*}\cr &= \left(h_{1} \left(- \tau F_{n} + F_{n + 1}\right) + h_{2} \left(- \tau F_{n + 1} + F_{n + 2}\right)\right) a^{*}}]with[\lim\limits_{n \rightarrow \infty} \left(- \tau F_{n} + F_{n + 1}\right) = 0 \hbox{ and } \lim\limits_{n \rightarrow \infty} \left(- \tau F_{n + 1} + F_{n + 2}\right) = 0,]yielding eventually [\lim\limits_{n \rightarrow \infty} \left(\pi^{\perp}(S^{n}{\bf H})\right) = 0 \hbox{ and } \lim\limits_{n \rightarrow \infty} \left(F(S^{n}{\bf H})\right) = F({\bf 0}).]The scaling of the diffraction vectors H by [S^{n}] corresponds to a hyperbolic rotation (Janner, 1992[link]) with angle [n\varphi], where [\sinh \varphi = 1/2] (Fig. 4.6.3.11)[link]: [\eqalign{\left(\matrix{0 &1\cr 1 &1\cr}\right)^{2n} &= \left(\matrix{\cosh 2n\varphi &\sinh 2n\varphi\cr \sinh 2n\varphi &\cosh 2n\varphi\cr}\right),\cr \left(\matrix{0 &1\cr 1 &1\cr}\right)^{2n + 1} &= \left(\matrix{\sinh \left[\left(2n + 1\right)\varphi\right] &\cosh \left[\left(2n + 1\right)\varphi\right]\cr \cosh \left[\left(2n + 1\right)\varphi\right] &\sinh \left[\left(2n + 1\right)\varphi\right]\cr}\right).}]

[Figure 4.6.3.11]

Figure 4.6.3.11 | top | pdf |

Scaling operations of the Fibonacci sequence. The scaling operation S acts six times on the diffraction vector [{\bf H} = (4\bar{2})] yielding the sequence [(4\bar{2}) \rightarrow (\bar{2}2) \rightarrow (20) \rightarrow (02) \rightarrow (22) \rightarrow (24) \rightarrow (46)].

4.6.3.3.2. Decagonal phases

| top | pdf |

A structure quasiperiodic in two dimensions, periodic in the third dimension and with decagonal diffraction symmetry is called a decagonal phase. Its holohedral Laue symmetry group is [K = 10/mmm]. All reciprocal-space vectors [{\bf H} \in M^{*}] can be represented on a basis (V basis) [{\bf a}_{i}^{*} = a_{i}^{*}\;(\cos 2\pi i/5, \sin 2\pi i/5, 0),] [i = 1, \ldots, 4] and [{\bf a}_{5}^{*} = a_{5}^{*} (0,0,1)] (Fig. 4.6.3.12)[link] as [{\bf H} = {\textstyle\sum_{i = 1}^{5}} h_{i}{\bf a}_{i}^{*}]. The vector components refer to a Cartesian coordinate system in physical (parallel) space. Thus, from the number of independent reciprocal-basis vectors necessary to index the Bragg reflections with integer numbers, the dimension of the embedding space has to be at least five. This can also be shown in a different way (Hermann, 1949[link]).

[Figure 4.6.3.12]

Figure 4.6.3.12 | top | pdf |

Reciprocal basis of the decagonal phase in the 5D description projected upon [{\bf V}^{\parallel}] (above left) and [{\bf V}^{\perp}] (above right). Below, a perspective physical-space view is shown.

The set [M^{*}] of all vectors H remains invariant under the action of the symmetry operators of the point group [10/mmm]. The symmetry-adapted matrix representations for the point-group generators, the tenfold rotation [\alpha = 10], the reflection plane [\beta = m_{2}] (normal of the reflection plane along the vectors [{\bf a}_{i}^{*}][{\bf a}_{i + 3}^{*}] with [i = 1, \ldots, 4] modulo 5) and the inversion operation [\Gamma (\gamma) = \bar{1}] may be written in the form [\openup-4pt\displaylines{\textstyle\Gamma (\alpha) = \pmatrix{0 &1 &\bar{1} &0 &0\cr 0 &1 &0 &\bar{1} &0\cr 0 &1 &0 &0 &0\cr \bar{1} &1 &0 &0 &0\cr 0 &0 &0 &0 &1\cr}_{D} \Gamma (\beta) = \pmatrix{0 &0 &0 &1 &0\cr 0 &0 &1 &0 &0\cr 0 &1 &0 &0 &0\cr 1 &0 &0 &0 &0\cr 0 &0 &0 &0 &1\cr}_{D}\cr\noalign{\vskip 6pt}\cr \Gamma (\gamma) = \pmatrix{\bar{1} &0 &0 &0 &0\cr 0 &\bar{1} &0 &0 &0\cr 0 &0 &\bar{1} &0 &0\cr 0 &0 &0 &\bar{1} &0\cr 0 &0 &0 &0 &\bar{1}\cr}_{D}.}]

By block-diagonalization, these reducible symmetry matrices can be decomposed into non-equivalent irreducible representations. These can be assigned to the two orthogonal subspaces forming the 5D embedding space [{\bf V} = {\bf V}^{\parallel} \oplus {\bf V}^{\perp}], the 3D parallel (physical) subspace [{\bf V}^{\parallel}] and the perpendicular 2D subspace [{\bf V}^{\perp}]. Thus, using [W\Gamma W^{-1} = \Gamma_{V} = \Gamma_{V}^{\parallel} \oplus \Gamma_{V}^{\perp}], we obtain[\eqalign{\Gamma_{V} (\alpha) &= \pmatrix{\cos (\pi/5) &-\sin (\pi/5) &0{\,\vrule height 8pt depth 6pt} &{\hskip -6pt}0 &0\cr\noalign{\vskip -4pt} \sin (\pi/5) &\cos (\pi/5) &0{\,\vrule height 8pt depth 6pt} &{\hskip -6pt}0 &0\cr\noalign{\vskip -4pt} 0 &0 &1{\,\vrule height 8pt depth 6pt} &{\hskip -6pt}0 &0\cr\noalign{\vskip -4pt} \noalign{\vskip 2pt\hrule\vskip 1pt}0 &0 &0{\,\vrule height 9pt depth 6pt} &{\hskip -6pt}\cos (3\pi/5) &-\sin (3\pi/5)\cr\noalign{\vskip -4pt} 0 &0 &0{\,\vrule height 9pt depth 6pt} &{\hskip -6pt}\sin (3\pi/5) &\cos (3\pi/5)\cr}_{V}\cr &= \pmatrix{\Gamma^{\parallel} (\alpha){\,\,\vrule height 8pt depth 6pt} &{\hskip -10pt}0\cr\noalign{\vskip -2pt} \noalign{\vskip 2pt\hrule}\quad0{\hskip 11pt\vrule height 12pt depth 5pt} &{\hskip -6pt}\Gamma^{\perp} (\alpha)\cr}_{V},\cr \Gamma_{V} (\beta) &= \pmatrix{1 &0 &0{\,\vrule height 8pt depth 6pt} &{\hskip -9pt}0 &0\cr\noalign{\vskip -6pt} 0 &\bar{1} &0{\,\vrule height 8pt depth 6pt} &{\hskip -9pt}0 &0\cr\noalign{\vskip -6pt} 0 &0 &1{\,\vrule height 8pt depth 6pt} &{\hskip -9pt}0 &0\cr\noalign{\vskip -2pt}\noalign{\hrule}\noalign{\vskip 1pt} 0 &0 &0{\,\vrule height 8pt depth 6pt} &{\hskip -9pt}1 &0\cr \noalign{\vskip -6pt}0 &0 &0{\,\vrule height 8pt depth 3pt} &{\hskip -9pt}0 &\bar{1}\cr}_{V}, \quad \Gamma_{V} (\gamma) \pmatrix{\bar1 &0 &0{\,\vrule height 8pt depth 6pt} &{\hskip -9pt}0 &0\cr\noalign{\vskip -6pt} 0 &\bar{1} &0{\,\vrule height 8pt depth 6pt} &{\hskip -9pt}0 &0\cr\noalign{\vskip -6pt} 0 &0 &\bar1{\,\vrule height 8pt depth 6pt} &{\hskip -9pt}0 &0\cr\noalign{\vskip -2pt}\noalign{\hrule}\noalign{\vskip 1pt} 0 &0 &0{\,\vrule height 8pt depth 6pt} &{\hskip -9pt}\bar{1} &0\cr \noalign{\vskip -6pt}0 &0 &0{\,\vrule height 8pt depth 3pt} &{\hskip -9pt}0 &\bar1\cr}_{V},}]where [W = \pmatrix{a_{1}^{*} \cos (2\pi/5) &a_{2}^{*} \cos (4\pi/5) &a_{3}^{*} \cos (6\pi/5) &a_{4}^{*} \cos (8\pi/5) &0\cr a_{1}^{*} \sin (2\pi/5) &a_{2}^{*} \sin (4\pi/5) &a_{3}^{*} \sin (6\pi/5) &a_{4}^{*} \sin (8\pi/5) &0\cr 0 &0 &0 &0 &a_{5}^{*}\cr\noalign{\vskip 2pt}\noalign{\hrule}\noalign{\vskip 4pt} a_{1}^{*} \cos (6\pi/5) &a_{2}^{*} \cos (2\pi/5) &a_{3}^{*} \cos (8\pi/5) &a_{4}^{*} \cos (4\pi/5) &0\cr a_{1}^{*} \sin (6\pi/5) &a_{2}^{*} \sin (2\pi/5) &a_{3}^{*} \sin (8\pi/5) &a_{4}^{*} \sin (4\pi/5) &0\cr}.]The column vectors of the matrix W give the parallel- (above the partition line) and perpendicular-space components (below the partition line) of a reciprocal basis in V space. Thus, W can be rewritten using the physical-space reciprocal basis defined above as [W = \pmatrix{{\bf d}_{1}^{*}, &{\bf d}_{2}^{*}, &{\bf d}_{3}^{*}, &{\bf d}_{4}^{*}, &{\bf d}_{5}^{*}\cr},]yielding the reciprocal basis [{\bf d}_{i}^{*}, i = 1, \ldots, 5], in the 5D embedding space (D space): [{\bf d}_{i}^{*} = a_{i}^{*} \pmatrix{\cos (2\pi i/5)\cr \sin (2\pi i/5)\cr 0\cr \cos (6\pi i/5)\cr \sin (6\pi i/5)\cr}_{V}, \, i = 1, \ldots, 4 \hbox{ and } {\bf d}_{5}^{*} = a_{5}^{*} \pmatrix{0\cr 0\cr 1\cr 0\cr 0\cr}_{V}.]The [5 \times 5] symmetry matrices can each be decomposed into a [3 \times 3] matrix and a [2 \times 2] matrix. The first one, [\Gamma^{\parallel}], acts on the parallel-space component, the second one, [\Gamma^{\perp}], on the perpendicular-space component. In the case of [\Gamma (\alpha)], the coupling factor between a rotation in parallel and perpendicular space is 3. Thus, a [\pi/5] rotation in physical space is related to a [3\pi/5] rotation in perpendicular space (Fig. 4.6.3.12)[link].

With the condition [{\bf d}_{i} \cdot {\bf d}_{j}^{*} = \delta_{ij}], a basis in direct 5D space is obtained: [{\bf d}_{i} = {2 \over 5a_{i}^{*}} \pmatrix{\cos (2\pi i/5) - 1\cr \sin (2\pi i/5)\cr 0\cr \cos (6\pi i/5) - 1\cr \sin (6\pi i/5)\cr}, \, i = 1, \ldots, 4, \hbox{ and } {\bf d}_{5} = {1 \over a_{5}^{*}} \pmatrix{0\cr 0\cr 1\cr 0\cr 0\cr}.]The metric tensors G, [G^{*}] are of the type [\pmatrix{A &C &C &C &0\cr C &A &C &C &0\cr C &C &A &C &0\cr C &C &C &A &0\cr 0 &0 &0 &0 &B\cr}]with [A = 2a_{1}^{*2}, B = a_{5}^{*2}, C = - (1/2) a_{1}^{*2}] for the reciprocal space and [A = 4/5a_{1}^{*2}, B = 1/a_{5}^{*2}, C = 2/5a_{1}^{*2}] for the direct space. Thus, for the lattice parameters in reciprocal space we obtain [d_{i}^{*} = a_{i}^{*} (2)^{1/2}], [i = 1, \ldots, 4]; [d_{5}^{*} = a_{5}^{*}]; [\alpha_{ij}^{*} = 104.5^{\circ}], [i,j = 1, \ldots, 4]; [\alpha_{i5}^{*} = 90^{\circ}], [i = 1, \ldots, 4], and for those in direct space [d_{i} = 2/[a_{i}^{*} (5)^{1/2}]], [i = 1, \ldots, 4]; [d_{5} = 1/a_{5}^{*}]; [\alpha_{ij} = 60^{\circ}], [i,j =] [1, \ldots, 4]; [\alpha_{i5} = 90^{\circ}, i = 1, \ldots, 4]. The volume of the 5D unit cell can be calculated from the metric tensor G: [V = [\det (G)]^{1/2} = {4 \over 5 (5)^{1/2} (a_{1}^{*})^{4} a_{5}^{*}} = {(5)^{1/2} \over 4} (d_{1})^{4} d_{5}.]

Since decagonal phases are only quasiperiodic in two dimensions, it is sufficient to demonstrate their characteristics on a 2D example, the canonical Penrose tiling (Penrose, 1974[link]). It can be constructed from two unit tiles: a skinny (acute angle [\alpha_{\rm s} = \pi/5]) and a fat (acute angle [\alpha_{\rm f} = 2\pi/5]) rhomb with equal edge lengths [a_{\rm r}] and areas [A_{\rm S} = a_{\rm r}^{2} \sin (\pi/5)], [A_{\rm F} = a_{\rm r}^{2} \sin (2\pi/5)] (Fig. 4.6.3.13)[link]. The areas and frequencies of these two unit tiles in the Penrose tiling are both in a ratio 1 to τ. By replacing each skinny and fat rhomb according to the inflation rule, a τ-inflated tiling is obtained. Inflation (deflation) means that the number of tiles is inflated (deflated), their edge lengths are decreased (increased) by a factor τ. By infinite repetition of this inflation operation, an infinite Penrose tiling is generated. Consequently, this substitution operation leaves the tiling invariant.

[Figure 4.6.3.13]

Figure 4.6.3.13 | top | pdf |

A section of a Penrose tiling (thin lines) superposed by its τ-deflated tiling (above, thick lines) and by its [\tau^{2}]-deflated tiling (below, thick lines). In the middle, the inflation rule of the Penrose tiling is illustrated.

From Fig. 4.6.3.13[link] it can be seen that the sets of vertices of the deflated tilings are subsets of the set of vertices of the original tiling. The τ-deflated tiling is dual to the original tiling; a further deflation by a factor τ gives the original tiling again. However, the edge lengths of the tiles are increased by a factor [\tau^{2}], and the tiling is rotated around [36^{\circ}]. Only the fourth deflation of the original tiling yields the original tiling in its original orientation but with all lengths multiplied by a factor [\tau^{4}].

Contrary to the reciprocal-space scaling behaviour of [M^{*} = \{{\bf H}^{\parallel} = {\textstyle\sum_{i = 1}^{4}} h_{i} {\bf a}_{i}^{*} |h_{i} \in {\bb Z}\}], the set of vertices [M =] [\{{\bf r} = {\textstyle\sum_{i = 1}^{4}} n_{i} {\bf a}_{i} | n_{i} \in {\bb Z}\}] of the Penrose tiling is not invariant by scaling the length scale simply by a factor τ using the scaling matrix S: [S = \pmatrix{0 &1 &0 &\bar{1}\cr 0 &1 &1 &\bar{1}\cr \bar{1} &1 &1 &0\cr \bar{1} &0 &1 &0\cr}_{D} \hbox{ acting on vectors } {\bf r} = \pmatrix{n_{1}\cr n_{2}\cr n_{3}\cr n_{4}\cr}_{D}.]The square of S, however, maps all vertices of the Penrose tiling upon other ones: [S^{2} = \pmatrix{1 &1 &0 &\bar{1}\cr 0 &2 &1 &\bar{1}\cr \bar{1} &1 &2 &0\cr \bar{1} &0 &1 &1\cr}_{D}, \quad \Gamma (\alpha) S^{2} = \pmatrix{1 &1 &\bar{1} &\bar{1}\cr 1 &2 &0 &\bar{2}\cr 0 &2 &1 &\bar{1}\cr \bar{1} &1 &1 &0\cr}_{D}.][S^{2}] corresponds to a hyperbolic rotation with [\chi = \cosh^{-1}(3/2)] in superspace (Janner, 1992[link]). However, only operations of the type [S^{4n}], [n = 0, 1, 2 \ldots], scale the Penrose tiling in a way which is equivalent to the (4nth) substitutional operations discussed above. The rotoscaling operation [\Gamma (\alpha) S^{2}], also a symmetry operation of the Penrose tiling, leaves a pentagram invariant as demonstrated in Fig. 4.6.3.14[link] (Janner, 1992[link]). Block-diagonalization of the scaling matrix S decomposes it into two non-equivalent irreducible representations which give the scaling properties in the two orthogonal subspaces of the 4D embedding space, [{\bf V} = {\bf V}^{\parallel} \oplus {\bf V}^{\perp}], the 2D parallel (physical) subspace [{\bf V}^{\parallel}] and the perpendicular 2D subspace [{\bf V}^{\perp}]. Thus, using [WSW^{-1} = S_{V} = S_{V}^{\parallel} \oplus S_{V}^{\perp}], we obtain [S_{V} = \pmatrix{\tau &0{\hskip 3pt}{\vrule height 8pt depth 6pt} &0 &0\cr\noalign{\vskip -3pt}0 &\tau{{\hskip 3pt}\vrule height 8pt depth 6pt} &0 &0\cr\noalign{\vskip-3pt}\noalign{\hrule}\noalign{\vskip 2pt} 0 &0{{\hskip 3pt}\vrule height 8pt depth 6pt} &-1/\tau &0\cr\noalign{\vskip-3pt} 0 &0{\hskip 3pt}{\vrule height 8pt depth 6pt}&0 &-1/\tau\cr}_{V} = \pmatrix{S_{V}^{\parallel}{\hskip  3pt}{\vrule height 12pt depth 6pt}&{\hskip -12pt}0\cr\noalign{\vskip-2pt}\noalign{\hrule}\noalign{\vskip -2pt} 0{\hskip 8.8pt}{\vrule height 12pt depth 6pt} &{\hskip -6pt}S_{V}^{\perp}\cr}_{V},]where [ W = {\pmatrix{a_{1}^{*} \cos (2\pi /5) &a_{2}^{*} \cos (4\pi /5) &a_{3}^{*} \cos (6\pi /5) &a_{4}^{*} \cos (8\pi /5)\cr a_{1}^{*} \sin (2\pi /5) &a_{2}^{*} \sin (4\pi /5) &a_{3}^{*} \sin (6\pi /5) &a_{4}^{*} \sin (8\pi /5)\cr\noalign{\vskip 3pt}\noalign{\hrule}\noalign{\vskip 3pt} a_{1}^{*} \cos (4\pi /5) &a_{2}^{*} \cos (8\pi /5) &a_{3}^{*} \cos (2\pi /5) &a_{4}^{*} \cos (6\pi /5)\cr a_{1}^{*} \sin (4\pi /5) &a_{2}^{*} \sin (8\pi /5) &a_{3}^{*} \sin (2\pi /5) &a_{4}^{*} \sin (6\pi /5)\cr}}.]

[Figure 4.6.3.14]

Figure 4.6.3.14 | top | pdf |

Scaling symmetry of a pentagram superposed on the Penrose tiling. A vector pointing to a corner of a pentagon (star) is mapped by the rotoscaling operation (rotation around [\pi/5] and dilatation by a factor [\tau^{2}]) onto the next largest pentagon (star).

The 2D Penrose tiling can also be embedded canonically in the 5D space. Canonically means that the 5D lattice is hypercubic and that the projection of one unit cell upon the 3D perpendicular space [{\bf V}^{\perp}], giving a rhomb-icosahedron, defines the atomic surface. However, the parallel-space image [{\bf a}_{i}^{*}], [i = 1, \ldots, 4], with [{\bf a}_{0}^{*} = -({\bf a}_{1}^{*} + {\bf a}_{2}^{*} + {\bf a}_{3}^{*} + {\bf a}_{4}^{*})], of the 5D basis [{\bf d}_{i}^{*}, i = 1, \ldots, 4] is not linearly independent. Consequently, the atomic surface consists of only a subset of the points contained in the rhomb-icosahedron: five equidistant pentagons (one with diameter zero) resulting as sections of the rhomb-icosahedron with five equi­distant parallel planes (Fig. 4.6.3.15)[link]. The linear dependence of the 5D basis allows the embedding in the 4D space. The resulting hyper-rhombohedral hyperlattice is spanned by the basis [{\bf d}_{i}], [i = 1, \ldots, 4], discussed above. The atomic surfaces occupy the positions [p/5(1111)], [p = 1,\ldots, 4], on the body diagonal of the 4D unit cell. Neighbouring pentagons are in an anti position to each other (Fig. 4.6.3.16)[link]. Thus the 4D unit cell is decorated centrosymmetrically. The edge length [a_{\rm r}] of a Penrose rhomb is related to the length of physical-space basis vectors [a_{i}^{*}] by [a_{\rm r} = \tau S], with the smallest distance [S = (2\tau/5a_{i}^{*}), i = 1, \ldots, 4]. The point density (number of vertices per unit area) of a Penrose tiling with Penrose rhombs of edge length [a_{\rm r}] can be calculated from the ratio of the relative number of unit tiles in the tiling to their area: [\rho = {1 + \tau \over a_{\rm r}^{2}[\sin(\pi/5) + \tau \sin(2\pi/5)]} = (5/2)a_{i}^{*2} (2 - \tau)^{2} \tan(2\pi/5).]This is equivalent to the calculation from the 4D description, [\eqalign{\rho &= {{\textstyle\sum_{i = 1}^{4}} \Omega_{\rm AS}^{i} \over \Omega_{\rm UC}} = {{\textstyle\sum_{i = 1}^{4}} (5/2) \lambda^{2} \sin(2\pi/5) \over 4/[{5(5)}^{1/2}|a_{i}^{*}|^{4}]} \cr &= (5/2)a_{i}^{*2} (2 - \tau)^{2} \tan(2\pi/5),}]where [\Omega_{\rm AS}] and [\Omega_{\rm UC}] are the area of the atomic surface and the volume of the 4D unit cell, respectively. The pentagon radii are [\lambda_{1, \, 4} = 2(2 - \tau)/5a^{*}] and [\lambda_{2, \, 3} = 2(\tau - 1)/5a^{*}] for the atomic surfaces in [(p/5)(1111)] with [p = 1, 4] and [p = 2, 3]. A detailed discussion of the properties of Penrose tiling is given in the papers of Penrose (1974[link], 1979[link]), Jaric (1986)[link] and Pavlovitch & Kleman (1987)[link].

[Figure 4.6.3.15]

Figure 4.6.3.15 | top | pdf |

Atomic surface of the Penrose tiling in the 5D hypercubic description. The projection of the 5D hypercubic unit cell upon [{\bf V}^{\perp}] gives a rhomb-icosahedron (above). The Penrose tiling is generated by four equidistant pentagons (shaded) inscribed in the rhomb-icosahedron. Below is a perpendicular-space projection of the same pentagons, which are located on the [[1111]_{D}] diagonal of the 4D hyper-rhombohedral unit cell in the 4D description.

[Figure 4.6.3.16]

Figure 4.6.3.16 | top | pdf |

Projection of the 4D hyper-rhombohedral unit cell of the Penrose tiling in the 4D description upon the perpendicular space. In the upper drawing all edges between the 16 corners are shown. In the lower drawing the corners are indexed and the four pentagonal atomic surfaces of the Penrose tiling are shaded.

4.6.3.3.2.1. Indexing

| top | pdf |

The indexing of the submodule [M_{1}^{*}] of the diffraction pattern of a decagonal phase is not unique. Since [M_{1}^{*}] corresponds to a [{\bb Z}] module of rank 4 with decagonal point symmetry, it is invariant under scaling by [\tau^{n}, n \in {\bb Z}]: [S^{n}M^{*} = \tau^{n} M^{*}]. Nevertheless, an optimum basis (low indices are assigned to strong reflections) can be derived: not the metrics, as for regular periodic crystals, but the intensity distribution characterizes the best choice of indexing.

A correct set of reciprocal-basis vectors can be identified experimentally in the following way:

  • (1) Find directions of systematic absences or pseudo-absences determining the possible orientations of the reciprocal-basis vectors (see Rabson et al., 1991[link]).

  • (2) Find pairs of strong reflections whose physical-space diffraction vectors are related to each other by the factor τ.

  • (3) Index these reflections by assigning an appropriate value to [a^{*}]. This value should be derived from the shortest interatomic distance S and the edge length of the unit tiles expected in the structure.

  • (4) The reciprocal basis is correct if all observable Bragg reflections can be indexed with integer numbers.

4.6.3.3.2.2. Diffraction symmetry

| top | pdf |

The diffraction symmetry of decagonal phases can be described by the Laue groups [10/mmm] or [10/m]. The set of all vectors H forms a Fourier module [M^{*} = \{{\bf H}^{\parallel} = {\textstyle\sum_{i = 1}^{5}} h_{i}{\bf a}_{i}^{*} | h_{i} \in {\bb Z}\}] of rank 5 in physical space which can be decomposed into two submodules [M^{*} = M_{1}^{*} \oplus M_{2}^{*}]. [M_{1}^{*} = \{h_{1}{\bf a}_{1}^{*} + h_{2}{\bf a}_{2}^{*} + h_{3}{\bf a}_{3}^{*} + h_{4}{\bf a}_{4}^{*}\}] corresponds to a [{\bb Z}] module of rank 4 in a 2D subspace, [M_{2}^{*} = \{h_{5}{\bf a}_{5}^{*}\}] corresponds to a [{\bb Z}] module of rank 1 in a 1D subspace. Consequently, the first submodule can be considered as a projection from a 4D reciprocal lattice, [M_{1}^{*} = \pi^{\parallel} (\Sigma^{*})], while the second submodule is of the form of a regular 1D reciprocal lattice, [M_{2}^{*} = \Lambda^{*}]. The diffraction pattern of the Penrose tiling decorated with equal point scatterers on its vertices is shown in Fig. 4.6.3.17[link]. All Bragg reflections within [10^{-2} |F({\bf 0})|^{2} \,\lt\, |F({\bf H})|^{2} \,\lt\, |F({\bf 0})|^{2}] are depicted. Without intensity-truncation limit, the diffraction pattern would be densely filled with discrete Bragg reflections. To illustrate their spatial and intensity distribution, an enlarged section of Fig. 4.6.3.17[link] is shown in Fig. 4.6.3.18[link]. This picture shows all Bragg reflections within [10^{-4} |F({\bf 0})|^{2} \,\lt\, |F({\bf H})|^{2} \,\lt\, |F({\bf 0})|^{2}]. The projected 4D reciprocal-lattice unit cell is drawn and several reflections are indexed. All reflections are arranged along lines in five symmetry-equivalent orientations. The perpendicular-space diffraction patterns (Figs. 4.6.3.19[link] and 4.6.3.20[link]) show a characteristic star-like distribution of the Bragg reflections. This is a consequence of the pentagonal shape of the atomic surfaces: the Fourier transform of a pentagon has a star-like distribution of strong Fourier coefficients.

[Figure 4.6.3.17]

Figure 4.6.3.17 | top | pdf |

Schematic diffraction pattern of the Penrose tiling (edge length of the Penrose unit rhombs ar = 4.04 Å). All reflections are shown within [10^{-2} |F({\bf 0})|^{2} \lt |F({\bf H})|^{2} \lt |F({\bf 0})|^{2}] and [0 \leq |{\bf H}^{\parallel}| \leq] [2.5\;\hbox{\AA}^{-1}].

[Figure 4.6.3.18]

Figure 4.6.3.18 | top | pdf |

Enlarged section of Fig. 4.6.3.17[link]. All reflections shown are selected within the given limits from a data set within [10^{-4} |F({\bf 0})|^{2} \;\lt\; |F({\bf H})|^{2} \;\lt\; |F({\bf 0})|^{2}] and [0 \leq |{\bf H}^{\parallel}| \leq 2.5\;\hbox{\AA}^{-1}]. The projected 4D reciprocal-lattice unit cell is drawn and several reflections are indexed.

[Figure 4.6.3.19]

Figure 4.6.3.19 | top | pdf |

The perpendicular-space diffraction pattern of the Penrose tiling (edge length of the Penrose unit rhombs ar = 4.04 Å). All reflections are shown within [10^{-4} |F({\bf 0})|^{2} \;\lt\; |F({\bf H})|^{2} \;\lt\; |F({\bf 0})|^{2}] and [0 \leq |{\bf H}^{\parallel}| \leq 2.5\;\hbox{\AA}^{-1}].

[Figure 4.6.3.20]

Figure 4.6.3.20 | top | pdf |

Enlarged section of Fig. 4.6.3.19[link] showing the projected 4D reciprocal-lattice unit cell.

The 5D decagonal space groups that may be of relevance for the description of decagonal phases are listed in Table 4.6.3.1.[link] These space groups are a subset of all 5D decagonal space groups fulfilling the condition that the 5D point groups they are associated with are isomorphous to the 3D point groups describing the diffraction symmetry. Their structures are comparable to 3D hexagonal groups. Hence, only primitive lattices exist. The orientation of the symmetry elements in the 5D space is defined by the isomorphism of the 3D and 5D point groups. However, the action of the tenfold rotation is different in the subspaces [{\bf V}^{\parallel}] and [{\bf V}^{\perp}]: a rotation of [\pi/5] in [{\bf V}^{\parallel}] is correlated with a rotation of [3\pi/5] in [{\bf V}^{\perp}]. The reflection and inversion operations are equivalent in both subspaces.

Table 4.6.3.1| top | pdf |
3D point groups of order k describing the diffraction symmetry and corresponding 5D decagonal space groups with reflection conditions (see Rabson et al., 1991[link])

3D point groupk5D space groupReflection condition
[\displaystyle{{10 \over m} {2 \over m} {2 \over m}}] 40 [\displaystyle{P{10 \over m} {2 \over m} {2 \over m}}] No condition
[\displaystyle{P{10 \over m} {2 \over c} {2 \over c}}] [h_{1}h_{2}h_{2}h_{1}h_{5}]: [h_{5} = 2n]
[h_{1}h_{2}\overline{h_{2}}\overline{h_{1}}h_{5}]: [h_{5} = 2n]
[\displaystyle{P{10_{5} \over m} {2 \over m} {2 \over c}}] [h_{1}h_{2}\overline{h_{2}}\overline{h_{1}}h_{5}]: [h_{5} = 2n]
[\displaystyle{P{10_{5} \over m} {2 \over c} {2 \over m}}] [h_{1}h_{2}h_{2}h_{1}h_{5}]: [h_{5} = 2n]
[\displaystyle{{10 \over m}}] 20 [\displaystyle{P{10 \over m}}] No condition
[\displaystyle{P{10_{5} \over m}}] [0000h_{5}]: [h_{5} = 2n]
1022 20 P1022 No condition
[P10_{j}22] [0000h_{5}]: [jh_{5} = 10n]
10mm 20 P10mm No condition
P10cc [h_{1}h_{2}h_{2}h_{1}h_{5}]: [h_{5} = 2n]
  [h_{1}h_{2}\overline{h_{2}}\overline{h_{1}}h_{5}]: [h_{5} = 2n]
[P10_{5}mc] [h_{1}h_{2}\overline{h_{2}}\overline{h_{1}}h_{5}]: [h_{5} = 2n]
[P10_{5}cm] [h_{1}h_{2}h_{2}h_{1}h_{5}]: [h_{5} = 2n]
[\overline{10}m2] 20 [P\overline{10}m2] No condition
[P\overline{10}c2] [h_{1}h_{2}h_{2}h_{1}h_{5}]: [h_{5} = 2n]
[P\overline{10}2m] No condition
[P\overline{10}2c] [h_{1}h_{2}\overline{h_{2}}\overline{h_{1}}h_{5}]: [h_{5} = 2n]
[\overline{10}] 10 [P\overline{10}] No condition
10 10 P10 No condition
[P10_{j}] [0000h_{5}]: [jh_{5} = 10n]

4.6.3.3.2.3. Structure factor

| top | pdf |

The structure factor for the decagonal phase corresponds to the Fourier transform of the 5D unit cell, [F({\bf H}) = {\textstyle\sum\limits_{k = 1}^{N}} \,f_{k}({\bf H}^{\parallel})T_{k}({\bf H}^{\parallel}, {\bf H}^{\perp})g_{k}({\bf H}^{\perp}) \exp (2\pi i {\bf H}\cdot {\bf r}_{k}),]with 5D diffraction vectors [{\bf H} = {\textstyle\sum_{i = 1}^{5}} h_{i}{\bf d}_{i}^{*}], N hyperatoms, parallel-space atomic scattering factor [f_{k}({\bf H}^{\parallel})], temperature factor [T_{k}({\bf H}^{\parallel}, {\bf H}^{\perp})] and perpendicular-space geometric form factor [g_{k}({\bf H}^{\perp})]. [T_{k}({\bf H}^{\parallel}, {\bf 0})] is equivalent to the conventional Debye–Waller factor and [T_{k}({\bf 0},{\bf H}^{\perp})] describes random fluctuations along the perpendicular-space coordinate. These fluctuations cause characteristic jumps of vertices in physical space (phason flips). Even random phason flips map the vertices onto positions which can still be described by physical-space vectors of the type [{\bf r} = {\textstyle\sum_{i = 1}^{5}} n_{i}{\bf a}_{i}]. Consequently, the set [M = \{{\bf r} = {\textstyle\sum_{i = 1}^{5}} n_{i}{\bf a}_{i} |n_{i} \in {\bb Z}\}] of all possible vectors forms a [{\bb Z}] module. The shape of the atomic surfaces corresponds to a selection rule for the positions actually occupied. The geometric form factor [g_{k}({\bf H}^{\perp})] is equivalent to the Fourier transform of the atomic surface, i.e. the 2D perpendicular-space component of the 5D hyperatoms.

For example, the canonical Penrose tiling [g_{k}({\bf H}^{\perp})] corresponds to the Fourier transform of pentagonal atomic surfaces: [g_{k}({\bf H}^{\perp}) = (1/A_{\rm UC}^{\perp}) {\textstyle\int\limits_{A_{k}}} \exp (2\pi i {\bf H}^{\perp} \cdot {\bf r}) \hbox{ d}{\bf r},]where [A_{\rm UC}^{\perp}] is the area of the 5D unit cell projected upon [{\bf V}^{\perp}] and [A_{k}] is the area of the kth atomic surface. The area [A_{\rm UC}^{\perp}] can be calculated using the formula [A_{\rm UC}^{\perp} = (4/25a_{i}^{*2}) [(7 + \tau) \sin (2\pi/5) + (2 + \tau) \sin (4\pi/5)].]Evaluating the integral by decomposing the pentagons into triangles, one obtains [\eqalign{ g_{k}({\bf H}^{\perp}) &={1 \over A_{\rm UC}^{\perp}} \sin \left({2\pi \over 5}\right)\cr &\quad\times\sum\limits_{j = 0}^{4} {A_{j} [\exp (iA_{j + 1}\lambda_{k}) - 1] - A_{j + 1} [\exp (iA_{j}\lambda_{k}) - 1] \over A_{j}A_{j + 1} (A_{j} - A_{j + 1})}}]with [j = 0, \ldots, 4] running over the five triangles, where the radii of the pentagons are [\lambda_{j}], [A_{j} = 2\pi {\bf H}^{\perp}{\bf e}_{j}], [{\bf H}^{\perp} = \pi^{\perp}({\bf H}) = {\textstyle\sum\limits_{j = 0}^{4}} h_{j}a_{j}^{*} \pmatrix{0\cr 0\cr 0\cr \cos (6\pi j/5)\cr \sin (6\pi j/5)\cr}]and the vectors [{\bf e}_{j} = {1 \over a_{j}^{*}} \pmatrix{0\cr 0\cr 0\cr \cos (2\pi j/5)\cr \sin (2\pi j/5)\cr} \hbox{ with } j = 0, \ldots, 4.]

As shown by Ishihara & Yamamoto (1988)[link], the Penrose tiling can be considered to be a superstructure of a pentagonal tiling with only one type of pentagonal atomic surface in the nD unit cell. Thus, for the Penrose tiling, three special reflection classes can be distinguished: for [|{\textstyle\sum_{i = 1}^{4}} h_{i}| = m\hbox{ mod }5] and [m = 0] the class of strong main reflections is obtained, and for [m = \pm 1, \pm 2] the classes of weaker first- and second-order satellite reflections are obtained (see Fig. 4.6.3.18[link]).

4.6.3.3.2.4. Intensity statistics

| top | pdf |

This section deals with the reciprocal-space characteristics of the 2D quasiperiodic component of the 3D structure, namely the Fourier module [M_{1}^{*}]. The radial structure-factor distributions of the Penrose tiling decorated with point scatterers are plotted in Figs. 4.6.3.21[link] and 4.6.3.22[link] as a function of parallel and perpendicular space. The distribution of [|F({\bf H})|] as a function of their frequencies clearly resembles a centric distribution, as can be expected from the centrosymmetric 4D subunit cell. The shape of the distribution function depends on the radius of the limiting sphere in reciprocal space. The number of weak reflections increases to the power of four, that of strong reflections only quadratically (strong reflections always have small [{\bf H}^{\perp}] components). The radial distribution of the structure-factor amplitudes as a function of perpendicular space clearly shows three branches, corresponding to the reflection classes [{\textstyle\sum_{i = 1}^{4}} h_{i} = m\hbox{ mod }5] with [|m| = 0], [|m| = 1] and [|m| = 2] (Fig. 4.6.3.23)[link].

[Figure 4.6.3.21]

Figure 4.6.3.21 | top | pdf |

Radial distribution function of the structure factors [F({\bf H})] of the Penrose tiling (edge length of the Penrose unit rhombs ar = 4.04 Å) decorated with point atoms as a function of [{\bf H}^{\parallel}]. All structure factors within [10^{-4} |F({\bf 0})|^{2} \;\lt\; |F({\bf H})|^{2} \;\lt\; |F({\bf 0})|^{2}] and [0 \leq |{\bf H}^{\parallel}| \leq 2.5\;\hbox{\AA}^{-1}] have been used and normalized to [F(0000) = 1].

[Figure 4.6.3.22]

Figure 4.6.3.22 | top | pdf |

Radial distribution function of the structure factors [F({\bf H})] of the Penrose tiling (edge length of the Penrose unit rhombs ar = 4.04 Å) decorated with point atoms as a function of [{\bf H}^{\perp}]. All structure factors within [10^{-4} |F({\bf 0})|^{2} \;\lt\; |F({\bf H})|^{2} \;\lt\; |F({\bf 0})|^{2}] and [0 \leq |{\bf H}^{\parallel}| \leq 2.5\;\hbox{\AA}^{-1}] have been used and normalized to [F(0000) = 1].

[Figure 4.6.3.23]

Figure 4.6.3.23 | top | pdf |

Radial distribution function of the structure-factor magnitudes [|F({\bf H})|] of the Penrose tiling (edge length of the Penrose unit rhombs ar = 4.04 Å) decorated with point atoms as a function of [{\bf H}^{\perp}]. All structure factors within [10^{-4} |F({\bf 0})|^{2} \;\lt\; |F({\bf H})|^{2} \;\lt\; |F({\bf 0})|^{2}] and [0 \leq |{\bf H}^{\parallel}| \leq 2.5\;\hbox{\AA}^{-1}] have been used and normalized to [F(0000) = 1]. The branches with (a) [|{\textstyle\sum_{i = 1}^{4}} h_{i}| = 0\hbox{ mod }5], (b) [|{\textstyle\sum_{i = 1}^{4}} h_{i}| = 1\hbox{ mod }5] and (c) [|{\textstyle\sum_{i = 1}^{4}} h_{i}| = 2\hbox{ mod }5] are shown.

The weighted reciprocal space of the Penrose tiling contains an infinite number of Bragg reflections within a limited region of the physical space. Contrary to the diffraction pattern of a periodic structure consisting of point atoms on the lattice nodes, the Bragg reflections show intensities depending on the perpendicular-space components of their diffraction vectors (Figs. 4.6.3.19[link], 4.6.3.20[link] and 4.6.3.22[link]).

4.6.3.3.2.5. Relationships between structure factors at symmetry-related points of the Fourier image

| top | pdf |

Scaling the Penrose tiling by a factor [\tau^{-n}] by employing the matrix [S^{-n}] scales at the same time its reciprocal space by a factor [\tau^{n}]: [\eqalignno{ S{\bf H} &= \pmatrix{0 &1 &0 &\bar{1} &0\cr 0 &1 &1 &\bar{1} &0\cr \bar{1} &1 &1 &0 &0\cr \bar{1} &0 &1 &0 &0\cr 0 &0 &0 &0 &1\cr}_{D} \pmatrix{h_{1}\cr h_{2}\cr h_{3}\cr h_{4}\cr h_{5}\cr} = \pmatrix{h_{2} - h_{4}\cr h_{2} + h_{3} - h_{4}\cr - h_{1} + h_{2} + h_{3}\cr - h_{1} + h_{3}\cr h_{5}\cr}. &\cr}]Since this operation increases the lengths of the diffraction vectors by the factor τ in parallel space and decreases them by the factor [1/\tau] in perpendicular space, the following distribution of structure-factor magnitudes (for point atoms at rest) is obtained: [\displaylines{ \hfill \left|F(S^{n}{\bf H})\right| > \left|F(S^{n - 1}{\bf H})\right| > \ldots > \left|F(S^{1}{\bf H})\right| > \left|F({\bf H})\right|, \hfill \cr \hfill \left|F(\tau^{n}{\bf H}^{\|})\right| > \left|F(\tau^{n - 1}{\bf H}^{\|})\right| > \ldots > \left|F(\tau {\bf H}^{\|})\right| > \left|F({\bf H})\right|. \hfill \cr}]The scaling operations [S^{n}], [n \in {\bb Z}], the rotoscaling operations [(\Gamma (\alpha)S^{2})^{n}] (Fig. 4.6.3.14)[link] and the tenfold rotation [(\Gamma (\alpha))^{n}], where [(\Gamma (\alpha)S^{2})^{n} = \pmatrix{1 &1 &\bar{1} &\bar{1} &0\cr 1 &2 &0 &\bar{2} &0\cr 0 &2 &1 &\bar{1} &0\cr \bar{1} &1 &1 &0 &0\cr 0 &0 &0 &0 &1\cr}^{n}_{D},]connect all structure factors with diffraction vectors pointing to the nodes of an infinite series of pentagrams. The structure factors with positive signs are predominantly on the vertices of the pentagram while the ones with negative signs are arranged on circles around the vertices (Figs. 4.6.3.24[link][link][link] to 4.6.3.27[link]).

[Figure 4.6.3.24]

Figure 4.6.3.24 | top | pdf |

Pentagrammal relationships between scaling symmetry-related positive structure factors [F({\bf H})] of the Penrose tiling (edge length ar = 4.04 Å) in parallel space. The magnitudes of the structure factors are indicated by the diameters of the filled circles.

[Figure 4.6.3.25]

Figure 4.6.3.25 | top | pdf |

Pentagrammal relationships between scaling symmetry-related negative structure factors [F({\bf H})] of the Penrose tiling (edge length ar = 4.04 Å) in parallel space. The magnitudes of the structure factors are indicated by the diameters of the filled circles.

[Figure 4.6.3.26]

Figure 4.6.3.26 | top | pdf |

Pentagrammal relationships between scaling symmetry-related structure factors [F({\bf H})] of the Penrose tiling (edge length ar = 4.04 Å) in parallel space. Enlarged sections of Figs. 4.6.3.24[link] (above) and 4.6.3.25[link] (below) are shown.

[Figure 4.6.3.27]

Figure 4.6.3.27 | top | pdf |

Pentagrammal relationships between scaling symmetry-related structure factors [F({\bf H})] of the Penrose tiling (edge length ar = 4.04 Å) in perpendicular space. Enlarged sections of positive (above) and negative structure factors (below) are shown.

4.6.3.3.3. Icosahedral phases

| top | pdf |

A structure that is quasiperiodic in three dimensions and exhibits icosahedral diffraction symmetry is called an icosahedral phase. Its holohedral Laue symmetry group is [K = m\bar{3}\bar{5}]. All reciprocal-space vectors [{\bf H} = \sum_{i = 1}^{6} h_{i}{\bf a}_{i}^{*} \in M^{*}] can be represented on a basis [{\bf a}_{1}^{*} = a^{*}(0,0,1)], [{\bf a}_{i}^{*} =] [a^{*} [\sin \theta \cos (2\pi i/5), \sin \theta \sin (2\pi i/5), \cos \theta]], [i = 2, \ldots, 6] where [\sin \theta = 2/(5)^{1/2}], [\cos \theta = 1/(5)^{1/2}] and [\theta \simeq 63.44^{\circ}], the angle between two neighbouring fivefold axes (Fig. 4.6.3.28)[link]. This can be rewritten as [\pmatrix{{\bf a}_{1}^{*}\cr {\bf a}_{2}^{*}\cr {\bf a}_{3}^{*}\cr {\bf a}_{4}^{*}\cr {\bf a}_{5}^{*}\cr {\bf a}_{6}^{*}\cr} = a^{*} \pmatrix{0 &0 &1\cr \sin \theta \cos (4\pi/5) &\sin \theta \sin (4\pi/5) &\cos \theta\cr \sin \theta \cos (6\pi/5) &\sin \theta \sin (6\pi/5) &\cos \theta\cr \sin \theta \cos (8\pi/5) &\sin \theta \sin (8\pi/5) &\cos \theta\cr \sin \theta &0 &\cos \theta\cr \sin \theta \cos (2\pi/5) &\sin \theta \sin (2\pi/5) &\cos \theta\cr} \openup3pt{\pmatrix{{\bf e}_{1}^{V}\cr {\bf e}_{2}^{V}\cr {\bf e}_{3}^{V}\cr}},]where [{\bf e}_{i}^{V}] are Cartesian basis vectors. Thus, from the number of independent reciprocal-basis vectors needed to index the Bragg reflections with integer numbers, the dimension of the embedding space has to be six. The vector components refer to a Cartesian coordinate system (V basis) in the physical (parallel) space.

[Figure 4.6.3.28]

Figure 4.6.3.28 | top | pdf |

Perspective (a) parallel- and (b) perpendicular-space views of the reciprocal basis of the 3D Penrose tiling. The six rationally independent vectors [{\bf a}_{i}^{*}] point to the edges of an icosahedron.

The set [M^{*} = \{{\bf H}^{\parallel} = {\textstyle\sum_{i = 1}^{6}} h_{i} {\bf a}_{i}^{*} | h_{i} \in {\bb Z}\}] of all diffraction vectors remains invariant under the action of the symmetry operators of the icosahedral point group [K = m\bar{3}\bar{5}]. The symmetry-adapted matrix representations for the point-group generators, one fivefold rotation α, a threefold rotation β and the inversion operation γ, can be written in the form [\displaylines{\Gamma (\alpha) = \pmatrix{1 &0 &0 &0 &0 &0\cr 0 &0 &0 &0 &0 &1\cr 0 &1 &0 &0 &0 &0\cr 0 &0 &1 &0 &0 &0\cr 0 &0 &0 &1 &0 &0\cr 0 &0 &0 &0 &1 &0\cr}_{D}, \Gamma (\beta) = \pmatrix{0 &1 &0 &0 &0 &0\cr 0 &0 &0 &0 &0 &1\cr 0 &0 &0 &\bar{1} &0 &0\cr 0 &0 &0 &0 &\bar{1} &0\cr 0 &0 &1 &0 &0 &0\cr 1 &0 &0 &0 &0 &0\cr}_{D},\cr \Gamma (\gamma) = \pmatrix{\bar{1} &0 &0 &0 &0 &0\cr 0 &\bar{1} &0 &0 &0 &0\cr 0 &0 &\bar{1} &0 &0 &0\cr 0 &0 &0 &\bar{1} &0 &0\cr 0 &0 &0 &0 &\bar{1} &0\cr 0 &0 &0 &0 &0 &\bar{1}\cr}_{D}.}]

Block-diagonalization of these reducible symmetry matrices decomposes them into non-equivalent irreducible representations. These can be assigned to the two orthogonal subspaces forming the 6D embedding space [{\bf V} = {\bf V}^{\parallel} \oplus {\bf V}^{\perp}], the 3D parallel (physical) subspace [{\bf V}^{\parallel}] and the perpendicular 3D subspace [{\bf V}^{\perp}]. Thus, using [W \Gamma W^{-1} = \Gamma^{\rm red} = \Gamma^{\parallel} \oplus \Gamma^{\perp}], we obtain[\eqalign{\Gamma (\alpha) &= \pmatrix{\cos (2\pi/5) &- \sin (2\pi/5) &0\noalign\,{\vrule height 8pt depth 6pt} &0 &0 &0\cr\noalign{\vskip -4pt} \sin (2\pi/5) &\cos (2\pi/5) &0\noalign\,{\vrule height 8pt depth 6pt} &0 &0 &0\cr\noalign{\vskip -4pt} 0 &0 &1\noalign\,{\vrule height 8pt depth 6pt} &0 &0 &0\cr\noalign{\vskip -2pt}\noalign{\hrule}\noalign{\vskip -1pt} 0 &0 &0\noalign\,{\vrule height 12pt depth 6pt} &\cos (4\pi/5) &- \sin (4\pi/5) &0\cr\noalign{\vskip -4pt} 0 &0 &0\noalign\,{\vrule height 8pt depth 6pt} &\sin (4\pi/5) &\cos (4\pi/5) &0\cr\noalign{\vskip -4pt} 0 &0 &0\noalign\,{\vrule height 8pt depth 6pt} &0 &0 &1\cr}_{V}\cr &= \pmatrix{\Gamma^{\parallel}\noalign\,{\vrule height 8pt depth 6pt} &{\hskip -10pt}0\cr\noalign{\vskip -2pt}\noalign{\hrule}\noalign{\vskip -3pt} 0\noalign{\hskip 7pt}{\vrule height 12pt depth 6pt} &{\hskip-4pt}\Gamma^{\perp}\cr}_{V},}]where [W = a^{*} \pmatrix{0 &sc4 &sc6 &sc8 &s &sc2\cr 0 &ss4 &ss6 &ss8 &0 &ss2\cr 1 &c &c &c &c &c\cr\noalign{\vskip 3pt}\noalign{\hrule}\noalign{\vskip 3pt} 0 &- sc8 &- sc2 &- sc6 &- s &- sc4\cr 0 &- ss8 &- ss2 &- ss6 &0 &- ss4\cr 1 &- c &- c &- c &- c &- c\cr}_{V},][c = \cos \theta], [ s = \sin \theta], [scn = \sin \theta \cos (n\pi/5)], [ ssn = \sin \theta \sin (n\pi/5)]. The column vectors of the matrix W give the parallel- (above the partition line) and perpendicular-space components (below the partition line) of a reciprocal basis in V. Thus, W can be rewritten using the physical-space reciprocal basis defined above and an arbitrary constant c, [\eqalign{W &= \pmatrix{{\bf a}_{1}^{*} &{\bf a}_{2}^{*} &{\bf a}_{3}^{*} &{\bf a}_{4}^{*} &{\bf a}_{5}^{*} &{\bf a}_{6}^{*}\cr c{\bf a}_{1}^{*} &- c{\bf a}_{4}^{*} &- c{\bf a}_{6}^{*} &- c{\bf a}_{3}^{*} &- c{\bf a}_{5}^{*} &- c{\bf a}_{2}^{*}\cr}\cr\noalign{\vskip 9pt} &= \pmatrix{{\bf d}_{1}^{*} &{\bf d}_{2}^{*} &{\bf d}_{3}^{*} &{\bf d}_{4}^{*} &{\bf d}_{5}^{*} &{\bf d}_{6}^{*}\cr},}]yielding the reciprocal basis [{\bf d}_{i}^{*}, i = 1, \ldots, 6], in the 6D embedding space (D space) [{\bf d}_{1}^{*} = a^{*} \pmatrix{0\cr 0\cr 1\cr 0\cr 0\cr c\cr} \hbox{ and } {\bf d}_{i}^{*} = a^{*} \pmatrix{\sin \theta \cos (2\pi i/5)\cr \sin \theta \sin (2\pi i/5)\cr \cos \theta\cr - c \sin \theta \cos (4\pi i/5)\cr - c \sin \theta \sin (4\pi i/5)\cr - c \cos \theta\cr}, i = 2, \ldots, 6.]The [6\times 6] symmetry matrices can each be decomposed into two [3\times 3] matrices. The first one, [\Gamma^{\parallel}], acts on the parallel-space component, the second one, [\Gamma^{\perp}], on the perpendicular-space component. In the case of [\Gamma (\alpha)], the coupling factor between a rotation in parallel and perpendicular space is 2. Thus a [2\pi /5] rotation in physical space is related to a [4\pi /5] rotation in perpendicular space (Figs. 4.6.3.28[link] and 4.6.3.29[link]).

[Figure 4.6.3.29]

Figure 4.6.3.29 | top | pdf |

Schematic representation of a rotation in 6D space. The point P is rotated to P′. The component rotations in parallel and perpendicular space are illustrated.

With the condition [{\bf d}_{i} \cdot {\bf d}_{j}^{*} = \delta_{ij}] , the basis in direct 6D space is obtained: [\displaylines{{\bf d}_{1} = {1 \over 2a^{*}} \pmatrix{0\cr 0\cr 1\cr 0\cr 0\cr 1/c\cr} \hbox{ and } {\bf d}_{i} = {1 \over 2a^{*}} \pmatrix{\sin \theta \cos (2\pi i/5)\cr \sin \theta \sin (2\pi i/5)\cr \cos \theta\cr - (1/c) \sin \theta \cos (4\pi i/5)\cr - (1/c) \sin \theta \sin (4\pi i/5)\cr - (1/c) \cos \theta\cr},\hfill\cr\hfill i = 2, \ldots, 6.}]The metric tensors G, [G^{*}] are of the type [\pmatrix{A &\hfill B &\hfill B &\hfill B &\hfill B &\hfill B\cr B &\hfill A &\hfill B &-B &-B &\hfill B\cr B &\hfill B &\hfill A &\hfill B &-B &-B\cr B &-B &\hfill B &\hfill A &\hfill B &-B\cr B &-B &-B &\hfill B &\hfill A &\hfill B\cr B &\hfill B &-B &-B &\hfill B &\hfill A\cr},]with [A = (1 + c^{2}) a^{*2}, B = [(5)^{1/2}/5] (1 - c^{2}) a^{*2}] for the reciprocal space and [A = (1 + c^{2})/[4 (ca^{*})^{2}]], [B = [(5)^{1/2} (c^{2} - 1)]/[20 (ca^{*})^{2}]] for the direct space. For [c = 1] we obtain hypercubic direct and reciprocal 6D lattices.

The lattice parameters in reciprocal and direct space are [d_{i}^{*} = a^{*} (2)^{1/2}] and [d_{i} = 1/[(2)^{1/2} a^{*}]] with [i = 1, \ldots, 6], respectively. The volume of the 6D unit cell can be calculated from the metric tensor G. For [c = 1] it is simply [V = [\det (G)]^{1/2} = \{1/[(2)^{1/2} a^{*}]\}^{6}.]

The best known example of a 3D quasiperiodic structure is the canonical 3D Penrose tiling (see Janssen, 1986[link]). It can be constructed from two unit tiles: a prolate and an oblate rhombo­hedron with equal edge lengths [a_{\rm r}] (Fig. 4.6.3.30)[link]. Each face of the rhombohedra is a rhomb with acute angles [\alpha_{\rm r} = \hbox{arc } \cos [1/(5)^{1/2}] \simeq 63.44^{\circ}]. Their volumes are [V_{\rm p} =] [(4/5) a_{\rm r}^{3} \sin (2\pi /5), V_{\rm o} = (4/5) a_{\rm r}^{3} \sin (\pi /5) = V_{\rm p}/\tau], and their frequencies [\nu_{\rm p}]:[\nu_{\rm o} = \tau]:1. The resulting point density (number of vertices per unit volume) is [\rho_{\rm p} = (\tau + 1)/(\tau V_{\rm p} + V_{\rm o}) =] [(\tau /a_{\rm r}^{3}) \sin (2\pi /5)]. Ten prolate and ten oblate rhombohedra can be packed to form a rhombic triacontahedron. The icosahedral symmetry of this zonohedron is broken by the many possible decompositions into the rhombohedra. Removing one zone of the triacontahedron gives a rhomb-icosahedron consisting of five prolate and five oblate rhombohedra. Again, the singular fivefold axis of the rhomb-icosahedron is broken by the decomposition into rhombohedra. Removing one zone again gives a rhombic dodecahedron consisting of two prolate and two oblate rhombo­hedra. Removing the last remaining zone leads finally to a single prolate rhombohedron. Using these zonohedra as elementary clusters, a matching rule can be derived for the 3D construction of the 3D Penrose tiling (Levine & Steinhardt, 1986[link]; Socolar & Steinhardt, 1986[link]).

[Figure 4.6.3.30]

Figure 4.6.3.30 | top | pdf |

The two unit tiles of the 3D Penrose tiling: a prolate [[\alpha_{\rm p} = \hbox{arc} \cos (5^{-1/2}) \simeq 63.44^{\circ}]] and an oblate [(\alpha_{\rm o} = 180^{\circ} - \alpha_{\rm p})] rhombo­hedron with equal edge lengths [a_{\rm r}].

The 3D Penrose tiling can be embedded in the 6D space as shown above. The 6D hypercubic lattice is decorated on the lattice nodes with 3D triacontahedra obtained from the projection of a 6D unit cell upon the perpendicular space [{\bf V}^{\perp}] (Fig. 4.6.3.31)[link]. Thus the edge length of the rhombs covering the triacontahedron is equivalent to the length [\pi^{\perp} ({\bf d}_{i}) = 1/2a^{*}] of the perpendicular-space component of the vectors spanning the 6D hypercubic lattice [\Sigma = \{{\bf r} = {\textstyle\sum_{i = 1}^{6}} n_{i} {\bf d}_{i} | n_{i} \in {\bb Z}\}].

[Figure 4.6.3.31]

Figure 4.6.3.31 | top | pdf |

Atomic surface of the 3D Penrose tiling in the 6D hypercubic description. The projection of the 6D hypercubic unit cell upon [{\bf V}^{\perp}] gives a rhombic triacontahedron.

4.6.3.3.3.1. Indexing

| top | pdf |

There are several indexing schemes in use. The generic one uses a set of six rationally independent reciprocal-basis vectors pointing to the corners of an icosahedron, [{\bf a}_{1}^{*} = a^{*} (0,0,1)], [{\bf a}_{i}^{*} = a^{*} [\sin \theta \cos (2\pi i/5), \sin \theta \sin (2\pi i/5), \cos \theta]], [i = 2, \ldots, 6], [\sin \theta = 2/(5)^{1/2}, \cos \theta = 1/(5)^{1/2}], with [\theta \simeq 63.44^{\circ}], the angle between two neighbouring fivefold axes (setting 1) (Fig. 4.6.3.28)[link]. In this case, the physical-space basis corresponds to a simple projection of the 6D reciprocal basis [{\bf d}_{i}^{*}, i = 1, \ldots, 6]. Sometimes, the same set of six reciprocal-basis vectors is referred to a differently oriented Cartesian reference system (C basis, with basis vectors [{\bf e}_{i}] along the twofold axes) (Bancel et al., 1985[link]). The reciprocal basis is [\pmatrix{{\bf a}_{1}^{*}\cr {\bf a}_{2}^{*}\cr {\bf a}_{3}^{*}\cr {\bf a}_{4}^{*}\cr {\bf a}_{5}^{*}\cr {\bf a}_{6}^{*}\cr} = {a^{*} \over (1 + \tau^{2})^{1/2}} \pmatrix{0 &1 &\tau\cr -1 &\tau &0\cr -\tau &0 &1\cr 0 &-1 &\tau\cr \tau &0 &1\cr 1 &\tau &0\cr}_{C} \pmatrix{{\bf e}_{1}^{C}\cr {\bf e}_{2}^{C}\cr {\bf e}_{3}^{C}\cr}.]

An alternate way of indexing is based on a 3D set of cubic reciprocal-basis vectors [{\bf b}_{i}^{*}, i = 1, \ldots, 3] (setting 2) (Fig. 4.6.3.32)[link]: [\eqalign{\pmatrix{{\bf b}_{1}^{*}\cr {\bf b}_{2}^{*}\cr {\bf b}_{3}^{*}\cr} &= {\textstyle{1 \over 2}} \pmatrix{0 &\bar{1} &0 &0 &0 &1\cr 1 &0 &0 &\bar{1} &0 &0\cr 0 &0 &1 &0 &1 &0\cr}_{D} \pmatrix{{\bf a}_{1}^{*}\cr {\bf a}_{2}^{*}\cr {\bf a}_{3}^{*}\cr {\bf a}_{4}^{*}\cr {\bf a}_{5}^{*}\cr {\bf a}_{6}^{*}\cr}\cr &= {a^{*} \over (1 + \tau^{2})^{1/2}} \pmatrix{{\bf e}_{1}^{C}\cr {\bf e}_{2}^{C}\cr {\bf e}_{3}^{C}\cr}.}]The Cartesian C basis is related to the V basis by a [\theta /2] rotation around [[100]_{C}], yielding [[001]_{V}], followed by a [\pi /10] rotation around [[001]_{C}]: [\displaylines{ \pmatrix{{\bf e}_{1}^{C}\cr {\bf e}_{2}^{C}\cr {\bf e}_{3}^{C}\cr} = \pmatrix{\cos (\pi/10) &\sin (\pi/10) &0\cr - \cos (\theta/2)\sin (\pi/10) &\cos (\theta/2)\cos (\pi/10) &\sin (\theta/2)\cr \sin (\theta/2)\sin (\pi/10) &- \sin (\theta/2)\cos (\pi/10) &\cos (\theta/2)\cr}_{V} \pmatrix{{\bf e}_{1}^{V}\cr {\bf e}_{2}^{V}\cr {\bf e}_{3}^{V}\cr}.}]Thus, indexing the diffraction pattern of an icosahedral phase with integer indices, one obtains for setting 1 [{\bf H} = {\textstyle\sum_{i = 1}^{6}} h_{i}{\bf a}_{i}^{*}, h_{i} \in {\bb Z}]. These indices [(h_{1}\;h_{2}\;h_{3}\;h_{4}\;h_{5}\;h_{6})] transform into the second setting to [(h/h'\;k/k'\;l/l')] with the fractional cubic indices [h_{1}^{\rm c} = h + h' \tau], [h_{2}^{\rm c} = k + k'\tau], [h_{3}^{\rm c} = l + l' \tau]. The transformation matrix is [\pmatrix{h\cr h'\cr k\cr k'\cr l\cr l'\cr}_{C} = \pmatrix{0 &\bar{1} &0 &0 &0 &1\cr 0 &0 &\bar{1} &0 &1 &0\cr 1 &0 &0 &\bar{1} &0 &0\cr 0 &1 &0 &0 &0 &1\cr 0 &0 &1 &0 &1 &0\cr 1 &0 &0 &1 &0 &0\cr} \pmatrix{h_{1}\cr h_{2}\cr h_{3}\cr h_{4}\cr h_{5}\cr h_{6}\cr}_{D} = \pmatrix{h_{6} - h_{2}\cr h_{5} - h_{3}\cr h_{1} - h_{4}\cr h_{6} + h_{2}\cr h_{5} + h_{3}\cr h_{1} + h_{4}\cr}_{D}.]

[Figure 4.6.3.32]

Figure 4.6.3.32 | top | pdf |

Perspective parallel-space view of the two alternative reciprocal bases of the 3D Penrose tiling: the cubic and the icosahedral setting, represented by the bases [{\bf b}_{i}^{*}, i = 1,\ldots, 3], and [{\bf a}_{i}^{*}, i = 1,\ldots, 6], respectively.

4.6.3.3.3.2. Diffraction symmetry

| top | pdf |

The diffraction symmetry of icosahedral phases can be described by the Laue group [K = m\bar{3}\bar{5}]. The set of all vectors H forms a Fourier module [M^{*} = \{{\bf H}^{\parallel} = {\textstyle\sum_{i = 1}^{6}} h_{i}{\bf a}_{i}^{*} | h_{i} \in {\bb Z}\}] of rank 6 in physical space. Consequently, it can be considered as a projection from a 6D reciprocal lattice, [M^{*} = \pi^{\parallel}(\Sigma^{*})]. The parallel and perpendicular reciprocal-space sections of the 3D Penrose tiling decorated with equal point scatterers on its vertices are shown in Figs. 4.6.3.33[link] and 4.6.3.34[link]. The diffraction pattern in perpendicular space is the Fourier transform of the triacontahedron. All Bragg reflections within [10^{-4} |F({\bf 0})|^{2} \,\lt\, |F({\bf H})|^{2}\, \lt\, |F({\bf 0})|^{2}] are depicted. Without intensity-truncation limit, the diffraction pattern would be densely filled with discrete Bragg reflections.

[Figure 4.6.3.33]

Figure 4.6.3.33 | top | pdf |

Physical-space diffraction patterns of the 3D Penrose tiling decorated with point atoms (edge lengths of the Penrose unit rhombohedra ar = 5.0 Å). Sections with five-, three- and twofold symmetry are shown for the primitive 6D analogue of Bravais type P in (a, b, c), the body-centred 6D analogue to Bravais type I in (d, e, f) and the face-centred 6D analogue to Bravais type F in (g, h, i). All reflections are shown within [10^{-4} |F({\bf 0})|^{2} \;\lt\; |F({\bf H})|^{2} \;\lt\; |F({\bf 0})|^{2}] and [-6 \leq h_{i} \leq 6, i = 1,\ldots, 6].

[Figure 4.6.3.34]

Figure 4.6.3.34 | top | pdf |

Perpendicular-space diffraction patterns of the 3D Penrose tiling decorated with point atoms (edge lengths of the Penrose unit rhombohedra ar = 5.0 Å). Sections with (a) five-, (b) three- and (c) twofold symmetry are shown for the primitive 6D analogue of Bravais type P. All reflections are shown within [10^{-4} |F({\bf 0})|^{2} \;\lt\; |F({\bf H})|^{2} \;\lt\; |F({\bf 0})|^{2}] and [-6 \leq h_{i} \leq 6, i = 1,\ldots, 6].

The 6D icosahedral space groups that are relevant to the description of icosahedral phases (six symmorphous and five non-symmorphous groups) are listed in Table 4.6.3.2.[link] These space groups are a subset of all 6D icosahedral space groups fulfilling the condition that the 6D point groups they are associated with are isomorphous to the 3D point groups [{2 \over m}\bar{3}\bar{5}] and 235 describing the diffraction symmetry. From 826 6D (analogues to) Bravais groups (Levitov & Rhyner, 1988[link]), only three fulfil the condition that the projection of the 6D hypercubic lattice upon the 3D physical space is compatible with the icosahedral point groups [{2 \over m}\bar{3}\bar{5}, \;235]: the primitive hypercubic Bravais lattice P, the body-centred Bravais lattice I with translation 1/2(111111), and the face-centred Bravais lattice F with translations [1/2(110000) + 14] further cyclic permutations. Hence, the I lattice is twofold primitive (i.e. it contains two vertices per unit cell) and the F lattice is 16-fold primitive. The orientation of the symmetry elements in the 6D space is defined by the isomorphism of the 3D and 6D point groups. The action of the fivefold rotation, however, is different in the subspaces [{\bf V}^{\parallel}] and [{\bf V}^{\perp}]: a rotation of [2\pi/5] in [{\bf V}^{\parallel}] is correlated with a rotation of [4\pi/5] in [{\bf V}^{\perp}]. The reflection and inversion operations are equivalent in both subspaces.

Table 4.6.3.2| top | pdf |
3D point groups of order k describing the diffraction symmetry and corresponding 6D decagonal space groups with reflection conditions (see Levitov & Rhyner, 1988[link]; Rokhsar et al., 1988[link])

3D point groupk5D space groupReflection condition
[\displaystyle{{2 \over m}}\bar{3}\bar{5}] 120 [P\displaystyle{{2 \over m}}\bar{3}\bar{5}] No condition
[P\displaystyle{{2 \over n}}\bar{3}\bar{5}] [h_{1}h_{2}\overline{h_{1}h_{2}}h_{5}h_{6}]: [h_{5} - h_{6} = 2n]
[I\displaystyle{{2 \over m}}\bar{3}\bar{5}] [h_{1}h_{2}h_{3}h_{4}h_{5}h_{6}]: [{\textstyle\sum_{i = 1}^{6}} h_{i} = 2n]
[F\displaystyle{{2 \over m}}\bar{3}\bar{5}] [h_{1}h_{2}h_{3}h_{4}h_{5}h_{6}]: [{\textstyle\sum_{i \neq j = 1}^{6}} h_{i} + h_{j} = 2n]
[F\displaystyle{{2 \over n}}\bar{3}\bar{5}] [h_{1}h_{2}h_{3}h_{4}h_{5}h_{6}]: [{\textstyle\sum_{i \neq j = 1}^{6}} h_{i} + h_{j} = 2n]
[h_{1}h_{2}\overline{h_{1}h_{2}}h_{5}h_{6}]: [h_{5} - h_{6} = 4n]
235 60 P235 No condition
[P235_{1}] [h_{1}h_{2}h_{2}h_{2}h_{2}h_{2}]: [h_{1} = 5n]
I235 [h_{1}h_{2}h_{3}h_{4}h_{5}h_{6}]: [{\textstyle\sum_{i = 1}^{6}} h_{i} = 2n]
[I235_{1}] [h_{1}h_{2}h_{3}h_{4}h_{5}h_{6}]: [{\textstyle\sum_{i = 1}^{6}} h_{i} = 2n]
  [h_{1}h_{2}h_{2}h_{2}h_{2}h_{2}]: [h_{1} + 5h_{2} = 10n]
F235 [h_{1}h_{2}h_{3}h_{4}h_{5}h_{6}]: [{\textstyle\sum_{i \neq j = 1}^{6}} h_{i} + h_{j} = 2n]
[F235_{1}] [h_{1}h_{2}h_{3}h_{4}h_{5}h_{6}]: [{\textstyle\sum_{i \neq j = 1}^{6}} h_{i} + h_{j} = 2n]
  [h_{1}h_{2}h_{2}h_{2}h_{2}h_{2}]: [h_{1} + 5h_{2} = 10n]

4.6.3.3.3.3. Structure factor

| top | pdf |

The structure factor of the icosahedral phase corresponds to the Fourier transform of the 6D unit cell, [F({\bf H}) = {\textstyle\sum\limits_{k = 1}^{N}}\, f_{k} ({\bf H}^{\parallel})T_{k} ({\bf H}^{\parallel}, {\bf H}^{\perp})g_{k} ({\bf H}^{\perp}) \exp (2\pi i {\bf H}\cdot {\bf r}_{k}),]with 6D diffraction vectors [{\bf H} = {\textstyle\sum_{i = 1}^{6}} h_{i}{\bf d}_{i}^{*}], parallel-space atomic scattering factor [f_{k}(H^{\parallel})], temperature factor [T_{k}({\bf H}^{\parallel}, {\bf H}^{\perp})] and perpendicular-space geometric form factor [g_{k}({\bf H}^{\perp})]. [T_{k}({\bf H}^{\parallel}, {\bf 0})] is equivalent to the conventional Debye–Waller factor and [T_{k}({\bf 0},{\bf H}^{\perp})] describes random fluctuations in perpendicular space. These fluctuations cause characteristic jumps of vertices (phason flips) in the physical space. Even random phason flips map the vertices onto positions that can still be described by physical-space vectors of the type [{\bf r} = {\textstyle\sum_{i = 1}^{6}} n_{i}{\bf a}_{i}]. Consequently, the set [M = \{{\bf r} = {\textstyle\sum_{i = 1}^{6}} n_{i}{\bf a}_{i} \big| n_{i} \in {\bb Z}\}] of all possible vectors forms a [{\bb Z}] module. The shape of the atomic surfaces corresponds to a selection rule for the positions actually occupied. The geometric form factor [g_{k}({\bf H}^{\perp})] is equivalent to the Fourier transform of the atomic surface, i.e. the 3D perpendicular-space component of the 6D hyperatoms.

For the example of the canonical 3D Penrose tiling, [g_{k}({\bf H}^{\perp})] corresponds to the Fourier transform of a triacontahedron: [g_{k}({\bf H}^{\perp}) = (1/A_{\rm UC}^{\perp}) {\textstyle\int\limits_{A_{k}}} \exp (2\pi i {\bf H}^{\perp}\cdot {\bf r}) \hbox{ d}{\bf r},]where [A_{\rm UC}^{\perp}] is the volume of the 6D unit cell projected upon [{\bf V}^{\perp}] and [A_{k}] is the volume of the triacontahedron. [A_{\rm UC}^{\perp}] and [A_{k}] are equal in the present case and amount to the volumes of ten prolate and ten oblate rhombohedra: [A_{\rm UC}^{\perp} = 8a_{\rm r}^{3} \left[\sin (2\pi /5) + \sin (\pi /5)\right]]. Evaluating the integral by decomposing the triacontahedron into trigonal pyramids, each one directed from the centre of the triacontahedron to three of its corners given by the vectors [{\bf e}_{i}, i = 1, \ldots, 3], one obtains [g({\bf H}^{\perp}) = (1/A_{\rm UC}^{\perp}) {\textstyle\sum\limits_{R}} g_{k}(R^{T}{\bf H}^{\perp}),]with [k = 1, \ldots, 60] running over all site-symmetry operations R of the icosahedral group, [\eqalign{g_{k} ({\bf H}^{\perp}) &= -iV_{\rm r} [A_{2}A_{3}A_{4} \exp (iA_{1}) + A_{1}A_{3}A_{5} \exp (iA_{2}) \cr&\quad+ A_{1}A_{2}A_{6} \exp (iA_{3}) + A_{4}A_{5}A_{6}]\cr&\quad\times(A_{1}A_{2}A_{3}A_{4}A_{5}A_{6})^{-1},}][A_{j} = 2\pi {\bf H}^{\perp} \cdot {\bf e}_{j}], [ j = 1, \ldots, 3], [A_{4} = A_{2} - A_{3}], [A_{5} = A_{3} - A_{1}], [A_{6} = A_{1} - A_{2}] and [V_{\rm r} = {\bf e}_{1} \cdot ({\bf e}_{2} \times {\bf e}_{3})] the volume of the parallelepiped defined by the vectors [{\bf e}_{i}, i = 1, \ldots, 3] (Yamamoto, 1992b[link]).

4.6.3.3.3.4. Intensity statistics

| top | pdf |

The radial structure-factor distributions of the 3D Penrose tiling decorated with point scatterers are plotted in Fig. 4.6.3.35[link] as a function of parallel and perpendicular space. The distribution of [|F({\bf H})|] as a function of their frequencies clearly resembles a centric distribution, as can be expected from the centrosymmetric unit cell. The shape of the distribution function depends on the radius of the limiting sphere in reciprocal space. The number of weak reflections increases as the power 6, that of strong reflections only as the power 3 (strong reflections always have small [{\bf H}^{\perp}] components).

[Figure 4.6.3.35]

Figure 4.6.3.35 | top | pdf |

Radial distribution function of the structure factors [F({\bf H})] of the 3D Penrose tiling (edge lengths of the Penrose unit rhombohedra ar = 5.0 Å) decorated with point atoms as a function of [|{\bf H}^{\parallel}|] (above) and [|{\bf H}^{\perp}|] (below). All reflections are shown within [10^{-6} |F({\bf 0})|^{2} \;\lt\; |F({\bf H})|^{2} \;\lt\; |F({\bf 0})|^{2}] and [-6 \leq h_{i} \leq 6, i = 1,\ldots, 6].

The weighted reciprocal space of the 3D Penrose tiling contains an infinite number of Bragg reflections within a limited region of the physical space. Contrary to the diffraction pattern of a periodic structure consisting of point atoms on the lattice nodes, the Bragg reflections show intensities depending on the perpendicular-space components of their diffraction vectors.

4.6.3.3.3.5. Relationships between structure factors at symmetry-related points of the Fourier image

| top | pdf |

The weighted 3D reciprocal space [M^{*} = \{{\bf H}^{\parallel} =] [{\textstyle\sum_{i = 1}^{6}} h_{i} {\bf a}_{i}^{*} | h_{i} \in {\bb Z}\}] exhibits the icosahedral point symmetry [K = m \bar{3} \bar{5}]. It is invariant under the action of the scaling matrix [S^{3}]: [\displaylines{S = {\textstyle{1 \over 2}} \pmatrix{1\hfill &1\hfill &1\hfill &1\hfill &1\hfill &1\hfill\cr 1\hfill &1\hfill &1\hfill &-1\hfill &-1\hfill &1\hfill\cr 1\hfill &1\hfill &1\hfill &1\hfill &-1\hfill &-1\hfill\cr 1\hfill &-1\hfill &1\hfill &1\hfill &1\hfill &-1\hfill\cr 1\hfill &-1\hfill &-1\hfill &1\hfill &1\hfill &1\hfill\cr 1\hfill &1\hfill &-1\hfill &-1\hfill &1\hfill &1\hfill\cr}_{D}, S^{3} = \pmatrix{2\hfill &1\hfill &1\hfill &1\hfill &1\hfill &1\hfill\cr 1\hfill &2\hfill &1\hfill &-1\hfill &-1\hfill &1\hfill\cr 1\hfill &1\hfill &2\hfill &1\hfill &-1\hfill &-1\hfill\cr 1\hfill &-1\hfill &1\hfill &2\hfill &1\hfill &-1\hfill\cr 1\hfill &-1\hfill &-1\hfill &1\hfill &2\hfill &1\hfill\cr 1\hfill &1\hfill &-1\hfill &-1\hfill &1\hfill &2\hfill\cr}_{D},\cr \pmatrix{2\hfill &1\hfill &1\hfill &1\hfill &1\hfill &1\hfill\cr 1\hfill &2\hfill &1\hfill &-1\hfill &-1\hfill &1\hfill\cr 1\hfill &1\hfill &2\hfill &1\hfill &-1\hfill &-1\hfill\cr 1\hfill &-1\hfill &1\hfill &2\hfill &1\hfill &-1\hfill\cr 1\hfill &-1\hfill &-1\hfill &1\hfill &2\hfill &1\hfill\cr 1\hfill &1\hfill &-1\hfill &-1\hfill &1\hfill &2\hfill\cr}_{D} \pmatrix{{\bf a}_{1}^{*}\cr {\bf a}_{2}^{*}\cr {\bf a}_{3}^{*}\cr {\bf a}_{4}^{*}\cr {\bf a}_{5}^{*}\cr {\bf a}_{6}^{*}\cr} = \tau^{3} \pmatrix{{\bf a}_{1}^{*}\cr {\bf a}_{2}^{*}\cr {\bf a}_{3}^{*}\cr {\bf a}_{4}^{*}\cr {\bf a}_{5}^{*}\cr {\bf a}_{6}^{*}\cr}.}]The scaling transformation [(S^{-3})^{T}] leaves a primitive 6D reciprocal lattice invariant as can easily be seen from its application on the indices: [\openup2pt{\pmatrix{h^{'}_{1}\cr h^{'}_{2}\cr h^{'}_{3}\cr h^{'}_{4}\cr h^{'}_{5}\cr h^{'}_{6}\cr}} = \pmatrix{-2\hfill &1\hfill &1\hfill &1\hfill &1\hfill &1\hfill\cr 1\hfill &-2\hfill &1\hfill &-1\hfill &-1\hfill &1\hfill\cr 1\hfill &1\hfill &-2\hfill &1\hfill &-1\hfill &-1\hfill\cr 1\hfill &-1\hfill &1\hfill &-2\hfill &1\hfill &-1\hfill\cr 1\hfill &-1\hfill &-1\hfill &1\hfill &-2\hfill &1\hfill\cr 1\hfill &1\hfill &-1\hfill &-1\hfill &1\hfill &-2\hfill\cr}_{D} \pmatrix{h_{1}\cr h_{2}\cr h_{3}\cr h_{4}\cr h_{5}\cr h_{6}\cr}.]The matrix [(S^{-1})^{T}] leaves [M^{*} = \{{\bf H}^{\parallel} = {\textstyle\sum_{i = 1}^{6}} h_{i} {\bf a}_{i}^{*} | h_{i} \in {\bb Z}\}] invariant, [\openup2pt\pmatrix{h^{'}_{1}\cr h^{'}_{2}\cr h^{'}_{3}\cr h^{'}_{4}\cr h^{'}_{5}\cr h^{'}_{6}\cr} = {\textstyle{1 \over 2}} \pmatrix{-1\hfill &1\hfill &1\hfill &1\hfill &1\hfill &1\hfill\cr 1\hfill &-1\hfill &1\hfill &-1\hfill &-1\hfill &1\hfill\cr 1\hfill &1\hfill &-1\hfill &1\hfill &-1\hfill &-1\hfill\cr 1\hfill &-1\hfill &1\hfill &-1\hfill &1\hfill &-1\hfill\cr 1\hfill &-1\hfill &-1\hfill &1\hfill &-1\hfill &1\hfill\cr 1\hfill &1\hfill &-1\hfill &-1\hfill &1\hfill &-1\hfill\cr}_{D} \pmatrix{h_{1}\cr h_{2}\cr h_{3}\cr h_{4}\cr h_{5}\cr h_{6}\cr},]for any [{\bf H} = {\textstyle\sum_{i = 1}^{6}} h_{i} {\bf d}_{i}^{*}] with [h_{i}] all even or all odd, corresponding to a 6D face-centred hypercubic lattice. In a second case the sum [{\textstyle\sum_{i = 1}^{6}} h_{i}] is even, corresponding to a 6D body-centred hypercubic lattice. Block-diagonalization of the matrix S decomposes it into two irreducible representations. With [WSW^{-1} = S_{V} = S_{V}^{\parallel} \oplus S_{V}^{\perp}] we obtain [S_{V} = \pmatrix{\tau &0 &0{\,\vrule height 8pt depth 8pt} &0 &0 &0\cr\noalign{\vskip -5pt} 0 &\tau &0{\,\vrule height 8pt depth 8pt} &0 &0 &0\cr\noalign{\vskip -5pt} 0 &0 &\tau{\hskip 1.7pt}{\vrule height 8pt depth 8pt} &0 &0 &0\cr\noalign{\vskip -4pt}\noalign{\hrule}\noalign{\vskip 2pt} 0 &0 &0{\,\vrule height 8pt depth 8pt} &-1/\tau &0 &0\cr\noalign{\vskip -5pt} 0 &0 &0{\,\vrule height 8pt depth 8pt} &0 &-1/\tau &0\cr\noalign{\vskip -5pt} 0 &0 &0{\,\vrule height 8pt depth 2pt} &0 &0 &-1/\tau\cr}_{V} = \pmatrix{S^{\parallel}{\hskip 0.7pt}{\vrule height 8pt depth 6pt} &{\hskip -6pt}0\cr\noalign{\vskip -2pt}\noalign{\hrule}\noalign{\vskip 1pt}0{\hskip 5.0pt}{\vrule height 10pt depth 2pt} &S^{\perp}\cr}_{V},]the scaling properties in the two 3D subspaces: scaling by a factor τ in parallel space corresponds to a scaling by a factor [(-\tau)^{-1}] in perpendicular space. For the intensities of the scaled reflections analogous relationships are valid, as discussed for decagonal phases (Figs. 4.6.3.36[link] and 4.6.3.37[link], Section 4.6.3.3.2.5[link]).

[Figure 4.6.3.36]

Figure 4.6.3.36 | top | pdf |

Parallel-space distribution of (a) positive and (b) negative structure factors of the 3D Penrose tiling of the 6D P lattice type decorated with point atoms (edge lengths of the Penrose unit rhombohedra ar = 5.0 Å). The magnitudes of the structure factors are indicated by the diameters of the filled circles. All reflections are shown within [10^{-4} |F({\bf 0})|^{2} \;\lt\; |F({\bf H})|^{2} \;\lt\; |F({\bf 0})|^{2}] and [-6 \leq h_{i} \leq 6, i = 1,\ldots, 6].

[Figure 4.6.3.37]

Figure 4.6.3.37 | top | pdf |

Perpendicular-space distribution of (a) positive and (b) negative structure factors of the 3D Penrose tiling of the 6D P lattice type decorated with point atoms (edge lengths of the Penrose unit rhombohedra ar = 5.0 Å). The magnitudes of the structure factors are indicated by the diameters of the filled circles. All reflections are shown within [10^{-4} |F({\bf 0})|^{2} \;\lt\; |F({\bf H})|^{2} \;\lt\; |F({\bf 0})|^{2}] and [-6 \leq h_{i} \leq 6, i = 1,\ldots, 6].

4.6.4. Experimental aspects of the reciprocal-space analysis of aperiodic crystals

| top | pdf |

4.6.4.1. Data-collection strategies

| top | pdf |

Theoretically, aperiodic crystals show an infinite number of reflections within a given diffraction angle, contrary to periodic crystals. The number of reflections to be included in a structure analysis of a periodic crystal may be very high (one million for virus crystals, for instance) but there is no ambiguity in the selection of reflections to be collected: all Bragg reflections within a limiting sphere in reciprocal space, usually given by [0 \leq \sin \theta/\lambda \leq 0.7\;\hbox{\AA}^{-1}], are used. All reflections, observed and unobserved, are included to fit a reliable structure model.

However, for aperiodic crystals it is not possible to collect the infinite number of dense Bragg reflections within [0 \leq \sin \theta/\lambda \leq 0.7\;\hbox{\AA}^{-1}]. The number of observable reflections within this limiting sphere depends only on the spatial and intensity resolution.

What happens if not all reflections are included in a structure analysis? How important is the contribution of reflections with large perpendicular-space components of the diffraction vector which are weak but densely distributed? These problems are illustrated using the example of the Fibonacci sequence. An infinite model structure consisting of Al atoms with isotropic thermal parameter B = 1 Å2, and distances S = 2.5 Å and [{\rm L} = \tau{\rm S}], was used for the calculations (Table 4.6.4.1)[link].

Table 4.6.4.1| top | pdf |
Intensity statistics of the Fibonacci chain for a total of 161 322 reflections with [-1000 \leq h_{i} \leq 1000] and [0 \leq \sin \theta/\lambda \leq 2\;\hbox{\AA}^{-1}]

In the upper line, the number of reflections in the respective interval is given; in the lower line the partial sums [ {\textstyle\sum} {I}({\bf H})] of the intensities [I({\bf H})] are given as a percentage of the total diffracted intensity. The F(00) reflection is not included in the sums.

 [F({\bf H})/F({\bf H})_{\max} \geq 0.1][0.1\, \gt\, F({\bf H})/F({\bf H})_{\max} \geq 0.01][0.01 \,\gt\, F({\bf H})/F({\bf H})_{\max} \geq 0.001][F({\bf H})/F({\bf H})_{\max} \,\lt\, 0.001]
[0 \leq \sin \theta/\lambda \leq 0.2\;{\AA}^{-1}] 17 148 1505 14 511
[ {\textstyle\sum} {I}({\bf H})] 52.53% 2.56% 0.27% 0.03%
[0.2 \leq \sin \theta/\lambda \leq 0.4\;\hbox{\AA}^{-1}] 11 107 1066 14 998
[ {\textstyle\sum} {I}({\bf H})] 27.03% 2.03% 0.19% 0.02%
[0.4 \leq \sin \theta/\lambda \leq 0.6\;\hbox{\AA}^{-1}] 9 64 654 15 456
[ {\textstyle\sum} {I}({\bf H})] 9.84% 0.96% 0.12% 0.01%
[0.6 \leq \sin \theta/\lambda \leq 0.8\;\hbox{\AA}^{-1}] 6 27 326 15 823
[ {\textstyle\sum} {I}({\bf H})] 2.94% 0.34% 0.07% 0.01%
[0.8 \leq \sin \theta/\lambda \leq 2\;\hbox{\AA}^{-1}] 1 35 338 96 720
[ {\textstyle\sum} {I}({\bf H})] 0.23% 0.79% 0.06% 0.01%
Total sum 44 381 3389 157 508
92.57% 6.67% 0.70% 0.06%

It turns out that 92.6% of the total diffracted intensity of 161 322 reflections is included in the 44 strongest reflections and 99.2% in the strongest 425 reflections. It is remarkable, however, that in all the experimental data for icosahedral and decagonal quasicrystals collected so far, rarely more than 20 to 50 reflections along reciprocal-lattice lines corresponding to net planes with Fibonacci-sequence-like distances could be observed. The consequences for structure determinations with such truncated data sets are primarily a lower resolution in perpendicular space than in physical space. This corresponds to a smearing of the hyperatoms in the perpendicular space. For the derivation of the local structure-building elements (clusters) of aperiodic crystals this is only a minor problem: the smeared hyperatoms give rise to split atoms and a biased electron-density distribution. The information on the global aperiodic structure, however, which is contained in the detailed shape of the atomic surfaces, is severely reduced when using a low-resolution diffraction data set. A combination of high-resolution electron microscopy, lattice imaging and diffraction techniques allows a good characterization of the local and global order even in these cases. For a more detailed analysis of these problems see Steurer (1995)[link].

4.6.4.2. Commensurability versus incommensurability

| top | pdf |

The question whether an aperiodic crystal is really aperiodic or rather a high-order approximant is of different importance depending on the point of view. As far as real finite crystals are considered, definitions of periodic and aperiodic real crystals and of periodic and aperiodic perfect crystals have to be given first. Real crystals, despite periodicity or aperiodicity, are the actual samples under investigation. Partial information about their actual structure can be obtained today by imaging methods (scanning tunnelling microscopy, atomic force microscopy, high-resolution transmission electron microscopy, …). Basically, the real crystal structure can be determined using full diffraction information from Bragg and diffuse scattering. In practice, however, only `Bragg reflections' are included in a structure analysis. `Bragg reflections' result from the integration of diffraction intensities from extended volumes around a limited number of Bragg-reflection positions ([{\bb Z}] module). This process of intensity condensation at Bragg points corresponds in direct space to an averaging process. The real crystal structure is projected upon one unit cell in direct space defined by the [{\bb Z}] module in reciprocal space. Generally, the identification of appropriate reciprocal-space metrics is not a problem in the case of crystals. It can be problematic, however, in the case of aperiodic crystals, in particular quasicrystals (see Lancon et al., 1994[link]). The metrics, and to some extent the global order in the case of quasicrystals, are fixed by assigning the reciprocal basis. The spatial resolution of a diffraction experiment defines the accuracy of the resulting metrics. The decision whether the rational number obtained for the relative length of a satellite vector indicates a commensurate or an incommensurate modulation can only be made considering temperature- and pressure-dependent chemical and physical properties of the material. The same is valid for other types of aperiodic crystals.

4.6.4.3. Twinning and nanodomain structures

| top | pdf |

High-order approximants of quasicrystals often occur in orientationally twinned form or, on a smaller scale, as oriented nanodomain structures. These structures can be identified by electron microscopy, and, in certain cases, also by high-resolution X-ray diffractometry (Kalning et al., 1994[link]). If the intensity and spatial resolution is sufficient, characteristic reflection splitting and diffuse diffraction phenomena can be observed. It has been demonstrated that for the determination of the local structure (structure-building elements) it does not matter greatly whether one uses a data set for a real quasicrystal or one for a twinned approximant (Estermann et al., 1994[link]). Examples of reciprocal-space images of an approximant, a twinned approximant and the related decagonal phase are shown schematically in Fig. 4.6.4.1[link] and for real samples in Fig. 4.6.4.2[link].

[Figure 4.6.4.1]

Figure 4.6.4.1 | top | pdf |

Simulated diffraction patterns of (a) the [\simeq 52\;\hbox{\AA}] single-crystal approximant of decagonal Al–Co–Ni, (b) the fivefold twinned approximant, and (c) the decagonal phase itself (Estermann et al., 1994[link]).

[Figure 4.6.4.2]

Figure 4.6.4.2 | top | pdf |

Zero-layer X-ray diffraction patterns of decaprismatic Al73.5C21.7Ni4.8 crystals taken parallel and perpendicular to the crystal axis on an image-plate scanner (Mar Research) at different temperatures. In (a) and (b), room-temperature (RT) diffraction patterns from a sample quenched after annealing at 1073 K are shown. Reflections from both a crystalline approximant and a decagonal phase are visible. The period along the unique direction in the decagonal phase and the corresponding period in the approximant phase is [\simeq 8\;\hbox{\AA}] (b). At 1520 K, a single-phase decagonal quasicrystal is present with [\simeq 4\;\hbox{\AA}] fundamental structure (c, d). In (e, f), the RT diffraction patterns of the slowly cooled sample indicate a single-phase nanodomain structure with [\simeq 8\;\hbox{\AA}] periodicity along the unique direction.

References

Axel, F. & Gratias, D. (1995). Editors. Beyond Quasicrystals. Les Ulis: Les Editions de Physique and Berlin: Springer-Verlag.
Bancel, P. A., Heiney, P. A., Stephens, P. W., Goldman, A. I. & Horn, P. M. (1985). Structure of rapidly quenched Al–Mn. Phys. Rev. Lett. 54, 2422–2425.
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.
Dehlinger, U. (1927). Über die Verbreiterung der Debyelinien bei kaltbearbeiteten Metallen. Z. Kristallogr. 65, 615–631.
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.
Estermann, M., Haibach, T. & Steurer, W. (1994). Quasicrystal versus twinned approximant: a quantitative analysis with decagonal Al70Co15Ni15. Philos. Mag. Lett. 70, 379–384.
Goldman, A. I. & Kelton, K. F. (1993). Quasicrystals and crystalline approximants. Rev. Mod. Phys. 65, 213–230.
Gouyet, J. F. (1996). Physics and Fractal Structures. Paris: Masson and Berlin: Springer-Verlag.
Hausdorff, F. (1919). Dimension und äusseres Mass. Math. Ann. 79, 157–179.
Hermann, C. (1949). Kristallographie in Räumen beliebiger Dimensionszahl. I. Die Symmetrieoperationen. Acta Cryst. 2, 139–145.
Ishihara, K. N. & Yamamoto, A. (1988). Penrose patterns and related structures. I. Superstructure and generalized Penrose patterns. Acta Cryst. A44, 508–516.
Janner, A. (1992). Decagrammal symmetry of decagonal Al78Mn22 quasicrystal. Acta Cryst. A48, 884–901.
Janner, A. & Janssen, T. (1979). Superspace groups. Physica A, 99, 47–76.
Janner, A. & Janssen, T. (1980a). Symmetry of incommensurate crystal phases. I. Commensurate basic structures. Acta Cryst. A36, 399–408.
Janner, A. & Janssen, T. (1980b). Symmetry of incommensurate crystal phases. II. Incommensurate basic structures. Acta Cryst. A36, 408–415.
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.
Janot, Chr. (1994). Quasicrystals. A Primer. Oxford: Clarendon Press.
Janssen, T. (1986). Crystallography of quasicrystals. Acta Cryst. A42, 261–271.
Janssen, T. (1988). Aperiodic crystals: a contradictio in terminis? Phys. Rep. 168, 55–113.
Janssen, T. (1995). From quasiperiodic to more complex systems. In Beyond Quasicrystals, edited by F. Axel & D. Gratias, pp. 75–140. Les Ulis: Les Editions de Physique and Berlin: Springer-Verlag.
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.
Jaric, M. V. (1986). Diffraction from quasicrystals: geometric structure factor. Phys. Rev. B, 34, 4685–4698.
Kalning, M., Kek, S., Burandt, B., Press, W. & Steurer, W. (1994). Examination of a multiple twinned periodic approximant of the decagonal phase Al70Co15Ni15. J. Phys. Condens. Matter, 6, 6177–6187.
Kato, K. (1990). Strukturverfeinerung des Kompositkristalls im mehr­dimensionalen Raum. Acta Cryst. B46, 39–44.
Kelton, K. F. (1995). Quasicrystals and related structures. In Intermetallic Compounds. Principles and Practice, Vol. 1, edited by J. H. Westbrook & R. L. Fleischer, pp. 453–491. Chichester: John Wiley & Sons.
Korekawa, M. (1967). Theorie der Satellitenreflexe. Habilitation thesis, University of Munich.
Lancon, F., Billard, L., Burkov, S. & DeBoissieu, M. (1994). On choosing a proper basis for determining structures of quasicrystals. J. Phys. I France, 4, 283–301.
Levine, D. & Steinhardt, P. J. (1986). Quasicrystals. I. Definition and structure. Phys. Rev. B, 34, 596–616.
Levitov, L. S. & Rhyner, J. (1988). Crystallography of quasicrystals; application to icosahedral symmetry. J. Phys. France, 49, 1835–1849.
Luck, J. M., Godréche, C., Janner, A. & Janssen, T. (1993). The nature of the atomic surfaces of quasiperiodic self-similar structures. J. Phys. A Math. Gen. 26, 1951–1999.
Paciorek, W. A. & Chapuis, G. (1994). Generalized Bessel functions in incommensurate structure analysis. Acta Cryst. A50, 194–203.
Pavlovitch, A. & Kleman, M. (1987). Generalized 2D Penrose tilings: structural properties. J. Phys. A Math. Gen. 20, 687–702.
Penrose, R. (1974). The role of aesthetics in pure and applied mathematical research. Bull. Math. Appl. 10, 266–271.
Penrose, R. (1979). Pentaplexity. A class of non-periodic tilings of the plane. Math. Intell. 2, 32–37.
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.
Petricek, V., Maly, K., Coppens, P., Bu, X., Cisarova, I. & Frost-Jensen, A. (1991). The description and analysis of composite crystals. Acta Cryst. A47, 210–216.
Rabson, D. A., Mermin, N. D., Rokhsar, D. S. & Wright, D. C. (1991). The space groups of axial crystals and quasicrystals. Rev. Mod. Phys. 63, 699–733.
Rokhsar, D. S., Wright, D. C. & Mermin, N. D. (1988). Scale equivalence of quasicrystallographic space groups. Phys. Rev. B, 37, 8145–8149.
Senechal, M. (1995). Quasicrystals and Geometry. Cambridge University Press.
Shechtman, D., Blech, I., Gratias, D. & Cahn, J. W. (1984). Metallic phase with long-range orientational order and no translational symmetry. Phys. Rev. Lett. 53, 1951–1953.
Smaalen, S. van (1991). Symmetry of composite crystals. Phys. Rev. B, 43, 11330–11341.
Smaalen, S. van (1992). Superspace description of incommensurate intergrowth compounds and the application to inorganic misfit layer compounds. Mater. Sci. Forum, 100 & 101, 173–222.
Smaalen, S. van (1995). Incommensurate crystal structures. Crystallogr. Rev. 4, 79–202.
Socolar, J. E. S. & Steinhardt, P. J. (1986). Quasicrystals. II. Unit-cell configurations. Phys. Rev. B, 34, 617–647.
Steurer, W. (1990). The structure of quasicrystals. Z. Kristallogr. 190, 179–234.
Steurer, W. (1995). Experimental aspects of the structure analysis of aperiodic materials. In Beyond Quasicrystals, edited by F. Axel & D. Gratias, pp. 203–228. Les Ulis: Les Editions de Physique and Berlin: Springer-Verlag.
Steurer, W. (1996). The structure of quasicrystals. In Physical Metallurgy, Vol. I, edited by R. W. Cahn & P. Haasen, pp. 371–411. Amsterdam: Elsevier.
Willis, B. T. M. & Pryor, A. W. (1975). Thermal Vibrations in Crystallography. Cambridge University Press.
Wolff, P. M. de (1974). The pseudo-symmetry of modulated crystal structures. Acta Cryst. A30, 777–785.
Wolff, P. M. de (1977). Symmetry operations for displacively modulated structures. Acta Cryst. A33, 493–497.
Wolff, P. M. de (1984). Dualistic interpretation of the symmetry of incommensurate structures. Acta Cryst. A40, 34–42.
Wolff, P. M. de, Janssen, T. & Janner, A. (1981). The superspace groups for incommensurate crystal structures with a one-dimensional modulation. Acta Cryst. A37, 625–636.
Yamamoto, A. (1982). Structure factor of modulated crystal structures. Acta Cryst. A38, 87–92.
Yamamoto, A. (1992a). Unified setting and symbols of superspace groups for composite crystals. Acta Cryst. A48, 476–483.
Yamamoto, A. (1992b). Ideal structure of icosahedral Al–Cu–Li quasicrystals. Phys. Rev. B, 45, 5217–5227.
Zobetz, E. (1993). One-dimensional quasilattices: fractally shaped atomic surfaces and homometry. Acta Cryst. A49, 667–676.








































to end of page
to top of page