ARS MATHEMATICA CONTEMPORANEA ISSN 1855-3966 (printed edn.), ISSN 1855-3974 (electronic edn.) ARS MATHEMATICA CONTEMPORANEA 16 (2019) 157-171 https://doi.org/10.26493/1855-3974.997.7ef (Also available at http://amc-journal.eu) Regular polygonal systems Jurij Kovic * Institute for Mathematics, Physics and Mechanics, Jadranska 19, 1000 Ljubljana and FAMNIT, University of Primorska, Glagoljaska 8, 6000 Koper, Slovenia Received 19 December 2015, accepted 4 March 2018, published online 18 November 2018 Let M = M(Q) be any triangle-free tiling of a planar polygonal region Q with regular polygons. We prove that its face vector f (M) = (f3, f4, f5,...), its symmetry group S(M) and the tiling M itself are uniquely determined by its boundary angles code ca(M) = ca(Q) = (t1,... ,tr), a cyclical sequence of numbers ti describing the shape of Q. Keywords: Regular polygonal system, boundary code, face vector, symmetry group, reconstructibility from the boundary. Math. Subj. Class.: 05B40, 05B45 1 Introduction Systems of regular (planar or spherical) polygons joined edge to edge arise in various contexts (in tilings, in polyhedral maps, in nature, in chemistry, in art). A rich theory of such systems may be developed. Researchers usually focus on some particular class of such systems (defined by some conditions), try to determine all its elements and explore various questions related to their combinatorial description, parameters, enumeration, characterization, classification, coding, etc. To unify the investigations of such objects and to emphasize their common characteristics we propose a general concept of a regular polygonal system and make some first few steps towards a general theory of such systems. Here we give definitions, examples and remarks; the general reconstruction problem is presented in Section 2; results are gathered in Section 3. A polygonal system M is a (finite or infinite) incidence structure M = (V, E, F) whose elements are called vertices, edges and faces: faces are abstract polygons - cyclical sequences of vertices (v1,..., vn), and edges are pairs of vertices [vi, vi+1}. Two faces are *This work is supported in part by the Slovenian Research Agency (research program P1-0294 and research project N1-0032). E-mail address: jurij.kovic@siol.net (Jurij Kovic) Abstract ©® This work is licensed under https://creativecommons.Org/licenses/by/4.0/ 158 Ars Math. Contemp. 16 (2019) 157-171 incident if they share an edge, two edges are incident if they share a vertex. If M consists only of n-gons it is called a n-system; if the number of these n-gons is m, it is a nm-system. If all these faces are congruent polygons X, it is a monohedral system denoted Xm. If M consists only of ni-gons, n2-gons, ..., nk-gons, it is called a (ni,..., nk)-system. The face vector (or just the f-vector) of the polygonal system M is the sequence f (M) = (f3, f4, f5,...) where f (i) = fi(M) are the numbers of its faces with i edges. We use also the notation f (M) = (3f (3), 4f(4), 5f (5),...). A planar or spherical polygon P is called regular if there is a cyclic group G = (R | Rn = I) of rotations acting transitively on the vertices and edges of P, and if its boundary dP is simply connected (thus we exclude star polygons, as in Kepler solids). A regular polygonal system M = M(Q) consists of regular polygons joined edge to edge, covering a polygonal planar or spherical region Q. The symmetry group of M = M(Q) is the group of the rotations and reflections of Euclidean space E2 or E3 preserving the incidences in M. A code c(X) of a given mathematical object X is any (not necessarily discrete) structure from which X can be completely reconstructed (up to isomorphism). A boundary code of a given m-dimensional object X is any code c(dX) of the shape of its (m - 1)-dimensional boundary dX. The boundary angles code of the planar or spherical regular polygonal system M = M(Q) covering (tiling) a polygonal region Q is the cyclical sequence Ca(M) = Ca(Q) = (ti,t2,.. .,tr) of the angles ti = 180° - a G (-180°, 180°), where ai are the interior angles between adjacent edges Ai-1Ai and AiAi+1 of Q. The boundary faces-edges code Cf,e(M) is the cyclical sequence of symbols f (i)e(i), where f (i) is the number of edges of i-th boundary face and e(i) is the number of the boundary edges of this face (see examples in Appendix A). A few examples and remarks will help us better understand these definitions and the motivation for them. In chemistry, the mathematical model of benzenoid molecules are called benzenoid systems or polyhexes (composed of regular hexagons). Similar systems are polydiamonds (composed of equilateral triangles) andpolyominoes (composed of squares) (see Figure 1). The motivation for the boundary angles code ca comes from the "turtle geometry" [1]. The definition of boundary faces-edges code cf e is motivated by the boundary codes of various n-systems presented below. Several boundary codes exist for the planar regular hexagonal systems B. One of them is the boundary edges code [6]. This code ce(B) = (k1,k2,... ,kr) is a cyclical sequence counting the numbers ki of boundary edges in boundary hexagons; we travel around the boundary in the clockwise direction, starting at any hexagon, and having the interior of B always at our right (see [2, 7]). The code for the regular planar triangular systems (or 3-systems) T may be defined exactly in the same way - as the boundary edges code ce(T) counting the number of boundary edges in boundary triangles of T. Thus, for the 3-system T in Figure 1 its code is ¿aEbnra Figure 1: Planar regular n3-systems for n = 3,4,5,6. ce(T) = (2,1, 2). J. Kovic: Regular polygonal systems 159 The case of the planar square systems (or 4-systems) S is trickier. To get the right numbers in ce(S) we must count also the "zeros" of boundary edges in boundary vertices adjacent to no boundary edge. We can use also a simple boundary vertices-faces code cVff (S), a cyclical sequence counting for each boundary vertex how many squares (1, 2 or 3) are incident with that vertex. For the first 4-system S in Figure 1 we have ce(S) = (2,3,0, 3) and cVff (S) = (1,2,1,1, 3,1,1,2). There is also a code for the regular pentagonal systems P (described in [3]) that actually counts the numbers of vertices incident only with boundary edges. For the second 5-system from Figure 1 this code is (3,3, 2). One can easily imagine polygonal systems composed of regular polygons with different numbers of edges, and also non-planar generalizations of these concepts. Example 1.1. The face vector of the planar regular (3,4, 5)-system M in Figure 2 is f (M) = (33,4i, 52) = T3S1P2, since there are three triangles T (f3 = 3), one square S (f4 = 1) and two pentagons P (f5 = 2). Its symmetry group S(M) is generated by one reflection (over the vertical axis). Figure 2: A polygonal system M with the symmetry group S(M) = Z2. Example 1.2. The boundary faces-edges code of the planar system M in Figure 3 is Cf,e(M) = (54,40, 32,42). The boundary angles codes ca(M) of the spherical systems with the same Cf,e(M) depend on the size of the spherical triangle T. Increasing the size of T the interior angles a increase as well, hence the angles tj = 180° - a decrease. If T is big enough, the angle between the pentagon and the triangle vanishes and we get a region with Cf,e(M) = (53, 3i, 42). Figure 3: A planar regular polygonal (31, 41, 51 )-system M with the boundary angles code c0(M) = (120°, 30°, 90°, -18°, 72°, 72°, 72°, -78°). The maps of Platonic solids are examples of regular 3-systems (tetrahedron 34, octahedron 38, icosahedron 320), 4-systems (cube 46) and 5-systems (dodecahedron 512). The 160 Ars Math. Contemp. 16 (2019) 157-171 polygonal system of the square pyramid is (B, C, D, E), (A, B, C), (A, C, D), (A, D, E), (A, E, B), thus it is a (3,4)-system. The system (A, B, C, D), (A, B, D, C) represents the Mobius band tiled with two quadrilaterals. In a polygonal system more than two faces may share the same edge, as in the "3-page book" (A, B, C, D), (A, B, D, E), (A, B, F, G), which is a 3-system, too. 2 The reconstruction problem The first motivation for this research was the realization that it is possible to generalize the method that worked so well for benzenoid systems B: "Encode the boundary of a system as a cyclical sequence of numbers, and then obtain various information about B from its boundary code b(B)" to other planar polycyclic molecules. Thus b(B) was used to find various relations between the parameters of B, to determine its symmetry scheme, to calculate its face count, etc. (as in [5, 7, 8]). With this motivation in mind we have in the Introduction already defined two boundary codes cej(M) and ca(M). Now we can ask: Is it possible to reconstruct the face vector f (M), the symmetry group S(M) or even the whole planar regular polygonal system M = M(Q) from its boundary code Cf,e(M) or ca(M)? Obviously, some regular polygonal systems M are reconstructible from the chosen boundary code of the region Q covered by M. How to characterize such systems? This question may be stated in a more general form: Given some class C of objects X whose boundary is determined by some boundary code c(dX) characterize those X from C that can be reconstructed from its boundary code, in other words, find X for which c(dX) is also the code c(X) of X. The code is not necessarily a cyclical sequence of numbers. The boundary of any simply connected polycube P (a solid composed of cubes of unit length joined face to face) may be coded with a graph G(P) whose vertices are boundary square faces of P and whose edges are labeled with angles (0, n/2, -n/2) between adjacent boundary faces. The number of cubes in P, the symmetry group S(P) and the polycube P itself are obviously determined by its boundary. However, their actual reconstruction from the graph G(P) is complicated. The code is not always a discrete structure. In analysis, a differentiable function f is reconstructible from its derivative f' by the Newton-Leibniz formula up to an additive constant C. Likewise, a harmonic function f: Q ^ R is determined by its values on the boundary of Q. Not all the codes of the same object contain the same amount of information. For example, two similar planar triangles have different codes (lengths a, b, c and a',b', c'), but the same boundary code (consisting of angles a, p, y). In some cases the chosen boundary code c(dX) contains some additional information about the structure of X that cannot be deduced from dX alone. Thus it may happen that it is possible to reconstruct X from the code of dX, although it is not possible to reconstruct X only from dX. Example 2.1. Let Q be a planar polygonal region composed of a regular 12-gon and a square sharing one side. There are two tilings M1 and M2 of Q with 12 squares S and 20 equilateral triangles T, having the same boundary angles code (90°, 90°, 60°, 30°, 30°, 0°, 30°, 30°, 30°, 0°, 30°, 30°, 30°, 0°, 30°, 30°, 30°, -90°), but different boundary faces-ed- J. Kovic: Regular polygonal systems 161 ges codes (see Figure 4): ce,f (Mi) = (43, 4o, 3i, 4i, 3i, 3o, 3i, 4i, 3i, 4i, 4i, 3i, 4i, 3i, 3o, 3i, 4i, 3i, 4i, 4o), Ce,f (M2) = (43, 3o, 4i, 3i, 4i, 4i, 3i, 4i, 3i, 3o, 3i, 4i, 3i, 4i, 4i, 3i, 4i, 3i, 3o, 3o). Hence, it is not possible to reconstruct M from the shape of the boundary dQ, encoded by ca(Q). However, M is reconstructible from Cf,e(M), since there are only two possible tilings Mi and M2 of Q. In this sense, Cf,e(M) is a stronger code (containing more information), and ca(M) = ca(Q) is weaker (containing less information). Figure 4: Two different planar (3,4)-systems M1 and M2 with the same boundary and with the same face vector (320,412). Remark 2.2. Two closed planar polygonal regions Q and Q' may have the same boundary angles code ca = (t^ ... , tr), but if the ratio of the lengths of the corresponding edges AjAi+i and AjAj+1 are not all the same, then the shapes of these regions are different. Note that the lengths of the edges of a region Q tiled by regular polygons joined edge to edge are the same. This is the reason why the boundary angles code ca(Q) suffices to describe the shape of the boundary of Q, tiled by regular polygons.1 3 Results In this Section we show: if M is a regular triangle-free tiling of a planar polygonal region Q, then we can use ca(M) to find the face vector f (M) (Theorems 3.1 and 3.8), to determine the symmetry group S(M) (Theorem 3.12), or even to completely reconstruct M = M(Q) (Theorem 3.5). Theorem 3.1. Let M = M(Q) be a planar regular (3,4)-system covering the polygonal region Q. Then the boundary angles code ca(M) = ca(Q) determines the face vector f (M) of M. Proof. The area of Q can be calculated from the boundary angles code ca(Q) = (t 1,..., tn) as follows: fix the coordinates of two adjacent vertices A1 = (0,0) and A2 = (1,0) of Q, use vectors to find the coordinates of other vertices Aj, triangulate Q and sum the areas 1 However, as we see in Figure 4, the tiling M (Q) may not be uniquely determined by the shape of its boundary, although it is tiled by regular polygons. 162 Ars Math. Contemp. 16 (2019) 157-171 of all these triangles. The area of the regular n-gon Pn with side 1 is Area(Pn) — n ■ cot(360°/n). These numbers are incommensurable at least for n — 3,4 since Area(P3) — a/3/4, Area(P4) — 1. Now it is easy to see that the integer solutions of the equation xi(%/3)/4 + y1 — x2(\/3)/4 + y2 are possible only if x1 — x2, y1 — y2. Thus f (M) is determined by ca(M). Solving the equation x1^v/3)/4 + y1 — Area(M) is easy, since the calculated expression for Area(M) must appear in this form, from which we just read x1 and y1. □ Theorem 3.2. The sum of the angles in the boundary angles code ca(Q) — (t1, t2, ... ,tr) of a planar polygonal region Q is t1 +t2 + ■ ■ ■ + tr — 360°. Proof. This follows from the formula J2™=i a = (n — 2) • 180° for the sum of the interior angles of a n-gon: J2 "=1 ti = 2 "=1 (180° — a). Another proof in the context of the "turtle geometry" is given in [1, p. 175]. □ Regular (planar or spherical) n-gons are usually denoted by the symbol {n} or just n. The vertex type of the interior vertex of the regular polygonal system M is defined as the cyclical sequence (a.b.c.... ) of the faces a, b, c,... surrounding it. The planar vertex types in this notation are listed in [4]. It is easy to check the following very useful observation ([4, p. 60]). Theorem 3.3. If the planar regular n1-gon, ..., nr -gon surround a vertex without gaps and overlaps, then 3 < r < 6 and (n1 — 2)/n1 + • • • + (nr — 2)/nr = 2, hence there are 21 types of vertices surrounded by regular polygons in a plane without gaps or overlaps: 3.3.3.3.3.3, 3.3.3.6, 3.3.3.4.4, 3.3.4.3.4, 3.3.4.12, 3.4.3.12, 3.3.6.6, 3.6.3.6, 3.4.4.6, 3.4.6.4, 3.7.42, 3.9.18, 3.8.24, 3.10.15, 3.12.12, 4.4.4.4, 4.5.20, 4.6.12, 4.8.8, 5.5.10, 6.6.6. Theorem 3.4. The possible (interior or boundary) vertex types in spherical regular polygonal systems are (if we exclude spherical 2-gons): • 5 triangles: 3.3.3.3.3; • 4 triangles: 3.3.3.3.4, 3.3.3.3.5; • 3 triangles: 3.3.3, 3.3.3.4, 3.3.3.5; • 2 triangles: 3.3, 3.3.m, 3.m.3 (m > 4), 3.3.4.5, 3.4.3.5, 3.3.5.5; • 1 triangle: 3, 3.m (m > 4), 3.4.n, 3.n.4 (n > 4), 3.5.n, 3.n.5 (n > 5); • 0 triangles: m, 4.m, 4.4.m, 4.m.4 (m > 4), 4.5.n, 4.n.5 (19 > n > 4), 5.5.5. Proof. We just use the fact that for the interior angle an of the spherical n-gon it holds that 180° > an > 180° — 360°/n, and check all the possible cases. A similar treatment of faces is given in Appendix B. □ Theorem 3.5. Let M = M(Q) be a planar regular polygonal system covering the polygonal region Q. If M contains no triangles then it is reconstructible from its boundary angles code ca(M). Likewise, M is reconstructible from ca(M) also in the case M is without squares and hexagons. J. Kovic: Regular polygonal systems 163 Proof. The theorem is certainly true if M contains only one face. Now suppose it is true for the systems with m or less faces and let M = M(Q) be a system with m + 1 faces. Since the sum of the r inner angles a of Q is J2¿=1 a = (r - 2) • 180°, at least one a must be smaller than 180°. The vertex A of such a < 180° is incident with at most three regular polygons: with two or three triangles, with one triangle and one square or pentagon or with two squares, or with a single n-gon with the inner angle (n - 2) • 180°/n. Hence, if M is triangle-free, then the angle a < 180° is incident with only one n-gon Pn whose interior angle is (n - 2) • 180°/n = Oj. Removing this Pn from M we get a smaller system, whose tiling is unique, by the induction hypothesis. Hence the position of each polygon in M is uniquely determined. Similarly, if there are no squares and hexagons in M, then every interior angle a < 180° is incident either with a single polygon Pn (and then we proceed as before) or with a triangle and a pentagon. But there is no planar interior vertex type containing faces 3 and 5 (see Theorem 3.3). Therefore both the vertices of the edge AC shared by a boundary triangle ABC and a pentagon are boundary vertices (see Figure 5). If we interchange the Figure 5: Left: A (3, 5)-system with two boundary components and with the cyclical symmetry group C15. Right: An illustration of the fact that such system cannot have two different tilings. positions of two adjacent boundary faces 3 and 5 along the adjacent boundary edges then the interchanged face 5 covers a neighbourhood of C (as in Figure 5 right). So there is only one possible tiling of M in the neighbourhood of all boundary points incident with 3.5 or 5.3. Removing this triangle and pentagon we get a smaller region for which the theorem is true by the induction hypothesis. □ Remark 3.6. It may happen that the boundary of M - Pn is no longer a simple closed curve (see Figure 6); it may have crossings and it may have more than one component. However, by repeating the process of removing the polygons corresponding to interior angles a < 180° from the system and calculating the boundary angles code of the smaller system we may find the exact locations of each face in the system algorithmically. For example, removing a boundary square P4 from M changes the boundary code from ca (M) = (tl,t2, . . . , ti-i,tj = «4 = 90°, tj+1, . . . ,tr) into Ca(M - P4) = (¿1, ¿2,^-1 + «4, «4 -180°, a4 + ti+1,..., tr), as in Figure 6. For n > 5 we get a more complicated formula for ca(M - Pn), dependent on how many successive angles tj, tj+1,... ti+k in ca(M) are equal to (n - 2) • 180°/n. Just as in the planar case, there are polygonal regions on the sphere admitting more 164 Ars Math. Contemp. 16 (2019) 157-171 Figure 6: Removing of a boundary face may produce a system that is no more face-connected. than one regular tiling, too. A spherical pentagon, covered by five regular triangles sharing a vertex of the spherical icosahedron is an example of such a region. Theorem 3.7. Let M be a spherical regular polygonal system M without triangles and squares. If the boundary of M has only one component, then M is reconstructible from its boundary faces-edges code Cf,e(M). Proof. By Theorem 3.4 M has no interior points. Hence there can be no cycle composed of adjacent faces in M, and Cf,e(M) completely describes the system M. □ Theorem 3.8. The face vector f (M) = (f3, f4, f5,...) and consequently also the number of faces f = f3 + f4 + f5 + ■ ■ ■ of any regular planar polygonal system M = M (Q) without hexagons or triangles and with all sides of length 1 is uniquely determined by the area of Q, and this area is uniquely determined by its boundary angles code ca(Q) = ca(M). Proof. If M = M(Q) has at most two faces, then its face vector is obviously determined by its boundary code, except in the case of a hexagon, which may be decomposed into six triangles. Suppose the theorem is true for any system with n faces. Take a system M with n+2 faces. In every boundary vertex with a positive angle tj we can repeat the procedure of taking out the polygons of M as in the proof of Theorem 3.5. In the boundary vertices with the interior angles filled with 3.4, 4.3, 3.5, 5.3, we take away from M both combinations and get two smaller systems M* and M** composed of n faces which must have the same area, hence they have the same face vector by the induction hypothesis. Therefore M has the unique face vector, too. □ Theorem 3.9. The symmetry group S(M) of a planar or spherical regular polygonal system M = M(Q) is a subgroup of S(Q) = S(dQ). Proof. Any rotation or reflection preserving M preserves the region Q, tiled by the polygons of M. The symmetry group of the boundary of Q is obviously isomorphic to the symmetry group of Q. □ A symmetry of the boundary of Q does not necessarily induce a symmetry of M = M(Q). For example, if we glue together the two 12-gons from Figure 4 (without the added top square) along the vertical edge, we get a region with two reflection symmetries (over the vertical axis and over the horizontal axis), yet its tiling has only one reflection symmetry (over the horizontal axis). J. Kovic: Regular polygonal systems 165 Theorem 3.10. If there is only one regular polygonal system M = M(Q) covering the region Q then the groups S(M) and S(Q) are isomorphic. Proof. If the tiling M of Q with regular polygons is unique, then every symmetry of Q automatically induces an automorphism of M. □ Lemma 3.11. Let Q be a planar polygonal system with given lengths of its edges li = AiAi+1 and with the boundary angles code c(Q) = (t\,... ,tr). (i) If ti+k = tk and if li+k = lfc for some k > 2, then Q has a rotational symmetry for the angle 360o/(r/k). (ii) If ca(Q) = (t1,... ,tr) = (tr,... ,t1) and if the cyclical sequence ci = (l1, l2,..., lr) of the lengths li of the boundary edges AiAi+1 has a reflection symmetry, then Q has a reflection symmetry. Proof. (i) This is obviously true for k = 1, for in that case we have ca(Q) = (t,t,... ,t), hence all ti are the same, and all li are the same, hence Q is a regular polygon Pr invariant for the rotation for the angle 360o/r. In the case k = 2 we have ca(Q) = (t, s,t,s,...,t, s). Removing from Q the congruent triangles AA1A2A3, AA3A4A5,..., AA2n-1A2nA1 we get a region with r/k = r/2 boundary edges where all ti are the same and all li the same (case k = 1) and for which we know (i) is true; hence Q has the rotation for the angle 360o/(r/2), too. Similarly, for the k > 3 we remove r/k = n congruent triangles from Q to get a smaller region with r — r = nk — n = n(k — 1) boundary angles and with a period k — 1, hence by the induction hypothesis having a rotation Rn where n = k = "^l^. Hence Q has the rotation for the angle 360o/(r/k), too. (ii) Reflection symmetry of the sequence cl (Q) may have three forms2: (a) ci = (y,z,...,b,a,a,b,..., y, z), (b) ci = (x,y,z, .. . ,b,a,a,b,.. ., y, z), (c) ci = (x,y,z, .. . ,b,a,a,b,.. ., y, z, w), and the same holds for the cyclical sequences of the lengths li. Obviously (ii) is true if Q has 3,4, 5 or 6 sides, since the only possible cases of cl (Q) are: (x, y, y), (y, z, z, y), (x, y, y, z), (x, y, z, z, y), (x, y, z, z, y, w) and (x, y, z, z, y, x). Now we can use a simple induction argument to see that if Q has more sides than 6, then we can cut it in two pieces with the boundary angles code of the types (c) or (b), and each of these two pieces has the same reflection symmetry (whose axis is the symmetral of the same line XY), hence Q has the same reflection symmetry (see Figure 7 middle). This is clear for regions of the type (a); the regions of the types (b) and (c) are obtained from a region of the type (a) by glueing one or two triangles to it. □ Theorem 3.12. Let M = M(Q) be a regular planar polygonal system without triangles covering the region Q. Then the symmetry group S(M) of M and the center of the rotation R" are determined by its boundary angles code ca(M) as follows: (i) if ca(M) has a reflection symmetry ca(M) = (t1,... ,tr) = (tr,... ,t1) then M has a reflection symmetry; 2The letters x,y, z,a,b,... are used here to denote lengths li of boundary edges of Q. 166 Ars Math. Contemp. 16 (2019) 157-171 Figure 7: Left: a region with rotational symmetry; middle: a region with reflection symmetry; right: a region whose all boundary angles are identical (equal to 120°, as in a regular hexagon), yet it has no symmetry. (ii) if ca(M) = (tl,...,tr) is a periodical cyclical sequence: ti+k = ti for some k > 2 then M has a rotational symmetry for the angle 360° /(r/k); (iii) if f = 1 (mod n) then the center of the rotation Rn is in a face; (iv) if f = 0 (mod n) and n > 2 then the center of the rotation is in the vertex (if n > 2); (v) if f = 0 (mod n) and n = 2 then the center of the rotation is in an edge. Proof. By Theorem 3.5, M = M(Q) is uniquely determined by its boundary angles code ca(M) = ca(Q). By Theorem 3.9, the group of symmetries of M is isomorphic to the group of symmetries of the boundary of Q. Now (i) and (ii) follow directly from Lemma 3.11. Indeed, if the boundary angles code ca(M) = ca(Q) of length r remains the same after the shift k (mod m), then r = kn is divisible by k and Q has a rotational symmetry Rn (see Lemma 3.11). Likewise, the reflection symmetry of ca (Q) implies the reflection symmetry of Q, since all the edges of Q are of the same size (since they are tiled by the regular polygons of M joined edge to edge) and we can apply Lemma 3.11. Now we use again S(Q) = S(M), implied by Theorem 3.5. The center of the rotational symmetry Rn of any planar regular polygonal system can be in a vertex (this is possible only in the cases when n = 3,4,6), in an edge center (this is possible only if n = 2) or in a center of a m-gon, where n divides m. If there are no triangles in the system, and if n is odd, then the center of the rotation Rn is in a face center. If the rotational symmetry Rn has the center in a m-gon, then fm = 1 (mod n) while the number of all other i-gons fi is divisible by n, hence f = 1 (mod n). If the center of the rotation is in a vertex or in an edge center, then fi = 0 (mod n) for every i, hence f = 0 (mod n). By Theorem 3.8 the face vector of M without triangles may be obtained from the boundary angles code ca(M) so the numbers fi are known, hence we can easily distinguish between the two possible cases with rotational symmetry: f = 1 (mod n) and f = 0 (mod n). If f = 0 (mod n) and n > 2 then the center of the rotation must be in a vertex. If f = 0 (mod n) and n = 2 then the center of the rotation must be in an edge. □ Remark 3.13. So the information hidden in the boundary angles code suffices to find out whether the center of rotation is located in a face, in a vertex or in the middle of an edge.3 3Another possible approach to this question is to calculate the center of rotation directly via orbit barycenters, and compare its coordinates with the coordinates of all the vertices, face centers and edge centers of the tiling M. J. Kovic: Regular polygonal systems 167 4 Summary and open questions We proposed a concept of a regular polygonal system, a mathematical model of chemical (planar or non-planar) molecules composed of chains of atoms forming regular n-gonal cycles of equal or various length (n = 3,4, 5, 6,... ). We proved that the structure (and hence the symmetry group) of any regular planar polygonal system M without triangles or without squares and hexagons is reconstructible from its boundary angles code ca(M ). The proof of Theorem 3.1 implicitly uses a simple incommensurability argument. This type of argument can be generalized as in the following definition and lemma: Definition 4.1. The quantities q1,q2,... ,qn are called incommensurable, if they satisfy the following condition: if two linear combinations of these quantities with integer coefficients are the same, i.e. J2n=i aiqi = n=i Mi, then their corresponding coefficients must be the same, i.e. ai = bi. From this definition immediately follows: Lemma 4.2. If any quantity q can be expressed as an integer linear combination of incommensurable quantities, then the coefficients ai are uniquely determined by q. Conjecture 4.3. The areas of all planar regular polygons with the same side length (except the triangle or hexagon) are incommensurable quantities. Conjecture 4.4. The volumes of all Platonic and Archimedean solids with the same side are incommensurable quantities. If one could prove Conjecture 4.3, this would automatically prove also Theorem 3.8. However, the converse is not true. If one could prove Conjecture 4.4 this would prove another conjecture: Conjecture 4.5. If we glue together copies of Platonic and Archimedean solids (with unit sides) face to face then we can find the number of each of them just by knowing the volume ofsuch composed solid. References [1] H. Abelson and A. diSessa, Turtle Geometry: The Computer as a Medium for Exploring Mathematics, MIT Press Series in Artificial Intelligence, MIT Press, Cambridge, Massachusetts, 1981. [2] A. T. Balaban, Chemical graphs-VII: Proposed nomenclature of branched cata-condensed benzenoid polycyclic hydrocarbons, Tetrahedron 25 (1969), 2949-2956, doi:10.1016/ s0040-4020(01)82827-1. [3] M. Deza, P. W. Fowler and V. P. Grishukhin, Allowed boundary sequences for fused polycyclic patches and related algorithmic problems, J. Chem. Inf. Comput. Sci. 41 (2001), 300-308, doi: 10.1021/ci000060o. [4] B. Grunbaum and G. C. Shephard, Tilings and Patterns, W. H. Freeman and Company, New York, 1987. n i= 1 168 Ars Math. Contemp. 16 (2019) 157-171 [5] I. Gutman and S. J. Cyvin, Introduction to the Theory of Benzenoid Hydrocarbons, SpringerVerlag, Berlin, 1989, doi:10.1007/978-3-642-87143-6. [6] P. Hansen, C. Lebatteux and M. Zheng, The boundary-edges code for polyhexes, J. Mol. Struct. THEOCHEM 363 (1996), 237-247, doi:10.1016/0166-1280(95)04139-7. [7] J. Kovic, How to obtain the number of hexagons in a benzenoid system from its boundary edges code, MATCH Commun. Math. Comput. Chem. 72 (2014), 27-38, http://match.pmf.kg. ac.rs/electronic_versions/Match72/n1/match72n1_2 7-38.pdf. [8] J. Kovic, T. Pisanski, A. T. Balaban and P. W. Fowler, On symmetries of benzenoid systems, MATCH Commun. Math. Comput. Chem. 72 (2014), 3-26, http://match.pmf.kg.ac. rs/electronic_versions/Match72/n1/match7 2n1_3-2 6.pdf. [9] N. J. A. Sloane (ed.), The On-Line Encyclopedia of Integer Sequences, published electronically at https://oeis.org. J. Kovic: Regular polygonal systems 169 Appendix A Planar regular (3,4)-systems These systems may be classified by two parameters: the numbers /3 and /4 of triangular and square faces, as in Figure 8. As we see, the number of these systems grows very quickly with the total number of faces / = /3 + /4. There are 16 such systems with 1, 2 or 3 faces. In Figure 9 they are arranged by the increasing total number of faces. We give also their boundary faces-edges codes cf,e and boundary angles codes ca. Observe that in all these cases all the faces are boundary faces. Since for the planar regular polygonal ( 3,4)-systems all their boundary angles are multiples of 30°, we give the boundary angles codes also in the form of multiples of 30°. The corresponding boundary faces-edges code and boundary angles codes of these ( 3,4)-systems, ordered as in Figure 9, are: cf,e (M ) c0(M ) (33 ) (120,120,120)° = 30° 4,4,4) (44 ) (90, 90, 90, 9o)° = 30° (3, 3, 3, 3) (32 32) (60,120, 60,120)° = 30° (2, 4, 2, 4) (43 32) (120, 30, 90, 90, 30)° = 30° 4,1, 3, 3,1) (43 43) (90, 0, 90, 90, 0, 90)° = 30° (3, 0, 3, 3, 0, 3) (32 3i, 32 3o) (60, 60,120,0,120)° = 30° (2, 2, 4, 0, 4) (32 3o, 43 3i) (-30, 90, 90, 30, 60,120)° = 30° -1, 3, 3,1, 2,4) 32 , 3i, 43, 3o) (90, 90, -30,120, 60, 30)° = 30° 3, 3,-1,4, 2,1) (32, 4i, 32, 4i) (120, 30, 30,120, 30, 30)° = 30° 4,1,1, 4,1,1) (32, 4o, 32, 42) (120, -30,120, 30, 90, 30)° = 30° 4,-1, 4,1, 3,1) (43, 4o, 32, 42) (90, -60,120, 30, 90, 0, 90)° = 30° 3, -2, 4,1, 3,0, 3) (42, 32, 4o, 43) (30,120, -60, 90, 90, 0, 90)° = 30° 1,4, -2, 3, 3,0, 3) (43, 4i, 32, 4i) (90, 0, 30,120, 30, 0, 90)° = 30° 3,0,1, 4,1,0, 3) (43, 3o, 43, 3i) (90, 90, -60, 90, 90, 30, 30)° = 30° 3, 3, -2, 3, 3,1,1) (43, 4i, 43, 4i) (90, 90, 0, 0, 90, 90,0,0)° = 30° 3, 3, 0, 0, 3, 3, 0, 0) (43, 4o, 43, 42) (90, 90, -90, 90, 90, 0, 90,0)° = 30° 3, 3, -3, 3, 3,0, 3, 0) /4 =0 /4 = 1 /4=2 /4=3 /3 = 0 □ m irr: r /3 = 1 A ù àitû m> etc. Figure 8: Regular planar ( 3,4)-systems classified by two parameters. 170 Ars Math. Contemp. 16 (2019) 157-171 au & fatn & & q rf) fVi m> -OO* rm:!; <6 Figure 9: Planar regular polygonal (3,4)-systems with f = f3 + f4 < 3 faces. Appendix B Types of faces The idea to classify the types of hexagons in a benzenoid system with respect to their contacts to adjacent faces in a system [2, 7] may be generalized to any n-gonal faces in any polygonal system as follows: Definition B.1. The type of the n-gonal face f in a polygonal system M is the cyclical binary sequence (bi,..., bn) where b = 0 if the ¿-th edge of f is incident to at least one other face of M and b = 1 if it is a boundary edge of f. The number of types of n-gons having k entries 1 in the binary code is denoted T(n, k). The number of possible types of a n-gon is denoted T(n). Theorem B.2. The following relations hold: (i) T (n) equals the number of binary cyclical sequences of length n. (ii) T(3) = 4, T(4) = 6, T(5) = 8, T(6) = 13. Proof. (i) This is obvious, and trivially implies also relations T (n, k) = T (n, n - k) and T (n) = T (n, 0) + T (n, 1) + • • • + T (n, n). To each type t of a n-gon exists the opposite type t* with the entries 6* = 1 - 6j, hence t** = t. For (ii) see Figure 10. Summing the adjacent entries 1 and ignoring the entries 0 we can, at least for the triangles T, squares S and pentagons P, denote their possible boundary J. Kovic: Regular polygonal systems 171 types like this: To = (0,0,0), Ti = (1,0,0), T2 = (1,1,0) = T1Î, T3 = (1,1,1)= T0, So = (0,0, 0, 0), S1= (1, 0,0,0), S2 = (1,1,0,0) = S2, Si+i = (1,0,1, 0) = SÎ+1, S3 = (1,1,1,0) = SÎ, S4 = (1,1,1,1), Po = (0,0, 0, 0,0), P1= (1, 0,0,0,0), P2 = (1,1,0,0, 0), P1+1 = (1,0,1,0,0), P1+2 = (1, 0,1,1,0) = Pj P3 = (1,1,1,0, 0) = P2*, P4 = (1,1,1,1,0) = Pj\ P5 = (1,1,1,1,1). The inner faces are T0 = (0,0,0), S0 = (0,0, 0, 0), P0 = (0,0,0,0, 0). The 13 types of boundary hexagons H with 1,2,3,4 or 5 boundary edges, e.g. Hi+i+i = (1,0,1,0,1,0) or H2 = (1,1,0,0,0,0) are given in [2] and more precisely classified in [7]. □ Remark B.3. The sequence T(n) (although defined only for n > 3) corresponds to the sequence A000029 in OEIS [9].