International
Tables for
Crystallography
Volume D
Physical properties of crystals
Edited by A. Authier

International Tables for Crystallography (2006). Vol. D. ch. 1.10, pp. 243-264
https://doi.org/10.1107/97809553602060000637

Chapter 1.10. Tensors in quasiperiodic structures

T. Janssena*

a Institute for Theoretical Physics, University of Nijmegen, 6524 ED Nijmegen, The Netherlands
Correspondence e-mail: ted@sci.kun.nl

This chapter is devoted to the symmetry-related physical properties of quasiperiodic crystals. In the first part the symmetry properties are described: point groups, superspace groups and the action of symmetry groups. The second part concerns the properties of tensors in higher-dimensional spaces, with emphasis on the particular cases of the piezoelectric, elastic and electric field gradient tensors. The last section gives tables of characters of some point groups for quasicrystals and of the matrices of the corresponding irreducible representations.

Keywords: Fourier module; compensating gauge transformations; elastic constants; electric field gradient; icosahedral quasicrystals; incommensurate structures; irreducible representations; metric tensor; modulation wavevector; octagonal quasicrystals; piezoelectric tensor; quasicrystals; quasiperiodic structures; superspace; superspace groups; tensors.

1.10.1. Quasiperiodic structures

| top | pdf |

1.10.1.1. Introduction

| top | pdf |

Many materials are known which show a well ordered state without lattice translation symmetry, often in a restricted temperature or composition range. This can be seen in the diffraction pattern from the appearance of sharp spots that cannot be labelled in the usual way with three integer indices. The widths of the peaks are comparable with those of perfect lattice periodic crystals, and this is a sign that the coherence length is comparable as well.

A typical example is K2SeO4, which has a normal lattice periodic structure above 128 K with space group Pcmn, but below this temperature shows satellites at positions [\gamma {\bf c}^{*}], where [\gamma] is an irrational number, which in addition depends on temperature. These satellites cannot be labelled with integer indices with respect to the reciprocal basis [{\bf a}^{*}], [{\bf b}^{*}], [{\bf c}^{*}] of the structure above the transition temperature. Therefore, the corresponding structure cannot be lattice periodic.

The diffraction pattern of K2SeO4 arises because the original lattice periodic basic structure is deformed below 128 K. The atoms are displaced from their positions in the basic structure such that the displacement itself is again periodic, but with a period that is incommensurate with respect to the lattice of the basic structure.

Such a modulated structure is just a special case of a more general type of structure. These structures are characterized by the fact that the diffraction pattern has sharp Bragg peaks at positions [{\bf H}] that are linear combinations of a finite number of basic vectors: [{\bf H} = \textstyle\sum\limits_{i=1}^{n}h_{i}{\bf a}_{i}^{*} \quad({\rm integer} \;\; h_{i}).\eqno(1.10.1.1)]Structures that have this property are called quasiperiodic. The minimal number n of basis vectors such that all [h_{i}] are integers is called the rank of the structure. If the rank is three and the vectors [{\bf a}_i] do not all fall on a line or in a plane, the structure is just lattice periodic. Lattice periodic structures form special cases of quasiperiodic structures. The collection of vectors [{\bf H}] forms the Fourier module of the structure. For rank three, this is just the reciprocal lattice of the lattice periodic structure.

The definition given above results in some important practical difficulties. In the first place, it is not possible to show experimentally that a wavevector has irrational components instead of rational ones, because an irrational number can be approximated by a rational number arbitrarily well. Very often the wavevector of the satellite changes with temperature. It has been reported that in some compounds the variation shows plateaux, but even when the change seems to be continuous and smooth one can not be sure about the irrationality. On the other hand, if the wavevector jumps from one rational position to another, the structure would always be lattice periodic, but the unit cell of this structure would vary wildly with temperature. This means that, if one wishes to describe the incommensurate phases in a unified fashion, it is more convenient to treat the wavevector as generically irrational. This experimental situation is by no means dramatic. It is similar to the way in which one can never be sure that the angles between the basis vectors of an orthorhombic lattice are really 90°, although this is a concept that no-one has problems understanding.

A second problem stems from the fact that the wavevectors of the Fourier module are dense. For example, in the case of K2SeO4 the linear combinations of [{\bf c}^{*}] and [\gamma{\bf c}^{*}] cover the c axis uniformly. To pick out a basis here could be problematic, but the intensity of the spots is usually such that choosing a basis is not a problem. In fact, one only observes peaks with an intensity above a certain threshold, and these form a discrete set. At most, the occurrence of scale symmetry may make the choice less obvious.

1.10.1.2. Types of quasiperiodic crystals

| top | pdf |

One may distinguish various families of quasiperiodic systems. [Sometimes these are also called incommensurate systems (Janssen & Janner, 1987[link]).] It is not a strict classification, because one may have intermediate cases belonging to more than one family as well. Here we shall consider a number of pure cases.

An incommensurately modulated structure or incommensurate crystal (IC) phase is a periodically modified structure that without the modification would be lattice periodic. Hence there is a basic structure with space-group symmetry. The periodicity of the modification should be incommensurate with respect to the basic structure. The position of the jth atom in the unit cell with origin at the lattice point [{\bf n}] is [{\bf n} + {\bf r}_{j}] ([j = 1, 2, \ldots, s]).

For a displacive modulation, the positions of the atoms are shifted from a lattice periodic basic structure. A simple example is a structure that can be derived from the positions of the basic structure with a simple displacement wave. The positions of the atoms in the IC phase are then [{\bf n} + {\bf r}_{j} + {\bf f}_{j}({\bf Q}\cdot{\bf n})\quad [{\bf f}_{j}(x) = {\bf f}_{j}(x + 1)]. \eqno(1.10.1.2)]Here the modulation wavevector [{\bf Q}] has irrational components with respect to the reciprocal lattice of the basic structure. One has[{\bf Q} = \alpha{\bf a}^{*} + \beta{\bf b}^{*} + \gamma{\bf c}^{*},\eqno(1.10.1.3)]where at least one of α, β or γ is irrational. A simple example is the function [{\bf f}_{j}(x) = {\bf A}_{j}\cos(2\pi{x} + \varphi_{j})], where [{\bf A}_{j}] is the polarization vector and [\varphi_{j}] is the phase of the modulation. The diffraction pattern of the structure (1.10.1.2)[link] shows spots at positions[{\bf H} = h_{1}{\bf a}^{*} + h_{2}{\bf b}^{*} + h_{3}{\bf c}^{*} + h_{4}{\bf Q}. \eqno(1.10.1.4)]Therefore, the rank is four and [{\bf a}_{4}^{*} = {\bf Q}]. In a more general situation, the components of the atom positions in the IC phase are given by[{\bf n}^{\alpha} + {\bf r}_{j}^{\alpha} + \textstyle\sum\limits_{m}{\bf A}_{j}^{\alpha} ({\bf Q}_{m})\cos(2\pi {\bf Q}_{m}\cdot{\bf n} + \varphi_{jm\alpha}), \quad\alpha = x, y, z. \eqno(1.10.1.5)]Here the vectors [{\bf Q}_{m}] belong to the Fourier module of the structure. Then there are vectors [{\bf Q}_{j}] such that any spot in the diffraction pattern can be written as[{\bf H} = \textstyle\sum\limits_{i=1}^{3}h_{i}{\bf a}_{i}^{*} + \textstyle\sum\limits_{j=1}^{d} h_{3+j}{\bf Q}_{j}\eqno(1.10.1.6)]and the rank is [3+d]. The peaks corresponding to the basic structure [the combinations of the three reciprocal-lattice vectors [{\bf a}_{i}^*] ([i= 1, 2, 3])] are called the main reflections, the other peaks are satellites. For the latter, at least one of the [h_{4},\ldots, h_{n}] is different from zero.

A second type of modulation is the occupation or composition modulation. Here the structure can again be described on the basis of a basic structure with space-group symmetry. The basic structure positions are occupied with a certain probability by different atom species, or by molecules in different orientations. In CuAu(II), the two lattice positions in a b.c.c. structure are occupied by either Cu and Au or by Au and Cu with a certain probability. This probability function is periodic in one direction with a period that is not a multiple of the lattice constant. In NaNO2, the NO2 molecules are situated at the centre of the orthorhombic unit cell. There are two possible orientations for the V-shaped molecule, and the probability for one of the orientations is a periodic function with periodicity along the a axis. In this case, the modulation wavevector [\alpha{\bf a}^{*}] has a component [\alpha] that strongly depends on temperature in a very narrow temperature range.

If the probability of finding species A in position [{\bf n} + {\bf r}_{j}] or of finding one orientation of a molecule in that point is given by [P_{j}({\bf Q}\cdot{\bf n})], the probability for species B or the other orientation is of course [1-P_{j}({\bf Q}\cdot{\bf n})]. In the diffraction pattern, the spots belong to the Fourier module with basic vectors [{\bf a}^{*}], [{\bf b}^{*}], [{\bf c}^{*}] and [{\bf Q}]. The analogous expression for a more general situation with more modulation wavevectors, or with more species or orientations, is a straightforward generalization.

The first examples of IC phases were found in magnetic systems (see Section 1.5.1.2.3[link] ). For example, holmium has a spiral spin arrangement with a periodicity of the spiral that does not fit with the underlying lattice. For the [\alpha] component ([\alpha = x,y,z]) of the magnetic moment at position [{\bf n}+{\bf r}_{j}] one has in an incommensurate magnetic system a superposition of waves[S_{\alpha}({\bf n}j) = \textstyle\sum\limits_{m}M_{m\alpha j}\cos (2\pi {\bf Q}_{m}\cdot{\bf n} + \varphi_{m\alpha}).\eqno(1.10.1.7)]The most general expression is[S_{\alpha}({\bf n}j) = \textstyle\sum\limits_{{\bf H}\in M^*}M_{\alpha j}({\bf H})\exp (i{\bf H}\cdot{\bf n}),\eqno(1.10.1.8)]where [M^*] is the Fourier module (1.10.1.1)[link].

A following class of quasiperiodic materials is formed by incommensurate composite structures. To this belong misfit structures, intercalates and incommensurate adsorbed layers. An example is Hg3−xAsF6. This consists of a subsystem of AsF6 octahedra forming a (modulated) tetragonal system and two other subsystems consisting of Hg chains, one system of chains in the x direction and one in the y direction. Because the average spacing between the Hg atoms is irrational with respect to the lattice constant of the host AsF6 system in the same direction, the total structure does not have lattice periodicity in the a or b direction.

In general, there are two or more subsystems, labelled by [\nu], and the atomic positions are given by[{\bf n}_{\nu} + {\bf r}_{\nu j} + {\rm modulation},\eqno(1.10.1.9)]where [{\bf n}_{\nu}] belongs to the [\nu]th lattice, and where the modulation is a quasiperiodic displacement from the basic structure. The diffraction pattern has wavevectors[{\bf H} = \textstyle\sum\limits_{\nu}\textstyle\sum\limits_{i_{\nu}=1}^{3}h_{i_{\nu}}{\bf a}_{\nu i_{\nu}}^{*} = \textstyle\sum\limits_{i=1}^{n}h_{i}{\bf a}_{i}. \eqno(1.10.1.10)]Each of the reciprocal-lattice vectors [{\bf a}_{\nu j}^{*}] belongs to the Fourier module [M^*] and can be expressed as a linear combination with integer coefficients of the n basis vectors [{\bf a}_{i}^{*}].

Very often, composite structures consist of a host system in the channels of which another material diffuses with a different, and incommensurate, lattice constant. Examples are layer systems in which foreign atoms intercalate. Another type of structure that belongs to this class is formed by adsorbed monolayers, for example a noble gas on a substrate of graphite. If the natural lattice constant of the adsorbed material is incommensurate with the lattice constant of the substrate, the layer as a whole will be quasiperiodic.

In general, the subsystems can not exist as such. They form idealized lattice periodic structures. Because of the interaction between the subsystems the latter will, generally, become modulated, and even incommensurately modulated because of the mutual incommensurability of the subsystems. The displacive modulation will, generally, contain wavevectors that belong to the Fourier module (1.10.1.10)[link]. However, in principle, additional satellites may occur due to other mechanisms, and this increases the rank of the Fourier module.

The last class to be discussed here is that of quasicrystals. In 1984 it was found (Shechtman et al., 1984[link]) that in the diffraction pattern of a rapidly cooled AlMn alloy the spots were relatively sharp and the point-group symmetry was that of an icosahedron, a group with 120 elements and one that can not occur as point group of a three-dimensional space group. Later, ternary alloys were found with the same symmetry of the diffraction pattern, but with spots as sharp as those in ordinary crystals. These structures were called quasicrystals. Others have been found with eight-, ten- or 12-fold rotation symmetry of the diffraction pattern. Such symmetries are also noncrystallographic symmetries in three dimensions. Sometimes this noncrystallographic symmetry is considered as characteristic of quasicrystals.

