Polygons and Chaos

Year: 2001 Authors: Jay Kappraff; Gary W. Adamson

Core claim

Star n-gon edge cycles encode periodic dynamics of the logistic equation and two polynomial families, with each odd n-gon having a characteristic cycle length.

Topics

chaos, star polygons, periodic trajectories, number theory

Domains

dynamical systems, logistic equation, Chebyshev polynomials, Euler phi-function, geometry, music theory, visualization

Methods

semi-conjugacy, modular arithmetic, polynomial mapping, cyclic edge construction

Media

star polygons, cyclotomic n-gons, figures, tone circle

Paper text

The text below is the locally extracted OCR/Markdown version of the paper. Raw PDF files remain local and are not published here.

BRIDGES Mathematical Connections in Art, Music, and Science

Polygons and Chaos

Jay Kappraff Department of Mathematics New Jersey Institute of Technology Kappraff@aol.com

Gary W. Adamson P.O. Box 124571 San Diego, CA 92112-4571

Abstract

An infinite number of periodic trajectories are derived for the logistic equation of dynamic systems theory at a value of the parameter corresponding to the extreme point on the real axis of the Mandelbrot set. Beginning with the edge of a family of star -gons as the seed, the trajectory of the logistic map cycles through a sequence of edges of other star -gons. Each -gon for odd is shown to have its own characteristic cycle length. The logistic map is shown to be the first of two infinite families of maps, all exhibiting periodic trajectories, derived from two families of polynomials, the Chebyshev polynomials and another related to the Lucas sequence. These dynamics are shown to be closely related to properties of number.

1. Introduction

This paper describes a remarkable connection between the edges of star polygons and dynamical systems in the state of chaos. A sequence of dynamical maps are derived from the Chebyshev polynomials and another family of polynomials related to the Lucas series, and these maps exhibit periodic trajectories of all lengths with each polygon having its own characteristic cycle length. The first polynomial of the family is the logistic equation at a value of its parameter corresponding to the extreme left-hand point on the real axis of the Mandelbrot set. This leads to new connections between chaotic dynamics and both Euclidean geometry and the theory of numbers.

2. Star Polygons

Consider the polynomial equation for complex. The solutions are,

referred to as the -th roots of unity. The points in a cartesian coordinate system,

lie at the vertices of a regular -gon with unit radius that we shall refer to as a cyclotomic -gon. If the point is distinguished, the other points satisfy the equation,

68 Jay Kappraff and Gary W. Adamson

referred to as the th cyclotomic polynomial when is prime. We shall continue to use this terminology for non-prime but odd.

Consider the transformation:

This transformation has the effect of rotating clockwise by degrees, i.e., . If is taken as the seed of this transformation, then, where the sequence of points, for is the trajectory of the transformation.

Let for and relatively prime integers, and take the seed value to be ,

i.e., , then the trajectory forms a modular system with indices and principal values . The vectors represent the directed edges from vertex to and the system of edges, for are the edges of a regular star polygon symbolized by . The star polygons are denoted simply as and represent regular -gons. For 0 < m \leq \left\lfloor \frac{n-1}{2} \right\rfloor the edges in the cycle are clockwise and have positive lengths. For \left\lfloor \frac{n-1}{2} \right\rfloor \leq m < n, the stars are counterclockwise or retrograde and have negative edge lengths. If the star -gon has an edge cycle oriented clockwise then has retrograde edges of the same length, i.e., . The three positively oriented star 7-gons are shown in Figure 1.

img-0.jpeg Figure 1: The three positively oriented 7-gons.

If and are not relatively prime, then in lowest terms and reduces to the star -gon . The case of the six positively oriented star polygons for are shown in Figure 2. Only and 5 correspond to star 12-gons; is a square while an equilateral triangle, and a digon (a degenerate case of a polygon).

Polygons and Chaos 69

img-1.jpeg

img-2.jpeg

img-3.jpeg

img-4.jpeg

img-5.jpeg

img-6.jpeg Figure 2: The family of star 12-gons. Only and are star 12-gons. The star polygons is a 4-gon, is a 3-gon, and is a digon.

The star polygons corresponding to and all of its factors constitute a family of star -gons. In general, star polygons are related to any -gon if the digon and the polygon with a single vertex are included. This follows from the property of numbers that,

