Discriminating different types of chaos is still a very challenging topic, even for dissipative three-dimensional systems for which the most advanced tool is the template. Nevertheless, getting a template is, by definition, limited to three-dimensional objects based on knot theory. To deal with higher-dimensional chaos, we recently introduced the templex combining a flow-oriented BraMAH cell complex and a directed graph (a digraph). There is no dimensional limitation in the concept of templex. Here, we show that a templex can be automatically reduced into a “minimal” form to provide a comprehensive and synthetic view of the main properties of chaotic attractors. This reduction allows for the development of a taxonomy of chaos in terms of two elementary units: the oscillating unit (O-unit) and the switching unit (S-unit). We apply this approach to various well-known attractors (Rössler, Lorenz, and Burke–Shaw) as well as a non-trivial four-dimensional attractor. A case of toroidal chaos (Deng) is also treated.

The nonlinear essence of many (if not most) of the dynamical processes in nature unavoidably led to the discovery of chaos, a deterministic dynamics for which long term prediction is not possible. With the emergence of desktop computers, many different types of chaos were reported, but a complete taxonomy is still under development. One of the limitations of the topological characterization of chaotic attractors is the use of knot theory, whose applicability is limited to three-dimensional space. Templexes are objects allowing for a definition of an algebraic topology on a cell complex (which, unlike knots, has no dimensional restrictions) and a directed graph whose nodes are the locally highest-dimensional cells of the complex and whose edges indicate the way in which the flow visits such cells. This work shows that a simplified templex can describe various types of chaotic attractors as different combinations of two elementary dynamical processes, namely, oscillating units and switching units. Such units pop up after the reduction procedure that merges cells or nodes, without altering the topological properties. The way in which the templex is reduced is crucial in order to identify the basic units underlying a given dynamics. Different types of chaotic attractors are discussed.

Understanding the structure of chaotic attractors has practical applications in various fields, including physics, biology, economics, and engineering. For instance, chaos theory has been applied in fields, such as weather prediction1–3 and control,4 fluid dynamics,5–7 optics,8 secure communication,9 and the modeling of ecological systems.10,11 Many topologically inequivalent chaotic attractors are now known in almost any field. There is no general relationship between the algebraic structure of the governing equations and the type of attractors they produce. Nevertheless, in the early 1970s, Otto E. Rössler developed complex chemical reactions as the combination of some basic circuits,12 an approach inherited from his background in electronics.13 He was following the works of Sugita who developed an analogy between chemical systems and logical circuits.14,15 Thus, Rössler introduced his eponym system as a combination between a two-dimensional oscillator and a nonlinear switch.16 He described the resulting chaotic attractor as the union of a “normal” strip and a Möbius strip, that is, using a branched manifold. The topological approach, thus, appeared to be the adequate way to describe chaotic attractors and to classify them. He then provided an early attempt to construct a “hierarchy of chaos” based on positive Lyapunov exponents.17,18 This led to the recent taxonomy of chaos19 by using the number p of positive Lyapunov exponents and the number q of null Lyapunov exponents to construct the marker C pT q for discriminating different classes of attractors, C p discriminating periodic or quasiperiodic ( p = 0) from chaotic ( p = 1) or hyperchaotic ( p = 2) behaviors: T q directly designates the torus T q on which the structure of the attractor is based (for instance, a limit cycle with q = 1 and a torus for q = 2).

Topological characterization of chaos provides a deep understanding of the intricate and complex behavior produced by chaotic dynamical systems.20–23 It offers analysis and description of the underlying structure of chaos in a systematic and mathematical way.19 Topology is a branch of mathematics that deals with the properties of a spatial object that remains unchanged under isotopy, that is, continuous deformations, such as stretching, bending, or twisting. In the context of chaos theory, topological methods are used to analyze and describe the structure of chaotic attractors produced by numerical simulations as well as measured in experiments.24–30 For many years, such a characterization was mainly based on knot theory and the construction of a template as a knot-holder.31 A template is a branched manifold sketched with a splitting chart, some strips with some local torsions and permutations. These strips are finally squeezed at the so-called joining chart. The way in which strips are squeezed can follow a “standard” convention as introduced by Tufillaro and co-workers.20 Characterizing higher-dimensional attractors was, therefore, a great challenge in spite of a few attempts.32,33 Indeed, building a template is strongly related to the relative organization of unstable periodic orbits extracted from the chaotic attractors, an organization that is quantified using topological invariants, which are only valid in three-dimensional space.

A very promising approach consisted of borrowing concepts from homology group theory, such as Betti numbers,34 to study a trajectory in the state space or in a reconstructed space. This approach was pushed further by Sciamarella and Mindlin30,35 who used cell complexes to describe branched manifolds in terms of homology and torsion groups. A cell complex is a layered set of cells of different dimensions designed to classify and describe the shape of a manifold up to topological equivalence—see  Appendix A for a few definitions and examples. This specific type of cell complex is called BraMAH,36 an acronym standing for Branched Manifold Analysis through Homologies. A BraMAH complex is built in a particular fashion so that its highest-dimensional cells locally approximate subsets of points lying on a branched manifold: it does not rely on unstable periodic orbits (knots), and consequently, there is no intrinsic limitation. Recent developments of this approach have led to the definition of an unprecedented mathematical object,37 called a templex.

The key feature is that the complex is endowed with a directed graph, whose nodes are identified with the locally highest-dimensional cells of the complex. The directed graph specifies how the flow visits the structure approximated by the complex. Of course, this only has sense for a complex produced by a dynamical system. A templex gathers in a single mathematical object a set of algebraic properties, which describes the topology of the cell complex and of the flow upon it. Here, we will extend the potential of a templex in showing that a chaotic system results from a combination between oscillators and nonlinear switches obtained after a reduction procedure, which decreases the number of cells in the complex and, consequently, of nodes in the digraph. Indeed, the reduced digraph provides a systematic way to describe chaotic attractors in terms of (at least two-dimensional) oscillators (O-units) and (at least one-dimensional) switches (S-units).

The subsequent part of this paper is organized as follows. Section II A introduces the concepts underlying a templex. Section II B describes the O- and S-units and how they are related to a taxonomy of chaos. Section III describes how a complex is reduced and the way to extract the units. In Sec. IV, a set of very different attractors (the spiral and funnel Rössler attractors, the Lorenz attractor, the Burke and Shaw attractor, a non-trivial 4D attractor, and the Deng toroidal chaos) are treated. Section V gives some conclusions.

A templex is based on a cell complex constructed from a cloud of points corresponding to a trajectory produced by a dynamical system. A BraMAH complex is built using this specificity as initiated by one of us30,35,38 and later improved.36,37,39 Cell complexes provide a combinatorial approach to topology.40 A cell complex K of dimension κ ( κ N) is a layered structure built up of a finite set of k-cells ( k = 0 , 1 , , κ): 0-cells are points, 1-cells are segments, 2-cells are filled polygons, 3-cells are solid polyhedra, and so forth. The dimension of the complex K is given by the dimension of its highest-dimensional cells. In a given attractor, it is possible that, depending on the location in the state space, the highest-dimensional cells are such that k < κ. It is, therefore, useful to introduce the local dimension κ l ( x ) of the highest-dimensional cell at point x: consequently, κ = max  κ l ( x ) | x. The topological properties of such a cell complex are independent of the number, shape, or distribution of the cells: what matters is how they are glued together. These properties are expressed by its homology groups H k ( K ) and its torsion coefficients. The k-generators of the homology groups are also called k-holes since they are non-trivial chains of k-cells that encircle empty areas ( k = 1), empty volumes ( k = 2), or empty hyper-volumes ( k 3 ). The 0-holes denote the number of connected components in the cell complex. In order to make a cell complex as simple as possible, a natural question is whether there is a minimal number of k-cells to approximate the attractor underlying the trajectory.

In a dynamical system, there is a flow underlying the cloud of points used to build the cell complex, and therefore, there is a deterministic sequence for visiting the cells. This can be schemed as a directed graph, here designated as a digraph.37 When the dynamics of the flow is not regular (where regular means periodic or quasiperiodic), the flow is structured by a branched manifold as required for chaotic behavior. A branched manifold is characterized by splitting and joining loci,37 which are necessary to describe the nonequivalent paths followed by the flow. The case of regular dynamics is briefly considered in  Appendix B.