Mathematical models for quasicrystals are quasiperiodic two- and three-dimensional tilings, plane or space coverings, without voids or overlaps, by copies of a finite number of `tiles'. Examples are the Penrose tiling or the standard octagonal tiling in two dimensions, and a three-dimensional version of Penrose tiling, a quasiperiodic space filling by means of two types of rhombohedra. For Penrose tiling, all spots of the diffraction pattern are linear combinations of the five basis vectors[{\bf a}_{m}^{*} = \{a\cos[2\pi(m-1)/5], a\sin[2\pi (m-1)/5]\} \quad (m=1,\ldots 5).\eqno(1.10.1.11)]Because the sum of these five vectors is zero, the rank of the spanned Fourier module is four. The Fourier module of the standard octagonal tiling is spanned by[{\bf a}_{m}^{*} = \{a\cos[(m-1)\pi/4], a\sin[(m-1)\pi/4]\} \quad (m=1,\ldots, 4).\eqno(1.10.1.12)]The rank of the Fourier module is four. The rank of the Fourier module of the three-dimensional Penrose tiling, consisting of two types of rhombohedra with a ratio of volumes of [(\sqrt{5}+1)/2], is six and basis vectors point to the faces of a regular dodecahedron.

An atomic model can be obtained by decorating the tiles with atoms, each type of tile in a specific way. Some quasicrystals can really be considered as decorated tilings.

1.10.1.3. Embedding in superspace

| top | pdf |

A simple example of a quasiperiodic function is obtained in the following way. Consider a function of n variables which is periodic with period one in each variable.[f(x_{1},\ldots,x_{n}) = f(x_{1}+1,x_{2},\ldots,x_{n}),\ldots.\eqno(1.10.1.13)]Now take n mutually irrational numbers [\alpha_{i}] and define the function [g(x)] with one variable as[g(x) = f(\alpha_{1}x,\alpha_{2}x,\ldots, \alpha_{n}x).\eqno(1.10.1.14)]Because of the irrationality, the function [g(x)] is not periodic. If we consider the Fourier transform of [f(x_{1},\ldots,x_{n})] we get[f(x_{1},\ldots, x_{n}) = \textstyle\sum\limits_{m_{1}}\ldots \textstyle\sum\limits_{m_{n}} A_{m_{1},\ldots, m_{n}}\exp[2\pi i(m_{1}x_{1} + \ldots + m_{n}x_{n})]\eqno(1.10.1.15)]and consequently[g(x) = \textstyle\sum\limits_{m_{1}\ldots m_{n}}A_{m_{1},\ldots,m_{n}} \exp\left[2\pi i\left(\textstyle\sum\limits_{i = 1}^{n}m_{i}\alpha_{i}\right)x\right],\eqno(1.10.1.16)]which proves that the function is quasiperiodic of rank n with n basis vectors [2\pi\alpha_{i}] in one dimension.

The quasiperiodic function [g(x)] is therefore the restriction to the line [(\alpha_{1}x,\ldots,\alpha_{n}x)] in n-dimensional space. This is a general situation. Each quasiperiodic function can be obtained as the restriction of a periodic function in n dimensions to a subspace that can be identified with the physical space. We denote the n-dimensional space in which one finds the lattice periodic structure (the superspace) by [V_{s}], the physical space by [V_{E}] and the additional space, called internal space, by [V_{I}], such that [V_{s}] is the direct sum of [V_{E}] and [V_{I}]. In the field of quasicrystals, one often uses the name parallel space for [V_{E}] and perpendicular space for [V_{I}].

On the other hand, one can embed the quasiperiodic function in superspace, which means that one constructs a lattice periodic function in n dimensions such that its restriction to physical space is the quasiperiodic function. Take as an example the displacively modulated structure of equation (1.10.1.2)[link]. Compare this three-dimensional structure with the array of lines[\left({\bf n} + {\bf r}_{j} + {\bf f}_{j}({\bf Q}\cdot{\bf n} + t), t \right) \quad({\rm real} \;\;t)\eqno(1.10.1.17)]in four-dimensional space. The restriction to the three-dimensional hyperplane [t = 0] gives exactly the structure (1.10.1.2)[link]. Moreover, the four-dimensional array of lines is lattice periodic. Because [{\bf f}_{j}] is periodic, the array is left invariant if one replaces t by [t+1], and for every lattice vector [{\bf m}] of the basic structure the array is left invariant if one replaces simultaneously t by [t-{\bf Q}\cdot{\bf m}]. This means that the array is left invariant by all four-dimensional lattice vectors of the lattice [\Sigma] with basis[{\bf a}_{si} = ({\bf a}_{i}, -{\bf Q}\cdot{\bf a}_{i})\quad (i = 1, 2, 3), \quad{\bf a}_{4} = (0,1). \eqno(1.10.1.18)]Indeed the quasiperiodic IC phase is the restriction to [V_{E}] ([t = 0]) of the lattice periodic function in four dimensions.

The reciprocal basis for (1.10.1.18)[link] consists of the basis vectors[{\bf a}_{si}^{*} = ({\bf a}_{i}^{*}, 0) \quad (i = 1, 2, 3), \quad {\bf a}_{4}^{*} = ({\bf Q},1). \eqno(1.10.1.19)]These span the reciprocal lattice [\Sigma^{*}]. The projection of this basis on [V_{E}] consists of the four vectors [{\bf a}_{i}^{*}] [(i = 1,2,3)] and [{\bf Q}], and these form the basis for the Fourier module of the quasiperiodic structure.

This is a well known situation. From the theory of Fourier transformation one knows that the projection of the Fourier transform of a function in n dimensions on a d-dimensional subspace is the Fourier transform of the restriction of that n-dimensional function to the same d-dimensional subspace. This gives a way to embed the quasiperiodic structure in a space with as many dimensions as the rank of the Fourier module. One considers the basis of the Fourier module as the projection of n linearly independent vectors in n-dimensional space. This means that for every vector of the Fourier module one has exactly one reciprocal-lattice vector in [V_{s}]. Suppose the quasiperiodic structure is given by some function, for example the density [\rho ({\bf r})]. Then[\rho ({\bf r}) = \textstyle\sum\limits_{{\bf H}\in M^*} \hat{\rho}({\bf H}) \exp (i{\bf H}\cdot{\bf r}).\eqno(1.10.1.20)]One may define a function in n-dimensional space by[\rho_{s}({\bf r}_{s}) = \textstyle\sum\limits_{{\bf H}_{s}\in \Sigma^*} \hat{\rho}({\bf H}) \exp (i{\bf H}_{s}\cdot{\bf r}_{s}),\eqno(1.10.1.21)]where [{\bf H}_{s}] is the unique reciprocal-lattice vector that is projected on the Fourier module vector [{\bf H}]. It is immediately clear that the restriction of [\rho_{s}] to physical space is exactly [\rho]. Moreover, the function [\rho_{s}] is lattice periodic with lattice [\Sigma], for which [\Sigma^{*}] is the reciprocal lattice.

This construction can be performed in the following equivalent way. Consider a point [{\bf r}] in physical space, where one has the quasiperiodic function [\rho ({\bf r})]. The Fourier module of this function is the projection on physical space of the n-dimensional reciprocal lattice [\Sigma^{*}] with basis vectors [{\bf a}_{si}^{*}] ([i = 1, 2,\ldots, n]). The reciprocal lattice [\Sigma^{*}] corresponds to the direct lattice [\Sigma]. A point r in [V_{E}] can also be considered as an element (r, 0) in n-dimensional space. By the translations of [\Sigma], this point is equivalent with a point [{\bf r}_{s}] with lattice coordinates[\xi_{i} = {\rm Frac}\left({\bf a}_{si}^{*}.({\bf r},0)\right) = {\rm Frac}({\bf a}_{i}^{*}\cdot{\bf r})\eqno(1.10.1.22)]in the unit cell of [\Sigma], where Frac(x) is x minus the largest integer smaller than x. If one puts [\rho_{s}({\bf r}_{s}) = \rho ({\bf r})], the function [\rho] determines the function [\rho_{s}] in the unit cell, and consequently in the whole n-dimensional space [V_{s}]. This means that all the information about the structure in [V_{E}] is mapped onto the information inside the n-dimensional unit cell. The information in three dimensions is exactly the same as that in superspace. Only the presentation is different.

In the case in which the crystal consists of point atoms, the corresponding points in d-dimensional physical space [V_{E}] are the intersection of ([n-d])-dimensional hypersurfaces with [V_{E}]. For displacively modulated IC phases in three dimensions with one modulation wavevector, one has [n = 4], [d = 3] and the hypersurfaces are just lines in superspace, as we have seen. For more independent modulation vectors the dimension of the hypersurfaces is larger than one. In this case, as often in the case of composite structures, the ([n-d])-dimensional surfaces do not have borders. This in contrast to quasicrystals, where they are bounded. All these hypersurfaces for which the intersection with physical space gives the atomic positions are called atomic surfaces.

1.10.2. Symmetry

| top | pdf |

1.10.2.1. Symmetry transformations

| top | pdf |

Because the embedded periodic structure in n dimensions has lattice periodicity, it has n-dimensional space-group symmetry as well. It is not a priori clear that such a symmetry group in the unphysical n-dimensional space is relevant for the physical structure, but we shall show here that the superspace description is indeed useful for the description of quasiperiodic systems. First we shall discuss some of the structures of these higher-dimensional space transformations.

Suppose the diffraction pattern has rotational symmetry. Consider for example an orthogonal transformation R that leaves the diffraction pattern invariant. In particular, any basis vector of the module is transformed into an element of the module, i.e. an integral linear combination of the basis vectors.[R{\bf a}_{i}^{*} = \textstyle\sum\limits_{j=1}^{n}M_{ji}{\bf a}_{j}^{*},\quad i = 1,2,\ldots, n. \eqno(1.10.2.1)]Because the matrix M depends on R and acts in reciprocal space, we denote it by [\Gamma^{*}(R) = M]. The matrix [\Gamma^{*}(R)] has integer entries. Because the intensity of the diffraction pattern is not constant on circles around the origin (that would imply that one can not distinguish separate peaks), the orthogonal transformation R is of finite order. Then a theorem from group theory states that [\Gamma^*(R)] is similar to an n-dimensional orthogonal transformation [R_{s}]. The latter certainly has an invariant subspace: the physical space. Therefore, one can find a basis transformation S such that the matrix [\Gamma^{*}(R)] is conjugated to the direct sum of an orthogonal transformation R in [V_{E}] and an additional orthogonal transformation in [V_{I}]:[S\Gamma^{*}(R)S^{-1} = \pmatrix{R&0\cr 0&R_{I} \cr}, \quad R\in O(d), \quad R_{I}\in O(n-d).\eqno(1.10.2.2)]We denote this orthogonal transformation in [V_{s}] as [R_{s}] or as a couple ([R,R_{I}]). Clearly, the transformation [R_{s}] leaves the embedded reciprocal lattice [\Sigma^{*}] invariant. Moreover, this transformation leaves the direct lattice [\Sigma] invariant as well. As always, the action of [R_{s}] on the basis of [\Sigma] for which [{\bf a}_{si}^{*}] form the reciprocal basis is then given by[R_{s}{\bf a}_{si} = \textstyle\sum\limits_{j=1}^{n}\Gamma (R)_{ji}{\bf a}_{sj}\;\;{\rm with}\;\;\Gamma (R)_{ij} = \Gamma^{*}(R^{-1})_{ji}.\eqno(1.10.2.3)]This is the usual relation between the action on a basis and the action on the reciprocal basis.

By construction, the orthogonal transformation [R_{s} = (R,R_{I})] leaves the lattice [\Sigma] invariant, and can therefore belong to the point group of a periodic structure with this lattice. In general, such a point-group element does not leave the periodic structure itself invariant, just as a point group in three dimensions does not leave a crystal with a nonsymmorphic space group invariant. One then has to combine the orthogonal transformation with a translation that in general does not belong to the lattice. Here a translation has components in physical as well as in internal space. A translation can be denoted by ([{\bf a}_{E},{\bf a}_{I}]). Then a general solid motion can be written as[g = \left\{ (R,R_{I}) \left| \right. ({\bf a}_{E},{\bf a}_{I}) \right\}. \eqno(1.10.2.4)]The action of such a transformation on a point [r_{s}] in superspace is given by[gr_{s} = g\left({\bf r}_{E},{\bf r}_{I}\right) = \left(R{\bf r}_{E}+{\bf a}_{E}, R_{I}{\bf r}_{I}+{\bf a}_{I} \right).\eqno(1.10.2.5)]If such a transformation leaves the periodic array of atomic surfaces in superspace invariant, it is a symmetry transformation. In particular, the elements ([{\bf a}_{E},{\bf a}_{I}]) of the translation group [\Sigma] are such symmetry transformations.

1.10.2.2. Point groups

| top | pdf |

The orthogonal transformations that leave the diffraction pattern invariant form a point group K, a finite subgroup of O(d), where d is the dimension of the physical space. All elements act on the basis of the Fourier module as in (1.10.2.1)[link] and the matrices [\Gamma^{*}(K)] form a representation of the group K, an integral representation because the matrices have all integer entries, and reducible because the physical space is an invariant subspace for [\Gamma^{*}(K)]. Because K is finite, this representation is equivalent with a representation in terms of orthogonal matrices. Moreover, by construction [\Gamma^{*}(K)] leaves the n-dimensional reciprocal lattice [\Sigma^{*}] invariant. It is an n-dimensional crystallographic point group. The components R of [R_{s}] form a d-dimensional point group [K_{E}], which is not necessarily crystallographic, and the components [R_{I}] form an ([n-d])-dimensional point group [K_{I}].

Consider as an example an IC phase with orthorhombic basic structure and one independent modulation wavevector [\gamma {\bf c}^{*}] along the c axis. Suppose that the Fourier module, which is of rank four, is invariant under the point group mmm. Then one has for the three generators [\matrix{{m}_x{\bf a}_1^* = -{\bf a}_1^*, \hfill&{m}_x{\bf a}_2^* = {\bf a}_2^*, \hfill&{m}_x{\bf a}_3^* = {\bf a}_3^*, \hfill&{m}_x{\bf a}_4^* = {\bf a}_4^*\hfill\cr {m}_y{\bf a}_1^* = {\bf a}_1^*, \hfill&{m}_y{\bf a}_2^* = -{\bf a}_2^*, \hfill&{m}_y{\bf a}_3^* = {\bf a}_3^*, \hfill&{m}_y{\bf a}_4^* = {\bf a}_4^*\hfill\cr{m}_z{\bf a}_1^* = {\bf a}_1^*, \hfill&{m}_z{\bf a}_2^* = {\bf a}_2^*, \hfill&{m}_z{\bf a}_3^* = -{\bf a}_3^*, \hfill&{m}_z{\bf a}_4^* = -{\bf a}_4^*.\hfill\cr}]Therefore, the corresponding matrices [\Gamma^{*}(R)] are[\displaylines{\pmatrix{ -1&0&0&0\cr 0&1&0&0\cr 0&0&1&0\cr 0&0&0&1}, \quad\pmatrix{ 1&0&0&0\cr 0&-1&0&0\cr 0&0&1&0\cr 0&0&0&1}, \cr\pmatrix{1&0&0&0\cr 0&1&0&0\cr 0&0&-1&0\cr 0&0&0&-1}, \cr\hfill(1.10.2.6)}]which implies that the three generators of the four-dimensional point group are ([m_{x},1]), ([m_{y},1]) and ([m_{z},\bar{1}]).

The diffraction pattern of the standard octagonal tiling has rank four, basis vectors of the Fourier module are[(1,0),\quad(\sqrt{1/2},\sqrt{1/2}),\quad(0,1),\quad(-\sqrt{1/2},\sqrt{1/2})]and the pattern is invariant under a rotation of [\pi /4] and a mirror symmetry. The action of these elements on the given basis of the Fourier module is[\Gamma^{*}(R_{1}) = \pmatrix{0&0&0&-1\cr 1&0&0&0\cr 0&1&0&0\cr 0&0&1&0 \cr},\quad \Gamma^{*}(R_{2}) = \pmatrix{0&0&0&1\cr 0&0&1&0\cr 0&1&0&0\cr 1&0&0&0\cr}.]By a basis transformation, one may bring these transformations into the form[\displaylines{\pmatrix{\cos(\pi{/4})&-\sin(\pi{/4})&0&0\cr \sin(\pi{/4})&\cos(\pi{/4})&0&0\cr 0&0&\cos(3\pi{/4})&-\sin(3\pi{/4})\cr 0&0&\sin(3\pi{/4}) &\cos(3\pi{/4})\cr},\cr \pmatrix{-1&0&0&0\cr 0&1&0&0\cr 0&0&-1&0\cr 0&0&0&1\cr}. }]Therefore, the [\pi /4] rotation in physical space is combined with a [3\pi/4] rotation in internal space in order to get a transformation that leaves a lattice invariant.

A three-dimensional example is the case of a quasicrystal with icosahedral symmetry. For the diffraction pattern all spots may be labelled with six indices with respect to a basis with basis vectors [\eqalign{{\bf a}_1^* &= (0,0,1)\cr{\bf a}_2^* &= (a,0,b)\cr{\bf a}_3^* &= \left(a\cos (2\pi /5), a\sin (2\pi /5),b\right)\cr{\bf a}_4^* &= \left(a\cos (4\pi /5), a\sin (4\pi /5),b\right)\cr{\bf a}_5^* &= \left(a\cos (4\pi /5), -a\sin (4\pi /5),b\right)\cr{\bf a}_6^* &= \left(a\cos (2\pi /5), -a\sin (2\pi /5),b\right),\cr}]with [a = 2/\sqrt{5}] and [b = 1/\sqrt{5}]. The rotation subgroup that leaves the Fourier module invariant is generated by[\displaylines{\Gamma^{*}(A) = \pmatrix{1&0&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 0&0&0&0&0&1\cr 0&1&0&0&0&0 \cr},\cr\Gamma^{*}(B) = \pmatrix{0&0&0&0&0&1\cr 1&0&0&0&0&0\cr 0&0&0&0&1&0\cr 0&0&-1&0&0&0\cr 0&0&0&-1&0&0\cr 0&1&0&0&0&0\cr}. \cr\hfill(1.10.2.7)}]Moreover, there is the central inversion [-E]. The six-dimensional representation of the symmetry group, which is the icosahedral group [\bar{5}\bar{3}m], is reducible into the sum of two nonequivalent three-dimensional irreducible representations. A basis for this representation in the six-dimensional space is then given by[\displaylines{({\bf a}_{1}^{*},c{\bf a}_{1}^{*}) \quad({\bf a}_{2}^{*},-c{\bf a}_{2}^{*}) \quad({\bf a}_{3}^{*},-c{\bf a}_{4}^{*})\cr ({\bf a}_{4}^{*},-c{\bf a}_{6}^{*}) \quad({\bf a}_{5}^{*},-c{\bf a}_{3}^{*})\quad ({\bf a}_{6}^{*},-c{\bf a}_{5}^{*}), \cr\hfill(1.10.2.8)}]which projects on the given basis in [V_{E}].

The point-group elements considered here are pairs of orthogonal transformations in physical and internal space. Orthogonal transformations that do not leave these two spaces invariant have not been considered. The reason for this is that the information about the reciprocal lattice comes from its projection on the Fourier module in physical space. By changing the length scale in internal space one does not change the projection but one would break a symmetry that mixes the two spaces. Nevertheless, quasicrystals are often described starting from an n-dimensional periodic structure with a lattice of higher symmetry. For example, the icosahedral 3D Penrose tiling can be obtained from a structure with a hypercubic six-dimensional lattice. Its reciprocal lattice is that spanned by the vectors (1.10.2.8)[link] where one puts c = 1. The symmetry of the periodic structure, however, is lower than that of the lattice and has a point group in reducible form. Therefore, we shall consider here only reducible point groups, subgroups of the orthogonal group O(n) which have a d-dimensional invariant subspace, identified with the physical space.

The fact that the spaces [V_E] and [V_I] are usually taken as mutually perpendicular does not have any physical relevance. One could as well consider oblique projections of a reciprocal lattice [\Sigma^*] on [V_E]. What is important is that the intersection of the periodic structure with the physical space should be the same in all descriptions. The metric in internal space [V_I] follows naturally from the fact that there is a finite group [K_I].

1.10.2.3. Superspace groups

| top | pdf |

The quasiperiodic function [\rho ({\bf r})] in d dimensions can be embedded as lattice periodic function [\rho_{s}({\bf r}_{s})] in n dimensions. The symmetry group of the latter is the group of all elements g (1.10.2.4)[link] for which[\rho_{s}({\bf r}_{s}) = \rho_{s}(g{\bf r}_{s}) = \rho_{s}(R_{E}{\bf r}+{\bf a}_{E},R_{I}{\bf r}_{I}+{\bf a}_{I}). \eqno(1.10.2.9)]This group is an n-dimensional space group G. It has an invariant subgroup of translations, which is formed by the lattice translations [\Sigma], and the quotient [G/\Sigma] is isomorphic to the n-dimensional point group K. However, not every n-dimensional space group can occur here because we made the restriction to reducible point groups. For example, the n-dimensional hypercubic groups do not occur in this way as symmetry groups of quasiperiodic systems.

The product of two superspace group elements is[\{R_{s1} | {\bf a}_{s1} \} \{R_{s2} | {\bf a}_{s2} \} = \{R_{s1} R_{s2} |{\bf a}_{s1} + R_{s1}{\bf a}_{s2}\}.\eqno(1.10.2.10)]On a lattice basis for [\Sigma], the orthogonal transformations [R_{s1}] and [ R_{s2}] are integer [n\times n] matrices and the translations [{\bf a}_{s1}] and [{\bf a}_{s2}] are column vectors. The orthogonal transformations [R_s] leave the origin invariant. The translations depend on the choice of this origin. For a symmorphic space group there is a choice of origin such that the translations a are lattice translations.

The point-group elements are reducible, which means that in the physical space one has the usual situation. If [d= 3] then the only intrinsic nonprimitive translations are those in screw axes or glide planes. An n-dimensional orthogonal transformation can always be written as the sum of a number r of two-dimensional rotations with rotation angle different from [\pi], a p-dimensional total inversion and a q-dimensional identity transformation. The integers r, p, q may be zero and [2r+p+q=n]. The possible intrinsic nonprimitive translations belong to the q-dimensional space in which the identity acts. For the three examples in the previous section, the internal component of the nonprimitive translation for [m_{x}] and [m_{y}] in the first example can be different from zero, but that for [m_{z}] in the same example is zero. For the octagonal case, only the second generator can have an intrinsic nonprimitive translation in the fourth direction, and for the icosahedral case the two generators have one two-dimensional invariant plane and one pointwise invariant line in [V_{I}].

In the diffraction pattern of an IC phase one can distinguish between main reflections and satellites. A symmetry operation cannot transform a main reflection into a satellite. This implies that for these structures the reciprocal lattice of the basic structure is left invariant by the point group, and consequently the latter must be a three-dimensional crystallographic point group. Therefore, the point groups for IC phases are the same as those for lattice periodic systems. They act in superspace as a representation of a three-dimensional crystallographic point group. This is not true for an arbitrary quasiperiodic structure. The restriction in the general case comes from the requirement that the three-dimensional point group must have a faithful integer matrix representation in superspace. There is a mathematical statement to the effect that the lowest dimension in which a p-fold rotation can be represented as an integer matrix is given by the Euler function, the number of integers smaller than p that do not divide p. For example, for a prime number p this number is p − 1. This implies that if one restricts the rank of the Fourier module (i.e. the dimension of the superspace) to six, only values 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14 and 18 are possible for p. The values 7, 9, 14 and 18 only occur for two-dimensional quasi­periodic structures of rank six. Therefore, the allowable three-dimensional point groups for systems up to rank six are limited to the groups given in Table 1.10.2.1[link]. The possible superspace groups for IC modulated phases of rank four are given in Chapter 9.8[link] of Volume C of International Tables (2004[link]). Superspace groups for quasicrystals of rank [n \leq 6] are given in Janssen (1988[link]).

Table 1.10.2.1| top | pdf |
Allowable three-dimensional point groups for systems up to rank six

Isomorphism classOrderThree-dimensional point groups
[C_1]11
[C_2]22, [\bar{1}], m
[C_3] 3 3
[C_4] 4 4, [\bar{4}]
[C_5] 5 5
[C_6] 6 6, [\bar{6}], [\bar{3}]
[C_8] 8 8, [\bar{8}]
[C_{10}] 10 10, [\overline{10}], [\bar{5}]
[C_{12}] 12 12, [\overline{12}]
[D_2] 4 222, [2/m], 2mm
[D_3] 6 32, 3m
[D_4] 8 422, 4mm, [\bar{4}2m]
[D_5] 10 52, 5m
[D_{6}] 12 622, [\bar{3}m], 6mm, [\bar{6}2m]
[D_8] 16 822, 8mm, [\bar{8}2m]
[D_{10}] 20 [10\,22], [10\,mm], [\overline{10}\,2m], [\bar{5}m]
[D_{12}] 24 [12\,22], [12\,mm], [\overline{12}\,2m]
[C_{4}\times C_2] 8 [4/m]
[C_{6}\times C_2] 12 [6/m]
[C_{8}\times C_2] 16 [8/m]
[C_{10}\times C_2] 20 [10/m]
[C_{12}\times C_2] 24 [12/m]
[D_{2}\times C_2] 8 [mmm]
[D_{4}\times C_2] 16 [4/mmm]
[D_{6}\times C_2] 24 [6/mmm]
[D_{8}\times C_2] 32 [8/mmm]
[D_{10}\times C_2] 40 [10/mmm]
[D_{12}\times C_2] 48 [12/mmm]
T1223
O24432, [\bar{4}3m]
I60532
[T\times C_2] 24 [m\bar{3}]
[O\times C_2] 48 [m\bar{3}m]
[I\times C_2] 120 [\bar{5}\bar{3}m]

The notation of higher-dimensional symmetry groups is discussed in two IUCr reports (Janssen et al., 1999[link], 2002[link]).

1.10.3. Action of the symmetry group

| top | pdf |

1.10.3.1. Action of superspace groups

| top | pdf |

The action of the symmetry group on the periodic density function [\rho_s] in n dimensions is given by (1.10.2.9)[link]. The real physical structure, however, lives in physical space. One can derive from the action of the superspace group on the periodic structure its action on the quasiperiodic d-dimensional one. One knows that the density function in [V_{E}] is just the restriction of that in [V_{s}]. The same holds for the transformed function.[g\rho_{s}({\bf r}_{s}) = \rho_{s}(g^{-1}{\bf r}_{s}) \rightarrow g\rho ({\bf r}) = \rho_{s}[R^{-1}({\bf r}-{\bf a}_{E}),-R_{I}^{-1} {\bf a}_{I}].\eqno(1.10.3.1)]This transformation property differs from that under an n-dimensional Euclidean transformation by the `phase shift' [-R_{I}^{-1}{\bf a}_{I}]. Take for example the IC phase with a sinusoidal modulation. If the positions of the atoms are given by[{\bf n} + {\bf r}_{j} + {\bf A}_{j}\cos (2\pi{\bf Q}\cdot{\bf n}+\varphi_{j}),]then the transformed positions are[R({\bf n} + {\bf r}_{j}) + R{\bf A}_{j}\cos (2\pi{\bf Q}\cdot{\bf n}+\varphi_{j} -R_{I}^{-1}{\bf a}_{I}) + {\bf a}_{E}.\eqno(1.10.3.2)]If the transformation g is a symmetry operation, this means that the original and the transformed positions are the same.[R({\bf n} + {\bf r}_{j}) + {\bf a}_{E} = {\bf n^\prime} + {\bf r}_{j^\prime}]and[R{\bf A}_{j}\cos (2\pi{\bf Q}\cdot{\bf n}+\varphi_{j}-R_{I}^{-1}{\bf a}_{I}) = {\bf A}_{j^\prime}\cos (2\pi{\bf Q}\cdot{\bf n^\prime}+\varphi_{j^\prime}).]This puts, in general, restrictions on the modulation.