where the summation is over all of the factors of . is the Euler phi-function which equals the number of integers relatively prime to and not greater than . For example, the family of star 12-gons is represented by the factor tree,

where , , , , , with a sum of 12.

70 Jay Kappraff and Gary W. Adamson

3. -Cycles

Consider the (complex) transformation restricted to the unit circle which can be denoted by . Observe that this map takes into . This is the defining property of the Chebyshev polynomials , viz.,

In dynamical system terms, there is a semi-conjugacy (see Appendix A) defined by,

from the map to the map , . It should be noted that much of what we discuss in the sequel can be inferred from this semi-conjugacy.

Beginning with the seed where ,

where the sequence for is the trajectory of .

If , then . Since , when is the smallest positive integer such that,

In this case the trajectory forms a cycle of length , or what we shall refer to as a -cycle. Since , Equation 6 is equivalent to

If is a prime number, then so that must be a factor of [1]. If is odd but not prime, then is replaced by where is the Euler phi-function. Table 1 shows the values of for and 25 and and 17.

Polygons and Chaos 71

Table 1. Exponents such that

m2n = 7n = 9n = 11n = 13n=15n=17
4335644
93538
16335322
2533528

4. Polygons and Chaos

For and relatively prime integers and odd, given a value of , , it can be shown that there exists a value of , such that,

for and relatively prime (e.g., cannot be even). The values, , are the edge lengths of either of the star -gons or depending on their signs [2]. Therefore we can view the -cycles of the previous section corresponding to an -cyclotomic polynomial as being a sequence of edges of star -gons.

This encourages us to consider the transformation of the circle of radius 2, , defined by,

This map takes into , which is also the defining property of the Lucas polynomials [3],

In other words, the maps restricted to are semi-conjugate (see Appendix A) to the Lucas polynomials for via the map

This has some interesting dynamical consequences.

72 Jay Kappraff and Gary W. Adamson

Letting , as the result of Equations 5b and 9b, and transformations 5a and 9a can be rewritten as,

where the first few Chebyshev and Lucas polynomials are listed in Table 2,

Table 2. The Lucas and Chebyshev Polynomials

L1x1T1x1
L23T23
L34T37
L47T417
L511T541
L618T699
L729T7239

Notice that, ignoring signs, a coefficient in row and column of Table 2 is given by the recursion relation, for the coefficients of the Lucas polynomials and for the coefficients of the Chebyshev polynomials. Also note that, ignoring signs, the sum of the coefficients of the Lucas polynomials form the Lucas sequence which satisfies the recursion relation: , i.e., it is a Fibonacci sequence beginning with 1,3. On the other hand the Chebyshev polynomials form the Chebyshev sequence (-sequence) which satisfies the recursion relation: , i.e., a Pell sequence beginning with 1,3. The standard Fibonacci sequence (-sequence) is: and while the standard Pell sequence (-sequence) is: and . Also,

where and . Also and . For this reason and are referred to as the golden and silver means respectively.

Consider the second Lucas polynomial , and its iterative map,

also known as the logistic map. It represents the extreme left-hand point on the real axis, i.e., , of the Mandelbrot set (see Figure 3) for the map,

Polygons and Chaos 73

where is a complex number. The arrow notation in Equation 11 means that we choose a seed value and place it into the polynomial to get . From we get , etc. and in this way generate the sequence , the trajectory of the map.

The map given by Equation 11a is a transformed version of the logistic map

which has been studied in great detail [4], [5]. The fact that in Equation 11b, means that this map is in a state of chaos. It can be shown that for values of \lambda < 4 (or c > -2) all points on the unit interval are “imprisoned” in the sense that their trajectories remain in the unit interval [1,0] for Equation 12 or for values of corresponding to in Equation 11a, the trajectories remain on the interval . However, beginning at (or ), orbits can escape; in fact the only imprisoned orbits lie on a Cantor set within the unit interval [0,1]. For any complex value of , the boundary in the complex plane of the prisoner set is what is called the Julia set. Therefore, the Julia sets for real values are what we refer to as “Cantor dusts.”

img-7.jpeg Figure 3: The Mandelbrot set.

