Category: GoV

Penrose’s aperiodic tilings

Published March 9, 2021 by lievenlb

Around 1975 Sir Roger Penrose discovered his aperiodic P2 tilings of the plane, using only two puzzle pieces: Kites (K) and Darts (D)

The inner angles of these pieces are all multiples of $36^o = \tfrac{180^o}{5}$, the short edges have length $1$, and the long edges have length $\tau = \tfrac{1+\sqrt{5}}{2}$, the golden ratio. These pieces must be joined together so that the colours match up (if we do not use this rule it is easy to get periodic tilings with these two pieces).

There is plenty of excellent online material available:

The two original Martin Gardner Scientific American articles on Penrose tiles have been made available by the MAA, and were reprinted in “Penrose Tiles to Trapdoor Ciphers”. They contain most of Conway’s early discoveries about these tilings, but without proofs.
A JavaScript application by Kevin Bertman to play around with these tilings. You can deflate and inflate tilings, find forced tiles and much more. Beneath the app-window there’s a detailed explanation of all the basics, including inflation and deflation of the P2-tiles, the seven types of local vertex configurations (naming by Conway, of course),

proofs of aperiodicity (similar to the one for Conway’s musical sequences), that every tile lies within an ace (similar to the LSL-subword in musical sequences) with application to local isomorphism (again similar to the $1$-dimensional case).
Course notes of an Oxford masterclass in geometry Lectures on Penrose tilings by Alexander Ritter, again with proofs of all of the above and a discussion about the Cartwheel tilings (similar to that in the post on musical sequences), giving an algorithm to decide whether or not a partial tiling can be extended to the entire plane, or not.

There’s no point copying this material here. Rather, I’d like to use some time in this GoV series of posts to talk about de Bruijn’s pentagrid results. For this reason, I now need to make the connection with Penrose’s ‘other’ tilings, the P3 tiles of ‘thin’ and ‘thick’ rhombi (sometimes called ‘skinny’ and ‘fat’ rhombi).

Every Penrose P2-tiling can be turned into a P3-rhombic tiling, and conversely.

From kites and darts to rhombi: divide every kite in two halves along its line of reflection. Then combine darts and half-kites into rhombi where a fat rhombus consists of a dart and two half-kites, joined at a long edge, and a skinny rhombus is made of two half-kites, joined at a short edge. This can always be done preserving the gluing conditions. It suffices to verify this for the deflated kite and dart (on the left below) and we see that the matching colour-conditions are those of the rhombi.

From rhombi to kites and darts: divide every fat rhombus into a dart (placed at the red acute angle) and two half-kites, joined at a long edge, and divide every skinny rhombus into two half-kites along its short diagonal.

All results holding for Kites and Darts tilings have therefore their versions for Rhombic tilings. For example, we have rhombic deflation and inflation rules

A Rhombic tiling can be seen as an intermediate step in the inflation process of a Penrose tiling $P$. Start by dividing all kites in two halves along their long diagonal and all darts in two halves along their short diagonal (the purple lines below).

If we consider all original black lines together with the new purple ones, we get a tiling of the plane by triangles and we call this the $A$-tiling. There are two triangles, s small triangle $S_A$ (the whiter ones, with two small edges ofd length $1$ and one edge of length $\tau$) and a large triangle $L_A$ (the greyer ones, with two edges of length $\tau$ and one edge of length $1$).

Next, we remove the black lines joining an $S_A$-triangle with an $L_A$-triangle, and get another tiling with triangles, the $B$-tiling, with two pieces, a large triangle $L_B$ with two sides of length $\tau$ and one side of length $1+\tau$ (the white ones) (the union of an $L_A$ and a $S_A$), and a small triangle $S_B$ which has the same form as $L_A$ (the greyer ones). If we now join two white triangles $L_B$ along their common longer edge, and two greyer triangles $S_B$ along their common smaller edge, we obtain the Rhombic tiling $R$ corresponding to $P$.