A joining locus exists in a BraMAH complex of dimension κ when there are at least three κ-cells sharing a ( κ 1)-cell. When κ = 2 as for a branched 2-manifold, the joining locus is a joining line. Let us designate as ingoing (outgoing) κ-cells, the κ-cells visited by the flow just before (after) crossing one of the joining ( κ 1)-cells belonging to a joining locus. As for strips in a template, a joining locus must have at least two ingoing and a single outgoing cell [Figs. 1(a) and 1(c)]; a splitting locus allows to have at least two κ-cells outgoing from a single ingoing κ-cell [Figs. 1(b) and 1(d)].

FIG. 1.

Three ingoing strips—the flow is from the top to the bottom—are joined into a single outgoing strip in (a). An ingoing strip is split into three outgoing strips according to the two critical points P 1 and P 2 in (b). In (a), boundaries are in solid thick lines. The thick line in (a) represents the joining locus, and the dashed thick line in (b) corresponds to the splitting locus. In (c), joining subgraph and in (d), splitting subgraph.

FIG. 1.

Three ingoing strips—the flow is from the top to the bottom—are joined into a single outgoing strip in (a). An ingoing strip is split into three outgoing strips according to the two critical points P 1 and P 2 in (b). In (a), boundaries are in solid thick lines. The thick line in (a) represents the joining locus, and the dashed thick line in (b) corresponds to the splitting locus. In (c), joining subgraph and in (d), splitting subgraph.

Close modal

To be able to compute homology and torsion groups, cell complexes must be directed or oriented; that is, each k-cell ( k > 0) must be assigned an arbitrary direction, regardless of whether there is an underlying flow. These orientations enable defining an arithmetic on the cells and, thus, an algebra of chains. In order to build a templex, there are special rules regarding how the highest dimensional cells should be oriented around the joining loci. This is performed by orienting the κ-cells that are outgoing from a joining locus according to the flow; the direction of its peripheral boundary is along the flow. The orientation is then propagated from a κ-cell to a κ-cell with an important caveat: orientation cannot be propagated across a joining locus without ambiguity. It is, therefore, sometimes necessary to restart the orientation process by repeating the procedure from a not yet oriented cell until completion of the orientation.

Each joining ( κ 1)-cell must be oriented according to the outgoing κ-cell from which it is a boundary. If there is a change in direction within a connected set of joining ( κ 1 )-cells, then there is a critical ( κ 2 )-cell, which divides this connected set into two different connected sets of joining ( κ 1 )-cells. The existence of such a critical ( κ 2 )-cell is induced by our convention for orienting the κ-cells according to the flow. This is strongly related to the convention for orienting the components of the Poincaré section from the center to the periphery to avoid different descriptions of the same dynamics.41 Once this is completed, each connected set of ( κ 1 )-cells sharing the same direction corresponds to one component of the Poincaré section. Each component is denoted by J i where i N. The concept of a multi-component Poincaré section was introduced for systems with symmetry42–44 and then generalized by Tsankov and Gilmore.45 When i = 2, there is a correspondence between a templex and a template in the dissipative three-dimensional limit, which was discussed in Ref. 37.

Once the κ-dimensional complex is built and oriented, the directed graph must encode the information regarding the way κ-cells are visited by the underlying flow. A digraph is made up of vertices or nodes, which are connected by directed boundaries, also called edges. We will use hereafter the terms nodes and edges. Two κ l-cells or nodes i and j are connected by a directed edge as i j when the flow switches from a κ l-cell γ i to a κ l-cell γ j. This naturally leads to the definition of a templex.37 

Definition 1

A templex T = ( K , G ) is made of a BraMAH complex K with dim ( K ) = κ and a digraph G = ( N , L ) whose underlying space is a branched κ-manifold associated with a dynamical system such that (i) the nodes N are κ l-cells and (ii) the edges L are the connections between the κ l-cells visited by the flow.

The advantage of a templex over a cell complex is that its algebraic descriptors include both the topology of the structure (complex) and of the organization of the flow on that structure (digraph). The topology of the structure is described in terms of homologies and torsions as in traditional topology textbooks, allowing for the distinction between, for instance, a torus and a Klein bottle. The topology of the flow on the cell complex is described through a set of definitions that enable computing the non-equivalent cell sequences or paths along the structure starting and ending in a joining locus.

Definition 2

A generatex G = ( K , C ) is a subtemplex in T = ( K , G ), where C is a cycle of the digraph G, and K K. A generatex is said to be of order θ with θ N, θ 1, if the cycle C has θ ingoing nodes. In this case, the cycle C is said to be of order θ too.

As strips in a template are always simple—in the sense that they are defined so that they do not pass through the set of joining loci more than once—we can introduce the definition of stripex as the generalization of a strip in a template.

Definition 3

Let G i be an order- θ generatex of T = ( K , G ). If θ = 1, the stripex S i = G i. If θ > 1, G i can be splitted into θ stripexes S i j with 1 j θ.

Last but not least, a stripex will be said to have a local twist if the free boundaries change their relative positions with respect to the orientation from the center to the periphery. Twisted stripexes correspond, in the dissipative 3D limit, to strips with uneven local torsions. From the point of view of an analysis in a Poincaré section, only the parity of the local torsion matters, and this intrinsic property can be directly extracted from the cell complex. If a finer characterization is required, it may be necessary to count the number of half-turns, as done with templates for three-dimensional cases.

In his first paper on chaos,16 Otto Rössler described chaos as the combination of a two-variable oscillator with a one-variable nonlinear switch. These two ingredients can be considered the fundamental units to build any type of chaos. As a fundamental unit, a periodic oscillator is characterized by a largest Lyapunov exponent, which is null. It can be furthermore either dissipative or conservative. An S-unit is a switch between two manifolds allowing stretching and squeezing, the process required for mixing trajectories and inducing the sensitivity to initial conditions.22 It is by combining an O-unit with an S-unit that Rössler got his first chaotic system,16 which was only slightly more complicated than the popular Rössler system.46 

After having proposed various systems for various types of chaos,17,47–50 Rössler proposed to construct a chaos hierarchy18,51 based on Lyapunov exponents. Such a hierarchy was then updated,52,53 before being transformed into a taxonomy of chaos,19 for defining and classifying classes of chaotic invariant sets based on shared characteristics. The taxonomic ranks, from the rough to the most refined, are as follows.19 

  • dynamical or not;

  • stochastic or deterministic;

  • positive Lyapunov exponents;

  • negative Lyapunov exponents;

  • bounding tori;

  • modality of the first-return map; and

  • a template or a templex.

As mentioned before, the numbers p of positive and q of negative Lyapunov exponents may be related to the existence of S- and O-units, respectively. This is used as guidelines for the present work, but further studies are required to predict the sign of the Lyapunov exponents from the O- and S-units.

Bounding tori are defined as permeable surfaces, which bound chaotic attractors.45 For instance, the spiral Rössler attractor is bounded by a genus-1 torus and the Lorenz attractor by a genus-3 torus. The genus g of bounding tori indicates the minimum number N c of components the Poincaré section must have: N c = g 1 for g 3 and N c = 1 for g = 1.19,45 Here, from a templex in a reduced form as we shall explain, we can work at the level of the first-return map (or the Poincaré section), that is, at the penultimate taxonomic rank, which corresponds already to a very refined level of description.

O-units will be associated with a two-dimensional oscillator (the Lotka system,54 the van der Pol equation,55 the Brusselator,56 etc.) producing a period-1 orbit. S-units—without which there is no sensitivity to initial conditions, that is, chaos—cannot exist without an O-unit: they must be bonded to an O-unit, as later described. The cycles computed from the reduced digraph will lead us to the decomposition of a reduced templex into such fundamental units, as we shall see in Sec. III.

The recipe to extract the fundamental units from a templex comprises two parts: reducing the templex and decomposing its reduced directed graph into O- and S-units. A graph decomposition is defined as follows.57 

Definition 4

A decomposition of a graph G = ( N , E ) is such that

  • G = 1 i k G i, where G i G, i.

  • E ( G i ) E ( G j ) = i , j.

We will here present how to proceed when the system is chaotic. The case of regular dynamics is discussed in  Appendix B. The starting point is a generating templex.

Definition 5

A generating templex T = ( K , G ) is a templex composed of a flow-oriented BraMAH complex associated with a digraph such that there is only one outgoing 2-cell per component of the joining locus.

Complex K is, thus, assumed to have a single cell outgoing from each component of the joining locus.37 The splitting locus in a templex T is the analog of the splitting chart in a template. It is formed by the cells that do not belong to the joining locus and from which the flow can visit at least two non-equivalent stripexes. We will also assume that a given κ l-cell cannot be bounded by a joining locus and a splitting locus so that there are at least two κ l-cells between a joining locus and a splitting locus. All the non-trivial dynamical processes will, thus, be captured between splitting and joining loci. The digraphs for our examples of joining and splitting loci are drawn in Figs. 1(c) and 1(d). We define joining and splitting nodes as follows.