The theory of dynamical system shows that as is increased to the Feigenbaum limit 3.569… the trajectories of the system go through period doubling bifurcations, i.e., cycles of length for . At the value 3.831… a trajectory with a cycle of length 3 appears, after which periods of every length are present according to the theorem of Sharkovskii. As a result of our analysis when is further increased to a value of 4, or alternatively is decreased to , the cycles can be characterized as edges of star -gons for odd in which each value of has its own characteristic cycle length. Therefore, in a sense, the edges can be thought to dance about on the grains of a Cantor dust.

74 Jay Kappraff and Gary W. Adamson

4. Results

For any value of the and transformations defined by Equation 10 result in -cycles, some of which are listed in Table 1, of semi-edges and edges, respectively, of star -gons. We shall study the dynamics of the logistic map of Equation 11a corresponding to , in some detail giving the results of the -cycles for cyclotomic -gons where and 17.

For the cyclotomic 7-gon, values of for are listed in Table 3 along with the values of for which (ignoring signs). In other words, a given value of corresponds to the edge of the star 14-gon. Also listed in Table 3 is the actual edge corresponding to taking into account its sign. For example, for which corresponds to for . Therefore, this represents the edge of the species of retrograde star 14-gon.

If is taken to be the seed of Transformation 11a, then we find the remarkable result that, the iterates are the sequence of edge lengths of different species of star 14-gons,

and since the sequence repeats with the 3-cycle,

As a result of the fact that , positive and negative values of the edges corresponding to (mod ) are identical. As a result, which corresponds in Sequence 13 to , and so the 3-cycle is represented by the -values in Table 3 and therefore by the sequences of edges . The orders of the edges in the cycles are also indicated in Table 3, beginning with the seed and in Figure 4 by the boxed integers. These edges are shown in Figure 4 drawn from vertex number 0 of the 14-gon to the darkened vertex . Notice the regular skip pattern of highlighted vertices: 4, 4, 6. We shall find this pattern to hold for all cycles. Also note that for each cycle, its mirror image, illustrated by open vertices, within the -gon is also a cycle, e.g.,

is another 3-cycle for the cyclotomic 7-gon. Each cycle of a cyclotomic -gon will have a corresponding mirror image cycle. Finally, Adamson has discovered that the product of the cycle values equals -1, e.g.,

and that the product of edges within any cycle will always equal for , a result that can be proven from a dynamic/number theoretic description.

Polygons and Chaos 75

img-8.jpeg Figure 4: The points shown with closed circles within the following 14-gon represent the edge-lengths from point 0 to point of the star 14-gon {14/k}. These are elements of the trajectories of the logistic equation corresponding to that 14-gon. The open circles represent the mirror image trajectory. The numbers in square boxes are the order of the points in the trajectory.

Table 3 Cycles for the Logistic Equation Corresponding to Cyclotomic Polygons.

Cyclotomic 7-gonCyclotomic 9-gon
K2cos 2πk/72sin πj/14{14/j}Order2cos 2πk/92sin πj/18{14/j}Order i
11.24696…3301.53208…550
2-0.44509…11310.34729…111
3-1.80189…592-0.500…
4-1.87938…7112
Cyclotomic 11-gonCyclotomic 17-gon
---------------------------
K2cos 2πk/112sin πj/22{14/j} jOrder I2cos 2πk/172sin πj/34{34/j} jOrder i
11.68250…7701.8649…13130
20.83082…3311.4780…991
3-0.28462…12130.89147550’

76 Jay Kappraff and Gary W. Adamson

4-1.30972…51720.1845…112
5-1.91898…9134-0.5473…3312’
6-1.2057271’
7-1.70043…11233’
8-1.9659…15193

Cyclotomic 13-gon

K2 cos 2 sin {26/j} jOrder I
11.77090…990
21.13612…551
30.24107…114
4-0.70920…3232
5-1.49702…7193
6-1.94188…11155

Next consider the cyclotomic 9-gon. Its values of are listed in Table 3 along with the corresponding values of for which . Using the value of as the seed in the logistic Map 11a, the corresponding 3-cycle is derived, as before, from the sequence of -values: 1, 2, 4, 8, … (mod 9) where which corresponds to the edge in Table 2. Therefore the 3-cycle is represented by values: . Notice that is not part of the cycle because 3 is not relatively prime to 9. The corresponding value gives rise to the edge , and this is the edge of a star 6-gon not an 18-gon. The sequence of edges is: (not shown). This time the skip pattern is: 4,8,6 (the pattern would have been 4,4,4,6 if were to be included) and the product of the edges in the cycle equals to .