Another view of the same transformation property is given by Fourier transforming (1.10.2.9)[link]. The result for the Fourier transform is[g\hat{\rho}_{s}({\bf k}_{s}) = \hat{\rho}_{s}(R_{s}^{-1}{\bf k}_{s}) \exp (-i{\bf k}_{s}\cdot{\bf a}_{s})\eqno(1.10.3.3)]and because there is a one-to-one correspondence between the vectors [{\bf k}_{s}] in the reciprocal lattice and the vectors k in the Fourier module one can rewrite this as[g\hat{\rho}({\bf k}) = \hat{\rho}(R^{-1}{\bf k})\exp (-i{\bf k}\cdot{\bf a}_{E} -{\bf k}_{I}\cdot{\bf a}_{I}).\eqno(1.10.3.4)]For a symmetry element one has [g\hat{\rho}({\bf k})=\hat{\rho}({\bf k})]. Therefore, the superspace group element g is a symmetry transformation of the quasiperiodic function [\rho] if[\hat{\rho}({\bf k}) = \hat{\rho}(R^{-1}{\bf k})\exp (-i{\bf k}\cdot{\bf a}_{E} -{\bf k}_{I}\cdot{\bf a}_{I}).\eqno(1.10.3.5)]This relation is at the basis of the systematic extinctions. If one has an orthogonal transformation R such that this in combination with a translation ([{\bf a}_{E},{\bf a}_{I}]) is a symmetry element and such that [R{\bf k} = {\bf k}], then[\hat{\rho}({\bf k}) = 0\;\;{\rm if}\;\;{\bf k}\cdot{\bf a}_{E}+{\bf k}_{I}\cdot{\bf a}_{I} \neq 2\pi \times\, {\rm integer}.\eqno(1.10.3.6)]Because the structure factor is the Fourier transform of a density function which consists of [\delta] functions on the positions of the atoms, for a quasiperiodic crystal it is the Fourier transform of a quasiperiodic function [\rho ({\bf r})]. Therefore, symmetry-determined absence of Fourier components leads to zero intensity of the corresponding diffraction peaks. Therefore, although there is no lattice periodicity for aperiodic crystals, systematic extinctions follow in the same way from the symmetry as in lattice periodic systems if one considers the n-dimensional space group as the symmetry group.

1.10.3.2. Compensating gauge transformations

| top | pdf |

The transformation property of the Fourier transform of the density given in the previous section can be formulated in another way. Consider a function [\rho ({\bf r})] which is invariant under a d-dimensional Euclidean transformation [\{R|{\bf a}\}] in physical space. Then its Fourier transform satisfies[\hat{\rho}({\bf k}) = \hat{\rho}(R^{-1}{\bf k})\exp (-i{\bf k}\cdot{\bf a}). \eqno(1.10.3.7)]Conversely, if the Fourier transform satisfies this relation, the Euclidean transformation is a symmetry operation for [\rho ({\bf r})]. The two equations (1.10.3.5)[link] and (1.10.3.7)[link] are closely related. One can also write (1.10.3.5)[link] as[\hat{\rho}({\bf k}) = \hat{\rho}(R^{-1}{\bf k})\exp (-i{\bf k}\cdot{\bf a}) \exp [i\Phi (R,{\bf k})], \eqno(1.10.3.8)]where [\Phi (R,{\bf k})] can be considered as a gauge transformation that compensates for the phase shift: it is a compensating gauge transformation. It is a function that is linear in k, [\Phi (R,{\bf k}+{\bf k}^\prime) = \Phi (R,{\bf k}) + \Phi (R,{\bf k}^\prime)\; ({\rm mod}\;2\pi), \eqno(1.10.3.9)]and satisfies a relation closely related to the one satisfied by nonprimitive translations.[\Phi (R,{\bf k}) + \Phi (S,R{\bf k}) = \Phi (RS, {\bf k}) \; ({\rm mod}\;2\pi).\eqno(1.10.3.10)][Recall that a system of nonprimitive translations [{\bf u}(R)] satisfies [{\bf u}(R)+R{\bf u}(S)={\bf u}(RS)] modulo lattice translations.] Therefore, the Euclidean transformation [\{R|{\bf a}\}] combined with the compensating gauge transformation with gauge function [\Phi (R,{\bf k})] is a symmetry transformation for [\rho ({\bf r})] if equation (1.10.3.8)[link] is satisfied. This is a three-dimensional formulation of the superspace group symmetry relation (1.10.3.5)[link].

1.10.3.3. Irreducible representations of three-dimensional space groups

| top | pdf |

A third way to describe the symmetry of a quasiperiodic function is by means of irreducible representations of a space group. For the theory of these representations we refer to Chapter 1.2[link] on representations of crystallographic groups.

Consider first a modulated IC phase. Suppose the positions of the atoms are given by[{\bf n} + {\bf r}_{j} + {\bf u}_{{\bf n}j},\eqno(1.10.3.11)]where n belongs to the lattice, [{\bf r}_{j}] is a position inside the unit cell and [{\bf u}_{{\bf n}j}] is a displacement. If the structure is quasiperiodic with Fourier module [M^{*}], the vectors [{\bf u}_{{\bf n}j}] can be written as a superposition of normal modes.[{\bf u}_{{\bf n}j} = \textstyle\sum\limits_{{\bf k}\in M^*,\nu} Q_{k\nu}{\boldvarepsilon}({\bf k}\nu |j) {\rm e}^{i{\bf k}\cdot{\bf n}} + c.c.,\eqno(1.10.3.12)]where the coefficient [Q_{{\bf k}\nu}] is a normal coordinate, [\nu] denotes the band index and [{\boldvarepsilon}({\bf k}\nu |j)] denotes the polarization of the normal mode. The normal coordinates transform under a space group according to one of its irreducible representations. The relevant space group here is that of the basic structure. For the simple case of a one-dimensional irreducible representation, for each [{\bf k}] the effect is simply multiplication by a factor of absolute value unity. For example, for the modulated phase with basic space group Pcmn and wavevector [{\bf k} = \gamma {\bf c}^{*}] there are four non-equivalent one-dimensional representations. It depends on the band index which representation occurs in the decomposition. The space-group element [\{R|{\bf a}\}], for which [R{\bf q}={\bf q}] (modulo reciprocal lattice) acts on [Q_{{\bf k}\nu}] according to[Q_{{\bf k}\nu} \rightarrow Q_{{\bf k}\nu}\exp (i{\bf k}.\cdot{\bf a})\chi_{\nu}(R),]where [\chi_{\nu}(R)] is the character of R in an irreducible representation associated with the branch [\nu]. Because the character of a one-dimensional representation is of absolute value unity, one may write it as [\exp[i\varphi_{\nu}(R,{\bf k})]]. Consequently, if the decomposition of the displacement contains only the vectors [\pm {\bf k}], the factor [\exp[i\varphi_{\nu}(R,{\bf k})]] describes a shift in the modulation function.

Consider again as an example a basic structure with space group Pcmn and a modulation wavevector [\gamma {\bf c}^{*}]. The point group [K_{{\bf k}}] that leaves the modulation wavevector invariant is generated by [m_{y}] and [m_{x}]. This point group mm2 has four elements and four irreducible representations, all one-dimensional. One of them has for the character [\chi(m_{x}) = +1], [\chi(m_{y}) = -1]. If the displacements of the atoms are described by a normal mode belonging to this irreducible representation, then the compensating phase shifts for [c_{x}] and [m_{y}] are, respectively, 0 and [\pi]. In the notation for superspace groups, this is the group Pcmn(00[\gamma])1s[\bar{1}]. The same structure can be described by the irreducible representation characterized as [\Delta_{3}], because the modulation wavevector is the point [\Delta] in the Brillouin zone and the irreducible representation [\Gamma_{3}] has the character mentioned above.

In this way there is a correspondence between superspace groups for (3 + 1)-dimensional modulated structures and two-dimensional irreducible representations of three-dimensional space groups.

1.10.4. Tensors

| top | pdf |

1.10.4.1. Tensors in higher-dimensional spaces

| top | pdf |

A vector in an n-dimensional space V transforms under an element of a point group as [{\bf r}\rightarrow R{\bf r}]. With respect to a basis [{{\bf a}_i}], the coordinates and basis vectors transform according to [\eqalign{{\bf a}_{i}^\prime &= \textstyle\sum\limits_{j = 1}^n R_{ji}{\bf a}_j\cr {\bf r} &= \textstyle\sum\limits_{i = 1}^{n}x_{i}{\bf a}_{i} \rightarrow {\bf r}^\prime = \textstyle\sum\limits_{i=1}^{n}x_{i}^\prime{\bf a}_{i}, \quad x_{i}^\prime = \textstyle\sum\limits_{j = 1}^{n}R_{ij}x_{j}}]and the reciprocal basis vectors and coordinates in reciprocal space according to [\eqalign{{\bf a}_i^{*}&=\textstyle\sum\limits_{j=1}^{n}R_{ij}{\bf a}_{j}^{*^\prime}\cr {\bf k} &= \textstyle\sum\limits_{i = 1}^{n}\kappa_{i}{\bf a}_{i}^{*} \rightarrow {\bf k}^\prime = \textstyle\sum\limits_{i = 1}^{n}\kappa_{i}^\prime{\bf a}_{i}^{*},\quad \kappa_{i}^\prime = \textstyle\sum\limits_{j = 1}^{n}R^{-1}_{ji}\kappa_{j}.}]

With respect to an orthonormal basis in V the transformations are represented by orthogonal matrices. For orthogonal matrices [R^{-1}=R^{{T}}], the vectors in reciprocal space transform in exactly the same way as in direct space: [\eqalign{{\bf r} &= \textstyle\sum\limits_{i = 1}^n x_i{\bf e}_i \rightarrow {\bf r}^\prime = \textstyle\sum\limits_{i = 1}^n x_i^\prime{\bf e}_i \quad x_i^\prime = \textstyle\sum\limits_{i = 1}^n R_{ij}x_i\cr {\bf k} &= \textstyle\sum\limits_{i = 1}^n \kappa_i {\bf e}_i^* \rightarrow {\bf k}^\prime = \textstyle\sum\limits_{i = 1}^n \kappa_i^\prime{\bf e}_i^* \quad \kappa_i^\prime = \textstyle\sum\limits_{i = 1}^n R_{ij}\kappa_j.}]As discussed in Section 1.2.4[link] , a tensor is a multilinear function of vectors and reciprocal vectors. Consider for example a tensor of rank two, the metric tensor g. It is a function of two vectors [{\bf r}_{1}] and [{\bf r}_{2}] which results in the scalar product of the two. [g({\bf r}_{1},{\bf r}_{2}) = {\bf r}_{1}\cdot{\bf r}_{2}.]It clearly is a symmetric function because [g({\bf r}_{1},{\bf r}_{2}) = g({\bf r}_{2},{\bf r}_{1})]. It is a function that is linear in each of its arguments and therefore [g({\bf r}_{1},{\bf r}_{2}) = g\left(\textstyle\sum\limits_{i=1}^{n}x_{i}{\bf e}_{i}, \textstyle\sum\limits_{j=1}^{n}y_{j}{\bf e}_{j}\right) = \textstyle\sum\limits_{ij}x_{i}y_{j}\delta_{ij} = \textstyle\sum\limits_{i}x_{i}y_{j}]if [x_{i}] and [y_{j}] are Cartesian coordinates of [{\bf r}_{1}] and [{\bf r}_{2}], respectively. For another basis, for example a lattice basis, one has coordinates [\xi_{i}] and [\eta_{j}], and the same function becomes [g({\bf r}_{1},{\bf r}_{2}) = g\left(\textstyle\sum\limits_{i=1}^{n}\xi_{i}{\bf a}_{i}, \textstyle\sum\limits_{j=1}^{n}\eta_{j}{\bf a}_{j}\right) = \textstyle\sum\limits_{ij}\xi_{i}\eta_{j}g_{ij}\eqno(1.10.4.1)]with [g_{ij}=g({\bf a}_{i},{\bf a}_{j})]. The relation between the Cartesian tensor components and the lattice tensor components follows from the basis transformation from orthonormal to a lattice basis. If [{\bf a}_{j} = \textstyle\sum\limits_{k}S_{kj}{\bf e}_{k},\eqno(1.10.4.2)] then the lattice tensor components are [g_{ij} = \textstyle\sum\limits_{k}S_{ki}S_{kj}.]For example, in the two-dimensional plane a lattice spanned by [a(1,0)] and [a(-{\textstyle{1\over2}},{\textstyle{1\over2}}\sqrt{3})] has a basis obtained from an orthonormal basis by the basis transformation [S = a\pmatrix{1&-{1\over2} \cr 0&{1\over2}\sqrt{3} \cr}]and consequently the tensor components in lattice coordinates are [g_{ij} = \pmatrix{a^{2}&-{1\over2} a^{2}\cr -{1\over2} a^{2}&a^{2} \cr}. ]

The transformation of the tensor g under an orthogonal transformation follows from its definition. The transformation of the Cartesian tensor under the orthogonal transformation R is [g_{ij}^\prime = \textstyle\sum\limits_{kl}R_{ki}R_{lj}g_{kl} = \textstyle\sum\limits_{kl}R_{ki}R_{lj}\delta_{kl} = \delta_{ij}]because of the fact that the matrix [R_{ij}] is orthogonal. The transformation of the tensor components with respect to the lattice basis, on which R is given by [\Gamma (R)], is [g_{ij}^\prime = \textstyle\sum\limits_{kl}\Gamma (R)_{ki}\Gamma (R)_{lj}g_{kl},\eqno(1.10.4.3)] or in matrix form [g^\prime=\Gamma (R)^{T}g\Gamma (R)].

The metric tensor is invariant under a point group K if [g_{ij} = \textstyle\sum\limits_{kl}\Gamma (R)_{ki}\Gamma (R)_{lj}g_{kl}\quad \forall\;R\in K.\eqno(1.10.4.4)]On the one hand this formula can be used to determine the symmetry of a lattice with metric tensor g and on the other hand one may use it to determine the general form of a metric tensor invariant under a given point group. This comes down to the determination of the free parameters in g for given group of matrices [\Gamma (K)]. These are the coordinates in the space of invariant tensors.

1.10.4.2. Tensors in superspace

| top | pdf |

The tensors occurring for quasiperiodic structures are defined in a higher-dimensional space, but this space contains as privileged subspace the physical space. Since physical properties are measured in this physical space, the coordinates are not all on the same footing. This implies that sometimes one has to make a distinction between the various tensor elements as well.

The distinction between physical and internal (or perpendicular) coordinates can be made explicit by using a split basis. This is a basis for the superspace such that the first d basis vectors span the physical subspace and the other nd basis vectors the internal space. A lattice basis is, generally, not a split basis.

Let us consider again the metric tensor which is used to characterize higher-dimensional lattices as well, and in particular those corresponding to quasiperiodic structures. The elements [g_{ij}=g({\bf a}_i,{\bf a}_j)] transform according to [g^\prime_{ij} = g({\bf a}^\prime_i,{\bf a}^\prime_j) = \textstyle\sum\limits_{kl}R_{ki}R_{lj}g_{kl}.]The symmetry of an n-dimensional lattice with metric tensor g is the group of nonsingular [n\times n] integer matrices S satisfying [g = S^{T} g S,\eqno(1.10.4.5)]where T means the transpose. For a lattice corresponding to a quasiperiodic structure, this group is reducible into a d- and an (nd)-dimensional component, where d is the dimension of physical space. This means that the d-dimensional component, which forms a finite group, is equivalent with a d-dimensional group of orthogonal transformations. In general, however, this does not leave a lattice in physical space invariant. However, it leaves the Fourier module of the quasiperiodic structure invariant. The basis vectors, for which the metric tensor determines the mutual relation, belong to the higher-dimensional superspace. Therefore, in this case the external and internal components of the basis vectors do not need to be treated differently. For the metric tensor g on a split basis one has [g_{ij}=0 \quad{\rm if}\;\;i\leq d, j> d \;\;{\rm or}\;\;i>d,j\leq d.]

A quasiperiodic structure has an n-dimensional lattice embedding such that the intersection of [\Sigma] with the physical space [V_{E}] does not contain a d-dimensional lattice. Because of the incommensurability, however, there are lattice points of [\Sigma] arbitrarily close to [V_{E}]. This means that by an arbitrarily small shear deformation one may get a lattice in the physical space. The deformed quasiperiodic structure then becomes periodic. In general, the symmetry of the lattice then changes. This is certainly the case if the point group of the quasiperiodic structure is noncrystallographic, because then there cannot be a lattice in physical space left invariant by such a point group. For a given lattice [\Sigma] with symmetry group K one may ask which subgroups allow a deformation of the lattice that gives periodicity in [V_{E}].

Physical tensors give often relations between vectorial or tensorial properties. Then they are multilinear functions of p vectors (and possibly q reciprocal vectors). An example is the dielectric tensor [\varepsilon] that gives the relation between E and D fields. This relation and the corresponding expression for the free energy F are [D_{i} = \textstyle\sum\limits_{j}\varepsilon_{ij}E_{j}\;\;{\rm or}\;\;F = \textstyle\sum\limits_{ij} E_{i}\varepsilon_{ij}E_{j} = \varepsilon ({\bf E},{\bf E}).\eqno(1.10.4.6)]Therefore, the [\varepsilon] tensor is a bilinear function of vectors. The difference from the metric tensor is that here the vectors E and D are physical quantities which have d components and lie in physical space. The transformation properties therefore only depend on the physical-space components [R_{E}] of the superspace point group, and not on the full transformations R.