Definition 6

A joining node (a splitting node) with an indegree (outdegree) 2 in a digraph G = ( N , E ) corresponds to a cell that is adjacent to a joining (splitting) locus in the complex.

Joining and splitting nodes cannot be removed during a merging process. Notice that joining and splitting nodes must be necessarily distinct; that is, a node cannot be simultaneously joining and splitting. By convention, a joining (splitting) node will be underlined (overlined). An implementation in Python of the reduction procedure detailed below can be found in https://git.cima.fcen.uba.ar/caterina.mosto/templex_properties_python.

The rationale of the reduction procedure consists of merging as many cells as possible in complex K (and its corresponding nodes in digraph G) without altering the topological properties (homology and torsion groups, joining and splitting loci, set of generatexes or stripexes). In other words, we can merge cells (nodes) as long as the operation does not create or eliminate stripexes, torsions, or holes. Once the complex is reduced ( K ¯), it remains to decompose the reduced digraph G ¯ into O- and S-units. This procedure requires identifying the set of bonds B in a templex.

Definition 7

Let T ¯ = ( K ¯ , G ¯ ) be a reduced templex with G ¯ = ( N ¯ , E ¯ ). A bond B i B is such that B i E ¯ and connects a joining node to a splitting node.

Bonds correspond to a part of the templex where the dynamics is trivial. In other words, the specificity of the attractor occurs between the splitting and joining loci, that is, in the complementary part of the bond. Bonds can be thought of as the handle(s) of a template, holding together the stripexes and marking the paths where the dynamics is irrelevant. When there are m bonds ( m N 2) belonging to a single cycle in G ¯, the bond is said to be degenerated and of order m.

Once the bonds are identified and assembled in B, we are ready to decompose the graph into units. We will start by selecting the O-units among the cycles of the reduced digraph. By construction, a cycle i in G ¯ surrounds a center of rotation c i and is related to a bond (or a set of degenerated bonds) B i.

Definition 8

The O-unit associated with a bond (or with a set of degenerate bonds) is the cycle in G ¯ (i) having the minimum number of nodes and (ii) lying closest to a center of rotation c i.

There is just one O-unit per bond (or per set of degenerated bonds). If a bond is degenerated, the associated O-unit will also be considered degenerated.

Once the O-units are identified, the S-units must be selected from the remaining cycles in G ¯ so that the set of units forms a decomposition of G ¯.

Definition 9

An S-unit is a path obtained by subtracting the edges of the O-units from the remaining cycles in G ¯.

Removing the edges corresponding to the O-units from the remaining cycles associated with the S-units exhibits that S-units cannot exist without the related O-unit. Notice that, in the absence of degenerated bonds, the total number of units (of both types) in the decomposition of G ¯ will typically coincide with the number of elements in the stripex set.

In this section, we provide the reduced templexes and their decomposition into O- and S-units for all the attractors, which were investigated in Ref. 37. The Deng toroidal chaos is also treated.

Let us now consider the Lorenz system58,
(1)
with parameter values as selected by Lorenz; that is, R = 28, σ = 10, and b = 8 3. The system is equivariant under a rotation symmetry R z ( π )43 and produces the Lorenz attractor, which is bounded by a genus-3 torus.45 

Figure 2(a) shows the generating complex K L leading to the templex T L = ( K L , G L ) for the Lorenz attractor previously provided in Ref. 37. Complex K L has a single 0-generator since the attractor has no disconnected components: H 0 [ K L ] Z. Since all the attractors here have a single connected component, we will not provide the 0th-order homology group in the subsequent part of this paper. Complex K L has two 1-generators, one around the center of rotation in each wing of the Lorenz attractor ( H 1 [ K L ] Z 2) and no enclosed cavities ( H 2 [ K L ] 0). The Lorenz attractor has a joining locus made of two components, thus corresponding to a two-component Poincaré section.

FIG. 2.

Semi-planar representation of original generating complex K L (a) from Ref. 37 with oriented components of the joining locus as thick arrows and components of the splitting locus as dashed lines; semi-planar representation of reduced cell complex K ¯ L (b); the template of the Lorenz attractor in which cells of K ¯ L are reported (c); and reduced digraph G ¯ L with bonds in color (d).

FIG. 2.

Semi-planar representation of original generating complex K L (a) from Ref. 37 with oriented components of the joining locus as thick arrows and components of the splitting locus as dashed lines; semi-planar representation of reduced cell complex K ¯ L (b); the template of the Lorenz attractor in which cells of K ¯ L are reported (c); and reduced digraph G ¯ L with bonds in color (d).

Close modal

Let us undertake the templex reduction. In K L [Fig. 2(a)], cells γ 1 and γ 8 are each bounded simultaneously by the joining locus and the splitting locus. To prevent this, cells γ 2 and γ 3 ( γ 9 and γ 10) are merged into one cell γ ¯ 3 ( γ ¯ 10). The templex can be further reduced without altering its properties: cells γ 5 and γ 7 ( γ 4 and γ 6) can be merged into one cell γ ¯ 5 ( γ ¯ 4); the same merging can be applied in the right wing, leading to the reduced complex K ¯ L [Fig. 2(b)]. The reduced templex T ¯ L has eight 2-cells instead of 14 but retains the same homologies, and the same generatex set is made of three elements, one being of order 2. The stripex set has four elements: two of them have a local twist. Merging more cells would alter the templex properties. The reduction is, thus, accomplished.

The reduced complex K ¯ L can be drawn in a semi-planar representation [Fig. 2(b)], as well as in a template-like representation. It is possible to report in a template the cells of the reduced cell complex as drawn in Fig. 2(c) so that the correspondence between the template and templex is clear. The template representation serves the purpose of displaying the genus- 3 surface of the bounding torus. Notice that the two focus holes in the bounding torus correspond to the two 1-holes exhibited in the cell complex, while the saddle hole around the rotation axis of the Lorenz symmetry in the bounding torus does not have the property required by a 1-hole as considered with homologies.

The reduced digraph G ¯ L is shown in Fig. 2(d) with the point representing the rotation axis.43 A reduced digraph is typically formed by joining nodes, splitting nodes, and a set of ‘regular” nodes, the latter ones being found in the path from a splitting node to a joining node. The eight-node reduced digraph for the Lorenz attractor has a joining locus made of two components (exhibited by the two corresponding cells, which are underlined) and a two-component splitting locus (the corresponding cells are overlined). The digraph has four regular cells. In the reduced digraph, we dressed regular nodes 5 and 12 with the superscript “ +” to indicate the presence of a positive half-twist, corresponding to the local twist in the associated stripexes.

There are two bonds: B L = { B 1 , B 2 }, where B 1 = ( 1 _ , 3 ¯ ) and B 2 = ( 8 _ , 10 ¯ ) in G ¯ L. They are highlighted in colors in Fig. 2(d).

We are now in a position to decompose G ¯ L into units. We have three cycles:
The first two cycles are of order- 1, and the third one is of order- 2. Cycles C 1 ( G ¯ L ) and C 2 ( G ¯ L ) share bond B 1; and cycles C 2 ( G ¯ L ) and C 3 ( G ¯ L ) share bond B 2, one being the symmetric of the other under the rotation symmetry.

C 1 ( G ¯ L ) associated with B 1 and C 2 ( G ¯ L ) associated with B 2 are cycles having the minimum number of nodes (three) and lying closer to the rotation centers c 1 and c 2 of each wing of the attractor. They both qualify to be O-units.

Now that the O-units have been determined, let us remove the edges already present in the O-units from the remaining cycle C 3 ( G ¯ L ). This leads to the following disjoint paths:
defining the two paths that are qualified to S-units. Digraph G ¯ L has, thus, been decomposed into four dynamical units: two O-units and two S-units.
Note that the reduced digraph is decomposed in a way, which is left invariant under the rotation symmetry. Note also that the O-units correspond to the stripexes, respectively, associated with the left and right wings of the attractor,
while the S-units are parts of the two stripexes forming the generatex of order 2,

The result is consistent with the Lorenz attractor being a C 1T 1 chaos with a structure O-S S-O where designates the transition from one wing to the other; this dynamical feature is produced by the S-units.

