ON THE SPECTRUM OF THE PENROSE LAPLACIAN
MICHAEL DAIRYKO, CHRISTINE HOFFMAN, JULIE PATTYSON, AND HAILEE PECK
Abstract. Since the early 1960’s, aperiodic tilings have been a topic of par- ticular interest to mathematicians. The Penrose tiling is one example of such atiling.Inthispaper,wepresentbasicsubdivisionrulesforgeneratingthe Penrose tiling using a set of four Robinson tiles. We introduce code written in MatLab that creates and stores the Penrose tiling as a set of data struc- tures and produces the Laplacian matrix of the tiling. In doing so, we provide insight into the spectrum of the Penrose tiling.
Contents 1. Introduction 1 2. Generating a Penrose Tiling 2 2.1. Substitution Method 3 2.2. Rules for Substitution with Robinson Triangles 4 2.3. Computational Methods for Penrose Tiling 7 3. The Laplacian and Spectrum 7 3.1. Computational Methods of the Laplacian and Spectrum 9 4. The Hausdor↵Dimension 10 4.1. Computational Methods for The Hausdor↵Dimension 11 5. Results 11 6. Conclusion and Future Work 13 7. Acknowledgments 14 Appendix A: MatLab Code 15 A.1. genPenTiling.m 15 A.2. tilingLaplacian.m 22 A.3. plotMat.m 25 A.4. cdf.m 25 A.5. ndsTable.m 26 References 28
1. Introduction Before 1982, crystalline structures were defined as having periodic lattice struc- ture. However, on April 8, 1982, Dan Shechtman discovered an aluminum-manganese quasicrystalline structure that had long-range order but did not exhibit the peri- odic patterns that characterized crystals [12]. Instead, the structure had icosahedral
Date:August1,2013. This work was supported by NSF grant DMS-0739338. 1 2 M.DAIRYKO,C.HOFFMAN,J.PATTYSON,ANDH.PECK symmetry that was not possible, by definition, in crystal form. In 2011, Shecht- man was awarded the Nobel Prize in Chemistry for his discovery of and work with quasicrystals. Since his groundbreaking discovery, many examples of quasicrystals have been examined. However, mathematical interest in periodicity began before these discoveries of the physical applications occurred. The following discussion of this history comes from [8]. Questions relating to tiling and periodicity began with Hao Wang’s work in 1960 on tiling the plane with Wang dominoes, which are sets of unit squares with colored edges. This work began because of its relation to questions in symbolic logic. Wang conjectured that there was a way to tile the plane periodically given any set of tiles so that bordering edges share the same color. In 1964, Robert Berger provided a counterexample to the Wang conjecture. Berger showed that there is a set (of more than 20,000 Wang dominoes) that tiles only nonperiodically, and later much smaller sets (of 104 and 92 Wang dominoes) with the same property were found.
Figure 1. Penrose rhombus and diamond prototiles divided into Robinson triangles
Berger’s findings inspired further exploration with polygonal tiles. Raphael Robinson and Robert Ammann, in 1971 and 1977, respectively, constructed two di↵erent sets of 6 “square” tiles that tile only nonperiodically. These tiles are like jigsaw puzzle pieces which can be rotated and reflected. The projections and indentations on these tiles force specific tiles to match together. In 1973, Roger Penrose found a set of six “non-square” tiles that tile only nonpe- riodically. He later reduced it to a set of four tiles, and then two tiles. Each of these sets of Penrose tiles generate nonperiodic tilings, which are two-dimensional models of quasicrystals. Robinson took a set of two prototiles, a rhombus and diamond, and divided them into a set of four triangles as shown in Figures 1 and 2. These tiles are known as Robinson triangles. Our work is based on these Penrose tilings, which can be generated from any of the sets of prototiles that Penrose found. In Section 2, we discuss the method we use to construct the Penrose tiling and give an explanation of MatLab code that constructs finite tilings. In Section 3, we define the Laplacian matrix and spectrum and explain the corresponding MatLab code. In Section 4, we examine the Hausdor↵dimension and describe the MatLab code developed to estimate the Hausdor↵dimension of the finite Penrose tilings.
2. Generating a Penrose Tiling Let us first present relevant definitions about tilings. ON THE SPECTRUM OF THE PENROSE LAPLACIAN 3
Definition 1. [7] Prototiles are a set of finite inequivalent tiles (i.e. are not equivalent under rigid motions, expansions, or contractions). Definition 2. [7] A tiling is an arrangement of tiles, such that their union covers and packs R2 so that distinct tiles have non-intersecting interiors. Definition 3. [9] A nonperiodic tiling does not have a period. Definition 4. [9] An aperiodic set of prototiles admits infinitely many tilings of the plane, none of which are periodic. Remark 1. An iteration is one application of the substitution method to a set of tiles. There are three main methods that can be used to construct the Penrose tiling. The matching method assigns matching rules to the edges of the prototiles, only allowing pairings of certain edges. These rules guarantee a nonperiodic tiling. The substitution method starts with a finite subset of the tiling and inflates and sub- divides existing tiles to expand the tiling to cover the plane. The cut and project method takes a higher dimensional lattice and then reduces the order, projecting it to a lower dimension. For a more detailed explanation, see [5] and [1]. 2.1. Substitution Method. In this paper, we construct the Penrose tiling using the substitution method. The substitution we use is based on the set of two Penrose prototiles, the rhombus and diamond (see Figure 1). The two tiles are each sliced temporarily in half into two triangles (called Robinson triangles), resulting in a substitution with four base prototiles: B, P , Y , and W (see Figure 2). As the substitution method for the rhombus and diamond allows for overlap, we implement the Robinson triangle substitution to avoid this problem. This method inflates 1+p5 each triangle prototile by 2 and replaces the inflated tile with a predetermined subdivision.
Figure 2. How to subdivide Robinson triangle prototiles (not drawn to scale)
The rhombus is composed of two triangles, which we label B and P . As shown in Figure 1, Triangle B is inflated and subdivided into three triangles: P 1,Y1, and B1. Triangle P is also inflated and subdivided into three triangles: B2,W2, and 4 M.DAIRYKO,C.HOFFMAN,J.PATTYSON,ANDH.PECK
P 2. Similarly, the diamond is composed of two triangles, which we label Y and W . The triangle Y is inflated and subdivided into two triangles: P 3 and Y 3. The triangle W is also inflated and subdivided into two triangles: W 4 and B4. Figure 2 shows this inflation and subdivision of the four prototiles.
Figure 3. First three iterations with Robinson triangles
We begin with ten prototiles (five triangle B tiles and five triangle P tiles) and use this substitution to generate larger iterations of the Penrose tiling. See Figure 3 for the first two iterations of the tiling using the substitution. For a more detailed explanation see [7] and [11]. Figures 1, 2, and 3 were inspired by and modified from [7].
2.2. Rules for Substitution with Robinson Triangles. The sides of each tri- angle are marked with either one arrow, two arrows, or no arrows. In the tiling, a one arrow side can only be adjacent to another one arrow side, a two arrow side can only be adjacent to a two arrow side, and a side with no arrows can only be adjacent to another side with no arrows (we will refer to this as the third side). In addition, the arrows on the shared sides of adjacent tiles must point in the same direction. The nine rules in Table 1 are the only possible ways that tiles can share edges. For example, the second side of a W tile can only be matched to the second side of a P tile or to the second side of a Y tile (following rules 7 and 8, respectively). Initially, one might think that the second side of a W tile could be matched to the second side of a B tile. However, if you were to match these two sides, you would see that arrows of the second side of each of these tiles do not point in the same direction. It is not possible to rotate these tiles in any way and match them on their second side so that the arrows on the second side of these tiles point in the same direction on their shared edge. Therefore these two tiles will not meet on their respective second sides. In Table 1, the “Shared Edge” column contains a tree diagram for each of the nine rules mentioned above. The two vertices connected by the solid black line represent which two parent prototiles share an edge. The side that they share is the number in superscript. The dashed arrows denote the children of each adjacent ON THE SPECTRUM OF THE PENROSE LAPLACIAN 5 parent. Each child is directly connected only to its parent tile, where the subscripts index the children. The “External Adjacencies” column contains the external adjacencies formed be- tween the new child prototiles after subdivision of the two parent prototiles which shared an edge. We say an adjacency is external if it connects children of dis- tinct prototiles. The adjacencies are separated by commas. Again, the superscript represents the side on which the tiles are adjacent and the subscript denotes the corresponding child. For example, Rule 5 describes P and Y prototiles that are adjacent on the first side (superscript 1). The children of the P parent prototile are B2, W2, and P2. The children of the Y parent prototile are P3 and Y3. For the external adjacencies, the child Y3 is adjacent to the child W2 on the third side (superscript 3). Also, the child P3 is adjacent to the child B2 on the second side (superscript 2). Because we split the original two prototiles (the rhombus and the diamond) into two triangles, each of whose third side cuts through the interior of the rhombus or diamond, the third side of each of these triangles will always only be matched with the third side of the other triangle that made up the original prototile. In this way, the structure of the original two prototile Penrose tiling is preserved. We need to return to the original set of two prototiles instead of four Robinson triangles to obtain an end result of a Penrose tiling. To reconstruct the rhombus, if triangle B is adjacent to triangle P on the third side with no arrows (see Figure 2), we delete this side that they share. Similarly, to reconstruct the diamond, if triangle Y and triangle W are adjacent on their respective third sides, we delete this side that they share (see Figure 2). 6 M.DAIRYKO,C.HOFFMAN,J.PATTYSON,ANDH.PECK
Table 1. Rules for Shared Edges of Adjacent Tiles
Rule Shared Edge External Adjacencies adj P 3 B3 1 P 3 B3, B1 P 1 1 2 1 2 B2 W2 P2 P1 Y1 B1
adj B1 P 1 2 P 2 B2, Y 3 W 3 1 2 1 2 P1 Y1 B1 B2 W2 P2
adj Y 3 W 3 3 Y 2 W 2 3 4 P3 Y3 W4 B4
adj P 2 B2 4 B3 P 3 1 2 B2 W2 P2 P1 Y1 B1
adj P 1 Y 1 5 Y 3 W 3, P 2 B2 3 2 3 2 B2 W2 P2 P3 Y3
adj B1 W 1 6 Y 3 W 3, P 2 B2 1 4 1 4 P1 Y1 B1 W4 B4
adj P 2 W 2 7 P 3 B3 2 4 B2 W2 P2 W4 B4
adj Y 2 W 2 8 P 3 B3 3 4 Y3 P3 W4 B4
adj B2 Y 2 9 B3 P 3 1 3 P1 Y1 B1 P3 Y3 ON THE SPECTRUM OF THE PENROSE LAPLACIAN 7
2.3. Computational Methods for Penrose Tiling.
2.3.1. Explanation of the MatLab Code. This section discusses the MatLab code written to generate finite iterations of the Penrose tiling. See Appendix A.1 for the code and more detailed comments. The goal of the function genPenTiling.m is to construct iterations of the Penrose tiling. The function takes input of oldTiles and outputs newTiles, which are the tiles in the next iteration. We store tiles as structures which each have four fields. The first is a positive integer, denoted as ‘id’ that is unique to each tile in a given iteration. The next is a representation of the color of the tile: ‘p’, ‘b’, ‘y’ and ‘w’, denoting pink, blue, yellow, and white. The next field is an array that contains the id numbers of the tile’s children. Our last field is a list of the tile’s neighbors on each side: first, second, and third, which we denote as ‘n1’, ‘n2’, and ‘n3’ respectively. For example, we see in Figure 3 that the first iteration of the tiling contains five P and five B tiles. So the input for oldTiles is: tile1 = struct ( ‘color’, ‘p’, ‘id’, 1, ‘n3’, 2, ‘n1’ , 10); tile2 = struct ( ‘color’, ‘b’, ‘id’, 2, ‘n3’, 1, ‘n1’ ,3 ); tile3 = struct ( ‘color’, ‘p’, ‘id’, 3, ‘n3’, 4, ‘n1’ ,2); tile4 = struct ( ‘color’, ‘b’, ‘id’, 4, ‘n3’, 3, ‘n1’ ,5); tile5 = struct ( ‘color’, ‘p’, ‘id’, 5, ‘n3’, 6, ‘n1’ ,4); tile6 = struct ( ‘color’, ‘b’, ‘id’, 6, ‘n3’, 5, ‘n1’ ,7); tile7 = struct ( ‘color’, ‘p’, ‘id’, 7, ‘n3’, 8, ‘n1’ ,6); tile8 = struct ( ‘color’, ‘b’, ‘id’, 8, ‘n3’, 7, ‘n1’ ,9); tile9 = struct ( ‘color’, ‘p’, ‘id’, 9, ‘n3’, 10, ‘n1’ ,8); tile10 = struct ( ‘color’, ‘b’, ‘id’, 10, ‘n3’, 9 , ‘n1’ , 1); This sets each P and B tile as neighbors on the third side, as shown in Figure 2. The function genPenTiling.m first creates the new tiles of the next iteration, which are the children of each of the oldTiles. Then the function sets the internal adjacencies (adjacencies between children of the same prototile) shown in Figure 2, and the external adjacencies based on the rules defined in Table 1. This is done by determining the neighbors of the tile and following the appropriate rule. Now genPenTiling.m outputs the set of newTiles which make up the next iteration of the tiling.
3. The Laplacian and Spectrum Let us first present some relevant definitions and background information related to the Laplacian matrix and its spectrum. Let G =(V,E)beasimple(undirected edges with no looped vertices) graph with vertex set V = v1,v2 ...,vn and edge set E = e ,e ,...,e , where e is an unordered element{ of V V. } { 1 2 m} i ⇥ Definition 5. The adjacency matrix, A(G), is an n n matrix such that the ⇥ entry aij is the number of edges connecting vertex i to vertex j.
Definition 6. The degree matrix, D(G),isann n matrix such that the diagonal ⇥ entry aii is the degree, the number of edges of the vertex vi.
Definition 7. [10] The Laplacian matrix of G, (G), is D(G) A(G). 8 M.DAIRYKO,C.HOFFMAN,J.PATTYSON,ANDH.PECK
Definition 8. [3] In a graph G, let u and v be vertices and dv the degree of vertex v. An equivalent definition of the Laplacian matrix (G) for u, v V is 2 dv if u = v, (u, v)= 1 if u and v are adjacent, 8 <>0 otherwise.
Definition 9. Let f : V (G) :>R be a function that assigns a real value to each vertex of the graph G.TheLaplacian! operator is a function acting on f, defined by f(v)= f(v) f(w) w:d(Xv,w)=1 which sums the di↵erence of the real values of adjacent nodes. Remark 2. The above definition does not require that G has finitely many vertices.
Definition 10. [3] Let 0 1 ... n 1 be the set of eigenvalues of (G). This set is denoted as ({ (G)), and is called the} spectrum of (G). In the following example, we compute (G) and ( (G)) for a given graph G. Example 1. Let G be a simple graph on six vertices and seven edges, as shown below.
6 4 5
3 2 1
Using Definition 7 of the Laplacian, we find the degree matrix and the adjacency matrix to compute (G). This is equivalent to using the piecewise function as shown above. Once (G) is constructed, we can compute ( (G)).
200000 010010 030000 101010 0 0020001 0 0101001 D(G)= A(G)= B 000300C B 001011C B C B C B 000030C B 110100C B C B C B 000001C B 000100C B C B C @ A @ A 2 10010 13 10 10 0 0 12 100 1 (G)= B 00131 1 C B C B 1 10130C B C B 000101C B C Using MatLab, we find@ the spectrum A ( (G)) = 0.0, 0.7216, 1.6826, 3.0, 3.7046, 4.8912 . { } ON THE SPECTRUM OF THE PENROSE LAPLACIAN 9
Many of the physical properties of quasicrystals can be revealed through study of the integrated density of states function, which we approach by finding the cumulative distribution function of the spectrum of finite iterations of the Penrose tiling. Definition 11. The cumulative distribution function (CDF) of a random variable X is given by F (x)=P [X x] for all x R 2 where P denotes the probability function. Definition 12. The support of a CDF is the smallest closed set whose complement has probability zero. The CDF gives the graph representing the probability that a random is less than or equal to a certain real-valued number, labeled along the x-axis. 3.1. Computational Methods of the Laplacian and Spectrum. This section discusses the MatLab code which we use to generate the spectrum of the Laplacian matrix and the CDF of finite iterations of the spectrum. See Appendix A.2, A.3, and A.4 for the code and more detailed comments. The function tilingLaplacian.m takes input newTiles (which is a vector of tile structures) from the genPenTiling.m function and outputs the eigenvalues of the Laplacian matrix of that iteration of the tiling. First, we construct the adjacency matrix for the tiling based on the substitution with four prototiles. We search through all of the new tiles and check if each tile has a neighbor on the first, second, and third side. If the tile does have a neighbor, we put a ‘1’ in the row and column corresponding to the tile and its neighbor tile. This builds the adjacency matrix for the tiling based o↵of the set of four Robinson triangles. Next, the function finds pairs of tiles that are neighbors on the third side and records their id numbers. If two tiles are neighbors on the third side, we combine them into one tile to return to the original two Penrose prototiles as shown in Figure 1. To do this in the adjacency matrix, we add corresponding rows and columns and delete the extra rows and columns. For example, if i and j (with i