# Symmetries, Cluster Synchronization, and Isolated Desynchronization in Complex Networks

###### Abstract

Synchronization is of central importance in power distribution, telecommunication, neuronal, and biological networks. Many networks are observed to produce patterns of synchronized clusters, but it has been difficult to predict these clusters or understand the conditions under which they form, except for in the simplest of networks. In this article, we shed light on the intimate connection between network symmetry and cluster synchronization. We introduce general techniques that use network symmetries to reveal the patterns of synchronized clusters and determine the conditions under which they persist. The connection between symmetry and cluster synchronization is experimentally explored using an electro-optic network. We experimentally observe and theoretically predict a surprising phenomenon in which some clusters lose synchrony while leaving others synchronized. The results could guide the design of new power grid systems or lead to new understanding of the dynamical behavior of networks ranging from neural to social.

Synchronization in complex networks is essential to the proper functioning of a wide variety of natural and engineered systems, ranging from electric power grids to neural networks Motter2013 (). Global synchronization, in which all nodes evolve in unison, is a well-studied effect, the conditions for which are related to the network structure through the master stability function Pecora1998 (). Equally important, and perhaps more commonplace, is partial, or cluster-synchronization (CS), in which patterns or sets of synchronized elements emerge Allefeld2007 (); Ji2013 (); Zhou2006 (). Recent work on cluster synchronization has been restricted to networks where the synchronization pattern is induced either by tailoring the network geometry or by the intentional introduction of heterogeneity in the time delays or node dynamics Do2012 (); Dahms2012 (); Fu2013 (); Kanter2011 (); Rosin2013 (); Sorrentino2007 (); Williams2013 (); Belykh2008 (). These anecdotal studies illustrate the interesting types of cluster synchronization that can occur, and suggest a broader relationship between the network structure and synchronization patterns. Recent studies have begun to draw a connection between network symmetry and cluster synchronization, although all have considered simple networks where the symmetries are apparent by inspection D'Huys2008 (); Nicosia2013 (); Russo2011 (). More in-depth studies have been done involving bifurcation phenomena and synchronization in ring and point-symmetry networks GolubitskyBOOKII (); GolubitskySewartBOOK (). Here we address the more common case where the intrinsic network symmetries are neither intentionally produced nor easily discerned. We present a comprehensive treatment of cluster synchronization, which uses the tools of computational group theory to reveal the hidden symmetries of networks and predict the patterns of synchronization that can arise. We use irreducible group representations to find a block-diagonalization of the variational equations that can predict the stability of the clusters. We further establish and observe a generic symmetry-breaking bifurcation termed isolated desynchronization, in which one or more clusters lose synchrony while the remaining clusters stay synchronized. The analytical results are confirmed through experimental measurements in a spatio-temporal electro-optic network. By statistically analyzing the symmetries of several types of networks, as well as electric power distribution networks, we argue that symmetries, clusters, and isolated desynchronization are commonplace and important in many complex networks.

The general dynamical equations to describe a network of coupled identical oscillators are

(1) |

where is the -dimensional state vector of the -th oscillator, describes the dynamics of each oscillator, is a symmetric matrix of 1’s and 0’s that describes the connectivity of the network, is the coupling strength, and is the output function of each oscillator. Eq. (1) can be extended to discrete-time systems or more general coupling schemesFink2000 ().

The symmetries of the network form a (mathematical) group . Each symmetry of the group can be described by a permutation matrix that re-orders the nodes in a way that leaves the dynamical equations unchanged (i.e., each commutes with ). The set of symmetries (or automorphisms) GolubitskyBOOKII (); TinkhamBOOK () of a network can be quite large, even for small networks, but they can be calculated from using widely available discrete algebra routines Stein (); GAP4 (). Figure 1a shows three graphs generated by randomly removing 6 edges from an otherwise fully connected 11-node network. Although the graphs appear similar and exhibit no obvious symmetries, the first instance has no symmetries (other than the identity permutation), while the others have 32 and 5,760 symmetries, respectively. So for even a moderate number of nodes (11) finding the symmetries can become impossible by inspection.