The Rössler attractors are produced by one of the most common chaotic systems proposed by Otto E. Rössler, read as46 
(2)
Setting b = 2 and c = 4, a spiral attractor is obtained for a = 0.432 and a funnel attractor for a = 0.492.21 A spiral Rössler attractor is characterized by p = 1 positive and q = 1 negative Lyapunov exponents: it is a C 1T 1 attractor.

Figure 3(a) shows the generating complex K R s for the spiral Rössler templex taken from Fig. 12 in Ref. 37. The corresponding reduced complex K ¯ R s is drawn in Fig. 3(b). In terms of homologies, the complex K R s has one 1-generator ( H 1 [ K R s ] Z) surrounding the saddle-focus singular point: there are no enclosed cavities ( H 2 [ K R s ] 0). It has a one-component joining locus: the thick black arrows in Figs. 3(a) and 3(b). The joining locus appears duplicated because K R s is drawn as an unfolded planar object. In planar representations of cell complexes, repeated labels indicate gluing prescriptions. There are two stripexes for the spiral Rössler templex, one of which has a local twist, corresponding to the negative half-twist in the template of Fig. 3(c).21 This is indicated with the superscript “ ” in node 6.

FIG. 3.

Generating complex K R s from Fig. 13 in Ref. 37 for the spiral Rössler attractor and its reduced version K ¯ R s (b). The splitting lines are dashed and the joining lines are thick. The corresponding reduced digraph G ¯ R s with its dynamical units is drawn in (d). The bond 1 _ 2 ¯ is colored in blue. (a) Complex K R s. (b) Reduced complex K ¯ R s. (c) Template with the cells of K ¯ R s. (d) Reduced digraph G ¯ R s with the bond as a blue arrow.

FIG. 3.

Generating complex K R s from Fig. 13 in Ref. 37 for the spiral Rössler attractor and its reduced version K ¯ R s (b). The splitting lines are dashed and the joining lines are thick. The corresponding reduced digraph G ¯ R s with its dynamical units is drawn in (d). The bond 1 _ 2 ¯ is colored in blue. (a) Complex K R s. (b) Reduced complex K ¯ R s. (c) Template with the cells of K ¯ R s. (d) Reduced digraph G ¯ R s with the bond as a blue arrow.

Close modal
The reduced digraph for the spiral Rössler attractor G ¯ R s is shown in Fig. 3(d). It has two order-1 cycles,
There is a single bond B R s = { ( 1 _ , 2 ¯ ) }, drawn as a blue arrow in Fig. 3(d). The two cycles have the same number of nodes, but the one lying closer to the center of rotation is C 1 ( G ¯ R s ), identified hence as the O-Unit. Removing the bond from the remaining cycle C 2 ( G ¯ R s ) leads to the S-unit, described by the path 2 ¯ 6 1 _.

As before, the number of units corresponds to the number of stripexes (or strips in the dissipative 3D limit), that is, to the number of monotone branches in the first-return map (not shown). The reduced digraph actually works at the same taxonomic rank as the first-return map.

When a = 0.492, the first-return map to a Poincaré section has three branches, meaning that there is a third stripex (strip) in the templex (template).21,37 From Fig. 16 in Ref. 37, a reduced complex is obtained as shown in Fig. 4(a). The cells are reported in the template as initially extracted by Letellier and co-workers21 as shown in Fig. 4(b). As expected, the complex K R f has the same homology groups as K R f: a single 1-generator around the focus-type hole of the attractor ( H 1 [ K R f ] Z) and no enclosed cavities ( H 2 [ K R f ] 0). There is a one-component joining locus in the 1-cell 0 , 1 , which appears thrice (and not twice as before) in the planar diagram, showing that the gluing prescription is more complex than for the spiral attractor. In the reduced digraph [Fig. 4(c)], the three order-1 cycles are
There is a single bond B R f = { ( 1 _ , 2 ¯ ) }. The three cycles C 1 , C 2 , C 3 have the same number of nodes, but C 1 is closest to the rotation center. It is, therefore, identified with the O-unit. The S-units are obtained by removing the bond (which is already in the O-unit) from the two remaining cycles, leading to ( 2 ¯ 5 1 _ ) and ( 2 ¯ 6 1 _ ).
FIG. 4.

Reduced templex for the funnel Rössler attractor T ¯ R f = ( K ¯ R f , G ¯ R f ) (corresponding to Fig. 15 in Ref. 37). (a) Planar diagram of reduced complex K ¯ R f. (b) Template with the cells of K ¯ R f. (c) Reduced digraph G ¯ R f with the bond as a blue arrow.

FIG. 4.

Reduced templex for the funnel Rössler attractor T ¯ R f = ( K ¯ R f , G ¯ R f ) (corresponding to Fig. 15 in Ref. 37). (a) Planar diagram of reduced complex K ¯ R f. (b) Template with the cells of K ¯ R f. (c) Reduced digraph G ¯ R f with the bond as a blue arrow.

Close modal

With respect to the spiral Rössler case, the funnel Rössler (Fig. 4) presents an additional S-unit. The two S-units are in “parallel” (as in an electric circuit); that is, they are bonded to the O-unit with the same path ( 1 _ 2 ¯ ). From Ref. 37, we know that stripex S 2 (related to the cycle C 2) has a local twist. The template strip corresponding to stripex S 3 has a negative 2 π-twist, which is indicated with superscript “ 2 ” attributed to node 6 in Fig. 4(c).

Another non-trivial example is provided by the Burke and Shaw attractor produced by the Burke and Shaw system,59 
(3)
with parameter values as a = 10 and b = 4.271. Integrating this system, a particular attractor is obtained, with a rotation symmetry R z ( π ) and whose direct template is drawn in Ref. 60.

The cell complex K ¯ BS of the Burke and Shaw attractor has a single 1-generator surrounding the focus-type hole in the attractor ( H 1 [ K BS ] Z) and no enclosed cavities ( H 2 [ K BS ] 0). There is a joining locus with the two components J 1 ( BS ) = 1 , 0 and J 2 ( BS ) = 3 , 2 , corresponding to the two components of the Poincaré section and associated with the two outgoing cells γ 1 and γ 7. As usual, joining loci in the complex are marked with thick arrows.

The reduced complex K ¯ BS for this attractor, obtained from the complex drawn in Fig. 19(b) from Ref. 37, is shown in a semi-planar representation in Fig. 5(a). The direct template representation shown in Fig. 5(b) has the two centers of rotation merged into a single one, illustrating how this attractor is actually bounded by a genus-one torus.61 The reduced digraph G ¯ BS is shown in Fig. 5(c). There are four cycles of order- 2,
Strictly speaking, only the first three cycles are needed for all the nodes to be visited once. There are two bonds, namely, B BS = { B 1 , B 2 } = { ( 1 _ , 2 ¯ ) , ( 7 _ , 8 ¯ ) }. The four cycles have the same number of nodes, but the closest to the rotation center is the first one. Cycle C 1 is, thus, an O-unit, and it is degenerated since it crosses the joining locus twice. The degeneracy of this O-unit comes from the order-2 symmetry viewed as the union of two copies of an image as discussed in Refs. 43, 44, and 62 (a complete treatment of such an approach is currently under development). From the remaining cycles, we remove all the edges already present in the O-unit, and we are left with two S-units formed by the paths
Let us recall that the four order-2 cycles lead to the following four stripexes:
where the union of the first two stripexes forms the degenerated O-unit. If the Burke and Shaw attractor appears as the dual Lorenz attractor, its structure is quite different. Note that, from the theory of bounding tori, a single-component Poincaré section is suggested. The templex approach clearly exhibits the two-components of the joining locus, and so for the Poincaré section, in agreement with early developments for computing first-return maps to this attractor.60,61
FIG. 5.

Reduced templex T ¯ BS = ( K ¯ BS , G ¯ BS ) for the Burke and Shaw attractor. Semi-planar reduced K ¯ BS is shown in (a), followed by the Burke–Shaw template representation with the cell numbers on the strips (b). The reduced digraph G ¯ BS with the colored bond, and the fundamental units are shown in (c).

FIG. 5.

Reduced templex T ¯ BS = ( K ¯ BS , G ¯ BS ) for the Burke and Shaw attractor. Semi-planar reduced K ¯ BS is shown in (a), followed by the Burke–Shaw template representation with the cell numbers on the strips (b). The reduced digraph G ¯ BS with the colored bond, and the fundamental units are shown in (c).

Close modal
Reference 37 computes the templex for a four-dimensional dynamical system, obtained as an extension of a three-dimensional one proposed by Deng.63 The system reads30 
(4)
For the parameter values reported in the caption of Fig. 6, this system produces a quite complicated attractor. Its cell complex K S has five 1-generators, and there are no enclosed cavities. This reads H 1 [ K S ] Z 5 , H 2 [ K S ] 0. The joining locus has two components.
FIG. 6.

