Randomness and Topological Invariants in Pentagonal Tiling Spaces

We analyze substitution tiling spaces with ﬁvefold symmetry. In the substitution process, the introduction of randomness can be done by means of two methods which may be combined: composition of inﬂation rules for a given prototile set and tile rearrangements. The conﬁgurational entropy of the random substitution process is computed in the case of prototile subdivision followed by tile rearrangement. When aperiodic tilings are studied from the point of view of dynamical systems, rather than treating a single one, a collection of them is considered. Tiling spaces are deﬁned for deterministic substitutions, which can be seen as the set of tilings that locally look like translates of a given tiling. ˇCech cohomology groups are the simplest topological invariants of such spaces. The cohomologies of two deterministic pentagonal tiling spaces are studied.


Introduction
Aperiodic tilings of the plane appeared in the literature in the works of Wang 1 and Penrose 2 .Substitution tilings with noncrystallographic planar symmetries have been intensively studied in the last decades, mainly since the discovery of quasicrystals.The structures are meaningful in several areas like the study of quasicrystalline materials and artificially fabricated macroscopic structures that can be used as photonic or phononic devices.
A property of quasicrystal structures is the appearance of sharp peaks in their diffraction patterns and recent results in this direction use methods familiar from statistical mechanics and from the long-range aperiodic order of tilings 3 .A suitable approach is to work with translation invariant families of mathematical quasicrystals, instead of dealing with a single one.The action of the dynamical systems, which is usually on time, is now translation on space.For recent advances in the mathematics of diffraction in the context of dynamical systems and stochastic spatial point processes, see 4 .On the other hand, the atoms in a material modeled with a quasicrystal tiling are distributed in such a way that they determine a quasiperiodic potential.The spectrum of the associated Schr ödinger Hamiltonian has infinitely many gaps, and its distribution is related to the integer Čech cohomology of the corresponding space of tilings 5, 6 .Other questions tied to the Čech cohomology are connected to the derivation of the internal structure of a material from diffraction data or the type of deformations of the molecular structure that are consistent with the combinatorics of the molecular bonds 7, 8 .
Pentagonal, octagonal, decagonal, and dodecagonal quasicrystalline materials have been found in experiments, and tilings with the corresponding symmetries in the diffraction patterns are candidates of their structural models.In this paper, we first study the derivation of both deterministic and random tilings and then we discuss the cohomology groups of some tilings generated by an inflation process.In Section 2, we review a general geometric construction for substitutions and we show how it is applied to the case d 5.
Then we analyze the introduction of randomness in the tilings generation and we compute the configurational entropy for certain cases.To conclude, we consider tiling spaces and we study the cohomology groups of the deterministic pentagonal patterns in Section 4.

The Substitution Tilings
In 9 , the authors studied the systems of tangents of the deltoid and derived a method for the construction of substitution tilings for d-fold symmetries, with d being odd and not divisible by 3.They are particular cases of subsystems of lines included in the constructions described below and in 10, 11 .By a system of lines, we mean a set of d straight lines d-system appearing in d orientations.In what follows, E n denotes the integer part of n and s ν ≡ sin νπ/d .The following sets of lines in the xy plane define for each d, up to mirror images, finite patterns made with triangles having edge lengths s ν .The patterns contain not only the prototiles but also, when arrows are added as explained below, the first inflation step level-1 supertiles and therefore the necessary information for the derivation of the substitution or inflation rules.For d 2m, m 3, 4, . . . the systems of lines are formed by x 0, y 0, and The inflation factors for simple tilings are ν d s ν /s 1 , ν 2, 3, . . .d, and composite tilings, where the inflation factors are products of the form ν d μ d • • • λ d , can be obtained also; although in order to include all the inflation factors, additional constructions must be used 11 .In Section 3, we will study the case d 5 which may be generated with the constructions considered in this section.
In general, arrows must be added to the triangle edges in order to generate tilings by means of substitution rules.The arrows are not necessary for the d 2m tilings with inflation factors s m /s 1 because the edge substitution rules are palindrome 10, 11 .A 2dsystem contains a d-system inheriting a structure on each edge.The d-system prototile edges appear subdivided into two segments.We adopt the criterion of assigning an arrow in the direction going from the longest to the shortest segment, and no arrow when the edge has a barycentric subdivision.When the tiles are juxtaposed along an edge, the arrows on the edge match.
The system of lines for d 10 are Figure 1 x 0, y 0,

