The map proposed by R. Lozi over 40 years ago is one of the simplest two dimensional maps that exhibits chaotic behavior and generates a wide range of interesting structures. The map may be defined thus:
where are real parameters.
We have discussed our explorations of this map at some length before. At that point we also discussed our discovery of certain generalizations of this class of maps. In course of the quest for the generalized Lozi-like maps, we discovered another generalization remarkable for the diversity of forms it produced. We held back from presenting it because we wished to investigate its geometry further. As a result we are now able to explain the basic geometric determinants of the map; however, there are other aspects that still remain mysterious to us. Hence, we are presenting this map with whatever understanding we could arrive at along with the issues that are still open questions for us. They are maps of great beauty; hence, they are also worth beholding for the sake of their aesthetics.
This map is defined thus:
Alternatively, we can also define it as:
Here are real parameters: is in the range (see below for explanation), while we explored primarily in range as this is the range where there is higher tendency for the map to yield chaotic attractors. If we take the first definition then the map is bilaterally symmetric about the axial line . It we take the second definition it is similarly symmetric about the axis . For the below discussion we refer only to the first definition because the second one shows comparable behavior with change in signs of some of the parameters.
To explore the basic geometric principles behind the attractors generated by this map we started with a set of 200 equally spaced starting points all of which lie on the circle . With each starting point the map was iterated 2000 times. Then the evolution of each starting point under the map was plotted in one of 11 different colors. We were able to determine the following rules for the structure of these attractors:
1) The most important determinant of the geometry of the attractor is sum of the parameters . When are closer to 1 and where are mutually prime integers, we get a p-ad structure in the attractor (e.g. Figure 1).
In Figure 1, , and takes successively, by row, the following 4 pairs of values (-0.5450419, 0.1); (-0.2450419; -0.2); (-0.7450419, 0.3); (-0.9450419; 0.5) . In each case . This determines the central heptad structure but the rest of the attractor shows considerable variability that cannot be explained in a straightforward way from these parameters.
If and then there is greater tolerance for the range of values can take.
With we get a nonad structure even though is relatively low.
When and the p-ad structure in the attractor tends to be retained (Figure 3).
In this case is close to ; thus, we get a core with a heptad structure.
Irrespective of , if are close to each other in magnitude but opposite in sign then we get square structure in the core of the attractor (Figure 4). This is because , which is
2) For low values of both or at least , including when they are below 0 and we see a loss of the p-ad structure or a flattening of the attractor along the axis of symmetry from the top left side (Figure 5).
In both panels ; . In the left panel and the expected pentad structure is retained. In the right panel ; we notice that the pentad structure has been lost with flattening along the axis of symmetry.
These are the basic rules by which the core structure of the attractors can be accounted for. The parameters can be seen as scaling parameters whereas can be seen as rotational parameters. The latter pair sets up the rotational structure based on how close their sum approaches . This is the cosine principle that occurs as the primary structure-determinant in other chaotic maps, including maps which we discovered and discussed before (1, 2). As in the previously discussed cases, this sets the limits of as (-2, 2); however, here there is greater scope for diversity because the cosine principle is distributed over the sum of two separate parameters. Beyond the cosine principle,n these attractors exhibit additional features whose origins remain mysterious to us:
1) For example, In Figure 3 in addition to the core heptad one can see additional harmonics like a 11-ad, 18-ad and 29-ad structures. What is the explanation for them? In the previous examples of chaotic maps we could explain higher harmonics based on other whose might lie close to the primary cosine. However, in this case the 11-ad, the 18-ad and 29-ad correspond to no such cosines. Hence, their emergence remains a mystery.
2) Further, there is an effect of the starting , which can result in the appearance of a n-ad structure independently of the cosine principle (Figure 6).
In both the above cases the four parameters are respectively, -0.567115129642189, -0.761931293738961, -0.347582505706949, -0.069076271019876. However, the left panel was initialized with a circle of radius 0.9 and the right panel was initialized with a circle of radius . Only in the second case we see the emergence of a heptad structure, which cannot be accounted for by the sum of . The role of the starting in generating n-ad structures is another open question.
To further explore the diversity of chaotic attractors within the above-stated parameter space for the 4 parameters of the map we set up our code to search for such maps thus:
1) We set as the starting point.
2) We then randomly generated a set of parameters in the range .
3) For each such set we allowed to evolve for 1000 iterations under the map. If the evolution resulted in the divergence to or convergence to one or few point attractors, then we discarded those parameters.
4) A heuristic for chaos is that very small differences in the initial conditions can result in very different end results upon evolution under the map after a certain number of iterations. Hence, for those parameters which survived the above filters we tested if a second , which differed from the first by in each coordinate diverged from the evolutionary path of the original after the same number of iterations. If it did so, the parameter set was retained as it was the sign of being on a chaotic attractor.
5) These surviving parameter sets were then explored more fully by studying the evolution of 200 equally spaced points on the circle for 2000 iterations of the map.
We present below few examples of attractors emerging from the above procedure with parameters listed in order from to .
Figure 7. -0.21947270759847, 0.33858182718046, 0.610763950282708, 0.315013298424892
This attractor assumes of the form reminiscent in some ways of the gingerbread man seen in the classical Lozi attractor. The gingerbread man has a “heart” in the form a five-lobed structure with a ellipse within it.
Figure 11. -0.482874695337377, -0.819933905000798, -0.381899853958748, 0.100596646997146
This is a representative of a prevalent type of attractor generated by this map that may be termed the “Sombrero hat” type. The central heptad structure in this attractor is a mystery because it is not explicable by the cosine conditions above-presented .
Figure 12. -0.442830848232843, 0.0637841782858595, 0.495140142193995, 0.24925118080806
This form of the attractor resembles the butterfly attractor generated by the square-root modification of the Lozi map that we had described earlier.
Figure 14. -0.517300257436, 0.795868384265341, 1.00516145080794, 0.493695715968497
For this attractor the radius of the circle was changed from 0.3 to 1, though it is stable and similar at the former value too. It vaguely resembles some Cambrian animal.
Figure 15. -0.307423420087434, -0.0870326082082465, -0.337442451966926, 0.766096443301067
This attractor assumes a bun-shaped morph that is seen quite often under these maps. An explanation for the multiple harmonics of this attractor remains as yet mysterious to us.
Figure 16. -0.0823114865785465, -0.766779959597625, 0.868352874149568, -0.0100218075420707
Another Sombrero hat type form. The radius of the circle was changed from 0.3 to .4 for this attractor. How one accounts for the central pentagonal zone of restriction and the triad of octagonal restriction zones around it remains mysterious.
Figure 17. -0.16997885087505, -0.0931982280220837, -0.631216614013538, -0.248738911962137
Here the radius of the circle was changed from 0.3 to 0.31. It shows a central pentad surrounded by 9 further pentads. An explanation for the emergence of these pentads and the 9 fold harmonic however remains elusive.