Cyclotomic Aperiodic Substitution Tilings

The class of Cyclotomic Aperiodic Substitution Tilings (CAST) is introduced. Its vertices are supported on the 2n-th cyclotomic field. It covers a wide range of known aperiodic substitution tilings of the plane with finite rotations. Substitution matrices and minimal inflation multipliers of CASTs are discussed as well as practical use cases to identify specimen with individual dihedral symmetry Dn or D2n, i.e. the tiling contains an infinite number of patches of any size with dihedral symmetry Dn or D2n only by iteration of substitution rules on a single tile.


Introduction
Tilings have been subject of wide research. Many of their properties are investigated and discussed in view of their application in physics and chemistry, in detail the research of crystals and quasicrystals, such as D. Shechtman et al.'s renowned Al-Mn-alloy with icosahedral point group symmetry [SBGC84]. Tilings are also of mathmatical interest on their own [GS87,BG13]. Without any doubt they have great aesthetic qualities as well. Some of the most impressing examples are the tilings in M. C. Escher's art works [EBL82], H. Voderberg's spiral tiling [Vod36,Vod37] and the pentagonal tilings of A. Dürer, J. Kepler and R. Penrose [GS87,Lü00], just to name a few. In contrast to a scientist, a designer may have different requirements. Nevertheless the result may have interesting mathematical properties. Due to the lack of a general criteria for subjective matters of taste, we consider the following properties as preferable: • The tiling shall be aperiodic and repetitive (locally indistinguishable) to have an interesting (psychedelic) appearance. • The tiling shall have a small inflation multiplier for reasons of economy. Large inflation multipliers either require large areas to be covered or many tiles of a small size to be used. • The tiling shall yield "individual dihedral symmetry" D n or D 2n with n ≥ 4. I.e.
it shall contain an infinite number of patches of any size with dihedral symmetry only by iteration of substitution rules on a single tile. Similar to G. Maloney we demand symmetry of individual tilings and not only symmetry of tiling spaces [Mal14,Mal15].
The most common methods to generate aperiodic tilings are: 1 • Matching rules, as introduced in the very first publication of an aperiodic tiling, the Wang-tiling by R. Berger [Ber66]. See [Soc89,Soc90,Gä93] for more details. • Cut-and-project scheme, first described by de Bruijn for the Penrose tiling in [dB81], later extended to a general method, see [Lag96] and [Moo97] and references therein for more details, notably the earlier general method described in [Mey72]. • Duals of multi grids as introduced by de Bruijn [dB86] with equidistant spacings are equivalent to the cut-and-project scheme [GR86]. Duals of multi grids with aperiodic spacings as introduced by Ingalls [Ing92,Ing93] are close related to Ammann bars, see [GS87,Soc89,Lü93] and [Sch98,Ch. 5] for details and examples. • The idea of substitution rules or substitutions in general is a rather old concept, e.g. Koch's snowflake [vK04,vK06,Man77] or Rep-Tilings [Gar63]. However, it seems its first consequent application to tile the whole Euclidean plane aperiodically appeared with the Penrose tiling [Pen74, Gar77, Pen79, GS87, BG13].
Among those methods substitution rules may be the easiest approach to construct aperiodic tilings. Additionally they have some other advantages: • The inflation multipliers of tilings obtained by the cut-and-project scheme are limited to PV-numbers. According to [Lag96] and [Moo97] this was first noted by [Mey72]. This limitation does not apply to substitution tilings. • Matching rules tend to be more complex than substitution rules. See [GS98] and [GS03] for examples.
In view of the preferable properties we will introduce the class of Cyclotomic Aperiodic Substitution Tilings (CAST). Its vertices are supported on the 2n-th cyclotomic field. It covers a wide range of known aperiodic substitution tilings of the plane with finite rotations. Its properties, in detail substitution matrices, minimal inflation multiplier and aperiodicity are discussed in Section 2. CASTs with minimal or at least small inflation multiplier are presented in Sections 3 and 4, which also includes a generalization of the Lançon-Billard tiling. Section 5 focuses on several cases of rhombic CASTs and their minimal inflation multiplier. In Section 6 the "Gaps to Prototiles" algorithm is introduced, which allows to identify large numbers of new CASTs. Finally, examples of Girih CASTs with n ∈ {4, 5, 7} are presented in Section 7. Except the generalized Lançon-Billard tiling all CASTs in this article yield local dihedral symmetry D n or D 2n .
For terms and definitions we stay close to [BG13] and [FGH]: • A "tile" in R d is defined as a nonempty compact subset of R d which is the closure of its interior. • A "tiling" in R d is a countable set of tiles, which is a covering as well as a packing of R d . The union of all tiles is R d . The intersection of the interior of two different tiles is empty. • A "patch" is a finite subset of a tiling. • A tiling is called "aperiodic" if no translation maps the tiling to itself.
• "Prototiles" serve as building blocks for a tiling.
• Within this article the term "substitution" means, that a tile is expanded with a linear map -the "inflation multiplier" -and dissected into copies of prototiles in original size -the "substitution rule". • A "supertile" is the result of one or more substitutions, applied to a single tile.
Within this article we use the term for one substitutions only. • We use ζ k n to denote the n-th roots of unity so that ζ k n = e 2ikπ n and its complex conjugate ζ k n = e − 2ikπ n . • Q (ζ n ) denotes the n-th cyclotomic field. Please note that Q (ζ n ) = Q (ζ 2n ) for odd n. • The maximal real subfield of Q (ζ n ) is Q ζ n + ζ n . • Z [ζ n ] denotes the the ring of algebraic integers in Q (ζ n ).
• Z ζ n + ζ n denotes the the ring of algebraic integers (which are real numbers) in Q ζ n + ζ n .