Once the symmetries are identified, the nodes of the network can be partitioned into clusters by finding the “orbits” of the symmetry group: the disjoint sets of nodes that, when all of the symmetry operations are applied, permute among one another GolubitskyBOOKII (). Because Eq. 1 is essentially unchanged by the by the permutations the dynamics of the nodes in each cluster can be equal, which is exact synchronization. Hence, there are M synchronized motions , one for each cluster. In Fig. 1a, the nodes have been colored to show the clusters. For the first example, which has no symmetries, the network divides into trivial clusters with one node in each. The other instances have 5 and 3 clusters, respectively. Once the clusters are identified, Eq. (1) can be linearized about a state where synchronization is assumed among all of the nodes within each cluster. This linearized equation is the variational equation and it determines the stability of the clusters.

Equation (1) is expressed in the “node” coordinate system, where the subscripts and are identified with enumerated nodes of the network. Beyond identifying the symmetries and clusters, group theory also provides a powerful way to transform the variational equations to a new coordinate system in which the transformed coupling matrix has a block-diagonal form that matches the cluster structure. The transformation matrix is not a simple node re-ordering, nor is it an eigendecomposition of . The process for computing is non-trivial, and involves finding the irreducible representations (IRR) of the symmetry group. We call this new coordinate system the “IRR coordinate system.” A more detailed description of this process is given later in this article.

Figure 1c shows the coupling matrix in the IRR coordinate system for the three example networks. The upper-left block is an matrix that describes the dynamics within the synchronization manifold. The remaining diagonal blocks describe motion transverse to this manifold and so are associated with loss of synchronization. Thus, the diagonalization completely decouples the transverse variations from the synchronization block, and partially decouples the variations among the transverse directions. In this way the stability of the synchronized clusters can be calculated using the separate, simpler, lower dimensional ODEs of the transverse blocks to see if the non-synchronous transverse behavior decays to zero.

The general form of the transformed variational equations for clusters is,

(2) |

where we have linearized about synchronized cluster states , is the vector of variations of all nodes transformed to the IRR coordinates, and are the Jacobians of the nodes’ vector field and coupling function, respectively, and is the block diagonalization of the coupling matrix . Further details are given in a later section. We note that this analysis holds for any node dynamics, steady-state, periodic, chaotic, etc.

Figure 2a shows the optical system used to study cluster synchronization. Light from a 1550 nm light emitting diode (LED) passes through a polarizing beamsplitter (PBS) and quarter wave plate (QWP), so that it is circularly polarized when it reaches the spatial light modulator (SLM). The SLM surface imparts a programmable spatially-dependent phase shift between the polarization components of the reflected signal, which is then imaged, through the polarizer, onto an infrared cameraHagerstrom2012 (). The relationship between the phase shift applied by the SLM and the normalized intensity recorded by the camera is . The resulting image is then fed back through a computer to control the SLM.

The dynamical oscillators that form the network are realized as square patches of pixels selected from a tiling of the SLM array. Figure 2b shows an experimentally measured camera frame captured for one of the 11-node networks considered earlier in Fig. 1. The patches have been falsely colored to show the cluster structure, and the links of the network are overlaid to illustrate the connectivity. The phase shift of the -th region, , is updated iteratively according to:

(3) |

where is the self-feedback strength, and the offset is introduced to suppress the trivial solution . Eq. (3) is a discrete-time equivalent of Eq. (1). Depending on the values of , and , Eq. (3) can show constant, periodic or chaotic dynamics. There are no experimentally-imposed constraints on the adjacency matrix , which makes this system an ideal platform to explore synchronization in complex networks.

Figure 3 plots the time-averaged root-mean square (RMS) synchronization error for all four of the non-trivial clusters shown in Fig. 2b, as a function of the feedback strength . The RMS synchronization error was calculated for each cluster as where indicates an average over a time interval (here taken to be 500 iterations) and denotes a spatial average over the nodes within the cluster. In Fig. 3c-e, we plot the observed intra-cluster deviations for three specific values of indicated by the vertical lines in Fig. 3a-b, showing different degrees of partial synchronization that can occur, depending on the parameters.

Together, Fig. 3a and Figs. 3c-e illustrate two examples of a bifurcation commonly seen in experiment and simulation: isolated desynchronization, where one or more clusters lose stability, while all others remain synchronized. At (Fig. 3c), all four of the clusters synchronize. At (Fig. 3d), the magenta cluster, which contains four nodes, has split into two smaller clusters of 2 nodes each, while the other two clusters remain synchronized.