As a final example consider the cyclotomic 17-gon. Its and values are listed in Table 3 and its 4-cycle can be generated from the seed value corresponding to as the sequence, once again, of -values: 1, 2, 4, 8, 16,… (mod 17) where which corresponds in Table 2 to the edge . This leads to the sequence of star 34-gon edges corresponding to -values: … whose edge lengths are found in Table 3, and to the cycle of edges:

However, since 17 is a prime number there are eight distinct star 34-gons (not considering orientation) and we have accounted for only four of them where 4 is a factor of 8 according to the results of Section 1. The other four can be obtained by beginning with a seed value corresponding to resulting in the sequence: 3, 6, 12, 24, 48, … (mod 17) where which corresponds in Table 2 to the edge . But since corresponding to the edge and , this leads to the sequence of edges corresponding to the values: found in Table 3 and to the edge sequence: . This 4-cycle is distinguished from the other in Table 3 by order numbers denoted by primes (’). If the vertices corresponding to both of these edge cycles are highlighted (not shown), they lead to the skip pattern: 4, 4, 4, 6, 4, 4, 4, 4. Once again, the product of the edges equals -1 for the first 4-cycle and +1 for the second.

Polygons and Chaos 77

5. A Cycle Algorithm

These examples lead to a simple algorithm for any cyclotomic -gon for odd:

  1. Compute for , .
  2. The k-values for the cycle related to the cyclotomic -gon is the sequence: , disregarding signs, i.e., for . This sequence repeats when (mod ) for in which case the cycle has length and the values of give the order of terms in the sequence beginning with a seed value of corresponding to . Record the order numbers for the sequence of iterates to the logistic Equation 3 as .
  3. Once the -values of the cycle are determined, compute the values of such that , ignoring signs. The values of come from the sequence: where is relatively prime to .
  4. If the sign of the cosine is positive, then the species of star polygon is ; if the sign is negative then the star polygon is the retrograde version, .
  5. The process ends when all values from step 3 are accounted for. If all of the sine values have not been used, then a second sequence of -values gives another cycle of length given by: for where for and steps 2-4 are repeated. Record the order numbers as . If there are still additional unaccounted for -values, then another -cycle can be determined from the sequence , etc.
  6. The period length is a divisor (or factor) of the number of integers , relatively prime to .
  7. Observe that the product of edges from a -cycle satisfies the equation,

Try applying this algorithm to . Note that among the integers only 1,2,4, and 7 are relatively prime to 15, so that according to step 6 of the Algorithm, the length of the cycle can be either 2 or 4. However, since (mod 15), according to step 4 of the Algorithm the period is 4.

The results of applying this algorithm to the cyclotomic 11- and 13-gons are found in Table 3 and in Figures 5a and b. The cyclotomic 11-gon results in a 5-cycle of edges of the star 22-gon family while the cyclotomic 13-gon results in a 6-cycle within a 26-gon. The equal-tempered chromatic scale can be represented by a tone circle with 12 tones to the octave, or a 12-gon with each tone distant from the next by a semitone. Therefore a 24-sided polygon can be thought of as a tone circle in which each tone represents an interval of a quarter-tone. This means that the 5- and 6-cycles of the cyclotomic 11- and 13-gon, along with their symmetric opposites can be viewed as tonal subscales of almost quarter-tone intervals, one with tones slightly greater than quarter-tones and the other with tones slightly less.

Jay Kappraff and Gary W. Adamson

img-9.jpeg

img-10.jpeg Figure 5: Trajectories of the logistic equation corresponding to the edges of (a) a star 22-gon; (b) a star 26-gon (see figure caption of Figure 4).

Polygons and Chaos 79

6. Other Lucas Polynomials

What we have discovered for the second Lucas polynomial, appears to hold true for all of the other Lucas polynomials from Table 2. For example, the third Lucas polynomial with alternating signs leads to the recursive map,