Properties of Cyclotomic Aperiodic Substitution Tilings
We define Cyclotomic Aperiodic Substitution Tilings (CASTs) similar to the concept of Cyclotomic Model Sets as described in [BG13,Ch. 7.3].
Definition 2.1. A (substitution) tiling T in the complex plane is cyclotomic if the coordinates of all vertices are algebraic integers in Z [ζ 2n ], i.e. an integer sum of the 2n-th roots of unity. As a result all vertices of all substituted prototiles and the inflation multiplier are algebraic integers as well.
Theorem 2.1. For the case that all prototiles P k of a CAST T have areas equal to A k = A(P k ) = R · sin kπ n , R ∈ R >0 , k, n ∈ N , 0 < k < n we can use a given inflation multiplier η to calculate the substitution matrix if |η 2 | can be written as and the conditions |η 2 | / ∈ N and min((max(c k ), odd k), (max(c k ), even k)) ≥ 1, (even n > 4) are met. (For simplification, all prototiles with the same area are combined.) The inflation multiplier η of such a substitution tiling can be written as a sum of 2n-th roots of unity: Please note that there are multiple ways to describe η. We assume a k to be chosen so that the sum is irreducible, i.e. n−1 k=0 |a k | is minimal. A substitution tiling with l ≥ 2, l ∈ N prototiles and substitution rules is partially characterized by its substitution matrix M ∈ N l×l 0 with an eigenvalue λ and an eigenvector x A .
The elements of the right eigenvector x A contain the areas of the prototiles A k = A(P k ).
Since M ∈ N l×l 0 , the elements of x A generate a ring of algebraic integers which are real numbers.
The elements of the left eigenvector x f represent the frequencies of the prototiles f k = f (P k ), so that: The eigenvalue λ can be interpreted as inflation multiplier regarding the areas during a substitution. In the complex plane we can conclude: With (2.1) and (2.4) the eigenvalue λ can also be written as a sum of 2n-th roots of unity. In other words, the elements of the eigenvector x A span a ring of algebraic integers which are real numbers. This ring is isomorphic to the ring of algebraic integers Z ζ 2n + ζ 2n in Q ζ 2n + ζ 2n , which is the maximal real subfield of the cyclotomic field Q (ζ 2n ). (2.5) Because of the conditions regarding a k we ensure that no combination of roots of unity in the right part of the Equation (2.5) can sum up to a integer which is a real number. Since we need at least two roots of unity to have an inflation multiplier η > 1, it follows that The length of the k-th diagonal µ n,k of a regular n-gon with side length 1 can be written as absolute value of a sum of 2n-th roots of unity. Note that µ n,1 = 1 because it refers to a single side of the regular n-gon.
As a result the diagonals of a n-gon span a ring of diagonals Z [µ n ]. With Equations (2.8), (2.9) and (2.10) can be shown that: Because of Equations (2.11), (2.12), (2.13) and (2.14), we choose the substitution matrix as M ∈ N (n−1)×(n−1) 0 and the eigenvector as: With Equation (2.14) and (2.19), we can find a matrix M n,k , n/2 ≥ k ≥ 1 with eigenvalue µ n,k for the given eigenvector x A : (2.20) 1 0 · · · · · · · · · · · · 0 0 1 0 . . . . . . 0 1 0 0 · · · · · · · · · · · · 0 1 To get the substitution matrix for a given eigenvalue λ as defined in Equation (2.11) we just need to sum up the matrices M n,k with the coefficients c k : We recall that a substitution matrix M n must be primitive and real positive. To ensure the latter property its eigenvalue λ ∈ R must be a positive sum of elements of the eigenvector x A . For this reason, we have to modify Equation (2.11) accordingly: With Equation (2.14) it can be shown that the λ with positive coefficients c k in Equation (2.23) span a commutative semiring N 0 [µ n ] which is a subset of Z [µ n ]: Since all M n,k are symmetric, this is also true for M n . In this case, left and right eigenvector of M n are equal so that x A = x f and x A represents the frequencies of prototiles as well as their areas.
Because of Equations (2.15) and (2.16), the matrices and the eigenvector can be reduced so that M n ; M n,k ∈ N n/2 × n/2 0 and x A ∈ R n/2 . In detail we omit the redundant elements of x A and rows of M n and replace the orphaned entries in M n with valid ones. This also enforces n/2 ≥ k ≥ 1. The results of the reduction are shown in Equations (2.26) and (2.28).
For x A we can write: For odd n the substitution matrices M n can be described by the following scheme: The scheme can be continued for any odd n. The matrices are still symmetric, so that and x A still represents the frequencies of prototiles.
For even n, the substitution matrices M n be described by the following scheme: The scheme can be continued for any even n, but the matrices are not symmetric anymore. With Equations (2.3), (2.20) and (2.28), we can derive: Remark 2.1. For a given matrix M n of a CAST T as defined above the eigenvalue λ can be calculated easily with Equation (2.23). From Equations (2.26) and (2.28) follows that: In other words, the coefficients c k can be read directly from the bottom line of matrix M n .
We also can find a vector notation for Equation (2.23): For n ∈ {2, 3} the M n = c 1 M n,1 ∈ N; λ ∈ N and x A ; x F ∈ R >0 . As a result M n is always primitive.
For the following discussion of cases with n ≥ 4 we define some sets of eigenvalues (and its powers) based on Equation (2.23): (2.33) Equations (2.26) and (2.28) show that substitution matrix M n is strictly positive (M n 0) and so also primitive if the corresponding eigenvalue λ ∈ S f ull . Substitution matrix M n is also primitive if some finite power of eigenvalue λ p ∈ S f ull (p ∈ N).
For odd n ≥ 5 Equations (2.14) and (2.15) show that some finite power of diagonals of a regular n-gon exist so that: It is also possible to show that some finite power of an eigenvalue λ ∈ N 0 [µ n ] \ N 0 exist so that: As a result any substitution matrix M n in Equation (2.26) with an eigenvalue λ ∈ N 0 [µ n ] \ N 0 is primitive. For even n ≥ 4 Equations (2.14) and (2.15) show that the powers of diagonals of a regular n-gon have the following porperties: (2.41) µ p n,k ∈ S odd (odd k) (2.42) µ p n,k ∈ S odd (even k; even p) (2.43) µ p n,k ∈ S even (even k; odd p) With the same equations it is also possible to show that: It is also possible to show that some finite power of an eigenvalue λ ∈ S mixed \ N 0 exist so that: As a result any substitution matrix M n in Equation (2.28) with an eigenvalue λ ∈ S mixed \N 0 is primitive. That follows that the coefficients c k have to be chosen so that: (2.51) min ((max (c k ) , odd k) , (max (c k ) , even k)) ≥ 1 (even n ≥ 4) Remark 2.2. Every substitution tiling defines a substitution matrix. But not for every substitution matrix M n exists a substitution tiling.
Theorem 2.2. CASTs T in Theorem 2.1 with n ≥ 4 are aperiodic under the condition that at least one substituion rule contains a pair of copies of a prototile whose rotational orientations differs in kπ n ; n > k > 0; k ∈ N.
Proof. We recall that any substitution matrix M n of CASTs T in Theorem 2.1 is primitive and that every CAST T yields finite rotations.
Under this conditions [Fre08, Theorem 2.3] applies, so that all orientations of a prototile appear in the same frequency.
We furthermore recall areas and frequencies of prototiles of CASTs T in Theorem 2.1 are given by right and left eigenvector x A and x f .
With Equations (2.27) and (2.29) the ratio between the frequencies of prototiles P 2 and P 1 is given by: Since f (P 2 ) f (P 1 ) / ∈ Q and all orientations of a prototile appear in the same frequency, no periodic configuration of prototiles P 1 and P 2 is possible.
As a result the CASTs in Theorem 2.1 with n ≥ 4 are aperiodic.
Remark 2.3. The inflation multiplier makes no statement whether a CAST meets the condition in Theorem 2.2 or not. As a consequence possible solutions have to be checked case by case. However, all examples within this article fulfill the condition in Theorem 2.2.
Remark 2.4. All prototiles of a CASTs T in Theorem 2.1 with n ∈ {2, 3} have identical areas. For these CASTs aperiodicity has to be proven with other methods.
With the following inequalities we can identify all eigenvalues which are smaller than λ min . As noted in Equation (2.23) all c k ≥ 0, therefore every eigenvalue λ must be a sum of diagonals µ n,k .
Theorem 2.4. The definition of CASTs with n ≥ 4 can be extended to all CASTs with other prototiles than in Theorem 2.1.
Proof. Let a CAST T * exist with l ≥ l min prototiles P * k , l ≥ k ≥ 1 with a real positive, primitive substitution matrix M * n ∈ N l×l 0 and a real positive eigenvector x * A ∈ R l whose elements represent the relative areas A(P * k ) A(P * 1 ) of the prototiles: . . .
Since T * is a CAST by definition and all coordinates of all vertices are algebraic integers in Z [ζ 2n ] the areas of all pototiles are real numbers so that As a consequence a transformation matrix T ∈ Z l× n/2 and its inverse T −1 ∈ Q n/2 ×l exist such that: In other words, if there is a CAST T * then also a corresponding CAST T (as defined in Theorem 2.1) with the same eigenvalue λ does exist. That includes but may not be limited to all cases where the CASTs T and T * are mutually locally derivable.
A similar approach is possible for the frequencies of the corresponding CASTs T and T * . As a result Theorem (2.2) also applies to CAST T * .
The corresponding CASTs T and T * have the same eigenvalue, so Theorem (2.3) also applies to CAST T * .
Remark 2.5. l min is given by the algebraic degree of λ which depends on n. In other words, λ is a solution of an irreducible polynomial with integer coefficients, and λ is of at least l min -th degree. l min can be described by Euler's totient function denoted by Φ, in detail l min = Φ(n) 2 , odd n and l min = Φ (n), even n. Remark 2.6. The eigenvalue λ as noted in Equation (2.23) and (2.31) is unique if l min = n/2 . This is the case if n is a prime. If n is not a prime, λ may have more than one corresponding vector c. E.g. for n = 9 we can show that µ 9,4 = µ 9,2 + µ 9,1 . An eigenvalue λ = µ 9,4 − µ 9,1 = µ 9,2 implies that the corresponding vector c ∈ Only the latter one is real positive. As a consequence we will consider λ ∈ N 0 [µ n ] as noted in Equation (2.24) if at least one real positive c exist.