Between and , two clusters, shown in Fig. 1 as red and blue, undergo isolated desynchronization together. In Fig. 3a, the synchronization error curves for these two clusters are visually indistinguishable. The synchronization of these two clusters is intertwined: they will always either synchronize together or not at all. While it is not obvious from a visual inspection of the network that the red and blue clusters should form at all, their intertwined synchronization properties can be understood intuitively by examining the connectivity of the network. Each of the two nodes in the blue cluster is coupled to exactly one node in the red cluster. If the blue cluster is not synchronized, the red cluster cannot synchronize because its two nodes are receiving different input. The group analysis treats this automatically and yields a transverse block in Fig. 1c.

The isolated desynchronization bifurcations we observe are predicted by computation of the maximum Lyapunov exponent (MLE) of the transverse blocks of Eq. (2), shown in Fig. 3. The region of stability of each cluster is predicted by a negative MLE. While there are four clusters in this network, there are only three MLEs: the two intertwined clusters are described by a 2-dimensional block in the block-diagonalized coupling matrix . These stability calculations reveal the same bifurcations as seen in experiment.

The existence of isolated desynchronizations in the network experiments raises several questions. Since the network is connected why doesn’t the desynchronization pull other clusters out of sync? What is the relation of ID to cluster structure and network symmetry? Is ID a phenomenon that is common to many networks? We provide answers to all these questions using geometric decomposition of a group which was developed in MacArthur2009 (); MacArthur2008 (). This technique enables a finite group to be written as a direct product of subgroups where is the number of subgroups and all the elements in one subgroup commute with all the elements in any other subgroup. This means that the set of nodes permuted by one subgroup is disjoint from the set of nodes permuted by any other subgroup. Then each cluster (say, ) is permuted only by one of the subgroups (say, ), but not by any others. There can be several clusters permuted by one subgroup. This is the case of the red and blue clusters in the 32 symmetry network in Fig.1, because the associated cannot have a geometric decomposition, but may have a more structured decomposition such as a wreath product Dionne1996 ().

One can show (see the Supplementary Information) that the above decomposition guarantees that the nodes associated with different subgroups all receive the same total input from the other subgroups’ nodes. Hence, nodes of each cluster do not see the effects of individual behavior of the other clusters associated with different subgroups. This enables the clusters to have the same synchronized dynamics even when another cluster desynchronizes. If that state is stable we have ID.

How common is such an ID situation we outlined above? We have examined statistics for some classes of random and semi-random graph types that suggest that when symmetries are present the opportunity for ID dynamics will be common although the stability for such will depend on the dynamical systems of the network nodes.

We examined 10,000 realizations of three random and semirandom networks: (1) randomly connected nodes (random graphs) similar to Erdos-Renyi graphs BollobasBOOK (), (2) scalefree tree graphs following Barabasi and Albert Barabasi1999 (); Albert2002 (), and (3) scalefree using the construction to give predetermined degree distributions characterized by exponent with same number of connections as in (1) goh2001universal (). In the Supplementary Information we detail how we generated the realizations and tested for duplicates and statistical relevancy.

Figure 4a shows the cumulative distribution of symmetries for each type of network. The scalefree graphs generally have fewer symmetries than the other two types. The Barabasi and Albert scalefree graphs often have many orders of magnitude more symmetries than the others which is a result of their hub and tree structure. All have similar distributions overall, but on different scales of symmetries. Cases of no symmetry are extremely rare for all graphs in these parameter ranges.

As shown in Fig. 4b,c, almost all graphs for each type have several nontrivial clusters and more than 1 subgroup with the Barabasi and Albert distribution skewing toward somewhat larger numbers. The median numbers of clusters for the random, Barabasi and Albert, and networks are 3, 5, and 3, respectively. The median numbers of subgroups are 3, 4, and 3, respectively. The percents of cases where the number of subgroups is less than the number of clusters (intertwined cases) are 33%, 59%, and 33%, respectively. Thus, the scenario is present for almost all of these networks to experience ID.

