ON THE SPECTRUM OF THE PENROSE LAPLACIAN MICHAEL DAIRYKO, CHRISTINE HOFFMAN, JULIE PATTYSON, AND HAILEE PECK

Transcription

1 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 particular 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 structures 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 Substitution Method Rules for Substitution with Robinson Triangles Computational Methods for Penrose Tiling 7 3. The Laplacian and Spectrum Computational Methods of the Laplacian and Spectrum 9 4. The Hausdor Dimension Computational Methods for The Hausdor Dimension Results Conclusion and Future Work 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 Introduction Before 1982, crystalline structures were defined as having periodic lattice structure. However, on April 8, 1982, Dan Shechtman discovered an aluminum-manganese quasicrystalline structure that had long-range order but did not exhibit the periodic patterns that characterized crystals [12]. Instead, the structure had icosahedral Date: August1,2013. This work was supported by NSF grant DMS

2 2 M. DAIRYKO, C. HOFFMAN, J. PATTYSON, AND H. PECK symmetry that was not possible, by definition, in crystal form. In 2011, Shechtman 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 nonperiodically. 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.

3 ON THE SPECTRUM OF THE PENROSE LAPLACIAN 3 Definition 1. [7] Prototiles are a set of finite inequivalent tiles (i.e. equivalent under rigid motions, expansions, or contractions). are not Definition 2. [7] A tiling is an arrangement of tiles, such that their union covers and packs R 2 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 subdivides 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] 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 each triangle prototile by 1+p 5 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 4 M. DAIRYKO, C. HOFFMAN, J. PATTYSON, AND H. 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] Rules for Substitution with Robinson Triangles. The sides of each triangle 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

5 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 between 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 distinct 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 B 2, W 2, and P 2. The children of the Y parent prototile are P 3 and Y 3. For the external adjacencies, the child Y 3 is adjacent to the child W 2 on the third side (superscript 3). Also, the child P 3 is adjacent to the child B 2 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 6 M. DAIRYKO, C. HOFFMAN, J. PATTYSON, AND H. PECK Table 1. Rules for Shared Edges of Adjacent Tiles Rule Shared Edge External Adjacencies adj P 3 B 3 1 P1 3 B2, 3 B1 1 P2 1 B 2 W 2 P 2 P 1 Y 1 B 1 B 1 adj P 1 2 P 2 1 B 2 2, Y 3 1 W 3 2 P 1 Y 1 B 1 B 2 W 2 P 2 Y 3 adj W 3 3 Y 2 3 W 2 4 P 3 Y 3 W 4 B 4 P 2 adj B 2 4 B 3 1 P 3 2 B 2 W 2 P 2 P 1 Y 1 B 1 P 1 adj Y 1 5 Y 3 3 W 3 2, P 2 3 B 2 2 B 2 W 2 P 2 P 3 Y 3 B 1 adj W 1 6 Y 3 1 W 3 4, P 2 1 B 2 4 P 1 Y 1 B 1 W 4 B 4 P 2 adj W 2 7 P 3 2 B 3 4 B 2 W 2 P 2 W 4 B 4 Y 2 adj W 2 8 P 3 3 B 3 4 Y 3 P 3 W 4 B 4 B 2 adj Y 2 9 B 3 1 P 3 3 P 1 Y 1 B 1 P 3 Y 3

7 ON THE SPECTRUM OF THE PENROSE LAPLACIAN Computational Methods for Penrose Tiling 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 = {v 1,v 2...,v n } and edge set E = {e 1,e 2,...,e m }, where e i is an unordered element of V V. Definition 5. The adjacency matrix, A(G), is an n n matrix such that the entry a ij is the number of edges connecting vertex i to vertex j. Definition 6. The degree matrix, D(G), is an n n matrix such that the diagonal entry a ii is the degree, the number of edges of the vertex v i. Definition 7. [10] The Laplacian matrix of G, (G), is D(G) A(G).

8 8 M. DAIRYKO, C. HOFFMAN, J. PATTYSON, AND H. PECK Definition 8. [3] In a graph G, let u and v be vertices and d v the degree of vertex v. An equivalent definition of the Laplacian matrix (G) for u, v 2 V is 8 >< d v if u = v, (u, v) = 1 if u and v are adjacent, >: 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 X f(v) = f(v) f(w) w:d(v,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 apple 1 apple...apple 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 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)). 0 D(G) = (G) = 1 C A Using MatLab, we find the spectrum 0 A(G) = C A ( (G)) = {0.0, , , 3.0, , }. 1 C A