An intermediate case occurs for the strain. The strain tensor S gives the relation between a displacement and its origin: the point [{\bf r}] is displaced to [{\bf r}+\Delta{\bf r}] with [\Delta{\bf r}] linear in [{\bf r}]: [\Delta{\bf r}_{i} = \textstyle\sum\limits_{j}S_{ij}{\bf r}_{j}.]In ordinary elasticity, both [{\bf r}] and [\Delta{\bf r}] belong to the physical space, and the relevant tensor is the symmetric part of S: [{\textstyle{1\over2}}(\partial_i \Delta{\bf r}_j + \partial_j \Delta{\bf r}_i).]For a quasiperiodic structure, [\Delta{\bf r}] may be either a vector in physical space or in superspace and may depend both on physical and internal coordinates. That means that the matrix [\sigma] is either [d\times d], or [n\times d] or [n\times n]. Displacements in physical space are said to affect the phonon degrees of freedom, those in internal space the phason degrees of freedom. The phonon and phason displacements are functions of the physical-space coordinates. The transformation of the strain tensor under an element of a superspace group is [\eqalign{S_{ij}^\prime &= \textstyle\sum\limits_{k=1}^{d}\textstyle\sum\limits_{l=1}^{d}R_{Eki}R_{Elj}S_{kl} \hbox{ for phonon degrees},\cr S_{ij}^\prime &= \textstyle\sum\limits_{k=d+1}^{n}\textstyle\sum\limits_{l=1}^{d}R_{Iki}R_{Elj}S_{kl} \hbox{ for phason degrees},\cr S_{ij}^\prime &= \textstyle\sum\limits_{k=1}^{n}\textstyle\sum\limits_{l=1}^{n}R_{ski}R_{slj}S_{kl} \hbox{ for the general case.} \cr }]The first two of these expressions apply only to a split basis, but the third can be written on a lattice basis. [\textstyle\sum\limits_{k,l=1}^{n}\Gamma (R)_{ki}\Gamma (R)_{lj}S_{kl}.\eqno(1.10.4.7)]

The tensor of elastic stiffnesses c gives the relation between stress T and strain S. The stress tensor is a physical tensor of rank two and dimension three. For the phonon strain one has [S_{ij} = \textstyle\sum\limits_{kl}^{3}c_{ijkl}T_{kl},\quad (i,j=1,\ldots, 3).\eqno(1.10.4.8)]The phonon part of the elasticity tensor is symmetric under interchange of ij and kl, i and j, and k and l. It can be written in the usual notation [c_{\mu\nu}] with [\mu, \nu = 1,2,\ldots, 6] with 1 = (11), 2 = (22), 3 = (33), 4 = (23), 5 = (13), 6 = (12). Its transformation property under a three-dimensional orthogonal transformation is [c_{ijkl}^{\prime} = \textstyle\sum\limits_{i{^\prime}j{^\prime}k{^\prime}l{^\prime}}R_{i{^\prime}i}R_{j{^\prime}j}R_{k{^\prime}k}R_{l^\prime l} c_{i{^\prime}j{^\prime}k{^\prime}l^\prime}.]For the phason part a similar elasticity tensor is defined. This and the third elastic contribution, the coupling between phonons and phasons, will be discussed in Section 1.10.4.5.[link]

1.10.4.3. Inhomogeneous tensors

| top | pdf |

A vector field in d-dimensional space assigns a vector to each point of the space. This vector-valued function may, for a quasiperiodic system, have values in physical space or in superspace. In both cases one has the transformation property[g{\bf f}_{i}({\bf r}) = \textstyle\sum\limits_{j}R_{ji}{\bf f}_{j}(g^{-1}{\bf r}). \eqno(1.10.4.9)]For a vector field in physical space, i and j run over the values 1, 2, 3. This vector field may, however, be quasiperiodic. This means that it may be embedded in superspace. Then[g{\bf f}_{i}({\bf r}_{s}) = \textstyle\sum\limits_{j=1}^{3}R_{Eji}{\bf f}_{j} \left [R_{E}^{-1}({\bf r}-{\bf a}),R^{-1}_{I}({\bf r}_{I}-{\bf a}_{I})\right]. \eqno(1.10.4.10)]Here i = 1, 2, 3. If the vector field has values in superspace, as one can have for a displacement, one has[g{\bf f}_{i}({\bf r}_{s}) = \textstyle\sum\limits_{j=1}^{n}R_{sji}{\bf f}_{j} \left[R_{E}^{-1}({\bf r}-{\bf a}),R^{-1}_I({\bf r}_{I}-{\bf a}_{I}) \right]. \eqno(1.10.4.11)]Here [i = 1,\ldots, n]. For Cartesian coordinates with respect to a split basis, [R_{s}] acts separately on physical and internal space and one has[g{\bf f}_{i}({\bf r}_{s}) = \textstyle\sum\limits_{j = d+1}^{n}R_{Iji}{\bf f}_{j} \left[R_{E}^{-1}({\bf r}-{\bf a}),R^{-1}_I({\bf r}_{I}-{\bf a}_{I})\right]\eqno(1.10.4.12)]for [i = d+1,\ldots, n].

Just as for homogeneous tensors, inhomogeneous tensors may be divided into physical tensors with components in physical space only and others that have components with respect to an n-dimensional lattice. A physical tensor of rank two transforms under a space-group element [g=\{R_{s}|{\bf a}_{s}\}] as[(gT)_{ij}({\bf r}_{s}) = \textstyle\sum\limits_{k=1}^{d}\textstyle\sum\limits_{l=1}^{d}R_{Eki}R_{Elj}T_{kl} [R_{E}^{-1}({\bf r}_{E}-{\bf a}_{E}),R_{I}^{-1}({\bf r}_{I}-{\bf a}_{I})].\eqno(1.10.4.13)]This implies the following transformation property for the Fourier components:[(g\hat{T})_{ij}({\bf k}) = \textstyle\sum\limits_{k=1}^{d}\textstyle\sum\limits_{l=1}^{d}R_{Eki}R_{Elj}T_{kl} (R_{E}^{-1}{\bf k})\exp (iR_{E}{\bf k}.{\bf a}_{E}+iR_{I}{\bf k}_{I}. {\bf a}_{I}).\eqno(1.10.4.14)]This gives relations between various Fourier components and restrictions for wavevectors [{\bf k}] for which [R_{E}{\bf k} = {\bf k}]:[\hat{T}_{ij}({\bf k}) = \textstyle\sum\limits_{k=1}^{d}\textstyle\sum\limits_{l=1}^{d}R_{Eki}R_{Elj} T_{kl}({\bf k})\exp (iR_{E}{\bf k}.{\bf a}_{E}+iR_{I}{\bf k}_{I}. {\bf a}_{I}). \eqno(1.10.4.15)]

For tensors with superspace components, the summation over the indices runs from 1 to n. An invariant tensor then satisfies[\hat{T}_{ij}({\bf k}) = \textstyle\sum\limits_{k=1}^{n}\textstyle\sum\limits_{l=1}^{n}R_{ski}R_{slj} T_{kl}({\bf k})\exp (iR_{E}{\bf k}.{\bf a}_{E}+iR_{I}{\bf k}_{I}. {\bf a}_{I}).\eqno(1.10.4.16)]The generalization to higher-rank tensors is straightforward.

1.10.4.4. Irreducible representations

| top | pdf |

For the characterization of vectors and tensors one needs the irreducible and vector representations of the point groups. If the point group is crystallographic in three dimensions, these can be found in Chapter 1.2[link] . All point groups for IC phases or composite structures belong to this category. Exceptions are the point groups for quasicrystals. For the finite point groups for structures up to rank six these are given in Table 1.10.5.1[link]. This table presents:

  • (1) The character tables for the point groups[\eqalign{&5, \bar{5}, 5m, 52, \bar{5}m\cr &10,\overline{10},10/m,10\,mm,10\, 22,\overline{10}\,2m,10/mmm\cr & 8,\bar{8},8/m,8mm,822,\bar{8}2m,8/mmm\cr &12,\overline{12},12/m,12\,mm,12 \,22,\overline{12}\,2m,12/mmm\cr & 532,\bar{5}\bar{3}m.}]

  • (2) Matrices for the generators in the irreducible representations of the groups[\bar{5}m, 10/mmm, 8/mmm, 12/mmm, \bar{5}\bar{3}m.]

  • (3) The vector representations and some tensor representations for the groups in the systems[\bar{5}m, 10/mmm, 8/mmm, 12/mmm, \bar{5}\bar{3}m.]

The character tables can be used to determine the number of independent tensor elements. This is the dimension of subspace of tensors transforming with the identity representation. Tensors transform according to (properly symmetrized or antisymmetrized) tensor products of vector representations. The number of times the identity representation occurs in the decomposition of the tensor product into irreducible components is equal to the number of independent tensor elements and can be calculated with the multiplicity formula. A number of examples are given in the following section.

1.10.4.5. Determining the number of independent tensor elements

| top | pdf |

1.10.4.5.1. Piezoelectric tensor

| top | pdf |

(See Sections 1.1.4.4.3[link] and 1.1.4.10.1[link] .) The strain in a crystal is determined by its displacement field. For a quasiperiodic crystal, this displacement can have components in the physical space [V_E] as well as in the internal space [V_I]. The first implies a local displacement of the material, the latter corresponds to a local deformation because of the shift in the internal coordinate, which is, for example, the phase of a modulation wave or a phason jump for a quasicrystal. The displacement in the point [{\bf r}] is [u = [u_{E}({\bf r}),u_{I}({\bf r})]]. Denote [u_E] by [{\bf v}] and [u_I] by [{\bf w}]. The strain tensor then is given by [\partial_i v_j] and by [\partial_i w_k]. Here i and j run from 1 to the physical dimension d, and k from 1 to the internal dimension nd. The antisymmetric part of [\partial_i v_j] corresponds to a global rotation, which does not lead to an energy change. Therefore, the relevant tensors are [\displaylines{e_{ij} = (\partial_i v_j + \partial_j v_i)/2,\quad f_{ik} = \partial_i w_k, \cr (i,j=1,\ldots, d; k=1,\ldots, n-d).\cr\hfill(1.10.4.17)}]Both the phonon part e and the phason part f may be coupled to an external electric field E. A linear coupling is given by the piezoelectric tensor [p_{ijk}]. The free energy is given by [F = \textstyle\int {\rm d}r \;\left(\textstyle\sum\limits_{ijk}p_{ijk}^{e}e_{ij}E_{k}+\textstyle\sum\limits_{ijk}p_{ijk}^f f_{ij}E_{k} \right).]The tensor e transforms with the symmetrized square of the vector representation in physical space, the tensor f according to the product of the vector representations in physical and internal space. Then [p^e] and [p^f] transform according to the product of these two representations with the vector representation in physical space, because E is a physical vector.

As an example, consider the decagonal phase with point group 10 mm(103mm). The physical space is three-dimensional and carries a (2 + 1)-reducible representation ([\Gamma_1 \oplus \Gamma_5]), the internal space an irreducible two-dimensional representation ([\Gamma_7]). The symmetrized square of the first is six-dimensional, and the product of first and second is also six-dimensional. The products of these two with the three-dimensional vector representation in physical space are both 18-dimensional. The first contains the identity representation three times, the other does not contain the identity representation. This implies that the piezoelectric tensor has three independent tensor elements, all belonging to [p^e]. The tensor [p^f] is zero.

1.10.4.5.2. Elasticity tensor

| top | pdf |

(See Section 1.3.3.2[link] .) As an example of a fourth-rank tensor, we consider the elasticity tensor. The lowest-order elastic energy is a bilinear expression in e and f: [F = \textstyle\int {\rm d}{\bf r}\; \left({\textstyle{1 \over2}} \textstyle\sum\limits_{ijkl} c_{ijkl}^{E}e_{ij}e_{kl}+ {\textstyle{1 \over2}}\textstyle\sum\limits_{ijkl}c_{ijkl}^{I}f_{ij}f_{kl}+\textstyle\sum\limits_{ijkl}c_{ijkl}^{EI} e_{ij}f_{kl} \right). \eqno(1.10.4.18)]

The elastic free energy is a scalar function. The integrand must be invariant under the operations of the symmetry group. When [\Gamma_E (K)] is the vector representation of K in the physical space (i.e. the vectors in [V_E] transform according to this representation) and [\Gamma_I (K)] the vector representation in [V_I], the tensor [e_{ij}] transforms according to the symmetrized square of [\Gamma_E] and the tensor [f_{ij}] transforms according to the product [\Gamma_{E}\otimes\Gamma_{I}]. Let us call these representations [\Gamma_e] and [\Gamma_f], respectively. This implies that the term that is bilinear in e transforms according to the symmetrized square of [\Gamma_e], that the term bilinear in f transforms according to the symmetrized square of [\Gamma_f], and that the mixed term transforms according to [\Gamma_e\otimes\Gamma_f]. The number of elastic constants follows from their transformation properties. If d = 3 and n = 3 + p, the number of constants [c^E] is 21, the number of constants [c^I] is [3p(3p+1)/2] and the number of [c^{IE}] is 18p. Therefore, without symmetry conditions, there are altogether [3(2+p)(7+3p)/2] elastic constants. For arbitrary dimension d of the physical space and dimension n of the superspace this number is [\eqalign{&d(d+1)(d^2 +d+2)/8+pd(pd+1)/2+d^2 (d+1)p/2\cr&\quad = d(2p+d+1)(2+d+d^2+2pd)/8.}]The number of independent elastic constants is the number of independent coefficients in F, and this is given by the number of invariants, i.e. the number of times the identity representation occurs as irreducible component of, respectively, the symmetrized square of [\Gamma_e], the symmetrized square of [\Gamma_f], and of [\Gamma_e\otimes\Gamma_f]. The first number is the number of elastic constants in classical theory. The other elastic constants involve the phason degrees of freedom, which exist for quasiperiodic structures. The theory of the generalized elasticity theory for quasiperiodic crystals has been given by Bak (1985[link]), Lubensky et al. (1985[link]), Socolar et al. (1986[link]) and Ding et al. (1993[link]).

As an example, we consider an icosahedral quasicrystal. The symmetry group 532 has five classes, which are given in Table 1.10.5.1[link]. The vector representation is [\Gamma_{2}]. It has character [\chi (R) =] [3,1+\tau, -\tau, 0,-1]. The character of its symmetrized square is 6, 1, 1, 0, 2. Then the character of the representation with which the elasticity tensor transforms is 21, 1, 1, 0, 5. This representation contains the trivial representation twice. Therefore, there are two free parameters ([c_{1111}] and [c_{1122}]) in the elasticity tensor for the phonon degrees of freedom.

For the phason degrees of freedom, the displacements transform with the representation [\Gamma_{3}]. In this case, the phason elasticity tensor transforms with the symmetrized square of the product of [\Gamma_{2}] and [\Gamma_{3}]. Its character is 45, 0, 0, 0, 5. This representation contains the identity representation twice. This implies that this tensor also has two free parameters.

Finally, the coupling term transforms with the product of the symmetrized square of [\Gamma_2], [\Gamma_2] and [\Gamma_3]. This representation has character 54, −1, −1, 0, 2 and consequently contains the identity representation once. In total, the number of independent elastic constants is five for icosahedral tensors. The fact that we have only used the rotation subgroup 532, instead of the full group [\bar{5}\bar{3}m], does not change this number. The additional central inversion makes the irreducible representations either even or odd. The elasticity tensors should be even, and there are exactly as many even irreducible representations as odd ones. This is shown in Table 1.10.4.1[link] (cf. Table. 1.10.5.1[link] for the character table of the group 532).

Table 1.10.4.1| top | pdf |
Characters of the point group [10\,mm(10^3\,mm)] for representations relevant for elasticity

[\tau = (\sqrt{5}-1)/2].

RepresentationClassesReduction
EAA2BAB
1 12 122015
[\Gamma_E] 3 [1+\tau] [-\tau] 0 [-1] [\Gamma_2]
[\Gamma_I] 3 [-\tau] [1+\tau]0 [-1] [\Gamma_3]
[\Gamma_{E}^{2}]9[2+\tau] [1-\tau] 0 1  
[\Gamma_{E}(g^2)]3 [-\tau] [1+\tau] 0 3  
[\Gamma_{e}=(\Gamma_E)_{s}^2] 6 1 1 0 2 [\Gamma_{1}+\Gamma_{5}]
[\Gamma_{e}^{2}] 36 1 1 0 4  
[\Gamma_{e}(g^{2})] 6 1 1 0 6  
[(\Gamma_{e})_{s}^{2}] 21 1 1 0 5 [2\Gamma_{1}+\Gamma_{4}+3\Gamma_{5}]
[\Gamma_{f}=\Gamma_E\times\Gamma_I] 9 [-1] [-1] 0 1 [\Gamma_4 +\Gamma_5]
[(\Gamma_{f})_{s}^{2}] 45 0 0 0 5 [2\Gamma_{1}] [+] [\Gamma_{2}] [+] [\Gamma_{3}] [+] [3\Gamma_{4}] [+] [5\Gamma_{5}]
[\Gamma_{e}\times\Gamma_{f}] 54 [-1] [-1] 0 2 [\Gamma_{1}] [+] [2\Gamma_{2}] [+] [2\Gamma_{3}] [+] [4\Gamma_{4}] [+] [5\Gamma_{5}]

1.10.4.5.3. Electric field gradient tensor

| top | pdf |

As an example, we consider a symmetric rank-two tensor, e.g. an electric field gradient tensor, in a system with superspace group symmetry Pcmn(00γ)1s[\bar{1}]. The Fourier transform of the tensor [T_{ij}] is nonzero only for multiples of the vector [\gamma {\bf c}^{*}]. The symmetry element consisting of a mirror operation [M_y] and a shift [{\textstyle{1 \over 2}}{\bf a}_{s4}] in [V_I] then has [R_{E}{\bf k}={\bf k},\quad{\bf k}\cdot{\bf a}_{E}=0,\quad R_{I}=+1,\quad{\bf k}_{I}\cdot{\bf a}_{I}=\pi. ] Then equation (1.10.4.15)[link] leads to the relation [\displaylines{\hat{T}(m\gamma{\bf c}^{*}) = \pmatrix{a_{11}&a_{12}&a_{13}\cr a_{21}&a_{22}&a_{23}\cr a_{31}&a_{32}&a_{33}\cr}\cr = (-1)^{m}\pmatrix{1&0&0\cr 0&-1&0\cr 0&0&1\cr} \pmatrix{a_{11}&a_{12}&a_{13}\cr a_{21}&a_{22}&a_{23}\cr a_{31}&a_{32}&a_{33}\cr} \pmatrix{1&0&0\cr 0&-1&0\cr 0&0&1\cr}}]with solution [\displaylines{\hat{T} = \pmatrix{a_{11}&0&a_{13}\cr 0&a_{22}&0\cr a_{13}&0&a_{33}\cr} \,\,\,(m\,\,{\rm even}),\cr \hat{T} = \pmatrix{0&a_{12}&0\cr a_{12}&0&a_{23}\cr 0&a_{23}&0\cr} \,\,\,(m\,\,{\rm odd}).}]This symmetry of the tensor can, for example, be checked by NMR (van Beest et al., 1983[link]).