2.2
This 10-system contains the subsystem thick lines in Figure 1 x 0, The lines of the 10-system intersect the 5-system on points that determine the arrow directions Figure 1 b .The 2D tilings can be described in terms of word sequences in D0L-systems as in 12 .A 0L system is a triple G {Σ, h, ω} where Σ {x 1 , x 2 , . . ., x n } is an alphabet, h is a finite substitution on Σ into the set of subsets of Σ * , and ω ∈ Σ * is the axiom or starting symbol.G is called a D0L system if # h x i 1, for every x i ∈ Σ, that is to say, there is only one possible substitution for each tile.Now we consider the tilings generated with the 5-system.The triangular tile T m x, y, z has edges x, y, z placed anticlockwise, and the index m ∈ Z 10 denotes relative orientation.The letters A m , B m represent the prototiles T m α, α, β , T m α, β, β , respectively, where α, β have lengths s 2 , s 1 , with s k ≡ sin kπ/5 .The tiles of type A then have edges with relative lengths 1, 1, θ, and the tiles of type B have sides θ, θ, 1, where θ s 2 /s 1 is the golden number.The alphabet is Σ {A m , A m , B m , B m , , } with m ∈ Z 10 .It contains two brackets and letters of type X i and X i which represent reflected tiles.The set of production rules h for the tiling T − is Figure 2 a :

2.4
Any letter representing an oriented tile can be used as axiom.By iterating the production rules applied to an axiom, we get word sequences that describe the tiling growth.We take the axiom A m for a given m.In the word B m 4 A m 6 A m 3 , if two letters follow one another inside a bracket, the corresponding oriented triangles are glued face to face in a unique way.In the next derivation step, which gives Φ − B m 4 Φ − A m 6 Φ − A m 3 , two oriented triangles represented by consecutive words enclosed by brackets like Φ − B m 4 and Φ − A m 6 are glued face to face and again the prescription is unique.
A different type of tilings T can be generated with substitution rules characterized by Φ X m Φ − X m for all the prototiles represented by X m ∈ Σ with X m X m .Penrose tilings 2, 13 can be obtained also with the same system of lines but the arrowing is changed for some edges Figure 2 b .Having in mind the different arrow decorations for the prototiles, the substitution rules for the Penrose pattern can be obtained formally from Φ − if we make the replacements A m ↔ A m .Fragments of T − , T and Penrose tilings can be seen in Figures 3  and 4.

Random Substitutions
Random tiling problems are of interest from a mathematical and also physical point of view.Quasicrystal order may appear even for random tilings, and it is not yet known if deterministic tilings are better candidates than random tilings for the description of materials with noncrystallographic symmetries.A method for the derivation of random substitutions with arbitrarily high symmetries has been introduced recently 11 .We analyze two types of nondeterministic tilings, one is related to composition of inflation rules and the other to tile rearrangements.

Composition of Inflation Rules
The composition of inflation rules gives tilings that are edge to edge.This is a property that is not in general fulfilled by other random tilings studied in the literature  14, 15 .Random substitutions can be described in terms of stochastic L-systems.A stochastic 0L-system is a four-tuple G π {Σ, h, ω, π}.The alphabet, the set of productions, and the axiom are defined as in a 0L-system.The function π : h → 0, 1 , called the probability distribution, maps the set of productions into the set of production probabilities.We define a stochastic DT0L-system as a four-tuple G π {Σ, H, ω, π} where H is a set of homomorphisms and π : H → 0, 1 .The map is defined for every element in H in contrast to the previous definition where it is defined for every production rule in h.For the case, we are considering H {Φ − , Φ } and a fragment of one of the patterns can be seen in Figure 5. New vertex configurations appear which are not present in T − , T or Penrose tilings.

Tile Rearrangements in the Inflation Rules
Another way to get nondeterministic tilings is to introduce tile rearrangements in the inflation rules.First we derive a tiling with the deterministic inflation rules corresponding to an m-system followed by prototile subdivision by the d-system which contains it d 2m, 4m, . . ., and then we apply tile rearrangements.In Figure 6 a we can see the 10-system with possible edge flips forming edges of a pentagon.The tiles corresponding to the 5system are represented in Figure 6 b .We begin by applying the substitutions Φ − or Φ , and then we get two possibilities for subdivision of one of the tiles and just one for the other Figure 6 c .The configurational entropy S for a sequence of patterns obtained by a random substitution process is defined as the logarithm of the number of patterns of a given size and shape or level-n supertiles divided by the number of tiles N n in the thermodynamic limit The frequencies of the tiles F A , F B in the tiling are given by the elements of the normalized eigenvector corresponding to the eigenvalue with largest modulus, or Perron-Frobenius eigenvalue, of the 2D prototile substitution matrix which in this case is The Perron-Frobenius eigenvalue is θ 2 with algebraic conjugate θ −2 , and the characteristic polynomial is x 2 − 3x 1.The frequencies are then