9 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 where P denotes the probability function. F (x) =P [X apple x] for all x 2 R 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 than or equal to a certain real-valued number, labeled along the x-axis. is less 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<j) are two tiles that are neighbors on the third side, the function finds (i, j) as a pair. In the adjacency matrix created above, add row j to row i and add column j to column i. Then we delete row j and column j. Now row i and column i correspond to (i, j). After doing this for each pair, we change the diagonal of the adjacency matrix to all zeros. The function then constructs the degree matrix by adding the entries in each row from the adjacency matrix and entering the sum on the diagonal of the degree matrix. Now the diagonal entry (i, i) for each vertex i contains its degree. The rest of the degree matrix is filled with zeros. The Laplacian matrix is constructed by taking D(G) A(G) (by Definition 7). Then tilinglaplacian.m outputs the set of eigenvalues of the Laplacian. Next, the plotmat.m function takes input of a cell array of eigenvalues (which are output of tilinglaplacian.m) and outputs the eigenvalues of each iteration on a plot. This displays the spectrum for the finite iterations of the Penrose Laplacian (see Figure 5). The cdf.m function inputs a vector of eigenvalues found from the Laplacian and constructs a CDF. At each point on the horizontal axis that is an eigenvalue, the

10 10 M. DAIRYKO, C. HOFFMAN, J. PATTYSON, AND H. PECK function will jump vertically by 1 n where n is the number of eigenvalues in the input vector (see Figure 6). 4. The Hausdorff Dimension We begin with some preliminary definitions. Definition 13. Let X R n be nonempty, and define U := sup { x We say {U i } is a -cover of X if S (1) X 1 U i i=1 (2) 0 apple U i apple y : x, y 2 U}. Definition 14. [6, Equations 2.1, 2.2] Let X R n be nonempty and s be any nonnegative real number. For >0 define ( 1 ) X H s (X) :=inf U i s : U i is a -cover of X. i=1 The s-dimensional Hausdor measure of X, denoted H s (X) is lim!0 H s (X). Observe that for a fixed <1, H s (X) is non-increasing with respect to s. We see that in Figure 2, H s (X) is always either 0 or 1, except for at most one point. In the following example, this point will be (1, 0). We call this point the Hausdor dimension. Definition 15. Let X R n be nonempty and s be any nonnegative real number. For >0 define dim H (X) :=inf{s : H s (X) =0} =sup{s : H s (X) =1}. We call dim H (X) the Hausdor dimension of X. Therefore, we will have a single point, the Hausdor dimension, at which the measure is not zero or infinity (see Example 2). Example 2. Let us divide the interval A =[0, 1] into n closed subintervals of equal length. Let U i := i n, i+1 n, where i = {0, 1,...,n 1}. Therefore {Ui } is a -cover of A. Let X s s Ui s = n = n n be an approximation of H s (A). As! 0, that is n!1we have 8 >< 0 if s>1, H s (A) = 1 if s =1 >: 1 otherwise. This indicates that the Hausdor dimension of A is 1 which agrees with our intuition about the topological dimension of A.

11 ON THE SPECTRUM OF THE PENROSE LAPLACIAN 11 H s (X) Figure 4. Hausdor Dimension for [0, 1] with = 1 n s gleaned from [6] 4.1. Computational Methods for The Hausdor Dimension. This section discusses the MatLab code which generates an estimate of the Hausdor Dimension. See Appendix A.5 for the code and more detailed comments. The ndstable.m is a function that takes input of a set of eigenvalues of the Laplacian and an estimate s of the Hausdor dimension and creates a column of values. These columns are calculated by n s,where = 1 2 for i 2{1, 2,...,7}, n i is the number of -intervals needed to contain all eigenvalues, and s is an estimate for the Hausdor dimension (a value between 0 and 1). This column of values is then used to estimate the Hausdor dimension of the spectrum. As all of the eigenvalues of the spectrum of the Laplacian were found to be between the finite interval [0, 8], the function first divides this interval into partitions of length. Each row in the output 10 1 vector corresponds to a unique value of. The function ndstable.m then finds the number of eigenvalues that lie within each of these partitions and the number of -intervals needed to contain all eigenvalues, which is the number of -partitions containing at least one eigenvalue. Each row output of this function is n( 1 2 ) s,wherenisthenumber of -intervals i needed for the corresponding = 1 2. If the output 10 1 column vector increases i (respectively decreases), the Hausdor dimension of the spectrum is greater than (respectively less than) the s value. 5. Results In this section, we highlight the results of the spectrum and probability functions for the tiling. Using tilinglaplacian.m we find the eigenvalues of the Laplacian matrices of seven iterations of the Penrose tiling. Then with plotmat.m, we plot the eigenvalues for the first seven iterations as shown in Figure 5. From [13] we know that the spectrum of the Penrose Laplacian should be bounded by twice the highest degree of any node. As we are working with the Penrose rhombus and diamond which are four-sided tiles, the highest possible degree of a node is four. Thus, the spectrum should be bounded by eight. As shown in Figure 5, this is illustrated by our results on the first seven iterations.