Finally, we examined two existing networks: the Nepal power grid NepalElectReport2011 () and the Mesa del Sol electrical grid abdollahy2012pnm (). We show the Nepal grid since its small size is easier to display in Fig. 5. Also shown is the block diagonalization of the coupling matrix. Here we treat the grid analogous to Motter2013 () in which all power stations are identical with the same bidirectional coupling along each edge. This man-made network has 86,400 symmetries, three nontrivial clusters (plus two trivial ones), and three subgroups (one for each nontrivial cluster). This implies it is possible for this network to split into three sets of synchronized clusters and one of those could lose stability while the others remain synchronized which is ID. The Mesa del Sol grid has 4096 symmetries, 20 nontrivial clusters, and 10 subgroups. The network has three intertwined clusters, two with 4 clusters and one with 5 clusters, making ID a possibility.

Many other networks were studied for symmetries in MacArthur2008 (); MacArthur2009 () for the purpose of finding motifs and redundancies, but not dynamics. Those networks were Human B Cell Genetic Interactions, C. Elegans genetic interactions, BioGRID data sets (Human, S. cerevisiae Drosphila, and Mus musculus), the internet (Autonomous Systems Level), and the US Power Grid. All the networks had many symmetries ranging in number from on the order of to , and could be decomposed into many subgroups (from 3 to more than 50). The subgroups were 90% or more made up of basic factors (not intertwined) consisting of various orders of the symmetric group . Hence, viewed as dynamical networks, all could show ID in the right situations.

The phenomena of symmetry-induced cluster synchronization and ID appear to be possible in many model, man-made, and natural networks. We’ve show that ID is explained generally as a manifestation of clusters and subgroup decompositions. Furthermore, computational group theory can greatly aid in identifying cluster synchronization in complex networks where symmetries are not obvious or far too numerous for visual identification. It also enables explanation of types of desynchronization patterns, and transformation of dynamic equations into more tractable forms. This leads to an encompassing of or overlap with other phenomena which are usually presented as separate. This list includes (1) remote synchronization Nicosia2013 () in which nodes not directly connected by edges can synchronize (this is just a version of cluster synchronization), (2) some types of chimera states Abrams2004 (); Hagerstrom2012 () which can appear when the number of trivial clusters is large and the number of nontrivial clusters is small, but the clusters are big (see Laing2009 () for some simple examples), (3) partial synchronization where only part of the network is synchronized (shown for some special cases in Belykh2000 ()). We note that although we have concentrated mostly on the maximal symmetry case, we can also examine the cases of lower symmetry induced by bifurcations that break the original symmetry and the same group theory techniques will apply to those cases. Some of this is developed for simple situations (rings or simple networks with point group symmetry) in Ref. GolubitskyBOOKII (), but we now have the ability to extend this to arbitrary complex networks. Finally, we note that it is possible to extend this approach to systems with nonidentical oscillators and weighted and/or directional coupling or to hypernetworks Irving2012 ().

## Symmetries, Synchronization Clusters, and Block Diagonalization of the Variational Equations.

Here we outline the steps necessary to determine the symmetries of the network, obtain the clusters, find the irreducible representations (IRRs), and the most crucial part, calculate the transformation from the node coordinates to the IRR coordinates that will block-diagonalize , since commutes with all symmetries of the group SaganBOOK ().

Using the discrete algebra software it is straightforward to,

(1) Determine the group of symmetries of .

(2) Extract the orbits which give the nodes in each cluster and extract the permutation matrices

(3) Using the character table of the group and the traces of the ’s determine which IRRs are present in the node-space representation of the group. Remark: This step is discussed in any book on representations of finite groups (e.g. Ref. TinkhamBOOK ())

(4) Put each into its appropriate conjugacy class.

The next steps are to generate the transformaion from the group information and they require writing code on top of the discrete algebra software.

(5) For each IRR present construct the projection operator TinkhamBOOK () from the node coordinates onto the subspace of that IRR, where indexes the set of IRRs present. Thus,

(4) |

where is a conjugacy class, is the character of that class for the th IRR, is the dimension of the th IRR and is the order (size) of the group. Remark: The trivial representation (all IRR matrices=1 and ) is always present and is associated with the synchronization manifold. All other IRRs are associated with transverse directions.

(6) Use singular value decomposition on to find the basis for the projection subspace for the th IRR.

(7) Construct by stacking the row basis vectors of all the IRRs which will form an matrix.

Once we have we can transform the variational equations as follows. Let be the set of nodes in the th cluster with synchronous motion . Then the original variational equations about the synchronized solutions are (in vectorial form),

(5) |

where the -dimensional vector and is an -dimensional diagonal matrix such that

(6) |