3.3
The general solution to the difference equation The number of A and B tiles after n iterations is , respectively, and the total number of tiles is The number of patterns after iterating n times level-n supertiles is where ω Φ X is the number of tile rearrangements in the first inflation step of X A, B. By taking into account 3.9 The case d 5 included in d 20 can be seen in Figure 7.We have now ω Φ A 5, ω Φ B 3 Figure 8 and S θ Log 2 ≈ 1.12.

3.10
If we randomize T − or T by considering the 5-system included in a d-system with higher d, we can obtain random tilings with increasing entropy values.Also the combination of the two methods discussed in this section is possible, and they give always edge to edge tilings.
The relation between diffuse scattering and randomness has been studied in 4, 16 .By comparing diffraction of point sets ranging from deterministic to fully stochastic, the authors show that diffraction is in some cases insensitive to the degree of order.The following is an open question: what is the role of entropy in the diffraction patterns of the structures described above.

Cohomology Groups of the Deterministic Tiling Spaces
A tiling space Ω T can be seen as the set of tilings that locally look like translates of a tiling T.
Anderson and Putnam studied the cohomology of substitution tiling spaces as inverse limits of branched manifolds 17 .They proved that the cohomology can be computed by means of a CW complex Γ 1 on collared tiles which are formed to get a substitution that "forces its border", a concept introduced in 6 .The cell complex Γ 0 contains one copy of every kind of tile that is allowed with some edges identified for the 2D cases, and the result is a branched surface.If somewhere in the tiling, a tile shares an edge with another tile, then those two edges are identified.Tiles labelled by the pattern of their nearest neighbors are called collared tiles, and the cell complex Γ 1 is obtained by stitching one copy of each collared tile.If σ is the map representing the substitution rule, we denote by γ k k 0, 1 the map induced on the cell complex Γ k by σ.A substitution is said to force the border if there is a positive integer n such that any two level-n supertiles of the same type have the same pattern of neighboring tiles.A method for describing an arbitrary substitution tiling space by a substitution that forces the border was introduced in 17 .It is obtained by rewriting the substitution in terms of collared tiles.If H k Γ 1 , Z denotes the Čech cohomology with integer coefficients of the complex Γ 1 , then it is shown in 17 that the cohomology Ȟk Ω σ of the tiling space associated with σ is isomorphic to the direct limit of the system of abelian groups for k 0, 1, 2, . ... If the substitution forces its border, then the same conclusions hold replacing Γ 1 and γ 1 by Γ 0 and γ 0 .

Penrose Tiling Space
Now we study the cohomology of Penrose tiling spaces Ω P along the lines of 8, 18 but with different inflation rules.The rotation group Z 10 acts freely on edges and tiles.We have two edges α, γ with length s 2 , 2 edges β, δ with length s 1 , and four tile types A T r 5 δ, r 8 γ, r 7 α , B T r 5 γ, r 9 δ, rβ , A T δ, r 8 α, r 7 γ , and B T γ, r 4 β, r 6 δ with r 10 1 and each appears in 10 orientations.The vertices of A and B can be seen in Figure 2  The four irreducible representations of Z 10 over the integers are the 1-dimensional scalar r 1 and pseudoscalar r −1 representations and two 4-dimensional representations.The vector and the pseudovector representations have r acting by multiplication on the rings R 1 Z r / r 4 − r 3 r 2 − r 1 and R 2 Z r / r 4 r 3 r 2 r 1 , respectively.In this case the vertices appear in the scalar and pseudoscalar representations, while the edges and faces appear in all representations.If the cochain groups are denoted by C k , then, for k 0, 1, the coboundary maps δ k ∂ T k 1 : C k → C k 1 are obtained from the boundary maps: in all representations and in the scalar and pseudoscalar representations.

Pseudoscalar Representation r −1
In the pseudoscalar representation, C 0 Z 2 , C 1 Z 4 , and C 2 Z 4 also.But now both δ 0 and δ 1 have rank 2 and

Vector Representation
There are no vertices in the vector representation, while C 1 , C 2 are free modules of dimensions 4 over the ring R 1 .The matrix δ 1 has rank 3 over R 1 and, as abelian groups, we have

Pseudovector Representation
In the pseudovector representation, C 1 , C 2 are free modules of dimensions 4 over the ring R 2 .The map δ 1 is an isomorphism and H 0 H 1 H 2 0. We have obtained the cohomology of the uncollared complex Γ 0 which is enough because the Penrose tiling forces the border 17 , in four steps Figure 9 .To get the cohomology of the tiling space Ω P , we need to compute the direct limit of the cohomologies under the substitution.But the substitution is invertible and the direct limit of each H k is simply H k .By taking into account all the irreducible representations, we get the well-known result:

