A Stochastic Cellular Automaton for Three-Coloring Penrose Tiles Mark McClure Department of Mathematics University of North Carolina at Asheville Asheville, North Carolina 28804 email: [email protected] August 24, 2001 Abstract

We present a three state, stochastic cellular automaton that runs on Penrose tilings and seems to evolve to a three-colored equilibrium.

1 Introduction In 1973 and 1974, Roger Penrose discovered three sets of polygons each of which tiles the plane aperiodically and (if certain matching conditions are enforced) only aperiodically. Later, John H. Conway asked if such tilings can be three-colored, where adjacent tiles are to receive di erent colors. This question has been answered armatively for two types of Penrose tilings, but appears to be open for the remaining type. In this paper, we present an algorithm that seems to three-color nite parts of Penrose tilings of all types. The algorithm works by running a particular three-state, stochastic cellular automaton on a given Penrose tiling. The cellular automaton is chosen so that three-colorings are stable and it seems to generally evolve to such an equilibrium.

1

2 Penrose tilings There are three types of Penrose tilings: tilings by kites and darts, tilings by rhombs, and tilings by pentacles. We describe them brie y here. More detailed references are [2] and [3].

2.1 Tilings by kites and darts

Figure 1 shows the kite and dart. The sides have length either 1 or  , the golden ratio, and the angles are all integer multiples of 5 . The lled and un lled disks at the vertices are used to enforce a matching condition. When tiling the plane with kites and darts, we demand that lled disks meet lled disks and un lled disks meet un lled disks. This matching condition guarantees that any tiling by kites and darts will be aperiodic, i.e. no translation of the tiling maps each tile to another tile. Figure 2 shows part of such a tiling.

2.2 Tilings by rhombs

Figure 3 shows the fat and skinny rhombs. The sides all have length 1 and the angles are all integer multiples of 5 . The matching condition is slightly more complicated. We demand that lled disks meet lled disks, un lled disks meet un lled disks, and oriented edges meet with the correct orientation. Again, this matching condition guarantees that any tiling by rhombs will be aperiodic. Figure 4 shows part of such a tiling.

2.3 Tilings by pentacles

Figure 5 shows the pentacles. As with the rhombs, all of the sides have length one and the angles are all integer multiples of 5 . The labels indicate a matching condition, which again assures aperiodicity. The edges labeled 0 must t against edges labeled 0, 1 against 1, and 2 against 2. Note that the three pentagons are congruent, but have di erent matching conditions. A portion of a tiling by pentacles is shown in gure 6.

2

3 Coloring the tiles A tiling is called three-colorable if we may assign one of three distinct colors to each tile such that adjacent tiles have di erent colors. Tiles are said to be adjacent if their intersection is a line segment. Figures 7, 8, and 9 show three-colored tilings by kites and darts, rhombs, and pentacles respectively. Sibley and Wagon [5] proved that tilings by rhombs are three-colorable and Babilon [1] proved that tilings by kites and darts are three-colorable. The equivalent question for the pentacles seems to be open. Our pictures are the nal stage of a three-state, stochastic cellular automaton that can run on any tiling. The cellular automaton works as follows. First, assign one of three possible colors to each tile randomly. Then, allow the cellular automaton to evolve according to the following set of rules:  If the value of a cell (or tile) equals the value of a bordering cell that is closer to the origin (as measured by some arbitrary point chosen within each tile), then with 90% probability, the cell changes value randomly to one of the other two colors.  If the value of a cell does not equal the value of a bordering cell that is closer to the origin, but does equal the value of a cell farther away from the origin, then with 10% probability, the cell changes value.  If the value of the cell does not equal the value of any bordering cell, the cell does not change value. Note that three-colorings are stable under these rules. The hope is that threecolorings are attractive equilibria. Figure 10 demonstrates the algorithm on a small piece of a kite and dart tiling. The origin is located at the lower left corner of the triangle. The probabilities are chosen so that tiles close to the origin are generally colored correctly before tiles farther away from the origin. The values of 90% and 10% were chosen experimentally. Larger dynamic images are available on the author's web page: http://www.unca.edu/~mcmcclur/mathematicaGraphics/PenroseColoring/

4 Related coloring schemes Clearly the basic idea of this paper is applicable in other situations. A change in the number of states yields a class of algorithms for n-coloring 3