We can repeat this process starting with the Rhombic tiling $R$. We divide all fat rhombi in two along their long diagonals, and the skinny rhombi in two along their short diagonals (the purple lines). We obtain a tiling of the plane by triangles, which is of course just the $B$-tiling above.

Remove the edge joining a small $S_B$ with a large $L_B$ triangle, then we get a new tiling by triangles, the $\tau A$-tiling consisting of large triangles $L_{\tau A}$ (the white ones) with two long edges of length $1+\tau$ and one short edge of length $\tau$ (note that $L_{\tau A}$ is $\tau$ times the triangle $L_A$), and a smaller one $S_{\tau A}$ (the grey ones) having two short edges of length $\tau$ and one long edge of length $1+\tau$ (again $S_{\tau A}$ is $\tau$ times the triangle $S_A$). If we join two $L_{\tau A}$-triangles sharing a common long edge we obtain a Kite ($\tau$-times larger than the original Kite) and if we joint two $S_{\tau A}$-triangles along their common small edge we get a Dart ($\tau$ times larger than the original Dart). The Penrose tiling we obtain is the inflation $inf(P)$ of the original $P$.

If we repeat the whole procedure starting from $inf(P)$ instead of $P$ we get, in turn, triangle tilings of the plane, subsequently the $\tau B$-tiling, the $\tau^2 A$-tiling, the $\tau^2 B$-tiling and so on, the triangles in a $\tau^n A$-tiling being of size $\tau^2$-times those of the $A$-tiling, and those in the $\tau^n B$-tiling $\tau^n$-times as large as those of the $B$-tiling.

The upshot of this is that we can associate to a Penrose tiling $P$ of the plane a sequence of $0$’s and $1$’s. Starting from $P$ we have an infinite sequence of associated triangle tilings
\[
A,~B,~\tau A,~\tau B,~\tau^2 A, \dots,~\tau^n A,~\tau^n B,~\tau^{n+1} A, \dots \]
with larger and larger triangle tiles. Let $p$ be a point of the plane lying in the interior of an $A$-tile, then we define its index sequence
\[
i(P,p) = (x_0,x_1,x_2,\dots ) \quad x_{2n} = \begin{cases} 0~\text{if $p \in L_{\tau^n A}$} \\ 1~\text{if $p \in S_{\tau^n A}$} \end{cases}~ x_{2n+1} = \begin{cases} 0~\text{if $p \in L_{\tau^n B}$} \\ 1~\text{if $p \in S_{\tau^n B}$} \end{cases} \]
That is, $x_0=1$ if $p$ lies in a small triangle of the $A$-tiling and $x_0=1$ if it lies in a large triangle, $x_1=1$ if $p$ lies in a small triangle of the $B$-tiling and is $0$ if it lies in a large $B$-triangle, and so on.

The beauty of this is that every infinite sequence $(x_0,x_1,x_2,\dots )$ whose terms are $0$ or $1$, and in which no two consecutive terms are equal to $1$, is the index sequence $i(P,p)$ of some Penrose tiling $P$ in some point $p$ of the plane.

From the construction of the sequence of triangle-tilings it follows that a small triangle is part of a large triangle in the next tiling. For this reason an index sequence can never have two consecutive ones. An index sequence gives explicit instructions as to how the Penrose tiling is constructed. In each step add a large or small triangle as to fit the sequence, together with the matching triangle (the other half of a Penrose or Rhombic tile). Next, look at the patches $P_i$ of Kites and Darts (in the $\tau^i A$ tiling), for $i=0,1,\dots$, then $P_{i-k}$ is contained in the $k$-times inflated $P_i$, $i^k(P_i)$ (without rescaling), for each $i$ and all $1 \leq k \leq i$. But then, we have a concentric series of patches
\[
P_0 \subset i(P_1) \subset i^2(P_2) \subset \dots \]
filling the entire plane.

The index sequence depends on the choice of point $p$ in the plane. If $p’$ is another point, then as the triangle tiles increase is size at each step, $p$ and $p’$ will lie in the same triangle-tile at some stage, and after that moment their index sequences will be the same.