The Tiling Space Ω Ξ
The tiling space Ω Ξ is associated with Φ described in Section 2. Iteration of the inflation rules shows that the tiles appear in 5 different orientations.We modify the inflation rules in such a way that, in any inflation step, the prototiles appear in the same 5 orientations Figure 10 .The possible vertex configurations can be seen in Figure 11.The analysis of level-6 supertiles shows that the substitution forces the border.In fact, the level-6 superedges with relative length s 2 have the following pattern of vertex configurations Figure 12 5327329253273292532925, 4.5 Therefore, it is enough to consider, up to rotation, 2 edges α, β with lengths s 2 , s 1 , respectively, and 4 tiles:

4.7
The rotation group Z 5 acts freely on edges and tiles.We have four tile types and two edges, each in five orientations.All the vertices are identified to one vertex a with a ra.A different case with the same property is the Ammann-Beenker octagonal tiling 19 .The Anderson-Putnam complex Γ 0 of Ω Ξ has Euler characteristic χ 11.
The two irreducible representations of Z 5 over the integers are the 1-dimensional scalar r 1 representation and the 4-dimensional vector representation which has r acting by multiplication on the ring Z r / r 4 r 3 r 2 r 1 .The vertex appears only in the scalar representation, while the edges and faces appear in all representations.For k 0, 1, the coboundary maps δ k ∂ T k 1 : C k → C k 1 are deduced from the boundary maps: and ∂ 1 0.

Scalar Representation r 1
In the scalar representation, the cochain groups are C 0 Z, C 1 Z 2 , and C 2 Z 4 .The rank of δ 1 is 2, and then we have H 0 Ker δ 0 Z, H 1 Ker δ 1 /Im δ 0 0, and Now we look at the direct limit of the cohomologies.The substitution acts as the identity on the vertex a.The substitution matrix on 2-chains is, according to 4.7 , 4.9 while the matrix on 1-chains is The induced matrix on 2-cochains is σ * 2 σ T 2 which is an isomorphism, and the direct limit gives a contribution of Z 2 to Ȟ2 .

Vector Representation
In the vector representation, C 0 is trivial while C 1 , C 2 are free modules of dimensions 2 and 4 over the ring R Z r / r 4 r 3 r 2 r 1 .The matrix δ 1 has rank 1 over R. In this case, H 0 0,H 1 has one copy of the ring R, and H 2 has three copies of R; therefore, as abelian groups, we have H 1 Z 4 , H 2 Z 12 .The ranks of σ * 1 , σ * 2 over R are 2 and 4, respectively, and the direct limit gives a contribution of Z 4 to Ȟ1 and Z 12 to Ȟ2 .
Adding up the contributions of each representation, we get

4.11
The group Ȟ1 has one copy of the vector representation for both Ω Ξ and Ω P and one of the scalar representation only for Ω P .We also see that Ȟ2 contains two copies of the scalar representation for both Ω Ξ and Ω P , and the additional terms are two copies of the pseudoscalar representation for Ω P , while the vector representation gives a contribution of Z 4 for Ω P and Z 12 for Ω Ξ .The tiling spaces Ω T that we consider in this paper are formed by the closure of the translational orbit of one tiling.The finite rotation group Z n acts on Ω T , and the quotient Ω T /Z n yields the space Ω T of tilings modulo rotation about the origin 8 .The cohomology of Ω T is the rotationally invariant part of the cohomology of Ω T , and we have for T P, Ξ , but Ȟ1 Ω P Z whereas Ȟ1 Ω Ξ 0.

The Tiling Space Ω Ξ −
In this case, connected with Φ − in Section 2, the substitution does not force the border and we have to use collared tiles.We have now eight vertices 1, 2, 3, 4, 5, 5, 6, 6.In Figure 13, where we have represented the uncollared vertex configurations, the vertices without bilateral symmetry correspond to 5, 6, 7, 11.Disregarding the vertices in the borders, there are, in the collared tiles, eight edge types: α, β, γ, γ with length s 2 and δ, η, , with length s 1 associated with the following vertex sequences Figure 14 in the level-6 superedges commas are used when a vertex is represented by two digits :  D , α, β , E γ , , δ , F β, , δ , and G γ, , , respectively.Observe that with this notation A, B, C, D have Section 2A-shape and E, F, G the B-shape.Now by taking into account the 8 vertices, we have, up to rotation, 86 edges and 106 prototiles.The edges and tiles appear in five different orientations and the complex Γ 1 of Ω Ξ − has Euler characteristic χ 108.The inflation rules for the prototiles, up to mirror reflection, can be seen in Table 1, where each letter representing a tile type is followed by their vertices and by the prototiles that appear in the substitution rules.In Figure 15, it is shown a fragment of a pattern obtained by assigning colors to the collared prototiles.We have treated this case without the use of representations, along the lines of 17 see also 20 for recent results applicable when the tiling does not force the border .The matrices are to large to list here but may be read off from Table 1 The substitution matrix on vertices is 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0