planar maps or graphs in general. For example, taking the number of states to be two, we can use the algorithm to two-color a checker board. While twocoloring a checker board is very simple, this gives us a rudimentary way of measuring the eciency of the algorithm. Experiments indicate that a 8  8 checkerboard is two-colored in about 32 generations on average. A 16  16 checkerboard is two-colored in about 82 generations on average. This seems ecient considering that an n  n checkerboard has 2n2 distinct colorings by two colors, only two of which are two-colorings. The algorithm may also be used to four color maps. A map of the United States is four colored in about 114 generations on average. While the algorithm is broadly applicable, it does not seem quite as good as more speci cally designed algorithms. The recursive algorithm based on Kempe chains described in chapter 24 of [6] four-colors most maps a little faster, for example. Furthermore, all of the maps described here so far have adjacency graphs all of whose vertices have small degree. More complicated maps (like Martin Gardner's April fools map described in [6]) can stump the algorithm completely.

5 Implementation notes All the images for this paper were generated with Mathematica. The tilings were generated using the DigraphFractals Mathematica package by the author as described in [4] These images were then converted to PlanarMap and PlanarGraph objects as de ned in the GraphColoring Mathematica package by Stan Wagon [6]. Code to run the cellular automaton on the PlanarGraph objects was written by the author. Final three-colored images were rendered by the ShowMap function de ned in the GraphColoring package. All code, dynamic images, and more examples are available on the author's web page: http://www.unca.edu/~mcmcclur/mathematicaGraphics/PenroseColoring/

References [1] Babilon, R., 3-colourability of Penrose kite-and-dart tilings. Discrete Mathematics, 2000, 235, 137-143. 4

[2] Gardner, M., Penrose Tiles to Trapdoor Ciphers, W. H. Freeman, New York, 1989. [3] Grunbaum, B. and Shepard, G. C., Tilings and Patterns, W. H. Freeman, New York, 1987. [4] McClure, M., Digraph self-similar sets and aperiodic tilings. Submitted to The Mathematical Intelligencer. [5] Sibley, T. and Wagon, S., Rhombic Penrose tilings can be 3-colored. The American Mathematical Monthly, 2000, 107, 251-253. [6] Wagon, S., Mathematica in Action, 2 ed., Springer-Verlag, New York, 1999.

5

The Kite

The Dart



1



1



1



1

Figure 1: The kite and dart

6

Figure 2: Part of a tiling by kites and darts

7

The Skinny Rhomb

The Fat Rhomb

Figure 3: The fat and skinny rhombs

8

Figure 4: Part of a tiling by rhombs

9

1

1

1

1

1

1

1

1 1

1

1

2

2

1

0

0 0 0

0

1

1

0

1

0

1

2

1 1

1

1 1

1

0

0 2

0 Figure 5: The pentacles 10

Figure 6: Part of a tiling by pentacles

11

Figure 7: Part of a three-colored tiling by kites and darts 12

Figure 8: Part of a three-colored tiling by rhombs

13

Figure 9: Part of a three-colored tiling by pentacles

14

Figure 10: Evolution to a three-coloring

15

A Stochastic Cellular Automaton for Three-Coloring ... - Mark McClure

Aug 24, 2001 - of the tiling maps each tile to another tile. Figure 2 shows part of such a tiling. ... dynamic images are available on the author's web page:.

178KB Sizes 3 Downloads 239 Views

Recommend Documents

Generating Google™ maps - Mark McClure
complete information is presented at the Google Maps™ API reference [2]. ... and display a Google map: an HTML file for the webpage and a javascript file that ...

Generating Google™ maps - Mark McClure
If you view the HTML file in a web browser, you should see the ..... We can illustrate the algorithm using the first 360 points of the Mt. Mitchell path. dataXML ...

Boundary scanning and complex dynamics - Mark McClure
with this notebook. .... [2] J. Glynn and T. Gray, The Beginner's Guide to Mathematica Version 4. Cambridge University Press, NY, 2000. BoundaryScanPP.nb. 9.

the prevalent dimension of graphs - Mark McClure
The extension of the various notions of \almost every" in Rn to infinite dimen- sional spaces is an interesting and difficult problem. Perhaps the simplest and most successful generalization has been through the use of category. Banach's application

Vibration of the Koch drum - Mark McClure
A preprint version of a “Mathematical graphics” column from .... More precisely, there are seven functions f0 , f1 , …, f6 that map the snow- flake onto the ...

development and performance of cellular automaton ...
critical load of a network if queueing costs are taken ... three cost functions on the critical load and through- ... for the development of our CA model of PSNs. The.

The connected locus for complex cubic iteration - Mark McClure
SupplementaryFiles directory, which should come with this notebook. The initialization .... 9 » a »2 +2 M, then the orbit of z0 will diverge to ¶. (See [4], page 266.) ...

The connected locus for complex cubic iteration - Mark McClure
SupplementaryFiles directory, which should come with this notebook. The initialization .... 9 » a »2 +2 M, then the orbit of z0 will diverge to ¶. (See [4], page 266.) ...

Boundary scanning and complex dynamics - Mark McClure
with this notebook. .... [2] J. Glynn and T. Gray, The Beginner's Guide to Mathematica Version 4. Cambridge University Press, NY, 2000. BoundaryScanPP.nb. 9.

the prevalent dimension of graphs - Mark McClure
An easy but important property of is that it respects closure. That is. (E) = (E). Another ( F] p. 41) is that the limsup need only be taken along any sequence fcng1n=1 where c 2 (01) and we still obtain the same value. One problem with is that it is