As a result, there exists an uncountable infinity of distinct Penrose tilings of the plane.

From the discussions we see that a Penrose tiling is determined by the index-sequence in any point in the plane $(x_0,x_1,\dots )$ consisting of $0$’s and $1$’s having no two consecutive $1$’s, and that another such sequence $(x’_0,x’_1,\dots )$ determines the same Penrose tilings if $x_n=x’_n$ for all $n \geq N$ for some number $N$. That is, the Penrose tiling is determined by the the equivalence class of such sequences, and it is easy to see that there are uncountably many such equivalence classes.

If you want to play a bit with Penrose tiles, you can order the P2 tiles or the P3 tiles from
Cherry Arbor Design.

Conway’s musical sequences (2)

Published March 3, 2021 by lievenlb

A Conway musical sequence is an infinite word in $L$ and $S$, containing no two consecutive $S$’s nor three consecutive $L$’s, such that all its inflations remain musical sequences.

We’ve seen that such musical sequences encode an aperiodic tiling of the line in short ($S$) and long ($L$) intervals, and that such tilings are all finite locally isomorphic.

But, apart from the middle $C$-sequences (the one-dimensional cartwheel tilings) we gave no examples of such tilings (or musical sequences). Let’s remedy this!

Take any real number $c$ as long as it is not an integral combination of $1$ and $\tfrac{1}{\tau}$ (with $\tau$ the golden ratio) and assign to any integer $a \in \mathbb{Z}$ a tile
\[
P_c(a) = \begin{cases} S \\ L \end{cases} ~\text{iff}~\lceil c+(a+1)\frac{1}{\tau} \rceil – \lceil c+a \frac{1}{\tau} \rceil = \begin{cases} 0 \\ 1 \end{cases} \]
(instead of ceilings we might have taken floors, because of the restriction on $c$).

With a little bit of work we see that the deflated word determined by $P_c$ is again of this type, more precisely $def(P_c) = P_{-(c-\lfloor c \rfloor)\frac{1}{\tau}}$. But then it also follows that inflated words are of this type, meaning that all $P_c$ define a musical sequence.

Let’s just check that these sequences satisfy the gluing restrictions. If there is no integer between $c+a\tfrac{1}{\tau}$ and $c+(a+1)\tfrac{1}{\tau}$, because $2 \tfrac{1}{\tau} \approx 1.236$ there must be an interval in the preceding and the following $\tfrac{1}{\tau}$-interval, showing that an $S$ in the sequence has an $L$ on its left and right, so there are no two consecutive $S$’s in the sequences.

Similarly, if two consecutive $\tfrac{1}{\tau}$-intervals have an integer in them, the next interval cannot contain an integer as $3 \tfrac{1}{\tau} \approx 1.854 < 2$.

Now we come to the essential point: these sequences can be obtained by the cut-and-project method.

Take the line $L$ through the origin with slope $\tfrac{1}{\tau}$ and $L^{\perp}$ the line perpendicular it.

Consider the unit square $H$ and $H_{\vec{\gamma}}=H + \vec{\gamma}$ its translation under a shift vector $\vec{\gamma}=(\gamma_x,\gamma_y)$ and let $\pi$ (or $\pi^{\perp}$) be the orthogonal projection of the plane onto $L$ (or onto $L^{\perp}$). One quickly computes that
\[
\pi(a,b) = (\frac{\tau^2 a + \tau b}{1+\tau^2},\frac{\tau a + b}{1+\tau^2}) \quad \text{and} \quad
\pi^{\perp}(a,b) = (\frac{a-\tau b}{1+\tau^2},\frac{\tau^2b-\tau a}{1+\tau^2}) \]
In the picture, we take $\vec{\gamma}=(c,-\tau c)$.

The window $W$ will be the strip, parallel with $L$ with basis $\pi^{\perp}(H_{\vec{\gamma}})$.

We cut the standard lattice $\mathbb{Z}^2$, of all points with integer coordinates in the plane, by retricting to the window $\mathcal{P}=\mathbb{Z}^2 \cap W$.

