---
abstract: 'A new section of databases and programs devoted to double crystallographic groups (point and space groups) has been implemented in the Bilbao Crystallographic Server (http://www.cryst.ehu.es). The double crystallographic groups are required in the study of physical systems whose Hamiltonian includes spin-dependent terms. In the symmetry analysis of such systems, instead of the irreducible representations of the space groups, it is necessary to consider the single- and double-valued irreducible representations of the double space groups. The new section includes databases of symmetry operations ([DGENPOS]{}) and of irreducible representations of the double (point and space) groups ([REPRESENTATIONS DPG]{} and [REPRESENTATIONS DSG]{}). The tool [DCOMPREL]{} provides compatibility relations between the irreducible representations of double space groups at different $\mathbf{k}$-vectors of the Brillouin zone when there is a group-subgroup relation between the corresponding little groups. The program [DSITESYM]{} implements the so-called site-symmetry approach, which establishes symmetry relations between localized and extended crystal states, using representations of the double groups. As an application of this approach, the program [BANDREP]{} calculates the band representations and the elementary band representations induced from any Wyckoff position of any of the 230 double space groups, giving information about the properties of these bands. Recently, the results of [BANDREP]{} have been extensively applied in the description and the search of topological insulators.'
author:
- Luis Elcoro
- Barry Bradlyn
- Zhijun Wang
- 'Maia G. Vergniory'
- Jennifer Cano
- Claudia Felser
- 'B. Andrei Bernevig'
- Danel Orobengoa
- Gemma de la Flor
- 'Mois I. Aroyo'
bibliography:
- 'elcoro.bib'
title: Double crystallographic groups and their representations on the Bilbao Crystallographic Server
---

Introduction
============