Concluding Remarks
Several types of deterministic and random substitutions have been analyzed.We have shown that, apart from the vertex configurations and relative orientations, an essential difference between the spaces of pentagonal tilings Ω Ξ , Ω Ξ − and the very well-known space of Penrose tilings Ω P is the cohomology, in spite of the fact that the first inflation step in the tiling growth seems to be very similar.Both tiling spaces have finitely generated torsion-free cohomology groups see 21 for cartesian product tiling spaces related to the constructions in Section 1 .
For projection tilings, it has been shown in 22 that in the first cohomology group there is at least a subgroup isomorphic to the reciprocal lattice of the tiling.The minimal dimension of the lattice for an N-fold symmetry tiling in 2D obtained by projection is given by the Eulers totient function φ N , which counts the number of positive integers less than N that are coprime to N. For the Penrose tiling φ 10 4 and its description by projection can be done with the root-lattice A 4 23 .As far as the author knows, Ω Ξ , Ω Ξ − have the simplest first cohomology group and hence K 1 -group of a tiling space in 2D with noncrystallographic symmetries, and its rank coincides with φ 5 .The complexity of the tiling spaces, if measured with the set of possibly collared vertex configurations, is reflected in the rank of the second cohomology group and as a consequence in the K 0 -group.For certain cases 8 , only the rotationally invariant part of the second cohomology group contributes to the gap labeling group, and therefore the Penrose and T tilings have isomorphic groups.However, in contrast to the Penrose tiling, T has a vanishing first cohomology group on the space of tilings modulo a rotation.Another tiling with the same prototiles was derived by the T übingen group 23 .For a study of its cohomology, which is computationally demanding, see 24 .
We have obtained the configurational entropy of sequences of finite patterns generated by particular random substitution processes.In general, the determination of the entropy of the tiling spaces associated to the random tiling ensembles defined possibly as direct sums of uniquely ergodic spaces and its physical meaning are open questions.
The constructions studied in Section 2 have been used for the generation of nonperiodic ordered tilings in multiconnected flat manifolds.For instance, the d-systems corresponding to d 8, 12 were used in 10 to get random triangulations of the torus and the Klein bottle.In Figure 6 a , it can be seen how a pentagon appears in the modification of the 10-system, the decagon is contained in the 20-system modification Figure 6 and in general we have an m-gone included in the 2m-system.This fact can be used to get recursive random simplicial structures for all the closed 2-manifolds, by using two polygons with certain identifications.By considering also higher dimensional manifolds this type of constructions could be of interest in different fields.
Some open problems that should also be considered in the future are the relationship with projection methods, matching rules 25 , and the analysis of the topological invariants for both deterministic and random substitution tiling spaces with N-fold symmetries generated with the constructions given in Section 2 and in 11 .

Figure 1 :
Figure 1: System of lines for d 5, 10 a and induced arrows for d 5 b .

Figure 2 :
Figure 2: Arrowed prototiles and inflation rules of a T − and b Penrose tilings.

Figure 3 :
Figure 3: Fragments of a T − , b T and c Penrose tilings.

Figure 6 :
Figure 6: a The 10-system with thick lines indicating flipped edges.b The first inflation step corresponding to d 5 is indicated by thick lines; prototile subdivision by the 10-system is represented with thin lines and the thicker line indicates the possible tile flip.c First inflation step followed by prototile subdivision and the two possibilities of subdividing one of the prototiles.
2 are constants determined by the initial conditions about the prototiles content.

Figure 7 :Figure 8 :
Figure 7: The 20-system with thick lines indicating flipped edges and the d 5 system.
b , while the vertices of A and B can be obtained by replacing {a, b, c} → {b, d, a} and {c, d, b} → {d, a, c} on A and B. They satisfy a rb, b ra, c rd, and d rc.The uncollared Anderson-Putnam complex Γ 0 has Euler characteristic χ 4.

Figure 13 :
Figure 13: Atlas of vertex configurations of T − .

Figure 15 :
Figure 15:Fragment of a pattern obtained by assigning colors to the collared tiles, after iteration of the inflation rules in Table1for Ω Ξ− .