CASTs with Minimal Inflation Multiplier
In the this section, we discuss CASTs with minimal inflation multiplier as noted in Theorem 2.3.
3.1. The case odd n. The substitution matrix M n,min for CASTs with minimal eigenvalue λ min = µ n,3 + 1, odd n is given by the following scheme: The case n = 5 describes the Penrose tiling with rhombs or Robinson triangles and individual dihedral symmetry D 5 . The cases with n > 5 are more difficult and require additional prototiles. For n = 7 we give an example of a CAST T * with individual dihedral symmetry D 7 as shown in Fig. 3.1 and the following properties. A(P k ) stands for the areas of a prototile P k .
A further example with n = 9, minimal inflation multiplier but without individual dihedral symmetry D 9 has been found but is not included here.
Another CAST with n = 7, individual dihedral symmetry D 7 and the same minimal inflation multiplier is given in [ND96, Fig. 1, Sec. 3, 2nd matrix] and shown in Fig. 3.2.
The status of Conjecture 3.1 is subject to further research.
Formally all preferable conditions are met. However, for large n, this type of CAST tends to be complex and the density of patches with individual dihedral symmetry D n tends to be small. 3.2. The case even n. The substitution matrices M n,min for CASTs with minimal eigenvalue λ A,min = µ n,2 + 2 is given by the following scheme:   The Lançon-Billard tiling (also known as binary tiling) [LB88, GL92, BG13] is a rhomb tiling with inflation multiplier |η| = |1 + ζ 1 2n | for the case n = 5. It is possible to generalize it for all n ≥ 4. The approach is very similar to the approach in [Har05]. The generalized substitution rules are shown in Fig. 3.3.
The eigenvector x B and substitution matrix M * n are given by: . . .
Remark 3.2. For odd n, the scheme and the substitution matrix M * n can be separated into two independent parts. We choose the part which relays on prototiles P k , n > k ≥ 1, odd k only. The eigenvector is given by the areas of the rhombic prototiles with side length 1 and area A(P k ) = sin kπ n . With Equation (2.15) and (2.16) we can write: . . .
As a result, we can use Equation (2.26) to describe the substitution matrix: (3.9) M n = M n,2 + 2E (odd n) Regarding generalized Lançon-Billard tilings, we can confirm the results sketched in [Hib15].
The generalized Lançon-Billard tiling does not contain any patches with individual dihedral symmetry D n or D 2n . However, for n = 4 and n = 5, other CASTs with individual symmetry D 8 and D 10 have been derived and submitted to [FGH]. An example for n = 4 is mentioned in Section 6 and shown in Fig. 6.3.

