Generalized connection graph method for synchronization in asymmetrical networks

https://doi.org/10.1016/j.physd.2006.09.014Get rights and content

Abstract

We present a general framework for studying global complete synchronization in networks of dynamical systems with asymmetrical connections. We extend the connection graph stability method, originally developed for symmetrically coupled networks, to the general asymmetrical case. The principal new component of the method is the transformation of the directed connection graph into an undirected graph. In our method for symmetrically coupled networks we have to choose a path between each pair of nodes. The extension of the method to asymmetrical coupling consists in symmetrizing the graph and associating a weight to each path. This weight involves the “node unbalance” of the two nodes. This quantity is defined to be the difference between the sum of connection coefficients of the outgoing edges and the sum of the connection coefficients of the incoming edges to the node. The synchronization condition for this symmetrized-and-weighted network then also guarantees synchronization in the original asymmetrical network.

Introduction

The increasing interest in synchronization in limit-cycle and chaotic dynamical systems [1], [2], [3] has led many researchers to consider the phenomenon of synchronization in large complex networks of coupled oscillators (see, e.g. [4], [5] for a sampling of this large field).

Much of this research has been inspired by technological and biological examples, including coupled synchronized lasers [6], [7], networks of computer clocks [8], and synchronized neuronal firing [9], [10]. Networks of identical or slightly non-identical oscillators often synchronize, or else form synchronous patterns that depend on the symmetry of the underlying network [11], [12].

The strongest form of synchrony in oscillator networks is complete synchronization when all oscillators do the same thing at the same time. An important problem in the study of complete synchrony is how the stability of a synchronized behavior, where the behavior could be a fixed point, a limit cycle or a chaotic attractor, is influenced by the network topology and kind of interaction. This problem was intensively studied for networks of biological oscillators [13], [14], [15], [16], and more generally of limit-cycle and chaotic oscillators[17], [18], [19], [21], [22], [23], [24], [25], [26], [27], [28], [29], [30], [31], [32], [33], [34], [35], [36], [37], [38], [39], [40], [41], [42], [43], [44].

Most methods for determining stability of synchronization in linearly coupled networks of chaotic systems are based on the calculation of the eigenvalues of the connection matrix and a term depending mainly on the dynamics of the individual oscillators (see, e.g., [17], [18], [19], [20], [21], [22], [23], [24], [25], [26], [27]). Pecora and Carroll [22] developed a general approach to the local stability of complete synchronization for any linear coupling network architecture. This approach, called the Master Stability function, is based on the calculation of the maximum Lyapunov exponent for the least stable transversal mode of the synchronous manifold and the eigenvalues of the connection matrix. This powerful method is widely used in local stability studies of synchronization in complex oscillator networks [28], [29], [30], [31], [32], [33]. Global stability results based on the calculation of the connection matrix eigenvalues were also derived for oscillator networks coupled via undirected [34], [35] and directed graphs [36]. These studies show that both local and global stabilities of complete synchronization depend on the eigenvalues of the Laplacian connection matrix.

We have previously developed an alternate way to establish synchrony which does not depend on explicit knowledge of the spectrum of the connection matrix [40]. This connection graph method combines the Lyapunov function approach with graph theoretical reasoning. It guarantees complete synchronization from arbitrary initial conditions and not just local stability of the synchronization manifold. It is also applicable to time-dependent networks. This approach was originally developed for undirected graphs and applied to global synchronization in complex networks [41], [42]. More recently, we showed that the method can be directly applied to asymmetrically coupled networks with node balance [43], [44]. Node balance means that the sum of the coupling coefficients of all edges directed to a node equals the sum of the coupling coefficients of all the edges directed outward from the node. We proved that for node balanced networks it is sufficient to symmetrize all connections by replacing a unidirectional coupling with a bidirectional coupling of half the coupling strength. The bound for global synchronization in this undirected network then holds also for the original directed network.

In this paper we extend our approach to networks with arbitrary asymmetrical connections. The connection graph of such a network is directed and the coupling coefficient from node i to node j is in general different from the coupling coefficient for the reverse direction. The new ingredient of the method is the transformation of the directed connection graph into an undirected weighted graph. This is done by symmetrizing the graph and associating a weight to each edge of the undirected graph and to each path between any two nodes. This weight involves the “node unbalance” of the two nodes. This quantity is defined to be the difference between the sum of connection coefficients of the outgoing edges and the sum of the connection coefficients of the incoming edges to the node. As in the case of node-balanced networks, the synchronization criterion derived for this symmetrical network then guarantees synchronization in the asymmetrical directed network.

The layout of this paper is as follows. First, in Section 2, we state the problem under consideration. Then, in Section 3, we derive a graph-based criterion for global synchronization in asymmetrically coupled networks and formulate the main theorem of our method. In Section 4, we show how to apply the generalized connection graph method to several examples of concrete networks. We start with the simplest network of two unidirectionally coupled oscillators, then we continue with a star-configuration and a directed network with an irregular topology. In Section 5, a brief discussion of the obtained results is given.

Section snippets

Systems under study