12 12 M. DAIRYKO, C. HOFFMAN, J. PATTYSON, AND H. PECK Iteration of Penrose Tiling 8 Spectrum of Finite Iterations of Penrose Tiling Eigenvalues Figure 5. Eigenvalues of the first seven iterations The program cdf.m constructs the CDF of eigenvalues from each iteration. Since we have a small set of eigenvalues for the first few iterations of the tiling, this function is step-like and jumps at each discrete eigenvalue. Once we found the CDF of finite iterations of the Penrose tiling, Theorem 3.2 from [4] tells us that, in the infinity norm, these functions approach a limiting function. This limiting function will be the integrated density of states of the Penrose Laplacian, and the support of this limiting function is the spectrum of the Penrose Laplacian. In the finite approximations, the support of the CDF is the set of jumps in the graph. These jumps occur at the set of eigenvalues, which for the finite approximations is the Penrose Laplacian spectrum. Thus, in the infinite case, the interval where the function is increasing is the support of the function, and this support is the spectrum. Table 2 will help to estimate the Hausdor dimension of the Penrose Laplacian spectrum. Theory says that, as approaches 0, H s will be either zero or infinity. Moving down through the columns of values in Table 2 is representative of taking to 0. Therefore, if the values increase down the column, we can definitively say the estimated s is not the Hausdor dimension, as the values are going to infinity. If the values decrease down the column, the chosen s is again not the Hausdor dimension, as the values are going to zero. However, the various values of s in which the columns neither clearly decrease nor increase provide an interval in which the actual Hausdor dimension can be found. Therefore, we estimate the Hausdor dimension of the Penrose Laplacian spectrum to be between.85 and.99, which implies fractal structure.

13 ON THE SPECTRUM OF THE PENROSE LAPLACIAN 13 1 Cumulative Distribution Function 0.8 Probability Real Numbers Figure 6. Cumulative distribution plot for the first seven iterations Table 2. Values of n s as Functions of i and s s i Conclusion and Future Work Shechtman s 1982 discovery of a quasicrystal sparked research interest in the field of quasicrystalline structures as they relate to aperiodic tilings. We are interested in the spectrum of the Penrose tiling as it reveals certain properties of quasicrystals. We created an algorithm to generate finite iterations of the Penrose tiling. Then we were able to create the Laplacian matrix for each respective finite tiling. With this, we calculated the spectrum of each Laplacian matrix. Using the spectrum, we estimated the Hausdor dimension of the Penrose tiling. For future work, we would like to explore ways to generate more iterations of the Penrose tiling to improve the estimate of the Hausdor dimension. We hope to use the methods established here to generalize these results to other aperiodic tilings.

14 14 M. DAIRYKO, C. HOFFMAN, J. PATTYSON, AND H. PECK 7. Acknowledgments Figures 1, 2, and 3 related to the substitution method were inspired by and modified from images in [7]. The presentation of Table 2 was inspired by [2]. We would like to thank Dr. May Mei and Drew Zemke for their help with this project. We would also like to thank the Summer Math Institute and the Mathematics Department at Cornell University for the use of their resources. This work was supported by NSF grant DMS