Next, we project $\mathcal{P}$ onto the line $L$, and we get a set of endpoints of intervals which divide the line $L$ into short intervals of length $\tfrac{1}{\sqrt{1+\tau^2}}$ and long intervals of length $\tfrac{\tau}{\sqrt{1+\tau^2}}$.

For $(a,b) \in W$, the interval will be short if $(a,b+1) \in W$ and long if $(a+1,b) \in W$.

Because these intervals differ by a factor $\tau$ in length, we get a tiling of the line by short intervals $S$ and long intervals $L$. It is easy to see that they satisfy the gluing restrictions (remember, no two consecutive short intervals and no three consecutive long intervals): the horizontal width of the window $W$ is $1+\tau \approx 2.618$ (so there cannot be three consecutive long intervals in the projection) and the vertical width of the window $W$ is $1+\tfrac{1}{\tau} = \tau \approx 1.618$ so there cannot be two consecutive short intervals in the projection.

Ready for the punchline?

The sequence obtained from projecting $\mathcal{P}$ is equal to the sequence $P_{(1+\tau^2)c}$. So, we get all musical sequences of this form from the cut-and-project method!

On $L^{\perp}$ the two end-points of the window are
\[
\begin{cases}
\pi^{\perp}(c+1,-\tau c) = (\frac{(1+\tau^2)c+1}{1+\tau^2},- \tau \frac{(1+\tau^2)c+1}{1+\tau^2}) \\
\pi^{\perp}(c,-\tau c+1) = (\frac{(1+\tau^2)c-\tau}{1+\tau^2},-\tau \frac{(1+\tau^2)c-\tau}{1+\tau^2})
\end{cases} \]
Therefore, a point $(a,b) \in \mathbb{Z}^2$ lies in the window $W$ if and only if
\[
(1+\tau^2)c-\tau < a-\tau b < (1+\tau^2)c+1 \] or equivalently, if \[ (1+\tau^2)c+(b-1)\tau < a < (1+\tau^2)c+b \tau + 1 \] Observe that \[ \lceil (1+\tau^2)c + b\tau \rceil - \lceil (1+\tau^2)c+(b-1)\tau \rceil = P_{(1+\tau^2)c}(b-1) + 1 \in \{ 1,2 \} \] We separate the two cases: (1) : If $\lceil (1+\tau^2)c + (b+1)\tau \rceil - \lceil (1+\tau^2)c+b \tau \rceil =1$, then there must be an integer $a$ such that $(1+\tau^2)c +(b+1) \tau -1 < a < (1+\tau^2) b+1$, and this forces $\lceil (1+\tau^2)c + (b+2)\tau \rceil - \lceil (1+\tau^2)c+(b+1)\tau \rceil =2$. With $b_i = (1+\tau^2)c+(b+i)\tau$ and $d_i = b_i+1$ we have the situation

and from the inequalities above this implies that both $(a+1,b+1)$ and $(a+1,b+2)$ are in $W$, giving a short interval $S$ in the projection.

(2) : If $\lceil (1+\tau^2)c + (b+1)\tau \rceil – \lceil (1+\tau^2)c+b \tau \rceil =1$, then there must be an integer $a$ such that $(1+\tau^2)c+b \tau < a < (1+\tau^2)cv + (b+1)\tau -1$, giving the situation

giving from the inequalities that both $(a+1,b+1)$ and $(a+2,b+1)$ are in $W$, giving a long interval $L$ in the projection, finishing the proof.

Conway’s musical sequences

Published March 2, 2021 by lievenlb

Before we’ll come to applications of quasicrystals to viruses it is perhaps useful to illustrate essential topics such as deflation, inflation, aperiodicity, local isomorphism and the cut-and project method in the simplest of cases, that of $1$-dimensional tilings.

We want to tile the line $\mathbb{R}^1$ with two kinds of tiles, short ($S$) and ($L$) long intervals, differing by a golden ratio factor $\tau=\tfrac{1}{2}(1+\sqrt{5}) \approx 1.618$.