Vibration of the Koch drum - Mark McClure
We begin by setting up the boundary of the snowflake. The level .... Norm@interiorGrid@@#DD - KochVertices@@nDDD § stepSize к 2 &D@@1DD;.

The Read-Bajraktarevic Operator - Mark McClure
0.4. 0.6. 0.8. 1. References. [1] Massopust, Peter R. Fractal functions, fractal surfaces, and wavelets. Academic Press, Inc., San Diego, CA, 1994. ReadBajPP.nb.

Vibration of the Koch drum - Mark McClure
The fundamental modes of vibration of this drum can be modelled by the eigenfunctions of the .... We begin by setting up the boundary of the snowflake.

lattice gas cellular automaton modeling of surface ... - Semantic Scholar
now equal to Lpl, and we define a uniform random vari- able that determines where the adatom, if created, will land. Thus, we require only two random number gen- erations. Once the adatom lands on a particular lattice site r, we increase the height a

Multiagent Coordination by Stochastic Cellular ... - Semantic Scholar
work from engineering, computer science, and mathemat- ics. Examples ..... ing serves to smooth out differences between connected cells. However, if this ...

Generating self-affine tiles and their boundaries - Mark McClure
Now for each pair Ha, bL where a and b are chosen from , we want MHa, bL to denote the set of pairs of digits. Hd, d'L so that b = A a ..... used by the program.

"Decremental tag systems and random trees". - Mark McClure
1. Introduction. We fix a positive natural number m and consider sequences of the form xn ... tree construction we call the use it or lose it construction. In fact, our ...

"Decremental tag systems and random trees". - Mark McClure
We fix a positive natural number m and consider sequences of the form xn. = (x1,...,xn), where each ..... Sequences, http://www.research.att.com/~njas/sequences/.

Generating self-affine tiles and their boundaries - Mark McClure
For example, the image in figure 1 is a self-affine four-tile (i.e. it consists of four parts) ... Self-similarity and iterated function systems are, by now, fairly well known concepts. .... In figure 5, we see the image of figure 4 under the mapping

Parametric L-Systems and borderline fractals - Mark McClure
parameter and their application to borderline fractals. Note: To reduce the size of the ..... Maeder, R., The Mathematica programmer (Chapter 8). AP Professional ...

Multiagent Coordination by Stochastic Cellular Automata
A common thread in all multi- agent systems is the issue of .... r t fe¡ , p§uni becomes unstable and the other equilibria, p§#$!%!%!%!$ p§Q , become stable.

Parallel Automaton
parallel/synchronization automata table. ... changing conditions, by just changing the automata table. 2. ..... condition is reached (according to a reference-table.

The Borel Structure of the Collections of Sub-Self ... - Mark McClure
|T(x) -T(y)| = r|x-y| Vx y GRd: If r < 1, then T is called contractive. A fundamental result ( Ed], Thm. 4.1.3) states that if Ti : Rd → Rd is a contractive similarity for ...