1.10.4.6. Determining the independent tensor elements

| top | pdf |

In the previous sections some physical tensors have been studied, for which in a number of cases the number of the independent tensor elements has been determined. In this section the problem of determining the invariant tensor elements themselves will be addressed.

Consider an orthogonal transformation R acting on the vector space V. Its action on basis vectors is given by [{\bf e}_i^\prime = \textstyle\sum\limits_j R_{ji} {\bf e}_j. \eqno(1.10.4.19)]If the basis is orthonormal, the matrix [R_{ij}] is orthogonal ([RR^{T}=E]). For a point group in superspace the action of R in [V_E] differs, in general, from that on [V_I]. [{\bf e}_{Ei}^\prime=\textstyle\sum\limits_j R_{Eji}{\bf e}_{Ej};\quad{\bf e}_{Ii}^\prime=\textstyle\sum\limits_j R_{Iji} {\bf e}_{Ij}. \eqno(1.10.4.20)]The action of R on the tensor product space [V_1 \otimes V_2], with [V_i] either [V_E] or [V_I], is given by [{\bf e}_{1i}^\prime \otimes {\bf e}_{2j}^\prime = \textstyle\sum\limits_k \textstyle\sum\limits_l R_{ki}^1R_{lj}^2{\bf e}_{1k}\otimes{\bf e}_{2l}. \eqno(1.10.4.21)]If both [R_i] are orthogonal matrices, the tensor product is also orthogonal. For the symmetrized tensor square [(V \otimes V)_{\rm sym}] the basis formed by [{\bf e}_{i}\otimes {\bf e}_i] ([i=1,\ldots ]) and [({\bf e}_{i}\otimes {\bf e}_j + {\bf e}_j \otimes {\bf e}_i)/\sqrt{2}] ([i \,\lt\, j]) is orthogonal.

A vector [\textstyle\sum\limits_{ij}c_{ij}{\bf e}_i\otimes {\bf e}_j] in the tensor product space is invariant if [RcR^{T}=c. \eqno(1.10.4.22)]A tensor as a (possibly symmetric or antisymmetric) bilinear function with coefficients [f_{ij}=f({\bf e}_i,{\bf e}_j)] is invariant if the matrix [f_{ij}] satisfies [R^{T}fR=f. \eqno(1.10.4.23)]For orthogonal bases the equations (1.10.4.22)[link] and (1.10.4.23)[link] are equivalent.

Which spaces have to be chosen for [V_i] depends on the physical tensor property. The algorithm for determining invariant tensors starts from the transformation of the basis vectors [{\bf e}_i], from which the basis transformation in tensor space follows after due orthogonalization in the case of (anti)symmetric tensors. This procedure can be continued to obtain higher-rank tensors. For orthogonal bases the invariant subspace is spanned by vectors corresponding to the independent tensor elements. We give a number of examples below.

1.10.4.6.1. Metric tensor for an octagonal three-dimensional quasicrystal

| top | pdf |

From the Fourier module for an octagonal quasicrystal in 3D the generators of the point group can be expressed as 5D integer matrices. They are [\displaylines{A=8(8^3)=\pmatrix{0&0&0&-1&0\cr 1&0&0&0&0\cr 0&1&0&0&0\cr 0&0&1&0&0\cr 0&0&0&0&1\cr}\cr B=m_z (1)=\pmatrix{1&0&0&0&0\cr 0&1&0&0&0\cr 0&0&1&0&0\cr 0&0&0&1&0\cr 0&0&0&0&-1\cr}}] and [C=m(m)=\pmatrix{-1&0&0&0&0\cr 0&0&0&1&0\cr 0&0&1&0&0\cr 0&1&0&0&0\cr 0&0&0&0&1\cr},]and span an integer representation of the point group 8/mmm(831mm). Solution of the three simultaneous equations [S^{T}g S=g] is equivalent with the determination of the subspace of the 15D symmetric tensor space that is invariant under the point group. The space has as basis the elements [e_{ij}] with [i\leq j]. The solution is given by [g = \pmatrix{g_{11}&g_{12}&0&-g_{12}&0\cr g_{12}&g_{11}&g_{12}&0&0\cr 0&g_{12}&g_{11}&g_{12}&0\cr -g_{12}&0&g_{12}&g_{11}&0\cr 0&0&0&0&g_{55}}.]If [{\bf e}_i\otimes{\bf e}_j] is denoted by ij, the solution follows because 55 is left invariant by A, B and C, whereas the orbits of 11 and 12 are [11\rightarrow 22 \rightarrow 33 \rightarrow 44 \rightarrow 11] and [12 \rightarrow 23 \rightarrow 34 \rightarrow -14 \rightarrow 12], respectively.

1.10.4.6.2. EFG tensor for Pcmn

| top | pdf |

The electric field gradient tensor transforms as the product of a reciprocal vector and a vector. In Cartesian coordinates the transformation properties are the same. The point group for the basic structure of many IC phases of the family of A2BX4 compounds is mmm, and the point group for the modulated phase is the 4D group [mmm(11\bar{1})], with generators [\displaylines{A = \pmatrix{-1&0&0&0\cr0&1&0&0\cr0&0&1&0\cr0&0&0&1\cr},\quad B = \pmatrix{1&0&0&0\cr 0&-1&0&0\cr 0&0&1&0 \cr 0&0&0&1\cr},\cr C = \pmatrix{1&0&0&0\cr0&1&0&0\cr 0&0&-1&0\cr0&0&0&-1\cr}.}]The tensor elements [\partial_i E_j] being indicated by ij, the transformation under the generators gives a factor [\pm 1] as shown in Table 1.10.4.2[link].

Table 1.10.4.2| top | pdf |
Sign change of [\partial_iE_j] under the generators A, B, C

 ABC
11+++
12+
13+
21+
22+++
23+
31+
32+
33+++
41+
42+
43+++

From this, it follows that the four independent tensor elements are [\partial_1 E_1], [\partial_2 E_2], [\partial_3 E_3] and the phason part [\partial_4 E_3].

1.10.4.6.3. Elasticity tensor for a two-dimensional octagonal quasicrystal

| top | pdf |

The point group of the standard octagonal tiling is generated by the 2D orthogonal matrices [A = \pmatrix{\sqrt{1/2}&-\sqrt{1/2}\cr \sqrt{1/2}&\sqrt{1/2}\cr},\quad B = \pmatrix{-1&0\cr 0&1\cr}. ]In the tensor space one has the following transformations of the basis vectors; they are denoted by ij for [{\bf e}_i\otimes{\bf e}_j]: [\eqalign{11 &\rightarrow {\textstyle{1\over2}}(11+12+21+22)\cr 12 &\rightarrow {\textstyle{1\over2}}(-11+12-21+22)\cr 21 &\rightarrow {\textstyle{1\over2}}(-11 - 12+21+22)\cr 22 &\rightarrow {\textstyle{1\over2}}(11 -12-21+22).\cr}]In the space spanned by a = 11, [b=\sqrt{1/2}(12+21)] and c = 22, the eightfold rotation is represented by the matrix [S_E = \pmatrix{{1\over2}&-\sqrt{1/2}&{1\over2} \cr \sqrt{1/2}&0&-\sqrt{1/2} \cr {1\over2}&\sqrt{1/2}&{1\over2} \cr}.]In the six-dimensional space with basis aa, [\sqrt{1/2}(ab+ba)], [\sqrt{1/2}(ac+ca)], bb, [\sqrt{1/2}(bc+cb)] and cc, the rotation gives the transformation [\pmatrix{{1\over4}&{1\over2}&\sqrt{2}/4&{1\over2}&{1\over2}&{1\over4}\cr -{1\over2} &-{1\over2}&0&0&{1\over2}&{1\over2}\cr \sqrt{2}/4&0&{1\over2}&-\sqrt{1/2}&0&\sqrt{2}/4\cr {1\over2}&0&-\sqrt{1/2}&0&0&{1\over2}\cr -{1\over2}&{1\over2}&0&0&-{1\over2}&{1\over2}\cr {1\over4}&-{1\over2}&\sqrt{2}/4&{1\over2}&-{1\over2}&{1\over4}\cr}.]The vector [{\bf v}] such that [{\bf S}\cdot {\bf v}={\bf v}] then is of the form [ {\bf v} = (v_1, 0, (v_1 -v_4)\sqrt{2}, v_4, 0, v_1)^{T}. ]This vector is also invariant under the mirror B. This means that there are two independent phonon elastic constants [c_{1111}^{E}] and [c_{1212}^{E}], whereas the other tensor elements satisfy the relations [\displaylines{c_{1112}^{E} = c_{1222}^{E} = 0,\quad c_{2222}^{E} = c_{1111}^{E}, \cr c_{1122}^{E} = (c_{1111}^{E} + c_{1212}^{E})\sqrt{1/2}.}]The internal component of the eightfold rotation is A3, that of the mirror B is B itself. The phason strain tensor transforms with the tensor product of external and internal components. This implies that the basis vectors, denoted by ij (i = 1, 2; j = 3, 4), transform under the eightfold rotation according to [\eqalign{13 &\rightarrow (-13+14 - 23+24)/2\cr 14 &\rightarrow (-13 - 14-23-24)/2\cr 23 &\rightarrow (13-14-23+24)/2\cr 24 &\rightarrow (13+14-23-24)/2.\cr}]The symmetrized tensor square of this matrix gives the transformation in the space of phason–phason elasticity tensors, the direct product of the transformations in the 3D phonon strain space and the 4D phason strain space gives the transformation in the space of phonon–phason elasticity tensors. The first matrix is given by [{1\over16}\pmatrix{1&\sqrt{2}&-\sqrt{2}&-\sqrt{2}&1&-\sqrt{2}&-\sqrt{2}&1&\sqrt{2}&1\cr -\sqrt{2}&0&2&0&\sqrt{2}&0&-2&-\sqrt{2}&0&\sqrt{2}\cr \sqrt{2}&2&0&0&\sqrt{2}&0&0&-\sqrt{2}&-2&-\sqrt{2}\cr -\sqrt{2}&0&0&2&\sqrt{2}&-2&0&\sqrt{2}&0&-\sqrt{2}\cr 1&-\sqrt{2}&-\sqrt{2}&\sqrt{2}&1&\sqrt{2}&-\sqrt{2}&1&-\sqrt{2}&1\cr -\sqrt{2}&0&0&-2&\sqrt{2}&2&0&\sqrt{2}&0&-\sqrt{2}\cr \sqrt{2}&-2&0&0&\sqrt{2}&0&0&-\sqrt{2}&2&-\sqrt{2}\cr 1&\sqrt{2}&\sqrt{2}&\sqrt{2}&1&\sqrt{2}&\sqrt{2}&1&\sqrt{2}&1\cr-\sqrt{2}&0&-2&0&\sqrt{2}&0&2&-\sqrt{2}&0&\sqrt{2}\cr 1&-\sqrt{2}&\sqrt{2}&-\sqrt{2}&1&-\sqrt{2}&\sqrt{2}&1&-\sqrt{2}&1\cr}. ]Vectors invariant under this operation and the transformation corresponding to the mirror B correspond to invariant elasticity tensors. For the transformation B, all tensor elements with an odd number of indices 1 or 3 are zero. In the space of phason–phason tensors the general invariant vector is [ (x_1,0,0,-x_6+(x_5 -x_1)\sqrt{2},x_5,x_6,0,x_5,0,x_1). ]

There are three independent elastic constants, [x_1=c_{1313}], [x_5=c_{1414}] and [x_6=c_{1423}]. For the phonon–phason elastic constants the corresponding invariant vector is [(x,0,0,x,0,x/\sqrt{2},-x/\sqrt{2},0,-x,0,0,-x). ] The independent elastic constant is x = [c_{1113}] = [c_{1124}] = [c_{1214}\sqrt{2}] = [-c_{1223}\sqrt{2}] = [-c_{2213}] = [-c_{2224}].

1.10.4.6.4. Piezoelectric tensor for a three-dimensional octagonal quasicrystal

| top | pdf |

A quasicrystal with octagonal point group 8/mmm(831mm) will not show a piezoelectric effect because the point group contains the central inversion. We consider here the point group 8mm(83mm) which is a subgroup without central inversion. It is generated by the matrices [\displaylines{A=\pmatrix{\alpha&-\alpha&0&0&0\cr \alpha&\alpha&0&0&0\cr 0&0&1&0&0 \cr 0&0&0&-\alpha&-\alpha \cr0&0&0&\alpha&-\alpha\cr},\cr B=\pmatrix{-1&0&0&0&0 \cr 0&1&0&0&0 \cr 0&0&1&0&0 \cr 0&0&0&-1&0 \cr 0&0&0&0&1\cr}.}]Here [\alpha =\sqrt{2}/2]. There are two components for the strain, a phonon component e and a phason component f. The phonon strain tensors form a 6D space, the phason strain tensors also a 6D space. The phonon strain space transforms with the symmetrized square of the physical parts of the operations, the phason strain space with the product of physical and internal parts. For the eightfold rotation the corresponding matrices are [\displaylines{S_e = \pmatrix{{1\over2}&-\alpha&0&{1\over2}&0&0\cr \alpha&0&0&-\alpha&0&0\cr 0&0&\alpha&0&-\alpha&0\cr {1\over2}&\alpha&0&{1\over2}&0&0\cr 0&0&\alpha&0&\alpha&0 \cr 0&0&0&0&0&1\cr}, \cr S_f = \pmatrix{-{1\over2}&-{1\over2}&{1\over2}&{1\over2}&0&0 \cr {1\over2}&-{1\over2}&-{1\over2}&{1\over2}&0&0\cr -{1\over2}&-{1\over2}&-{1\over2}&-{1\over2}&0&0\cr {1\over2}&-{1\over2}&{1\over2}&-{1\over2}&0&0\cr 0&0&0&0&-\alpha&-\alpha\cr 0&0&0&0&\alpha&-\alpha\cr}.}]The action on the space of piezoelectric tensors is given, for the phonon and the phason part, by taking the product of these matrices with the physical part [A_E]. The invariant vectors under these matrices give the invariant tensors. If the second generator is taken into account, which requires that the number of indices 1 or 4 is even, this results in the independent tensor elements [x_3 = c_{113}, \quad x_{16} = c_{322}, \quad x_{18} = c_{333}]with relation [c_{223} = c_{113}], whereas all other elements are zero for the coupling between the electric field and phonon strain. There is no nontrivial invariant vector in the second case. Therefore, all tensor elements for the coupling between the electric field and phason strain are zero.

1.10.4.6.5. Elasticity tensor for an icosahedral quasicrystal

| top | pdf |

The point group of an icosahedral quasicrystal is 532(5232) with generators having components [\displaylines{A_E =\pmatrix{1&\tau&-1-\tau\cr\tau&1+\tau&1\cr1+\tau&-1&\tau\cr} /2, \cr B_E = \pmatrix{-\tau&1+\tau&-1\cr1+\tau&1&\tau\cr 1&-\tau&-1-\tau\cr}/2 }]in physical space and components [A_I =\pmatrix{-1&\tau&-1-\tau\cr -\tau&\tau^{-1}&1\cr1+\tau&1&-\tau\cr} /2,\quad B_I = \pmatrix{0&0&-1\cr-1&0&0 \cr0&1&0\cr}]in internal space (see Table 1.10.5.2[link]). The phonon and phason strain tensors form a 6D, respectively 9D, vector space, in which the point group acts with matrices[\displaylines{S_e ={1\over4}\pmatrix{1&\tau\sqrt{2}&-\varphi\sqrt{2}&1-\tau&-\sqrt{2}&2+\tau\cr \tau\sqrt{2}&2&0&\sqrt{2}&-2&-\varphi\sqrt{2}\cr \varphi\sqrt{2}&0&-2&-\tau\sqrt{2}&2&-\sqrt{2}\cr 1-\tau&\sqrt{2}&\tau\sqrt{2}&2+\tau&\varphi\sqrt{2}&1\cr \sqrt{2}&2&2&-\varphi\sqrt{2}&0&\tau\sqrt{2}\cr 2+\tau &-\varphi\sqrt{2} &\sqrt{2}&1& -\tau\sqrt{2}& 1-\tau\cr},\cr S_f ={1\over4}\pmatrix{-\tau&\varphi&-1&\tau -1&1&-\tau &1&-2-\tau&\varphi\cr \varphi &1&\tau&1&\tau&1-\tau&-2-\tau&-\varphi&-1\cr\ 1&-\tau&-\varphi&\tau &\tau -1&-1&-\varphi&1&2+\tau\cr \tau -1&1&-\tau &-1&2+\tau&-\varphi&-\tau&\varphi&-1\cr 1&\tau&1-\tau& 2+\tau&\varphi&1&\varphi&1&\tau\cr \tau&\tau -1&-1&\varphi &-1&-2-\tau & 1&-\tau&-\varphi\cr -1&2+\tau&-\varphi&\tau&-\varphi&1&\tau -1&1&-\tau\cr 2+\tau &\varphi &1 & -\varphi& -1& -\tau&1 &\tau&1-\tau\cr \varphi& -1&-2-\tau&-1&\tau&\varphi& \tau & \tau -1& -1\cr}}]and[\displaylines{T_e ={1\over4}\pmatrix{1&-\tau\sqrt{2}&\varphi\sqrt{2}&1-\tau&-\sqrt{2}&2+\tau\cr \tau\sqrt{2}&-2&0&\sqrt{2} &-2&-\varphi\sqrt{2}\cr -\varphi\sqrt{2} & 0& -2& \tau\sqrt{2}& -2& \sqrt{2}\cr 1-\tau& -\sqrt{2}& -\tau\sqrt{2} & 2+\tau& \varphi\sqrt{2}&1\cr -\sqrt{2}& 2&2&\varphi\sqrt{2}&0&-\tau\sqrt{2}\cr 2+\tau&\varphi\sqrt{2}&-\sqrt{2}&1& -\tau\sqrt{2}&1-\tau\cr},\cr T_f ={1\over2}\pmatrix{0&0&1&0&0&-\tau&0&0&\varphi\cr 1&0&0&-\tau&0&0&\varphi&0&0\cr 0&-1&0&0&\tau&0&0&-\varphi&0\cr 0&0&\tau&0&0&-\varphi&0&0&-1\cr \tau&0&0&-\varphi&0&0&-1&0&0\cr 0&-\tau&0&0&\varphi&0&0&1&0\cr 0&0&-\varphi&0&0&-1&0&0&\tau\cr -\varphi&0&0&-1&0&0&\tau&0&0\cr 0&\varphi&0&0&1&0&0&-\tau&0\cr}.}]This implies that the phonon elasticity tensors form a 21D space, the phason elasticity tensors a 45D space and the phonon–phason coupling a 54D space. The invariant vectors under these orthogonal transformations correspond to invariant elastic tensors. Their coordinates are the elastic constants. For the given presentation of the point group, these are given in Table 1.10.4.3[link]. The tensor elements are expressed in parameters x and y where there are two independent tensor elements. The tensor elements that are not given are zero or equal to that given by the permutation symmetry. If bases for the phonon and phason strain are introduced by [[1]=11, [2]=12, [3]=13, [4]=22, [5]=23, [6]=33]for the phonon part and [\displaylines{[1]=14, [2]=15, [3]=16, [4]=24, [5]=25, [6]=26, \cr [7]=34, [8]=35, [9]=36}]for the phason part, the elastic tensors may be given in matrix form as[\displaylines{c^{ee}=\pmatrix{x+y&0&0&x&0&x\cr 0&y&0&0&0&0\cr 0&0&y&0&0&0\cr x&0&0&x+y&0&x\cr 0&0&0&0&y&0\cr x&0&0&x&0&x+y\cr}, \cr c^{ff}=\pmatrix{z&\tau^2 u&- \tau u& - \tau u& \tau u&- \tau u&-u&0& \tau u \cr \tau^2 u&z- 2\tau u&u& \tau u&u&0&0&\tau^2 u& \tau u \cr - \tau u &u& z&- \tau u&0&-\tau^2 u& \tau u& \tau u& \tau u\cr - \tau u& \tau u&-\tau u&z& \tau u&u&-\tau^2 u& \tau u&0\cr \tau u&u&0& \tau u&z&-\tau^2 u& \tau u&- \tau u&- \tau u\cr - \tau u&0&-\tau^2 u&u&-\tau^2 u&z-2 \tau u&0&- \tau u&u\cr -u&0& \tau u&-\tau^2 u& \tau u&0&z-2 \tau u&-u&-\tau^2 u\cr 0&\tau^2 u& \tau u& \tau u&- \tau u&- \tau u&-u&z&- \tau u\cr \tau u& \tau u& \tau u&0&- \tau u&u&-\tau^2 u&- \tau u &z\cr}, \cr c^{ef}=\pmatrix{-v&- \tau v&-\tau^2 v&-\tau^3v& \tau v&-\tau^2 v&v&- \tau v&-\tau^{-1}v\cr -\tau^3v& \tau v&-\tau^2 v&\tau^2v&v& \tau v&0&0&0\cr v&- \tau v&-\tau^{-1}v&0&0&0&- \tau v&-\tau^2 v&-\tau^3 v\cr \tau^{-1}v&v& \tau v&-\tau^2 v&v&- \tau v&-\tau^2v&-\tau^3 v& \tau v\cr\ 0&0&0&-\tau^2 v&\tau^3 v& \tau v& \tau v&-\tau^{-1} v&v\cr -\tau v&-\tau^2v&-\tau^3v& \tau v&-\tau^{-1} v&v&- \tau v&\tau^2 v&v \cr}.}]