Clearly, no two tiles may overlap and we impose a gluing restriction: there can be no two consecutive $S$-intervals in the tiling, and no three consecutive $L$-intervals.

The code of a tiling is a doubly infinite word in $S$ and $L$ such that there are no two consecutive $S$’s nor three consecutive $L$’s. For example
\[
\sigma = \dots LSLLSLS\underline{L}LSLLSL \dots \]
We underline one tile to distinguish the sequence from shifts of it.

Conway’s musical sequences will be special codes (or tilings), allowing for the inverse operations of inflation and deflation, terms coined by John Conway in relation to Penrose tilings. The musical sequences are important to understand Conway’s worms (sometimes called “wormholes”) which are strings of Long and Short bow ties in a Penrose tiling, and to measure the distances between Amman bars. In fact, many of the properties of Penrose tilings and $3$-dimensional quasicrystals (for example, local isomorphism) have their counterparts for tilings having a musical sequence as code.

Conway’s investigations of Penrose tiles held up work on the ATLAS-project and caused some problems at home:

“In pursuing his investigations, he unsurped some of his wife Eileen’s territory, covering the dining table with an infinite nuisance of tiles. He cut them out himself, causing his right hand to hurt with cramps for days. To Eileen’s dismay, he studied the dining table mosaic for a year, relegating family meals to the kitchen and prohibiting dinner parties.”

From “Genius at Play – The curious mind of John Horton Conway” by Siobhan Roberts

Let’s investigate inflation and deflation of these tilings.

The point of the golden factor $\tau$ is to allow for deflation. That is, we can replace a tiling by another one with tiles $S$ and $L$ both a factor $\tfrac{1}{\tau}=\tau-1 \approx 0.618$ smaller than the original tiles. If the original $S$-tile has length $a$ (and the $L$-tile length $\tau a$), then the new tile $S$ wil have length $\tfrac{1}{\tau}a$ and the new $L$-tile length $a$.
We do this by replacing each old $S$-tile by a new $L$-tile, and to break up any old $L$-tile in a new $L$ and new $S$-tile, as $\tau a = a + \tfrac{1}{\tau}a$ (note that $\tau^2=\tau+1$)

To get the code of the deflated tiling we replace each letter $S$ by a letter $L$ and each $L$ by $LS$. The underlined letter will be the first letter of the deflated underlined letter in the original sequence. The deflated sequence of the one above is
\[
def(\sigma) = \dots LSLLSLSLLSL\underline{L}SLSLLSLSLLS \dots \]
and it is easy to see that the deflated tiling satisfies again the gluing condition.

Certain of these tilings (not all!) allow for an inverse to deflation, called inflation, increasing the size of the tiles by a factor $\tau$.

Starting from a tiling we divide each $L$-tile in half and these mid-points will be end-points of the tiles in the new tiling, erasing all endpoints of the original one. The inflated tiling will have two sorts of tiles, a new short one $S$ of length $\tau a$ obtained from the end-half of an original $L$-tile, followed b the start-half of an original $L$-tile, and a new long tile $L$ of length $\tau^2 a = (\tau+1) a$, made of the end-half of an original $L$, followed by an original $S$, followed by the start-half of an original $L$.

We get the code of the inflated tiling by replacing first each $L$ by $ll$ and subsequently replace each word $lSl$ by a letter $L$ and each $ll$ by $S$. An example,
\[
\sigma = \dots LSLLSLSLLSL\underline{L}SLSLLSLSLLS \dots \\
\dots llSllllSllSllllSlll\underline{l}SllSllllSllSllllS \dots \\
inf(\sigma) = \dots (l)LSLLSLS\underline{L}LSLLS(lS) \dots \]

But, the inflated tiling may no longer satisfy the gluing condition. An example
\[
\dots LSLSLSL \dots \mapsto \dots llSllSllSll \dots \mapsto \dots (l)LLL(l) \dots \]

A Conway musical sequence is the code of a tiling $\sigma$ such that all its consecutive inflations $inf^n(\sigma)$ satisfy the gluing condition. For the corresponding Conway tiling $\sigma$ we have that
\[
def(inf(\sigma))=\sigma=inf(def(\sigma)) \]