CASTs with Inflation Multiplier Equal to the Longest Diagonal of a Regular Odd n-Gon
Another interesting approach to identify CASTs with preferred properties is to choose a CAST as described in Theorem 2.1 and a relative small inflation multiplier, in detail the longest diagonal of a regular n-gon with odd n: (4.1) η = µ n, n/2 (odd n) Because of Equations 2.4, 2.14 and 2.26, we can write: Conjecture 4.1. CASTs as described in Theorem 2.1 with inflation multiplier η = µ n, n/2 , substitution matrix M n = M n, n/2 2 and exactly n/2 prototiles and n/2 corresponding substitution rules exist for every odd n ≥ 5.
The substitution rules have no dihedral symmetry, so different combinations of the prototiles chirality within the substitution rules define different CASTs with identical inflation multiplier. For the cases 7 ≤ n ≤ 11, solutions with individual dihedral symmetry D n have been found by trial and error. For the case n = 7, see  Formally all preferable conditions are met. However, for CASTs with large n as described in Conjecture 4.2 the density of patches with individual dihedral symmetry D n tend to be small.

Rhombic CASTs with Symmetric Edges and Substitution Rules
For many cases of CASTs with given n and rhombic prototiles, the minimal inflation multiplier can be described if some additional preconditions are met. In the following section, we will focus on cases where all edges of all substitution rules are congruent and the inner angles of the rhomb prototiles are integer multiples of π n . Definition 5.1. "Edge" means here a segment of the boundary of a supertile including the set of tiles which are crossed by it.  Despite these restrictions we still have several options left: • There are two ways to place rhombs on the edge of substitution rules. We recall that the inner angles of the rhombs are integer multiples of π n . We can place all rhombs on the edge so that the inner angles either with even or odd multiples of π n are bisected by the boundary of the supertile. We will call these two cases "even" and "odd edge configuration", for details see Fig. 5.1. A "mixed" configuration is not allowed, because it would force the existence of rhombs with inner angle equal to k + 1 2 π n .
• We can choose the symmetry of the substitution rules and their edges. Possible choices are dihedral symmetry D 1 and D 2 . Edges with dihedral symmetry D 1 can have the boundary of the supertile or its perpendicular bisector as line of symmetry. The smallest nontrivial solution for the latter case is the generalized Goodman-Strauss tiling [Har05]. Since this example does not provide individual dihedral symmetry D n or D 2n in general, we will focus on the other case. Substitution rules of rhombs which appear on the edge of a substitution rule are forced to have the appropriate dihedral symmetry D 1 as well. This is also true for substitution rules of prototiles which lie on the diagonal i.e. a line of symmetry of a substitution rule. The orientations of the edges have to be considered as well. These three conditions may force the introduction of additional rhomb prototiles and substitution rules. Additionally, the existence of edges with orientations may require additional preconditions.
To avoid this problem, a general dihedral symmetry D 2 can be chosen for the substitution rules and their edges. • Parity of the chosen n may require different approaches in some cases, similar to the example of the generalized Lançon-Billard tiling in Section 3.2 and Fig. 3.3.
The options above can be combined without restrictions, so in total we have eight cases:  Within this section, we will denote a rhomb with side length d 0 ≡ 1 and an inner angle kπ n , n > k > 0 which is bisected by the boundary of a supertile as R k and its diagonal which lies on the boundary of a supertile as d k . The line segments which appear on the edges of substitution rules in the cases with even edge configuration are denoted by R 0 and their length by d 0 . Under those conditions, we can write the inflation multiplier η as sum of the d k : N 0 , cases 2 and 4) The inflation multiplier η can also be written as a sum of 2n-th roots of unity, where the roots occur in pairs ζ i 2n + ζ i 2n so that η ∈ R: cases 1 and 3) (5.4) We borrow a remark regarding worms from [FH13]: "By definition, a parallelogram has two pairs of parallel edges. In a parallelogram tiling therefore there are natural lines of tiles linked each sharing a parallel edge with the next. These lines are called worms and were used by Conway in studying the Penrose tilings [GS87]. We will follow Conway to call these lines worms." Every rhomb R k with an inner angle kπ n , n > k > 0 which is bisected by the boundary of a supertile work as entry-exit-node of two worms. We also can say that such a rhomb on the edge "reflects" a worm back into the substitution rule. The line segments R 0 work as entry-exit-node of one worm only.
The Kannan-Soroker-Kenyon (KSK) criterion as defined in [KS92,Ken93] decides whether a polygon can be tiled by parallelograms. We use a simplified phrasing for the criterion: All line segments R 0 and all inner line segments of rhombs R k , n > k > 0 on the edge of the substitution rule serve as nodes. All corresponding nodes of the substitution rules are connected by lines. The KSK criterion is fulfilled if for every intersection of lines the inner angles between the corresponding nodes are larger than 0 and smaller than π.
We start with observations of the substitution rule for rhomb prototile R 1 or R n−1 . In detail we will focus on those both edges which enclose the inner angle equal π n and which we will denote as "corresponding edges". We will denote the area where the both edges meet as the "tip" of the substitution rule.
(1) The rhombs on the edges must not overlap. For this reason, for the tip of the substitution rule, only three configurations are possible as shown in Fig. 5.2. Obviously, a tip as shown in Fig. 5.2 (c) is compliant to the cases 2 and 4 with odd edge configuration and Fig. 5.2 (b) to cases 1 and 3 with even edge configuration. Fig. 5.2 (a) requires the even edge configuration as well. Since all edges are congruent, it must be the start and the end of the same edge, which meet on that vertex. Since start and end of the edge are different, it can not have dihedral symmetry D 2 . For this reason, the tip in Fig. 5.2 (a) is not compliant to case 3. (2) Any rhomb R k with n − 3 ≥ k ≥ 3 on one edge implies the existence of a rhomb R k−2 on the corresponding edge. In turn, rhomb R k on one edge is implied by a rhomb R k+2 on the corresponding edge or a rhomb R k on the opposite edge. An example is shown in Fig. 5.3. (3) Any rhomb R 2 on one edge implies the existence of a line segment R 0 on the corresponding edge. In turn, rhomb R 2 on one edge is implied by a rhomb R 4 on the corresponding edge or a rhomb R 2 on the opposite edge.
(4) Any rhomb R 1 on one edge implies the existence of a rhomb R 1 on the corresponding edge. In turn, rhomb R 1 on one edge is implied by a rhomb R 3 on the corresponding edge or a rhomb R 1 on the opposite edge. (5) Any line segment R 0 on one edge is implied by a rhomb R 2 on the corresponding edge or a line segment R 0 on the opposite edge. (6) Any rhomb R n−2 on one edge implies the existence of a rhomb R n−4 on the corresponding edge. In turn, rhomb R n−2 on one edge is implied by a rhomb R n−2 on the opposite edge. (Rhomb R n does not exist, because the inner angle would be zero.) (7) Any rhomb R n−1 on one edge implies the existence of a rhomb R n−3 on the corresponding edge. In turn, rhomb R n−1 on one edge is implied by a rhomb R n−1 on the opposite edge. (Rhomb R n+1 does not exist, it would have an inner angle greater than π or smaller the 0.) (8) If for a rhomb R k , n > k > 0 on one edge two related elements (rhomb or line segment) R |k−2| and R k+2 exist on the corresponding edge, R |k−2| is closer to the tip than R k+2 Since all edges are congruent we can derive the following inequalities for α k : (5.5) α 0 ≥ α 2 ≥ α 4 ≥ ...α 2 (n−1)/2 ≥ 0 (α k ∈ N 0 , cases 1 and 3)
For case 1b and 3b: • Any line segment R 0 on the edge implies the existence of a rhomb R n−1 on the correspondent edge or a line segment R 0 on the opposite edge. As shown in Fig. 5.4 (a) and (b), the existence of rhomb R n−1 on the edge is not required to meet the KSK criterion. • Any rhomb R 2 on the edge implies the existence of a rhomb R n−3 on the correspondent edge or a rhomb R 2 on the opposite edge. As shown in Fig. 5.4 (c) and (d), the KSK criterion is only met if at least one R n−3 exists on the edge.
(5.9) α n−3 > α n−1 ≥ 0 (cases 1b and 3b) For case 1a and 3a: • The line segment R 0 on the edge implies the existence of a line segment R 0 on the opposite edge only (rhomb R n does not exist). • Any rhomb R 2 on the edge implies the existence of a rhomb R n−2 on the correspondent edge or a rhomb R 2 on the opposite edge. So the KSK criterion is only met if at least one R n−2 exists on the edge.
(5.10) α n−2 ≥ 1 (cases 1a and 3a) For case 2b and 4b: • Any rhomb R 1 on the edge implies the existence of a rhomb R n−2 on the correspondent edge or a rhomb R 1 on the opposite edge. So the KSK criterion is only met if at least one R n−2 exists on the edge.
(5.11) α n−2 ≥ 1 (cases 2b and 4b) For case 2a and 4a: • Any rhomb R 1 on the edge implies the existence of a rhomb R n−1 on the correspondent edge or a rhomb R 1 on the opposite edge. So the KSK criterion is only met if at least one R n−1 exists on the edge.
(5.12) α n−1 ≥ 1 (cases 2a and 4a) In the next part we will determine the minimal inflation multiplier η min of a rhombic CAST. We will discuss only case 1 in detail to show the general concept.
We continue our observations of the substitution rule R 1 or R n−1 . There are four possible combinations of the corresponding edge orientations as shown in Fig. 5.5. Figure 5.6. Substitution rule for rhomb R 1 or R n−1 : KSK criterion for corresponding edges with different orientations with α 0 = α 2 and α 2k = α 2k+2 ; k > 1.
We recall Equations (5.5) and (5.9). To get a small inflation multiplier, the α k with even k have to be chosen as small as possible. We use the KSK criterion to check whether substitution rules can exist for different orientations of the corresponding edges and equal amounts of rhombs R k and R k+2 so that α k = α k+2 . The checks are illustrated in Fig. 5.6. Actually, the figure just shows the case α k = α k+2 = 1, while the result is true for α k = α k+2 > 1 as well.
The results for α 0 = α 2 in Fig. 5.6 (a), (b), (c) show that the KSK criterion is fulfilled for more than one combination.
The results for α k = α k+2 ; k > 0 in Fig. 5.6 (d), (e), (f) show that the KSK criterion is fulfilled for one orientation only.
At this point we have to introduce an additional condition. For cases 1 and 2, we consider only rhombic CASTs whose prototiles R 1 or R n−1 yield at least two different orientations of the corresponding edges.
Remark 5.1. The author is only aware of one example of a rhombic CAST whose prototiles R 1 or R n−1 show only one orientation of the corresponding edges, namely the Ammann-Beenker tiling with n = 4. For n > 4, the existence of such rhombic CASTs seem to be unlikely. However, a proof to rule out this possibility is not available yet.
With the results and conditions above we can conclude: With Equations (2.9), (5.3), (5.5) and (5.7), we can describe the minimal inflation multiplier: For case 1b with odd n, we can use the same concept of orientations as shown in [ND96, Fig. 2] which leads to rhombs with orientations as shown in Fig. 5.7. As a result, the orientations of the edges are globally defined. In detail, the edges which enclose an inner angle π n have orientations as shown in Fig. 5.5 (c) and (d) and a tip as shown in Fig. 5.2 (a). For case 1a with even n this simplification can not be used. The edges which enclose an inner angle π n have not only orientations as shown in Fig. 5.5 (c) and (d). Depending on how exactly the orientations are defined one of the orientations in Fig. 5.5 (a) or (b) may appear. Similar approaches are possible for cases 2, 3 and 4. For cases 3 and 4 the symmetry conditions of the substitution rules and its edges enforce α k even, k > 1, just α 0 is not required to be even.
With the minimal inflation multiplier η min as shown in Table 2 -5, we can start the search of rhombic CASTs by trial and error. The minimal rhomb edge sequences therein were also obtained in this way, so additional solutions with varied sequences obtained by permutation may exist. Examples for all eight cases are shown in Fig. 5.8 -5.11.
Conjecture 5.1. For all cases as listed in Table 1, there exist rhombic CASTs for any n with edge configuration and inflation multiplier η min as shown in Tables 2 -5 which yield individual dihedral symmetry D n or D 2n .
A general proof of Conjecture 5.1 is subject to further research. However, a proof for case 4b already exists. The results in [KR15,KR16] for odd n match very well with our results in case 4b. The case even n in the same publication is almost identical to our case 3a, up to a difference of one line segment R 0 in the center of the edge of the substitution rules. Remark 5.3. Socolar's 7-fold tiling [FGH, credited to J. Socolar] can be derived from case 1b with n = 7 with an modified minimal rhomb edge sequence.
Remark 5.4. Two substitution steps of the Penrose rhomb tiling are equivalent to one substitution step in case 1b with n = 5. Table 2. Rhombic CAST substitution rule edge configuration for case 1.