. Note that , where is the -dimensional identity matrix.

Applying to Eq. (5) we arrive at the variational matrix equation shown in Eq. (2),where , is the transformed , and is the block diagonalization of the coupling matrix . We can write the block diagonal as a direct sum , where is a (generally complex) matrix with the multiplicity of the th IRR in the permutation representation , the number of IRRs present, and the dimension of the th IRR, so that Vallintin2009 (); GoodmanBOOK (). For many transverse blocks is a scalar, i.e. . However, the trivial representation which is associated with the motion in the synchronization manifold has . The form of the variational equation for the first examples is shown in Fig. 1c. Each block in Fig. 1c is governed by a separate variational ODE as given in Eq. (2). Note that the vector field can contain a self-feedback term as in the experiment and other feedbacks are possible, e.g. row sums of , as long as those terms commute with the .

Acknowledgements We acknowledge help and guidance with computational group theory from Prof. David Joyner of the US Naval Academy and information and computer code from Ben D. MacArthur and Rubén J. Sánchez-García both of the University of Southampton to address the problem of finding group decompositions. Thanks to Vivek Bhandari for providing us a copy of the Nepal Electricity Authority Annual Report NepalElectReport2011 (). Thanks to Shahin Abdollahy and Andrea Mammoli for providing us data on the Mesa Del Sol electric network abdollahy2012pnm ().

## References