Table 1.10.4.3| top | pdf |
Elastic constants for icosahedral quasicrystals

TypeFree parametersRelations
Phonon–phonon2[ c_{1111}=c_{2222}=c_{3333}=x]
  [c_{1122}=c_{1133}=c_{2233}=y]
  [c_{1212}=c_{1313}=c_{2323}=x-y]
Phason–phason2[c_{1414}=c_{1616}=c_{2424}=c_{2525}=c_{3535}=c_{3636}=x]
  [c_{1416}=c_{1424}=-c_{1425}=c_{1426}=-c_{1436}=y]
  [c_{1524}=c_{1536}=-c_{1624}=c_{1634}=c_{1635}=-y]
  [c_{1636}=c_{2425}=c_{2435}=c_{2534}=-c_{2535}=-y]
  [c_{2536}=c_{2635}=c_{3536}=y, c_{1515}=c_{2626}=c_{3434}=x+2y]
  [c_{1415}=c_{1535}=-c_{1626}=-c_{2434}=-c_{2526}=-c_{3436}=y\tau]
  [c_{1434}=-c_{11516}=-c_{1525}=-c_{2426}=-c_{2636}=c_{3435}=y/\tau]
Phonon–phason1[-c_{1114}=c_{1134}=c_{1225}=c_{2225}=c_{2336}=c_{3326}=c_{3336}]
  [\quad=-c_{1115}/\tau =c_{1125}/\tau =-c_{1135}\tau=c_{1215}\tau =c_{1226}\tau]
  [\quad=-c_{1334}\tau =-c_{2226}\tau =c_{2236}\tau =c_{2326}\tau =c_{2334}\tau ]
  [\quad=-c_{3334}\tau =-c_{1116}\tau^2 =-c_{1126}\tau^2 =-c_{1216}\tau^2 ]
  [\quad=-c_{1335}\tau^2 =-c_{2224}\tau^2 =-c_{2234}\tau^2 =c_{3335}\tau^2 ]
  [\quad=-c_{1136}\tau =c_{1224}\tau =-c_{1316}\tau =-c_{2335}\tau ]
  [\quad=-c_{1124}/\tau^3 =-c_{1336}/\tau^3 =c_{2235}/\tau^3 =c_{2325}/\tau^3 ]

The parameters x, y, z, u, v are the five independent elastic constants.

1.10.5. Tables

| top | pdf |

In this section are presented the irreducible representations of point groups of quasiperiodic structures up to rank six that do not occur as three-dimensional crystallographic point groups.

Table 1.10.5.1[link] gives the characters of the point groups [C_n] with n = 5, 8, 10, 12, [D_n] with n = 5, 8, 10, 12, and the icosahedral group I. The direct products with [{\bb Z}_2] then follow easily. Although these direct products of a group K with [{\bb Z}_2] do not belong to the isomorphism class of K, their irreducible representations are nevertheless given in the table for K because these irreducible representations have the same labels as those for K apart from an additional subindex u. The reresentations of the subgroup K of [K\times {\bb Z}_2] are the same as for K itself, those for the cosets get an additional minus sign. In the tables, the characters for the groups [K\times {\bb Z}_2] are separated from those for K by a horizontal rule. In addition to the characters are given the realizations of crystallographic point groups, and the irreducible components of the vector representations in direct space [V_E] and internal space [V_I] for these realizations. The vector representation in [V_I] is called the perpendicular representation.

Table 1.10.5.1| top | pdf |
Character tables of some point groups for quasicrystals

(a) [C_5] [[\omega =\exp (2\pi i/5)]].

[C_{5}][\varepsilon][\alpha][\alpha^{2}][\alpha^{3}][\alpha^{4}]
n11111
Order 15555
[\Gamma_{1}]11111
[\Gamma_{2}]1[\omega][\omega^{2}][\omega^{3}][\omega^{4}]
[\Gamma_{3}]1[\omega^{2}][\omega^{4}][\omega][\omega^{3}]
[\Gamma_{4}]1[\omega^{3}][\omega][\omega^{4}][\omega^{2}]
[\Gamma_{5}]1[\omega^{4}][\omega^{3}][\omega^{2}][\omega]

 Generators Vector representationPerpendicular representation
5[\alpha = C_{5z}] [\Gamma_{1}\oplus\Gamma_{2}\oplus\Gamma_{5}][\Gamma_{3}\oplus\Gamma_{4}]

(b) [D_5] [[\tau = (\sqrt{5}-1)/2]].

[D_{5}][\varepsilon][\alpha][\alpha^{2}] [\beta]
n1225
Order 1 5 5 2
[\Gamma_{1}]1 1 1 1
[\Gamma_{2}]1 1 1 [-1]
[\Gamma_{3}]2 [\tau] [-1-\tau] 0
[\Gamma_{4}]2 [-1-\tau] [\tau] 0

 Generators Vector representationPerpendicular representation
52 [\alpha = C_{5z}] [\Gamma_{2}\oplus \Gamma_{3}][\Gamma_{4}]
 [\beta = C_{2x}]  
5m [\alpha = C_{5z}] [\Gamma_{1}\oplus \Gamma_{3}][\Gamma_{4}]
 [\beta = m_{x}]   
[\bar{5}m] [\sim 52 \times {\bb Z}_{2}] [\Gamma_{1u}\oplus\Gamma_{3u}][\Gamma_{4u}]

(c) [C_8] [[\omega =\exp (\pi i/4) = (1+i)/\sqrt{2}]].

[C_{8}][\varepsilon][\alpha][\alpha^{2}][\alpha^{3}][\alpha^{4}][\alpha^{5}][\alpha^{6}][\alpha^{7}]
[ n]11111111
Order 18482868
[\Gamma_{1}]11111111
[\Gamma_{2}]1[\omega]i[\omega^{3}][-1] [\omega^{5}][-i][\omega^{7}]
[\Gamma_{3}]1i[-1][-i][1] i[-1][-i]
[\Gamma_{4}]1[\omega^{3}][-i][\omega][-1] [\omega^{7}]i[\omega^{5}]
[\Gamma_{5}]1[-1][1][-1][1] [-1][1][-1]
[\Gamma_{6}]1[\omega^{5}]i[\omega^{7}][-1] [\omega][-i][\omega^{3}]
[\Gamma_{7}]1[-i][-1]i[1] [-i][-1]i
[\Gamma_{8}]1[\omega^{7}][-i][\omega^{5}][-1] [\omega^{3}]i[\omega]

 Generators Vector representationPerpendicular representation
8[\alpha = C_{8z}] [\Gamma_{1}\oplus\Gamma_{2}\oplus\Gamma_{8}] [\Gamma_{4}\oplus\Gamma_{6}]
[\bar{8}] [\alpha =S_{8z}] [\Gamma_4\oplus\Gamma_5\oplus\Gamma_6] [\Gamma_2\oplus\Gamma_8]
[8/m] [\sim 8\times {\bb Z}_2] [\Gamma_{1u}\oplus\Gamma_{2u}\oplus\Gamma_{8u}] [\Gamma_{4u}\oplus\Gamma_{6u}]

(d) [D_8]

[D_{8}][\varepsilon][\alpha] [\alpha^{2}][\alpha^{3}][\alpha^{4}][\beta][\alpha \beta]
[ n]1222144
Order 1 8 4 8 2 2 2
[\Gamma_{1}]1 1 1 1 1 1 1
[\Gamma_{2}]1 1 1 1 1 [-1] [-1]
[\Gamma_{3}]1 [-1] 1 [-1 ]1 1 [-1]
[\Gamma_{4}]1 [-1 ]1 [-1] 1 [-1] 1
[\Gamma_{5}]2 [\sqrt{2}] 0 [-\sqrt{2}] [-2] 0 0
[\Gamma_{6}]2 0 [-2 ]0 2 0 0
[\Gamma_{7}]2 [-\sqrt{2}] 0 [\sqrt{2}] [-2] 0 0

 GeneratorsVector representationPerpendicular representation
822[\alpha = C_{8z}] [\Gamma_{2}\oplus\Gamma_{5}][\Gamma_{7}]
 [\beta = C_{2x}]  
[8mm][\alpha = C_{8z}] [\Gamma_{1}\oplus\Gamma_{5}][\Gamma_{7}]
 [\beta = m_{x}]  
[\bar{8}2m][\alpha = S_{8z}] [\Gamma_{3}\oplus\Gamma_{7}][\Gamma_{5}]
 [\beta = C_{2x}]  
[8/mmm][\sim 822\times{\bb Z}_{2}][\Gamma_{2u}\oplus\Gamma_{5u}][\Gamma_{7u}]

(e) [C_{10}] [[\omega=exp(2\pi i/5) ]].

[C_{10}][ \varepsilon][\alpha^2][\alpha^4][\alpha^6][\alpha^8]
n11111
Order 15555
[\Gamma_1]11111
[\Gamma_2]1[\omega][\omega^2][\omega^3][\omega^4]
[\Gamma_3]1[\omega^2][\omega^4][\omega][\omega^3]
[\Gamma_4]1[\omega^3][\omega][\omega^4][\omega^2]
[\Gamma_5]1[\omega^4][\omega^3][\omega^2][\omega]
[\Gamma_6]11111
[\Gamma_7]1[\omega][\omega^2][\omega^3][\omega^4]
[\Gamma_8]1[\omega^2][\omega^4][\omega][\omega^3]
[\Gamma_9]1[\omega^3][\omega][\omega^4][\omega^2]
[\Gamma_{10}]1[\omega^4][\omega^3][\omega^2][\omega]

[C_{10}][\alpha^5][\alpha^7][\alpha^9][\alpha][\alpha^3]
n11111
Order 210101010
[\Gamma_1]11111
[\Gamma_2]1[\omega][\omega^2][\omega^3][\omega^4]
[\Gamma_3]1[\omega^2][\omega^4][\omega][\omega^3]
[\Gamma_4]1[\omega^3][\omega][\omega^4][\omega^2]
[\Gamma_5]1[\omega^4][\omega^3][\omega^2][\omega]
[\Gamma_6][-1][-1][-1][-1][-1]
[\Gamma_7][-1][-\omega][-\omega^2][-\omega^3][-\omega^4]
[\Gamma_8][-1][-\omega^2][-\omega^4][-\omega][-\omega^3]
[\Gamma_9][-1][-\omega^3][-\omega][-\omega^4][-\omega^2]
[\Gamma_{10}][-1][-\omega^4][-\omega^3][-\omega^2][-\omega]

 Generators Vector representationPerpendicular representation
10[\alpha = C_{10z}] [\Gamma_{1}\oplus\Gamma_{7}\oplus\Gamma_{10}][\Gamma_{8}\oplus\Gamma_{9}]
[{\bar{5}}][\alpha = S_{5z}] [\Gamma_6\oplus\Gamma_8\oplus\Gamma_9] [\Gamma_7\oplus\Gamma_{10}]
[\overline{10}][\alpha =S_{10z}] [\Gamma_2\oplus\Gamma_4\oplus\Gamma_6] [\Gamma_3\oplus\Gamma_5]
[10/m ][\sim 10\times {\bb Z}_2] [\Gamma_{1u}\oplus\Gamma_{7u}\oplus\Gamma_{10u}][\Gamma_{8u}\oplus\Gamma_{9u}]

(f) [D_{10}] [[\tau = (\sqrt{5}-1)/2]].

[D_{10}] [\varepsilon] [\alpha] [\alpha^{2}] [\alpha^{3}]
n1222
Order 1 10 5 10
[\Gamma_{1}] 1 1 1 1
[\Gamma_{2}] 1 1 1 1
[\Gamma_{3}] 1 [-1 ]1 [-1]
[\Gamma_{4}] 1 [-1] 1 [-1]
[\Gamma_{5}]2 [1+\tau] [\tau] [-\tau]
[\Gamma_{6}]2 [\tau] [-1-\tau] [-1-\tau]
[\Gamma_{7}]2 [-\tau] [-1-\tau] [1+\tau]
[\Gamma_{8}]2 [-1-\tau] [\tau] [\tau]

[D_{10}][\alpha^{4}] [\alpha^{5}] [\beta] [\alpha \beta]
n2155
Order 5 2 2 2
[\Gamma_{1}] 1 1 1 1
[\Gamma_{2}] 1 1 [-1 ][-1]
[\Gamma_{3}] 1 [-1] 1 [-1]
[\Gamma_{4}] 1 [-1 ][-1] 1
[\Gamma_{5}][-1-\tau] [-2 ]0 0
[\Gamma_{6}][\tau] 2 0 0
[\Gamma_{7}][\tau ] [-2 ]0 0
[\Gamma_{8}][-1-\tau] 2 0 0

 GeneratorsVector representationPerpendicular representation
[10\,22][\alpha = C_{10z}] [\Gamma_{2}\oplus\Gamma_{5}][\Gamma_7]
 [\beta = C_{2x}]  
[10\,mm][\alpha = C_{10z}] [\Gamma_{1}\oplus\Gamma_{5}][\Gamma_7]
 [\beta = m_{x}]  
[\overline{10}\,2m] [\alpha = S_{10z}] [\Gamma_{4}\oplus\Gamma_{8}][\Gamma_6]
 [\beta = C_{2x}]  
[\bar{5}m] [\alpha = S_{5z}] [\Gamma_{3}\oplus\Gamma_{7}][\Gamma_5]
 [\beta = C_{2x}]   
[10/mmm][\sim 10\,22\times{\bb Z}_{2}][\Gamma_{2u}\oplus\Gamma_{5u}][\Gamma_{7u}]

(g) [C_{12}] [[\omega = \exp (\pi i/6)]].

[C_{12}][\varepsilon] [\alpha][\alpha^{2}][\alpha^{3}][\alpha^4][\alpha^{5}]
n111111
Order11264312
[\Gamma_{1}]111111
[\Gamma_{2}]1[\omega][\omega^2]i[\omega^4][\omega^5]
[\Gamma_{3}]1[\omega^2][\omega^4][-1][-\omega^2][-\omega^4]
[\Gamma_4]1i[-1][-i]1i
[\Gamma_{5}]1[\omega^4][-\omega^2][1][\omega^4][-\omega^2]
[\Gamma_{6}]1[\omega^5][-\omega^4]i[-\omega^2][\omega]
[\Gamma_{7}]1[-1]1[-1]1[-1]
[\Gamma_{8}]1[-\omega][\omega^2][-i][\omega^4][-\omega^5]
[\Gamma_{9}]1[-\omega^2][\omega^4][1][-\omega^2][\omega^4]
[\Gamma_{10}]1[-i][-1]i1[-i]
[\Gamma_{11}]1[-\omega^4][-\omega^2][-1][\omega^4][\omega^2]
[\Gamma_{12}]1[-\omega^5][-\omega^4][-i][-\omega^2][-\omega]

[C_{12}][\alpha^{6}][\alpha^{7}][\alpha^{8}][\alpha^{9}][\alpha^{10}][\alpha^{11}]
n111111
Order21234612
[\Gamma_{1}]11 1111
[\Gamma_{2}][-1][-\omega][-\omega^2][-i][-\omega^4][-\omega^5]
[\Gamma_{3}]1[\omega^2][\omega^4][-1][-\omega^2][-\omega^4]
[\Gamma_4][-1][-i]1i[-1][-i]
[\Gamma_{5}]1[\omega^4][-\omega^2][1][\omega^4][-\omega^2]
[\Gamma_{6}][-1][-\omega^5][\omega^4][-i][\omega^2][-\omega]
[\Gamma_{7}]1[-1] 1[-1]1[-1]
[\Gamma_{8}][-1][\omega][-\omega^2]i[-\omega^4][\omega^5]
[\Gamma_{9}]1[-\omega^2][\omega^4][1][-\omega^2][\omega^4]
[\Gamma_{10}][-1]i1[-i][-1]i
[\Gamma_{11}]1[-\omega^4][-\omega^2][-1][\omega^4][\omega^2]
[\Gamma_{12}][-1][\omega^5][\omega^4]i[\omega^2][\omega]

 Generators Vector representationPerpendicular representation
12[\alpha=C_{12z}] [\Gamma_{1}\oplus\Gamma_{2}\oplus\Gamma_{12}][\Gamma_{6}\oplus\Gamma_{8}]
[\overline{12}] [\alpha=S_{12z}] [\Gamma_6\oplus\Gamma_7\oplus\Gamma_8] [\Gamma_2\oplus\Gamma_{12}]
[12/m ][\sim 12\times {\bb Z}_2] [\Gamma_{1u}\oplus\Gamma_{2u}\oplus\Gamma_{12u}][\Gamma_{6u}\oplus\Gamma_{8u}]

(h) [D_{12}]

[D_{12}][\varepsilon] [\alpha] [\alpha^{2}] [\alpha^{3}]
n1222
Order 1 12 6 4
[\Gamma_{1}]1 1 1 1
[\Gamma_{2}]1 1 1 1
[\Gamma_{3}]1 [-1] 1 [-1 ]
[\Gamma_{4}]1 [-1] 1 [-1]
[\Gamma_{5}]2 [\sqrt{3}] 1 0
[\Gamma_{6}]2 1 [-1 ][-2]
[\Gamma_{7}]2 0 [-2] 0
[\Gamma_{8}]2 [-1] [-1 ]2
[\Gamma_{9}]2 [-\sqrt{3}] 1 0