The *Bilbao Crystallographic Server* (http://www.cryst.ehu.es) website offers crystallographic databases and programs [@aroyo2006]. It can be used free of charge from any computer with a web browser via Internet. The applications on the server are organized into different sections depending on their degree of complexity, in such a way that the most complex tools make use of the results obtained by the simpler ones. The server is built on a core of databases that includes data from the *International Tables for Crystallography*, Vol. A: *Space- group symmetry* (, henceforth abbreviated as [[*IT*]{}A ]{}), Vol. A1: *Symmetry Relations between Space Groups* [@ita1] and Vol. E: *Subperiodic groups* [@ite]. A **k**-vector database with Brillouin-zone figures and classification tables of all the wave vectors for all 230 space groups is also available. Databases of magnetic space groups and magnetic structures have recently been implemented, alongside with a set of computational tools that facilitate the systematic application of symmetry arguments in the study of magnetic materials [@perezmato2015]. The database of incommensurate structures, hosted by the server, contains both single-modulated structures and composites. Parallel to the databases and the crystallographic software there are a number of programs facilitating the study of specific problems related to solid-state physics, structural chemistry and crystallography involving crystallographic groups, their group-subgroup relations and irreducible representations.

In a number of physical applications it is necessary to include spin-dependent terms in the Hamiltonian of a crystal system: for example, in taking into account relativistic effects in band structure calculations or spin-orbit coupling in crystal-field theory. Then, instead of the crystallographic groups and their representations, it is necessary to consider the so called *double crystallographic groups* and the related *single-valued* or *vector* and *double-valued* or *spinor* representations. Since their introduction during the first half of the last century [@bethe1929; @opechowski1940], the double groups and their representations have been discussed in detail in the literature (see for example, the account given by and the references therein). Several reference compilations of the character tables of double point groups (*cf.* ), and of the irreducible representations of double space groups [@cracknell1979], abbreviated as CDML exist. However, we are not aware of any online available databases of double crystallographic groups nor of any computing tools for the calculation of their representations. Here, we report on the recently completed development of such programs and their implementation on the Bilbao Crystallographic Server (BCS).

For what follows, it will be convenient for us to recall briefly some essential features of the double crystallographic groups. Suppose that $\overline{{\mbox{$\mathcal{G}$}}}$ is a point group which consists of pure rotations, *i.e.* a subgroup of the group $SO(3)$. The two-to-one homomorphism $\varphi$ between the group $SU(2)$ of all $(2\times2)$ unitary unimodular matrices onto $SO(3)$ can be used for the formal definition of the double groups $\overline{^d{\mbox{$\mathcal{G}$}}}$ [@opechowski1940]: the double group $\overline{^d{\mbox{$\mathcal{G}$}}}$ of a group $\overline{{\mbox{$\mathcal{G}$}}}$ of order $|\overline{{\mbox{$\mathcal{G}$}}}|$ which is a subgroup of the group $SO(3)$, is the abstract group of order 2$|\overline{{\mbox{$\mathcal{G}$}}}|$ having the same multiplication table as the 2$|\overline{{\mbox{$\mathcal{G}$}}}|$ matrices of $SU(2)$ which correspond under $\varphi$ to the elements of the group $\overline{{\mbox{$\mathcal{G}$}}}$. The assignment of the $SU(2)$ matrices to the crystallographic symmetry operations used in the databases and programs of the server follows the choice by .

The kernel of the homomorphism $\varphi:SU(2)\longrightarrow SO(3)$ consists of two elements: $\{1=\begin{pmatrix}1&0\\0&1\end{pmatrix}, {^d1}=\begin{pmatrix}-1&0\\0&-1\end{pmatrix}\}$. As a result, the preimage of each element $R$ of $\overline{{\mbox{$\mathcal{G}$}}}$ consist of two elements of $\overline{^d{\mbox{$\mathcal{G}$}}}$, namely $R$ and $^dR={^d1}R$. Formally, one can write $$\label{def_double_group}
\overline{^d{\mbox{$\mathcal{G}$}}}=\{R\}\cup\{^dR\}.$$ It is important to note that the subset of elements {R} does not form a subgroup of $\overline{^d{\mbox{$\mathcal{G}$}}}$ as it is not closed under the binary group operation. For instance, if $R$ represents a 2-fold rotation or a mirror plane, $R^2=\,^d1$. In that sense it is wrong to refer to eq. (\[def\_double\_group\]) as a coset decomposition of $\overline{^d{\mbox{$\mathcal{G}$}}}$ with respect to $\overline{{\mbox{$\mathcal{G}$}}}$. Physically, it is considered that the $SU(2)$ matrices act on the spinors of a 1/2 spin space. Then, the operation $^d1$ is of order 2 and it is often interpreted as a $2\pi$ rotation.

The above statements about double point groups of pure rotations can be generalized in a straightforward way to the case of double point groups that contain improper rotations. For this purpose it is sufficient to indicate that the two $SU(2)$ matrices $\overline{1}=\begin{pmatrix}1&0\\0&1\end{pmatrix}$ and $^d\overline{1}={^d1}\overline{1}=\begin{pmatrix}-1&0\\0&-1\end{pmatrix}$ are assigned to the symmetry operation of inversion.

Consider a space group and its decomposition ${\mbox{$\mathcal{G}$}}:{\mbox{$\mathcal{T}$}}$ into cosets with respect to its translation subgroup , i.e. ${\mbox{$\mathcal{G}$}}={\mbox{$\mathcal{T}$}}\cup {\mbox{$\mathcal{T}$}}\{R_2|\mathbf{v}_2\} \cup \dots \cup {{\mbox{$\mathcal{T}$}}}\{R_n|\mathbf{v}_n\}$. In a similar way, the symmetry operations of the double space group $^d{\mbox{$\mathcal{G}$}}$ can be conveniently represented using its coset decomposition with respect to as: $$\label{coset_double}
^d{\mbox{$\mathcal{G}$}}={{\mbox{$\mathcal{T}$}}}\cup {{\mbox{$\mathcal{T}$}}}\{R_2|\mathbf{v}_2\} \cup \dots \cup {{\mbox{$\mathcal{T}$}}}\{R_n|\mathbf{v}_n\}\cup {{\mbox{$\mathcal{T}$}}}\{^d1|\mathbf{o}\}\cup\cdots\cup {{\mbox{$\mathcal{T}$}}}\{{^dR}_n|\mathbf{v}_n\},$$ where $\{R_1|\mathbf{v}_1\}=\{1|\mathbf{o}\}$ is omitted. Here, $R_i$ and ${^dR}_i$ are the elements of the double point group $\overline{^d{\mbox{$\mathcal{G}$}}}$ of $^d{\mbox{$\mathcal{G}$}}$. Consequently, there are two elements $\{R_i|\mathbf{v}_i\}$ and $\{^dR_i|\mathbf{v}_i\}$ of the double space group $^d{\mbox{$\mathcal{G}$}}$, that correspond to every element of the space group ${\mbox{$\mathcal{G}$}}$. The translation subgroup is an invariant subgroup of the double space group $^d{\mbox{$\mathcal{G}$}}$.

In the following, we shall discuss the development and implementation of databases and programs involving the double crystallographic groups for the BCS. We start with the presentation of the double space group database and the retrieval tools that access the stored crystallographic symmetry information (Section \[sec:doublegroups\]). The introduction to the basic programs available on BCS for the computation of representations of double crystallographic groups is given in Section \[sec:representations\]. The last sections of the article are devoted to the presentation of the accompanying applications of the representations of double crystallographic groups, such as their compatibility relations (Section \[sec:compatibilityrelations\]), the site-symmetry approach (Section \[sec:sitesymmetry\]), and finally, the determination of the band representations and the elementary band representations (Section \[sec:bandrepresentations\]). In the Appendix we briefly describe the normal-subgroup induction procedure.

\[sec:doublegroups\]Double crystallographic groups
==================================================

The double space groups are infinite groups, *i.e.* they contain an infinite number of symmetry operations generated by the set of all translations of the space group. As already noted, a practical way to represent the symmetry operations of the double space group $^d{\mbox{$\mathcal{G}$}}$ is based on the coset decomposition of $^d{\mbox{$\mathcal{G}$}}$ with respect to its translation subgroup ${\mbox{$\mathcal{T}$}}$, *cf.* eq.(\[coset\_double\]). The set of coset representatives $\{\{R_i|\mathbf{v}_i\},\{^dR_i|\mathbf{v}_i\}, i=1,\dots, n\}$ of the decomposition $^d{\mbox{$\mathcal{G}$}}:{\mbox{$\mathcal{T}$}}$ (often referred to as *General positions* of $^d{\mbox{$\mathcal{G}$}}$) represents, in a clear and compact way, the infinite number of symmetry operations of the double space group $^d{\mbox{$\mathcal{G}$}}$. The infinite symmetry operations in a coset have the same linear part $R_i$ (or $^dR_i$) while their translation parts differ by lattice translations. The translations $\{1|\mathbf{t}\}\in{\mbox{$\mathcal{T}$}}$ form the first coset with the identity $\{1|\mathbf{o}\}$ as a coset representative. The number of cosets is always finite and is equal to the order of the double point group $\overline{^d{{\mbox{$\mathcal{G}$}}}}$ of the double space group $^d{\mbox{$\mathcal{G}$}}$.

The database of double space groups, available on BCS, includes the lists of the representatives of the general positions of each double space group $^d{\mbox{$\mathcal{G}$}}$. These data can be retrieved using the [DGENPOS]{} tool (http://www.cryst.ehu.es/cryst/dgenpos), by specifying the sequential [[*IT*]{}A ]{}number or the Hermann-Mauguin symbol of the space group corresponding to $^d{\mbox{$\mathcal{G}$}}$. The symmetry operations are specified by their matrix presentations, *i.e* by the $(3\times4)$ matrix-column pairs and the $(2\times2)$ matrices of $SU(2)$, and by Seitz symbols.

-   Matrix presentations

    With reference to a coordinate system (O, $\mathbf{a}_1,\mathbf{a}_2, \mathbf{a}_3$), consisting of an origin O and a basis $\mathbf{a}_k$, the symmetry operations of the space group $^d{\mbox{$\mathcal{G}$}}$ are described by matrix-column pairs $({\mbox{\bfseries{\itshape{W}}}},\,{\mbox{\bfseries{\itshape{w}}}})$. The set of translations are represented by the matrix-column pairs (, ${\mbox{\bfseries{\itshape{t}}}}_i$), where ${\mbox{\bfseries{\itshape{I}}}}$ is the ($3\times 3$) unit matrix and ${\mbox{\bfseries{\itshape{t}}}}_i$ is the column of coefficients of the translation vector $\mathbf{t}_i$ that belongs to the [*vector lattice*]{} $\mathbf{L}$ of $^d{\mbox{$\mathcal{G}$}}$. The programs and databases of the double space groups, as well the rest of the computer tools on BCS, use specific settings of space groups (hereafter referred to as standard, or default, settings) that coincide with the conventional space-group descriptions found in [[*IT*]{}A ]{}. For space groups with more than one description in [[*IT*]{}A ]{}, the following settings are chosen as standard: *unique axis b* setting, *cell choice 1* for monoclinic groups, *hexagonal axes* setting for rhombohedral groups, and *origin choice 2* (origin on $\overline{1}$ ) for the centrosymmetric groups listed with respect to two origins in [[*IT*]{}A ]{}.

    The pair of symmetry operations $\{\{R_i|\mathbf{v}_i\},\{^dR_i|\mathbf{v}_i\}\}$ related by $\{^d1|\mathbf{o}\}$ have the same matrix-column presentation $({\mbox{\bfseries{\itshape{W}}}}_i,\,{\mbox{\bfseries{\itshape{w}}}}_i)$ while their $(2\times2)$ matrices differ by a sign. As already noted, the chosen correspondence scheme between the rotations and the $(2\times2)$ matrices follows closely the choice made by . The shorthand descriptions of the $(3\times4)$ matrix-column pairs and of the $(2\times2)$ matrices of the symmetry operations are also shown in the general-position table.

-   Seitz symbols

    The Seitz symbols $\{R_i|\mathbf{v}_i\}$ of space-group symmetry operations consist of two parts: a rotation (or linear) part $R$ and a translation part $\mathbf{v}$ [@glazer2014]. The Seitz symbol is specified between braces and the rotational and the translational parts are separated by a vertical line. The translation parts **v** correspond exactly to the columns ${\mbox{\bfseries{\itshape{w}}}}$ of the matrix-column presentation $({\mbox{\bfseries{\itshape{W}}}},{\mbox{\bfseries{\itshape{w}}}})$ of the symmetry operations. The rotation parts $R$ consist of symbols that specify (i) the type and the order of the symmetry operation, and (ii) the orientation of the corresponding symmetry element with respect to the space-group basis. The symbols 1 and $\overline{1}$ are used for the identity and the inversion, $m$ for reflections, the symbols 2, 3, 4 and 6 denote rotations and $\overline{3}$ , $\overline{4}$ and $\overline{6}$ rotoinversions. For rotations and rotoinversions of order higher than 2, a superscript + or - is used to indicate the sense of the rotation. The orientation of the symmetry element is denoted by the direction of the axis for rotations or rotoinversions, or by the direction of the normal to reflection planes. Subscripts of the symbols specify the characteristic direction of the operation: for example, the subscripts $100$, $010$ and $1\overline{1} 0$ refer to the directions $[100]$, $[010]$ and $[1\overline{1} 0]$, respectively.

    The symmetry operations of the double space groups are denoted by the *modified* Seitz symbols. The *modified* Seitz symbols include a superscript $d$ added to the rotational part of the symmetry operations to distinguish between the symmetry operations $\{R|\mathbf{v}\}$ and those obtained by their combinations with $\{^d1|\mathbf{o}\}$.

As an example, consider the general-position table of the double space group ${\mathit{P}2_12_12_1}$ (No. 19) shown in Fig. \[fig:dgenpos\]. The listed eight symmetry operations are the chosen coset representatives of the decomposition of ${\mathit{P}2_12_12_1}$ with respect to its translation subgroup. The pairs of elements related by $\{^d1|\mathbf{o}\}$ are clearly distinguished by their Seitz symbols $\{R_i|\mathbf{v}_i\}$ and $\{^dR_i|\mathbf{v}_i\}$, *e.g.* $(2) \{2_{001}|1/2,0,1/2\}$ and $(6) \{^d2_{001}|1/2,0,1/2\}$. As explained, the pairs of symmetry operations $\{\{R_i|\mathbf{v}_i\},\ \{^dR_i|\mathbf{v}_i\}\}$ are represented by the same $3\times4$ matrix-column pair while their $2\times2$ matrices differ by a sign.

![Figure 1: Screenshot of the output of the program [DGENPOS]{}. The figure shows the general-position table of the double space group ${\mathit{P}2_12_12_1}$ (No. 19). The symmetry operations are specified by their matrix representations, shorthand notation (coordinate triplets and spin components) and Seitz symbols.[]{data-label="fig:dgenpos"}](./Figure01.png){width="60.00000%"}

\[sec:representations\]Representations of the double crystallographic groups
============================================================================

The problem
-----------

There are two programs on the Bilbao Crystallographic Server that compute the irreducible representations (irreps) of space groups explicitly, namely [REPRES]{} (http://www.cryst.ehu.es/cryst/repres) and [Representations SG]{} (http://www.cryst.ehu.es/cryst/ representationsSG). Given a space group $\cal{G}$ and a [**k**]{}-vector, both programs calculate the irreps of space groups following the algorithm based on a normal-subgroup induction method, *i.e.* the irreps of a group $\cal{G}$ are calculated starting from those of a normal subgroup $\cal{H}\lhd\cal{G}$. The main steps of the procedure involve the construction of all irreps of $\cal{H}$ and their distribution into orbits under $\cal{G}$, determination of the corresponding little groups and the allowed (small) irreps and finally, construction of the irreps of $\cal{G}$ by induction from the allowed irreps. The labels assigned to the irreps calculated by [REPRES]{} and [Representations SG]{} correspond to those used by CDML. The correct assignment of the labels to the irreps calculated by [REPRES]{} is achieved with the help of the database of space-group irreps recently developed by .

Here we report on the development of the program [Representations DSG]{} (http://www.cryst.ehu.es/cryst/ representationsDSG), available on BCS, for the computation of the irreps of double crystallographic groups. The method for calculating of the irreps of double crystallographic groups is based on a generalization of the normal-subgroup induction procedure for the calculation of space-group irreps implemented in the programs [REPRES]{} and [Representations SG]{}.

To make the exposition self consistent, we include in the Appendix the main concepts and statements of the normal-subgroup induction procedure (for a detailed presentation, *cf*. ).

\[irrep\_method\]The method
---------------------------

The normal-subgroup induction procedure for the calculation of space-group representations (see the Appendix) can be generalized in a straightforward way for the calculations of the representations of the double space groups (*cf.* Miller & Love, 1967). Consider the coset decomposition of the double group $^d{\mbox{$\mathcal{G}$}}$ with respect to its translation subgroup (*cf.* eq. (\[coset\_double\])):

$$^d{\mbox{$\mathcal{G}$}}=\bigcup_i{\mbox{$\mathcal{T}$}} {\mbox{\textit{\textsf{q}}}}_i \cup {\mbox{$\mathcal{T}$}}\,{^d{\mbox{\textit{\textsf{q}}}}}_i.$$

It will be convenient to rewrite the coset decomposition of $^d{\mbox{$\mathcal{G}$}}$ with respect to the group $^d{\mbox{$\mathcal{T}$}}={\mbox{$\mathcal{T}$}}\otimes\{\{1|000\},\{^d1|000\}\}$: $$^d{\mbox{$\mathcal{G}$}}={^d{\mbox{$\mathcal{T}$}}}\cup {^d{\mbox{$\mathcal{T}$}}}{\mbox{\textit{\textsf{q}}}}_2 \cup \cdots \cup {^d{\mbox{$\mathcal{T}$}}}{\mbox{\textit{\textsf{q}}}}_m.$$

The group ${^d{\mbox{$\mathcal{T}$}}}$ is a trivial central extension of ${\mbox{$\mathcal{T}$}}$ by the group generated by $\{^d1|000\}$; it is an abelian group, and it is a normal subgroup of $^d{\mbox{$\mathcal{G}$}}$, *i.e.* $^d{\mbox{$\mathcal{G}$}}\triangleright {^d{\mbox{$\mathcal{T}$}}}$. Each irrep $\mathbf{\Gamma}^\textbf{k}$ of generates two irreps $\mathbf{\Gamma}^{\mathbf{k}}$ and $\overline{\mathbf{\Gamma}}^{\mathbf{k}}$ of ${^d{\mbox{$\mathcal{T}$}}}$:

$$\begin{aligned}
\mathbf{\Gamma}^{\mathbf{k}}(\{1|\mathbf{t}\})=\mathrm{exp}(-i\textbf{k}\cdot\textbf{t})\textrm{,}&\hspace{0.15cm}\mathbf{\Gamma}^{\mathbf{k}}(\{^d1|\mathbf{t}\})=\mathrm{exp}(-i\textbf{k}\cdot\textbf{t}),\\
\overline{\mathbf{\Gamma}}^{\mathbf{k}}(\{1|\mathbf{t}\})=\mathrm{exp}(-i\textbf{k}\cdot\textbf{t})\textrm{,}&\hspace{0.15cm}\overline{\mathbf{\Gamma}}^{\mathbf{k}}(\{^d1|\mathbf{t}\})=-\mathrm{exp}(-i\textbf{k}\cdot\textbf{t}).\end{aligned}$$

The irrep $\mathbf{\Gamma}^{\mathbf{k}}$ is known as a *single-valued* irrep: the same (one-dimensional) matrix represents the elements $\{1|\mathbf{t}\}$ and $\{^d1|\mathbf{t}\}$. On the contrary, the (one-dimensional) matrices of the elements $\{1|\mathbf{t}\}$ and $\{{^d1}|\mathbf{t}\}$ differ by a sign in the *double-valued* irrep $\overline{\mathbf{\Gamma}}^{\mathbf{k}}$.

Note that the term single-valued (double-valued) representation is kept for any representation $\mathbf{D}^{\mathbf{k}} (\overline{\mathbf{D}}^{\mathbf{k}})$ induced from $\mathbf{\Gamma}^{\mathbf{k}} (\overline{\mathbf{\Gamma}}^{\mathbf{k}})$ as all such representations are characterized by the same relationship between the matrices representing the elements related by the operation ${^d{\mbox{\textit{\textsf{1}}}}}$, namely:

$$\begin{aligned}
\mathbf{D}^{\mathbf{k}}(\{R|\mathbf{t}\})=\mathbf{D}^{\mathbf{k}}(\{^dR|\mathbf{t}\}),\\
\overline{\mathbf{D}}^{\mathbf{k}}(\{R|\mathbf{t}\})=-\overline{\mathbf{D}}^{\mathbf{k}}(\{^dR|\mathbf{t}\}).\end{aligned}$$

As the wave vector is left invariant under the action of ${^d{\mbox{\textit{\textsf{1}}}}}$, the double little cogroup $^d\overline{{\mbox{$\mathcal{G}$}}}\,^{\mathbf{k}}$ of the wave vector $\mathbf{k}$ consists of the sets of elements of $\{R^{\mathbf{k}}\}$and $\{^dR^{\mathbf{k}}\}$ that correspond to the elements of $\overline{{\mbox{$\mathcal{G}$}}}\,^{\mathbf{k}}$, *cf.* eq. (\[def\_little\_cogroup\]): $$\label{def_double_little_cogroup}
^d\overline{{\mbox{$\mathcal{G}$}}}\,^{\mathbf{k}}=\{R\}\cup\{^dR\}.$$ The group $^d\overline{{\mbox{$\mathcal{G}$}}}\,^{\mathbf{k}}$ is a subgroup of the double point group $^d\overline{{\mbox{$\mathcal{G}$}}}$ of the space group $^d{\mbox{$\mathcal{G}$}}$. The index $\mid{^d\overline{{\mbox{$\mathcal{G}$}}}}/{^d\overline{{\mbox{$\mathcal{G}$}}}}\,^{\mathbf{k}}\mid$ equals the index $\overline{{\mbox{$\mathcal{G}$}}}\,^{\mathbf{k}}$ in $\overline{{\mbox{$\mathcal{G}$}}}$, and for the coset decomposition of $^d\overline{{\mbox{$\mathcal{G}$}}}$ with respect to $^d\overline{{\mbox{$\mathcal{G}$}}}\,^{\mathbf{k}}$ one can choose the same set of coset representatives as of the decomposition of $\overline{{\mbox{$\mathcal{G}$}}}$ with respect to $\overline{{\mbox{$\mathcal{G}$}}}\,^{\mathbf{k}}$. As a result, the star $*\mathbf{k}$ of the wave vector $\mathbf{k}$ in $^d{{\mbox{$\mathcal{G}$}}}$ coincides with $*\mathbf{k}$ in .

In analogy to the relationship between the little co-group $\overline{{\mbox{$\mathcal{G}$}}}\,^{\mathbf{k}}$ and the little group ${{\mbox{$\mathcal{G}$}}}\,^{\mathbf{k}}$ (*cf.* eq. (\[litgr\])), the little group $^d{\mbox{$\mathcal{G}$}}\,^{\mathbf{k}}$ can be defined as:

$$^d{\mbox{$\mathcal{G}$}}^{\mathbf{k}}=\{\{R^{\mathbf{k}}|\mathbf{v}^{\mathbf{k}}\} \in 
{\mbox{$\mathcal{G}$}}|R^{\mathbf{k}} \in \,^d\overline{{\mbox{$\mathcal{G}$}}}\,^{\mathbf{k}}\}.\label{litgr_double}$$

The allowed single-valued $\mathbf{D}^{\mathbf{k},i}$ and double-valued irreps $\overline{\mathbf{D}}^{\mathbf{k},j}$ of $^d{\mbox{$\mathcal{G}$}}\,^{\mathbf{k}}$ can be determined from the allowed single-valued $\mathbf{\Gamma}^{\mathbf{k}}$ and double-valued irreps $\overline{\mathbf{\Gamma}}^{\mathbf{k}}$ of $^d{\mbox{$\mathcal{T}$}}$ by stepwise induction along the composition series of $^d{\mbox{$\mathcal{G}$}}\,^{\mathbf{k}}$:

$$\label{comp_little_double}
^d{\mbox{$\mathcal{G}$}}^{\mathbf{k}}\rhd{^d{\mbox{$\mathcal{H}$}}}_1^{\mathbf{k}}\rhd\ldots\rhd{^d{\mbox{$\mathcal{H}$}}}_{m-1}^{\mathbf{k}}
\rhd{^d{\mbox{$\mathcal{H}$}}}_{m}^{\mathbf{k}} \rhd\ldots\rhd{^d{\mbox{$\mathcal{H}$}}}_n^{\mathbf{k}}   ={^d{\mbox{$\mathcal{T}$}}}$$

where $\mid{^d{\mbox{$\mathcal{H}$}}}_{m-1}^{\mathbf{k}}/{^d{\mbox{$\mathcal{H}$}}}_{m}^{\mathbf{k}}\mid=2\mbox{ or }3$.

The full single-valued irreps $\mathbf{D}^{\mathbf{*k},i}$ and full double-valued irreps $\overline{\mathbf{D}}^{\mathbf{*k},j}$ of the double space group $^d{\mbox{$\mathcal{G}$}}$ are induced from the allowed single-valued and double-valued irreps of the little group $^d{\mbox{$\mathcal{G}$}}^{\mathbf{k}}$. The full irreps are of dimension $r\times s$ where $s$ is the number of arms of the star $*\mathbf{k}$ and $r$ is the dimension of the corresponding allowed irrep of $^d{\mbox{$\mathcal{G}$}}^{\mathbf{k}}$. The matrices of the full irreps have a block structure with $s\times s$ blocks, each of dimension $r\times r$.

\[subsec:complexconjugation\]Complex conjugation
------------------------------------------------

For applications involving time-reversal symmetry, it is necessary to recall briefly the classification of irreps with respect to complex conjugation or, as it is commonly referred to, with respect to their *reality* (for more details, *cf.* ). Note that the concepts and results formulated in the following for ordinary representations of groups, are equally valid for single- and double-valued representations of double groups.

An irrep $\mathbf{D}$ of the group ${\mbox{$\mathcal{G}$}}$ is of:

-   the *first kind* (or *real*) if $\mathbf{D}$ is equivalent to a group of real matrices;

-   the *second kind* (or *pseudoreal*) if $\mathbf{D}$ is equivalent to $\mathbf{D}^*$ but not to any group of real matrices;

-   the *third kind* (or *complex*) if $\mathbf{D}$ is not equivalent to $\mathbf{D}^*$.

One can show that [@herring1937] $$\frac{1}{\left| \cal{G} \right|}\sum_{j=1}^{|{\mbox{$\mathcal{G}$}}|}\chi\left(g_j^2\right)=\left\{
\begin{array}{rl}
1&\textrm{iff $\mathbf{D}$ is of the first kind},\\
-1&\textrm{iff $\mathbf{D}$ is of the second kind},\\
0&\textrm{iff $\mathbf{D}$ is of the third kind},\\
\end{array}\right.
\label{test_general}$$ where $\chi$ is the character of the irrep $\mathbf{D}$ of . The obvious difficulties in the direct application of the *reality* test (eq. \[test\_general\]) to space-group irreps can be overcome if the sum over all operations of can be replaced by a sum over a relatively small number of space-group elements. This can be achieved using the relationship between the space-group irreps and the allowed irreps of the little group from which they are induced. As a result, the reality test for space group irreps can be written in the following form (see for details of the proof): the irrep $\mathbf{D}^{*\mathbf{k},j}$ of the space group induced from the allowed irrep $\mathbf{D}^{\mathbf{k},j}$ of the little group ${\mbox{$\mathcal{G}$}}^\mathbf{k}$, $\mathbf{D}^{*\mathbf{k},j}=\mathbf{D}^{\mathbf{k},j}\uparrow{\mbox{$\mathcal{G}$}}$, is of the first, second or third kind according to:

$$\frac{1}{|\mathbf{H}_{\mathbf{k}}|}\sum_{\{R_s|\mathbf{v_s}\}\in \mathbf{H}_{\mathbf{k}}}\chi^{\mathbf{k},j}\left(\{R_s|\mathbf{v_s}\}^2\right)=1, -1, \textrm{or} \ 0,$$

where $\chi^{\mathbf{k},j}$ is the character of the allowed irrep $\mathbf{D}^{\mathbf{k},j}$ of the little group ${\mbox{$\mathcal{G}$}}^\mathbf{k}$ and the sum is restricted to the set $\mathbf{H}_{\mathbf{k}}$ of coset representatives $\{R_s|\mathbf{v_s}\}$ of with respect to whose rotation parts transform $\mathbf{k}$ into a vector equivalent to $-\mathbf{k}$. Obviously, the element $\{R_s|\mathbf{v_s}\}^2$ leaves $\mathbf{k}$ invariant and therefore it belongs to ${\mbox{$\mathcal{G}$}}^\mathbf{k}$.

The irreps of the first and of the second kind are also known as *self-conjugate* while the irreps of the third kind form pairs of conjugated irreps $(\mathbf{D}, \mathbf{D}^*)$ which, in general, may be induced from allowed irreps of wave vectors belonging to different stars.

The concepts of *physically irreducible* representations or *time-reversal invariant* representations used in some physical applications are closely related to the above-discussed *reality* of the representations. Once the reality of a space-group irrep has been determined and, for the complex irreps, the pairs of conjugated irreps have been identified, the *time-reversal (TR)-invariant* irreps (single- and double-valued irreps of the double space groups) can be constructed according to:

-   If the irrep $\mathbf{D}$ is (a) single-valued and real, or (b) double-valued and pseudo-real, it is TR-invariant.

-   If the irrep $\mathbf{D}$ is (a) single-valued and pseudo-real or (b) double-valued and real, the TR-invariant irrep is the direct sum of $\mathbf{D}$ with itself. The dimension of the TR-invariant irrep doubles the dimension of $\mathbf{D}$. The label of the TR-invariant irrep consists of two copies of the label of $\mathbf{D}$.

-   If $\mathbf{D}_1$ and $\mathbf{D}_2$ form a pair of mutually conjugated irreps, the direct sum of both irreps is TR-invariant. The label of the TR-invariant representation is the union of the labels of the two irreps.

The program [Representations DSG]{}\[progDRepres\]
--------------------------------------------------

An algorithm based on the normal-subgroup procedure for the calculation of the irreps of the double space groups (Section \[irrep\_method\]) is implemented in the program [Representations DSG]{}.

As **Input**, the program needs the specification of the double space group $^d{\mbox{$\mathcal{G}$}}$ by the sequential [*IT*]{}A number of . (As already noted, the programs of BCS use the standard [[*IT*]{}A ]{}settings for the description of the space groups.) The [**k**]{}-vector data can be introduced by choosing the [**k**]{}-vector directly from a table provided by the program, where CDML labels are used to designate the symmetry [**k**]{}-vector types. The listed [**k**]{}-vector coefficients (called *conventional* **k**-vector coefficients) refer to a basis $\mathbf{a^*},\mathbf{b^*},\mathbf{c^*}$ which is dual to the conventional *IT*A settings of the space groups.

The **Output** of the program consists of two tables showing the matrices of the little-group irreps and the matrices of the irreps of the space group (*i.e* the full representations). In detail, the program produces:

-   The table of little-group representations (illustrated by the allowed irreps of the point $X:(0,1/2,0)$ of the double space group P$4/ncc$ (No. 130) shown in Fig. \[fig:representationsdsg\]).

    ![Figure 2: Screenshot of the output of the program [Representations DSG]{}, which shows the allowed irreducible representations of the little group of the point $X\,(0,1/2,0)$ in the double space group P$4/ncc$ (No. 130). Only part of the longer output has been included in the figure.[]{data-label="fig:representationsdsg"}](./Figure02.png){width="\textwidth"}

    The rows of the table are labeled by the symmetry operations of the little group while the columns are specified by the allowed little-group irreps. The symmetry operations are represented by the $(3\times 4)$ matrix-column pairs and $(2\times 2)$ $SU(2)$ matrices, and by the (modified) Seitz symbol. The first row corresponds to a general translation $\{1|\mathbf{t}\}$, while the subsequent rows show the data of the representatives $\{R|\mathbf{v}\}$ of the coset decomposition of the little group with respect to its translation subgroup. (From the fact that any element of the group $\{R|\mathbf{v}+\mathbf{t}\}$ can be expressed as a combination of a translation $\{1|\mathbf{t}\}$ and a coset representative $\{R|\mathbf{v}\}$ follows that the irrep matrices of $\{R|\mathbf{v}+\mathbf{t}\}$ are equal to products of the type $\mathbf{D}^{\mathbf{k}}(\{R|\mathbf{v}+\mathbf{t}\})=\mathbf{D}^{\mathbf{k}}(\{1|\mathbf{t}\})\mathbf{D}^{\mathbf{k}}(\{R|\mathbf{v}\})$.) The columns of single-valued irreps are followed by those of the double-valued irreps that are ordered according to their dimension. The labels of the single-valued irreps follow the notation of CDML: the labels consist of a $\mathbf{k}$-vector letter(s) and a sequential index. The labels of the double-valued irreps are constructed in a similar way: a bar over the $\mathbf{k}$-vector letter(s) distinguishes the double-valued irreps from the single-valued ones. In the example of the space group P$4/ncc$ (No. 130) and $\mathbf{k}$ vector $X$, $\mathbf{k}=(0,1/2,0)$ (Fig. \[fig:representationsdsg\]), the single-valued irreps have labels $X_1$ and $X_2$, and the two double-valued irreps are labelled $\overline{X}_3$ and $\overline{X}_4$. The irrep matrices are shown explicitly if the dimension of the irrep is $\leq 4$. When the dimension of the irrep is larger than 4, the output shows only the non-zero elements of the matrix, in the format: $(i;j):r$, signifying that the $(i,j)$ element of the matrix has the value $r$.

-   [The full-group irreps (illustrated by the full-group irreps of the double space group P$4/ncc$ at the point $X:(0,1/2,0)$ shown in Fig. \[fig:representationsdsg2\]).]{}

    ![Figure 3: Screenshot of the output of the program [Representations DSG]{}, which shows the irreducible representations of the double space group P$4/ncc$ (No. 130) at the point $X\,(0,1/2,0)$ (full-group). The integer in parenthesis after the irrep label indicates the *reality* of the representation; the corresponding pair of complex-conjugate irreps is also specified. Only part of the longer output has been included in the figure.[]{data-label="fig:representationsdsg2"}](./Figure03.png){width="\textwidth"}

    The indication of the arms of the star $*\mathbf{k}$ precedes the table of the matrices of the full-group irreps $\mathbf{D}^{*k,j}$ of $^d{\mbox{$\mathcal{G}$}}$ induced from the allowed irreps $\mathbf{D}^{\mathbf{k},j}$ of the little group $^d{\mbox{$\mathcal{G}$}}^\mathbf{k}$. The structure and the organization of the data of the full-group irrep table is very similar to that of the little-group irreps: the coset representatives of the decomposition $^d{\mbox{$\mathcal{G}$}}:{\mbox{$\mathcal{T}$}}$ of the double space group with respect to the translation subgroup specify the rows of the table, while the columns correspond to the full space-group irreps. The symmetry operations are described by their matrix-column pairs, $SU(2)$ matrices and Seitz symbols. The labels of the space-group irrep follow the labels of the allowed little-group irreps from which they are induced: *e.g.* as in the case of little-group irreps, a bar over the $\mathbf{k}$ vector letter indicates a double-valued irrep.

    The number in brackets after the irrep label specifies the reality of the irrep per eq. (\[test\_general\]): (1) indicates an irrep of the first kind, [*i.e.*]{} [*real*]{}; (-1) - an irrep of the second kind, or [*pseudoreal*]{}; and (0) - an irrep of the third kind, or [*complex*]{}. The pairs of the complex conjugated irreps are also indicated in the output. For example, the double-valued irreps $^*\overline{X}_3$ and $^*\overline{X}_4$ form a complex-conjugated pair as indicated immediately before the table of full irreps (Fig. \[fig:representationsdsg2\]).\
    As another example, the double-valued irrep $\overline{P}_7$ in the cubic double space group I$a\overline{3}$ (No. 206) is real, so that it doubles when time-reversal is considered. On the contrary, the double-valued irrep $\overline{P}_7$ in the cubic double space group I$4_132$ (No. 214) is pseudoreal and it does not double with time-reversal.

\[point\_irreps\]Representations of the double point groups
-----------------------------------------------------------

The crystallographic double point groups and their representations have been extensively discussed in the literature. Sets of irrep compilations can be found, for example, in , and . The data on the irreps of the crystallographic double point groups are now also online accessible via BCS. The irrep data was obtained by the program [Representations DSG]{} (*cf.* Section \[progDRepres\]) for the particular case of $\mathbf{k} = \Gamma (0,0,0)$.

For each of the 32 crystallographic double point groups, specified by their international (Hermann-Mauguin) and Schoenflies symbols, the retrieval tool [Representations DPG]{} (http://www.cryst.ehu.es/cryst/representationsDPG) displays the following set of tables:

1.  *Character table*

    As usual, the table entries correspond to characters of the irreps (rows) for the conjugacy classes of symmetry operations (columns) of the chosen double point group. The irreps are labelled according to the notation of , by the $\Gamma$ labels given by and by the labels generated by the program [Representations DSG]{}. A bar over the irrep label distinguishes the double-valued (known also as [*spinor*]{}) representations. The distribution into conjugacy classes of the symmetry operations of the double point group (designated by (modified) Seitz symbols), is also indicated.

2.  *Table of irrep matrices*

    The matrices of the irreps, as calculated by [Representations DSG]{} for the particular case of $\mathbf{k}=\Gamma\,(0,0,0)$ are also provided by [Representations DPG]{}. The matrices are explicitly listed for each symmetry operation (rows) and irrep (columns). The symmetry operations of the double point group are specified by the pair of $(3\times 3)$ rotation and $(2\times 2)$ $SU(2)$ matrices, and by the (modified) Seitz symbol. The classification of irreps with respect to complex conjugation is revealed by the number in brackets after the irrep label: as already noted, (1) indicates an irrep of the first kind, [*i.e.*]{} [*real*]{}; (-1) - an irrep of the second kind, or or [*pseudoreal*]{}; and (0) - an irrep of the third kind, or [*complex*]{}. The pairs of complex conjugated irreps are also listed.

Screenshots of the character table and the table of irrep matrices for the double point group ${\mathit{}mm2} (C_{2v})$ are shown in Fig. \[fig:character\_point\] and Fig. \[fig:matrices\_point\].

![Figure 4: Screenshot of the output of the program [Representations DPG]{} which shows the character table of the double point group ${\mathit{}mm2} (C_{2v})$. The irreps are labelled according to the notation of (1) , (2) and (3) by the labels generated by the program. The distribution into conjugacy classes of the symmetry operations is also indicated. Only the first part of the output has been included in the figure.[]{data-label="fig:character_point"}](./Figure04.png){width="85.00000%"}

![Figure 5: Screenshot of the output of the program [Representations DPG]{} which shows the irrep matrices of the double point group ${\mathit{}mm2} (C_{2v})$. The symmetry operations are specified by the matrix presentations and Seitz symbols. The integer in parenthesis after the irrep label indicates the *reality* of the representation. Only the last part of the output has been included in the figure.[]{data-label="fig:matrices_point"}](./Figure05.png){width="90.00000%"}

\[sec:compatibilityrelations\]Compatibility relations
=====================================================

The problem
-----------

The so-called *compatibility relations* have different applications in solid-state physics as, for example, in the analysis of the electronic band structures or phonon dispersion curves. The compatibility relations are essential in the study of connectivity of energy bands as we move in a continuous way from one **k**-vector point to a neighbouring one with a different symmetry, or in crystal-field splitting problems that arise when a high-symmetry crystal field is perturbed by a field of lower symmetry.

From a group-theoretical point of view, the compatibility relations correspond to the so-called *subduction* relationships between the little-group representations of different **k**-vectors of the same space group $\cal{G}$ whose little groups form a group-subgroup pair. Let $\cal{G}$ and $\cal{H}$ be two groups with group-subgroup relation $\cal{G} > \cal{H}$ with $ n = | \cal{G} | / | \cal{H} |$ being the index of $\cal{H}$ in $\cal{G}$. Consider an *irrep* $\mathbf{D}_{\beta}= \{\mathbf{D}_{\beta}(g), g\in\cal{G} \}$ of a group $\cal{G}$: The subduction of $\mathbf{D}_{\beta}$ to the subgroup $\cal{H}$ results in a representation of the subgroup, the so-called *subduced* representation, formed by the matrices of those elements of $\cal{G}$ that also belong to the subgroup $\cal{H}$, that is, $\mathbf{D}_{\beta}(h)$ with $h\in\cal{H}<\cal{G}$. This subduced representation $(\mathbf{D}_{\beta} \downarrow \cal{H})=\mathbf{D}^S$ is in general reducible, and is decomposable into irreps $\mathbf{d}_{\alpha}$ of $\cal{H}$:

$$(\mathbf{D}_\beta \downarrow {\cal{H}})(h) = \mathbf{D}^S(h) \sim \bigoplus_\alpha n_\alpha^{(\beta)}\mathbf{d}_\alpha(h),\qquad h\in\cal{H}.$$

The multiplicities $n_\alpha^{(\beta)}$ of the *irreps* $\mathbf{d}_{\alpha}$ of $\cal{H}$ in the subduced representation $\mathbf{D}_\beta \downarrow \cal{H}$ can be calculated by the *reduction formula* (known also as the Schur orthogonality relation or *‘magic’ formula*):

$$\label{eq:magic}
n_\alpha^{(\beta)} = \frac{1}{|\cal{H}|}\sum_h\chi^S(h)\chi_\alpha(h)^*,$$

where $\chi^S(h)$ is the character of the subduced representation and $\chi_\alpha(h)$ is the character of the irrep $\mathbf{d}_\alpha$ for the same element $h \in \cal{H}$.\
Consider two **k**-vectors **k**$_1$ and **k**$_2$ where **k**$_2$ = **k**$_1$ + $\boldsymbol{\kappa}$ with $\boldsymbol{\kappa}$ - an infinitesimal **k**-vector; *e.g.* **k**$_1$ could be the wave vector of a symmetry point (line) ’sitting’ on a symmetry line (plane) **k**$_2$. The little group ${\cal{G}}^{\textbf{k}_2}$ is in general a subgroup of the little group ${\cal{G}}^{\textbf{k}_1}$ and the compatibility relations between the *irreps* at **k**$_1$ and **k**$_2$ (in the limit $\kappa \rightarrow 0$) can be established by the subduction of the representations **D**$^{\textbf{k}_1,i}$ of the little group ${\cal{G}}^{\textbf{k}_1}$ onto the little subgroup ${\cal{G}}^{\textbf{k}_2}$. In the following we say that these two $\mathbf{k}$-vectors are *connected*.

Tables of compatibility relations of double space groups can be found in Miller & Love (1967) but the compiled tables provide only the relationships between the irreps of **k**-vector points and lines; neither the relations between the **k**-vector lines and planes nor the relations between the **k**-vector points and planes are made available. The program [DCOMPREL]{} (http://www.cryst.ehu.es/cryst/dcomprel) calculates the compatibility relations between the little-group irreps of the double space groups for any pair of symmetry-related wave vectors. Compatibility relations between little-group irreps of ordinary space groups can be calculated by the program [COMPATIBILITY RELATIONS]{}, also accessible in BCS (http://www.cryst.ehu.es/rep/rep$\_$correlation$\_$relations).

The method
----------

Essentially, the same algorithm is implemented in the programs [COMPATIBILITY RELATIONS]{} and [DCOMPREL]{}: it follows closely the subduction procedure explained above applied to the special case of little-group representations. Consider two symmetry-related wave vectors **k**$_1$ and **k**$_2$, where **k**$_2$ = **k**$_1$ + $\boldsymbol{\kappa}$ with ${\cal{G}}^{\textbf{k}_2}<{\cal{G}}^{\textbf{k}_1}$. The matrices of an irrep of the little group ${\cal{G}}^{\textbf{k}_1}$ associated to the symmetry operations that belong to ${\cal{G}}^{\textbf{k}_2}$ form a representation of the little group of $\mathbf{k}_2$, which in general, is reducible. The compatibility relations are extracted from the decomposition of the subduced representation ${\cal{G}}^{\textbf{k}_1} \downarrow {\cal{G}}^{\textbf{k}_2}$ into irreps of ${\cal{G}}^{\textbf{k}_2}$:

$$\mathbf{D}^{\textbf{k}_1,i} \downarrow {\cal{G}}^{\textbf{k}_2}\sim \displaystyle\bigoplus_{j=1}^s m_{\textbf{k}_2,j} \mathbf{D}^{\textbf{k}_2,j},$$

where $m_{\textbf{k}_2,j}$ represents the multiplicity of the irrep $\mathbf{D}^{\textbf{k}_2,j}$ in the subduced representation. Given the characters $\chi_{\mathbf{D}^{\textbf{k}_1,i}}$ of $\mathbf{D}^{\textbf{k}_1,i}\downarrow {\cal{G}}^{\textbf{k}_2}$ and $\chi_{\mathbf{D}^{\textbf{k}_2,j}}$ of $\mathbf{D}^{\textbf{k}_2,j}$ for all $g\in{\cal{G}}^{\textbf{k}_2}$ (*e.g.* calculated by the program [Representations DSG]{}), the multiplicities can be obtained by a variation of the reduction formula:

$$\label{magic_inf}
m_{\textbf{k}_2,j} = \frac{1}{|{\cal{G}}^{\textbf{k}_2}|}  \displaystyle\sum_{g} \chi_{\mathbf{D}^{\textbf{k}_2,j}}(g) \chi_{\mathbf{D}^{\textbf{k}_1,i}}(g)^*,   \quad \quad \forall g \in {\cal{G}}^{\textbf{k}_2}.$$

The expression (\[magic\_inf\]) is not convenient for practical calculations due to the infinite order of the little group ${\cal{G}}^{\textbf{k}_2}$. However, the summation in eq. (\[magic\_inf\]) can be first performed over all translations in ${\cal{G}}^{\textbf{k}_2}$ which will reduce the infinite sum to a finite one over the set $\mathbf{G}^{\textbf{k}_2}_o$ of the representatives $g_k$ of the decomposition of ${\cal{G}}^{\textbf{k}_2}$ in cosets with respect to its translation subgroup:

$$m_{\textbf{k}_2,j} = \frac{1}{n}  \displaystyle\sum_{k=1}^n \chi_{\mathbf{D}^{\textbf{k}_2,j}}(g_k) \chi_{\mathbf{D}^{\textbf{k}_1,i}}(g_k)^*;   \quad \quad \forall g_k \in {\mathbf{G}}^{\textbf{k}_2}_o.$$

The number *n* of coset representatives $g_k$ is equal to the order of the little co-group $\overline{\cal{G}}^{\textbf{k}_2}$.

The program [DCOMPREL]{}
------------------------

The procedure described above for the calculation of the compatibility relations between the little-group irreps of space groups is implemented in the program `DCOMPREL`. The program is available on BCS, and for a number of cases the results have been successfully checked against the compatibility-relation data listed in .

#### Input:

The specification of the double space group by the [[*IT*]{}A ]{}number leads to a menu with the list of the different symmetry $\mathbf{k}$-vectors for the group chosen. The choice of a $\mathbf{k}$-vector, $\mathbf{k}_1$, produces a second output with the list of all $\mathbf{k}$-vectors that can be connected to it. The user can then ask for the compatibility relations between the chosen $\mathbf{k}_1$-vector and a single $\mathbf{k}$-vector or between $\mathbf{k}_1$ and all the $\mathbf{k}$-vectors in the list.

#### Output:

After a summary of the input data, a table with the compatibility relations between the little-group *irreps* of the selected **k**-vectors is shown.

Fig. \[fig:compatibility\] shows a partial output of the program [DCOMPREL]{} of the compatibility relations of the $\mathbf{k}$-vector line $W(0,1/2,w)$ with all symmetry-related wave vectors in the double space group ${\mathit{P}4/ncc}$. Note that in some cases, $W$ is the high-symmetry point of the pair (as in $W\to B,F,GP$) while in others $W$ is the point of lower symmetry (as in $R\to W$). The number in parenthesis indicates the dimension of the representation.

![Figure 6: Screenshot of the output of the program [DCOMPREL]{} which shows all the compatibility relations which involve the irreps of the $\mathbf{k}$-vector $W(0,1/2,w)$ in the double space group ${\mathit{P}4/ncc}$ (No. 130). The first column gives the list of $\mathbf{k}$-vectors symmetry-related to $W$, and the second column shows the corresponding compatibility relations. In some cases, $W$ is the high-symmetry point of the pair (as in $W\to B,F,GP$) while in others $W$ is the point of lower symmetry (as in $R\to W$). Only part of the longer output has been included in the figure.[]{data-label="fig:compatibility"}](./Figure06.png){width="55.00000%"}

Example: electronic bands of germanium
--------------------------------------

As an application of the use of the representations of the double space groups and the compatibility relations we consider the electronic bands of germanium, with space group F$d\bar{3}m$ (No. 227) and occupied atomic position $8a$ $(1/8,1/8,1/8)$. Figure \[fig:dresselhaus\] reproduces the figures 12.10 and 14.41 in . The labels of the irreducible representations have been changed to the labels used in [REPRESENTATIONS DSG]{}. Figure \[fig:dresselhaus\] (a) shows the band structure when the spin-orbit coupling is not considered and the relevant irreps are the single-valued irreducible representations, and Figure \[fig:dresselhaus\] (b) reproduces the band structure when the spin-orbit coupling is considered and the relevant irreps are the double-valued irreducible representations. Figure \[fig:dresselhauscomp\] shows the compatibility relations between the $\Gamma(0,0,0)$ point and the $\Delta(0,v,0)$ line and between the $X(0,1/2,0)$ point and the $\Delta$ line in the double space group F$d\bar{3}m$. It can be checked that these compatibility relations agree with the paths $\Gamma \leftrightarrow\Delta\leftrightarrow X$ in Figures \[fig:dresselhaus\](a) and \[fig:dresselhaus\](b). In some cases, these relations could be useful to identify the irreducible representations at the $\vec{k}$ points of high symmetry. The compatibility relations give the connectivities between the branches of the electronic band along the Brillouin zone and the dimensions of the irreps at each point give the degeneracies of the energy.

![Figure 7: Electronic bands of germanium (a) without and (b) with spin-orbit coupling (adapted from Dresselhaus et al. (2008). The labels of the electronic bands correspond to the irreps labels used by [REPRESENTATIONS DSG]{}.[]{data-label="fig:dresselhaus"}](./Figure07.png){width="\textwidth"}

![Figure 8: Compatibility relations between the point $\Gamma$ and the line $\Delta$ (left) and between the point $X$ and the line $\Delta$ (right) of the space group F$d\bar{3}m$, given by the program [DCOMPREL]{}.[]{data-label="fig:dresselhauscomp"}](./Figure08.png){width="60.00000%"}

\[sec:sitesymmetry\]Site-symmetry approach
==========================================

The problem
-----------

The symmetry of the *extended states* (phonons, soft modes, electronic bands, *etc.*) of crystal structures (described by the irreps of the space group of the crystal) over the entire Brillouin zone is related to the symmetry of *localized states* (local atomic displacements, atomic orbitals, *etc.*) of the constituent structural units (described by the irreps of the local symmetry groups). A procedure for the determination of such a relationship is very useful, as it allows the prediction of the symmetry of the possible extended states starting from the crystal structural data. In particular, it is useful in distinguishing topological materials [@NaturePaper]. Formally, the procedure relating localized and extended crystalline states can be described by induction of a representation of a space group $\cal{G}$ from a finite subgroup $\cal{H}$, followed by a reduction into irreps. In other words, the induction method permits the calculation of the symmetry of the compatible extended states transforming according to irreps of crystal space group $\cal{G}$ induced by a localized state described by an irrep of the local or site symmetry group $\cal{H}=\cal{S}$. However, the calculation of the space-group irreps induced by the irreps of the site-symmetry group is not an easy task: the fact that the site-symmetry group $\cal{S}$ (isomorphic to a point group) is a subgroup of infinite index of $\cal{G}$ implies that the representation of $\cal{G}$ induced by an irrep of $\cal{S}$ must be of infinite dimension, and therefore difficult to calculate directly. The so-called *site-symmetry approach* resolves this problem by applying the *Frobenius reciprocity theorem*, which states that the multiplicities of the irreps of a group $\cal{G}$ in the induced representation from an irrep of a subgroup $\cal{H}$ of $\cal{G}$ can be determined from the multiplicities of the irreps of in the representations subduced from $\cal{G}$ to . (For a detailed presentation of the method, its discussion and applications, *cf.* ). In this way, the knowledge of the subduced representation of the irreps of the space group onto the site-symmetry group, which are relatively easy to compute once the space group irreps are known, is enough to obtain all the necessary information about the representations induced by the irreps of the site-symmetry group into the space group.

The site-symmetry method is implemented in two programs of BCS: the program [SITESYM]{} for ordinary representations of space groups, and the program [DSITESYM]{} for the double space groups. The following explanations of the site-symmetry method are equally valid for the cases of space groups and double space groups.

The method
----------

The objective of the site-symmetry programs is to find the symmetry of the crystal extended states induced by localized states of some of the constituent structural units. Examples of the application of this method can be found in the supplementary material of . In group-theoretical terms, this task requires the derivation of the irreps of a space group $\cal{G}$ at any point in the reciprocal space (which classify the extended states of the structure) induced by the irreps of the site-symmetry group of a Wyckoff position (according to which localized states are classified). Two basic concepts of representation theory, namely *subduction* and *induction* are essential for the site-symmetry approach. As already explained, the *subduction coefficients* specify the decomposition of a representation of a group $\cal{G}$ into irreps of one of its subgroups $\cal{H}< \cal{G}$ while the induction procedure permits the construction of a representation of $\cal{G}$ starting from a representation of $\cal{H}$. Consider the decomposition of $\cal{G}$ in cosets with respect to $\cal{H}$ with coset representatives $\{{\mbox{\textit{\textsf{q}}}}_m, m=1,\dots,n\}$, *cf.* eq. (\[coset\]). If $\mathbf{d}_\alpha=\{\mathbf{d}_\alpha(h),h\in\cal{H}\}$ is an irrep of $\cal{H} < \cal{G}$, then the matrices of the induced representation $(\mathbf{d}_\alpha \uparrow \cal{G}) = \mathbf{D}^I$ of $\cal{G}$ are constructed as follows:

$$\label{ind-1}
\mathbf{D}^I(g)_{kt,js} =\left\{
\begin{array}{l l}
\mathbf{d}_\alpha(q_k^{-1} g q_j)_{ts} & \text{if }  q_k^{-1} g q_j \in \cal{H}\\
\rule{0pt}{4ex}0 & \text{if } q_k^{-1} g q_j \not \in \cal{H}
\end{array}\right.,$$

where $k,j = 1,\dots,n$ and $t,s=1,\dots,m$ with $m$ the dimension of the irrep $\mathbf{d}_\alpha$ of $\cal{H}$. Its characters are given by: $$\chi^I_\gamma(g) =\sum_j \chi^\alpha(q_j^{-1}g q_j),$$

where $\chi^\alpha(q_j^{-1}g q_j)$ is the trace of the $j$-th diagonal block of $\mathbf{D}^I(g)$. In general, the induced representations are reducible, and as such it is possible to decompose them into irreps $\mathbf{D}_\beta$ of $\cal{G}$: $$\label{ind-dec}
(\mathbf{d}_\alpha \uparrow{\cal{G}}) = \mathbf{D}^I \sim \bigoplus_\beta n^{(\alpha)}_\beta \mathbf{D}_\beta,$$ where $n^{(\alpha)}_\beta$ are the multiplicities of the irreps of $\cal{G}$ in the induced representation, and it is possible to calculate them applying the reduction formula (\[eq:magic\]):

$$\label{ind-multi}
n^{(\alpha)}_\beta = \frac{1}{|\cal{G}|}\sum_g\chi^I(g)\chi^*_\beta(g),$$

where $\chi_\beta(g)$ is the character of the element $g$ of the irrep $\mathbf{D}_\beta$ of $\cal{G}$.

The dimension of the induced representation can be directly read off the equation for its construction, eq.(\[ind-1\]): $$\label{eq:ind-dim}
dim(\mathbf{d}_\alpha \uparrow {\cal{G}})= dim (\mathbf{d}_\alpha) \frac{|\cal{G}|}{|\cal{H}|}.$$

This result points to the difficulties for the direct calculation of a representation of a space group $\cal{G}$ induced from an irrep of a finite subgroup $\cal{H}$ of $\cal{G}$: due to the infinite index of $\cal{H}$ in $\cal{G}$, the dimension of the representation of $\cal{G}$ induced from an irrep of $\cal{H}$ is infinite. By means of the site-symmetry approach it is possible to determine the multiplicities $n^{(\alpha)}_\beta$ of an irrep of $\cal{G}$ in the induced representation without the necessity of constructing the infinite-dimensional representation. As stated above, the method is based on the Frobenius reciprocity theorem, according to which the multiplicity of an irreducible irrep $\mathbf{D}_\beta(g)$ of $\cal{G}$ in a representation $\mathbf{d}_\alpha \uparrow\cal{G}$ of $\cal{G}$ induced by an irrep $\mathbf{d}_\alpha$ of $\cal{H}<\cal{G}$ is equal to the multiplicity of the irrep $\mathbf{d}_\alpha$ of $\cal{H}$ in the representation $\mathbf{D}_\beta \downarrow \cal{H}$ subduced by $\mathbf{D}_\beta$ of $\cal{G}$ to $\cal{H}$, *i.e.* $n_\alpha^{(\beta)}=n_\beta^{(\alpha)}$. In other words, in order to calculate the frequencies $n_\beta^{(\alpha)}$, eq.(\[ind-dec\]), it is sufficient to calculate the multiplicities $n_\alpha^{(\beta)}$ of $\mathbf{d}_\alpha$ in the subduced representation $\mathbf{D}_\beta \downarrow \cal{H}\sim\bigoplus_\alpha$ $n_{\alpha}^{(\beta)}d_{\alpha}$.

As already discussed, the irreps of a space group $\cal{G}$ are classified according to the reciprocal space wave vectors and for each vector $\bf{k}$ there is a finite set of irreps of $\cal{G}$. The main idea of the site-symmetry approach is that although the representation induced by an irrep of a site-symmetry group in a space group $\cal{G}$ has infinite dimension, it is possible to know the part of this induced representation that corresponds to any set (finite) of irreps of $\cal{G}$ for a given wave vector $\bf{k}$ (calculated in advance). As in most applications only the irreps of $\cal{G}$ related to a few wave vectors are of interest, this partial knowledge of the induced representation proves to be sufficient.

The algorithmic procedure, based on the site-symmetry approach and implemented in the site-symmetry programs [SITESYM]{} and [DSITESYM]{}, is the following:

1.  Consider an occupied Wyckoff position and determine its site-symmetry group $\cal{S}$. Note that it is not sufficient to determine the point group isomorphic to the site-symmetry group: it is necessary to obtain the set of space-group symmetry operations, with their rotation and translations parts, that belong to the site-symmetry group, since the representation matrices of symmetry operations related by translations are, in general, different.

2.  Calculate the irreps of the space group $\cal{G}$ for the wave vectors $\bf{k}$ of interest. The BCS program [REPRES]{} for ordinary space-group representations (or [Representations DSG]{} for single- and double-valued irreps of double space groups) provides the irrep matrices for any element of the space group (or of the double space group).

3.  From the obtained space-group irreps, calculate the representations subduced to the site-symmetry group $\cal{S}$ obtained in the first step.

4.  From the irreps of the site-symmetry group (tabulated by the BCS programs [POINT]{} [@aroyo2006] for ordinary point-group representations or [Representations DPG]{} for single- and double-valued irreps of double point groups), and making use of the reduction formula (\[eq:magic\]), calculate the multiplicities of the site-symmetry irreps in the subduced representations.

5.  Apply the Frobenius reciprocity theorem to obtain the multiplicities of the irreps in the induced representation of $\cal{G}$ from the multiplicities of irreps of $\cal{S}$.

The programs [SITESYM]{} and [DSITESYM]{}
-----------------------------------------

This algorithm for the calculation of the multiplicities of space-group irreps in the representations induced by the irreps of a site-symmetry group ($\cal{S}<\cal{G}$) has been implemented in the site-symmetry programs of BCS. The description that follows refers to the program [DSITESYM]{} (http://www.cryst.ehu.es/cryst/dsitesym) but similar explanations are also applicable to the program [SITESYM]{} (http://www.cryst.ehu.es/cryst/sitesym).

The site symmetry approach allows the determination of the symmetry relationships between the extended (Bloch) and localised (Wannier-type) electronic states in crystals. As an example of the utility of [DSITESYM]{} we will consider the determination of the symmetry of electronic states in (GaAs)$_m$(AlAs)$_n$ semiconductor superlattices grown along the \[001\] direction, a problem discussed in detail in . The output of the program will be illustrated by the specific calculations of the symmetry relationships between the electronic band states at the point $\vec{k}=M(1/2,1/2,0)$ of the double space group P$\overline{4}m2$ (No. 115) (the structure-symmetry group of (GaAs)$_m$(AlAs)$_n$) and the localised states of atomic orbitals at the Wyckoff position $2g$ $(0,1/2,z)$. In other words, the multiplicities of the irreps of the double space group P$\overline{4}m2$ at the point $\vec{k}=M(1/2,1/2,0)$ shown in the output of [DSITESYM]{}, describe the transformation properties of the extended electronic states induced by the irreps of the double site-symmetry group of the Wyckoff position $2g$ $(0,1/2,z)$ according to which the localised one-electron wave functions transform.

#### Input:

The required input data include the specification of the double space group, of the occupied Wyckoff position and of the wave-vector label of the space-group irreps whose induced-representation multiplicities are to be calculated. The information is entered using three forms: in the first one, the double space group is specified by its [[*IT*]{}A ]{}sequential number; the wave vector of interest and the occupied Wyckoff position can be selected from the corresponding lists produced by the program.

#### Output:

After a header that reproduces the input data, the program displays the following tables with results of the intermediate steps of the procedure (the induction table, that lists the final results of the site-symmetry calculations, is at the bottom of the screen (v)):

1.  *List of operations of the double site-symmetry group $\cal{S}$*: each of the symmetry operations of the double space group that leaves the Wyckoff position representative point invariant, is represented in a $x,y,z$ and in a matrix notation. Labels ($g_1,...,g_n$) necessary for later referencing are assigned to each element of $\cal{S}$.

    The double site-symmetry group of the position 2$g \ (0,\frac12,z)$ of the double space group ${\mathit{P}\bar 4m2}$ is formed by 8 symmetry operations as shown in the screenshot (Fig. \[fig:symmetryoperations\]) of the program [DSITESYM]{}.

    ![Figure 9: Screenshot of the partial output of the program [DSITESYM]{} which shows the 8 symmetry operations of the double site-symmetry group of the Wyckoff position 2$g \ (0,\frac12,z)$ of the double space group ${\mathit{P}\bar 4m2}$ (No. 115). The symmetry operations are specified by the matrix presentations and shorthand notation.[]{data-label="fig:symmetryoperations"}](./Figure09.png){width="70.00000%"}

2.  *Character table of the double point group*. The table reproduces the character table of the irreps of the double point group isomorphic to the site-symmetry group, as provided by the program [Representations DPG]{} of BCS. The labels of the irreps are given in notations by and as generated by the program. The site symmetry group of $(0,\frac12,z)$ is isomorphic to the double point group $mm2$ and has five irreps. Its character table is shown in the Fig. \[fig:charactertable\].

    ![Figure 10: Screenshot of the partial output of the program [DSITESYM]{} which shows the chatacter table of the double point group $mm2$, isomorphic to the site-symmetry group of the Wyckoff position $(0,1/2,z)$ of the space group ${\mathit{P}\bar 4m2}$ (No. 115). The irreps are labelled according to the notation of (1) and (2) the labels generated by the program.[]{data-label="fig:charactertable"}](./Figure10.png){width="\textwidth"}

3.  *Table of characters of the subduced representations.* The program [Representations DSG]{} calculates the characters of the elements of the site-symmetry group (obtained in the first step) for each of the irreps $\mathbf{D}_\beta$ of ${^d\cal{G}}$ of the selected wave vector. In this way, we obtain the characters of the subduced representations $\mathbf{D}_\beta \downarrow \cal{S}$ of $\cal{S}$.

    The double space group ${\mathit{P}\bar 4m2}$ has seven irreps for the wave-vector $M$ which subduce seven representations of the double site-symmetry group $\cal{S}$=$mm2$ with the characters shown in Fig. \[fig:characterssubducedirreps\].

    ![Figure 11: Screenshot of the partial output of the program [DSITESYM]{} which shows the characters of the subduced representations of the site-symmetry group $\cal{S}$=$mm2$ by the irreps of the double space group ${\mathit{P}\bar 4m2}$ (No. 115) at the $\mathbf{k}$-point $M$. The columns are specified by the symmetry operations of the site-symmetry group, *cf.* Fig.\[fig:symmetryoperations\]. []{data-label="fig:characterssubducedirreps"}](./Figure11.png){width="75.00000%"}

4.  *Table of the decompositions of the subduced representations*. The multiplicities $n_\alpha^{(\beta)}$ of the double point-group irreps $\mathbf{d}_{\alpha}$ of $\cal{S}$ in the subduced representations $\mathbf{D}_{\beta} \downarrow \cal{S}$ are obtained by the application of the reduction formula, *cf.* eq. (\[eq:magic\]). In the example, the decompositions of the representations $^{*}M_i \downarrow mm2$, $i=1,\ldots,7$, into irreps of $mm2$ are shown in Fig. \[fig:decomposition\]:

    ![Figure 12: Screenshot of the partial output of the program [DSITESYM]{} which shows the decompositions of the subduced representations $M_i \downarrow mm2$ (see Fig. \[fig:characterssubducedirreps\]) into irreps of $mm2$, *cf.* Fig. \[fig:charactertable\].[]{data-label="fig:decomposition"}](./Figure12.png){width="55.00000%"}

5.  *Table of induced representations.* According to Frobenius reciprocity theorem, the multiplicities $n_\beta^{(\alpha)}$ of the irreps of ${^d\cal{G}}$ for a given **k**-vector in the representations $\mathbf{d}_\alpha \uparrow{^d\cal{G}}$ (induced from the irreps $\mathbf{d}_\alpha$ of the site-symmetry group $\cal{S}$) are obtained by “transposing” the table of the decompositions of the subduced representations $\mathbf{D}_{\beta} \downarrow \cal{S}$.

    The table of representations of the double space group ${\mathit{P}\bar 4m2}$ at the point $M$ induced by the irreps of the site-symmetry group $mm2$ of the Wyckoff position 2$g$ is shown in Fig. \[fig:representationsbar4m2\]. The rows of the table correspond to the irreps $\mathbf{d}_\alpha$ of the site-symmetry group $mm2$ (cf. Fig. \[fig:charactertable\]); the entries in each row indicate the multiplicities of the $M$-irreps of ${\mathit{P}\bar 4m2}$ in the (infinite) induced representation $\mathbf{d}_\alpha\uparrow {\mathit{P}\bar 4m2}$:

    ![Figure 13: Screenshot of the partial output of the program [DSITESYM]{} which shows the representations of the double space group ${\mathit{P}\bar 4m2}$ (No. 115) at the point $M$ induced by the irreps of the site-symmetry group $mm2$ of the Wyckoff position 2$g$.[]{data-label="fig:representationsbar4m2"}](./Figure13.png){width="55.00000%"}

      -------------------------------------------------------------------------------------------------------
      $A_1 \uparrow {\mathit{P}\bar 4m2} \sim {{M_5}} \oplus \cdots$
      $A_2 \uparrow {\mathit{P}\bar 4m2} \sim {M_5}\oplus \cdots$
      $B_2 \uparrow {\mathit{P}\bar 4m2} \sim {M_1}\oplus{M_2}\oplus \cdots$
      $B_1 \uparrow {\mathit{P}\bar 4m2} \sim {M_3}\oplus{M_4}\oplus \cdots$
      $\overline{E} \uparrow {\mathit{P}\bar 4m2} \sim {\overline{M}_6}\oplus{\overline{M}_7}\oplus \cdots$
      -------------------------------------------------------------------------------------------------------

    The obtained results coincide exactly with the corresponding data of Table 2 of ref. .

\[sec:bandrepresentations\]Band representations and Elementary band representations
===================================================================================

\[bandrepproblem\]The problem
-----------------------------

The concept of a band representation (BR) was introduced by , as a set of extended energy states over the entire reciprocal space, $E_n(\mathbf{k})$, related to the symmetry of (exponentially) localized states (Wannier orbitals). The basic structure description of a crystal includes the assignation of a space group and the determination of the lattice parameters and of the atomic coordinates (occupied Wyckoff positions) of a minimal set of atoms that belong to the asymmetric unit (known also as the set of independent atoms). The atomic positions of all the atoms in the crystal are obtained through the action of the symmetry operations of the space group onto the coordinates of the independent atoms. As explained in Section \[sec:sitesymmetry\], the finite set of symmetry operations of the space group that keeps a point of a Wyckoff position invariant is its site-symmetry group ${\mbox{$\mathcal{S}$}}$, and it is isomorphic to a point group. It is important to note that the site-symmetry groups ${\mbox{$\mathcal{S}$}}_q$ of all points $q$ belonging to a Wyckoff position $\mathbf{Q}=\{q\}$ are isomorphic to the same point-group type, and in that sense one speaks of a site-symmetry group of a Wyckoff position ${\mbox{$\mathcal{S}$}}_{\mathbf{Q}}$. The localized states (the atomic orbitals, for example) of an atom that occupy a given Wyckoff position transform according to a representation $\mathbf{d}_{\alpha}$ of its site-symmetry group. When the spin-orbit coupling is also considered, the localized states transform according to a representation of the double site-symmetry group, isomorphic to a double point group. A *band representation* can be defined as the induced representation $\mathbf{d}_{\alpha}\uparrow\mathcal{G}$, being a particular case of the general description of the site-symmetry approach discussed in Section \[sec:sitesymmetry\] If the representation $\mathbf{d}_{\alpha}$ of the site-symmetry group ${\mbox{$\mathcal{S}$}}$ is reducible, then it decomposes into irreps of the site-symmetry group, $$\label{ind-decband}
\mathbf{d}_\alpha= \bigoplus_\beta n_\beta^{(\alpha)} \mathbf{d}_\beta,$$ and the induced BR decomposes into BRs induced from the irreducible representations $\mathbf{d}_\beta$ $$\label{ind-decbandirrep}
\mathbf{d}_\alpha\uparrow\mathcal{G}= \bigoplus_\beta  n_\beta^{(\alpha)}(\mathbf{d}_\beta\uparrow\mathcal{G}).$$ Therefore, in order to have a complete information of the BRs induced from a given Wyckoff position, it is necessary to calculate only the BRs induced from the *irreps* of its site-symmetry group.

The BRs induced from different irreps of the same site-symmetry group are not equivalent, but BRs induced from different Wyckoff positions could be equivalent. The definition of equivalence of BRs is different from the definition of equivalence of representations, where we say that two representations are equivalent if all the multiplicities of their decomposition into irreps are the same. However, for certain applications, as for example, the study of topological phases, we need a stronger form of equivalence [@michel1992; @zeiner2000; @NaturePaper; @EBRtheory].

**Definition 1:** Two band representations $\rho_{{\mbox{$\mathcal{G}$}}}^{\mathbf{k}}$ and $\sigma_{{\mbox{$\mathcal{G}$}}}^{\mathbf{k}}$ are equivalent iff there exists a unitary matrix-valued function $S(\mathbf{k},s,{\mbox{\textit{\textsf{g}}}})$ smooth in $\mathbf{k}$ and continuous in $s$ such that for all ${\mbox{\textit{\textsf{g}}}}\in\mathcal{G}$

-   $S(\mathbf{k},s,{\mbox{\textit{\textsf{g}}}})$ is a band representation for all $s\in[0,1]$,

-   $S(\mathbf{k},0,{\mbox{\textit{\textsf{g}}}})=\rho_{{\mbox{$\mathcal{G}$}}}^{\mathbf{k}}({\mbox{\textit{\textsf{g}}}})$ and

-   $S(\mathbf{k},1,{\mbox{\textit{\textsf{g}}}})=\sigma_{{\mbox{$\mathcal{G}$}}}^{\mathbf{k}}({\mbox{\textit{\textsf{g}}}})$

This definition implies that the BRs $\rho_{{\mbox{$\mathcal{G}$}}}^{\mathbf{k}}({\mbox{\textit{\textsf{g}}}})$ and $\sigma_{{\mbox{$\mathcal{G}$}}}^{\mathbf{k}}({\mbox{\textit{\textsf{g}}}})$ subduce into the same little group representations at all points in the Brillouin zone. This definition also implies the following: consider two Wyckoff positions $\mathbf{Q}$ and $\mathbf{Q}\,'$ and let $\rho_{{\mbox{$\mathcal{G}$}}}^{\mathbf{k}}({\mbox{\textit{\textsf{g}}}})$ and $\sigma_{{\mbox{$\mathcal{G}$}}}^{\mathbf{k}}({\mbox{\textit{\textsf{g}}}})$ be two BRs induced from an irrep $\rho$ of the site-symmetry group $\mathcal{S}_{\mathbf{Q}}$ of $\mathbf{Q}$ and an irrep $\sigma$ of the site-symmetry group $\mathcal{S}_{\mathbf{Q}\,'}$ of $\mathbf{Q}\,'$. The intersection $\mathcal{S}_{\mathbf{Q}_0}=\mathcal{S}_{\mathbf{Q}}\cap\mathcal{S}_{\mathbf{Q}\,'}$ of the two site-symmetry groups is the site-symmetry group of another Wyckoff position $\mathbf{Q}_0$ that can be identified relatively easy: some of the point coordinates of $\mathbf{Q}_0$ should be represented by variable parameters that interpolate between the point coordinates of the Wyckoff positions $\mathbf{Q}$ and $\mathbf{Q}\,'$. If for a given irrep $\rho_0$ of $\mathcal{S}_{\mathbf{Q}_0}$, the induced representations into the site-symmetry groups of $\mathbf{Q}$ and $\mathbf{Q}\,'$ satisfy, $\rho_0\uparrow S_{\mathbf{Q}}=\rho$ and $\rho_0\uparrow S_{\mathbf{Q}\,'}=\sigma$,then the two BRs $(\rho\uparrow G)\downarrow G^{\mathbf{k}}=\rho_G^{\mathbf{k}}$ and $(\sigma\uparrow G)\downarrow G^{\mathbf{k}}=\sigma_G^{\mathbf{k}}$ are equivalent.

Once a criterion of equivalence of BR is established, we define:

**Definition 2:** A band representation is called *composite* if it is equivalent to the direct sum of other band representations. A band representation that is not composite is called *elementary*.

All the elementary band representations (EBR) are induced from the so-called *Wyckoff positions of maximal symmetry* (*cf.* Section \[wyckoffmaximal\] for the definition of the concept) but the opposite is not true. This fact was already pointed out for single-valued BRs in , where also the complete list of exceptions was given. The list of exceptions for BRs induced from double-valued irreducible representations of the double space groups can be found in the supplementary material of .

The generalization of the procedure for the calculation of EBRs for systems with time-reversal symmetry is straightforward. In this case, it is enough to consider the BRs induced from the physically irreducible representations of the site-symmetry groups of the Wyckoff positions of maximal symmetry. The physically irreducible or TR-invariant representations of the site-symmetry groups can be constructed following the steps explained in Section \[subsec:complexconjugation\]. The TR-invariant irreps induce TR-invariant BRs in $\mathcal{G}$, and using the definitions 1 and 2 above the TR-invariant EBRs can be calculated. Not all the single-valued TR-invariant BRs induced from the double site-symmetry groups of Wyckoff positions are elementary: the list of exceptions is given in [@NaturePaper]. However, all the double-valued TR-invariant BRs induced from the double site-symmetry groups of Wyckoff positions of maximal symmetry are elementary, as a consequence of Kramer’s theorem [@NaturePaper].

In a series of concomitant articles [@EBRtheory; @NaturePaper; @GraphDataPaper], in relation to the problem of topological insulators, we have studied the *connectivities* of energy bands of the EBRs. The energy bands are continuous functions defined in the Brillouin zone. The possible different ways in which the Bloch states at bands of $\mathbf{k}$-*vectors of maximal symmetry* (*cf.* Section \[kmax\]) are connected through intermediate bands of $\mathbf{k}$-vectors of lower symmetry, are restricted by the corresponding compatibility relations (*cf.* Section \[sec:compatibilityrelations\]). In most cases the EBRs are fully connected, which means that we can continuously ’travel’ along all the points that form the energy band. However, in some cases, depending on the specific values of the electronic band energies at each $\mathbf{k}$-vector, the compatibility relations allow the decomposition of the whole group of bands into disconnected branches, separated by an energy gap. Exactly such EBRs describe topological insulators. In a companion paper [@GraphDataPaper] we describe in detail the algorithms used to determine if an EBR (with and without TR) is decomposable or not.

\[sec:methodbandrep\]The method
-------------------------------

The algorithms implemented in the program [BANDREP]{} to calculate the BRs follows the site-symmetry procedure explained in Section \[sec:sitesymmetry\] and makes use of different tools of the BCS. Given a space group , a BR is fully identified by a Wyckoff position, an irrep of its site-symmetry group and the irreps of the little groups of any $\mathbf{k}$-vector of the reciprocal space. In fact, it is only necessary to consider the $\mathbf{k}$-vectors of maximal symmetry (*cf.* Section \[kmax\]), because the multiplicities of the irreps of the little groups of any other $\mathbf{k}$-vector can be calculated from the former ones using the corresponding compatibility relations described in Section \[sec:compatibilityrelations\]. The multiplicities of the irreps of the little groups of any $\mathbf{k}$-vector in BRs are determined applying the site-symmetry approach.

Given a space group $\mathcal{G}$ and a Wyckoff position $\mathbf{Q}=\{q\}$, first the program identifies the symmetry elements included in the site-symmetry group ${\mbox{$\mathcal{S}$}}_q$: for each of the coset representatives $\{R|\mathbf{v}\}$ of ${\mbox{$\mathcal{G}$}}:{\mbox{$\mathcal{T}$}}$ the program calculates $\mathbf{t}$ that satisfies the equation, $$\label{sitesymchek}
\mathbf{t}=Rq+\mathbf{v}-q.$$ If $\mathbf{t}$ belongs to the set of lattice translations of $\mathcal{G}$, then $\{R|\mathbf{v}-\mathbf{t}\}$ belongs to the site-symmetry group $\mathcal{S}_{q}$. Once all elements of $\mathcal{S}_{q}$ are determined, the tool [IDENTIFY GROUP]{} (www.cryst.ehu.es/cryst/identify\_group) of the BCS is used to identify the point (or double point) group type isomorphic to $\mathcal{S}_{q}$ and to establish the corresponding isomorphic mapping between the group elements. As a next step, for each element ${\mbox{\textit{\textsf{g}}}}\in\mathcal{S}_q$, we calculate: (i) the characters $\chi_\alpha({\mbox{\textit{\textsf{g}}}})$ of each irrep $\mathbf{d}_\alpha$ of the double point group (using [Representations DPG]{}), and (ii) the characters $\chi_\beta^{^{*}\mathbf{k}}({\mbox{\textit{\textsf{g}}}})$ of the irreps $D_\beta^{*\mathbf{k}}$ of $\mathcal{G}$ for each of the tabulated $\mathbf{k}$-vectors in the Brillouin zone, induced from the allowed irreps (or TR-invariant irreps) $\mathbf{d}_\beta^{\mathbf{k}}$ of the little group $\mathcal{G}^{\mathbf{k}}$ of $\mathbf{k}$ (applying the tool [Representations DSG]{}).

The decomposition of the subduced representations $\mathbf{D}_\beta^{*\mathbf{k}}\downarrow\mathcal{S}_{q}$ into irreps $\mathbf{d}_\alpha$ of the site-symmetry group ${\mbox{$\mathcal{S}$}}_q$ is performed using eqs. (\[eq1\]) and (\[eq2\]):

$$\label{eq1}
\mathbf{D}_\beta^{*\mathbf{k}}\downarrow\mathcal{S}_{q}\sim\bigoplus n_\alpha^{\mathbf{k},(\beta)}\mathbf{d}_\alpha,$$

where the multiplicities are, $$\label{eq2}
n_\alpha^{\mathbf{k},(\beta)}=\frac{1}{\left|\mathcal{S}_q\right|}\sum_{{\mbox{\textit{\textsf{g}}}}\in\mathcal{S}_q}\chi_\alpha({\mbox{\textit{\textsf{g}}}})^*\chi_\beta^{*\mathbf{k}}({\mbox{\textit{\textsf{g}}}}).$$

As explained in Section \[sec:sitesymmetry\], according to the Frobenius reciprocity theorem, these multiplicities coincide with the multiplicities of the irreps of the little group $\mathcal{G}^{\mathbf{k}}$ in the BR induced from the irrep $\mathbf{d}_\alpha$ of the double point group $\mathcal{S}_{q}$:

$$\label{eq:multiplicities}
(\mathbf{d}_\alpha\uparrow \mathcal{G})\downarrow\mathcal{G}^{\mathbf{k}}\sim\bigoplus n_{\beta}^{\mathbf{k},(\alpha)}   \mathbf{d}_\beta^{\mathbf{k}}.$$

We thus fully identify the BR induced from the irrep $\mathbf{d}_\alpha$ of the site-symmetry group of the Wyckoff position $\mathbf{Q}$ giving the tabulated irreps $\mathbf{d}_\beta^{\mathbf{k}}$ of the little group for every $\mathbf{k}$ of the space group $\mathcal{G}$ (obtained by [REPRESENTATIONS DSG]{}) and calculate the multiplicities $n_{\alpha}^{\mathbf{k},(\beta)}$.

Finally, we check if the BR induced from the Wyckoff position $\mathbf{Q}$ is elementary or not. As all EBRs of a space group are induced from Wyckoff positions of maximal symmetry, the check is reduced to BRs induced from such Wyckoff positions. Each BR is characterized by a set of multiplicities of the little-group irreps but only for $\mathbf{k}$-vectors of maximal symmetry. Using these lists of multiplicities we can easily check if the given BR decomposes at every **k** as the direct sum of two or more irreps. If such decompositions exist, the BR is a candidate to be a composite BR.

However, as explained in Section \[sec:bandrepresentations\], this condition is not sufficient to consider the BR as equivalent to the direct sum of BRs induced from a different Wyckoff position $\mathbf{Q}\,'$. In addition, it is necessary to calculate the intersection of the site-symmetry groups ${\mbox{$\mathcal{S}$}}_{\bf{Q}}$ and ${\mbox{$\mathcal{S}$}}_{\bf{Q}}\,'$ of $\mathbf{Q}$ and $\mathbf{Q}\,'$, respectively, and identify the Wyckoff position $\mathbf{Q}_0$ (of non-maximal symmetry) which has this intersection ${\mbox{$\mathcal{S}$}}_{\bf{Q}_0}$ as its site-symmetry group. Then, we induce representations of ${\mbox{$\mathcal{S}$}}_{\bf{Q}}$ and ${\mbox{$\mathcal{S}$}}_{\bf{Q}}\,'$ from every irrep of ${\mbox{$\mathcal{S}$}}_{\bf{Q}_0}$. If for some irrep of ${\mbox{$\mathcal{S}$}}_{\bf{Q}_0}$ we get $\rho$, the irrep of ${\mbox{$\mathcal{S}$}}_{\bf{Q}}$ that induces the candidate BR to be composite, and a reducible representation $\sigma$ of ${\mbox{$\mathcal{S}$}}_{\bf{Q}}\,'$, that induces a composite BR that decomposes at every $\mathbf{k}$ as $\rho\uparrow{\mbox{$\mathcal{G}$}}$ does, the candidate BR is composite. Otherwise, it is elementary.

\[wyckoffmaximal\]Wyckoff positions of maximal symmetry
-------------------------------------------------------

The so-called *Wyckoff positions of maximal symmetry* are of importance for certain applications, and in particular, for the calculation of the EBRs. In the following, we comment briefly on the definition of Wyckoff positions of maximal symmetry The intersection $\mathcal{S}_{\mathbf{Q}_0}=\mathcal{S}_{\mathbf{Q}}\cap\mathcal{S}_{\mathbf{Q}\,'}$ of the two site-symmetry groups is the site-symmetry group of another Wyckoff position $\mathbf{Q}_0$ that can be identified relatively easy: some of the point coordinates of $\mathbf{Q}_0$ should be represented by variable parameters that interpolate between the point coordinates of the Wyckoff positions $\mathbf{Q}$ and $\mathbf{Q}\,'$. If for a given irrep $\rho_0$ of $\mathcal{S}_{\mathbf{Q}_0}$, the induced representations into the site-symmetry groups of $\mathbf{Q}$ and $\mathbf{Q}\,'$ satisfy, $\rho_0\uparrow S_{\mathbf{Q}}=\rho$ and $\rho_0\uparrow S_{\mathbf{Q}\,'}=\sigma$, then the two BRs $(\rho\uparrow G)\downarrow G^{\mathbf{k}}=\rho_G^{\mathbf{k}}$ and $(\sigma\uparrow G)\downarrow G^{\mathbf{k}}=\sigma_G^{\mathbf{k}}$ are equivalent.and on the procedure we used for their calculation. Note that the sets of Wyckoff positions of double space groups and those of the corresponding space groups are closely related. In fact, the essential difference concerns the site-symmetry groups which, in the former case, are isomorphic to double point groups, while in the latter, only to point groups. As a consequence, practically the same procedure for the determination of the Wyckoff positions of maximal symmetry can be used for space groups and for double space groups.

It is common to describe a Wyckoff position by its multiplicity, Wyckoff letter, symbol of the site-symmetry group and a set of coordinate triplets of the points in the unit cell that belong to the Wyckoff position, possibly depending on one or two variable parameters (three for the general position) (for a detailed introduction to Wyckoff positions of space groups, *cf.* [[*IT*]{}A ]{}). A Wyckoff position $\bf{Q}$ with a site-symmetry group ${\mbox{$\mathcal{S}$}}_{\bf{Q}}$ has maximal symmetry if it is not *connected* to another Wyckoff position $\bf{Q}\,'$ whose site-symmetry group ${\mbox{$\mathcal{S}$}}_{\bf{Q}\,'}$ is a supergroup of ${\mbox{$\mathcal{S}$}}_{\bf{Q}}$. We say that two Wyckoff positions are *connected* if (i) the coordinate triplet of at least one of them depends on one or more variable parameters, and (ii) if for specific values of the variable parameters the coordinate triplets of the two Wyckoff positions coincide. For instance, in the space group P$2/m$ (No. 10), the site-symmetry group of the Wyckoff position $2i: (0,y,0)$ is isomorphic to the point group $2$. It is not a Wyckoff position of maximal symmetry because it is *connected* to the Wyckoff position $2a: (0,0,0)$: the coordinate-triplet description of $2a$ is obtained from that of $2i$ by the substitution $y=0$, and the site-symmetry group of the position $2a$ is isomorphic to $2/m$ which is supergroup of $2$.

The algorithm to identify the Wyckoff positions of maximal symmetry is straightforward. In the space groups with no points in a special position (the so-called *fixed-point-free* or *Bieberbach* groups) the general Wyckoff position is the only Wyckoff position of maximal symmetry. For the rest of space groups, we distribute the special Wyckoff positions into three subsets, according to the number $(0,1,2)$ of variable parameters of their coordinate triplets. All Wyckoff positions with $0$ variable parameters are Wyckoff positions of maximal symmetry. Those Wyckoff positions of the subset with one variable parameter which are *connected* to at least one Wyckoff position of the subset with no variable parameter, in the sense explained above, are not maximal. The rest of Wyckoff positions with one variable parameter are maximal. Finally, we repeat the check for the subset of Wyckoff positions with two variable parameters, trying to find if they are *connected* to at least one of the Wyckoff positions of the subsets with 0 and 1 variable parameter. Those that have no *connection* are Wyckoff positions of maximal symmetry.

As an example, consider the Wyckoff positions of the (double) space group P$4/ncc$ (No. 130) shown in Table \[table:wyckoffmax\]. The Wyckoff positions $4a,\, 4b,\,4c$ and $8d$ are of maximal symmetry while $8e,\, 8f$ and $16g$ are not. If a Wyckoff position is not maximal, then the *corresponding* Wyckoff position of maximal symmetry (*i.e.* the one to which it is *connected*) together with the specific values of the variable parameters for which the two coordinate-triplet descriptions coincide, are indicated in the last column of the table.

  ---------- ---------------- ----------------- -----------------
   Wyckoff    Site-symmetry      coordinate          maximal
   position       group            triplet          symmetry
     $4a$       $222(D_2)$     $(3/4,1/4,1/4)$         yes
     $4b$     $\bar{4}(S_4)$    $(3/4,1/4,0)$          yes
     $4c$        $4(C_4)$       $(1/4,1/4,z)$          yes
     $8d$     $\bar{1}(C_i)$      $(0,0,0)$            yes
     $8e$        $2(C_2)$       $(3/4,1/4,z)$          no
                                                  $z=1/4\to 4a$
     $8f$        $2(C_2)$      $(x+1/2,x,1/4)$         no
                                                  $x=1/4\to 4a$
    $16g$        $1(C_1)$         $(x,y,z)$            no
                                                 $x=y=z=0\to 8d$
  ---------- ---------------- ----------------- -----------------

  : Table 1: Wyckoff positions of the (double) space group P$4/ncc$ (No. 130). Each Wyckoff position is specified by its multiplicity and Wyckoff letter (first column), the Hermann-Mauguin (Schönflies) symbol of its site-symmetry group (second column), the coordinate-triplet description of a representative of the orbits of the Wyckoff positions (third column). In the last column it is indicated if the Wyckoff position is of maximal symmetry or not. For the Wyckoff positions of non-maximal symmetry, the *corresponding* Wyckoff position of maximal symmetry, together with the specific values of the variable parameters for which the two coordinate-triplet descriptions coincide, are also specified.[]{data-label="table:wyckoffmax"}

\[kmax\]$\mathbf{k}$-vectors of maximal symmetry
------------------------------------------------

In analogy to the Wyckoff positions of maximal symmetry introduced in the previous section, for each space group we can define a set of $\mathbf{k}$-*vectors of maximal symmetry*. These vectors play also an important role, for example, in the analysis of the connectivities of the BRs (*cf.* [@GraphTheoryPaper; @GraphDataPaper; @EBRtheory]).

Similar to the distribution of points of direct space into Wyckoff positions, the set of all $\mathbf{k}$-vectors can be distributed into the so-called $\mathbf{k}$-*vector types*. A $\mathbf{k}$-vector type consists of complete orbits of $\mathbf{k}$-vectors and thus of full stars of $\mathbf{k}$-vectors. The $\mathbf{k}$-vectors belonging to a $\mathbf{k}$-vector type are represented by a $\mathbf{k}$-vector letter (here we follow CDML notation), by the point-group type of the little co-groups of $\mathbf{k}$-vectors, and by a set of $\mathbf{k}$-vectors coefficients: the zero, one or two variable parameters in the $\mathbf{k}$-vector coefficients correspond to special $\mathbf{k}$-vector types, *i.e.*, they define symmetry points, symmetry lines or symmetry planes in the Brillouin zone. Three variable parameters indicate a general $\mathbf{k}$-vector type. We say that a $\mathbf{k}$-vector type (or just a $\mathbf{k}$-vector, for short) is *of maximal symmetry* if its little co-group is not a subgroup of the little co-group of another $\mathbf{k}\,'$-vector type *connected* to $\mathbf{k}$ in the same sense as the *connected* Wyckoff positions discussed in Section \[wyckoffmaximal\]. The procedure to identify the $\mathbf{k}$-vectors of maximal symmetry is analogous to the procedure for the determination of the Wyckoff positions of maximal symmetry.

In general, the set of $\mathbf{k}$-vectors of maximal symmetry for non-centrosymmetric space groups (*i.e* space groups whose point groups $\overline{{\mbox{$\mathcal{G}$}}}$ do not include the operation of (space) inversion $\overline{1}$) could be modified when time-reversal symmetry is taken into account. One can show that in such cases wave vectors of the same $\mathbf{k}$-vector type (with respect to the spacial symmetry) could behave differently under the action of time reversal: some wave vectors are time-reversal invariant, the so-called *Time Reversal Invariant Momentum* (TRIM) points, while others are not. In other words, when TR symmetry is taken into account, the TRIM points are $\mathbf{k}$-vectors of maximal symmetry and the corresponding physically irreducible representations determine the transformation properties of the eigenfunctions of the Hamiltonian of the system. As an example, consider the polar space group ${\mathit{P}4}$ (No. 75) and the symmetry line $\Lambda\,(0,0,w)$ which is a $\mathbf{k}$-vector line of maximal symmetry as its little co-group is the point group 4 ($C_4$). There are 4 one-dimensional single-valued irreps at the $\Lambda\,(0,0,w)$ point: $\Lambda_i$, $i=1,\ldots,4$. If however, time-reversal symmetry is taken into account, the $\mathbf{k}$-vector type $\Lambda\,(0,0,w)$ ’splits’ into three $\mathbf{k}$-vector types: the points $\Gamma\, (0,0,0)$ and $Z\,(0,0,1/2)$ are TRIM points and become $\mathbf{k}$-vector types of maximal symmetry while the rest of the points of the line $\Lambda$ form a $\mathbf{k}$-vector type of *non*-maximal symmetry. Two of the four single-valued irreducible representations at these points, $\Gamma_3$, $\Gamma_4$ and $Z_3$, $Z_4$, form a pair of complex conjugated irreps and become doubly-degenerate. Tables \[table:kvecsnoTR\] and \[table:kvecsTR\] show the lists of $\mathbf{k}$-vectors of maximal symmetry of ${\mathit{P}4}$ (No. 75), without and with TR, respectively.

The above example indicates the important consequences of the time-reversal symmetry in the derivation of the compatibility relations, and in particular, in the study of EBR connectivities between pairs of $\mathbf{k}$-vectors of maximal symmetry.

     label     little co-group   coefficients
  ----------- ----------------- ---------------
   $\Lambda$      $4(C_4)$         $(0,0,w)$
      $W$         $4(C_4)$        $(0,1/2,w)$
      $V$         $4(C_4)$       $(1/2,1/2,w)$

  : Table 2: $\mathbf{k}$-vectors of maximal symmetry of the (double) space group P$4$ (No. 75) if time-reversal symmetry is not considered. The $\mathbf{k}$-vector labels follow the notation of CDML. The second column gives the little co-group of the $\mathbf{k}$-vector. The third column shows the coefficients of a representative of the $\vec{k}$-vector star.[]{data-label="table:kvecsnoTR"}

    label     little co-group    coefficients
  ---------- ----------------- -----------------
   $\Gamma$      $4(C_4)$          $(0,0,0)$
     $Z$         $4(C_4)$         $(0,0,1/2)$
     $X$         $4(C_4)$         $(0,1/2,0)$
     $R$         $4(C_4)$        $(0,1/2,1/2)$
     $M$         $4(C_4)$        $(1/2,1/2,0)$
     $A$         $4(C_4)$       $(1/2,1/2,1/2)$

  : Table 3: $\mathbf{k}$-vectors of maximal symmetry of the (double) space group P$4$ (No. 75) with time-reversal symmetry taken into account. The $\mathbf{k}$-vector labels in the first column follow the notation of CDML. The second column gives the little co-group of the $\mathbf{k}$-vector. The third column shows the coefficients of a representative of the $\vec{k}$-vector star.[]{data-label="table:kvecsTR"}

\[bandrep\]The program [BANDREP]{}
----------------------------------

The program [BANDREP]{}, recently added to BCS, applies the method explained in Section \[sec:methodbandrep\] to calculate the BRs of any of the 230 double space groups. The fact that the list of single-valued BRs identified as elementary by the program, coincides exactly with the list obtained by , alongside with the many checks performed and the agreement between the program and the results of , can be considered a proof test of the program.

#### Input:

The main page of the input requires the specification of double space group (by its [[*IT*]{}A ]{}sequential number) and offers four options: calculation either of EBRs of the double space group with or without TR symmetry, or of the BRs with or without TR symmetry. For the last two options the program produces the list of Wyckoff positions of the selected space group, separated into those of maximal and non-maximal symmetry. Figures \[fig:bandrepinput1\] and \[fig:bandrepinput2\] show screenshots of the input pages.

![Figure 14: Screenshot of the first input form of the program [BANDREP]{}. The input data include (i) specification of the space group, and (ii) choice among the four options: (1) elementary band representations (EBR) without time-reversal (TR); (2) EBRs with TR; (3) band representations (BR) without TR; and (4) BRs with TR.[]{data-label="fig:bandrepinput1"}](./Figure14.png){width="\textwidth"}

![Figure 15: Screenshot of the second input form of the program [BANDREP]{} which shows the set of Wyckoff positions of the specified space group divided into positions of maximal and non-maximal symmetry. This input form is shown only for the options (3) and (4) of the main input (see figure \[fig:bandrepinput1\]). []{data-label="fig:bandrepinput2"}](./Figure15.png){width="70.00000%"}

#### Output:

The option selected by the user determines the specific output produced by the program:

1.  Option: *Elementary band representations without TR symmetry*

    A screenshot of the table output of the program [BANDREP]{} for the double space group ${\mathit{P}4/ncc}$ (No. 130) is shown on Figure \[fig:bandrepoutput1\]. The EBRs are listed in columns specified by (i) the Wyckoff positions of ${\mathit{P}4/ncc}$ (No. 130) followed by the symbol of the double point group ${\mbox{$\mathcal{S}$}}_{\mathbf{Q}}$ isomorphic to the site-symmetry group (first row of the header), and (ii) the irrep $\mathbf{d}_\alpha$ of the site-symmetry group from which the EBR $\mathbf{d}_\alpha\uparrow {\mbox{$\mathcal{G}$}}$ is induced (second row of the header); the dimension of $\mathbf{d}_\alpha\uparrow {\mbox{$\mathcal{G}$}}$ is shown in brackets after the EBR symbol. The entries of the third row of the header indicate if the EBR is decomposable or not (*cf.* Section \[bandrepproblem\]). The entries of the output table show the decompositions of EBRs into irreps of the little groups of the $\mathbf{k}$-vectors of maximal symmetry which denote the rows of the table. The dimensions of the little-group irreps are given in brackets after their symbols.

    ![Figure 16: Screenshot of the output given by the program [BANDREP]{} for $P4/ncc$ (No. 130) and for the option ’elementary band representations without time-reversal symmetry’ (*cf.* Fig. \[fig:bandrepinput1\]). For a detailed description of the displayed data, see Section \[bandrep\]. Only part of the wider output has been included in the figure.[]{data-label="fig:bandrepoutput1"}](./Figure16.png){width="\textwidth"}

    Clicking on the button “Show all types of $\mathbf{k}$-vectors” generates a table with the decompositions of each EBR into irreps of the little groups of maximal and non-maximal $\mathbf{k}$-vectors.

    Clicking on the button “Minimal set of paths and compatibility relations to analyse the connectivity”, produces the set of paths between pairs of $\mathbf{k}$-vectors of maximal symmetry to be considered in the analysis of the connectivities of the EBRs in the given space group, once the known redundancies have been removed [@GraphDataPaper], and the corresponding independent compatibility relations. Figure \[fig:paths1\] shows the screenshot of the minimal set of paths and compatibility relations for the space group ${\mathit{P}4/ncc}$. For more details about the problem of connectivity of EBRs see .

    ![Figure 17: Screenshot of the output given by the program [BANDREP]{} which shows the independent sets of compatibility relations along intermediate paths between the $\mathbf{k}$-vectors of maximal symmetry of $P4/ncc$ (No. 130), necessary for the analysis of the connectivity structure of elementary band representations. Only part of the longer output has been included in the figure.[]{data-label="fig:paths1"}](./Figure17.png){width="\textwidth"}

2.  Option: *Elementary band representations with TR symmetry*

    The essential difference in comparison with the output of Option 1 is that the shown results refer to *physically-irreducible* representations (for a discussion on physically-irreducible representations, *cf.* Section \[subsec:complexconjugation\]).

3.  Option: *Band representations from a Wyckoff position without TR symmetry*

    Under this option the program shows the band representations for a specific Wyckoff position chosen from a list of Wyckoff positions of maximal or of non-maximal symmetry. Fig. \[fig:bandrepoutput2\] shows the output produced by the program for the case of the double space group ${\mathit{P}4/ncc}$, and the Wyckoff position of maximal symmetry $8d$. The structure of the output table is similar to those of elementary bands: the essential difference is that the program shows not only the elementary but also the composite band representations (specified in the second row of the table header). Whether the EBRs are decomposable or not is also indicated. In the partial output shown in Fig. \[fig:bandrepoutput2\] all the EBRs are decomposable. By clicking on the button “Decomposable”, we obtain all different ways to decompose the EBR into disconnected branches, which correspond to topological insulators. The four possible ways to decompose the EBR $A_g\uparrow {\mbox{$\mathcal{G}$}}$ induced from the $A_g$ irreducible representation of the site-symmetry group of the $8d$ Wyckoff position of the ${\mathit{P}4/ncc}$ space group are shown in Fig. \[fig:bandrepoutputdesc\]. Each column gives the arrangements of the set of the little-group irreps at each $\mathbf{k}$-vector of maximal symmetry in a branch. The band-graph (not yet available in the BCS) that illustrates a disconnected EBR is shown in Fig. \[fig:bandgraph\]. The arrangements of irreps at each $\mathbf{k}$-vector in Fig. \[fig:bandgraph\] correspond to the decomposition shown on the first row of Fig. \[fig:bandrepoutputdesc\].

    ![Figure 18: Screenshot of the output given by the program [BANDREP]{} which shows the band representations of $P4/ncc$ (No. 130) and Wyckoff position $8d$. For a detailed description of the displayed data, see Section \[bandrep\][]{data-label="fig:bandrepoutput2"}](./Figure18.png){width="\textwidth"}

    ![Figure 19: Screenshot of the output given by the program [BANDREP]{} which shows all (four) different ways to decompose the elementary band representation induced from the $A_g$ representation of the site-symmetry group of the Wyckoff position $8d$ of $P4/ncc$ (No. 130). The elementary band representation can be decomposed into two further indecomposable branches (branch 1 and branch 2 in the figure). Each column of the table gives the arrangements of the set of the little-group irreps at each $\mathbf{k}$-vector of maximal symmetry in a branch.[]{data-label="fig:bandrepoutputdesc"}](./Figure19.png){width="60.00000%"}

    ![Figure 20: Band-graph that illustrates a disconnected elementary band representation induced from the $A_g$ representation of the site-symmetry group of the Wyckoff position $8d$ of $P4/ncc$ (No. 130). The arrangements of irreps at each $\mathbf{k}$-vector correspond to the decomposition shown in the first row of Fig. \[fig:bandrepoutputdesc\]. The interconnections between the little-group irreps of neighboring $\mathbf{k}$-vectors are drawn in accordance to the corresponding compatibility relations.[]{data-label="fig:bandgraph"}](./Figure20.png){width="\textwidth"}

4.  Option: *Band representations from a Wyckoff position with TR symmetry*. This provides similar information as Option 3 but for the *physically-irreducible* representations.

Conclusions
===========

The group-theoretical description of physical systems where the Hamiltonian depends on spin components, would require the use of the so-called *double crystallographic groups* and their single- and double-valued representations. In this paper, we describe a set of databases and programs of double crystallographic groups that recently have been implemented in the Bilbao Crystallographic Server (http://www.cryst.ehu.es). As the rest of the programs on BCS, the new tools are freely available and can be accessed via user-friendly web interfaces. Some of the algorithms applied in the new programs are extensions of the algorithms used in the BCS for ordinary space groups.

The tool [DGENPOS]{} provides in different formats the symmetry operations of the 230 double space groups in the standard or conventional setting. The program [REPRESENTATIONS DPG]{} gives access to the irreducible representations of the 32 crystallographic double point groups while [REPRESENTATIONS DSG]{} calculates the irreducible representations of the double space groups and analyses their *reality* indicating also the pairs of conjugated irreducible representations. The program [DSITESYM]{} applies the site-symmetry approach to the double space groups. It can be considered as a bridge between a local description of the atomic orbitals on-site and a global description through extended states along the Brillouin zone in a structure. The program [DCOMPREL]{} calculates the compatibility relations between the irreducible representations of double space groups at high- and low-symmetry points in the Brillouin zone, necessary in the analysis of the connectivity of the functions defined in the reciprocal space (*e.g.* in the analysis of the structure of the electronic bands and their connectivity through the Brillouin zone). As an application of the site-symmetry approach, we have developed the program [BANDREP]{} that provides the band representations and the subset of elementary band representations induced from any Wyckoff position of any double space group. The tool also identifies the subset of decomposable elementary band representations and the different ways of their decomposition, together with large amount of additional data necessary for their study. Concomitantly with the current paper, the results provided by [BANDREP]{} have been successfully applied in a novel method for the description, search and prediction of topological insulators [@NaturePaper; @GraphTheoryPaper; @GraphDataPaper; @EBRtheory].

Normal-subgroup induction procedure
===================================

The irreps of a space group ${\cal{G}}$ are obtained by induction from the irreps of its translation group ${\cal{T}}$. Assuming the Born-von Karman (periodic) boundary conditions $(\textbf{I}, \textbf{t}_i)^{N_i} = (\textbf{I}, \textbf{o})$ to hold, where **t**$_i$ = (1,0,0), (0,1,0) or (0,0,1) and $N_i$ is a large integer for *i* = 1, 2 or 3, respectively. Then, the irreps $\mathbf{\Gamma}^{\mathbf{k}}$ of the translation group ${\cal{T}}$ are given by: $$\mathbf{\Gamma}^\textbf{k} [ (\textbf{I},\textbf{t})] = \mathrm{exp}(-i\textbf{k}\cdot\textbf{t}).$$ There are $N_1N_2N_3$ different irreps of ${\mbox{$\mathcal{T}$}}$ which are distinguished by the wave vectors: $$\label{eqkve}
\mathbf{k} = \sum_{i=1}^3 k_{i}\,{\bf a}_{i}^{*},$$ where $ k_{i} = q_i/N_i;\ q_i=0,\,1,\,2,\,\ldots\,,\,N_i-1$. The basis **a**$_1^{*}$, **a**$_2^{*}$, **a**$_3^{*}$ is called the basis of the *reciprocal lattice* **L**$^*$ and it is the *dual basis* of $\mathbf{a}_1$, **a**$_2$, **a**$_3$ of **L**. The vectors $\mathbf{a}_i^*$ are defined by the relations $\mathbf{a}_i\cdot\mathbf{a}_j^*=2\pi\delta_{ij}$ where $\delta_{ij}$ is the Kronecker symbol.

The wave vectors **k** and **k’** = **k** + **K**, where **K** is a vector of the reciprocal lattice **L**$^*$, describe the same irreps of $\cal{T}$. Therefore, in order to determine all the irreps of ${\cal{T}}$, it is necessary to consider only the **k**-vectors of the first Brillouin-zone.

The *little co-group* of **k** is the point group consisting of all the rotational parts $R^\textbf{k}$ of the symmetry operations of the space group ${\cal{G}}$ that either leave the **k**-vector invariant, or map it to an equivalent vector, *i.e.*: $${\bf k} = {\bf k}{\mbox{\bfseries{\itshape{W}}}}^{\mathbf{k}} + {\bf K,~~ K} \in {\bf L}^*. \label{def_little_cogroup}$$ Here, ${\mbox{\bfseries{\itshape{W}}}}^{\mathbf{k}}$ is the $(3\times3)$ matrix representation of $R^\textbf{k}$. The little co-group $\overline{{\cal{G}}}^\textbf{k}$ is a subgroup of the point group $\overline{{\cal{G}}}$ of the space group ${\cal{G}}$. The vector [**k**]{} is called a *general* **k** *vector* if the little co-group contains the identity operation only, *i.e.* $\overline{\cal{G}}\,^{\mathbf{k}} = \{\cal{I}\}$; otherwise $\overline
{\cal{G}}\,^{\mathbf{k}} > \{\cal{I}\}$, and [**k**]{} is called a *special* **k** *vector*.

Consider the coset decomposition of $\overline{{\mbox{$\mathcal{G}$}}}$ relative to $\overline{{\mbox{$\mathcal{G}$}}}\,^{\mathbf{k}}$. If $\{R_m\}$ is the corresponding set of coset representatives, then the set $\mathbf{*k}=\{{\bf k}{\mbox{\bfseries{\itshape{W}}}}_m+\mathbf{K}\}$, with $\mathbf{K}\in$**L**$^*$, is called the *star of [**k**]{}* and the vectors $\mathbf{k}{\mbox{\bfseries{\itshape{W}}}}_m+\mathbf{K}$ are called the *arms of* $\mathbf{*k}$.

An orbit of $\mathbf{\Gamma}^{\mathbf{k}}({\mbox{$\mathcal{T}$}})$ under conjugation by comprises all irreps $\mathbf{\Gamma}^{\mathbf{k}'}({\mbox{$\mathcal{T}$}})$ with $\mathbf{k}'$ belonging to $\mathbf{*k}$. Irreps of belonging to the same orbit give rise to equivalent irreps of , *i.e.* in order to obtain each irrep of exactly once it is necessary to consider one **k** vector per star. (A simply connected subset of the Brillouin zone which contains exactly one **k** vector per $\mathbf{*k}$, is called a *representation domain*.)

Given a space group , its translation subgroup , and an irrep $\mathbf{\Gamma}^{\mathbf{k}}({\mbox{$\mathcal{T}$}})$, one can define the *little group* ${\mbox{$\mathcal{G}$}}^{\mathbf{k}}$ of the wave vector **k**: it is a space group that consists of all those elements of whose rotation parts $R^{\mathbf{k}}$ leave either **k** unchanged or transform it into an equivalent vector,

$${\mbox{$\mathcal{G}$}}^{\mathbf{k}}=\{\{R^{\mathbf{k}}|\mathbf{v}^{\mathbf{k}}\} \in 
{\mbox{$\mathcal{G}$}}|R^{\mathbf{k}} \in \overline{{\mbox{$\mathcal{G}$}}}\,^{\mathbf{k}}\}.\label{litgr}$$

The irreps of space groups are obtained by induction from the so-called *allowed* irreps of the little groups ${\mbox{$\mathcal{G}$}}^{\mathbf{k}}$ of **k**. If ${\mbox{\bfseries{\itshape{D}}}}^{\mathbf{k},i}$ is an allowed irrep of ${\mbox{$\mathcal{G}$}}^{\mathbf{k}}$, then ${\mbox{\bfseries{\itshape{D}}}}^{\mathbf{k},i}
(\{1|\,\mathbf{t}\})=\exp{(-i\,{\bf k\cdot\,t})}\,{\mbox{\bfseries{\itshape{I}}}}$. (The matrix is the identity matrix with $\dim{\mbox{\bfseries{\itshape{I}}}}=\dim{\mbox{\bfseries{\itshape{D}}}}^{\mathbf{k},\,i}$).

The allowed irreps of the little group ${\mbox{$\mathcal{G}$}}^{\mathbf{k}}$ are determined by an induction procedure [@zak1960] which is based on the fact that the little groups (as all crystallographic groups) are *solvable* groups, *i.e.* for each group ${\mbox{$\mathcal{G}$}}^{\mathbf{k}}$ there exists a series of subgroups ${\mbox{$\mathcal{H}$}}_i^{\mathbf{k}}$ (the so-called *composition series*), such that:

$$\label{comp_little}
{\mbox{$\mathcal{G}$}}^{\mathbf{k}}\rhd{\mbox{$\mathcal{H}$}}_1^{\mathbf{k}}\rhd\ldots\rhd{\mbox{$\mathcal{H}$}}_{m-1}^{\mathbf{k}}
\rhd{\mbox{$\mathcal{H}$}}_{m}^{\mathbf{k}}\rhd\ldots\rhd{\mbox{$\mathcal{H}$}}_n^{\mathbf{k}}={{\mbox{$\mathcal{T}$}}}$$

and that the factor groups ${\mbox{$\mathcal{H}$}}_{m-1}^{\mathbf{k}}/{\mbox{$\mathcal{H}$}}_{m}^{\mathbf{k}}$ are cyclic groups of prime order. The (allowed) irreps of ${\mbox{$\mathcal{G}$}}^{\mathbf{k}}$ can be obtained from the (allowed) irreps of ${{\mbox{$\mathcal{T}$}}}$ by applying several times the general induction procedure ’climbing up’ the chain of normal subgroups (eq. \[comp\_little\]). Important for the irrep calculation is the observation that the factor groups in the composition series of crystallographic groups have orders 2 or 3 which simplifies considerably the induction procedure. The corresponding induction formulae and a detailed example of application of the induction method in the case of the space group $P4bm$ and $\mathbf{k}=\mathbf{X}(0,1/2,0)$ can be found, for example, in Aroyo *et al.* (2006).

Finally, following the normal-subgroup induction procedure, the irreps of a space group (the so-called *full* irreps) for a given **k** vector are obtained by induction from the allowed irreps $\mathbf{D}^{\mathbf{k},\,i}$ of the corresponding little group ${\mbox{$\mathcal{G}$}}^{\mathbf{k}}$. Let the elements ${\mbox{\textit{\textsf{q}}}}_m=\{R_m|\,\mathbf{v}_m\},\ \ m=1,\ \ldots\ ,s$ be the representatives of the cosets in the decomposition of relative to ${\mbox{$\mathcal{G}$}}^{\mathbf{k}}$:

$$\label{coset}
{\cal{G}} = {\cal{G}}^{\mathbf{k}} \cup {\mbox{\textit{\textsf{q}}}}_2{\cal{G}}^{\mathbf{k}}\cup \dots \cup {\mbox{\textit{\textsf{q}}}}_s{\cal{G}}^{\mathbf{k}}.$$

If $dim\,{\mbox{\bfseries{\itshape{D}}}}^{\mathbf{k},\,i}=r$, and if $s$ is the number of arms in (the *order* of) the star of **k**, then the induced irrep ${\mbox{\bfseries{\itshape{D}}}}^{\mathbf{*k},\,i}({\mbox{$\mathcal{G}$}})$ has the dimension $r\times s$ and its matrices can be written in the form: $$\label{sgirrep}
{\mbox{\bfseries{\itshape{D}}}}^{\mathbf{*k},\,i}(\{R|\,\mathbf{v}\}_{mp,nq})=
{\mbox{\bfseries{\itshape{M}}}}(\{R|\,\mathbf{v}\}_{m,n}){\mbox{\bfseries{\itshape{D}}}}^{\mathbf{k},\,i}
(\{R^{\mathbf{k}}|\,\mathbf{v}^{\mathbf{k}}\}_{p,q}) \,,$$ where $\{R^{\mathbf{k}}|\,\mathbf{v}^{\mathbf{k}}\}=
({\mbox{\textit{\textsf{q}}}}_m)^{-1}\,\{R|\,\mathbf{v}\}\,{\mbox{\textit{\textsf{q}}}}_n$ is an element of the little group ${\mbox{$\mathcal{G}$}}^{\mathbf{k}}$, with $n,m=1,\ldots,s$. Because the $s \times s$ *induction matrix* ${\mbox{\bfseries{\itshape{M}}}}(\{R|\,\mathbf{v}\})$ is a monomial matrix, the matrices ${\mbox{\bfseries{\itshape{D}}}}^{\mathbf{*k},\,i}(\{R|\,\mathbf{v}\})$ have a block structure with exactly one non-zero $(r \times r)$ block in every column and every row; the block is the matrix ${\mbox{\bfseries{\itshape{D}}}}^{\mathbf{k},\,i}(\{R^{\mathbf{k}}|\,\mathbf{v}^{\mathbf{k}}\})$, and $\{R^{\mathbf{k}}|\,\mathbf{v}^{\mathbf{k}}\}$ is fixed by the condition $\{R^{\mathbf{k}}|\,\mathbf{v}^{\mathbf{k}}\}=({\mbox{\textit{\textsf{q}}}}_m)^{-1}\,\{R|\,\mathbf{v}\}\,{\mbox{\textit{\textsf{q}}}}_n\,\in{\mbox{$\mathcal{G}$}}^{\mathbf{k}}$.

The work of LE, GF and MIA was supported by the Government of the Basque Country (project IT779-13) and the Spanish Ministry of Economy and Competitiveness and FEDER funds (project MAT2015-66441-P). The work of MVG was supported by FIS2016- 75862-P and FIS2013-48286-C2-1-P national projects of the Spanish MINECO. ZW and BAB, as well as part of the development of the initial theory and further ab-initio work, were supported by the NSF EAGER Grant No. DMR-1643312, ONR - N00014-14-1-0330, ARO MURI W911NF-12-1-0461, and NSF-MRSEC DMR-1420541. The development of the practical part of the theory, tables, some of the code development, and ab-initio work was funded by Department of Energy de-sc0016239, Simons Investigator Award, the Packard Foundation, and the Schmidt Fund for Innovative Research. BB, JC, ZW, and BAB acknowledge the hospitality of the Donostia International Physics Center, where parts of this work were carried out. JC acknowledges the hospitality of the Kavli Institute for Theoretical Physics, and BAB acknowledges the hospitality and support of the École Normale Supérieure and Laboratoire de Physique Théorique et Hautes Energies.
