research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoFOUNDATIONS
ADVANCES
ISSN: 2053-2733

Layer groups: Brillouin-zone and crystallographic databases on the Bilbao Crystallographic Server

crossmark logo

aInstitute of Applied Geosciences, Karlsruhe Institute of Technology, Karlsruhe, Germany, bInstitute for Mathematics, Astrophysics and Particle Physics, Radboud University Nijmegen, The Netherlands, and cDepartamento de Física, Universidad del País Vasco UPV/EHU, Spain
*Correspondence e-mail: gemma.delaflor@kit.edu

Edited by U. Grimm, The Open University, United Kingdom (Received 7 July 2021; accepted 30 July 2021; online 24 September 2021)

The section of the Bilbao Crystallographic Server (https://www.cryst.ehu.es/) dedicated to subperiodic groups contains crystallographic and Brillouin-zone databases for the layer groups. The crystallographic databases include the generators/general positions (GENPOS), Wyckoff positions (WYCKPOS) and maximal subgroups (MAXSUB). The Brillouin-zone database (LKVEC) offers k-vector tables and Brillouin-zone figures of all 80 layer groups which form the background of the classification of their irreducible representations. The symmetry properties of the wavevectors are described applying the so-called reciprocal-space-group approach and this classification scheme is compared with that of Litvin & Wike [(1991), Character Tables and Compatibility Relations of the Eighty Layer Groups and Seventeen Plane Groups. New York: Plenum Press]. The specification of independent parameter ranges of k vectors in the representation domains of the Brillouin zones provides a solution to the problems of uniqueness and completeness of layer-group representations. The Brillouin-zone figures and k-vector tables are described in detail and illustrated by several examples.

1. Introduction

The Bilbao Crystallographic Server (https://www.cryst.ehu.es/) (Aroyo et al., 2006[Aroyo, M. I., Perez-Mato, J. M., Capillas, C., Kroumova, E., Ivantchev, S., Madariaga, G., Kirov, A. & Wondratschek, H. (2006). Z. Kristallogr. 221, 15-27.], 2011[Aroyo, M. I., Perez-Mato, J. M., Orobengoa, D., Tasci, E. S., de la Flor, G. & Kirov, A. (2011). Bulg. Chem. Commun. 43, 183-197.]; Tasci et al., 2012[Tasci, E. S., de la Flor, G., Orobengoa, D., Capillas, C., Perez-Mato, J. M. & Aroyo, M. I. (2012). EJP Web Conf. 22, 00009.]) is a free web site that grants access to specialized databases and tools for the resolution of different types of crystallographic, structural chemistry and solid-state physics problems. The server has been operating for more than 20 years, and is in constant improvement and development, offering free-of-charge tools to study an increasing number of crystallographic systems (Elcoro et al., 2017[Elcoro, L., Bradlyn, B., Wang, Z., Vergniory, M. G., Cano, J., Felser, C., Bernevig, B. A., Orobengoa, D., de la Flor, G. & Aroyo, M. I. (2017). J. Appl. Cryst. 50, 1457-1477.]; Gallego et al., 2019[Gallego, S. V., Etxebarria, J., Elcoro, L., Tasci, E. S. & Perez-Mato, J. M. (2019). Acta Cryst. A75, 438-447.]; de la Flor et al., 2019[Flor, G. de la, Orobengoa, D., Evarestov, R. A., Kitaev, Y. E., Tasci, E. & Aroyo, M. I. (2019). J. Appl. Cryst. 52, 1214-1221.]). The programs on the server are organized into different sections depending on their degree of complexity, in such a way that the more complex tools make use of the results obtained by the simpler ones. The Bilbao Crystallographic Server (hereafter referred to as BCS) is built on a core of databases that include data from International Tables for Crystallography, Vol. A, Space-group Symmetry (Aroyo, 2016[Aroyo, M. I. (2016). Editor. International Tables for Crystallography, Vol. A, Space-group Symmetry, 6th ed. Chichester: Wiley. ]; henceforth abbreviated as ITA), Vol. A1, Symmetry Relations between Space Groups (Wondratschek & Müller, 2010[Wondratschek, H. & Müller, U. (2010). Editors. International Tables for Crystallography, Vol. A1, Symmetry Relations between Space Groups, 2nd ed. Chichester: John Wiley & Sons.]), and Vol. E, Subperiodic Groups (Kopský & Litvin, 2010[Kopský, V. & Litvin, D. (2010). Editors. International Tables for Crystallography, Vol. E, Subperiodic Groups, 2nd ed. Chichester: Wiley. ]; henceforth referred to as ITE). A k-vector database with Brillouin-zone figures and classification tables of the wavevectors for all 230 space groups (Aroyo et al., 2014[Aroyo, M. I., Orobengoa, D., de la Flor, G., Tasci, E. S., Perez-Mato, J. M. & Wondratschek, H. (2014). Acta Cryst. A70, 126-137.]) is also available in the server, together with the magnetic and double-space-group databases. The magnetic and the incommensurate structure databases are also hosted on the server.

Aside from the subgroups of space groups with three-dimensional lattices which are again space groups, there also exist subgroups called subperiodic groups with translation lattices of dimensions one or two. These groups are essential to describe polymers, nanotubes, nanowires, layered and multilayered materials. The interest in materials with sub­periodic symmetry is constantly growing due to their outstanding properties and possible technological applications. There are three types of subperiodic groups: frieze groups (two-dimensional groups with one-dimensional translation lattices), rod groups (three-dimensional groups with one-dimensional translation lattices) and layer groups (three-dimensional groups with two-dimensional translation lattices). The crystallographic data for subperiodic groups are compiled in ITE and also now offered online in the BCS section dedicated to subperiodic groups. This section includes programs which give access to the generators/general positions (GENPOS), Wyckoff positions (WYCKPOS) and maximal subgroups (MAXSUB) databases.

In addition, we have developed the Brillouin-zone database of layer groups which contains k-vector tables and Brillouin-zone figures that form the background of the classification of the irreducible representations of layer groups. The Brillouin-zone figures and the wavevector data for space groups are well established and many authors have contributed to the standardization of the data (see e.g. Miller & Love, 1967[Miller, S. C. & Love, W. F. (1967). Tables of Irreducible Representations of Space Groups and Co-Representations of Magnetic Space Groups. Boulder: Pruett Press.]; Bradley & Cracknell, 1972[Bradley, C. J. & Cracknell, A. P. (1972). The Mathematical Theory of Symmetry in Solids: Representation Theory for Point Groups and Space Groups. Oxford: Clarendon Press.]; Cracknell et al., 1979[Cracknell, A. P., Davies, B. L., Miller, S. C. & Love, W. F. (1979). Kronecker Product Tables, Vol. 1, General Introduction and Tables of Irreducible Representations of Space Groups. New York: IFI/Plenum.]; Stokes & Hatch, 1988[Stokes, H. T. & Hatch, D. M. (1988). Isotropy Subgroups of the 230 Crystallographic Space Groups. Singapore: World Scientific.]). For layer groups, however, the Brillouin-zone and wavevector descriptions proposed by several authors (Ipatova & Kitaev, 1985[Ipatova, I. P. & Kitaev, Y. E. (1985). Prog. Surf. Sci. 18, 189-246.]; Hatch & Stokes, 1986[Hatch, D. M. & Stokes, H. T. (1986). Phase Transit. 7, 87-279.]; Milošević et al., 2008[Milošević, I., Nikolić, B., Damnjanović, M. & Krčmar, M. (2008). J. Phys. A, 31, 3625-3648.]) are incomplete and difficult to compare due to the lack of standards in the classification and nomenclature of the k vectors. To the best of our knowledge, the only complete compilation of layer-group k vectors together with the Brillouin-zone diagrams is found in the work of Litvin & Wike (1991[Litvin, D. B. & Wike, T. R. (1991). Character Tables and Compatibility Relations of the Eighty Layer Groups and Seventeen Plane Groups. New York: Plenum Press.]) (in the following, referred to as LW). In this description the k vectors are labelled following the classification scheme of space-group k vectors used by Cracknell et al. (1979[Cracknell, A. P., Davies, B. L., Miller, S. C. & Love, W. F. (1979). Kronecker Product Tables, Vol. 1, General Introduction and Tables of Irreducible Representations of Space Groups. New York: IFI/Plenum.]). Based on the group–subgroup relations between space and layer groups and using the so-called reciprocal-space-group approach (cf. Wintgen, 1941[Wintgen, G. (1941). Math. Ann. 118, 195-215.]; Aroyo & Wondratschek, 1995[Aroyo, M. I. & Wondratschek, H. (1995). Z. Kristallogr. 210, 243-254.]; Aroyo et al., 2014[Aroyo, M. I., Orobengoa, D., de la Flor, G., Tasci, E. S., Perez-Mato, J. M. & Wondratschek, H. (2014). Acta Cryst. A70, 126-137.]), we have derived the k-vector data and Brillouin-zone figures for all 80 layer groups and compared them with the classification of LW. An auxiliary tool for the complete characterization of wavevectors is also available: given the wavevector coefficients referred to a primitive or conventional dual basis, the program assigns the corresponding wavevector symmetry type, specifies its LW label, determines the layer little co-group of the wavevector and generates the arms of the wavevector star. The BCS Brillouin-zone database of layer groups is accessed by the retrieval tool LKVEC.

The aim of this contribution is to present the crystallographic databases and the Brillouin-zone database for layer groups available in the BCS. In the following, the retrieval tools GENPOS, WYCKPOS, MAXSUB and LKVEC, and the procedure to derive and classify the k vectors of layer groups are described in detail. The utility of the programs is demonstrated by several illustrative examples.

2. Crystallographic databases for layer groups

The BCS section Subperiodic Groups: Layer, Rod and Frieze Groups hosts the layer-group crystallographic databases. The structure of these databases is similar to that of the space groups – they include information on generators, general positions, Wyckoff positions and maximal subgroups for the 80 layer groups. Apart from the data shown in ITE the server offers additional information and computer tools that allow the generation of data not available in ITE. The BCS programs and databases use the so-called standard or default settings of the layer groups. These are the specific settings of layer groups that coincide with the conventional layer-group descriptions found in ITE. For layer groups with more than one description in ITE, the following settings are chosen as standard: (i) cell-choice 1 description for the two monoclinic/oblique layer groups p11a (No. 5) and p112/a (No. 7) described with respect to three cell choices in ITE, and (ii) origin choice 2 descriptions (i.e. origin at an inversion centre) for the three layer groups p4/n (No. 52), p4/nbm (No. 62) and p4/nmm (No. 64) listed with respect to two origins in ITE.

Note that, following the conventions of ITE, the ab plane is the plane of periodicity for the layer groups and therefore the translation vectors are of the form

[\left(\matrix{t_{1}\cr t_{2}\cr 0}\right).]

2.1. Generators and general positions

The BCS database of layer groups includes the list of generators/general positions of each layer group. These data can be retrieved using the program GENPOS, by specifying the sequential ITE layer-group number (which, if unknown, can be determined by choosing the corresponding Hermann–Mauguin symbol from a table with the layer-group symbols). The generators and/or general positions of layer groups are specified by their coordinate triplets, the matrix-column representations of the corresponding symmetry operations and their geometric interpretation:

(a) The list of coordinate triplets (x, y, z) reproduces the data of the General Positions blocks of layer groups found in ITE. The coordinate triplets may also be interpreted as shorthand descriptions of the matrix presentation of the corresponding symmetry operations.

(b) Matrix-column presentation of symmetry operations. With reference to a coordinate system, consisting of an origin O and a basis (a1, a2, a3), the symmetry operations of layer groups are described by (3 × 4) matrix-column pairs.

(c) Geometric interpretation. The geometric interpretation of symmetry operations is given (i) following the conventions in ITE [including the symbol of the symmetry operation, its glide or screw component (if relevant) and the location of the related symmetry element], and (ii) using the Seitz notation. It is worth pointing out that, in contrast to the online and printed editions of ITE, the programs of BCS use the standard International Union of Crystallography Seitz notation. For example, a twofold rotation around the c axis is denoted by 2001 instead of 2z (for details, cf. Litvin & Kopský, 2014[Litvin, D. B. & Kopský, V. (2014). Acta Cryst. A70, 677-678.]).

Fig. 1[link] shows the general position for the layer group p112/a (No. 7) in the default setting (cell-choice 1).

[Figure 1]
Figure 1
General-position table of the layer group p112/a (No. 7), provided by the program GENPOS.

The program GENPOS lists the generators and/or general positions of layer groups in the standard/default setting as well as in conventional non-standard settings of monoclinic/rectangular and orthorhombic/rectangular layer groups described in Table 1.2.6.1 of ITE (option ITE settings). In addition, the program can produce the data in any non-conventional setting if the transformation relating the non-conventional setting to the standard one is specified (option Non-conventional setting). The matrix-column pair (P, p) of the transformation relating the two settings consists of two parts: a linear part P defined by a (3 × 3) matrix, which describes the change of direction and/or length of the basis vectors (a1, a2, a3)non-conv = (a1, a2, a3)stanP, and an origin shift p = (p1,p2,p3) defined as a (3 × 1) column, whose coefficients describe the position of the non-conventional origin with respect to the standard one.

The URL of the program GENPOS is https://www.cryst.ehu.es/subperiodic/get_sub_gen.html.

2.2. Wyckoff positions

The BCS Wyckoff-positions database for layer groups is accessible under the WYCKPOS program. The data on Wyckoff positions can be retrieved by specifying the ITE number of the layer-group type. As a result, the program WYCKPOS shows a table with the Wyckoff positions [see Fig. 2[link](a)]. Following ITE, each Wyckoff position is characterized by its multiplicity, Wyckoff letter, site-symmetry group and a set of coordinate triplets of the Wyckoff-position points in the unit cell. For centred subperiodic groups, the centring translations are listed above the coordinate triplets. The site-symmetry groups [see column three of Fig. 2[link](a)] are described by oriented symbols which display the same sequence of symmetry directions as the layer-group symbols (cf. Table 1.2.4.1 of ITE). An explicit listing of the symmetry operations of the site-symmetry group of a point is obtained by clicking directly on its coordinate triplet. A recently implemented auxiliary tool permits the identification of the Wyckoff position and the site-symmetry operations of a point specified by its coordinate triplet [see Fig. 2[link](b)]. The program accepts as input relative point coordinates in fractions, decimals or variable parameters (indicating a generic value). Fig. 3[link] shows the list of the symmetry operations of the site-symmetry group of the layer group cmmm (No. 47) for the points (5/4,3/4,z). The points belong to the Wyckoff position 8h with site-symmetry group ..2.

[Figure 2]
Figure 2
A screenshot of the output of the program WYCKPOS, showing (a) the Wyckoff positions of the layer group cmmm (No. 47), and (b) the interface of the auxiliary tool for the determination of the Wyckoff position and the site-symmetry group of a specific point.
[Figure 3]
Figure 3
Output of the auxiliary tool in WYCKPOS that shows the symmetry operations of the site-symmetry group for points with coordinates (5/4,3/4,z) of the layer group cmmm (No. 47). The points belong to the Wyckoff position 8h with site-symmetry group ..2.

Apart from the standard/default setting option the program is also able to calculate the Wyckoff positions in different ITE (conventional) settings (option ITE settings) or with respect to a non-conventional setting if the corresponding coordinate transformation (P, p) is defined (option Non-conventional setting).

The URL of the program WYCKPOS is https://www.cryst.ehu.es/subperiodic/get_sub_wp.html.

2.3. Maximal subgroups

The listing of maximal subgroups of layer groups available in ITE is incomplete and lacks additional information, such as, for example, possible unit-cell transformations and/or origin shifts involved. In contrast, the BCS database of maximal subgroups of layer groups provides the complete listing (not just by type but individually) of (i) all maximal non-isotypic subgroups for each layer group, and (ii) all maximal isotypic subgroups of indices 2, 3 and 4. The list of maximal subgroups is retrieved by the access tool MAXSUB. Each subgroup is specified by its ITE number, Hermann–Mauguin symbol, index, subgroup type (t for translationengleiche or k for klassengleiche) and transformation matrix-column pair (P, p) that relates the standard setting of the group with that of the subgroup (see Fig. 4[link]). The different maximal subgroups are distributed into conjugacy classes. The identification of the subgroup symmetry operations as a subset of the elements of the group is achieved by an optional tool of MAXSUB that transforms the general-position representatives of the subgroup to the coordinate system of the group.

[Figure 4]
Figure 4
Maximal subgroups of the layer group cmmm (No. 47) provided by the program MAXSUB.

The URL of the program MAXSUB is https://www.cryst.ehu.es/subperiodic/get_sub_maxsub.html.

3. Classification of wavevectors

A layer group [{\cal L}] is defined as a three-dimensional crystallographic group with periodicity restricted to a two-dimensional subspace. Just as for space groups, the general strategy for determining the irreducible representations, called irreps for short, of such a group is to exploit the fact that it contains the translation subgroup [{\cal T}] as a normal abelian subgroup and that irreps of abelian groups are one-dimensional. The elements of [{\cal T}] have matrix-column pairs of the form [({\bf I},{\bf t})] where [{\bf t}] lies in a two-dimensional lattice [{\bf L}] with basis [{\bf a}_{1},{\bf a}_{2}]. Each one-dimensional irrep of [{\cal T}] is then of the form [\Gamma^{{\bf k}}({\bf I},{\bf t}) = \exp(-2\pi i\,{\bf k}\cdot {\bf t})] for a [{\bf k}] vector [{\bf k} = k_{1}{\bf a}_{1}^{*}+k_{2}{\bf a}_{2}^{*}] in the reciprocal space spanned by the reciprocal basis [{\bf a}_{1}^{*},{\bf a}_{2}^{*}] defined by [{\bf a}_{i}\cdot{\bf a}_{j}^{*} = \delta_{ij}].

Starting from an irrep [\Gamma^{{\bf k}}] of [{\cal T}], an irrep of the full layer group [{\cal L}] is obtained by first extending [\Gamma^{{\bf k}}] to the stabilizer of [\Gamma^{{\bf k}}] under the conjugation action of [{\cal L}] and then inducing to the full group [{\cal L}]. When the induced irrep of [{\cal L}] is restricted to [{\cal T}], it is the sum of those one-dimensional irreps [\Gamma^{{\bf k}^{\prime}}] of [{\cal T}] which lie in the orbit of [\Gamma^{{\bf k}}] under the conjugation action of [{\cal L}] and these are precisely the irreps of [{\cal T}] that yield irreps of [{\cal L}] equivalent to the one obtained from [\Gamma^{{\bf k}}]. For the action of [{\cal L}] on the [{\bf k}] vectors only the part acting on the plane of periodicity of [{\cal L}] is relevant and the restriction of [{\cal L}] to this plane gives rise to a plane group [{\cal P}] associated to [{\cal L}]. In the convention of ITE, the matrix-column pairs of [{\cal L}] are expressed with respect to a basis [{\bf a},{\bf b},{\bf c}] such that the first two basis vectors span the plane of periodicity. Therefore, the matrices of the linear parts of the plane group [{\cal P}] are simply the upper 2×2 diagonal blocks of the linear parts of the layer group [{\cal L}].

An explicit computation shows that the [{\bf k}^{\prime}] in the orbit of [{\bf k}] are of the form [{\bf k}^{\prime} = {\bf k}{\bf W}+{\bf K}] for [{\bf W}] an element of the point group [\overline{{\cal P}}] of [{\cal P}] and [{\bf K}] an element from the reciprocal lattice [{\bf L}^{*}] of [{\bf L}], i.e. for [{\bf K} = k_{1}{\bf a}_{1}^{*}+k_{2}{\bf a}_{2}^{*}] with k1,k2 integers. The set of elements [{\bf W}] for which [{\bf k} = {\bf k}{\bf W}+{\bf K}] is called the little co-group [\overline{{\cal P}}^{k}] of [{\bf k}]. The k vector is called general if [\overline{{\cal P}}^{k}] contains only the identity element of [\overline{{\cal P}}], i.e. [\overline{{\cal P}}^{k}] = {I}; otherwise [\overline{{\cal P}}^{k}] > {I} and k is special (Bradley & Cracknell, 1972[Bradley, C. J. & Cracknell, A. P. (1972). The Mathematical Theory of Symmetry in Solids: Representation Theory for Point Groups and Space Groups. Oxford: Clarendon Press.]; Dresselhaus et al., 2008[Dresselhaus, M. S., Dresselhaus, G. & Jorio, A. (2008). Group Theory Application to the Physics of Condensed Matter. Berlin: Springer.]). If {[{\bf W}_{m}]} is the set of coset representatives of the decomposition of [\overline{{\cal P}}] with respect to [\overline{{\cal P}}^{k}], then the set of [{\bf k}] vectors {[{\bf k}{\bf W}_{m}]} is called the star of [{\bf k}] while the vectors [{\bf k}{\bf W}_{m}] are called the arms of the star. Using the relation between the elements of the layer group [{\cal L}] and the plane group [{\cal P}] associated to [{\cal L}], in analogy to [\overline{{\cal P}}^{k}], one defines the little layer co-group [\overline{{\cal L}}^{k}] which is essential for the derivation of the irreps of [{\cal L}].

The preceding discussion indicates the importance of the so-called reciprocal plane group [({\cal P})^{*}], defined as the symmorphic plane group having the reciprocal lattice [{\bf L}^{*}] as translation vectors and [\overline{{\cal P}}] as point group. It is crucial that the operations in [({\cal P})^{*}] act on reciprocal space, i.e. they act on rows, as opposed to [{\cal P}] which acts on columns. In order to identify [({\cal P})^{*}] with one of the symmorphic plane groups in their standard setting, the action on reciprocal space (on rows) has to be transformed into one on direct space (on columns).

3.1. Transformation between reciprocal space and direct space

After fixing a basis [{\bf a}_{1},{\bf a}_{2}] for direct space, the action of the point group [\overline{{\cal P}}] of [{\cal P}] on [{\bf L}] is given by 2×2 matrices [{\bf W}] acting on columns. We now have to relate the action on the reciprocal lattice [{\bf L}^{*}] to that on [{\bf L}]. Defining

[{\bf b}_{i} = \left(\matrix{b_{i,1}\cr b_{i,2}}\right) = {\bf W}\left(\matrix{a_{i,1}\cr a_{i,2}}\right) = {\bf W}{\bf a}_{i},]

the basis [{\bf a}_{1},{\bf a}_{2}] is mapped by [{\bf W}] to the new basis [{\bf b}_{1},{\bf b}_{2}]. The corresponding matrix [{\bf W}^{*}] for the action on reciprocal space (on rows) must then map the reciprocal basis [{\bf a}_{1}^{*},{\bf a}_{2}^{*}] to the reciprocal basis [{\bf b}_{1}^{*},{\bf b}_{2}^{*}] of [{\bf b}_{1},{\bf b}_{2}]. From

[{\bf a}_{j}^{*}\cdot{\bf a}_{i} = \delta_{ij} = {\bf b}_{j}^{*}\cdot {\bf b}_{i} = ({\bf a}_{j}^{*}{\bf W}^{*})({\bf W}{\bf a}_{i}) = {\bf a}_{j}^{*}({\bf W}^{*}{\bf W}){\bf a}_{i}]

it follows that [{\bf W}^{*}{\bf W}] must be the identity matrix [{\bf I}]. We thus have [{\bf W}^{*} = {\bf W}^{-1}], i.e. the action on the rows of reciprocal space is given by the inverse matrices acting on the columns of direct space. Since [\overline{{\cal P}}] is a group, it contains with every matrix [{\bf W}] of course also the inverse matrix [{\bf W}^{-1}]; hence the set of matrices remains the same, but they now act from the right on rows.

As an example, we take the layer group [p\bar{6}m2] (No. 78). Restricting the action to the plane of periodicity gives the plane group p3m1 (No. 14). If we look at the [{\bf k}] vector [(-1.21,1)], applying the threefold rotation with matrix

[\left(\matrix{0&-1\cr 1&-1}\right)]

to the row [{\bf k}] gives [{\bf k}^{\prime} = (1,0.21)] and applying it again gives [{\bf k}^{\prime\prime} = (0.21,-1.21)]. On the other hand, applying the reflection with reflection line 2x,x and matrix

[\left(\matrix{1&0\cr 1&-1}\right),]

we see that the row [{\bf k}] is mapped to [(-0.21,-1)], which is equal to [{\bf k}] up to a reciprocal-lattice vector. The star of [{\bf k}] consisting of the [{\bf k}] vectors in the orbit under [{\cal L}] is therefore [{}^{*}{\bf k} = \{(-1.21,1),(1,0.21),(0.21,-1.21)\}].

As we have just demonstrated, it is natural to consider the action of the reciprocal plane group [({\cal P})^{*}] on rows, but on the other hand [({\cal P})^{*}] is by construction isomorphic to a symmorphic plane group [{\cal P}_{0}]. In order to identify [({\cal P})^{*}] with the proper symmorphic plane group, the action on rows has to be transformed to an action on columns. This is of course achieved by simply transposing the matrices of the point group [\overline{{\cal P}}], since [({\bf k}{\bf W})^{\rm tr} = {\bf W}^{\rm tr}{\bf k}^{\rm tr}] and the transpose [{\bf k}^{\rm tr}] of a row is a column. The symmorphic plane group [{\cal P}_{0}] isomorphic to [({\cal P})^{*}] is therefore the group

[{\cal P}_{0} = \{({\bf W}^{\rm tr},{\bf t}_{{\bf K}}^{\rm tr})\mid{\bf W}\in \overline{{\cal P}},{\bf t}_{{\bf K}}\in{\bf L}^{*}\}.]

The only problem that may occur is that [{\cal P}_{0}] is not given in its standard setting. This can be seen from different (strongly interrelated) perspectives:

(i) The reciprocal basis [{\bf a}_{1}^{*},{\bf a}_{2}^{*}] may not be a conventional basis of the lattice [{\bf L}^{*}]. For example, the reciprocal lattice of a hexagonal lattice is also a hexagonal lattice, but while the vectors [{\bf a}_{1},{\bf a}_{2}] in the conventional basis have an angle of 120°, the reciprocal basis vectors [{\bf a}_{1}^{*},{\bf a}_{2}^{*}] enclose an angle of 60° and are therefore not a conventional basis of a hexagonal lattice.

(ii) The matrices [{\bf W}] of the point group fix the metric tensor [{\bf G}] of the lattice [{\bf L}], i.e. one has [{\bf W}^{\rm tr}{\bf G}{\bf W} = {\bf G}]. The transposed matrices fix the inverse metric tensor [{\bf G}^{-1}] as can be seen from inverting the above equation: [{\bf W}^{\rm tr}{\bf G}{\bf W}] = [ {\bf G}\,\Longrightarrow\,{\bf W}^{-1} {\bf G}^{-1}{\bf W}^{-{\rm tr}}] = [{\bf G}^{-1}\,\Longrightarrow\,{\bf W} {\bf G}^{-1}{\bf W}^{\rm tr} = {\bf G}^{-1}]. This corresponds to the fact that the reciprocal lattice [{\bf L}^{*}] has metric tensor [{\bf G}^{-1}]. But [{\bf G}^{-1}] may not be the metric tensor of a lattice with respect to its conventional basis. For example, a hexagonal lattice has a metric tensor of the form

[{\bf G} = a^{2}\left(\matrix{2&-1\cr -1&2}\right)]

which has inverse

[{\bf G}^{-1} = {{1} \over {3a^{2}}}\left(\matrix{2&1\cr 1&2}\right).]

This is the metric tensor of a hexagonal lattice with basis vectors having an angle of 60°, i.e. in a non-conventional setting.

(iii) The transposed matrices [{\bf W}^{\rm tr}] simply do not occur as matrices of one of the point groups of plane groups in their standard setting. For example, the transposed matrix

[\left(\matrix{0&1\cr -1&-1}\right)]

of the matrix

[\left(\matrix{0&-1\cr 1&-1}\right)]

of a threefold rotation 3+ does not belong to the point group of any of the plane groups.

For plane groups, the hexagonal lattice is actually the only case in which the transposed matrices do not belong to a point group in its standard setting. For all other lattices, the inverse of the metric tensor still belongs to a conventional basis or, in other words, the set of matrices of the point group does not change by transposing.

In order to identify the symmorphic plane group [{\cal P}_{0}] in the case of a hexagonal lattice, the group has to be transformed to a basis in the conventional setting, i.e. which has a metric tensor equal to a multiple of

[\left(\matrix{2&-1\cr -1&2}\right).]

Such a transformation is

[{\bf P} = \left(\matrix{1&-1\cr 0&1}\right)]

(or its compositions with powers of a sixfold rotation) and the transposed matrices of the point group have to be conjugated by this. This interchanges the two types of reflections (in p6mm), those with normals along the two basis vectors and their sum and those with reflection lines along the basis vectors (and their sum). For example, for

[{\bf W} = \left(\matrix{1&0\cr 1&-1}\right),]

the reflection with reflection line 2x,x, one gets

[{\bf P}^{-1}{\bf W}^{\rm tr}{\bf P} = \left(\matrix{1&1\cr 0&1}\right) \left(\matrix{1&1\cr 0&-1}\right) \left(\matrix{1&-1\cr 0&1}\right) = \left(\matrix{1&-1\cr 0&-1}\right)]

which is the reflection with reflection line along the a axis.

As a result, for a layer group with plane group of type p3m1, the symmorphic plane group [{\cal P}_{0}] isomorphic to the reciprocal plane group [({\cal P})^{*}] is of type p31m and vice versa.

Based on the identification of the reciprocal plane group [({\cal P})^{*}] with a symmorphic plane group [{\cal P}_{0}], the special [{\bf k}] vectors for a layer group [{\cal L}] can be directly read off from the Wyckoff positions of [{\cal P}_{0}].

3.2. Crystallographic conventions in the classification of layer-group irreps

The isomorphism between a reciprocal plane group [({\cal P})^{*}] and a symmorphic plane group [{\cal P}_{0}] allows the application of crystallographic conventions in the classification of the wavevectors (and henceforth of the irreps) of the layer groups [{\cal L}]:

(i) The unit cells of the symmorphic plane groups listed in ITA can replace the Brillouin zones as unit cells of the reciprocal space. To find all irreps of [{\cal L}], it is necessary to consider only the wavevectors of the so-called representation domain. It is defined as a simply connected part of the (first) Brillouin zone (a unit cell of the reciprocal space) which contains exactly one k vector of each orbit of k. The asymmetric units of plane groups can serve as representation domains. The advantage of choosing the crystallographic unit cells and their asymmetric units becomes evident in layer groups where the Brillouin zones may belong to different topological types depending on the ratios of the lattice parameters. Lines on the Brillouin zone may appear or disappear or change their relative sizes depending on the lattice parameters. In contrast to that, the unit cells and their asymmetric units of ITA are independent of the ratios of the lattice parameters.

(ii) The action of the reciprocal plane group [({\cal P})^{*}] on the wavevectors results in their distribution into orbits of symmetry-equivalent k vectors with respect to [({\cal P})^{*}]. Thanks to the isomorphism of ([{\cal P}]) * with the symmorphic plane group [{\cal P}_{0}], the different types of k vectors correspond to the different kinds of point orbits (Wyckoff positions) of [{\cal P}_{0}]. In this way, a complete list of the special sites in the Brillouin zone of ([{\cal P}]) * is provided by the Wyckoff positions of [{\cal P}_{0}] found in ITA. The site-symmetry groups of ITA correspond to the little co-groups of the wavevectors and the number of arms of the star of a wavevector follows from the multiplicity of the Wyckoff position. The Wyckoff positions with zero, one and two variable parameters correspond to special k-vector points, special k-vector lines and special (or general) k-vector planes, respectively. A k-vector type, i.e. the set of all k vectors corresponding to a Wyckoff position, consists of complete orbits of k vectors and thus of full stars of k vectors. The different orbits (and stars) of a k-vector type are obtained by varying the free parameters. Correspondingly, the irreps of k vectors of a k-vector type are interrelated by parameter variation and are said to belong to the same type of irreps (Boyle, 1986[Boyle, L. L. (1986). Proceedings of the 14th International Colloquium on Group Theoretical Methods in Physics, pp. 405-408. Singapore: World Scientific.]). In this way all wavevector stars giving rise to the same type of irreps are related to the same Wyckoff position and designated by the same Wyckoff letter.

(iii) A complete set of irreps of [{\cal L}] is derived by considering exactly one k-vector representative per k-vector orbit. To achieve that, it is necessary to specify the exact parameter ranges of the independent k-vector regions within the representation domain (or the asymmetric unit). While such data are not available in the literature, the Brillouin-zone database of BCS offers the listing of the exact parameter ranges for the k vectors which are absolutely necessary for the solution of the problems of uniqueness and completeness of layer-group irreps. For this purpose it is advantageous to describe the different k-vector stars belonging to a Wyckoff position applying the so-called uni-arm description. Two k vectors of a Wyckoff position are called uni-arm if one can be obtained from the other by parameter variation. The description of k-vector stars of a Wyckoff position is called uni-arm if the k vectors representing these stars are uni-arm. Frequently, in order to achieve a uni-arm description, it is necessary to transform k vectors to equivalent ones. In addition, to enable a uni-arm description, symmetry lines outside the asymmetric unit may be selected as orbit representatives. Such a segment of a line is called a flagpole. Examples of flagpoles are displayed in Fig. 9, and Figs. 10 and 11 for the layer group cm2m (No. 35) [cf. Section 5.2[link] for a detailed account of the use of flagpoles in the uni-arm description of k-vector types that belong to Wyckoff position 2a of the reciprocal plane group (c1m1)*].

4. The Brillouin-zone database for layer groups

The k-vector data of the Brillouin-zone database of the BCS are accessed by the retrieval tool LKVEC which uses as input the ITE number of the layer group. The output consists essentially of wavevector tables and figures. There are several sets of figures and tables for the same layer group when its Brillouin-zone shape depends on the lattice parameters of the reciprocal lattice. The k-vector data are the same for layer groups of the same arithmetic crystal class.

In the k-vector tables, the wavevector data of LW are compared with the Wyckoff-position data of ITA. In the figures, the Brillouin zones and the representation domains of LW, and the asymmetric units, chosen often in analogy to those of ITA, are displayed.

LW describe the monoclinic/rectangular layer groups with respect to a setting that is different from the conventional one found in ITE. Using the relationship between the two settings, we have transformed the k-vector data of LW to the conventional setting of ITE for all six monoclinic/rectangular arithmetic crystal classes: 211p, 211c, m11p, m11c, 2/m11p and 2/m11c. The transformed special k-vector points, lines and planes keep the LW labels of the k vectors from which they were derived.

The URL of the program LKVEC is https://www.cryst.ehu.es/subperiodic/get_layer_kvec.html.

4.1. Guide to the tables

Each k-vector table is headed by the corresponding Hermann–Mauguin symbol of the layer group, its ITE number and the symbol of the arithmetic crystal class to which the layer group belongs. If there is more than one table for an arithmetic crystal class, then these tables refer to different geometric conditions for the lattice parameters that are indicated after the symbol of the arithmetic crystal class. The set of layer groups of the arithmetic crystal class are also indicated in the headline block. They are followed by the symbol of the corresponding reciprocal plane-group type together with the conditions for the lattice parameters of the reciprocal lattice, if any (asterisks denote reciprocal-space quantities). From the k-vector table there is a link to the corresponding Brillouin-zone figure.

The k-vector tables consist of two parts: (i) `Litvin & Wike' description and (ii) `Plane-group description'. The first three columns under the heading `Litvin & Wike' refer to the description of the k-vectors found in Tables 24 and 25 of LW. It consists of labels of k-vectors (column 1), their parameter descriptions (column 2) and their layer little co-group (column 3 for primitive lattices and column 4 for c-centred). Note that LW substitutes the Greek-character labels for the symmetry points and lines inside the Brillouin zone by a symbol consisting of two Roman characters, e.g. GM instead of Γ, LD instead of Λ etc. In order to enable the uni-arm description new k-vector types, equivalent to those of LW, are added to the k-vector lists. Equivalent k vectors (related by the sign ∼) are designated by the same labels; additional indices distinguish the new k vectors.

Different k vectors with the same LW label always belong to the same k-vector type, i.e. they correspond to the same Wyckoff position. k Vectors with different LW labels may either belong to the same or to different types of k vectors. When k vectors with different LW labels belong to the same k-vector type, the corresponding parameter descriptions are followed by the letters `ex' (from Latin, with the meaning of `from' or `out of'). Symmetry points or lines of symmetry of LW, related to the same Wyckoff position, are grouped together in a block. In the k-vector tables, neighbouring Wyckoff-position blocks are distinguished by a slight difference in the background colour. The parameter description of the uni-arm region of a k-vector type is shown in the last row of the corresponding Wyckoff-position block.

The wavevector coefficients of LW (column 2 of the k-vector tables) refer always to a primitive basis irrespective of whether the conventional description of the group in ITE is with respect to a centred or primitive basis. For that reason, for layer groups with centred lattices, the wavevector coefficients with respect to the usual conventional reciprocal basis, i.e. dual to the conventional centred basis, are listed in the column under the heading `Conventional' of the k-vector tables. The relations between the conventional coefficients (k1, k2) and the primitive coefficients ([{k}_{{\rm p}1}], [{k}_{{\rm p}2}]) are summarized in Table 1[link]. (For layer groups with primitive lattices, the wavevector coefficients referred to a primitive basis coincide with those referred to the basis dual to the conventional one of ITA.)

Table 1
`Conventional' k-vector coefficients kj (i.e. with respect to a basis dual to the conventional basis of ITA) expressed by the `primitive' k-vector coefficients [{k}_{{\rm p}j}] (i.e. referred to a primitive basis) for the different Bravais lattices in direct space

Bravais lattice k1 k2
mp, op, tp, hp [{k}_{{\rm p}1}] [{k}_{{\rm p}2}]
oc [{k}_{{\rm p}1}+ {k}_{{\rm p}2}] [-{k}_{{\rm p}1}+{k}_{{\rm p}2}]

The layer little co-group data of each k vector are listed under the heading `Layer little co-group' of the k-vector tables. The layer little co-groups are subgroups of the point group of the layer group and are described by oriented point-group symbols (as is customary for site-symmetry groups of Wyckoff positions).

The data for the crystallographic classification scheme of the wavevectors are listed under the heading `Plane-group description' in the k-vector tables. The columns `Wyckoff positions' show the `multiplicity', `Wyckoff letter' and `site-symmetry' of the Wyckoff positions of the corresponding symmorphic plane group [{\cal P}_{0}] of ITA which is isomorphic to the reciprocal plane group ([{\cal P}]) *. The multiplicity of a Wyckoff position divided by the number of lattice points in the conventional unit cell of ITA is equal to the number of arms of the star of the corresponding k vector. The alphabetical sequence of the Wyckoff positions determines the sequence of the LW labels. Unlike in ITA, the tables start with the Wyckoff letter `a' for the Wyckoff position of the highest site symmetry. Site symmetries are described by means of oriented point-group symbols which are also links to more detailed information on the symmetry operations of the site-symmetry group. Besides the shorthand description, the matrix-column representation and the geometric representation of the symmetry operations of the site-symmetry group, the program also provides a table with the relationship between the symmetry operations of the site-symmetry group and the layer little co-group (cf. Fig. 5[link] and Section 3.1[link] for detailed explanations).

[Figure 5]
Figure 5
Relation between the symmetry operations of the site-symmetry group .m. of the plane group c1m1 (No. 5) and the layer little co-group m2m of the layer group cm2m (No. 35) for the k vector DT = (−u, u).

The parameter description of the Wyckoff positions is shown in the last column of the wavevector tables under the heading `Coordinates'. It consists of a representative coordinate doublet of the Wyckoff position and algebraic statements for the description of the independent parameter range. In some cases, the algebraic expressions are substituted by the designation of the parameter region in order to avoid clumsy notation. Because of the dependence of the shape of the Brillouin zone on the lattice parameter relations there may be vertices of the Brillouin zone with a variable coordinate. If such a point is displayed and designated in the tables and figures by an upper-case letter, then the label of its variable coefficient used in the parameter-range descriptions is the same letter but typed in lower case.

Because of the isomorphism between [{\cal P}_{0}] and ([{\cal P}]) * the coordinate doublets of the Wyckoff positions of [{\cal P}_{0}] can be interpreted as k-vector coefficients ([{k}_{{\rm a}1}], [{k}_{{\rm a}2}]) determined with respect to the conventional ITA basis of [{\cal P}_{0}]. The relation between the ITA coefficients ([{k}_{{\rm a}1}], [{k}_{{\rm a}2}]) and the conventional coefficients ( k1, k2) is shown in Table 2[link].

Table 2
`Conventional' k-vector coefficients kj (i.e. with respect to a basis dual to the conventional basis of ITA) expressed by the ITA k-vector coefficients [{k}_{{\rm a}j}] (i.e. referred to the conventional ITA basis of [{\cal P}_{0}]) for the different Bravais lattices in direct space

Bravais lattice k1 k2
mp, op, tp [{k}_{{\rm a}1}] [{k}_{{\rm a}2}]
oc [{2k}_{{\rm a}1}] [{2k}_{{\rm a}2}]
hp [{k}_{{\rm a}1}-{k}_{{\rm a}2}] [{k}_{{\rm a}2}]

At the bottom of the web page with the k-vector table one finds an auxiliary tool which allows the complete characterization of a wavevector of the reciprocal space (not restricted to the first Brillouin zone): given the k-vector coefficients referred either to a primitive (LW) or to a conventional basis, the program assigns the k vector to the corresponding wavevector symmetry type, specifies its LW label, and calculates the layer little co-group and the arms of the k-vector star. Consider again the example of the k vector with coefficients [(-1.21,1)] for the layer group [p\bar{6}m2] (No. 78), cf. Fig. 6[link] and Fig. 8. It is a vector outside the first Brillouin zone and its coefficients do not correspond to any of the parameter descriptions of the k-vector representatives listed in Fig. 6[link]. The output of the auxiliary tool indicates that k [(-1.21,1)] is a point of a special k-vector line of type SM and belongs to the Wyckoff-position block 3c. As already commented in Section 3.1[link], its star consists of three k vectors, [{}^{*}{\bf k}] = {(−1.21, 1), (1, 0.21), (0.21, −1.21)}. The site-symmetry group ..m is generated by a reflection plane that can be identified by direct inspection among the symmetry operations of (p31m) *. The layer little co-group, however, is mm2, due to the additional reflection in the layer plane.

[Figure 6]
Figure 6
k-Vector table of the layer group [p\bar{6}m2] (No. 78) as shown on the BCS. The Brillouin-zone diagram is shown in Fig. 8.

4.2. Guide to the figures

The headline of each Brillouin-zone figure includes the same information as the k-vector tables: the Hermann–Mauguin symbol of the layer group, the ITE number and the symbol of the arithmetic crystal class to which the layer group belongs. Different figures for the same arithmetic crystal class are distinguished by the geometric conditions for the lattice. The corresponding conditions for the lattice parameters of the reciprocal lattice are indicated after the symbol of the reciprocal plane group.

The Brillouin zones are two-dimensional objects in the reciprocal space. The coordinate axes are designated by kx and ky, and the origin with coefficients (0, 0) always coincides with the centre of the Brillouin zone and is called Γ (indicated as GM in the k-vector tables). In the Brillouin-zone figures the representation domains of LW are compared with the asymmetric units of ITA. A statement of whether the representation domain of LW and the asymmetric unit are identical or not is given below the k-vector table. The asymmetric units are often not fully contained in the Brillouin zone, but protrude from it, in particular by flagpoles.

The representatives of the orbits of k-vector symmetry points or k-vector symmetry lines, as well as the edges of the representation domains of LW and of the asymmetric units are brought out in colour (see Fig. 7[link]):

[Figure 7]
Figure 7
Colour coding of points and lines applied in the Brillouin-zone diagrams.

(a) Symmetry points. A representative point of each orbit of special k-vector points is designated by a circle filled in red with its label also in red. Note that a point is coloured red only if it is really a special point, i.e. a point whose layer little co-group is a supergroup of the little co-groups of the points in its neighbourhood. In the figures, a point is marked by its label and an empty circle if it is listed in the corresponding k-vector table but is not a point of special symmetry. For example, points listed by LW are not coloured if they form part of a symmetry line or a symmetry plane. The same designation is used for the auxiliary points that have been added in order to facilitate the comparison between the traditional and the reciprocal plane-group descriptions of the k-vector types.

(b) Symmetry lines. A line is coloured in red with its label also in red only if it is a special k-vector line, i.e. the layer little co-groups of the points on the line are supergroups of the little co-groups of the points in its neighbourhood. The colour of the line is pink for an edge of the asymmetric unit which is not a symmetry line. The colour of the line is brown with the name in red for a line which is a symmetry line as well as an edge of the asymmetric unit. The edges of the representation domains are coloured light blue if the representation domain of LW does not coincide with the asymmetric unit. Edges of the representation domain and their labels are coloured dark blue if they are symmetry lines. Flagpoles are always coloured in red. Coordinate axes, edges of the Brillouin zone or auxiliary lines are displayed by thin solid black lines.

5. Examples

The relation between the traditional and the reciprocal group descriptions of the wavevector types is illustrated by the following examples. The figures and tables included here form part of the output of the access tool LKVEC.

5.1. k-Vector table and Brillouin zone for the layer group p62m (No. 78)

The k-vector table and the Brillouin-zone diagram of the hexagonal layer group [p\bar{6}m2] (No. 78) are shown in Fig. 6[link] and Fig. 8[link], respectively. The reciprocal lattice of a hexagonal p lattice is also a hexagonal p lattice and the Brillouin zone is a hexagon. The conventional basis for the reciprocal lattice has [\gamma^{*} = 60^{\circ}] while the ITA description of hexagonal layer groups is based on a basis aH, bH with [\gamma = 120^{\circ}]. In the Brillouin-zone diagrams, the axis kx is taken along aH while ky points out in the direction of aH + bH.

[Figure 8]
Figure 8
Brillouin zone and representation domain of LW, and asymmetric unit of the layer group [p\bar{6}m2] (No. 78). The reciprocal plane group is (p31m)* (No. 15). The representation domain of LW is different from the asymmetric unit.

The k vectors of [p\bar{6}m2] (No. 78) listed by LW (Fig. 6[link]) are distributed in four k-vector types: (i) the Wyckoff-position block 1a formed by the GM point, (ii) the block 2b formed by the K point, (iii) the k-vector point M and the k-vector lines SM and SN correspond to the block 3c, and (iv) the Wyckoff-position block 6d formed by the k-vector lines LD and T and the k-vector planes B and BB. The parameter description of a k-vector type is given in the last column of the table. Consider for example the line SM which, according to the ITA description, forms part of the k-vector type that is assigned to the Wyckoff position 3c with a site-symmetry group `..m'. Its parameter description x, 0: 0 < x < 1/2 indicates that the independent segment of the line x, 0 in the asymmetric unit is limited by the special k-vector points Γ (x = 0) and M (x = 1/2) with x varying between 0 and 1/2. The parameter descriptions of the uni-arm regions of the k-vector types are shown in the last row of the corresponding Wyckoff-position block. For example, in the block for position 3c, the k-vector line SN is equivalent (by a threefold rotation) to [x,0:-1/2\,\lt\,x\,\lt\,0] which in turn is equivalent (by a translation) to the line [x,0:1/2\,\lt\,x\,\lt\,1], denoted by SN1. This gives the uni-arm description M+SM+SN1 for the Wyckoff position 3c in Fig. 6[link].

As the asymmetric unit and the representation domain do not coincide, their edges are coloured in pink and light blue, respectively. It has already been pointed out that k-vector points and lines are brought out in red only if they are special k-vector points and lines. For example, the lines T and LD (indicated as Λ on Fig. 8[link]) are not coloured as special lines since they belong to the k-vector type of the Wyckoff-position block 6d, and their symmetry coincides with that of the neighbouring points of the symmetry plane B = [GM K M]. The points GM (indicated as Γ on Fig. 8[link]) and K, however, are represented by red circles as they are special k-vector points (cf. Fig. 6[link]). Likewise, the line SM (indicated as Σ on Fig. 8[link]) is coloured in brown because it is an edge of the asymmetric unit and at the same time is a special k-vector symmetry line. The k-vector line SN is coloured in dark blue as it is a special symmetry line along the edge of the representation domain. As already indicated, the k-vector line SN together with SM and the point M belong to the special k-vector type of the Wyckoff-position block 3c, i.e. all these different wavevectors belong to the same k-vector type. Although M is explicitly listed by LW as a special k-vector point, it is represented by an open circle in Fig. 8[link]: in fact, it joins the symmetry lines SM and SN1 to a continuous line as its little co-group type coincides with those of the points on the lines.

5.2. k-Vector table and Brillouin zone for the layer group cm2m (No. 35)

The layer group cm2m (No. 35) is an example of ortho­rhombic layer groups with a c-centred lattice. It belongs to the arithmetic crystal class m2mc which also includes the layer group cm2e (No. 36). The k-vector tables and Brillouin-zone figures for the layer groups belonging to the arithmetic crystal class m2mc are given in Fig. 9[link], and Figs. 10[link] and 11[link], respectively.

[Figure 9]
Figure 9
k-Vector tables of the layer groups with arithmetic crystal class m2mc, acute (a > b) and obtuse (a < b) cases, as shown on the BCS. The Brillouin-zone diagrams are shown in Figs. 10[link] and 11[link].
[Figure 10]
Figure 10
Brillouin zone and representation domain of LW, and asymmetric unit of the arithmetic crystal class m2mc: layer group cm2m (No. 35) and cm2e (No. 36) acute case (a > b). The reciprocal group is (c1m1)* (No. 5). The representation domain of LW does not coincide with the asymmetric unit.
[Figure 11]
Figure 11
Brillouin zone and representation domain of LW, and asymmetric unit of the arithmetic crystal class m2mc: layer group cm2m (No. 35) and cm2e (No. 36) obtuse case (a < b). The reciprocal group is (c1m1)* (No. 5). The representation domain of LW does not coincide with the asymmetric unit.

The k-vector tables show all special wavevectors with their coefficients and layer little co-groups as specified in Tables 24 and 25 of LW. The wavevector coefficients with respect to the conventional reciprocal basis, i.e. dual to the conventional centred basis, are listed in the column under the heading `Conventional' of the k-vector tables. For example, a k-vector point of the DT line (Fig. 9[link]) with primitive coefficients [(-1/4,1/4)] is described as (0,1/2) with respect to a basis dual to the conventional basis of cm2m.

The comparison of the wavevector list of LW and the reciprocal plane-group description indicates clearly the redundancy of most of the k vectors given by LW. In fact, for the derivation of a complete set of irreps it is necessary to consider just two k-vector types: a general one corresponding to the general Wyckoff position 4b, and the k-vector symmetry line related to the special Wyckoff position 2a. The large number of additional k vectors given in the tables of LW are due to two main reasons:

(a) The more symmetry a layer group has lost compared with its holosymmetric layer group (layer groups whose point groups are holohedral), the more k vectors are introduced in LW. In the case of cm2m, the holosymmetric layer group is cmmm (No. 47), and the lines DA and FA are examples of such additional k vectors. In most cases these additional k vectors can be avoided by extending the parameter range in the k-vector space.

(b) In the transition from a holosymmetric [{\cal{L}}_{\rm h}] to a non-holosymmetric layer group [{\cal L}], the order of the little co-group of a special k vector in [{\cal{L}}_{\rm h}] may be reduced in [{\cal L}] and, as a result, the special k vector in [{\cal{L}}_{\rm h}] may lose its `special nature' in [{\cal L}]. Such k vectors become part of a more general k-vector type (i.e. assigned to a Wyckoff position of lower site symmetry) and can be described by an extension of the corresponding parameter range. Consider, for example, the k-vector points Γ and Y in the tables of cm2m. They are special k-vector points lying on the reflection plane of m010 in the holosymmetric layer group cmmm, but in (the non-holosymmetric group) cm2m which does not contain m010 the two points form part of the special k-vector line with the uni-arm description: [0,y:-1/2\,\lt\,y\leq 1/2] (cf. the last row of the 2a Wyckoff-position block, Fig. 9[link]).

The Brillouin zone of the layer groups of the arithmetic crystal class m2mc is a non-regular hexagon (cf. Figs. 10[link] and 11[link]). Depending on the relation between the lattice parameters a and b, two topologically different Brillouin zones are to be distinguished: (i) the acute case with [a\,\gt\,b] and (ii) the obtuse case with [a\,\lt\,b]. Because of the reflection [-x,y] (with normal kx) of the reciprocal plane group (c1m1)*, the representation domain is only one half of the hexagon: for example, in the acute case (Fig. 10[link]) it is the trapezium with vertices [\Delta_{0},\ F_{0},\ F_{2},\ \Delta_{2}] (light blue boundary). The asymmetric unit is different from the representation domain: it is the rectangle with vertices J2J4V4V2 (pink boundary), i.e. the points [x,\ y:\ 0\leq x\leq{{1} / {2}},\ -{{1} / {4}}\leq y\leq{{1} / {4}}]. While the representation domains of the acute and obtuse unit cells have the more complicated form of a trapezium, the asymmetric units in both cases have the topologically identical and relatively simple shape of a rectangle.

As already indicated, the points Γ, Y2 and S (acute case), and Γ, Y and S (obtuse case) are not special k-vector points but form part of special lines and planes and in the diagrams they are represented by open circles. The line SM is not a symmetry line and is represented by a thin black line because it is located inside the Brillouin zone. The lines DT and DA are coloured in brown because they are symmetry lines and at the same time are edges of the asymmetric unit. Parts of DT and DA are also coloured in red because they correspond to flagpoles. The k-vector lines F and FA (acute case) are coloured in dark blue as they are symmetry lines along the edges of the representation domain.

Because of the special shape of the Brillouin zone and the representation domain for the acute case ([a\,\gt\,b]), the special k-vector line corresponding to the Wyckoff-position block 2a splits into several segments: the lines DT and DA, located inside the Brillouin zone, and the lines F and FA (coloured dark blue) at the border of the Brillouin zone (cf. Figs. 9[link] and 10[link]). For the description of the end points of the segments, it is necessary to introduce additional parameters as dt0 and f0 whose values depend on the specific relations between the lattice parameters. [In fact, the vertices [\Delta_{0}] and F0 of the representation domain have the coordinates [\Delta_{0}:] 0, 1 / 4+1 / 4(a* / b*) and F0: [0,\ {{1} / {4}}-{{1} / {4}}(a^{*} / b^{*})].] The use of flagpoles enables the uni-arm description: the flagpole [J2 Y2] [0,y:1/4\,\lt\,y\,\lt\,1/2] is equivalent to the segment [V4 Y] [1/2,y:-1/4\,\lt\,y\,\lt\,0] and the flagpole [V2 Y4] 0,y: [-1/2\,\lt\,y\,\lt\,-1/4] is equivalent to the segment [Y J4] [1/2,y:0\,\lt\,y\,\lt\,1/4]. The uni-arm description of the k-vector type of the Wyckoff position 2a is shown in the last row of the Wyckoff-position block and it is formed by the union of the points GM and Y2, the lines DT, DA, DT1(∼FA) and DA1 (∼F). Its parameter description (0, y) with y varying in the range (−1/2, 1/2) coincides with that of the acute case. The parameter description of the flagpole and its parameter range with respect to the basis of the reciprocal group are given below the k-vector table.

6. Conclusions

In this paper we have presented the layer-group crystallographic and wavevector databases of the BCS, together with the programs which give access to these data. Like the rest of the programs on the server, these tools are freely available and can be accessed via user-friendly web interfaces. The programs GENPOS, WYCKPOS and MAXSUB provide access to generators/general positions, Wyckoff positions and maximal subgroup information, respectively.

The wavevector database which contains the Brillouin-zone figures and wavevector tables for all 80 layer groups was recently implemented on the BCS. One can access the database through the program LKVEC. In this compilation, the representation domains and lists of special k vectors in the tables on layer-group representations by Litvin & Wike (1991[Litvin, D. B. & Wike, T. R. (1991). Character Tables and Compatibility Relations of the Eighty Layer Groups and Seventeen Plane Groups. New York: Plenum Press.]) are compared with the figures and wavevector data derived applying the reciprocal-space-group approach. This new database provides a solution to the completeness problem of layer-group representations by specifying the independent parameter ranges of general and special k vectors within the representation domains.

Acknowledgements

Open access funding enabled and organized by Projekt DEAL.

Funding information

This work has been supported by the Spanish Ministry of Science and Innovation (project MINECOG19/P21) and the Government of the Basque Country (project IT1301-19).

References

First citationAroyo, M. I. (2016). Editor. International Tables for Crystallography, Vol. A, Space-group Symmetry, 6th ed. Chichester: Wiley.   Google Scholar
First citationAroyo, M. I., Orobengoa, D., de la Flor, G., Tasci, E. S., Perez-Mato, J. M. & Wondratschek, H. (2014). Acta Cryst. A70, 126–137.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationAroyo, M. I., Perez-Mato, J. M., Capillas, C., Kroumova, E., Ivantchev, S., Madariaga, G., Kirov, A. & Wondratschek, H. (2006). Z. Kristallogr. 221, 15–27.  Web of Science CrossRef CAS Google Scholar
First citationAroyo, M. I., Perez-Mato, J. M., Orobengoa, D., Tasci, E. S., de la Flor, G. & Kirov, A. (2011). Bulg. Chem. Commun. 43, 183–197.  CAS Google Scholar
First citationAroyo, M. I. & Wondratschek, H. (1995). Z. Kristallogr. 210, 243–254.  CrossRef CAS Web of Science Google Scholar
First citationBoyle, L. L. (1986). Proceedings of the 14th International Colloquium on Group Theoretical Methods in Physics, pp. 405–408. Singapore: World Scientific.  Google Scholar
First citationBradley, C. J. & Cracknell, A. P. (1972). The Mathematical Theory of Symmetry in Solids: Representation Theory for Point Groups and Space Groups. Oxford: Clarendon Press.  Google Scholar
First citationCracknell, A. P., Davies, B. L., Miller, S. C. & Love, W. F. (1979). Kronecker Product Tables, Vol. 1, General Introduction and Tables of Irreducible Representations of Space Groups. New York: IFI/Plenum.  Google Scholar
First citationDresselhaus, M. S., Dresselhaus, G. & Jorio, A. (2008). Group Theory Application to the Physics of Condensed Matter. Berlin: Springer.  Google Scholar
First citationElcoro, L., Bradlyn, B., Wang, Z., Vergniory, M. G., Cano, J., Felser, C., Bernevig, B. A., Orobengoa, D., de la Flor, G. & Aroyo, M. I. (2017). J. Appl. Cryst. 50, 1457–1477.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationFlor, G. de la, Orobengoa, D., Evarestov, R. A., Kitaev, Y. E., Tasci, E. & Aroyo, M. I. (2019). J. Appl. Cryst. 52, 1214–1221.  Web of Science CrossRef IUCr Journals Google Scholar
First citationGallego, S. V., Etxebarria, J., Elcoro, L., Tasci, E. S. & Perez-Mato, J. M. (2019). Acta Cryst. A75, 438–447.  Web of Science CrossRef IUCr Journals Google Scholar
First citationHatch, D. M. & Stokes, H. T. (1986). Phase Transit. 7, 87–279.  CrossRef Google Scholar
First citationIpatova, I. P. & Kitaev, Y. E. (1985). Prog. Surf. Sci. 18, 189–246.  CrossRef CAS Web of Science Google Scholar
First citationKopský, V. & Litvin, D. (2010). Editors. International Tables for Crystallography, Vol. E, Subperiodic Groups, 2nd ed. Chichester: Wiley.   Google Scholar
First citationLitvin, D. B. & Kopský, V. (2014). Acta Cryst. A70, 677–678.  Web of Science CrossRef IUCr Journals Google Scholar
First citationLitvin, D. B. & Wike, T. R. (1991). Character Tables and Compatibility Relations of the Eighty Layer Groups and Seventeen Plane Groups. New York: Plenum Press.  Google Scholar
First citationMiller, S. C. & Love, W. F. (1967). Tables of Irreducible Representations of Space Groups and Co-Representations of Magnetic Space Groups. Boulder: Pruett Press.  Google Scholar
First citationMilošević, I., Nikolić, B., Damnjanović, M. & Krčmar, M. (2008). J. Phys. A, 31, 3625–3648.  Google Scholar
First citationStokes, H. T. & Hatch, D. M. (1988). Isotropy Subgroups of the 230 Crystallographic Space Groups. Singapore: World Scientific.  Google Scholar
First citationTasci, E. S., de la Flor, G., Orobengoa, D., Capillas, C., Perez-Mato, J. M. & Aroyo, M. I. (2012). EJP Web Conf. 22, 00009.  Google Scholar
First citationWintgen, G. (1941). Math. Ann. 118, 195–215.  CrossRef CAS Google Scholar
First citationWondratschek, H. & Müller, U. (2010). Editors. International Tables for Crystallography, Vol. A1, Symmetry Relations between Space Groups, 2nd ed. Chichester: John Wiley & Sons.  Google Scholar

This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.

Journal logoFOUNDATIONS
ADVANCES
ISSN: 2053-2733
Follow Acta Cryst. A
Sign up for e-alerts
Follow Acta Cryst. on Twitter
Follow us on facebook
Sign up for RSS feeds