With , taking the seed for odd , results in the sequence where this time (mod ), i.e., (mod ) for , a -cycle as predicted by Table 1. The sequence is based on powers of 3 since we are using the Lucas polynomial. The cycle lengths can be determined by the Algorithm with the only change being that -values of the cycle are now powers of 3 instead of 2. This generalizes to the -th Lucas polynomial with alternating signs in which case the iterates correspond to -values that are powers of . Presumably, these polynomial maps also represent dynamical systems in a state of chaos. However, the products of the edges in a cycle are not always equal .

7. Chaos and Number

There is an intimate relationship between chaos theory and number. We have shown that properties of number also lie at the basis of the polygon cycles. Let us once again consider the Lucas polynomial map. Expanding 1/7 in the base 4 (the square of 2), , a repeating decimal with a 3-cycle. Likewise . expanded in base 4, a 5-cycle (see Appendix B). Our conjecture is that for odd, expanded in base 4, has the identical cycle length as the cyclotomic -gon analyzed in the previous section. Furthermore, the identical cycle lengths occur for in base 9, 16, 25 or any base as for the cycle lengths of cyclotomic -gons corresponding to the -th Lucas polynomial maps as shown in Table 2. The validity of this claim and other parts of this analysis were computer checked by Malcolm Lichtenstein [6].

8. Conclusion

Many things have come together in this study. We have shown the close relationship between chaos and number. The well known theorem of Sharkovskii predicts that once a period 3-cycle appears in a dynamical system, periods of all lengths occur. We have shown that at critical point of the Mandelbrot set where orbits of the logistic equation begin to escape, each of these periods can be characterized by a sequence of edge lengths of a family of star -gons, each having a characteristic cycle length. Coxeter has shown star polygons to be related to the two-dimensional projections of higher-dimensional polyhedra or polytopes [7]. Geometry has shown itself once again to be the rich well-spring of mathematics. Rather than jettisoning these roots, the theory of dynamical systems and chaos has strongly embraced them.

Each of these star polygons can be looked at as a tone circle with the cycles represented by tones from the “octave.” In particular, the 5-cycle from the 22-gon and 6-cycle of the 26-gon are promising candidates for new musical scales. After all, the chromatic scale was built from the circle of fifths related to the star 12-gon.

80 Jay Kappraff and Gary W. Adamson

Bibliography

[1] L.Dickson. Introduction to the Theory of Numbers. New York: Dover (1957). [2] J. Kappraff. Beyond Measure: A Guided Tour through Nature, Myth, and Number. Singapore: World Scientific. In press. 2002. [3] Blackmore, D. Private communication [2001]. [4] H-O. Peitgen, H. Jurgens, and D. Saupe. Chaos and Fractals. New York: Springer-Verlag. 1992. [5] M. Schroeder. Fractals, Chaos, Power Laws. New York: W.H. Freeman Press. 1990. [6] M. Lichtenstein. Private communication. 2001. [7] H.S.M. Coxeter. Regular Polytopes. New York: Dover (1973).

Acknowledgement

I wish to acknowledge generous help that I have received from Denis Blackmore in strengthening the mathematical arguments of this paper and helping to place the results into the mainstream of research in dynamical systems, and to Robert Devaney for several of his helpful comments on this work.

Appendix A

Given a pair of transformations on set and on set , the transformation mapping to , such as the ones defined by Equations 5b and 9b, induce relationships between the transformations and via the following commuting diagram:

In this relationship where is not invertible. Such a transformation is known as a semi-conjugacy. When is invertible an isomorphism is induced between the and characterized by , and is called a conjugacy.

Appendix B

A decimal in base 10 can be written in any other base by the following procedure illustrated for converting in base 10 to the base 4.

  1. Multiply the decimal in base 10 by 4 and record a 0 if the result is less than 1, otherwise record the integer part. For example, so record a 0 as the in the first decimal place.
  2. Multiply the result again by 4 to get 0.2857142… and record a 2 as the decimal place.
  3. Multiply the decimal part of the preceding number by 4 to get: 1.1428568… and record a 1 as the decimal place.
  4. Again multiply the decimal part of the preceding number by 4, but since the decimal part repeats we have the repeating decimal in base 4: 0.021

In general, consider the rational fraction where is the decimal expansion of in base 10. Its decimal expansion in base is then: where equals the integer part of: for .

0 items under this folder.