[D_{12}][\alpha^{4}] [\alpha^{5}] [\alpha^{6}] [\beta] [\alpha \beta]
n22166
Order 3 12 2 2 2
[\Gamma_{1}]1 1 1 1 1
[\Gamma_{2}]1 1 1 [-1] [-1]
[\Gamma_{3}]1 [-1] 1 1 [-1]
[\Gamma_{4}]1 [-1] 1 [-1] 1
[\Gamma_{5}][-1] [-\sqrt{3}] [-2 ]0 0
[\Gamma_{6}][-1] 1 2 0 0
[\Gamma_{7}]2 0 [-2 ]0 0
[\Gamma_{8}][-1 ][-1 ]2 0 0
[\Gamma_{9}][-1 ][\sqrt{3}] [-2] 0 0

 GeneratorsVector representationPerpendicular representation
[12 \,22][\alpha = C_{12z}] [\Gamma_{2}\oplus\Gamma_{5}][\Gamma_9]
 [\beta = C_{2x}]  
[12\,mm][\alpha = C_{12z}] [\Gamma_{1}\oplus\Gamma_{5}][\Gamma_{9}]
 [\beta = m_{x}]  
[\overline{12}\,2m][\alpha = S_{12z}] [\Gamma_{4}\oplus\Gamma_{9}][\Gamma_5]
 [\beta = C_{2x}]  
[12/mmm][\sim 12 \,22\times{\bb Z}_{2}][\Gamma_{2u}\oplus\Gamma_{5u}][\Gamma_{9u}]

(i) I [[\tau = (\sqrt 5 -1)/2]].

I[\varepsilon] [\alpha] [\alpha^{2}] [\beta] [\alpha \beta]
n112122015
Order 1 5 5 3 2
[\Gamma_{1}]1 1 1 1 1
[\Gamma_{2}]3 [1+\tau] [-\tau] 0 [-1]
[\Gamma_{3}]3 [-\tau] [1+\tau] 0 [-1]
[\Gamma_{4}]4 [-1] [-1] 1 0
[\Gamma_{5}]5 0 0 [-1] 1

 GeneratorsVector representationPerpendicular representation
532 [\alpha =C_{5}][\Gamma_2][\Gamma_3]
 [\beta =C_{3d}]  
[\bar{5}\bar{3}m][\sim 532 \times{\bb Z}_{2}][\Gamma_{2u}][\Gamma_{3u}]

In Table 1.10.5.2[link] the representation matrices for the irreducible representations in more than one dimension are given (one-dimensional representations are just the characters). For the cyclic groups there are only one-dimensional representations, for the dihedral groups there are one- and two-dimensional irreducible representations. There are four irreducible representations of I of dimension larger than one. The four- and five-dimensional ones are given as integer representations. They form crystallographic groups in 4D and 5D. The two three-dimensional representations have the same matrices. The elements, however, are connected by an outer automorphism. That means that the ith element [R_i] is represented by [\Gamma_2 (R_i)] in the representation [\Gamma_2], and by [\Gamma_3 (R_i)=\Gamma_2 (\varphi R_i)] in [\Gamma_3]. The element [\varphi R_i] is another element [R_j]. The corresponding j for each i is given in Table 1.10.5.3[link].

Table 1.10.5.2| top | pdf |
Matrices of the irreducible representations of dimension [d \geq 2] corresponding to the irreps of Table 1.10.5.1[link]

(a) [{D}_{5}]

Representation [D(\alpha^{p})][D(\beta)]
[\Gamma_{3}][\pmatrix{\cos (2\pi p/5)& -\sin (2\pi p/5)\cr \sin (2\pi p/5)& \cos (2\pi p/5)\cr}][\pmatrix{0 & 1 \cr 1 & 0 \cr}]
[\Gamma_{4}][\pmatrix{\cos (4\pi p/5)& -\sin (4\pi p/5)\cr \sin (4\pi p/5)& \cos (4\pi p/5)\cr}][\pmatrix{ 0 & 1 \cr 1 & 0\cr}]

(b) [{D}_{8}]

Representation [D(\alpha^{p})][D(\beta)]
[\Gamma_{5}][\pmatrix{\cos (\pi p/4)& -\sin (\pi p/4)\cr \sin (\pi p/4)& \cos (\pi p/4)\cr}][\pmatrix{0 & 1 \cr 1 & 0 \cr}]
[\Gamma_{6}][\pmatrix{\cos (\pi p/2)& -\sin (\pi p/2)\cr \sin (\pi p/2)& \cos (\pi p/2) \cr}][\pmatrix{0 & 1 \cr 1 & 0\cr}]
[\Gamma_{7}][\pmatrix{\cos (3\pi p/4)& -\sin (3\pi p/4)\cr \sin (3\pi p/4)& \cos (3\pi p/4)\cr}][\pmatrix{ 0 & 1 \cr1 & 0 \cr}]

(c) [{D}_{10}]

Representation [D(\alpha^{p})][D(\beta)]
[\Gamma_{5}][\pmatrix{\cos (\pi p/5)& -\sin (\pi p/5)\cr \sin (\pi p/5)& \cos (\pi p/5)\cr}][\pmatrix{0 & 1 \cr 1 & 0\cr}]
[\Gamma_{6}][\pmatrix{\cos (2\pi p/5)& -\sin (2\pi p/5)\cr \sin (2\pi p/5)& \cos (2\pi p/5)\cr}][\pmatrix{0 & 1 \cr 1 & 0\cr}]
[\Gamma_{7}][\pmatrix{\cos (3\pi p/5)& -\sin (3\pi p/5)\cr \sin (3\pi p/5)& \cos (3\pi p/5)\cr}][\pmatrix{0 & 1 \cr 1 & 0 \cr}]
[\Gamma_{8}][\pmatrix{\cos (4\pi p/5)& -\sin (4\pi p/5)\cr \sin (4\pi p/5)& \cos (4\pi p/5)\cr}][\pmatrix{0 & 1 \cr 1 & 0 \cr}]

(d) [{D}_{12}]

Representation [D(\alpha^{p})][D(\beta)]
[\Gamma_{5}][\pmatrix{\cos (\pi p/6)& -\sin (\pi p/6)\cr \sin (\pi p/6)& \cos (\pi p/6)\cr} ][\pmatrix{0 & 1 \cr 1 & 0\cr}]
[\Gamma_{6}][\pmatrix{\cos (\pi p/3)& -\sin (\pi p/3)\cr \sin (\pi p/3)& \cos (\pi p/3)\cr} ][\pmatrix{0 & 1 \cr 1 & 0\cr}]
[\Gamma_{7}][\pmatrix{\cos (\pi p/2)& -\sin (\pi p/2)\cr \sin (\pi p/2)& \cos (\pi p/2)\cr} ][\pmatrix{0 & 1 \cr 1 & 0 \cr}]
[\Gamma_{8}][\pmatrix{\cos (2\pi p/3)& -\sin (2\pi p/3)\cr \sin (2\pi p/3)& \cos (2\pi p/3)\cr}][\pmatrix{0 & 1 \cr 1 & 0\cr}]
[\Gamma_{9}][\pmatrix{\cos (5\pi p/6)& -\sin (5\pi p/6)\cr \sin (5\pi p/6)& \cos (5\pi p/6)\cr}][\pmatrix{0 & 1 \cr 1 & 0 \cr}]

(e) I. First column: numbering of the elements. [f=(1+\sqrt{5})/2, t=(\sqrt{5}-1)/2]. Horizontal rules separate conjugation classes.