Two views of the chaotic attractor produced by the four-dimensional system (4). Parameter values: a = 7, b = 1.45, c = 1, d = 0.5, R = 6, α = 0.3, β = 7, ϵ 1 = 0.165, ϵ 2 = 0.01, ϵ 3 = 2, and μ = 1.543.

FIG. 6.

Two views of the chaotic attractor produced by the four-dimensional system (4). Parameter values: a = 7, b = 1.45, c = 1, d = 0.5, R = 6, α = 0.3, β = 7, ϵ 1 = 0.165, ϵ 2 = 0.01, ϵ 3 = 2, and μ = 1.543.

Close modal
The digraph G S T S is drawn in Fig. 7(a). The reduced digraph G ¯ S is shown in Fig. 7(b). We have five order-1 cycles,
and a single bond
Cycle C 1 ( G ¯ S ) is the shortest (in the number of nodes) and the closest to the rotation center: it corresponds, therefore, to the O-unit. Subtracting the edges of the O-unit from the remaining cycles, we obtain the four S-units shown in Fig. 7(b). There are as many units as stripexes in T S. Local twists are indicated with a × superscript in nodes 14, 34, and 40.
FIG. 7.

Digraph G S (a) and reduced digraph G ¯ S (b) for the four-dimensional attractor produced by the system (4). The bond is colored in blue. (a) G S in templex T S from Ref. 37. (b) G ¯ S with colored bond, O- and S-units.

FIG. 7.

Digraph G S (a) and reduced digraph G ¯ S (b) for the four-dimensional attractor produced by the system (4). The bond is colored in blue. (a) G S in templex T S from Ref. 37. (b) G ¯ S with colored bond, O- and S-units.

Close modal

Let us note that cells γ 1 γ 8 and γ 40 are 3-cells. The nodes corresponding to the three-dimensional cells are highlighted in bold in G S and G ¯ S. The O-unit is, thus, made of 3-cells, approximating the solid torus or ring (see Fig. 8). The S-unit associated with C 2 is particular since it is made of four nodes due to the unusual splitting locus in the outgoing 3-cell γ 40 (see Fig. 9). We need further investigations with other systems to determine whether this is a four-dimensional signature of the dynamics. The S-units associated with C 3 and C 4 contain five nodes, and the S-unit corresponding to C 5 contains six nodes. The complexity of this attractor is clearly exhibited from these low-order cycles.

FIG. 8.

Filled torus corresponding to the O-unit of the four-dimensional chaotic attractor produced by the system (4).

FIG. 8.

Filled torus corresponding to the O-unit of the four-dimensional chaotic attractor produced by the system (4).

Close modal
FIG. 9.

Direct template for the attractor produced by the four-dimensional system (4), with the cells labeled as in reduced complex K ¯ S. Splitting loci (dashed) and joining loci (thick) are drawn.

FIG. 9.

Direct template for the attractor produced by the four-dimensional system (4), with the cells labeled as in reduced complex K ¯ S. Splitting loci (dashed) and joining loci (thick) are drawn.

Close modal

We were able to draw a template for this attractor as shown in Fig. 9: according to the Birman and Williams theorem,31 this attractor must be C 1T 1 as confirmed by the four Lyapunov exponents (kindly computed for us by N. Stankevich): Λ 1 = 0.0469, Λ 2 = 0.0, Λ 3 = 1.006, and Λ 4 = 31.62. There is a specific foliated splitting locus issued from the 3-cell γ 40 revealing the uncommon structure of this attractor. Since the thick S-unit and this splitting locus are local features, they do not lead to a second positive Lyapunov exponent. The foliated splitting locus can be identified as follows. There is a single ingoing 3-cell, which shares a 1-cell with more than one outgoing 2-cells: there is, therefore, a splitting chart without a critical point as commonly observed in three-dimensional attractors, but with a critical line.

The Deng toroidal system63 is governed by the three ordinary differential equations,
(5)
For appropriate parameter values, this system produces the unimodal toroidal chaotic Deng attractor shown in Fig. 10(a).
FIG. 10.

Chaotic toroidal attractor produced by the Deng system (5) (a) and its cell complex K D (b). The two 1-generators of H 1 [ K D ] ( g 1 and g 2) are colored in purple and magenta on the cell complex and on the regular torus, respectively. The cell complex encloses a cavity: it has a single 2-generator in H 2 [ K D ]. Parameter values: a = 3, b = 0.8, c = 1, d = 0.1, e = 0.05, η = 3.312, R = 10, α = 2.8, β = 5, ϵ = 0.1, and λ = 2. Initial conditions: x 0 = y 0 = z 0 = 0.01. (a) Trajectory. (b) Original complex K D. (c) Sketch of a torus T 2.

FIG. 10.

Chaotic toroidal attractor produced by the Deng system (5) (a) and its cell complex K D (b). The two 1-generators of H 1 [ K D ] ( g 1 and g 2) are colored in purple and magenta on the cell complex and on the regular torus, respectively. The cell complex encloses a cavity: it has a single 2-generator in H 2 [ K D ]. Parameter values: a = 3, b = 0.8, c = 1, d = 0.1, e = 0.05, η = 3.312, R = 10, α = 2.8, β = 5, ϵ = 0.1, and λ = 2. Initial conditions: x 0 = y 0 = z 0 = 0.01. (a) Trajectory. (b) Original complex K D. (c) Sketch of a torus T 2.

Close modal

The templex for this attractor T D = ( K D , G D ) is here proposed for the first time. Complex K D is composed of twenty-eight 2-cells and has the homology groups of a torus, i.e., two 1-generators H 1 [ K D ] Z 2 [drawn in Figs. 10(b) and 10(c)] and a single 2-generator H 2 [ K D ] Z defining the enclosed toroidal cavity. K D presents a joining locus [Fig. 10(b)] as required for chaos.

Digraph G D is shown in Fig. 11(a): it can be reduced to a six-node digraph G ¯ D presented in Fig. 11(b). Node 1 is kept to prevent the digraph from having a single cell between a joining locus and a splitting locus. K ¯ D has a one-component joining locus associated with node 3 in the reduced digraph G ¯ D as shown in Fig. 11(b). There are three order-1 cycles,
corresponding to the three stripexes in T ¯ D labeled S 1 ( T ¯ D ), S 2 ( T ¯ D ), and S 3 ( T ¯ D ). There is only one bond,
C 1 has the minimum number of nodes and is the closest to the rotation center, qualifying it as the only O-unit of the templex. Subtracting the edges of the O-unit from the two remaining cycles leads to the two S-units ( 15 ¯ 28 ¯ 1 3 _ ) and ( 28 ¯ 18 3 _ ) in accordance with Fig. 11(c).
FIG. 11.

Digraphs for (a) the complex K D, (b) for the reduced complex K ¯ D, and the three units, which are identified in G ¯ D. The bond is the blue arrow. (a) Digraph G D. (b) Reduced digraph G ¯ D.

FIG. 11.

Digraphs for (a) the complex K D, (b) for the reduced complex K ¯ D, and the three units, which are identified in G ¯ D. The bond is the blue arrow. (a) Digraph G D. (b) Reduced digraph G ¯ D.

Close modal

The Deng attractor is structured around a torus T 2, and consequently, it is governed by two main frequencies as confirmed by a Fourier spectrum leading to f 1 = 1.022 Hz and f 2 = 2.046 Hz, with a ratio f 1 f 2 = 0.4997 0.5. Each of these frequencies is associated with one of the two 1-generators of the homology group H 1 ( K ¯ D ); that structure is a torus as sketched in Fig. 10(c). To better understand how stripexes S 1 ( T ¯ D ), S 2 ( T ¯ D ), and S 3 ( T ¯ D ) are organized, we extract three excerpts of a chaotic trajectory visiting each of the three stripexes. The duration associated with these excerpts is τ 1 = 0.926 non dimensional units, τ 2 = 1 , 980, and τ 3 = 2.916, respectively [see Fig. 12(a)].

FIG. 12.

Unstable periodic orbits extracted from the Deng toroidal chaotic attractor (a) and the parts of the trajectory associated with each cycle (b). The joining locus is drawn as a thick arrow in (b). (a) Trajectory visiting the cells of each cycle C i. (b) Unstable periodic orbits.

FIG. 12.

Unstable periodic orbits extracted from the Deng toroidal chaotic attractor (a) and the parts of the trajectory associated with each cycle (b). The joining locus is drawn as a thick arrow in (b). (a) Trajectory visiting the cells of each cycle C i. (b) Unstable periodic orbits.

