View source for Crystallographic restriction theorem

{{Short description|Theorem about admissible crystal symmetries}}
{{Use dmy dates|date=January 2022}}{{Missing information|a succinct statement of the theorem|date=December 2018}}
The '''crystallographic restriction theorem''' in its basic form was based on the observation that the [[rotational symmetry|rotational symmetries]] of a [[crystal]] are usually limited to 2-fold, 3-fold, 4-fold, and 6-fold. However, [[quasicrystal]]s can occur with other diffraction pattern symmetries, such as 5-fold; these were not discovered until 1982 by [[Dan Shechtman]].<ref>Shechtman et al (1982)</ref>

Crystals are modeled as discrete [[lattice (group)|lattice]]s, generated by a list of [[linear independence|independent]] finite [[translation (group theory)|translations]] {{Harv|Coxeter|1989}}. Because discreteness requires that the spacings between lattice points have a lower bound, the [[group (mathematics)|group]] of rotational symmetries of the lattice at any point must be a [[finite group]] (alternatively, the point is the only system allowing for infinite rotational symmetry). The strength of the theorem is that ''not all'' finite groups are compatible with a discrete lattice; in any dimension, we will have only a finite number of compatible groups.

==Dimensions 2 and 3==

The special cases of 2D ([[wallpaper group]]s) and 3D ([[space group]]s) are most heavily used in applications, and they can be treated together.

=== Lattice proof ===