Let’s construct at least two such Conway tilings (later we’ll see that there are uncountably many). Take $C_n=def^n(LS\underline{L})$ and write it in a special form to highlight symmetries.
\begin{eqnarray*}
C_0 =& (L.S)\underline{L} \\
C_1 =& L(S.L)\underline{L}S \\
C_2 =& LSL(L.S)\underline{L}SL \\
C_3 =& LSLLSL(S.L)\underline{L}SLLS \\
C_4 =& LSLLSLSLLSL(L.S)\underline{L}SLLSLSL
\end{eqnarray*}

The even terms have middle-part $(L.S)$ and the odd ones $(S.L)$. The remaing left and right parts are each others reflexion (or part of it). This is easily seen by induction as are the inclusions
\[
C_0 \subset C_2 \subset C_4 \subset \dots \subset C_{even} \quad \text{and} \quad C_1 \subset C_3 \subset C_5 \subset \dots \subset C_{odd} \]

$C_{even}$ and $C_{odd}$ are special Conway musical sequences, called the middle $C$-sequences, and are each others inflation and deflation. If you are familiar with Penrose tilings, these are the $1$-dimensional counterparts of the cartwheel Penrose tiling (here with the $10$ Conway worms emanating from the center, and with the borders of the first few cartwheels drawn).

A direct consequence of inflation on Conway’s musical sequences is that the corresponding tiling is aperiodic, that is, it has no translation symmetry.

For, inflation only depends on the local configuration of tiles, so if translation by $R$ is a symmetry of a musical sequence $\sigma$ then it is also a symmetry of $inf(\sigma)$, and so also of $inf^n(\sigma)$. But for large $n$ we will have that $R < \tau^n a$ (with $a$ the size of the tiles in $\sigma$). But then a tile in $inf^n(\sigma)$ and its translation by $R$ must overlap which is impossible if $+R$ is a translation symmetry of $inf^n(\sigma)$. Done!

Returning to the middle C-sequences, what was the point of starting with $C_0 = LSL$? Well, it follows directly from the gluing restrictions that any letter in a musical sequence is part of a subword $LSL$ of $\sigma$. But then, every finite subword $W$ of $\sigma$ is also a subword of $C_{2n}$ for some large $n$.

For, let $d$ be the length of the interval corresponding to $W$ and choose $n$ such that $d > \tau^{2n} a$ then the interval of the line corresponding to $W$ is contained in a single tile in $inf^{2n}(\sigma)$ and this tile belongs to a subword $LSL$ of $inf^{2n}(\sigma)$. But then $W$ will be a subword of the $2n$-th deflation of that interval $LSL \subset inf^{2n}(\sigma)$, which is $C_{2n}$.

Or, as Conway would phrase it with respect to Penrose tilings (quote again from Siobhan Roberts’ book)

“Every points is in the cartwheel somewhere. If you jab your finger anywhere, on any point anywhere on teh pattern, you are part of a cartwheel. The whole ting is overlapping cartwheels.”

An immediate consequence is the local isomorphism theorem: Every subword of a musical sequence $\sigma$ appears infinitely many times as subword of any other musical sequence. That is, one cannot distinguish two tilings of the line with musical sequence codes from each other by looking at finite intervals!

The argument is similar to the one above. The finite interval corresponding to the subword lies in a unique tile of $inf^n(\sigma)$ for $n$ large enough. Now, take another musical sequence $\mu$ and consider any of the infinitely many tiles of the same type in $inf^n(\mu)$, then $def^n$ of such a tile will contain the subword in $\phi$.

Another time, we’ll see that musical sequences can be produced by the ‘cut-and-project’-method (what I called the ‘windows’-method before).
This time we will project parts of the standard $2$-dimensional lattice $\mathbb{Z}^2$ onto the line, which is a lot easier to visualise than de Bruijn’s projection from $\mathbb{R}^5$ to produce Penrose tilings or the projection from six dimensional space to harvest quasicrystals.