<figure><img src="image_1.png" /></figure>

HAL
open science

# Bicrystallography and Beyond: Example of Group–Subgroup Phase Transformations

Denis Gratias, Marianne Quiquandon

► To cite this version:

Denis Gratias, Marianne Quiquandon. Bicrystallography and Beyond: Example of Group–Subgroup Phase Transformations. Crystals, 2020, 10, 10.3390/cryst10070560. [https://doi.org/10.3390/cryst10070560]

HAL Id: hal-03013167

https://hal.science/hal-03013167

Submitted on 18 Nov 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Bicrystallography and Beyond: Example of Group–Subgroup Phase Transformations

Denis Gratias * and Marianne Quiquandon

CNRS UMR 8247, Institut de Recherche de Chimie Paris, ENSCP, 11 rue Pierre et Marie Curie, F-75005 Paris, France; marianne.quiquandon@chimieparistech.psl.eu

* Correspondence: denis.gratias@chimieparistech.psl.eu

Received: 3 June 2020; Accepted: 28 June 2020; Published: 1 July 2020

**Abstract:** This paper presents the basic elementary tools for describing the global symmetry obtained by overlapping two or more crystal variants of the same structure, differently oriented and displaced one with respect to the other. It gives an explicit simple link between the concepts used in the symmetry studies on grain boundaries on one side and group-subgroup transformations on the other side. These questions are essentially of the same nature and boil down to the resolution of the same problem: identifying the permutation groups that are images of the corresponding applications. Examples are given from both domains, classical grain boundaries with coincidence lattices and group-subgroup phase transformations that illustrate the profound similarities between the two approaches.

**Keywords:** interfaces; bicrystal; N-crystal; group action theory; permutations homorphisms

## 1. Introduction

Understanding how geometric order can propagate in solids, for example with the generation of fascinating Moiré patterns in the superimposition of two or more identical twisted lattices, the specific relative orientations between twinned crystals or the distribution of ordered variants in the morphology of an initial single crystal after a phase transition, is a major part of describing the properties of coherent interfaces in crystals. The subject has met a particular revival since the discovery [1] of bilayers of twisted extended 2D structure like graphene and similar structures with very interesting electronic properties of the bilayer depending crucially on its symmetry. Our intention here is to give some general theoretical elements of an answer to this fascinating symmetry problem.

The notion of *bicrystal* is issued from the studies of grain boundaries in metals and minerals as first discussed by Bollman [2] after the seminal approach of Friedel on the symmetries of twinned crystals (see for instance [3,4]). First designated as dichromatic pattern (see Pond and Bollman [5] and Pond & Vlachavas [6]), the notion of bicrystal was implicitly based on the lattice symmetry mostly because examples at that time were specific studies of grain boundaries of simple metals the atomic structure of which can be described by symmorphic space groups with one unique atom per unit cell. This lack of generality has been discussed by Gratias & Portier [7] who used space groups and group action theory [8] to construct a general crystallographic framework for describing the geometry of homophase grain boundaries. In all these approaches, the basic idea was to decipher which symmetries appear in the abstract construction of the superimposition of two crystals of the same nature in different orientations for finding which elementary domain in space should be examined to discuss all the possible interfaces—grain boundaries or twins—that can actually exist between the two crystals.

The overall discussion addressed in the present paper takes its roots on the works of Guymont et al. [9] and Gratias & Portier [7] giving the basic tools for handling the symmetry properties of bicrystals using elementary group action theory [8]. Our present goal is to generalize this notion to
that of N-crystals, an abstract construct built of the superimposition of N identical crystals displaced and disoriented with respect to each other. This generalization allows the inclusion of the description of symmetry breaking in solids due to group-subgroup phase transformations that generate an assembly of crystalline domains of the low-symmetry phase out of a single crystal of the high-symmetry phase. Classical example of these phase transformation are the many chemical order-disorder transitions in metallic alloys. This leads to a complete geometric description of the interconnection between the variants generated by a group-subgroup phase transition.

The paper has three main sections:

* the first section “Standard bicrystallography” discusses the basic tools used to characterize interfaces in homogeneous crystals; a special attention is given to the building of the space group W$_{α}$ that define all possible space groups P$_{τ}$ of the bicrystal as function of the rigid-body translation τ; another key point here is to emphasize the fact that the symmetry group of the bicrystal results from only two kind of symmetry operations these that keep both crystals simultaneously invariant and those that exchange the crystals: there is an homomorphism between the group of the bicrystal and the permutation group of two elements;

* the second section shows how to generalize the bicrystal concept to the N-crystal one by constructing all sets of symmetry operations, if any, that realize one of the permutations of N objects;

* the third and last section proposes a direct application of the notion of N-crystal in the case of the group-subgroup transformations.

To make the reading easy, the following notations are used all over the paper:

* point groups are designated by uppercase letters as G or H$_{α}$; point symmetry elements are designated by lowercase letters as g or α;

* space groups, including translation groups, are noted by calligraphic letters like G or U$_{α}$ and space symmetry elements by $\hat{g}$ or $\hat{\alpha}$;

* for simplicity, the term lattice is indiscriminately used as a set of vectors defining the lattice nodes λ = na + mb + pc, n, m, p ∈ Z or as a set of symmetry elements (1|λ) belonging to a translation group;

* the space group of the crystal is G with point group (symmetry class) G and lattice Λ.

## 2. Standard Bicrystallography

The crystal is characterized by its space group G of lattice Λ as defined in the International Tables for Crystallography [10,11]. The space group is represented by the set of symmetry elements $\hat{g} = (g|t)$ defined by a point symmetry operation g (rotation, mirror, inversion, ...) associated with a translation t:

$$ \hat{g} r = (g|t) r = g r + t $$

The elements of G with the identity as point symmetry operation define the lattice $Λ = \{(1|λ) ∈ G\}$.

The geometric operation that transforms crystal I into crystal II is an isometric transformation that takes the same form as the symmetry elements of the crystal $\hat{\alpha} = (α|τ)$ where α is the point operation (rotation, mirror, inversion, ...) and τ a rigid-body translation: $\hat{\alpha}r = (α|τ)r = αr + τ$.

Because of the intrinsic internal symmetries of the crystal, a generic point $r_1$ has infinitely many equivalent points $Gr_1$, called its G-orbit, that can as well be chosen for defining $\hat{\alpha}$; the same applies in the second crystal. Thus, the chain of equivalences:

$$ r_1 \xrightarrow{G} G r_1 \xrightarrow{\hat{\alpha}} \hat{\alpha} G r_1 \xrightarrow{\hat{\alpha}G\hat{\alpha}^{-1}} \hat{\alpha} G \hat{\alpha}^{-1} \hat{\alpha} G r_1 \equiv \hat{\alpha} G r_1 $$

shows that the most general transformation from crystal I to II, is equivalently described by any element of the left coset $\hat{\alpha}G$ (see for instance Guymont et al. [9]):
$$r_I \xrightarrow{\hat{\alpha}\mathcal{G}} r_{II} = \hat{\alpha}\mathcal{G}r_I$$

The inverse transformation from II to I is characterized by the right coset $\mathcal{G}\hat{\alpha}^{-1}$ as shown in Figure 1.

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Passing from crystal I of space group $\mathcal{G}$ to crystal II of space group $\hat{\alpha}\mathcal{G}\hat{\alpha}^{-1}$ is achieved using any element of the left coset $\hat{\alpha}\mathcal{G}$; the inverse transformation from crystal II to I is described by using any element of the right coset $\mathcal{G}\hat{\alpha}^{-1}$.</figcaption></figure>

## 2.1. Translation Invariants

For a generic interface, the lattices $\Lambda$ and $\Lambda'$ of respectively crystal I and II have no translation in common in the sense that no translation of crystal II can be expressed as a linear combination with integer coefficients of the basic vectors of crystal I.

There are specific cases however, like coincidence grain boundaries and twins, where the two crystals share a same translation subset of rank 3 of $\Lambda$, called the coincidence lattice (see [12]) and noted $\mathcal{T}_{\alpha}$, which is defined by:

$$\mathcal{T}_{\alpha} = \{(1|T) \in \Lambda : \exists (1|T') \in \Lambda \text{ such that } T' = \alpha T\} \quad (1)$$

This coincidence lattice $\mathcal{T}_{\alpha}$ is a subgroup of order $\Sigma$ of $\Lambda$ and $\Lambda' = \alpha \Lambda \alpha^{-1}$ and is better defined as the intersection of these two lattices:

$$\mathcal{T}_{\alpha} = \Lambda \cap \alpha \Lambda \alpha^{-1}$$

The second important translation group is the group $\mathcal{U}_{\alpha}$ generated by the union:

$$\mathcal{U}_{\alpha} = \Lambda \cup \alpha \Lambda \alpha^{-1}$$

This group is a supergroup of order $\Sigma$ of both $\Lambda$ and $\Lambda'$ and is also independent of the rigid-body translation. It has often been called the *Displacement Shift Complete* lattice but is now better defined as *Displacement Symmetry Conserving* lattice (see [5]) with the same acronym DSC.

The translation group $\mathcal{U}_{\alpha}$ is the invariance translation group of the rigid-body translation $\tau$. The rigid-body translation is indeed defined up to any translation of crystal I and any translation of crystal II: for a given point operation $\alpha$, any rigid-body translation $\tau'$ deduced from $\tau$ by a translation of $\mathcal{U}_{\alpha}$ leads to an equivalent bicrystal.

The translation groups $\mathcal{T}_{\alpha}$ and $\mathcal{U}_{\alpha}$ have their image in reciprocal space as demonstrated by Grimmer [13]: the intersection of the reciprocal lattices $\Lambda^*$ and $\Lambda'^* = \alpha \Lambda^* \alpha^{-1}$ is the reciprocal lattice of $\mathcal{U}_{\alpha}: \mathcal{U}_{\alpha}^* = \Lambda^* \cap \alpha \Lambda^* \alpha^{-1}$ whereas the union group is the reciprocal lattice of $\mathcal{T}_{\alpha}$. These general group relations are shown in Figure 2.
$$ \begin{array}{ccc} \mathcal{U} = \Lambda \cup \Lambda' & \mathcal{T}^* = \Lambda^* \cup \Lambda'^* & \mathcal{T} = \Lambda \cap \Lambda' \\ \begin{array}{c} | \\ \Sigma \\ (\Lambda, \Lambda') \end{array} & \begin{array}{c} | \\ \Sigma \\ (\Lambda^*, \Lambda'^*) \end{array} & \begin{array}{c} | \\ \Sigma \\ \mathcal{U}^* = \Lambda^* \cap \Lambda'^* \end{array} \end{array} $$

<caption>**Figure 2.** Tree of the partial order of the translation groups $\mathcal{T}_\alpha$, $\Lambda, \Lambda' = \alpha \Lambda \alpha^{-1}$ and $\mathcal{U}_\alpha$ in direct and reciprocal spaces (see [13]). The order sequence goes from $\mathcal{T}_\alpha$ being a subgroup of $\Lambda$ ($\Lambda'$) of order $\Sigma$ and $\Lambda$ ($\Lambda'$) being a subgroup of $\mathcal{U}_\alpha$ of same order. This corresponds in reciprocal space to exchanging the roles of $\mathcal{T}_\alpha$ and $\mathcal{U}_\alpha$: $\mathcal{T}_\alpha^*$ is the union of the reciprocal lattices and $\mathcal{U}_\alpha^*$ is the intersection of the reciprocal lattices.</caption>

## 2.2. Symmetry and Rigid-Body Translation

Once a specific operation $\alpha$ generating a coincidence lattice is given, the question arises of which kind of space groups $\mathcal{P}_\tau$ can be found depending on the value of the associated rigid-body translation $\tau$. For that purpose, we will first determine the space group $\mathcal{W}_\alpha$ that generates the orbit of all equivalent of $\tau$ whatever its value: as will be discussed next, the little group of $\tau$ in $\mathcal{W}_\alpha$ will thus generate the space group $\mathcal{P}_\tau$.

We have seen in the preceding section that the invariant translation group of $\tau$, whatever its value, is the union group $\mathcal{U}_\alpha$ since $\tau$ is defined up to any translation of either crystals I or II. Concerning the orientational symmetries and since $\tau$ is a global translation of one of the crystals with respect to the other, only the point symmetry operations must be considered.

We use the general rule that *determining the global symmetry of the union of two equivalent objects (their normalizer) is made in searching for the symmetry elements that are common to both objects and the possible additional external elements, if any, that exchange these two objects.*

Let us designate by $G$ the point group of $\mathcal{G}$ and let $\alpha G$ be the (point) coset representing the transformation of I to II. We designate by:

$$ H_{\alpha} = G \cap \alpha G \alpha^{-1} \quad \text{or} \quad H(\alpha) = \{g \in G \text{ such that } \exists g' \in G \text{ with } \alpha g = g' \alpha\} \quad (2) $$

the point group intersection of the points groups of crystals I and II. This point group, keeping both I and II invariant, is part of the construction of $\mathcal{W}_\alpha(\alpha)$.

On the other hand, the elements common to $\alpha G$ and $G\alpha^{-1}$,

$$ E_{\alpha} = \alpha G \cap G\alpha^{-1} \quad \text{or} \quad E(\alpha) = \{\alpha g \in \alpha G \text{ such that } \exists g' \in G \text{ with } \alpha g = g' \alpha^{-1}\} \quad (3) $$

transform crystal I into II and vice-versa, generate other equivalent of $\tau$, but in flipping it in the opposite direction $\tau \to -\tau$. This set, if not empty, is also part of the construction of $\mathcal{W}_\alpha$.

The union:

$$ W_{\alpha} = H_{\alpha} \cup E_{\alpha} $$

is a (point) group with $H_\alpha \subsetneq W_\alpha$:

i. the right and left products of an element of $H_\alpha$ with an element of $E_\alpha$ is an element of $E_\alpha$; let $h_1 \in H_\alpha, \alpha h_2 \in E_\alpha$, then:

$$ h_1(\alpha h_2) = h_1 h_2' \alpha^{-1} \in G\alpha^{-1} \quad \text{and} \quad h_1(\alpha h_2) = \alpha h_1' h_2 \in \alpha G $$

$$ (\alpha h_2)h_1 = h_2' \alpha^{-1} \alpha h_1' \hat{\alpha}^{-1} \in G\alpha^{-1} \quad \text{and} \quad (\alpha h_2)h_1 = \alpha h_2 h_1 \in \alpha G $$
ii. the product of two elements of $E_\alpha$ is an element of $H_\alpha$; let $\alpha h_1, \alpha h_2 \in E_\alpha$, then:

$$(\alpha h_1)(\alpha h_2) = \alpha h_1 h_2' \alpha^{-1} \in \alpha G \alpha^{-1} \text{ and } (\alpha h_1)(\alpha h_2) = h_1' \alpha^{-1} \alpha h_2 = h_1' h_2 \in G$$

Therefore, the group $W_\alpha$ is either $H_\alpha$ itself if $E_\alpha$ is empty, or a supergroup of order 2 of $H_\alpha$ if $E_\alpha$ is not empty; in that later case, it can be written as:

$$W_\alpha = H_\alpha \cup \epsilon H_\alpha$$

with $\epsilon \in E_\alpha$.

It is now trivial to find the group $W_\alpha$ describing all the possible symmetries of the bicrystals generated by a given point operation $\alpha$ according to the rigid-body translation $\tau$; it is given by the direct product of $W_\alpha$ with $U_\alpha$:

$$W_\alpha = (H_\alpha \cup \epsilon H_\alpha) \times U_\alpha.$$

This space group $W_\alpha$ contains at once all the symmetry information for any value of $\tau$. It is enough to restrict $\tau$ inside the Dirichlet domain of $W_\alpha$ to find all the possible space groups $P_\tau$ depending on $\tau$ using the following properties:

* the number of different possible types of space groups $P_\tau$ of the bicrystal created by a point operator $\alpha$ and a rigid-body translation $\tau$ is the number of different position strata (see [14,15]) of the space group $W_\alpha$ that $\tau$ can take;

* the space group $P_\tau$ of the bicrystal is determined by the point group of the symmetry stratum of $\tau$ in $W_\alpha$; in particular, high-symmetry groups appear when $\tau$ points on special positions (strata of dimension 0) of $W_\alpha$; they correspond to the so-called symmetry dictated extrema as discussed by Cahn & Kalonji [16];

* in the construction of $P_\tau$, the generating elements belonging to the intersection group $H_\alpha$ are to be taken as they are, whereas those corresponding to the coset $\epsilon H_\alpha$ must be multiplied by the inversion before being injected in $P_\tau$ (for example, a mirror in the coset $\epsilon H_\alpha$ generates a two-fold axis perpendicular to the mirror as generator of $P_\tau$).

### 2.3. The Space Group of the Bicrystal

Calculating $P_\tau$ is similar to the previous derivation in using space groups instead of point groups. Here again, the space group $P_\tau$ is constituted of two kinds of space symmetry elements:

* those that leave invariant simultaneously I and II: they form a subgroup of $G$ noted $I_\alpha$ and defined by $I_\alpha = G \cap \hat{\alpha} G \hat{\alpha}^{-1}$;

* those that exchange I and II: $E_\alpha = \hat{\alpha} G \cap G \hat{\alpha}^{-1}$ that are outside $G$ and can reduce to the empty set.

The elements (For example, the subgroup $G_\alpha$ of $G$ defined by the elements of $G$ that commute with $\hat{\alpha}$, $G_\alpha = \{\hat{g} \in G \text{ such that } \hat{g} \hat{\alpha} = \hat{\alpha} \hat{g}\}$ is a subgroup of $I_\alpha$. This is the case when $\hat{\alpha}$ is a pure rotation sharing the same axis as a pure rotation $\hat{g}$ of $G$, then $\hat{g}$ is in $I_\alpha$) of $I_\alpha$ are defined by:

$$I_\alpha = \{\hat{g} \in G : \exists \hat{g}' \in G \text{ such that } \hat{g}' = \hat{\alpha} \hat{g} \hat{\alpha}^{-1} \text{ or } \hat{\alpha} \hat{g} = \hat{g}' \hat{\alpha}\}$$

The elements (For example, the binary elements of $\hat{\alpha} G$ modulo $G$ are elements of $E_\alpha$; let $(\hat{\alpha} \hat{g}_1)(\hat{\alpha} \hat{g}_1) = \hat{g}_2$, then:

$$\hat{\alpha} \hat{g}_1 = (\hat{\alpha} \hat{g}_1)(\hat{\alpha} \hat{g}_1) \hat{g}_1^{-1} \hat{\alpha}^{-1} = \hat{g}_2 \hat{g}_1^{-1} \hat{\alpha}^{-1} \in G \hat{\alpha}^{-1}.$$
As an example, the product $\hat{\alpha} \hat{m}$ of a rotation of angle $\alpha$ around an axis contained in the plane of a mirror $\hat{m}$ element of $\mathcal{G}$ is an element of $\mathcal{E}_{\alpha}$: it generates a mirror of $\mathcal{P}_{\tau}$ the plane of which is the mirror plane of $\hat{m}$ rotated by $\alpha/2$) of $\mathcal{E}_{\alpha}$ — if not empty — are defined by:

$$ \mathcal{E}_{\alpha} = \{ \hat{\alpha}\hat{g}, \hat{g} \in \mathcal{G} : \exists \hat{g}' \in \mathcal{G} \text{ such that } \hat{\alpha}\hat{g} = \hat{g}'\hat{\alpha}^{-1} \} $$

As in the preceding section, it is easily demonstrated that the union of the two sets $\mathcal{I}_{\alpha}$ and $\mathcal{E}_{\alpha}$ (if not empty) forms a group $\mathcal{P}_{\tau}$:

$$ \mathcal{P}_{\tau} = \mathcal{I}_{\alpha} \cup \mathcal{E}_{\alpha} = \mathcal{I}_{\alpha} \cup \hat{\mathcal{E}}_{\alpha}, \quad \hat{\mathcal{E}} \in \mathcal{E}_{\alpha}, \mathcal{I}_{\alpha} \triangleleft \mathcal{P}_{\tau} \qquad (4) $$

This space group of the bicrystal with translation group $\mathcal{I}_{\alpha}$ is also designated as the symmetry group of the *dichromatic pattern* (see [17]).

The explicit computation of $\mathcal{P}_{\tau}$ is made particularly simple in considering the little group of $\tau$ in $\mathcal{W}_{\alpha}$:

* the elements of $\mathcal{I}_{\alpha}$ are such that:

$$ (g|t)(\alpha|\tau) = (\alpha|\tau)(g'|t') \text{ or } (g\alpha|g\tau + t) = (\alpha g'|\alpha t' + \tau) $$

leading to $g\alpha = \alpha g'$ and $g\tau - \tau = \alpha t' - t$; this is achieved for $g$ belonging to the little group of $\tau$ ($g\tau = \tau$) where $t$ and $t'$ are translations of $\mathcal{I}_{\alpha}$ ($\alpha t' = t$);

* the elements of $\mathcal{E}_{\alpha}$ are such that:

$$ (g|t)(\alpha|\tau)^{-1} = (\alpha|\tau)(g'|t') \text{ or } (g\alpha^{-1}|t - g\alpha^{-1}\tau) = (\alpha g'|\alpha t' + \tau) $$

leading to $g\alpha^{-1} = \alpha g'$ and $g\alpha^{-1}\tau + \tau = -\alpha t' + t$; this is achieved for $g\alpha^{-1}$ being element of the little group of $\tau$ multiplied by the inversion ($g\alpha^{-1}\tau = -\tau$) where $t$ and $t'$ are translations of $\mathcal{I}_{\alpha}$ ($\alpha t' = t$).

This shows that the little (point) group of $\tau$ in $\mathcal{W}_{\alpha}$, say $W_{\tau} = H_{\tau} \cup \epsilon_{\tau}H_{\tau}$ is directly connected to the point group $P_{\tau}$ of the space group $\mathcal{P}_{\tau}$ that is built with $H_{\tau}$ and $\epsilon_{\tau}H_{\tau}$ multiplied by the inversion:

$$ P_{\tau} = H_{\tau} \cup \bar{1}\epsilon_{\tau}H_{\tau} \qquad (5) $$

The number of different strata of $\mathcal{W}_{\alpha}$ gives the number of various kind of space groups $\mathcal{P}_{\tau}$ depending on the values of $\tau$. In particular, the symmetry dictated values of $\tau$ are those corresponding to strata of dimension 0 of $\mathcal{W}_{\alpha}$.

## 2.4. A Simple 2D Example

To exemplify these rules, we consider the classical case of the bicrystal generated by the rotation of a 2D square structure (see for instance [18–21]) onto itself as shown in Figure 3. The actual structure belongs to either symmetry class G = 4 (space group $\mathcal{G}$ = p4) or G = 4mm that is the holohedry of the square system (space groups $\mathcal{G}$ = p4mm and p4gm). Coincidence lattices appear for rotations $\alpha = 2 \arctan(m/n)$ around the origin that superimpose the node $(n, -m)$ on top of the node $(n, m)$, with n, m coprime and $0 < m < n \in \mathbb{Z}$. This leads to the coincidence lattice defined by the vectors $T_1 = (n, m)$, $T_2 = (-m, n)$; the index is $\Sigma = n^2 + m^2$ and the union lattice $\mathcal{U}_{\alpha}$ is defined by $U_1 = T_1/\Sigma, U_2 = T_2/\Sigma$.
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Example of two square crystals of space groups p4mm (Λ) rotated by α = 53.1301°, drawn in red and blue, leading to a coincidence lattice T₁ = (2, 1), T₂ = (-1, 2), Σ = 5, in black with same square symmetry p4mm (T₁, T₂). The space group W_α = p4mm (U₁, U₂) is shown in turquoise with U₁ = T₁/5 and U₂ = T₂/5. The number of various possible space groups of the (red-blue) bicrystals according to the value of the rigid-body translation are the number of position strata of p4mm (U₁, U₂) as given in Table 1.</figcaption></figure>

The computation of W_α(α) is quite simple. The intersection H_α is the set

$$H(\alpha) = \{g \in G \text{ such that } g\alpha = \alpha g\} = 4$$

and the exchange set is the set of the mirrors in αG ∩ Gα⁻¹ corresponding to the initial mirrors, mₓ(mᵧ), mₓᵧ(mₓŷ) rotated by α/2 for the class 4m. This set is empty for the case of the class 4. Therefore:

$$\text{Class 4mm : } W_{\alpha}(\alpha) = p4mm(U_1, U_2), \quad \text{Class 4 : } W_{\alpha}(\alpha) = p4(U_1, U_2).$$

This demonstrated that whatever the value of the rotation α generating a coincidence lattice, the bicrystal of a p4 structure can have only four space groups: p4(×2), p2 and p1 corresponding to the four strata of the group p4, whereas structures of space groups p4mm and p4gm have seven possible space groups from p4mm to p1 corresponding to the seven strata of the group p4mm as shown on Table 1 and exemplified in Figure 4 for the values of τ corresponding to the special points.

Here, in the holohedral examples, all 4-fold and 2-fold rotations belong to H_α whereas the exchange elements are the original mirrors rotated by α/2. It is interesting to observe in the case of the strata of dimension 1 in Table 1, that, as expected, the actual mirrors of P_τ are perpendicular to those of the little group of τ as expected in relation (5); for instance a little group mₓ generates a space group pmᵧ or pgᵧ and a little group mₓᵧ generates a mirror mₓŷ.
<table><caption>**Table 1.** The complete set of symmetry groups $\mathcal{P}_\tau$ for general square structures of space groups $p4mm$, $p4gm$ and $p4$ as function of the rigid-body translation $\tau$. Structures of groups $p4mm$ and $p4gm$ share the same $\mathcal{W}_\alpha = p4mm$ and structures of group $p4$, generating no exchange set, have $\mathcal{W}_\alpha = p4$. The coordinates of the translation $\tau$ are given with respect to the unit cell of the union group $\mathcal{W}_\alpha$ defined by $U_1 = (n,m)/\Sigma$ and $U_2 = (-m,n)/\Sigma$ with $\Sigma = n^2 + m^2$. The (primitive) unit cell of the groups $\mathcal{P}_\tau$ are defined by $T_1 = (n,m)$ and $T_2 = (-m,n)$.</caption><thead><tr><th>τ</th><th>Stratum Dim</th><th>Little Group in W<sub>α</sub></th><th>P<sub>τ</sub> for p4mm</th><th>P<sub>τ</sub> for p4gm</th><th>P<sub>τ</sub> for p4</th></tr></thead><tbody><tr><td>(0,0)</td><td>0</td><td>4mm 4</td><td>p4mm</td><td>p4gm</td><td>p4</td></tr><tr><td>(1/2, 1/2)</td><td>0</td><td>4mm 4</td><td>p4gm</td><td>p4mm</td><td>p4</td></tr><tr><td>(1/2, 0)</td><td>0</td><td>2mm 2</td><td>p2mg</td><td>p2gm</td><td>p2</td></tr><tr><td>(x, y)</td><td>2</td><td>1</td><td>p1</td><td>p1</td><td>p1</td></tr><tr><td>(x, 0)</td><td>1</td><td>.m<sub>x</sub>.</td><td>pm<sub>y</sub></td><td>pg<sub>y</sub></td><td>—</td></tr><tr><td>(x, 1/2)</td><td>1</td><td>.m<sub>x</sub>.</td><td>pg<sub>y</sub></td><td>pm<sub>y</sub></td><td>—</td></tr><tr><td>(x, x)</td><td>1</td><td>.m<sub>x&#773;y</sub></td><td>cm<sub>xy</sub></td><td>cm<sub>xy</sub></td><td>—</td></tr></tbody></table>

<figure><img src="image_5.png" /><figcaption>**Figure 4.** The possible space groups $\mathcal{P}_\tau$ given in Table 1 for the bicrystal generated by the rotation $\alpha = 2 \arctan(m/n)$ (here for $n = 2, m = 1, \Sigma = 5, \alpha = 53.1301^\circ$) as a function of the rigid-body translation $\tau$ at the special points $(0,0)$, $(1/2, 1/2)$ and $(1/2, 0)$ for the three square symmetry space groups $p4mm$, $p4gm$ and $p4$. This scheme holds for any integer values of $n$ and $m$.</figcaption></figure>
## 2.5. The Special Case of General Quasi-Bicrystals

The physical meaning of the rigid-body translation $\tau$ deserves special attention for the case of general crystal superimposition with no coincidence lattice. This case corresponds to the limit $\Sigma \to \infty$ where the bicrystal becomes a quasiperiodic structure (see for instance [22,23]). The translation group $U_\alpha$ becomes thus a $\mathbb{Z}$-module and the size of the unit cell of $U_\alpha$ reduces to zero. We call it a quasi-bicrystal. This means that the general quasi-bicrystal is then independent of $\tau$: any two quasi-bicrystals built with the same I and II crystals that differ only by a change $\delta\tau$ in the rigid-body translation are locally isomorphic, i.e., any local configuration of finite size in one quasi-bicrystal is present in the other with the same frequency and vice-versa. The two quasi-bicrystals are thus undistinguishable in the sense that modifying the rigid-body translation between the two crystals has the only effect of redistributing in space the atomic configurations of finite size with the same frequencies without adding any new configurations nor removing any existing ones.

# 3. From Bicrystallography to N-Crystallography

We can extend the bicrystal concept to that of N-crystal where $N > 2$ crystals are superimposed in observing that, for $\mathcal{E}$ non-empty, the quotient group $\mathcal{P} / \mathcal{I}$ is isomorphic to the permutation group of two objects $\mathcal{S}_2$:

$$ \mathcal{P} / \mathcal{I} = \mathcal{S}_2. \qquad (6) $$

This basic isomorphism is the starting point for making an easy generalization to N-crystal.

## 3.1. The N-Crystal Generalization

The overall symmetry group of N equivalent objects for finite N, is isomorphic to a subgroup of the permutation group $\mathcal{S}_N$ of N objects. Thus, constructing the symmetry of the abstract N-crystal, consists of identifying which actual symmetry elements are representative of those of $\mathcal{S}_N$.

For notation coherency, we note the identity as $\hat{\alpha}_1 = \text{Id}$. The elements of $\mathcal{S}_N$, say $\pi_k$, transform the N-uplet of integers $(1, 2, \dots, N)$ into $(\pi_k(1), \pi_k(2), \dots, \pi_k(N))$. On the other hand, transforming crystal i into j is achieved by the elements of the coset (if not empty) $\hat{\alpha}_j \mathcal{G} \hat{\alpha}_i^{-1}$ as illustrated in Figure 5.

<figure><img src="image_6.png" /><figcaption>**Figure 5.** Transforming variant i into j in a N-crystal is achieved by transforming i back into the reference variant noted 1 using $\mathcal{G}\hat{\alpha}_i^{-1}$ and then into variant j by applying $\hat{\alpha}_j\mathcal{G}$. This leads to the global operation $\hat{\alpha}_j \mathcal{G} \hat{\alpha}_i^{-1}$. Then, the transposition $(i, j) \to (j, i)$ has the image $(\hat{\alpha}_j \mathcal{G} \hat{\alpha}_i^{-1}) \cap (\hat{\alpha}_i \mathcal{G} \hat{\alpha}_j^{-1})$ (possibly empty) in the standard crystallographic representation.</figcaption></figure>
At the substitution $j \to p_k(j)$, corresponds the coset $\hat{\alpha}_{p_k(j)} \mathscr{G} \hat{\alpha}_j^{-1}$. Thus, at the permutation $p_k$ of $N$ objects:

$$ (1, 2, \dots, j, \dots, N) \to (p_k(1), p_k(2), \dots, p_k(j), \dots, p_k(N)) $$

corresponds the intersection (if not empty) of the cosets $\hat{\alpha}_{p_k(j)} \mathscr{G} \hat{\alpha}_j^{-1}$ for all $j$ varying from 1 to $N$:

$$ (p_k(1), p_k(2), \dots, p_k(j), \dots, p_k(N)) \implies \bigcap_{j=1}^{N} \hat{\alpha}_{p_k(j)} \mathscr{G} \hat{\alpha}_j^{-1} $$

The symmetry group $\mathscr{P}$ of the $N$-crystal is the union of the $N!$ such intersections (many of them can be empty) associated with the $N!$ permutations of the group $\mathscr{S}_n$:

$$ \mathscr{P} = \bigcup_{k=1}^{N!} (\bigcap_{j=1}^{N} \hat{\alpha}_{p_k(j)} \mathscr{G} \hat{\alpha}_j^{-1}) \quad (7) $$

The mapping between each element of the permutation group and its corresponding coset intersections (possibly empty) insures $\mathscr{P}$ to be a group.

Introducing the invariant subgroup $\mathscr{I} \subset \mathscr{P}$ generated by the identity permutation $p_1(j) = j$:

$$ \mathscr{I} = \bigcap_{j=1}^{N} \hat{\alpha}_j \mathscr{G} \hat{\alpha}_j^{-1}, $$

we find the basic decomposition for the $N$-crystal symmetry group:

$$ \mathscr{P} = \bigcup_{k=1}^{N!} \varepsilon_k \mathscr{I}, \quad \text{with } \varepsilon_1 = \mathrm{Id} \quad (8) $$

and

$$ \mathscr{P} / \mathscr{I} \simeq \sigma_N \leq \mathscr{S}_N $$

where $\sigma_N$ is a subgroup of $\mathscr{S}_N$.

Expressions (7) and (8) are the explicit formula that define the symmetry group of the $N$-crystal as natural extensions of formula (4) of the bicrystal.

## 3.2. Application to Group–Subgroup Phase Transformations

The simplest and most obvious application of using relation (8) is the classical group–subgroup phase transformations as exemplified by the numerous cases encountered in solid state physics of alloys with the numerous examples of order–disorder transformations. These are characterized by a symmetry breaking at some critical temperature—or on reaching a two-phased region in the phase diagram—from a high-symmetry group $\mathscr{G}$ at high temperature to a low-symmetry subgroup $\mathscr{H}$ of order $N$ in $\mathscr{G}$ at low temperature:

$$ \mathscr{G} \to [\mathscr{G} : \mathscr{H}], \quad \mathscr{H} < \mathscr{G}, \quad \mathscr{G} = \bigcup_{i=1}^{N} \hat{\mathscr{g}}_i \mathscr{H}, \quad \hat{\mathscr{g}}_1 = \mathrm{Id} $$

The decomposition of $\mathscr{G}$ in cosets of $\mathscr{H}$ corresponds to a single crystal of space group $\mathscr{G}$ at high temperature transforming, at low temperature, into numerous crystallites of phase $\mathscr{H}$ that distribute into $N$ kind of differently oriented/translated crystals $\mathscr{H}$. Each family, called a variant, gather all crystallites described by the same space group $\mathscr{H}$. The coset $\hat{\mathscr{g}}_i \mathscr{H}$ is the set of all equivalent operations that transform a given chosen variant of space group $\mathscr{H}$ into another of space group $\hat{\mathscr{g}}_i \mathscr{H} \hat{\mathscr{g}}_i^{-1}$. Passing from the variant $i$ to the variant $j$ is achieved by the set of operations $\hat{\mathscr{g}}_j \mathscr{H} \hat{\mathscr{g}}_i^{-1}$. We recognize here the expressions discussed in the previous section where the operations $\hat{\alpha}_i$ are identified to $\hat{\mathscr{g}}_i$. Therefore, the symmetry group $\mathscr{P}_\tau$ is the group $\mathscr{G}$ itself since the $N$-crystal is the high temperature phase:

$$ \mathscr{G} = \mathscr{P} = \bigcup_{k=1}^{N!} (\bigcap_{j=1}^{N} \hat{\mathscr{g}}_{p_k(j)} \mathscr{H} \hat{\mathscr{g}}_{j}^{-1}), \quad \hat{\mathscr{g}}_{1} = \mathrm{Id} \quad (9) $$
The group $\mathcal{G}$ of the high temperature phase describes the way the $N$ variants of the low temperature phase transform into each other. We designate by $\mathcal{I}$ the invariant subgroup of order $K$ in $\mathcal{H}$ that is the intersection of the groups of the $N$ variants:

$$ \mathcal{I} = \bigcap_{i=1}^{N} \hat{\mathbf{g}}_i \mathcal{H} \hat{\mathbf{g}}_i^{-1} $$

This intersection group is the kernel of the application $\mathcal{G} \to [\mathcal{G} : \mathcal{H}]$ of $\mathcal{G}$ on its coset decomposition on $\mathcal{H}$. The pertinent group to be considered for discussing the symmetry breaking induced by the transformation is the image of the application. It is associated with the symmetry group $\sigma_N$ of order $N \times K$, isomorphic to a subgroup of $\mathcal{I}_N$ obtained by computing the coset decomposition $\mathcal{G} = \mathcal{P}$ of $\mathcal{I}_\alpha$ and extracting its quotient by $\mathcal{I}_\alpha$:

$$ \mathcal{G}(=\mathcal{P}) = \bigcup_{j=1}^{N \times K} \varepsilon_j \mathcal{I}_\alpha, \quad \sigma_n \simeq \mathcal{G} / \mathcal{I} $$

This image group $\sigma_N \simeq \mathcal{G} / \mathcal{I}$ gives the complete and faithful description of how the parent group $\mathcal{G}$ operates on the variants $\mathcal{H}$.

### 3.3. Simple Examples

#### 3.3.1. A 2D Toy-Model

We consider the hypothetical phase transition of blue atoms decomposing by chemical ordering into red, green and brown atoms as shown in Figure 6. This order-disorder transformation is characterized by the group-subgroup relation:

$$ \mathcal{G} = p3m1(a,b) \to \mathcal{H} = cm1(a,c = a + 2b) $$

with standard coset decomposition $[\mathcal{G} : \mathcal{H}]$:

$$ p3m1(a,b) = \{\hat{\mathbf{1}} + \hat{\mathbf{3}}^1 + \hat{\mathbf{3}}^2\}cm1(a,c) $$

corresponding to generate 3 variants deduced from each other by the 3-fold rotation lost during the transition.

<figure><img src="image_7.png" /><figcaption>**Figure 6.** Group-subgroup phase transition from (a) $p3m1(a,b)$ to (b) $cm1(a,c = a + 2b)$ generating three variants (c) deduced from each other by the permutation group $\mathcal{I}_3$ isomorphic to $3m$.</figcaption></figure>
We note $\hat{\mathbf{g}}_1 = \hat{\mathbf{1}}, \hat{\mathbf{g}}_2 = \hat{\mathbf{3}}^1 = (3^1|0,0)$ and $\hat{\mathbf{g}}_3 = \hat{\mathbf{g}}_2^{-1} = \hat{\mathbf{3}}^2 = (3^2|0,0)$. The permutation group (equivalent here to a three-color symmetry group) is obtained by computing the quotient group of $\mathcal{G}$ on the kernel $\mathcal{I}$:

$$ \mathcal{I} = \mathcal{H} \cap \hat{\mathbf{g}}_2 \mathcal{H} \hat{\mathbf{g}}_2^{-1} \cap \hat{\mathbf{g}}_3 \mathcal{H} \hat{\mathbf{g}}_3^{-1} = p1(a,b) $$

There are thus 6 cosets in the decomposition of $\mathcal{G}$ on $\mathcal{I}$ defined by the following variant permutations and corresponding cosets intersections:

$$
\begin{array}{l}
\mathcal{I}_{\alpha}: \quad \hat{\mathbf{1}} \quad \rightarrow \quad (1,2,3) \longleftrightarrow \hat{\mathbf{g}}_1 \mathcal{G} \hat{\mathbf{g}}_1^{-1} \cap \hat{\mathbf{g}}_2 \mathcal{G} \hat{\mathbf{g}}_2^{-1} \cap \hat{\mathbf{g}}_3 \mathcal{G} \hat{\mathbf{g}}_3^{-1} \\
\mathcal{E}_{\alpha}: \quad \left\{
\begin{array}{l}
\hat{\mathbf{3}}^1 \quad \rightarrow \quad (2,3,1) \longleftrightarrow \hat{\mathbf{g}}_2 \mathcal{G} \hat{\mathbf{g}}_1^{-1} \cap \hat{\mathbf{g}}_3 \mathcal{G} \hat{\mathbf{g}}_2^{-1} \cap \hat{\mathbf{g}}_1 \mathcal{G} \hat{\mathbf{g}}_3^{-1} \\
\hat{\mathbf{3}}^2 \quad \rightarrow \quad (3,1,2) \longleftrightarrow \hat{\mathbf{g}}_3 \mathcal{G} \hat{\mathbf{g}}_1^{-1} \cap \hat{\mathbf{g}}_2 \mathcal{G} \hat{\mathbf{g}}_3^{-1} \cap \hat{\mathbf{g}}_1 \mathcal{G} \hat{\mathbf{g}}_2^{-1} \\
\hat{\mathbf{m}}_1 \quad \rightarrow \quad (1,3,2) \longleftrightarrow \hat{\mathbf{g}}_1 \mathcal{G} \hat{\mathbf{g}}_1^{-1} \cap \hat{\mathbf{g}}_3 \mathcal{G} \hat{\mathbf{g}}_2^{-1} \cap \hat{\mathbf{g}}_2 \mathcal{G} \hat{\mathbf{g}}_3^{-1} \\
\hat{\mathbf{m}}_2 \quad \rightarrow \quad (3,2,1) \longleftrightarrow \hat{\mathbf{g}}_3 \mathcal{G} \hat{\mathbf{g}}_1^{-1} \cap \hat{\mathbf{g}}_2 \mathcal{G} \hat{\mathbf{g}}_2^{-1} \cap \hat{\mathbf{g}}_1 \mathcal{G} \hat{\mathbf{g}}_3^{-1} \\
\hat{\mathbf{m}}_3 \quad \rightarrow \quad (2,1,3) \longleftrightarrow \hat{\mathbf{g}}_2 \mathcal{G} \hat{\mathbf{g}}_1^{-1} \cap \hat{\mathbf{g}}_1 \mathcal{G} \hat{\mathbf{g}}_2^{-1} \cap \hat{\mathbf{g}}_3 \mathcal{G} \hat{\mathbf{g}}_3^{-1}
\end{array}
\right.
\end{array}
$$

and the explicit 3-crystal symmetry group can be written as:

$$ \mathcal{G} = \mathcal{P} = \{\hat{\mathbf{1}} + \hat{\mathbf{3}}^1 + \hat{\mathbf{3}}^2 + \hat{\mathbf{m}}_1 + \hat{\mathbf{m}}_2 + \hat{\mathbf{m}}_3\} p1(a,b) $$

the image of which is the group of permutation of three objects of order 6:

$$ \sigma_3 = \{(1,2,3), (2,3,1), (3,1,2), (1,3,2), (3,2,1), (2,1,3)\} $$

isomorphic to the point group 3m. This is shown in Figure 6c where the drawing of the three-fold axis and the three mirrors $m_1, m_2$ and $m_3$ explicitly demonstrates how the three variants connect to each other according to 3m.

A few other classical examples in the family of chemical ordering of metallic alloys are the following.

### 3.3.2. CuZn B$_2$:

$$ \text{Here we have } \mathcal{G} = \text{Im}3m(a), \mathcal{H} = \text{Pm}3m(a), N = 2, \mathcal{I} = \mathcal{H}, \sigma_2 \simeq 2 = \mathcal{I}_2 $$

$$ \text{Im}3m(a) = \{(1|000) + (1|1/2,1/2,1/2)\} \text{Pm}3m(a) $$

Here, we are in the simplest case of bicrystal induced by transformation where the translation $1/2(1,1,1)$ is the generator of the exchange set (color translation). The image group is $\mathcal{I}_2$.

### 3.3.3. Cu$_3$Au L1$_2$:

$$ \mathcal{G} = \text{Fm}3m(a), \mathcal{H} = \text{Pm}3m(a), N = 4, \mathcal{I} = \text{Pmmm}(a), \sigma_4 \simeq m3 = \mathcal{I}_4 $$

$$ \text{Fm}3m(a) = \{(1|000) + (1|1/2,1/2,0) + (1|0,1/2,1/2) + (1|1/2,0,1/2)\} \\ \times \{(1|000) + (C_3|0,0,0) + (C_3^2|0,0,0)\} \times \{(1|000) + (m_{x,y}|000)\} \text{Pmmm}(a) $$

Here the image group is isomorphic to m3 (regular tetrahedron). We can represent the four variants as the vertices of a regular tetrahedron. For example applying the translation $1/2(1,1,0)$ to the four variants corresponds to the permutation $(1,2,3,4) \rightarrow (2,1,4,3)$.

## 4. Conclusions

We have shown here that the symmetry group of the superimposition of two crystals related by the boundary operation $\hat{\mathbf{a}} = (a|\tau)$ is the combination of two kind of operations: those that
leave simultaneous invariant both crystals ($\mathcal{I}_\alpha = \mathcal{G} \cap \hat{\alpha}\mathcal{G}\hat{\alpha}^{-1}$) and those that exchange them ($\varepsilon\mathcal{I}_\alpha = \hat{\alpha}\mathcal{G} \cap \mathcal{G}\hat{\alpha}^{-1}$):

$$ \mathcal{P} = \mathcal{I}_\alpha \cup \varepsilon\mathcal{I}_\alpha $$

and, when the exchange set $\varepsilon\mathcal{I}_\alpha$ is non-empty, the quotient group $\mathcal{P}/\mathcal{I}_\alpha$ is isomorphic to $\mathcal{S}_2$, the permutation group of two objects.

The space group $\mathcal{P}$ depends on both the point operation $\alpha$ and the translation part $\tau$ of $\hat{\alpha}$. The number of different possible space groups associated with a given point operation $\alpha$ according to the $\tau$ value is equal to the number of different strata of the group:

$$ \mathcal{W}_\alpha = [(G \cup \alpha G\alpha^{-1}) \cap (\alpha G \cup G\alpha^{-1})] \times \mathcal{U}_\alpha $$

where $G$ is the point group of the structure and $\mathcal{U}_\alpha$ the union group of the lattices of crystals I and II.

The symmetry group of the superimposition of a finite number (if $N \to \infty$, as for cubic variants generated by iterative mirror twinning on the [1, 1, 1] plane, the present description fails in which case the analysis developed by Cayron (see for instance [24]) using groupoids that is to be considered), say $N$, of variants is the combination of all possible symmetry operation that permute the $N$ variants and can be written:

$$ \mathcal{P} = \bigcup_{k=1}^{N!} \left( \bigcap_{j=1}^{N} \hat{\alpha}_{p_k(j)} \mathcal{G} \hat{\alpha}_{j}^{-1} \right) $$

where the $N$ structurally identical variants are related by $\hat{\alpha}_j\mathcal{G}$ cosets with $j = 1, N$. Here again, the quotient group $\sigma_N \simeq \mathcal{P}/\mathcal{I}$, where $\mathcal{I} = \bigcap_{j=1}^{N} \hat{\mathbf{g}}_j \mathcal{H} \hat{\mathbf{g}}_j^{-1}$ is a subgroup of the permutation group $\mathcal{I}_N$ of $N$ objects.

This generalization allows the gathering in a unique scheme of the symmetry studies of $N$-crystallography on special grain boundaries of $N$ (finite) variants and the coset decomposition in $N$ variants due to a group-subgroup solid state transformation.

The only problem here was to determine the mapping between the actual symmetry crystallographic elements on one side and the abstract permutation operations on the other side. As always when using group theory, the pertinent mapping is obtained by extracting the image group of the application, say $\varphi$, as the quotient of the overall symmetry group ($\mathcal{P}$ or $\mathcal{G}$) by the kernel of the application given by the intersection group $\mathcal{I}$ that is invariant in the parent group:

$$ \mathcal{G} \xrightarrow{\varphi} [\mathcal{G} : \mathcal{H}], \quad \mathcal{I} = \bigcap_{j=1}^{N} \hat{\mathbf{g}}_j \mathcal{H} \hat{\mathbf{g}}_j^{-1} \\ \mathrm{Im} \, \varphi = \mathcal{G} / \mathcal{I} \simeq \sigma_N \leq \mathcal{I}_N. $$

or equivalently:

$$ \mathcal{G} \xrightarrow{\varphi} \mathcal{P} = \bigcup_{k=1}^{N!} \left( \bigcap_{j=1}^{N} \hat{\alpha}_{p_k(j)} \mathcal{G} \hat{\alpha}_{j}^{-1} \right), \quad \mathcal{I} = \bigcap_{j=1}^{N} \hat{\alpha}_{j} \mathcal{G} \hat{\alpha}_{j}^{-1} \\ \mathrm{Im} \, \varphi = \mathcal{P} / \mathcal{I} \simeq \sigma_N \leq \mathcal{I}_N. $$

In both cases, the mathematical skeleton is essentially the natural concept of permutations between objects: symmetry operations in crystallography do nothing but swapping equivalent objects with each others.

**Author Contributions:** Conceptual derivation, D.G. and M.Q.; writing—original draft preparation, D.G.; writing—review and editing, M.Q. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding beyond the salaries of the authors from CNRS-France.

**Acknowledgments:** The authors are pleased to acknowledge fruitful discussions on bicrystals and dichromatic patterns with Sylvie Lartigue-Korinelk and Richard Portier, and on the general notion of variants with Cyril Cayron.

**Conflicts of Interest:** The authors declare no conflict of interest.
References and Notes

1. De Laissardière, G.T.; Mayou, D.; Magaud, L. Localization of Dirac Electrons in Rotated Graphene Bilayers. *Nano Lett.* **2010** *10*, 804–808. [CrossRef]

2. Bollmann, W. *Crystal Defects and Crystalline Interfaces*: Springer: Berlin/Heidelberg, Germany, 1970.

3. Friedel, G. *Leçons de Cristallographie*; 2nd ed.; Paris, Berger-Lecrault Paris Blanchard: Paris, France, 1964.

4. Friedel, G. Etude sur les Groupements Cristallins. In *Bulletin de la Société de l'Industrie Minérale*; Quatrième Série, Tomes III et IV; Société de l'Imprimerie Thèolier J. Thomas et C: Saint-Etienne, France, 1904.

5. Pond, R.C.; Bollmann, W. The symmetry and interfacial structure of bicrystals. *Philos. Trans. R. Soc. Lond.* **1979**, *292*, 449–472.

6. Pond, R.C.; Vlachavas, D.S. Bicrystallography. *Proc. R. Soc. Lond.* **1983**, *386*, 95–143.

7. Gratias, D.; Portier, R. General geometrical models of grain boundaries. *J. Phys. Colloques* **1982**, *43*, 6–15. [CrossRef]

8. Bourbaki, N. Structures algébriques. In *Eléments de Mathématiques*; Editions Hermann: Paris, France, 1970.

9. Guymont, M.; Gratias, D.; Portier, R.; Fayard, M. Space group theoretical determination of translation, twin, and translation-twin boundaries in cell-preserving phase transitions. *Phys. Stat. Solidi* **1976**, *38*, 629–636. [CrossRef]

10. Hahn, T. *International Tables for Crystallography Volume A: Space-Group Symmetry*; Springer: Berlin/Heidelberg, Germany, 2005.

11. Grimmer, H.; Bollmann, W.; Warrington, H.D. Coincidence Site Lattices and Complete Pattern-Shift Lattices in Cubic Crystals. *Acta Cryst.* **1974**, A30 Pt 2, 197–205. [CrossRef]

12. The usual term coincidence site lattice [11] suggests a set of precise crystallographic sites, which are common to both crystals. In fact, it is a set of translation operations that act the same way in both crystals by transforming any point of either crystal into another equivalent one of the same crystal.

13. Grimmer, H. A reciprocity relation between the coincidence site lattice and the DSC lattice. *Scr. Met.* **1974**, *8*, 1221–1224. [CrossRef]

14. Each position *x* has a little group $G_x$ in $\mathcal{G}$ defined by $G_x = \{g \in \mathcal{G}, x = gx\}$; the set of all positions *x* with conjugated little groups in $\mathcal{G}$ form a stratum also called a Wyckoff [15] position by crystallographers.

15. Wyckoff, R.G. *Crystal Structures*, 2nd ed.; Interscience: New York, NY, USA, 1963.

16. Cahn, J.-W.; Kalonji, G. Symmetry in solid-state transformation morphologies. In Proceedings of the International Conference on Solid [to] Solid Phase Transformations, Pittsburgh, PA, USA, 10–14 August 1981.

17. To make the nature of the symmetry elements explicit, two different colors are usually attributed to each of the crystals: The elements of $\mathcal{I}_\alpha$ are grey elements since they do not change the crystals whereas those of $\mathcal{E}_\alpha$ are color elements since they exchange the two crystals.

18. Grimmer, H. The Generating Function for Coincidence Site Lattices in the Cubic System. *Acta Cryst.* **1984**, *40*, 108–112. [CrossRef]

19. Grimmer, H. Coincidence rotations for cubic lattices. *Scr. Met.* **1973**, *7*, 1295–1300. [CrossRef]

20. Ranganathan, S. On the Geometry of Coincidence-Site Lattices. *Acta Cryst.* **1966**, *21*, 197–199. [CrossRef]

21. Romeu, D.; Gómez, A. Recurrence Properties of O- Lattices and the Classification of Grain Boundaries. *Acta Cryst.* **2006**, *62*, 411–412. [CrossRef] [PubMed]

22. Shechtman, D.; Blech, I.; Gratias, D.; Cahn, J.W. Metallic Phase with Long-Range Orientational Order and No Translational Symmetry *Phys. Rev. Lett.* **1984**, *53*, 1951–1953. [CrossRef]

23. Levine, D.; Steinhardt, P.J. Quasicrystals: A New Class of Ordered Structures. *Phys. Rev. Lett.* **1984**, *53*, 2477–2480. [CrossRef]

24. Cayron, C. Groupoid of orientational variants. *Acta Cryst.* **2006**, A62, 21–40. [CrossRef] [PubMed]

<figure><img src="image_8.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).