15 ON THE SPECTRUM OF THE PENROSE LAPLACIAN 15 A.1. genpentiling.m. Appendix A: MatLab Code 1 % genpentiling.m 2 % Michael Dairyko, Christine Hoffman, Julie Pattyson, Hailee Peck 3 % 4 % Given the data for an iteration of the Penrose tiling, this function 5 % outputs the next iteration. 6 % A tiling is stored as a list of structures ; each has the following fields : 7 % id A positive integer that is unique to this tile in this 8 % iteration 9 % color p, b, y or w, based on whether the tile is pink, 10 % blue, yellow or white 11 % children An array that contains this tile s children s ids. 12 % n1, n2, n3 A list of this tile s neighbors on each side function [ newtiles ] = genpentiling ( oldtiles ) % First we create the new tiles without worrying 17 % about adjacencies. We will keep track of the old tiles children 18 % and each new tile s parent. 19 % Recall the rules : 1) each pink tile turns into a blue tile, a white tile 20 % and a pink tile, 2) each blue tile turns into a pink, a yellow, and a blue 21 % tile, 3) each white tile turns into a white tile and a blue tile, 22 % and 4) each yellow tile turns into a pink tile and a yellow tile 23 newtiles = {}; 24 idcounter = 1; % Setting up an empty temptiles set which will start empty but 27 % each time we create a new tile, will add that new tile to this set. 28 % This is a master list of all tiles created during the process of generating 29 % the tiling. 30 temptiles = {}; for i = 1: length ( oldtiles ) 33 % Go through old tiles one by one. 34 thistile = oldtiles{ i } ; % The following if statements set up the children of each of these 37 % tiles and their internal adjacencies % thistile is a pink tile. 40 % newb2, neww2, and newp2 are thistile s children. 41 if thistile. color == p % Tile 1 ( first child ) := blue 44 newb2 = struct ( id, idcounter, color, b ); 45 % Assign an index number so we can reference the tile by this number. 46 idcounter = idcounter + 1; % Tile 2 ( second child ) := white 49 neww2 = struct ( id, idcounter, color, w ); 50 idcounter = idcounter + 1; % Tile 3 ( third child ) := pink 53 newp2 = struct ( id, idcounter, color, p ); 54 idcounter = idcounter + 1; % Note that we labeled tiles coming from our pink 57 % tiles with a 2 ( because pink is triangle 2), in order 58 % to distinguish the pink, blue, white tiles that 59 % come from pink apart from the pink, blue, white that come from 60 % other parent tiles % Set these new tiles as children of thistile. 63 thistile. children = [newb2. id, neww2. id, newp2. id ]; % Internal adjacencies of p ( see Figure #5 in the paper for a 66 % visual of these ) : 67 % Blue connects to white on side newb2. n1 = neww2. id ; 69 % White connects to pink on side 2.

16 16 M. DAIRYKO, C. HOFFMAN, J. PATTYSON, AND H. PECK 70 neww2. n2 = newp2. id ; 71 % Likewise, pink connects to white on side newp2. n2 = neww2. id ; 73 % White connects to blue on side neww2. n1 = newb2. id ; % We must define each adjacency both ways so the 77 % adjacency is stored in each tile ( and not just in one tile ) 78 % because we want our adjacency matrix to reflect the symmetry 79 % of the adjacencies %Now add these children to our set of newtiles. 82 newtiles{ length ( newtiles )+1} = newb2 ; 83 newtiles{ length ( newtiles )+1} = neww2 ; 84 newtiles{ length ( newtiles )+1} = newp2 ; % thistile is a blue tile. 87 % newp1, newy1, and newb1 are thistile s children. 88 elseif thistile. color == b % As with the pink tile, we label these children with a 1 to distinguish 91 % that they came from the blue tile ( triangle 1) % Tile 1 ( first child ) := Pink 94 newp1 = struct ( id, idcounter, color, p ); 95 idcounter = idcounter + 1; % Tile 2 ( second child ) := Yellow 98 newy1 = struct ( id, idcounter, color, y ); 99 idcounter = idcounter + 1; % and Tile 3 ( third child ) := Blue 102 newb1 = struct ( id, idcounter, color, b ); 103 idcounter = idcounter + 1; % Set these new tiles as children of thistile. 106 thistile. children=[newp1. id, newy1. id, newb1. id ]; %Internal adjacencies of b ( see Figure 5) : 109 % Pink connects to yellow on side newp1. n1 = newy1. id ; 111 % Yellow connects to blue on side newy1. n2 = newb1. id ; 113 % Blue connects to yellow on side newb1. n2 = newy1. id ; 115 % Yellow connects to pink on side newy1. n1 = newp1. id ; % Now add these children to our set of newtiles. 119 newtiles{ length ( newtiles )+1} = newp1 ; 120 newtiles{ length ( newtiles )+1} = newy1 ; 121 newtiles{ length ( newtiles )+1} = newb1 ; % thistile is yellow a tile. 124 % newp3 and newy3 are thistile s children. 125 elseif thistile. color == y % As with the pink and blue tiles, we label these children with a 3 to distinguish 128 % that they came from the yellow tile ( triangle 3) % Tile 1 ( first child ) := Pink 131 newp3 = struct ( id, idcounter, color, p ); 132 idcounter = idcounter + 1; % Tile 2 ( second child ) := Yellow 135 newy3 = struct ( id, idcounter, color, y ); 136 idcounter = idcounter + 1; % Set these new tiles as children of thistile. 139 thistile. children = [newp3. id, newy3. id ];

17 ON THE SPECTRUM OF THE PENROSE LAPLACIAN % Internal adjacencies of y ( see Figure 5) : 142 % Pink connects to yellow on side newp3. n1 = newy3. id ; 144 % Yellow connects to pink on side newy3. n1 = newp3. id ; % Now add these children to our set of newtiles. 148 newtiles{ length ( newtiles )+1} = newp3 ; 149 newtiles{ length ( newtiles )+1} = newy3 ; % thistile is a white tile. 152 % newb4 and neww4 are thistile s children. 153 else % As with the previous tiles, we label these children with a 4 to distinguish 156 % that they came from the white tile ( triangle 4) % Tile 1 ( first child ) := Blue 159 newb4 = struct ( id, idcounter, color, b ); 160 idcounter = idcounter + 1; % Tile 2 ( second child ) := White 163 neww4 = struct ( id, idcounter, color, w ); 164 idcounter = idcounter + 1; % Set these new tiles as children of thistile. 167 thistile. children = [newb4. id, neww4. id ]; % Internal adjacencies of w ( see Figure 5) : 170 % Blue connects to white on side newb4. n1 = neww4. id ; 172 % White connects to blue on side neww4. n1 = newb4. id ; % Now add these children to our set of newtiles. 176 newtiles{ length ( newtiles )+1} = newb4 ; 177 newtiles{ length ( newtiles )+1} = neww4 ; 178 end % Store the tile that we just processed and changed. 181 temptiles{ length ( temptiles )+1 } = thistile ; 182 end % Now we set up the external adjacencies. We proceed 185 % by iterating through the temptiles which is where the old tiles that 186 % were changed in the previous if statement are stored. Note 187 % that when we reference an adjacency rule, we are referencing the 188 % numbered rules from the table in Appendix B of the paper % Go through the temptiles one by one. 191 for i = 1: length ( temptiles ) %The tile we are currently working with is thistile. 194 thistile = temptiles{ i } ; % thistile is pink. 197 if thistile. color == p % Checks if our pink tile has a neighbor on the 1 side (1 arrow ). 200 % If it does it will be blue or yellow, according to rule 2 or rule if isfield ( thistile, n1 ) % This is the neighbor on the 1 side of this pink tile, denoted by N. 204 N = temptiles{ thistile. n1 } ; % In Rule 2, the neighbor is a blue tile. 207 if N. color == b % Our neighbor N is a blue tile, so follow Rule 2.

18 18 M. DAIRYKO, C. HOFFMAN, J. PATTYSON, AND H. PECK 210 % Referencing the first child ( B2) of the pink tile. 211 Child1 = newtiles{ thistile. children (1) } ; 212 % B2 is connected on the second side to the first child ( P1 ) 213 % of the neighbor ( blue ) tile. 214 Child1. n2 = newtiles{n. children (1) }. id ; % Again, we define each adjacency both ways so the 217 % adjacency is stored in each tile ( and not just in one tile ) 218 % because we want our adjacency matrix to reflect the symmetry 219 % of the adjacencies. 220 newtiles{n. children (1) }. n2 = Child1. id ; % Referencing the second child (W2) of the pink tile. 223 Child2 = newtiles{ thistile. children (2) } ; 224 % W2 is connected on the third side to the second child (Y1) 225 % of the neighbor ( blue ) tile. 226 Child2. n3 = newtiles{n. children (2) }. id ; 227 newtiles{n. children (2) }. n3 = Child2. id ; % In Rule 5, the neighbor is a yellow tile. 230 elseif N. color == y 231 % Our neighbor N is a yellow tile, so follow Rule % Referencing the first child ( B2 ) of the pink tile. 233 Child1 = newtiles{ thistile. children (1) } ; 234 % B2 is connected on the second side to the first child ( P3 ) of 235 % the neighbor ( yellow ) tile. 236 Child1. n2 = newtiles{n. children (1) }. id ; 237 newtiles{n. children (1) }. n2 = Child1. id ; % Referencing the second child (W2) of the pink tile. 240 Child2 = newtiles{ thistile. children (2) } ; 241 % W2 connected on the third side to the second child (Y3) 242 % of the neighbor ( yellow ) tile 243 Child2. n3 = newtiles{n. children (2) }. id ; 244 newtiles{n. children (2) }. n3 = Child2. id ; 245 end 246 end % Checks if our pink tile has a neighbor on the 2 side (2 arrows ). 249 % If it does it will be blue or white, according to Rule 4 or Rule % In Rule 4, the neighbor is a blue tile. 251 % In Rule 7, the neighbor is a white tile. 252 if isfield ( thistile, n2 ) % This is the neighbor on the 2 side of this pink tile. 255 N = temptiles{ thistile. n2 } ; 256 % Our neighbor N is a blue tile, so follow Rule if N. color == b 258 % Referencing the third child ( P2 ) of the pink tile. 259 Child3 = newtiles{ thistile. children (3) } ; 260 % P2 is connected on the third side to the third child ( B1) 261 % of the neighbor ( blue) tile. 262 Child3. n3 = newtiles{n. children (3) }. id ; 263 newtiles{n. children (3) }. n3 = Child3. id ; % Our neighbor is a white tile, so follow Rule elseif N. color == w 267 % Referencing the third child ( P2 ) of the pink tile. 268 Child3 = newtiles{ thistile. children (3) } ; 269 % P2 is connected on the third side to the first child ( B4) 270 % of the neighbor ( blue ) tile. 271 Child3. n3 = newtiles{n. children (1) }. id ; 272 newtiles{n. children (1) }. n3 = Child3. id ; 273 end 274 end % Checks if our pink tile has a neighbor on the 3 side ( no arrows ). 277 % If it does it will be blue, according to Rule 1.

19 ON THE SPECTRUM OF THE PENROSE LAPLACIAN % In Rule 1, the neighbor is a blue tile. 279 if isfield ( thistile, n3 ) %This is the neighbor on the 3 side of this pink tile ( which will always be blue ). 282 N = temptiles{ thistile. n3 } ; 283 % Referencing the first child ( B2 ) of the pink tile. 284 Child4 = newtiles{ thistile. children (1) } ; 285 % B2 is connected on the third side to the first child ( P1 ) 286 % of the neighbor ( blue ) tile. 287 Child4. n3 = newtiles{n. children (1) }. id ; 288 newtiles{n. children (1) }. n3 = Child4. id ; % Referencing the third child ( P2 ) of the pink tile. 291 Child5 = newtiles{ thistile. children (3) } ; 292 % P2 is connected to on the first side to the third child ( B1 ) 293 % of the neighbor ( blue ) tile. 294 Child5. n1 = newtiles{n. children (3) }. id ; 295 newtiles{n. children (3) }. n1 = Child5. id ; 296 end 297 end %Now we run through the external adjacencies for tiles other than pink. 300 % thistile is blue. 301 if thistile. color == b % Checks if the blue tile has a neighbor on the first (1 arrow ) side. 304 if isfield ( thistile, n1 ) % This is the neighbor on the 1 side of the blue tile. 307 % Follow Rule 2 or Rule 6: 308 % In Rule 2, the neighbor is a pink tile. 309 % In Rule 6, the neighbor is a white tile. 310 N = temptiles{ thistile. n1 } ; 311 % Follow Rule if N. color == p 313 % Referencing the first child ( P1 ) of the blue tile. 314 Child6 = newtiles{ thistile. children (1) } ; 315 % P1 is connected on the second side to the first child ( B2) 316 % of the neighbor ( pink) tile. 317 Child6. n2 = newtiles{n. children (1) }. id ; 318 newtiles{n. children (1) }. n2 = Child6. id ; Child7 = newtiles{ thistile. children (2) } ; 321 % Referencing the second child (Y1) of the blue tile. 322 Child7. n3 = newtiles{n. children (2) }. id ; 323 % Y1 is conncetd on the third side to the second child ( W2) 324 % of the neighbor ( pink) tile. 325 newtiles{n. children (2) }. n3 = Child7. id ; % Follow Rule elseif N. color == w 329 % Referencing the first child ( P1 ) of the blue tile. 330 Child6 = newtiles{ thistile. children (1) } ; 331 % P1 is connected on the second side to the first child ( B4) 332 % of the neighbor ( white ) tile. 333 Child6. n2 = newtiles{n. children (1) }. id ; 334 newtiles{n. children (1) }. n2 = Child6. id ; 335 % Referencing the second child (Y1) of the blue tile. 336 Child7 = newtiles{ thistile. children (2) } ; 337 % Y1 is connected on the third side to the second child ( W4) 338 % of the neighbor ( white ) tile. 339 Child7. n3 = newtiles{n. children (2) }. id ; 340 newtiles{n. children (2) }. n3 = Child7. id ; 341 end 342 end

20 20 M. DAIRYKO, C. HOFFMAN, J. PATTYSON, AND H. PECK % Checks if the blue tile has a neighbor on the second (2 arrows ) side. 345 if isfield ( thistile, n2 ) % This is the neighbor on the 2 side of the blue tile. 348 % Follow Rule 4 or Rule 9: 349 % In Rule 4, the neighbor is a pink tile. 350 % In Rule 9, the neighbor is a yellow tile. 351 N = temptiles{ thistile. n2 } ; 352 % Follow Rule if N. color == p % Referencing the third child ( B1) of the blue tile. 356 Child8 = newtiles{ thistile. children (3) } ; 357 % B1 is connected on the third side to the third child ( P2 ) 358 % of the neighbor ( pink) tile. 359 Child8. n3 = newtiles{n. children (3) }. id ; 360 newtiles{n. children (3) }. n3 = Child8. id ; 361 % Follow Rule elseif N. color == y 363 % Referencing the third child ( B1) of the blue tile. 364 Child8 = newtiles{ thistile. children (3) } ; 365 % B1 is connected on the third side to the first child ( P3 ) 366 % of the neighbor ( yellow ) tile 367 Child8. n3 = newtiles{n. children (1) }. id ; 368 newtiles{n. children (1) }. n3 = Child8. id ; 369 end 370 end % Checks if the blue tile has a neighbor on the 3 ( no arrows ) side. 373 if isfield ( thistile, n3 ) % This is the neighbor on the 3 side of the blue tile ( will always be pink ). 376 % Follow Rule 1: In Rule 1, the neighbor is a pink tile. 377 N = temptiles{ thistile. n3 } ; % Referencing the first child ( P1 ) of the blue tile. 380 Child9 = newtiles{ thistile. children (1) } ; 381 % P1 is connected on the third side to the first child ( B2) 382 % of the neighbor ( pink) tile. 383 Child9. n3 = newtiles{n. children (1) }. id ; 384 newtiles{n. children (1) }. n3 = Child9. id ; 385 % Referencing the third child ( B1) of the blue tile. 386 Child10 = newtiles{ thistile. children (3) } ; 387 % B1 is connected on the first side to the third child ( P2 ) 388 % of the neighbor ( pink) tile. 389 Child10. n1 = newtiles{n. children (3) }. id ; 390 newtiles{n. children (3) }. n1 = Child10. id ; 391 end 392 end % thistile is white. 395 if thistile. color == w % Checks if the white tile has a neighbor on the first (1 arrow ) side. 398 if isfield ( thistile, n1 ) % This is the neighbor on the 1 side of the white tile ( will always be blue ). 401 % Follow Rule 6. In Rule 6, the neighbor is a blue tile. 402 N = temptiles{ thistile. n1 } ; 403 % Referencing the first child ( B4 ) of the white tile. 404 Child11 = newtiles{ thistile. children (1) } ; 405 % B4 is connected on the second side to the first child ( P1 )

21 ON THE SPECTRUM OF THE PENROSE LAPLACIAN % of the neighbor ( blue ) tile. 407 Child11. n2 = newtiles{n. children (1) }. id ; 408 newtiles{n. children (1) }. n2 = Child11. id ; 409 % Referencing the second child (W4) of the white tile. 410 Child12 = newtiles{ thistile. children (2) } ; 411 % W4 is connected on the third side to the second child ( Y1) 412 % of the neighbor ( blue ) tile. 413 Child12. n3 = newtiles{n. children (2) }. id ; 414 newtiles{n. children (2) }. n3 = Child12. id ; 415 end % Checks if the white tile has a neighbor on the second (2 arrow ) side. 418 if isfield ( thistile, n2 ) 419 % This is the neighbor on the 2 side of the white tile. 420 % Follow Rule 7 or Rule N = temptiles{ thistile. n2 } ; % In Rule 7, the neighbor is a pink tile. 424 if N. color == p 425 % Follow Rule % Referencing the first child ( B4 ) of the white tile. 427 Child13 = newtiles{ thistile. children (1) } ; 428 % B4 is connected on the third side to the third child ( P2 ) 429 % of the neighbor ( pink) tile 430 Child13. n3 = newtiles{n. children (3) }. id ; 431 newtiles{n. children (3) }. n3 = Child13. id ; % In Rule 8, the neighbor is a yellow tile. 434 elseif N. color == y 435 % Follow Rule % Referencing the first child ( B4 ) of the white tile. 437 Child13 = newtiles{ thistile. children (1) } ; 438 % B4 is connected on the third side to the first child ( P3 ) 439 % of the neighbor ( yellow ) tile. 440 Child13. n3 = newtiles{n. children (1) }. id ; 441 newtiles{n. children (1) }. n3 = Child13. id ; 442 end 443 end % Checks if the white tile has a neighbor on the third ( no arrows ) side. 446 if isfield ( thistile, n3 ) % This is the neighbor on the 3 side of the white tile ( will always be yellow ). 449 % Follow Rule 3. In Rule 3, the neighbor is a yellow tile. 450 N = temptiles{ thistile. n3 } ; % Referencing the second child (W4) of the white tile. 453 Child14 = newtiles{ thistile. children (2) } ; 454 % W4 is connected on the second side to the second child (Y3) 455 % of the neighbor ( yellow ) tile. 456 Child14. n2 = newtiles{n. children (2) }. id ; 457 newtiles{n. children (2) }. n2 = Child14. id ; 458 end 459 end % thistile is yellow. 462 if thistile. color == y % Checks if the yellow tile has a neighbor on the first (1 arrow ) side. 465 if isfield ( thistile, n1 ) % This is the neighbor on the first side of the yellow tile ( will always be pink ). 468 % Follow Rule 5. In Rule 5, the neighbor is a pink tile.

22 22 M. DAIRYKO, C. HOFFMAN, J. PATTYSON, AND H. PECK 469 N = temptiles{ thistile. n1 } ; % Referencing the first child ( P3 ) of the yellow tile. 472 Child15 = newtiles{ thistile. children (1) } ; 473 % P3 is connected on the second side to the first child ( B2) 474 % of the neighbor ( pink) tile. 475 Child15. n2 = newtiles{n. children (1) }. id ; 476 newtiles{n. children (1) }. n2 = Child15. id ; 477 %Referencing the second child (Y3) of the yellow tile. 478 Child16 = newtiles{ thistile. children (2) } ; 479 % Y3 is connected on the third side to the second child ( W2) 480 % of the neighbor ( pink) tile. 481 Child16. n3 = newtiles{n. children (2) }. id ; 482 newtiles{n. children (2) }. n3 = Child16. id ; 483 end % Checks if the yellow tile has a neighbor on the second (2 arrows ) side. 486 if isfield ( thistile, n2 ) % This is the neighbor on the second side of the yellow tile. 489 % Follow Rule 8 or Rule N = temptiles{ thistile. n2 } ; % In Rule 8, the neighbor is a white tile. 493 if N. color == w 494 % Follow Rule % Referencing the first child ( P3 ) of the yellow tile. 496 Child17 = newtiles{ thistile. children (1) } ; 497 % P3 is connected on the third side to the first child ( B4) 498 % of the neighbor ( white ) tile. 499 Child17. n3 = newtiles{n. children (1) }. id ; 500 newtiles{n. children (1) }. n3 = Child17. id ; % In Rule 9, the neighbor is a blue tile. 503 elseif N. color == b 504 % Follow Rule % Referencing the first child ( P3 ) of the yellow tile. 506 Child18 = newtiles{ thistile. children (1) } ; 507 % P3 is connected on the third side to the third child ( B1) 508 % of the neighbor ( blue ) tile. 509 Child18. n3 = newtiles{n. children (3) }. id ; 510 newtiles{n. children (3) }. n3 = Child18. id ; 511 end 512 end % Checks if the yellow tile has a neighbor on the third ( no arrows ) side. 515 if isfield ( thistile, n3 ) % This is the neighbor on the third side of the yellow tiley ( will always be white ). 518 N = temptiles{ thistile. n3 } ; 519 % Follow Rule 3. In Rule 3, the neighbor is a white tile. 520 % Referencing the second child (Y3) of the yellow tile. 521 Child19 = newtiles{ thistile. children (2) } ; 522 % Y3 is connected on the second side to the second child (W4) 523 % of the neighbor ( white ) tile. 524 Child19. n2 = newtiles{n. children (2) }. id ; 525 newtiles{n. children (2) }. n2 = Child19. id ; 526 end 527 end 528 end A.2. tilinglaplacian.m. 1 % tilinglaplacian.m

23 ON THE SPECTRUM OF THE PENROSE LAPLACIAN 23 2 % Michael Dairyko, Christine Hoffman, Julie Pattyson, Hailee Peck 3 % 4 % This function takes the numbered newtiles generated by genpentiling as input, combines 5 % tiles that share a third side, and computes the Laplacian matrix of the dual 6 % graph of iterations of the tiling. Two numbered tiles are considered to be adjacent 7 % if they share a side. The function outputs the eigenvalues of this matrix, 8 % which is the spectrum of the Laplacian function [ EIG ] = tilinglaplacian ( newtiles ) % Here we create the adjacency matrix. 13 % Creates a matrix of zeros ( an nxn matrix with n = # of tiles in newtiles ) 14 A = zeros ( length ( newtiles )); 15 % Go through newtiles one by one. 16 for i = 1: length ( newtiles ) 17 % thistile = ith entry of newtiles = ith row of the matrix 18 thistile = newtiles{ i } ; 19 % Checks if thistile has a neighbor on the first side. 20 if isfield ( thistile, n1 ) 21 % Then put a 1 in the matrix in the row and 22 % column corresponding to thistile and the neighbor on the first side. 23 A( thistile. id, thistile. n1 ) = 1; 24 end 25 % Checks if thistile has a neighbor on the second side. 26 if isfield ( thistile, n2 ) 27 % Put a 1 in the row/ column corresponding to thistile and the neighbor. 28 A( thistile. id, thistile. n2 ) = 1; 29 end 30 % Checks if thistile has a neighbor on the third side. 31 if isfield ( thistile, n3 ) 32 % Put a 1 in the row/ column corresponding to thistile and the neighbor. 33 A( thistile. id, thistile. n3 ) = 1; 34 end 35 end % Here we find the ordered pairs containing the id numbers of each set of 38 % two tiles which are adjacent on their third side. We need to 39 % then condense the rows/ columns that correspond to these tiles in the adjacency 40 % matrix into one row/ column ( transforming our triangles back into 41 % rhombi/ diamonds ). 42 % Start off with an empty array, which we will fill. 43 pairs = {}; 44 for i = 1: length ( newtiles ) 45 thistile = newtiles{ i } ; 46 L = length ( pairs ); 47 % Checks if thistile has a neighbor on the third side. 48 if isfield ( thistile, n3 ) 49 % Increase the empty array by L = L + 1; 51 % Put in the ordered pair of thistile and its neighbor on the third side. 52 pairs{l} = [ thistile. id, thistile. n3 ]; 53 end 54 end % We need to get rid of repeated ordered pairs. These repeats 57 % occur because if tile1 is connected to tile2 on side 3, then tile2 is connected to 58 % tile1 on side 3, so we count the same pair twice. 59 % Start off with an empty array which we will fill with the ordered pairs without repeats. 60 OP = {}; 61 % Start with an identification number of ID = 1; 63 L2 = length (OP); 64 % Go through each pair of tiles which are neighbors on the third side. 65 for j = 1: length ( pairs )