A rotation symmetry in dimension 2 or 3 must move a lattice point to a [[Group action (mathematics)#Orbits and stabilizers|succession]] of other lattice points in the same plane, generating a [[regular polygon]] of coplanar lattice points. We now confine our attention to the plane in which the symmetry acts {{Harv|Scherrer|1946}}, illustrated with lattice [[Vector (geometric)|vector]]s in the figure.

[[Image:Crystallographic restriction polygons.png|right|280px|thumb|Lattices restrict polygons<br />
&nbsp;''Compatible'': 6-fold (3-fold), 4-fold (2-fold)<br />
&nbsp;''Incompatible'': 8-fold, 5-fold]]

Now consider an 8-fold rotation, and the displacement vectors between adjacent points of the polygon. If a displacement exists between any two lattice points, then that same displacement is repeated everywhere in the lattice. So collect all the edge displacements to begin at a single lattice point. The [[edge vector]]s become radial vectors, and their 8-fold symmetry implies a regular octagon of lattice points around the collection point. But this is ''impossible'', because the new octagon is about 80% as large as the original. The significance of the shrinking is that it is unlimited. The same construction can be repeated with the new octagon, and again and again until the distance between lattice points is as small as we like; thus no ''discrete'' lattice can have 8-fold symmetry. The same argument applies to any ''k''-fold rotation, for ''k'' greater than 6.

A shrinking argument also eliminates 5-fold symmetry. Consider a regular pentagon of lattice points. If it exists, then we can take every ''other'' edge displacement and (head-to-tail) assemble a 5-point star, with the last edge returning to the starting point. The vertices of such a star are again vertices of a regular pentagon with 5-fold symmetry, but about 60% smaller than the original.

Thus the theorem is proved.

The existence of quasicrystals and [[Penrose tiling]]s shows that the assumption of a linear translation is necessary.  Penrose tilings may have 5-fold [[rotational symmetry]] and a discrete lattice, and any local neighborhood of the tiling is repeated infinitely many times, but there is no linear translation for the tiling as a whole.  And without the discrete lattice assumption, the above construction not only fails to reach a contradiction, but produces a (non-discrete) counterexample.  Thus 5-fold rotational symmetry cannot be eliminated by an argument missing either of those assumptions.  A Penrose tiling of the whole (infinite) plane can only have exact 5-fold rotational symmetry (of the whole tiling) about a single point, however, whereas the 4-fold and 6-fold lattices have infinitely many centres of rotational symmetry.

=== Trigonometry proof ===

Consider two lattice points A and B separated by a translation vector '''''r'''''. Consider an angle α such that a rotation of angle α about any lattice point is a symmetry of the lattice. Rotating about point B by α maps point A to a new point A'.  Similarly, rotating about point A by α maps B to a point B'.  Since both rotations mentioned are symmetry operations, A' and B' must both be lattice points.  Due to periodicity of the crystal, the new vector '''''r'''''' which connects them must be equal to an integer multiple of '''''r''''':

: <math> \mathbf{r}' = m\mathbf{r} </math>

with <math>m</math> integer. The four translation vectors, three of length <math>r=|\mathbf{r}|</math> and one, connecting A' and B', of length <math>r'=|\mathbf{r}'|</math>, form a trapezium.  Therefore, the length of '''''r'''''' is also given by:

: <math> r' = 2r\cos\alpha - r. </math>

Combining the two equations gives:

: <math> \cos\alpha = \frac{m+1}{2} = \frac{M}{2}  </math>

where <math>M=m+1</math> is also an integer. Bearing in mind that <math>|\cos\alpha|\le 1</math> we have allowed integers <math>M\in\{-2,-1,0,1,2\}</math>. Solving for possible values of <math>\alpha</math> reveals that the only values in the 0° to 180° range are 0°, 60°, 90°, 120°, and 180°.  In radians, the only allowed rotations consistent with lattice periodicity are given by 2π/''n'', where ''n'' = 1, 2, 3, 4, 6.  This corresponds to 1-, 2-, 3-, 4-, and 6-fold symmetry, respectively, and therefore excludes the possibility of 5-fold or greater than 6-fold symmetry.

=== Short trigonometry proof ===
[[Image:Crystallographic restriction.png|thumb|right]]
[[Image:Crystallographic restriction 2.png|thumb|right]]

Consider a line of atoms ''A-O-B'', separated by distance ''a''. Rotate the entire row by θ = +2π/''n'' and θ = −2π/''n'', with point ''O'' kept fixed. After the rotation by +2π/''n'', A is moved to the lattice point ''C'' and after the rotation by -2π/''n'', B is moved to the lattice point ''D''. Due to the assumed periodicity of the lattice, the two lattice points ''C'' and ''D'' will be also in a line directly below the initial row; moreover ''C'' and ''D'' will be separated by ''r'' = ''ma'', with ''m'' an integer. But by trigonometry, the separation between these points is:

: <math>2a\cos{\theta} = 2a\cos{\frac{2\pi}{n}}</math>.

Equating the two relations gives:

: <math>2\cos{\frac{2\pi}{n}}=m</math>

This is satisfied by only ''n'' = 1, 2, 3, 4, 6.

=== Matrix proof ===

For an alternative proof, consider [[matrix (mathematics)|matrix]] properties. The sum of the diagonal elements of a matrix is called the [[trace (matrix)|trace]] of the matrix. In 2D and 3D every rotation is a planar rotation, and the trace is a function of the angle alone. For a 2D rotation, the trace is 2 cos θ; for a 3D rotation, 1 + 2 cos θ.

'''Examples'''
*Consider a 60° (6-fold) [[rotation matrix]] with respect to an [[orthonormal basis]] in 2D.
::<math>\begin{bmatrix} {1/2} & -{\sqrt{3}/2} \\ {\sqrt{3}/2} & {1/2} \end{bmatrix}</math>
:The trace is precisely 1, an [[integer]].
*Consider a 45° (8-fold) rotation matrix.
::<math>\begin{bmatrix} {1/\sqrt{2}} & -{1/\sqrt{2}} \\ {1/\sqrt{2}} & {1/\sqrt{2}} \end{bmatrix}</math>
:The trace is 2/{{radic|2}}, not an integer.

Selecting a basis formed from vectors that spans the lattice, neither orthogonality nor unit length is guaranteed, only linear independence. However the trace of the rotation matrix is the same with respect to any basis. The trace is a [[similarity invariance|similarity invariant]] under linear transformations. In the lattice basis, the rotation operation must map every lattice point into an integer number of lattice vectors, so the entries of the rotation matrix in the lattice basis – and hence the trace – are necessarily integers. Similar as in other proofs, this implies that the only allowed rotational symmetries correspond to 1,2,3,4 or 6-fold invariance.  For example, wallpapers and crystals cannot be rotated by 45° and remain invariant, the only possible angles are: 360°, 180°, 120°, 90° or 60°.

'''Example'''
*Consider a 60° (360°/6) rotation matrix with respect to the [[oblique lattice]] basis for a [[tessellation|tiling]] by equilateral triangles.
::<math>\begin{bmatrix} 0 & -1 \\ 1 & 1 \end{bmatrix}</math>
:The trace is still 1. The [[determinant]] (always +1 for a rotation) is also preserved.

The general crystallographic restriction on rotations does ''not'' guarantee that a rotation will be compatible with a specific lattice. For example, a 60° rotation will not work with a square lattice; nor will a 90° rotation work with a rectangular lattice.

==Higher dimensions==
When the dimension of the lattice rises to four or more, rotations need no longer be planar; the 2D proof is inadequate. However, restrictions still apply, though more symmetries are permissible. For example, the [[hypercubic lattice]] has an eightfold rotational symmetry, corresponding to an eightfold rotational symmetry of the [[hypercube]]. This is of interest, not just for mathematics, but for the physics of quasicrystals under the [[Aperiodic tiling#Cut-and-project method|cut-and-project theory]]. In this view, a 3D quasicrystal with 8-fold rotation symmetry might be described as the projection of a slab cut from a 4D lattice.

The following 4D rotation matrix is the aforementioned eightfold symmetry of the [[hypercube]] (and the [[cross-polytope]]):
:<math>A = \begin{bmatrix} 0 & 0 & 0 & -1 \\ 1 & 0 & 0 & 0 \\ 0 & -1 & 0 & 0 \\ 0 & 0 & -1 & 0 \end{bmatrix}.</math>
Transforming this matrix to the new coordinates given by 
:<math>B = \begin{bmatrix} -1/2 & 0 & -1/2 & \sqrt 2/2 \\ 1/2 & \sqrt 2/2 & -1/2 & 0 \\ -1/2 & 0 & -1/2 & -\sqrt 2/2 \\ -1/2 & \sqrt 2/2 & 1/2 & 0 \end{bmatrix}</math> will produce:
:<math>B A B^{-1} = \begin{bmatrix} \sqrt 2/2 & \sqrt 2/2 & 0 & 0 \\ -\sqrt 2/2 & \sqrt 2/2 & 0 & 0 \\ 0 & 0 & -\sqrt 2/2 & \sqrt 2/2 \\ 0 & 0 & -\sqrt 2/2 & -\sqrt 2/2 \end{bmatrix}.</math>
This third matrix then corresponds to a rotation both by 45° (in the first two dimensions) and by 135° (in the last two). Projecting a slab of hypercubes along the first two dimensions of the new coordinates produces an [[Ammann–Beenker tiling]] (another such tiling is produced by projecting along the last two dimensions), which therefore also has 8-fold rotational symmetry on average.

The [[A4 lattice]] and [[F4 lattice]] have order 10 and order 12 rotational symmetries, respectively.

To state the restriction for all dimensions, it is convenient to shift attention away from rotations alone and concentrate on the integer matrices {{Harv|Bamberg|Cairns|Kilminster|2003}}. We say that a matrix A has '''[[order (group theory)|order]]''' ''k'' when its ''k''-th power (but no lower), A<sup>''k''</sup>, equals the identity. Thus a 6-fold rotation matrix in the equilateral triangle basis is an integer matrix with order 6. Let Ord<sub>''N''</sub> denote the set of integers that can be the order of an ''N''×''N'' integer matrix. For example, Ord<sub>2</sub> = {1, 2, 3, 4, 6}. We wish to state an explicit formula for Ord<sub>''N''</sub>.

Define a function ψ based on [[Euler's totient function]] φ; it will map positive integers to non-negative integers. For an odd [[prime number|prime]], ''p'', and a positive integer, ''k'', set ψ(''p''<sup>''k''</sup>) equal to the totient function value, 
φ(''p''<sup>''k''</sup>), which in this case is ''p''<sup>''k''</sup>−''p''<sup>''k−1''</sup>. Do the same for ψ(2<sup>''k''</sup>) when ''k'' > 1. Set ψ(2) and ψ(1) to 0. Using the [[fundamental theorem of arithmetic]], we can write any other positive integer uniquely as a product of prime powers, ''m'' = Π<sub>α</sub> ''p''<sub>α</sub><span><sup>''k''<sub> α</sub></sup></span>; set ψ(''m'') = Σ<sub>α</sub> ψ(''p''<sub>α</sub><sup><span>''k''<sub> α</sub></span></sup>). This differs from the totient itself, because it is a sum instead of a product.

The crystallographic restriction in general form states that Ord<sub>''N''</sub> consists of those positive integers ''m'' such that ψ(''m'') ≤ ''N''.

:{| class="wikitable"
|+ '''Smallest dimension for a given order'''{{OEIS2C|id=A080737}}
|- align="center"
| ''m'' || 1 || 2 || 3 || 4 || 5 || 6 || 7 || 8 || 9 || 10 || 11 || 12 || 13 || 14 || 15 || 16 || 17 || 18 || 19 || 20 || 21 || 22 || 23 || 24 || 25 || 26 || 27 || 28 || 29 || 30 || 31
|- align="center"
| ψ(''m'') || 0 || 0 || 2 || 2 || 4 || 2 || 6 || 4 || 6 ||  4 || 10 ||  4 || 12 ||  6 || 6 || 8 || 16 || 6 || 18 || 6 || 8 || 10 || 22 || 6 || 20 || 12 || 18 || 8 || 28 || 6 || 30
|}

For ''m''>2, the values of ψ(''m'') are equal to twice the [[Algebraic number#Properties|algebraic degree]] of cos(2&pi;/''m''); therefore, ψ(''m'') is strictly less than ''m'' and reaches this maximum value if and only if ''m'' is a [[prime number|prime]].

These additional symmetries do not allow a planar slice to have, say, 8-fold rotation symmetry. In the plane, the 2D restrictions still apply. Thus the cuts used to model quasicrystals necessarily have thickness.

Integer matrices are not limited to rotations; for example, a reflection is also a symmetry of order 2. But by insisting on determinant +1, we can restrict the matrices to [[proper rotation]]s.

==Formulation in terms of isometries==

The crystallographic restriction theorem can be formulated in terms of [[isometry|isometries]] of [[Euclidean space]]. A set of isometries can form a [[group (mathematics)|group]]. By a ''discrete isometry group'' we will mean an isometry group that maps each point to a discrete subset of '''R'''<sup>''N''</sup>, i.e. the [[Group action (mathematics)#Orbits and stabilizers|orbit]] of any point is a set of [[isolated point]]s. With this terminology, the crystallographic restriction theorem in two and three dimensions can be formulated as follows.

:For every discrete [[isometry group]] in two- and three-dimensional space which includes translations spanning the whole space, all isometries of finite [[order (group theory)|order]] are of order 1, 2, 3, 4 or 6.
Isometries of order ''n'' include, but are not restricted to, ''n''-fold rotations. The theorem also excludes ''S<sub>8</sub>'', ''S<sub>12</sub>'', ''D<sub>4d</sub>'', and ''D<sub>6d</sub>'' (see [[point groups in three dimensions]]), even though they have 4- and 6-fold rotational symmetry only.
Rotational symmetry of any order about an axis is compatible with translational symmetry along that axis.

The result in the table above implies that for every discrete isometry group in four- and five-dimensional space which includes translations spanning the whole space, all isometries of finite order are of order 1, 2, 3, 4, 5, 6, 8, 10, or 12.

All isometries of finite order in six- and seven-dimensional space are of order 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14, 15, 18, 20, 24 or 30 .

==See also==
*[[Crystallographic point group]]
*[[Crystallography]]

==Notes==
<references/>

== References ==
* {{Citation
  | last1 =Bamberg | first1 =John
  | last2 =Cairns | first2 =Grant
  | last3 =Kilminster | first3 =Devin
  | title =The crystallographic restriction, permutations, and Goldbach's conjecture
  | journal =[[American Mathematical Monthly]]
  |date=March 2003
  | volume =110  | issue=3 | pages =202–209
  | url =http://cage.ugent.be/~bamberg/Research_files/The%20crystallographic%20restriction,%20permutations%20and%20Goldbach%27s%20conjecture.pdf
  | doi =10.2307/3647934
  | jstor =3647934
| citeseerx =10.1.1.124.8582
  }}
* {{Citation
  | last =Elliott  | first =Stephen
  | author-link =Stephen Elliott (mathematician)
  | title =The Physics and Chemistry of Solids
  | publisher =Wiley
  | year =1998
  | isbn =978-0-471-98194-7
}}
* {{Citation
  | last =Coxeter  | first =H. S. M.
  | author-link =H. S. M. Coxeter
  | title =Introduction to Geometry
  | publisher =Wiley
  | year =1989
  | edition =2nd
  | isbn =978-0-471-50458-0
}}
* {{Citation
  | last =Scherrer | first =W.
  | title =Die Einlagerung eines regulären Vielecks in ein Gitter
  | journal =[[Elemente der Mathematik]]
  | year =1946
  | volume =1  | issue =6 | pages =97–98
  | url =http://www-gdz.sub.uni-goettingen.de/cgi-bin/digbib.cgi?PPN378850199_0001
}}
* {{Citation
  | last =Shechtman | first =D.
  | last2 =Blech | first2 =I.
  | last3 =Gratias | first3 =D.
  | last4 =Cahn | first4 =JW
  | year =1984
  | title =Metallic phase with long-range orientational order and no translational symmetry
  | journal =[[Physical Review Letters]]
  | volume =53  | issue =20  | pages =1951–1953
  | doi = 10.1103/PhysRevLett.53.1951 | bibcode=1984PhRvL..53.1951S
| doi-access =free}}

==External links==
*[http://www-history.mcs.st-and.ac.uk/~john/geometry/Lectures/A2.html The crystallographic restriction]
*[https://www.xtal.iqf.csic.es/Cristalografia/parte_03_1_1-en.html The crystallographic restriction theorem by CSIC]

[[Category:Crystallography]]
[[Category:Group theory]]
[[Category:Theorems in algebra]]
[[Category:Articles containing proofs]]