Close modal

By using a close-return technique applied to a common Poincaré section of the attractor, we find one period- 2 and two period- 3 orbits [Fig. 12(b)]. From G ¯ D, they correspond to cycles C 2, C 2 C 1, and C 3, respectively. Cycles C 2 and C 3 are of order- 1, and they should be, by definition, associated with period-1 orbits. Since cycle C 1 corresponds to the main frequency associated with a revolution around the main hole of the toroidal structure ( 1-generator g 1), the corresponding period-1 orbit should be associated with the 1-generator g 1, that is, with a trajectory in the center of the closed cavity defined by the torus: consequently, it cannot belong to the population of unstable periodic orbits embedded within the attractor. This means that the common Poincaré section is not the optimal one for investigating this toroidal chaotic flow. Using a Poincaré section based on the joining locus, the period- 3 orbit corresponding to C 2 C 1 is now a period- 2 orbit as expected for an orbit visiting two stripexes. By reducing the Poincaré section to the joining locus, the orbit periods are now equal to the number of times they visit a stripex. If we use the three symbols { 1 , 2 , 3 } to designate the three cycles C i, the extracted orbits are encoded by ( 2 ), ( 21 ), and ( 3 ), respectively. The period- 1 orbit ( 1 ) is missing, as explained above. This result is a strong confirmation that a unimodal toroidal chaos is associated with three stripexes and is, hence, described by the three-strip template in Fig. 13 as proposed in Ref. 64.

FIG. 13.

Direct template for the toroidal chaotic Deng attractor with some of the 2-cells appearing in the reduced digraph [Fig. 11(b)] (redrawn from Fig. 14 in Ref. 64). The splitting locus is dashed, and the joining locus is thick.

FIG. 13.

Direct template for the toroidal chaotic Deng attractor with some of the 2-cells appearing in the reduced digraph [Fig. 11(b)] (redrawn from Fig. 14 in Ref. 64). The splitting locus is dashed, and the joining locus is thick.

Close modal

A templex is a flow-based BraMAH complex approximating the structure on which the attractor lies, endowed with a digraph indicating how the flow visits the locally highest-dimensional cells. Reducing the templex by merging unnecessary cells and eliminating the corresponding nodes allows one to define classes of chaotic attractors at the penultimate taxonomic rank, that is, at the level of the first-return map. By considering various types of chaotic dynamics—unimodal or multimodal, with symmetry or not, three-dimensional or higher, with a toroidal underlying structure or not —, we here show that the reduced templex can be interpreted in terms of two elementary dynamical processes, namely, oscillating units (O-units) and switching units (S-units). The way these two types of units are combined can be directly extracted from the reduced templex. Our procedure has the advantage of not being limited to three-dimensional state space and not requiring the definition of a Poincaré section. As a result, splitting and joining loci are accurately identified. The toroidal chaos case exhibits the generality of these two loci types proving that the joining locus should be used to define the Poincaré section.

The number and nature of the units are derived from the properties of the stripex set, which are manipulated as cycles from the reduced digraph. When the dynamics can be synthesized by a template, the units can be recognized directly using the cell labels on the strips.

We, thus, show that it is indeed possible to decompose chaotic systems into elementary dynamical processes using a templex and a reduction procedure to algorithmically obtain the bonds between O-units and S-units. This is significant advancement for a complete taxonomy of chaos. Note that the reduced digraph does not contain all the information for ultimate discrimination between topologically non-equivalent attractors, but the templex does. The reduced digraph can be dressed with additional information in reporting the local twists of the stripexes, as well as additional (extrinsic) properties, such as the number of half-twists. The ultimate taxonomic rank is provided by the templex properties (or by the template in the dissipative three-dimensional limit). The characterization of other types of toroidal chaos with these theoretical tools is currently under study. The case of hyperchaos is also one of our next tasks. A paper devoted to the templex extraction from experimental data is currently underway to show that this method is not limited to theoretical cases. For particularly complex attractors where trajectories may remain trapped for a long time interval, the templex may vary depending on the time series being analyzed. In such cases, it would be beneficial to explore different domains of the attractor from well-selected time series to generate a templex for each part of the attractor, which can then be combined to provide a comprehensive description of the whole attractor.

This work is dedicated to Otto E. Rössler.

This work has received funding from the ANR project TeMPlex ANR-23-CE56-0002 (D.S.). G.D.C. gratefully acknowledges her postdoctoral scholarship within the ANR project. C.M. gratefully acknowledges her doctoral scholarship from CONICET. C. L. and co-workers thank N. Stankevich for the Lyapunov exponent computation.

The authors have no conflicts to disclose.

