• Abstract and Figures
• Public Full-text

Coupled Hypergraph Maps and Chaotic Cluster Synchronization Tobias B¨ohlea, Christian Kuehna,b, Raffaella Mulasc,d,e, and J¨urgen Jostc,f aFaculty of Mathematics, Technical University of Munich, Boltzmannstr. 3, 85748 Garching b. M¨unchen, Germany bComplexity Science Hub Vienna, Josefst¨adter Str. 39, 1080 Vienna, Austria cMax Planck Institute for Mathematics in the Sciences, Inselstr. 22, 04103 Leipzig, Germany dThe Alan Turing Institute, The British Library, London NW1 2DB, UK eUniversity of Southampton, University Rd, Southampton SO17 1BJ, UK f Santa Fe Institute for the Sciences of Complexity, 1399 Hyde Park Road Santa Fe, New Mexico 87501, USA Coupled map lattices (CMLs) are prototypical dynamical systems on networks/graphs. They exhibit complex patterns generated via the interplay of diffusive/Laplacian coupling and nonlinear reactions modelled by a single iterated map at each node; the maps are often taken as unimodal, e.g., logistic or tent maps. In this letter, we propose a class of higher-order coupled dynamical systems involving the hypergraph Laplacian, which we call coupled hypergraph maps (CHMs). By combining linearized (in-)stability analysis of synchronized states, hypergraph spectral theory, and numerical methods, we detect robust regions of chaotic cluster synchronization occurring in parameter space upon varying coupling strength and the main bifurcation parameter of the unimodal map. Furthermore, we find key differences between Laplacian and hypergraph Laplacian coupling and detect various other classes of periodic and quasi-periodic patterns. The results show the high complexity of coupled graph maps and indicate that they might be an excellent universal model class to understand the similarities and differences between dynamics on classical graphs and dynamics on hypergraphs. INTRODUCTION elling principle and we refer to [12] for a detailed survey of the area, which displays a high level of recent activ- ity [13–17] in the context of dynamics. In this work, There are many processes, which can be modeled el- we are interested in the effect that higher-order coupling egantly using network dynamical systems. For exam- between nodes can have on the interactions between cou- ple, this includes epidemic spreading [1], opinion for- pled maps. Coupled map lattices (CMLs) are a very mation [2], neural networks [3], synchronization [4], and classical simplified/universal class, which are quite well- game theory [5], just to name a few. The standard mod- studied on graphs. Yet, to the best of your knowledge, eling approach for network dynamics is that each ver- intertwining the dynamical properties of CMLs in the tex (or node) has one or more suitable state variables context of hypergraphs has not been considered. Here while edges (or links) provide the interaction between we propose a very general extension to chemical hyper- the vertices. Yet, in many applications, just considering graphs, which includes classical hypergraphs as a special a standard graph model is insufficient and it is necessary case. While classical hypergraphs are made by vertices to consider higher-order interactions, which can be mod- that are joined together by hyperedges, chemical hyper- elled via hypergraphs. For example, higher-order cou- graphs have the additional structure that each vertex– pling is of crucial importance in neuroscience [6], where hyperedge incidence is given either a plus sign, a minus it can encode co-firing of neurons, the joint activation sign, or both. They were introduced in [18] as a model for of several brain areas as well as more general functional chemical reaction networks, in which vertices represent relationships between neurons or brain regions. In epi- chemical elements and hyperedges represent chemical re- demiology [7], higher-order interactions can model the actions. The signs allow to distinguish reactants, prod- arXiv:2102.02272v2 [nlin.CD] 18 May 2021 different epidemiological interaction patterns for group ucts and catalysts within a chemical reaction. Moreover, dynamics. Similarly to epidemiology, contact processes in [18], also a normalized Laplace operator for chemical in opinion formation [8] clearly benefit from the inclu- hypergraphs was introduced. Such a operator is a natu- sion of group dynamics. In cell biology [9], classical graph ral generalization of the normalized Laplacian for graphs structures can also fall short to correctly represent the in- that was introduced by Fan Chung [19], and its spectrum, teractions between proteins, metabolites, or genes, where as in the graph case, encodes many qualitative properties chemical mechanisms actually dictate that certain bind- of the hypergraph to which it is associated. Conceptually, ing or catalysis patterns only form, when more than two there are two possible approaches for constructing hyper- biological components react. Historically, the need to graph Laplacians. The first approach considers a random model higher-order interactions was observed very early walk or a diffusion between the vertices of a hypergraph within ecology [10, 11], where the interaction between and identifies a Laplacian as the generator of such a pro- more than two species is very common in predator-prey cess. Thus, one looks at transition probabilities between systems due to intertwining of various competitive, mu- vertices. Since these are pairwise relations, however, they tualistic, or parasitic effects. Hence, higher-order inter- can be encoded as the weights of some ordinary graph. actions have recently become a common physical mod- The resulting Laplacian can then be represented as the 2 Laplacian of that underlying effective graph. Thus, we ber of neighbors of i. Classically, one has considered lat- would be back to graph theory. The second approach, tices G = L(d1, d2) (d1,2 ∈ N) with d1d2 nodes, or com- the one that we are adopting here, considers relations plete graphs G = Kd (d ∈ N) on d nodes; both classes al- between sets of vertices, for instance between the inputs ready display very surprising phenomena. Subsequently, and the outputs of reactions. Mathematically, one can triggered by the rise of network science, it was discovered define a boundary operator representing such relations, that new effects may arise in CMLs when the graph is not and when endowing the vertex and the hyperedge set complete or a lattice [21–23]. A commonly encountered with suitable scalar products, construct the correspond- theme in all classes of CMLs is synchronization [24]. Dy- ing Laplacian. Such a Laplacian, in contrast to the first namics is called synchronized if xn(i) = xn(j) for all approach, cannot be reduced to a graph Laplacian, but i, j and all times n ≥ n0 for some n0 ∈ N. Impor- rather genuinely represents and reflects the higher order tantly, synchronized dynamics need not be constant in relations encoded by the hypergraph. The essential point n, but could, for instance, show itself chaotic behavior. is that the operator compares the total input of an ori- In such a case, one speaks of the (complete) synchroniza- ented hyperedge to its total output, but this still allows tion of chaos [25]; we refer to [26] for the discovery of the individual contributions of the input nodes to vary, the general effect of chaotic synchronization in coupled and the same holds for the output nodes. When nodes are oscillators. An important object is the (complete) syn- contained in several hyperedges, they then are subjected chronization manifold M := {x(1) = x(2) = ··· = x(d)}, to the constraints coming from each of these hyperedges, and one is interested in the transverse stability of M. ∞ but it turns out that for many hypergraphs, there is still Consider an orbitγ ¯ = {x¯n}n=1 of the given map f. If enough freedom to make genuinely new dynamical phe- the CML is uncoupled (ǫ = 0), then the homogeneous nomena possible. solution γ = {xn(j)=¯xn} for all j is synchronized and In this work, we employ that hypergraph Laplacian to remains in M. Then one may linearize around γ, de- focus our study of the new class of coupled hypergraph rive a variational equation, and link stability for ǫ > 0 maps (CHMs) on chaotic dynamics and synchronization to the Lyapunov exponent of f and the eigenvalues of as these two phenomena are cornerstones in the study of ∆G. As shown in [21], (complete) synchronization on M more classical CMLs. Our work is structured as follows: is transversally linearly stable for (1) if First, we recall and develop the main mathematical tools µ0 we need for classical CMLs. Then we introduce CHMs |e (1 − ǫλk)| < 1 ∀k ∈{2,...,d}, (3) defined via iterated maps at each node and hypergraph Laplacian coupling between nodes. In the main dynamics where part, we study synchronization and chaos for the paradig- matic example of hyperflowers as well as some further N−1 1 ′ interesting hypergraphs. Last, we give a brief conclusion µ0 = lim log |f (¯xn)| N→∞ N and an outlook to future work. n=0 X is the Lyapunov exponent of f and λk are the nonzero COUPLED MAP LATTICES eigenvalues of the Laplacian (2), which can be ordered as To study dynamical systems induced by iterating cou- 0= λ1 < λ2 ≤···≤ λd. pled maps on networks is an established paradigm in non- linear dynamics originally termed coupled map lattices The eigenvalue 0 is simple because we assumed that G (CMLs) [20]. Each node/vertex i ∈ {1, 2,...,d} of a is connected. Obviously, it suffices to check (3) for the connected network/graph G evolves according to a time- eigenvalues λ2 and λd. Thus, in favorable cases, there discrete map f. Typical examples are the logistic map µ is a certain range of values of the coupling parameter ǫ f(x)= µx(1−x) or the tent map f(x)= 2 min{x, 1−x} for which (3) is satisfied. Note that the eigenfunctions each with parameter µ ∈ [0, 4]. The dynamics of the for the simple eigenvalue λ = 0 are the constants corre- R N 1 state xn(i) ∈ at node i at time n ∈ is defined via sponding to the tangential direction of the synchroniza- tion manifold.

Load more

  7

Copy Link