No. Order[\Gamma_{2}][\Gamma_{4}][\Gamma_{5}]
1 1 [\pmatrix{1&0&0 \cr 0&1&0\cr 0&0&1 \cr}][\pmatrix{1&0&0&0\cr 0&1&0&0 \cr 0&0&1&0 \cr 0&0&0&1 \cr}][\pmatrix{1&0&0&0&0\cr 0&1&0&0&0\cr 0&0&1&0&0 \cr 0&0&0&1&0\cr 0&0&0&0&1 \cr}]
25[\pmatrix{1/2&t/2&-f/2\cr t/2&f/2&1/2\cr f/2&-1/2&t/2 \cr}][\pmatrix{0&0&0&-1\cr1&0&0&-1\cr 0&1&0&-1\cr 0&0&1&-1\cr}][\pmatrix{1&0&0&0&-1\cr0&0&0&0&-1\cr0&1&0&0&-1\cr0&0&1&0&-1\cr 0&0&0&1&-1\cr}]
3 5 [\pmatrix{1/2&-t/2&f/2\cr-t/2&f/2&1/2\cr\-f/2&-1/2&t/2\cr}][\pmatrix{0&0&1&-1\cr 1&0&0&-1\cr0&0&0&-1\cr0&1&0&-1\cr}][\pmatrix{0&-1&1&0&0\cr0&-1&0&0&1\cr0&-1&0&0&0\cr 0&-1&0&1&0\cr 1&-1&0&0&0\cr}]
4 5 [\pmatrix{1/2&t/2&f/2\cr t/2&f/2&-1/2\cr -f/2&1/2&t/2\cr}][\pmatrix{-1&1&0&0\cr -1&0&1&0\cr -1&0&0&1\cr -1&0&0&0\cr}][\pmatrix{1&-1&0&0&0\cr 0&-1&1&0&0\cr 0&-1&0&1&0\cr 0&-1&0&0&1\cr 0&-1&0&0&0 \cr}]
5 5 [\pmatrix{t/2&-f/2&1/2\cr f/2&1/2&t/2\cr -1/2&t/2&f/2\cr}][\pmatrix{0&-1&0&0\cr 0&-1&0&1\cr 1&-1&0&0\cr 0&-1&1&0\cr}][\pmatrix{0&1&0&-1&0\cr 0&0&0&-1&1\cr 0&0&1&-1&0\cr 1&0&0&-1&0\cr 0&0&0&-1&0\cr}]
6 5 [\pmatrix{ f/2&-1/2&-t/2\cr 1/2&t/2&f/2\cr -t/2&-f/2&1/2\cr}][\pmatrix{0&1&-1&0\cr 0&0&-1&0\cr 0&0&-1&1\cr 1&0&-1&0\cr} ][\pmatrix{0&1&0&0&0\cr 0&0&0&1&0\cr 0&0&0&0&1\cr 0&0&1&0&0\cr 1&0&0&0&0\cr}]
7 5 [\pmatrix{ f/2&1/2&t/2\cr -1/2&t/2&f/2\cr t/2&-f/2&1/2\cr}][\pmatrix{0&-1&0&1\cr 0&-1&1&0\cr 1&-1&0&0\cr 0&-1&0&0\cr}][\pmatrix{-1&0&1&0&0\cr -1&1&0&0&0\cr -1&0&0&0&1\cr -1&0&0&0&0\cr -1&0&0&1&0\cr}]
8 5 [\pmatrix{ t/2&f/2&-1/2\cr -f/2&1/2&t/2\cr 1/2&t/2&f/2\cr}][\pmatrix{-1&0&1&0\cr -1&0&0&0\cr -1&0&0&1\cr -1&1&0&0 \cr}][\pmatrix{0&0&0&1&-1\cr 1&0&0&0&-1\cr 0&0&1&0&-1\cr 0&0&0&0&-1\cr 0&1&0&0&-1\cr}]
95[\pmatrix{ t/2&f/2&1/2\cr -f/2&1/2&-t/2\cr -1/2&-t/2&f/2\cr}][\pmatrix{0&0&-1&0\cr 0&0&-1&1\cr 0&1&-1&0\cr 1&0&-1&0\cr}][\pmatrix{-1&0&0&1&0\cr -1&0&1&0&0\cr -1&0&0&0&0\cr -1&1&0&0&0\cr -1&0&0&0&1\cr}]
10 5 [\pmatrix{ f/2&1/2&-t/2\cr -1/2&t/2&-f/2\cr -t/2&f/2&1/2\cr}][\pmatrix{0&-1&0&1\cr 1&-1&0&0\cr 0&-1&0&0\cr 0&-1&1&0\cr}][\pmatrix{0&0&0&0&1\cr 1&0&0&0&0\cr 0&0&0&1&0\cr 0&1&0&0&0\cr 0&0&1&0&0\cr} ]
115[\pmatrix{ 1/2&-t/2&-f/2\cr -t/2&f/2&-1/2\cr f/2&1/2&t/2\cr}][\pmatrix{0&1&-1&0\cr 0&0&-1&1\cr 1&0&-1&0\cr 0&0&-1&0\cr}][\pmatrix{0&0&-1&0&1\cr 0&0&-1&0&0\cr 1&0&-1&0&0\cr 0&0&-1&1&0\cr 0&1&-1&0&0\cr}]
125 [\pmatrix{ f/2&-1/2&t/1\cr 1/2&t/2&-f/2\cr t/2&f/2&1/2\cr}][\pmatrix{0&0&1&-1\cr 0&0&0&-1\cr 0&1&0&-1\cr 1&0&0&-1\cr}][\pmatrix{0&0&0&-1&0\cr 0&1&0&-1&0\cr 1&0&0&-1&0\cr 0&0&0&-1&1\cr 0&0&1&-1&0\cr}]
13 5 [\pmatrix{ t/2&-f/2&-1/2\cr f/2&1/2&-t/2\cr 1/2&-t/2&f/2\cr} ][\pmatrix{-1&0&0&1\cr -1&0&1&0\cr -1&0&0&0\cr -1&1&0&0\cr}][\pmatrix{0&0&-1&0&0\cr 0&0&-1&1&0\cr 0&1&-1&0&0\cr 1&0&-1&0&0\cr 0&0&-1&0&1\cr}]
14 5 [\pmatrix{ -t/2&f/2&-1/2\cr f/2&1/2&t/2\cr 1/2&-t/2&-f/2\cr} ][\pmatrix{0&0&-1&1\cr 0&0&-1&0\cr 1&0&-1&0\cr 0&1&-1&0\cr}][\pmatrix{1&0&0&-1&0\cr 0&0&0&-1&1\cr 0&0&0&-1&0\cr 0&1&0&-1&0\cr 0&0&1&-1&0\cr}]
15 5 [\pmatrix{ -t/2&f/2&1/2\cr f/2&1/2&-t/2\cr -1/2&t/2&-f/2\cr}][\pmatrix{0&-1&1&0\cr 0&-1&0&1\cr 0&-1&0&0\cr 1&-1&0&0\cr}][\pmatrix{1&0&-1&0&0\cr 0&0&-1&1&0\cr 0&0&-1&0&1\cr 0&0&-1&0&0\cr 0&1&-1&0&0\cr}]
16 5 [\pmatrix{ -f/2&1/2&-t/2\cr -1/2&-t/2&f/2\cr t/2&f/2&1/2\cr}][\pmatrix{0&0&-1&1\cr 1&0&-1&0\cr 0&1&-1&0\cr 0&0&-1&0\cr} ][\pmatrix{0&-1&0&0&0\cr 0&-1&0&1&0\cr 0&-1&1&0&0\cr 0&-1&0&0&1\cr 1&-1&0&0&0\cr} ]
17 5 [\pmatrix{ -t/2&-f/2&1/2\cr -f/2&1/2&t/2\cr -1/2&-t/2&-f/2\cr}][\pmatrix{ 0 & -1 & 0 & 0\cr 0 & -1 & 1 & 0\cr 0 & -1 & 0 & 1\cr 1 & -1 & 0 & 0\cr} ][\pmatrix{0&0&0&0&-1\cr 1&0&0&0&-1\cr 0&1&0&0&-1\cr 0&0&0&1&-1\cr 0&0&1&0&-1\cr}]
18 5 [\pmatrix{ -t/2&-f/2&-1/2\cr -f/2&1/2&-t/2\cr 1/2&t/2&-f/2\cr} ][\pmatrix{-1&0&0&1\cr -1&0&0&0\cr -1&1&0&0\cr -1&0&1&0\cr}][\pmatrix{-1&1&0&0&0\cr -1&0&1&0&0\cr -1&0&0&0&1\cr -1&0&0&1&0\cr -1&0&0&0&0\cr}]
19 5 [\pmatrix{ -f/2&-1/2&t/2\cr 1/2&-t/2&f/2\cr -t/2&f/2&1/2\cr}][\pmatrix{0&1&0&-1\cr 0&0&1&-1\cr 0&0&0&-1\cr 1&0&0&-1\cr}][\pmatrix{-1&0&0&0&1\cr -1&0&0&0&0\cr -1&0&1&0&0\cr -1&1&0&0&0\cr -1&0&0&1&0\cr}]
20 5 [\pmatrix{ -f/2&1/2&t/2\cr -1/2&-t/2&-f/2\cr -t/2&-f/2&1/2\cr}][\pmatrix{0&-1&1&0\cr 1&-1&0&0\cr 0&-1&0&1\cr 0&-1&0&0\cr}][\pmatrix{0&1&0&-1&0\cr 0&0&0&-1&0\cr 1&0&0&-1&0\cr 0&0&1&-1&0\cr 0&0&0&-1&1\cr}]
21 5 [\pmatrix{ 1/2&-t/2&f/2\cr t/2&-f/2&-1/2\cr f/2&1/2&-t/2\cr}][\pmatrix{-1&1&0&0\cr -1&0&0&1\cr -1&0&0&0\cr -1&0&1&0\cr}][\pmatrix{0&0&0&1&-1\cr 0&1&0&0&-1\cr 0&0&0&0&-1\cr 0&0&1&0&-1\cr 1&0&0&0&-1\cr}]
22 5 [\pmatrix{ 1/2&-t/2&-f/2\cr t/2&-f/2&1/2\cr -f/2&-1/2&-t/2\cr} ][\pmatrix{0&0&0&-1\cr 0&0&1&-1\cr 1&0&0&-1\cr 0&1&0&-1\cr}][\pmatrix{0&0&0&1&0\cr 0&0&1&0&0\cr 1&0&0&0&0\cr 0&0&0&0&1\cr 0&1&0&0&0\cr}]
235 [\pmatrix{ 1/2&t/2&f/2\cr -t/2&-f/2&1/2\cr f/2&-1/2&-t/2\cr} ][\pmatrix{0&0&-1&0\cr 1&0&-1&0\cr 0&0&-1&1\cr 0&1&-1&0\cr}][\pmatrix{0&0&-1&0&1\cr 0&1&-1&0&0\cr 0&0&-1&1&0\cr 1&0&-1&0&0\cr 0&0&-1&0&0\cr} ]
24 5 [\pmatrix{ 1/2&t/2&-f/2\cr -t/2&-f/2&-1/2\cr -f/2&1/2&-t/2\cr}][\pmatrix{-1&0&1&0\cr -1&0&0&1\cr -1&1&0&0\cr -1&0&0&0\cr}][\pmatrix{0&0&1&0&0\cr 0&0&0&0&1\cr 0&1&0&0&0\cr 1&0&0&0&0\cr 0&0&0&1&0\cr}]
25 5 [\pmatrix{ -f/2&-1/2&-t/2\cr 1/2&-t/2&-f/2\cr t/2&-f/2&1/2\cr}][\pmatrix{0&1&0&-1\cr 0&0&0&-1\cr 1&0&0&-1\cr 0&0&1&-1\cr}][\pmatrix{0&-1&1&0&0\cr 1&-1&0&0&0\cr 0&-1&0&1&0\cr 0&-1&0&0&0\cr 0&-1&0&0&1\cr}]
26 3 [\pmatrix{ -1/2&t/2&-f/2\cr -t/2&f/2&1/2\cr f/2&1/2&-t/2\cr}][\pmatrix{1&-1&0&0\cr 0&-1&1&0\cr 0&-1&0&0\cr 0&-1&0&1\cr}][\pmatrix{0&-1&0&1&0\cr 0&-1&0&0&1\cr 1&-1&0&0&0\cr 0&-1&1&0&0\cr 0&-1&0&0&0\cr} ]
27 3 [\pmatrix{ -1/2&-t/2&f/2\cr t/2&f/2&1/2\cr -f/2&1/2&-t/2\cr}][\pmatrix{1&0&-1&0\cr 0&0&-1&0\cr 0&1&-1&0\cr 0&0&-1&1\cr} ][\pmatrix{0&0&1&0&-1\cr 0&0&0&0&-1\cr 0&0&0&1&-1\cr 1&0&0&0&-1\cr 0&1&0&0&-1\cr}]
28 3 [\pmatrix{ -1/2&t/2&f/1\cr -t/2&f/2&-1/2\cr -f/2&-1/2&-t/2\cr} ][\pmatrix{0&0&0&1\cr 0&1&0&0\cr 1&0&0&0\cr 0&0&1&0\cr}][\pmatrix{0&0&-1&1&0\cr 0&0&-1&0&0\cr 0&1&-1&0&0\cr 0&0&-1&0&1\cr 1&0&-1&0&0\cr}]
29 3 [\pmatrix{ 0&0&1\cr 1&0&0\cr 0&1&0\cr}][\pmatrix{0&0&0&1\cr 1&0&0&0\cr 0&0&1&0\cr 0&1&0&0\cr}][\pmatrix{0&1&0&0&-1\cr 0&0&1&0&-1\cr 1&0&0&0&-1\cr 0&0&0&0&-1\cr 0&0&0&1&-1\cr}]
30 3 [\pmatrix{ -1/2&-t/2&-f/2\cr t/2&f/2&-1/2\cr f/2&-1/2&-t/2\cr} ][\pmatrix{0&0&1&0\cr 0&1&0&0\cr 0&0&0&1\cr 1&0&0&0\cr}][\pmatrix{0&-1&0&0&1\cr 0&-1&1&0&0\cr 0&-1&0&0&0\cr 1&-1&0&0&0\cr 0&-1&0&1&0\cr}]
31 3 [\pmatrix{ 0&0&-1\cr 1&0&0\cr 0&-1&0\cr} ][\pmatrix{1&-1&0&0\cr 0&-1&0&1\cr 0&-1&1&0\cr 0&-1&0&0\cr} ][\pmatrix{0&0&0&0&-1\cr 0&0&1&0&-1\cr 0&0&0&1&-1\cr 0&1&0&0&-1\cr 1&0&0&0&-1\cr}]
32 3 [\pmatrix{ f/2&1/2&-t/2\cr 1/2&-t/2&f/2\cr t/2&-f/2&-1/2\cr}][\pmatrix{-1&0&1&0\cr -1&1&0&0\cr -1&0&0&0\cr -1&0&0&1\cr}][\pmatrix{-1&1&0&0&0\cr -1&0&0&0&0\cr -1&0&0&1&0\cr -1&0&0&0&1\cr -1&0&1&0&0\cr}]
33 3 [\pmatrix{ 0&1&0\cr 0&0&1\cr 1&0&0\cr}][\pmatrix{0&1&0&0\cr 0&0&0&1\cr 0&0&1&0\cr 1&0&0&0\cr}][\pmatrix{0&0&1&-1&0\cr 1&0&0&-1&0\cr 0&1&0&-1&0\cr 0&0&0&-1&1\cr 0&0&0&-1&0\cr}]
34 3 [\pmatrix{ 0&0&-1\cr -1&0&0\cr 0&1&0\cr} ][\pmatrix{0&0&0&-1\cr 0&1&0&-1\cr 0&0&1&-1\cr 1&0&0&-1\cr}][\pmatrix{0&1&-1&0&0\cr 0&0&-1&0&1\cr 0&0&-1&1&0\cr 0&0&-1&0&0\cr 1&0&-1&0&0\cr} ]
35 3 [\pmatrix{ 0&-1&0\cr 0&0&1\cr -1&0&0\cr}][\pmatrix{-1&0&0&1\cr -1&1&0&0\cr -1&0&1&0\cr -1&0&0&0\cr} ][\pmatrix{0&0&0&-1&1\cr 1&0&0&-1&0\cr 0&0&0&-1&0\cr 0&0&1&-1&0\cr 0&1&0&-1&0\cr}]
36 3 [\pmatrix{ f/2&-1/2&t/2\cr -1/2&-t/2&f/2\cr -t/2&-f/2&-1/2\cr}][\pmatrix{1&0&0&0\cr 0&0&0&1\cr 0&1&0&0\cr 0&0&1&0\cr}][\pmatrix{0&-1&0&0&1\cr 0&-1&0&1&0\cr 1&-1&0&0&0\cr 0&-1&0&0&0\cr 0&-1&1&0&0\cr}]
37 3 [\pmatrix{ 0&0&1\cr -1&0&0\cr 0&-1&0\cr}][\pmatrix{-1&1&0&0\cr -1&0&0&0\cr -1&0&1&0\cr -1&0&0&1\cr}][\pmatrix{0&0&-1&0&0\cr 0&0&-1&0&1\cr 1&0&-1&0&0\cr 0&1&-1&0&0\cr 0&0&-1&1&0\cr}]
38 3 [\pmatrix{ 0&1&0\cr 0&0&-1\cr -1&0&0\cr} ][\pmatrix{1&0&0&-1\cr 0&0&0&-1\cr 0&0&1&-1\cr 0&1&0&-1\cr}][\pmatrix{-1&0&0&0&1\cr -1&0&0&1&0\cr -1&1&0&0&0\cr -1&0&1&0&0\cr -1&0&0&0&0\cr}]
39 3 [\pmatrix{ f/2&1/2&t/2\cr 1/2&-t/2&-f/2\cr -t/2&f/2&-1/2\cr}][\pmatrix{0&0&-1&0\cr 0&1&-1&0\cr 1&0&-1&0\cr 0&0&-1&1\cr}][\pmatrix{0&-1&0&0&0\cr 1&-1&0&0&0\cr 0&-1&0&0&1\cr 0&-1&1&0&0\cr 0&-1&0&1&0\cr}]
40 3 [\pmatrix{ -t/2&-f/2&1/2\cr f/2&-1/2&-t/2\cr 1/2&t/2&f/2\cr}][\pmatrix{1&0&-1&0\cr 0&1&-1&0\cr 0&0&-1&1\cr 0&0&-1&0\cr}][\pmatrix{-1&0&0&1&0\cr -1&0&0&0&1\cr -1&1&0&0&0\cr -1&0&0&0&0\cr -1&0&1&0&0\cr} ]
41 3 [\pmatrix{ -t/2&-f/2&-1/2\cr f/2&-1/2&t/2\cr -1/2&-t/2&f/2\cr} ][\pmatrix{0&0&1&0\cr 1&0&0&0\cr 0&1&0&0\cr 0&0&0&1\cr}][\pmatrix{0&0&-1&1&0\cr 1&0&-1&0&0\cr 0&0&-1&0&1\cr 0&1&-1&0&0\cr 0&0&-1&0&0\cr}]
42 3 [\pmatrix{ -t/2&f/2&1/2\cr -f/2&-1/2&t/2\cr 1/2&-t/2&f/2\cr}][\pmatrix{1&0&0&-1\cr 0&1&0&-1\cr 0&0&0&-1\cr 0&0&1&-1\cr}][\pmatrix{0&0&0&-1&0\cr 0&0&1&-1&0\cr 0&0&0&-1&1\cr 1&0&0&-1&0\cr 0&1&0&-1&0\cr}]
43 3 [\pmatrix{ -t/2&f/2&-1/2\cr -f/2&-1/2&-t/2\cr -1/2&t/2&f/2\cr}][\pmatrix{0&1&0&0\cr 0&0&1&0\cr 1&0&0&0\cr 0&0&0&1\cr}][\pmatrix{0&1&0&0&-1\cr 0&0&0&1&-1\cr 0&0&0&0&-1\cr 1&0&0&0&-1\cr 0&0&1&0&-1\cr} ]
44 3 [\pmatrix{ f/2&-1/2&-t/2\cr -1/2&-t/2&-f/2\cr t/2&f/2&-1/2\cr}][\pmatrix{1&0&0&0\cr 0&0&1&0\cr 0&0&0&1\cr 0&1&0&0\cr}][\pmatrix{0&0&1&-1&0\cr 0&0&0&-1&0\cr 0&0&0&-1&1\cr 0&1&0&-1&0\cr 1&0&0&-1&0\cr}]
453 [\pmatrix{ 0&-1&0\cr 0&0&-1\cr 1&0&0\cr}][\pmatrix{0&-1&0&0\cr 1&-1&0&0\cr 0&-1&1&0\cr 0&-1&0&1\cr}][\pmatrix{-1&0&1&0&0\cr -1&0&0&1&0\cr -1&0&0&0&0\cr -1&0&0&0&1\cr -1&1&0&0&0\cr}]
46 2 [\pmatrix{ -1&0&0\cr 0&1&0\cr 0&0&-1\cr}][\pmatrix{0&1&0&-1\cr 1&0&0&-1\cr 0&0&1&-1\cr 0&0&0&-1\cr}][\pmatrix{0&0&0&1&0\cr 0&1&0&0&0\cr 0&0&0&0&1\cr 1&0&0&0&0\cr 0&0&1&0&0\cr}]
47 2 [\pmatrix{ -f/2&1/2&t/2\cr 1/2&t/2&f/2\cr t/2&f/2&-1/2\cr}][\pmatrix{-1&0&0&0\cr -1&1&0&0\cr -1&0&0&1\cr -1&0&1&0\cr}][\pmatrix{0&0&1&0&0\cr 0&0&0&1&0\cr 1&0&0&0&0\cr 0&1&0&0&0\cr 0&0&0&0&1\cr}]
48 2 [\pmatrix{ -f/2&-1/2&-t/2\cr -1/2&t/2&f/2\cr -t/2&f/2&-1/2\cr}][\pmatrix{0&0&1&0\cr 0&0&0&1\cr 1&0&0&0\cr 0&1&0&0\cr}][\pmatrix{-1&0&0&0&0\cr -1&1&0&0&0\cr -1&0&0&1&0\cr -1&0&1&0&0\cr -1&0&0&0&1\cr}]
49 2 [\pmatrix{ -f/2&-1/2&t/2\cr -1/2&t/2&-f/2\cr t/2&-f/2&-1/2\cr}][\pmatrix{1&0&-1&0\cr 0&0&-1&1\cr 0&0&-1&0\cr 0&1&-1&0\cr}][\pmatrix{0&1&0&0&0\cr 1&0&0&0&0\cr 0&0&1&0&0\cr 0&0&0&0&1\cr 0&0&0&1&0\cr} ]
50 2 [\pmatrix{ -f/2&1/2&-t/2\cr 1/2&t/2&-f/2\cr -t/2&-f/2&-1/2\cr}][\pmatrix{-1&0&0&0\cr -1&0&1&0\cr -1&1&0&0\cr -1&0&0&1\cr} ][\pmatrix{0&0&0&-1&1\cr 0&1&0&-1&0\cr 0&0&1&-1&0\cr 0&0&0&-1&0\cr 1&0&0&-1&0\cr}]
51 2 [\pmatrix{ t/2&f/2&-1/2\cr f/2&-1/2&-t/2\cr -1/2&-t/2&-f/2\cr}][\pmatrix{0&1&0&0\cr 1&0&0&0\cr 0&0&0&1\cr 0&0&1&0\cr} ][\pmatrix{-1&0&0&0&0\cr -1&0&0&0&1\cr -1&0&1&0&0\cr -1&0&0&1&0\cr -1&1&0&0&0\cr}]
52 2 [\pmatrix{ t/2&f/2&1/2\cr f/2&-1/2&t/2\cr 1/2&t/2&-f/2\cr} ][\pmatrix{1&0&0&-1\cr 0&0&1&-1\cr 0&1&0&-1\cr 0&0&0&-1\cr}][\pmatrix{0&1&-1&0&0\cr 1&0&-1&0&0\cr 0&0&-1&0&0\cr 0&0&-1&1&0\cr 0&0&-1&0&1\cr}]
53 2 [\pmatrix{ -1/2&t/2&f/2\cr t/2&-f/2&1/2\cr f/2&1/2&t/2\cr}][\pmatrix{0&-1&1&0\cr 0&-1&0&0\cr 1&-1&0&0\cr 0&-1&0&1\cr}][\pmatrix{0&0&0&0&1\cr 0&0&1&0&0\cr 0&1&0&0&0\cr 0&0&0&1&0\cr 1&0&0&0&0\cr}]
54 2 [\pmatrix{ -1/2&t/2&-f/2\cr t/2&-f/2&-1/2\cr -f/2&-1/2&t/2\cr} ][\pmatrix{0&0&-1&1\cr 0&1&-1&0\cr 0&0&-1&0\cr 1&0&-1&0\cr}][\pmatrix{0&0&1&0&-1\cr 0&1&0&0&-1\cr 1&0&0&0&-1\cr 0&0&0&1&-1\cr 0&0&0&0&-1\cr} ]
55 2 [\pmatrix{ 1&0&0\cr 0&-1&0\cr 0&0&-1\cr}][\pmatrix{0&-1&0&1\cr 0&-1&0&0\cr 0&-1&1&0\cr 1&-1&0&0\cr}][\pmatrix{0&-1&0&1&0\cr 0&-1&0&0&0\cr 0&-1&1&0&0\cr 1&-1&0&0&0\cr 0&-1&0&0&1\cr}]
56 2 [\pmatrix{ -1&0&0\cr 0&-1&0\cr 0&0&1\cr}][\pmatrix{-1&0&0&0\cr -1&0&0&1\cr -1&0&1&0\cr -1&1&0&0\cr} ][\pmatrix{1&-1&0&0&0\cr 0&-1&0&0&0\cr 0&-1&0&0&1\cr 0&-1&0&1&0\cr 0&-1&1&0&0\cr}]
57 2 [\pmatrix{ -1/2&-t/2&-f/2\cr -t/2&-f/2&1/2\cr -f/2&1/2&t/2\cr} ][\pmatrix{1&-1&0&0\cr 0&-1&0&0\cr 0&-1&0&1\cr 0&-1&1&0\cr} ][\pmatrix{1&0&-1&0&0\cr 0&1&-1&0&0\cr 0&0&-1&0&0\cr 0&0&-1&0&1\cr 0&0&-1&1&0\cr}]
58 2 [\pmatrix{ t/2&-f/2&-1/2\cr -f/2&-1/2&t/2\cr -1/2&t/2&-f/2\cr}][\pmatrix{0&1&-1&0\cr 1&0&-1&0\cr 0&0&-1&0\cr 0&0&-1&1\cr}][\pmatrix{1&0&0&-1&0\cr 0&0&1&-1&0\cr 0&1&0&-1&0\cr 0&0&0&-1&0\cr 0&0&0&-1&1\cr}]
59 2 [\pmatrix{ t/2&-f/2&1/2\cr -f/2&-1/2&-t/2\cr 1/2&-t/2&-f/2\cr}][\pmatrix{0&0&1&-1\cr 0&1&0&-1\cr 1&0&0&-1\cr 0&0&0&-1\cr}][\pmatrix{1&0&0&0&-1\cr 0&0&0&1&-1\cr 0&0&1&0&-1\cr 0&1&0&0&-1\cr 0&0&0&0&-1\cr} ]
60 2 [\pmatrix{ -1/2&-t/2&f/2\cr -t/2&-f/2&-1/2\cr f/2&-1/2&t/2\cr}][\pmatrix{0&0&0&1\cr 0&0&1&0\cr 0&1&0&0\cr 1&0&0&0\cr}][\pmatrix{1&0&0&0&0\cr 0&0&0&0&1\cr 0&0&0&1&0\cr 0&0&1&0&0\cr 0&1&0&0&0\cr} ]

Table 1.10.5.3| top | pdf |
The representation matrices for [\Gamma_3]

The representation matrices for [\Gamma_3] are the same as for [\Gamma_2]. Correspondences are given as pairs i, j: [\Gamma_3(R_{i}) = \Gamma_2(R_{j})].

ijijijijijij
1 1 11 21 21 531 42 41 29 51 48
2 1412 16 22 6 32 45 42 39 52 54
3 23 13 17 23 8 33 36 43 33 53 46
4 15 14 4 24 10 34 27 44 30 54 50
5 25 15 2 25 11 35 26 45 38 55 52
6 24 16 13 26 34 36 28 46 49 56 57
7 19 17 12 27 35 37 31 47 53 57 59
8 20 18 7 28 43 38 40 48 51 58 56
9 18 19 9 29 44 39 37 49 47 59 58
10 22 20 3 30 41 40 32 50 55 60 60

References

First citation Bak, P. (1985). Symmetry, stability, and elastic properties of icosahedral incommensurate crystals. Phys. Rev. B, 32, 5764–5772.Google Scholar
First citation Beest, B. W. van, Janner, A. & Blinc, R. (1983). 87Rb electric field gradient tensors and the symmetry of the incommensurate phase in Rb2ZnBr4 and Rb2ZnCl4. J. Phys. C, 16, 5409–5416.Google Scholar
First citation Ding, D.-H., Yang, W.-G., Hu, C.-Z. & Wang, R.-H. (1993). Generalized elasticity theory of quasicrystals. Phys. Rev. B, 48, 7003–7010.Google Scholar
First citation International Tables for Crystallography (2004). Vol. C. Mathematical, physical and chemical tables, edited by E. Prince. Dordrecht: Kluwer Academic Publishers.Google Scholar
First citation Janssen, T. (1988). Aperiodic crystals: a contradictio in terminis? Phys. Rep. 168, 57–113.Google Scholar
First citation Janssen, T., Birman, J. L., Dénoyer, F., Koptsik, V. A., Verger-Gaugry, J. L., Weigel, D., Yamamoto, A., Abrahams, S. C. & Kopsky, V. (2002). Report of a subcommittee on the nomenclature of n-dimensional crystallography. II. Symbols for arithmetic crystal classes, Bravais classes and space groups. Acta Cryst. A58, 605–621.Google Scholar
First citation Janssen, T., Birman, J. L., Koptsik, V. A., Senechal, M., Weigel, D., Yamamoto, A., Abrahams, S. C. & Hahn, Th. (1999). Report of a subcommittee on the nomenclature of n-dimensional crystallography. I. Symbols for point-group transformations, families, systems and geometric crystal classes. Acta Cryst. A55, 761–782.Google Scholar
First citation Janssen, T. & Janner, A. (1987). Incommensurability in crystals. Adv. Phys. 36, 519–624.Google Scholar
First citation Lubensky, T. C., Ramaswamy, S. & Toner, J. (1985). Hydrodynamics of icosahedral quasicrystals. Phys. Rev. B, 32, 7444–7452.Google Scholar
First citation 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.Google Scholar
First citation Socolar, J. E. S., Lubensky, T. C. & Steinhardt, P. J. (1986). Phonons, phasons and dislocations in quasicrystals. Phys. Rev. B, 34, 3345–3360.Google Scholar








































to end of page
to top of page