We consider a network of n interacting nonlinear l-dimensional dynamical systems (oscillators). We assume that the individual oscillators are all identical, even though our results can be generalized to slightly non-identical systems. The composed dynamical system is described by the n×l ordinary differential equations ẋi=F(xi)+k=1ndik(t)Pxk,i=1,,n, where xi=(xi1,,xil) is the l-vector containing the coordinates of the ith oscillator, the function F:RlRl is nonlinear and capable of

Stability system for the difference variables

Since we are interested in complete synchronization, we introduce the difference variables Xij=xjxi for any i and j. Similarly to our previous works [40], [43], we can write the stability system for the difference variables Ẋij=F(xj)F(xi)+k=1n{djkPXjkdikPXik},i,j=1,,n. The function difference F(xj)F(xi) can be rewritten in a compact vector form F(xj)F(xi)=[01DF(βxj+(1β)xi)dβ]Xij, where DF is an l×l Jacobi matrix of F. Hence, we obtain Ẋij=[01DF(βxj+(1β)xi)dβ]Xij+k=1n{djkPXjkdikPXi

Examples: Application of the method

To find an upper bound for the synchronization threshold in concrete networks, we should follow the steps of the above study.

Conclusions

We have given a sufficient condition for global complete synchronization in an arbitrary network of diffusively coupled identical dynamical systems. The condition is composed of a set of inequalities which have to be satisfied, one inequality for each edge of the connection graph. Each inequality involves a term that depends only on the individual dynamical systems, namely the coupling strength that guarantees global synchronizing of two systems. The other terms of the inequality depend only on

Acknowledgments

I.B. acknowledges the financial support of the Georgia State University Research Initiation Program (Grant FY07) and a Cariplo Foundation fellowship. V.B. acknowledges the support from the RFBR (grant No. 05-01-00509), NWO-RFBR (grant No. 047-017-018), and RFBR-MF (grant No. 05-02-19815). M.H. acknowledges the support from the SNSF (grant No. 200021-112081) and EU Commission (FP6-NEST project N 517133).

References (45)

  • N. Kopell et al.

    Math. Biosci.

    (1988)
  • D. Somers et al.

    Physica D

    (1995)
  • V.S. Afraimovich et al.

    Physica D

    (1997)
  • V.N. Belykh et al.

    Physica D

    (2004)
  • I.V. Belykh et al.

    Physica D

    (2004)
  • H. Fujisaka et al.

    Prog. Theor. Phys.

    (1983)
  • V.S. Afraimovich et al.

    Radiophys. Quantum Electron.

    (1986)
  • L.M. Pecora et al.

    Phys. Rev. Lett.

    (1990)
  • L. Fabiny et al.

    Phys. Rev. A

    (1993)
  • S.H. Strogatz

    Nature

    (2001)
  • D.L. Mills

    IEEE Trans. Communications

    (1991)
  • C.M. Gray et al.

    Proc. Natl. Acad. Sci. USA

    (1989)
  • R. Stoop et al.

    Nonlinearity

    (2000)
  • M. Golubitsky et al.

    SIAM J. Appl. Dyn. Syst.

    (2005)
  • M. Golubitsky et al.

    Bull. Amer. Math. Soc.

    (2006)
  • S.H. Strogatz et al.

    Physica D

    (1988)
  • I. Belykh et al.

    Phys. Rev. Lett.

    (2005)
  • C.W. Wu et al.

    IEEE Trans. Circ. Syst. -I: Fundam. Theory Appl.

    (1996)
  • J.F. Heagy et al.

    Phys. Rev. Lett.

    (1994)
  • L.M. Pecora et al.

    Chaos

    (1997)
  • Cited by (90)

    • Analysis of fixed-time outer synchronization for double-layered neuron-based networks with uncertain parameters and delays

      2020, Journal of the Franklin Institute
      Citation Excerpt :

      With the widespread application for synchronization, a set of theoretical results for synchronization of varieties of networks have been promoted. Most studies of synchronization come through the application of Lyapunov stableness theory, such as synchronization for asymmetrical networks [7], optimal synchronization of networks with unknown system dynamics [8], exponential synchronization of complex networks by impulsive control [9], cluster synchronization for networks with multiple networks [10,11] and synchronization of stochastic chaotic neural networks [12]. However, most of these synchronization behaviours are achieved in infinite time.

    • Robust distributed synchronization of networked linear systems with intermittent information

      2019, Automatica
      Citation Excerpt :

      The wide applicability of synchronization in science and engineering has promoted a rich set of theoretical results for a variety of classes of dynamical systems using diverse tools. The study of convergence and stability of synchronization come through the use of systems theory tools such as Lyapunov functions (Belykh, Belykh, & Hasler, 2006; Hui, Haddad, & Bhat, 2007), contraction theory (Slotine, Wang, & El-Rifai, 2004), and incremental input-to-state stability (Angeli, 2002; Cai, Zhou, & Liu, 2015). Results for asymptotic synchronization with continuous coupling between agents exist in both the continuous-time domain and the discrete-time domain; see, e.g., Scardovi and Sepulchre (2009), Moreau (2004) and Olfati-Saber, Fax, and Murray (2007), where the latter is a detailed survey of coordination and consensus for integrator dynamics, in both continuous-time and discrete-time.

    View all citing articles on Scopus
    View full text