- (1) A. E. Motter, S. A. Myers, M. Anghel, and T. Nishikawa, “Spontaneous synchrony in power-grid networks,” Nat. Phys. 9, 191–197 (2013)
- (2) L.M. Pecora and T.L. Carroll, “Master stability functions for synchronized coupled systems,” Physical Review Letters 80, 2109–2112 (1998)
- (3) C. Allefeld, M Müller, and J. Kurths, “Eigenvalue decomposition as a generalized synchronizatioun cluster analysis,” Int. J. Bif. Chaos 17, 3493–3497 (2007)
- (4) P. Ji, T.K. Peron, P.J. Menck, F.A. Rodrigues, and J. Kurths, “Cluster explosive synchronization in complex networks,” Physical Review Letters 110, 218701 (2013)
- (5) C. Zhou and J. Kurths, “Hierarchical synchronization in complex networks with heterogeneous degrees,” CHAOS 16, 015104 (2006)
- (6) A-L.Do, J. Hoefener, and T. Gross, “Engineering mesoscale structures with distinct dynamical implications,” New Journal of Physics 14, 115022 (2012)
- (7) T. Dahms, J. Lehnert, and E. Schöll, “Cluster and group synchronization in delay-coupled networks,” Physical Review E 86, 016202 (2012)
- (8) C. Fu, Z. Deng, L. Huang, and X. Wang, “Topological control of synchronous patterns in systems of networked chaotic oscillators,” Physical Review E 87, 032909 (2013)
- (9) I. Kanter, M. Zigzag, A. Englert, F. Geissler, and W. Kinzel, “Synchronization of unidirectional time delay chaotic networks and the greatest common divisor,” EPL 93, 6003–1–6 (2011)
- (10) D. P. Rosin, D. Rontani, D.J. Gauthier, and E. Schöll, “Control of synchronization patterns in neural-like boolean networks,” Physical Review Letters 110, 104102 (2013)
- (11) F. Sorrentino and E. Ott, “Network synchronization of groups,” physical Review E 76, 056114 (2007)
- (12) C. Williams, T. Murphy, R. Roy, F. Sorrentino, T. Dahms, and E. Schöll, “Experimental observations of group synchrony in a system of chaotic optoelectronic oscillators,” Physical Review Letters 110, 064104 (2013)
- (13) V. Belykh, G.V. Osipov, V.S. Petrov, J.K.A. Suykens, and J. Vandewalle, “Cluster synchronization in osccillatory networks,” CHAOS 18, 037106 (2008)
- (14) O. D’Huys, R. Vicente, T. Erneux, J. Danckaert, and I. Fischer, “Synchronization properties of network motifs: Influence of coupling delay and symmetry,” CHAOS 18, 037116 (2008)
- (15) V. Nicosia, M. Valencia, M. Chavez, A. Daz-Guilera, and V. Latora, “Remote synchronization reveals network symmetries and functional modules,” Physical Review Letters 110, 174102 (2013)
- (16) G. Russo and J-J. E. Slotine, “Symmetries, stability, and control in nonlinear systems and networks,” Physical Review E 84, 041929 (2011)
- (17) M. Golubitsky, I. Stewart, and D.G. Schaeffer, Singularities and groups in bifurcation theory, Vol. II (Springer-Verlag, New York, NY, 1985)
- (18) M. Golubitsky and I. Stewart, The Symmetry Perspective: From Equilibrium to Chaos in Phase Space and Physical Space, Vol. II (Berkhäuser-Verlag, Basel, 2002)
- (19) K. Fink, G. Johnson, D. Mar, T. Carroll, and L. Pecora, “Three-oscillator systems as universal probes of coupled oscillator stability,” Physical Review E 61, 5080–5090 (2000)
- (20) M. Tinkham, Group Theory and Quantum Mechanics (McGraw-Hill, New York, NY, 1964) the group analysis can be understood with only a basic knowledge of finite groups and group representation theory on the level of the first three chapters
- (21) William Stein, SAGE: Software for Algebra and Geometry Experimentation (http://www.sagemath.org/sage/, http://sage.scipy.org/, 2013)
- (22) The GAP Group, GAP: Groups, Algorithms, and Programming, Version 4.4 (http://www.gap-system.org, 2005)
- (23) A.M. Hagerstrom, T.E. Murphy, R. Roy, P. Hovel, I. Omelchenko, and E. Schöll, “Experimental observation of chimeras in coupled-map lattices,” Nature Physics 8, 658–61 (2012)
- (24) B.D. MacArthur and R.J. Sanchez-Garcia, “Spectral characteristics of network redundancy,” physical Review E 80, 026117 (2009)
- (25) B.D. MacArthur, R.J. Sanchez-Garcia, and J.W. Anderson, “On automorphism groups of networks,” Discrete Appl. Math. 156, 3525 (2008)
- (26) B. Dionne, M. Golubitsky, and I. Stewart, “Coupled cells with internal symmetry: I. wreath products,” Nonlinearity 9, 559Ð574 (1996)
- (27) B. Bollobas and O.M. Riordan, Handbook of Graphs and Networks (Wiley-VCH, Weinheim, 2003)
- (28) A.-L. Barabasi and R. Albert, “Emergence of scaling in random networks,” Science 286, 509–512 (1999)
- (29) R. Albert and A.-L. Barabasi, ‘‘Statistical mechanics of complex networks,” Reviews of Modern Physics 74, 48–96 (2002)
- (30) K-I Goh, B Kahng, and D Kim, “Universal behavior of load distribution in scale-free networks,” Physical Review Letters 87, 278701 (2001)
- (31) “Nepal electricity authority annual report 2011,” Available at www.nea.org.np
- (32) S Abdollahy, O Lavrova, A Mammoli, S Willard, and B Arellano, “PNM smart grid demonstration project from modeling to demonstration,” in Innovative Smart Grid Technologies (ISGT), 2012 IEEE PES (IEEE, 2012) pp. 1–6
- (33) D.M. Abrams and S.H. Strogatz, “Chimera states for coupled oscillators,” Physical Review Letters 93, 174102–1–4 (2004)
- (34) C.R. Laing, “The dynamics of chimera states in heterogeneous kuramoto networks,” Physica D 238, 1569–88 (2009)
- (35) V.N. Belykh, I.V. Belykh, and M. Hasler, “Hierarchy and stability of partially synchronous oscillations of diffusively coupled dynamical systems,” Physical Review E 62, 6332–45 (2000)
- (36) D. Irving and F. Sorrentino, “Synchronization of dynamical hypernetworks: dimensionality reduction through simultaneous block-diagonalization of matrices,” Physical Review E 86, 056102 (2012)
- (37) B.E. Sagan, The Symmetric Group (Wadsworth Brooks, Pacific Grove, CA, 1991)
- (38) F. Vallentin, “Symmetry in semidefinite programs,” Linear Algebra And Its Application 430, 360–369 (2009)
- (39) R. Goodman and N.R. Wallach, Representations and Invariants of the Classical Groups (Cambridge University Press, 1998)
- (40) Duncan J. Watts and Steven H. Strogatz, ‘‘Collective dynamics of ’smallworld’ networks,” Nature 393, 440–442 (1998)

## Supplementary Information

Adjacency, Transformation, and Block-Diagonal Matrices and orbits. Here are the adjacency matrices, the clusters (group orbits), the transformation, and block-diagonalized coupling matrices for Fig. 1 in more detail.

The 0-symmetry case:

There are no nontrivial clusters, the transformation matrix is just the identity matrix , and .

The 32-symmetry case:

The nontrivial clusters are the nodes [1, 7], [2, 3, 7, 9], [4, 6], [5, 10] (the numbering of nodes matches the row and column numbers of ). The transformation matrix is,

And the block diagonal coupling matrix is,

The 5670-symmetry case:

The nontrivial clusters are the nodes [2, 6, 1, 7, 8, 11], [3, 5, 9, 10]. The transformation matrix is,

And the block diagonal coupling matrix is,

Subgroup decomposition and cluster dynamics. To start let , a subgroup of , permute only cluster and be the permutation on the indices of nodes in for one permutation . Assume is not in so it is not permuted by and recall that commutes with all permutations in , then we have (just concentrating on the terms from ),

where is, in general, another node in and the sums over other clusters are unchanged. This shows that all nodes in are coupled into the th node in the same way (the same factor). Similarly, if we use a permution on the cluster containing we can show that all the nodes of are coupled in the same way to the nodes in . Hence, nodes of each receive the same input sum from the nodes of whether the nodes of are synchronized or not. This explains how the cluster can become desynchronized, but the nodes of can still be synchronized – they all have the same input despite the desynchronization, thus making the synchronous state flow invariant. If it is also stable, this is the case of ID. This argument is easily generalized to the case when permutes nodes of several clusters as this will just add other similar sums to Eq. (Supplementary Information). The latter case explains the intertwined desynchronization in the experiment and is a more general form of ID.

Statistics of random graphs. Random graphs were generated by starting with 25 nodes completely connected and randomly deleting 20 edges. Scalefree Barabasi and Albert graphs were based on the original Barabasi and Albert preference algorithm Barabasi1999 () using the SAGE routine RandomBarabasiAlbert. These had 25 nodes with 24 edges and a tree structure. Scalefree graphs with a specific power-law distribution were generated according to goh2001universal () using and . 10,000 realizations of each graph type were generated. We tested several 10,000 realizations and we see very little variation in statistics between realizations of the same class leading us to believe that we are sampling fairly and enough to trust our results. We also checked for equivalent (isomorphic) graphs to see how much repetition we had. The random systems yielded on average 1 equivalent pair per 10000 realizations. The scalefree cases yielded about 5 to 10% equivalent graphs. Apparently we are not near the maximum number of inequivalent graphs for any of the classes although the results suggest that the scalefree classes are much smaller than the random class. Even with just 100 realizations the main trends in number of symmetries and other statistics are evident although such small samples occasionally miss those symmetry cases that are not too common in the class.

The scalefree model. The model generates a scale free network with nodes and edges and a specified power law degree distribution exponent . Start with vertices, assign to each vertex , a weight , where the exponent lies in the range [0, 1). Assume that initially no edges are present among the network vertices, then edges are added one by one until connections are created. For each new edge, two vertices are randomly selected, each one with probability proportional to its weight, and they are connected unless a link already exists or the two selected nodes are the same. By following this procedure, the resulting network is scale free and the power law degree distribution , with exponent .

SmallWorld networks We also studied symmetries, clusters, and subgroup decompositions in smallworld graphs. Smallworld graphs Watts1998 (); BollobasBOOK () were generated by starting with a ring of nearest neighbor connected nodes, then adding a fixed number of edges to give the same number of edges as the random graphs in the text. We found we had to add many edges beyond the usual few used to generate the smallworld effect because adding only a few edges beyond the ring rarely resulted in any symmetries. As a result the smallworld examples approached being a like the random graphs so we do not display their results although the two systems each have symmetries that the other does not so they appear to not be exactly identical.

The Mesa del Sol power grid. In Supplementary Fig. 6 we show a circle plot of the Mesa del Sol network which, because of the network size (132 nodes), exposes the cluster structure much better. There are 20 nontrivial clusters and 10 subgroups in the decomposition. There are a large number of trivial clusters with only about 1/3 of the nodes being in a synchronizable cluster. However, in that subset of clustered nodes the subgroup decomposition shows that ID is dynamically possible.