Gaps to Prototiles Algorithm
In this section, we sketch a "Gaps to Prototiles" algorithm to identify CASTs for a given n and a selected edge of a substitution rule (as in Definition 5.1). Conditions: • All prototiles have inner angles equal kπ n . • All edges of all substitution rules are congruent and have dihedral symmetry D 2 .
• As discussed in Section 5 the tiles on the edge have to be placed, so that the inner angles either with even or odd multiples of π n are bisected by the boundary of the supertile.
• The tiles on the edge are bisected by one or two lines of symmetry of the edge. This implies dihedral symmetry D 1 or D 2 of the corresponding substitution rules. • The inflation multiplier η must fulfill the conditions in Theorem 2.1.
• The inflation multiplier η is defined by the sequence of tiles which are part of the edge. Algorithm: (1) We start with the prototiles which appear on the edge of the substitution rule.
(2) We start the construction of the substitution rules by placing the prototiles on the edge. (3) If the edge prototiles overlap the algorithm has failed. In this case, we may adjust the sequence of rhombs or other equilateral polygons on the edge and start another attempt. (4) We try to "fill up" the substitution rules with existing prototiles under consideration of the appropriate dihedral symmetry D 1 or D 2 . If gaps remain, they are defined as new prototiles and we go back to step (2). Please note, if a gap lies on one or two lines of symmetry, the substitution rule of the new prototile must also have the appropriate dihedral symmetry D 1 or D 2 . (5) If no gaps remain the algorithm was successful.
In some cases the "Gaps to Prototiles" algorithm delivers results with preferable properties as shown in Fig. 6.1, 6.2 and 6.3. However, some results contain a very large number of prototiles with different sizes and complex shapes as in Fig. 6.2. It is not known yet whether the algorithm always delivers solutions with a finite number of prototiles. Table 6. Inflation multipliers of CASTs identified by the "Gaps to Prototiles" algorithm n inflation multiplier 7 2µ 7,3 + 2µ 7,2 + 1 Fig. 6.1  7 µ 7,2 + 2