Caterina Mosto: Data curation (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Software (lead); Validation (equal); Writing – original draft (equal); Writing – review & editing (equal). Gisela D. Charó: Conceptualization (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Validation (equal); Writing – original draft (equal); Writing – review & editing (equal). Christophe Letellier: Conceptualization (equal); Formal analysis (equal); Funding acquisition (equal); Investigation (equal); Methodology (equal); Project administration (equal); Resources (equal); Validation (equal); Writing – original draft (equal); Writing – review & editing (equal). Denisse Sciamarella: Conceptualization (equal); Formal analysis (equal); Funding acquisition (equal); Investigation (equal); Methodology (equal); Project administration (equal); Resources (equal); Supervision (equal); Validation (equal); Writing – original draft (equal); Writing – review & editing (equal).

The data that support the findings of this study are available within the article as well as from the corresponding author upon reasonable request.

In algebraic topology, a k-cell within a cell complex K is a set that can be mapped continuously and invertibly to the interior of a k-disk, with its boundary divided into a finite number of lower-dimensional cells, called faces. As a whole, the cell complex K is a finite collection of cells where the faces of each cell are also elements of K, and the interiors of distinct cells do not overlap. The dimension κ of a cell complex K is defined by the dimension of its highest-dimensional cells. A cell complex is termed directed or oriented if each cell is assigned a direction. In a directed complex, a k-chain is a linear combination of k-cells, where the coefficients are integers. These chains capture the connectivity of the cells at each k level.

Let us consider the example of the Klein bottle in Fig. 14. A cell complex is provided, with four 2-cells σ i ( i = 1 , , 4), nine 1-cells ( a, b, c, d, e, f, g, h, i), and five 0-cells ( 1 , 2 , 3 , 4 , 5 ). Repeated labels are important in a planar diagram since they indicate gluing instructions. Notice that with the same set of cells and labels but changing some of their positions and directions, one can get a torus, as shown in Fig. 15.

FIG. 14.

A cell complex K with four 2-cells σ i ( i = 1 , , 4) for a Klein bottle. The repeated labels in the planar diagram indicate gluing instructions to form the Klein bottle in four dimensions.

FIG. 14.

A cell complex K with four 2-cells σ i ( i = 1 , , 4) for a Klein bottle. The repeated labels in the planar diagram indicate gluing instructions to form the Klein bottle in four dimensions.

Close modal
FIG. 15.

A cell complex K with four 2-cells σ i ( i = 1 , , 4) for a torus. The repeated labels in the planar diagram indicate gluing instructions to form the toroidal surface in three dimensions.

FIG. 15.

A cell complex K with four 2-cells σ i ( i = 1 , , 4) for a torus. The repeated labels in the planar diagram indicate gluing instructions to form the toroidal surface in three dimensions.

Close modal
Definition

Two k-chains are homologous if their difference is the border of some ( k + 1 )-chain.

In our example, c, f, and e are such that c + f and e are homologous ( c + f e) since c + f e = ( σ 2 ). The 2-cell σ 2 serves as a bridge between the edges c + f and e, implying that any information conveyed by c + f is also represented by e, and the reverse holds true as well. The homology equivalence relation leads to the definition of homology groups.

Definition

For a directed complex K, the k-dimensional homology group of K is the quotient group H k ( K ) = Z k ( K ) / B k ( K ), i.e., the group of equivalence classes of elements of Z k ( K ) with the homology relation.

The homologically independent k-dimensional cycles form the k-generators of H k, also called k-holes. The kth Betti number, denoted as β k, is the rank of H k, i.e., the number of elements in the basis. In simple terms, β 0 indicates the number of connected components in the complex, β 1 reflects the number of closed paths that encircle holes, and the number of enclosed cavities is captured by β 2. At each level, generators can be used to identify each of the connected components, holes, or cavities.

For the cell complex of the Klein bottle, there is a single 0-generator or a connected component, and therefore, β 0 ( K ) = 1, H 0 ( K ) Z with basis { 1 }. At level 1, we have H 1 ( K ) Z with basis { a }, and β 1 ( K ) = 1. Finally, at level 2, there are no enclosed cavities; i.e., H 2 ( K ) 0 and β 2 ( K ) = 0. In fact, since ( b + c d ) must be followed around twice to form a cycle, the surface is said to have a torsion coefficient of 2.

The torus has Betti numbers β 0 ( K ) = 1, β 1 ( K ) = 2, and β 2 ( K ) = 1. More precisely, H 0 ( K ) Z with basis { 1 }, H 1 ( K ) =∼ Z 2 with basis { a , b + c d } (there are two 1-holes) and H 2 ( K ) Z with basis { σ 1 + σ 2 + σ 3 + σ 4 }, enclosing the cavity. The torus has no torsion coefficients.

A torus is well represented by complex K in  Appendix A but can also be analyzed with a complex formed by a single 2-cell with at least two 1-cells, as shown in Figs. 4.13 and 4.14 in Kinsey (Ref. 40, p. 62). The fact that a single 2-cell is sufficient to build up the complex of a torus highlights the fact that a cell complex is not a simple mesh or tessellation of an object, and that each layered dimension k contains important information. When the regular surface of a torus is traveled by a quasiperiodic trajectory characterized by a templex without joining locus, the relevant information is already encoded in the cell complex since there are no distinct ways for flowing within its structure.

When the original templex has no joining locus, the underlying flow is periodic or quasiperiodic, and the manifold supporting the trajectories has no “branches.” In other words, there are no alternative paths for the flow along the underlying manifold, or equivalently, there are no stripexes. In this case, the templex can be simplified by reducing the number of highest-dimensional cells in the complex (and therefore the nodes in the digraph) as long as homologies and torsion groups do not change. The highest dimensional cells of the original complex can be merged until we get a complex with the minimal number of κ-cells, leading to the same homology and torsion groups. Thus, a single 2-cell is sufficient to describe a torus. Also, a regular flow on a torus can be represented by a reduced digraph with a self-edge representing the unique path of the flow.

An example of such a flow is provided by the Deng toroidal system63 formed by the following three ordinary differential equations:
(B1)

The symmetry of the system (B1) prevents it from developing chaos-prone folds63 so that the resulting system is quasiperiodic (Fig. 16). As expected, complex K QP has no joining locus. Its homology groups are those of a torus: K QP has two 1-generators, that is, H 1 [ K QP ] Z 2 and a 2-generator for its enclosed cavity H 2 [ K QP ] Z. The 2-cells covering the torus surface can all be merged into a single 2-cell, as shown in Fig. 16(b), without altering the homology groups. There are no torsion groups. The reduced digraph is a single node (corresponding to γ 1) with a single edge from it to itself, displaying the trivial flow organization on the torus. The digraph is shown in Fig. 17.

FIG. 16.

Original K QP and reduced K ¯ QP complexes for the quasiperiodic attractor produced by the system (B1). The reduced complex has a single 2-cell labeled γ 1.

FIG. 16.

Original K QP and reduced K ¯ QP complexes for the quasiperiodic attractor produced by the system (B1). The reduced complex has a single 2-cell labeled γ 1.

Close modal
FIG. 17.

Reduced digraph G ¯ QP for the quasiperiodic attractor produced by the system (B1) for the single 2-cell labeled γ 1.

FIG. 17.

Reduced digraph G ¯ QP for the quasiperiodic attractor produced by the system (B1) for the single 2-cell labeled γ 1.

Close modal
1.
E. N.
Lorenz
, “
A study of the predictability of a 28-variable atmospheric model
,”
Tellus
17
,
321
333
(
1965
).
2.
E. N.
Lorenz
, “On the range of atmospheric predictability,” in Proceedings of the First National Conference on Statistical Meteorology (American Meteorological Society, 1968), pp. 11–19.
3.
B.-W.
Shen
,
R. A.
Pielke
,
X.
Zeng
,
J.-J.
Baik
,
S.
Faghih-Naini
,
J.
Cui
, and
R.
Atlas
, “
Is weather chaotic?: Coexistence of chaos and order within a generalized Lorenz model
,”
Bull. Am. Meteorol. Soc.
102
,
E148
E158
(
2021
).
4.
Q.
Sun
,
T.
Miyoshi
, and
S.
Richard
, “
Control simulation experiments of extreme events with the Lorenz-96 model
,”
Nonlinear Process. Geophys.
30
,
117
128
(
2023
).
5.
E. N.
Lorenz
, “
The predictability of hydrodynamic flow
,”
Trans. N.Y. Acad. Sci.
25
,
409
432
(
1963
).
6.
J. A.
Yorke
and
E. D.
Yorke
, “Chaotic behavior and fluid dynamics,” in Hydrodynamic Instabilities and the Transition to Turbulence, edited by H. L. Swinney and J. P. Gollub (Springer, Berlin, 1985), pp. 77–95.
7.
H.
Aref
,
J. R.
Blake
,
M.
Budišić
,
S. S. S.
Cardoso
,
J. H. E.
Cartwright
,
H. J. H.
Clercx
,
K.
El Omari
,
U.
Feudel
,
R.
Golestanian
,
E.
Gouillart
,
G. F.
van Heijst
,
T. S.
Krasnopolskaya
,
Y.
Le Guer
,
R. S.
MacKay
,
V. V.
Meleshko
,
G.
Metcalfe
,
I.
Mezić
,
A. P. S.
de Moura
,
O.
Piro
,
M. F. M.
Speetjens
,
R.
Sturman
,
J.-L.
Thiffeault
, and
I.
Tuval
, “
Frontiers of chaotic advection
,”
Rev. Mod. Phys.
89
,
025007
(
2017
).
8.
S.
Sunada
,
T.
Harayama
,
K.
Arai
,
K.
Yoshimura
,
P.
Davis
,
K.
Tsuzuki
, and
A.
Uchida
, “
Chaos laser chips with delayed optical feedback using a passive ring waveguide
,”
Opt. Express
19
,
5713
5724
(
2011
).
9.
S.
Wang
,
J.
Kuang
,
J.
Li
,
Y.
Luo
,
H.
Lu
, and
G.
Hu
, “
Chaos-based secure communications in a large community
,”
Phys. Rev. E
66
,
065202
(
2002
).
10.
A.
Hastings
and
T.
Powell
, “
Chaos in a three-species food chain
,”
Ecology
72
,
896
903
(
1991
).
11.
A.
Hastings
,
C. L.
Hom
,
S.
Ellner
,
P.
Turchin
, and
H. C. J.
Godfray
, “
Chaos in ecology: Is mother nature a strange attractor?
,”
Annu. Rev. Ecol. Syst.
24
,
1
33
(
1993
).
12.
O. E.
Rössler
, “
Grundschaltungen von flüssigen Automaten und Relaxationssystemen
,”
Z. Naturforsch. B
27
,
333
343
(
1972
).
13.
C.
Letellier
and
V.
Messager
, “
Influences on Otto E. Rössler’s earliest paper on chaos
,”
Int. J. Bifurc. Chaos
20
,
3585
3616
(
2010
).
14.
M.
Sugita
, “
Functional analysis of chemical systems in vivo using a logical circuit equivalent. II. The idea of a molecular automaton
,”
J. Theor. Biol.
4
,
179
192
(
1963
).
15.
M.
Sugita
, “
Functional analysis of chemical systems in vivo using a logical circuit equivalent. III. Analysis using a digital circuit combined with an analogue computer
,”
J. Theor. Biol.
4
,
412
425
(
1963
).
16.
O. E.
Rössler
, “
Chaotic behavior in simple reaction system
,”
Z. Naturforsch. A
31
,
259
264
(
1976
).
17.
O. E.
Rössler
and
P. J.
Ortoleva
, “Strange attractors in 3-variable reaction systems,” in Theoretical Approaches to Complex Systems, Lecture Notes in Biomathematics Vol. 21 (Springer, 1978), pp. 67–73.
18.
O. E.
Rössler
, “
The chaotic hierarchy
,”
Z. Naturforsch. A
38
,
788
801
(
1983
).
19.
C.
Letellier
,
N.
Stankevich
, and
O. E.
Rössler
, “
Dynamical taxonomy: Some taxonomic ranks to systematically classify every chaotic attractor
,”
Int. J. Bifurc. Chaos
32
,
2230004
(
2022
).
20.
N. B.
Tufillaro
,
T.
Abbott
, and
J.
Reilly
,
An Experimental Approach to Nonlinear Dynamics and Chaos
(
Addison-Wesley
,
Redwood City, CA
,
1992
).
21.
C.
Letellier
,
P.
Dutertre
, and
B.
Maheu
, “
Unstable periodic orbits and templates of the Rössler system: Toward a systematic topological characterization
,”
Chaos
5
,
271
282
(
1995
).
22.
R.
Gilmore
, “
Topological analysis of chaotic dynamical systems
,”
Rev. Mod. Phys.
70
,
1455
1529
(
1998
).
23.
R.
Gilmore
and
M.
Lefranc
,
The Topology of Chaos
(
Wiley
,
2003
).
24.
H. G.
Solari
and
R.
Gilmore
, “
Relative rotation rates for driven dynamical systems
,”
Phys. Rev. A
37
,
3096
3109
(
1988
).
25.
G. B.
Mindlin
and
R.
Gilmore
, “
Topological analysis and synthesis of chaotic time series
,”
Physica D
58
,
229
242
(
1992
).
26.
M.
Lefranc
and
P.
Glorieux
, “
Topological analysis of chaotic signals from a CO 2 laser with modulated losses
,”
Int. J. Bifurc. Chaos
3
,
643
650
(
1993
).
27.
C.
Letellier
,
L. L.
Sceller
,
P.
Dutertre
,
G.
Gouesbet
,
Z.
Fei
, and
J. L.
Hudson
, “
Topological characterization and global vector field reconstruction from an experimental electrochemical system
,”
J. Phys. Chem.
99
,
7016
7027
(
1995
).
28.
R.
Gilmore
,
R.
Vilaseca
,
R.
Corbalán
, and
E.
Roldán
, “
Topological analysis of chaos in the optically pumped laser
,”
Phys. Rev. E
55
,
2479
2487
(
1997
).
29.
C.
Letellier
,
J.
Maquet
,
H.
Labro
,
L. L.
Sceller
,
G.
Gouesbet
,
F.
Argoul
, and
A.
Arnéodo
, “
Analyzing chaotic behaviour in a Belousov-Zhabotinskii reaction by using a global vector field reconstruction
,”
J. Phys. Chem. A
102
,
10265
10273
(
1998
).
30.
D.
Sciamarella
and
G. B.
Mindlin
, “
Unveiling the topological structure of chaotic flows from data
,”
Phys. Rev. E
64
,
036209
(
2001
).
31.
J.
Birman
and
R. F.
Williams
, “
Knotted periodic orbits in dynamical systems I. Lorenz’s equations
,”
Topology
22
,
47
82
(
1983
).
32.
M.
Lefranc
, “
Alternative determinism principle for topological analysis of chaos
,”
Phys. Rev. E
74
,
035202
(
2006
).
33.
M.
Lefranc
, “Reflections from the fourth dimension,” in Topology and Dynamics of Chaos, edited by C. Letellier and R. Gilmore (World Scientific Publishing, 2013), pp. 205–226.
34.
M. R.
Muldoon
,
R. S.
MacKay
,
J. P.
Huke
, and
D. S.
Broomhead
, “
Topology from time series
,”
Physica D
65
,
1
16
(
1993
).
35.
D.
Sciamarella
and
G. B.
Mindlin
, “
Topological structure of chaotic flows from human speech data
,”
Phys. Rev. Lett.
64
,
1450
1453
(
2001
).
36.
G. D.
Charó
,
G.
Artana
, and
D.
Sciamarella
, “
Topology of dynamical reconstructions from Lagrangian data
,”
Physica D
405
,
132371
(
2020
).
37.
G. D.
Charó
,
C.
Letellier
, and
D.
Sciamarella
, “
Templex: A bridge between homologies and templates for chaotic attractors
,”
Chaos
32
,
083108
(
2022
).
38.
D.
Sciamarella
, “Estructura topológica de flujos caóticos,” Ph.D. thesis (Universidad de Buenos Aires, 2001).
39.
G. D.
Charó
,
G.
Artana
, and
D.
Sciamarella
, “
Topological colouring of fluid particles unravels finite-time coherent sets
,”
J. Fluid Mech.
923
,
A17
(
2021
).
40.
L. C.
Kinsey
,
Topology of Surfaces
(
Springer-Verlag
,
New York
,
1993
).
41.
M.
Rosalie
and
C.
Letellier
, “
Systematic template extraction from chaotic attractors: I. Genus-one attractors with an inversion symmetry
,”
J. Phys. A
46
,
375101
(
2013
).
42.
C.
Letellier
,
P.
Dutertre
, and
G.
Gouesbet
, “
Characterization of the Lorenz system, taking into account the equivariance of the vector field
,”
Phys. Rev. E
49
,
3492
3495
(
1994
).
43.
C.
Letellier
and
R.
Gilmore
, “
Covering dynamical systems: Two-fold covers
,”
Phys. Rev. E
63
,
016206
(
2001
).
44.
C.
Letellier
, “Caractérisation topologique et reconstruction des attracteurs étranges,” Ph.D. thesis (University of Paris VII, Paris, France, 1994).
45.
T. D.
Tsankov
and
R.
Gilmore
, “
Topological aspects of the structure of chaotic attractors in R 3
,”
Phys. Rev. E
69
,
056206
(
2004
).
46.
O. E.
Rössler
, “
An equation for continuous chaos
,”
Phys. Lett. A
57
,
397
398
(
1976
).
47.
O. E.
Rössler
, “
Different types of chaos in two simple differential equations
,”
Z. Naturforsch. A
31
,
1664
1670
(
1976
).
48.
O. E.
Rössler
, “
Toroidal oscillation in a 3-variable abstract reaction system
,”
Z. Naturforsch. A
32
,
299
301
(
1977
).
49.
O. E.
Rössler
,
Continuous Chaos
(
Springer-Verlag
,
Berlin, Germany
,
1977
), pp.
174
183
.
50.
O.
Rössler
, “
Horseshoe-map chaos in the Lorenz equation
,”
Phys. Lett. A
60
,
392
394
(
1977
).
51.
O. E.
Rössler
, “The hierarchy of chaos,” in 10th Summer Seminar on Applied Mathematics (University of Utah, 1978), F. C. Hoppensteadt (organizer).
52.
M.
Klein
and
G.
Baier
, “Hierarchies of dynamical systems,” in A Chaotic Hierarchy (World Scientific Publishing, 1991), pp. 1–24.
53.
C.
Letellier
and
O. E.
Rössler
, “An updated hierarchy of chaos,” in Chaos: The World of Nonperiodic Oscillations (Springer, Cham, Switzerland, 2020), pp. 181–203.
54.
A. J.
Lotka
, “
Contribution to the theory of periodic reactions
,”
J. Phys. Chem.
14
,
271
274
(
1910
).
55.
B.
van der Pol
, “
On “relaxation-oscillations”
,”
Philos. Mag. VII
2
,
978
992
(
1926
).
56.
I.
Prigogine
and
R.
Lefever
, “
Symmetry breaking and instabilities in dissipative systems. II
,”
J. Chem. Phys.
48
,
1695
1700
(
1968
).
57.
J.
Kottarathil
,
S.
Naduvath
, and
J. V.
Kureethara
,
Graph Theory and Decomposition
(
CRC Press
,
2024
).
58.
E. N.
Lorenz
, “
Deterministic nonperiodic flow
,”
J. Atmos. Sci.
20
,
130
141
(
1963
).
59.
R.
Shaw
, “
Strange attractor, chaotic behavior and information flow
,”
Z. Naturforsch. A
36
,
80
112
(
1981
).
60.
C.
Letellier
,
P.
Dutertre
,
J.
Reizner
, and
G.
Gouesbet
, “
Evolution of multimodal map induced by an equivariant vector field
,”
J. Phys. A
29
,
5359
5373
(
1996
).
61.
G.
Byrne
,
R.
Gilmore
, and
C.
Letellier
, “
Distinguishing between folding and tearing mechanisms in strange attractors
,”
Phys. Rev. E
70
,
056214
(
2004
).
62.
R.
Miranda
and
E.
Stone
, “
The proto-Lorenz system
,”
Phys. Lett. A
178
,
105
113
(
1993
).
63.
B.
Deng
, “
Constructing homoclinic orbits and chaotic attractors
,”
Int. J. Bifurc. Chaos
4
,
823
841
(
1994
).
64.
S.
Mangiarotti
and
C.
Letellier
, “
Topological characterization of toroidal chaos: A branched manifold for the Deng toroidal attractor
,”
Chaos
31
,
013129
(
2021
).