Extended Girih CASTs
"Girih" is the Persian word for "knot" and stands for complex interlaced strap works of lines, which are a typical feature of Islamic architecture and design. A common definition is given in [All04]: "Geometric (often star-and-polygon) designs composed upon or generated from arrays of points from which construction lines radiate and at which they intersect." The oldest known examples of star pattern date back to the 8th century AD [Lee87]. Girih designs are known in many styles and symmetries, see [Bou73] for examples. A variant of Girih design relies on Girih tiles and tilings as shown in the reproduction of the Topkapi Scroll in [NAA95]. The decorations on the tiles consist of lines which run from tile to tile when they are joined together. So the borders between joined tiles seem to disappear.
According to [NAA95] the shapes of all girih tiles are equilateral polygons with the same side length and inner angles kπ 5 , k ∈ {2, 3, 4, 6} including: Because of these properties, Girih tilings are cyclotomic tilings as well.
A Girih (cyclotomic) aperiodic substitution tiling was derived from a mosaic at the Darb-i Imam Shrine, Isfahan, Iran which dates back from 1453. It relies on the regular decagon and two hexagons and has individual dihedral symmetry D 10 . It was published in [LS07a,LS07b]. However, the complete set of substitution rules can be found in [Ten08,Cro08]. Additional examples of Girih CASTs have been discovered and submitted to [FGH]. More examples have been published in [Lü90,LL94,Cro08].
That rises the question if Girih CASTs with other symmetries and relative small inflation multiplier exist. For this reasons we have to define the extended Girih CASTs. The following definition turned out to be useful for a given n ≥ 4.
• The substitution rule of the regular 2n-gon has dihedral symmetry D 2n for odd n and D n for even n.
The Girih CASTs in Fig. 7.1, 7.2 and 7.4 -7.6 have been obtained by a trial and error method under the following conditions: • In every corner of every substitution rule a regular 2n-gon is placed.
• Edge and inflation multiplier have been derived from a periodic pattern of regular 2n-gons and their inter space counterparts.
(For n ≥ 8 this approach might fail. As an alternative it is possible to reuse inflation multipliers from rhombic CASTs in Section 5, case 1 and 2. See Table 2 and 3 for details.) The Girih CAST in Fig. 7.3 has been derived from [Cro08, Fig. 14, 15]. In detail the nonconvex hexagons were replaced and the star shaped gap in the center of the substitution rule of the pentagon prototile [Cro08, Fig. 15 (b)] was eliminated. This was possible by changing the symmetry of the substitution rule from dihedral symmetry D 5 to D 1 . The tilings in [Cro08,Fig. 14,15] were derived by an analysis of patterns shown in the Topkapi Scroll, in detail [NAA95, Panels 28, 31, 32, 34].
Remark 7.1. Please note that the case even n requires special care to make sure that the regular 2n-gons with dihedral symmetry D n always match. It seems that the existence of a substitution rule of the regular 2n-gons with dihedral symmetry D 2n requires the existence of prototiles with inner angle π n which are forbidden due to our preconditions. As a result, additional substitution rules might be necessary for prototiles with the same shape but different orientations.
Remark 7.2. The decorations at the prototiles are related but not necessarily equivalent to Ammann bars.

Summary and Outlook
Although the motivation behind this article was mainly aesthetic, some significant results have been achieved. Cyclotomic Aperiodic Substitution Tilings (CASTs) cover a large number of new and well known aperiodic substitution tilings as shown in Table 8. The properties of CASTs, in detail their substitution matrices and their minimal inflation multipliers, can be used as practical starting point to identify previously unknown solutions. For many cases, such solutions yield individual dihedral symmetry D n or D 2n .
The different approaches to identify CASTs have their individual advantages and disadvantages. The preferable properties as listed in the introduction section may be complemented by a high frequency of patches with dihedral symmetry. A promising approach we do not discuss in this article may be the choice of inflation multipliers which are PV numbers.
The results in this paper focus strictly on the Euclidean plane so that Equation (2.4) applies. However, the methods described herein might be adapted for other cases as well.
Finally several conjectures have been made, which require further research.

Acknowledgment
The author dedicates this paper to his daughter Lili and his parents Marita and Herbert. He would like to thank M. Baake, D. Frettlöh, U. Grimm, R. Lück and C. Mayr for their support and encouragement.