1. Introduction
A complex network is a set of coupled elements with features that do not occur in simple networks such as lattices or random graphs, but often occur when trying to model real-world systems. Many human-made and natural systems are described by models of complex networks, such as the citation network of scientific papers, electrical networks, wireless communication networks, food chains, social networks, ecosystems, and so on [
1,
2,
3,
4,
5,
6].
The concept of complex networks introduced by Watts and Strogatz (WS) [
7] aims to describe the transition from a regular lattice to a random graph. These WS networks exhibit several distinctive traits such as a high degree of clustering and a small average path between any two nodes. Moreover, empirical results by Barabasi and Albert [
8] show that large WS networks display scale-free properties, such as the Internet, social networks, and some metabolic networks, to cite some examples. To understand the dynamical behaviors of the complex networks, we may extend the previous models by adding dynamic elements as nodes in the network. One of the most studied phenomena in this kind of network is synchronization, where the networked elements present the same behavior. These dynamical complex networks display various synchronization phenomena between the nodes [
9,
10,
11,
12,
13]. These synchronization phenomena have attracted increasing attention from various fields of science and engineering.
As an example, Wang and Chen introduced a simple uniform scale-free network model which reflects the complexity of the network structure but still, it refers only to a simple non-varying dynamical network—the coupling strength remains constant for all connections in the network [
9]. In the real world, most of the networks we found have different coupling and or varying coupling strengths for distinct nodes. Moreover, many networks we find in nature have varying coupling strengths and still exhibit an interesting synchronization phenomenon [
14]. For instance, Sorrentino and Ott [
15] used an adaptive strategy to guarantee synchronization in a network with non autonomous couplings that evolves under the effects of exogenous (unpredictable) factors. In Reference [
16] the authors deals with the problem of synchronization of systems in networks with fix and switching externally defined couplings and communication time delays, for the case of fluctuating structures, it changes or discretely switches the couplings, this is the network couplings changes and remains constant or fixed for a certain time, in our proposal the varying couplings are not constant.
Nevertheless, we observe that many of the existing research about synchronization has been devoted to completely regular networks; this is, networks where the structure and coupling strength remains fixed. Now, we consider possible changes in the links between nodes. In such networks, the changes can be due to the collective behavior (synchronization), or due to the dynamics in the interconnected nodes, or changes in the links between nodes, or a mixture of all of them.
In this note, the synchronization phenomenon in weighted networks is presented, this is, networks where coupling strengths connect the nodes and have different values, moreover, these values are not fixed. The dynamics of these coupling strengths are divided into two; those who depend explicitly on time, and those who depend on the dynamics of the interconnected nodes. Therefore, we claim that the first case is a time-varying network and the latter is an evolving network.
These changes in couplings can be identified as changes in time and changes due to the behavior of the individuals. Changes along time refer to changes in the couplings between nodes which are defined as function purely of time. In this sense, the changes in the network are dictated by an external agent, and the nodes in the networks do not influence their couplings, for instance in Reference [
17] the authors consider varying couplings but a control action is implemented to achieve synchronization in a defined reference. Another example is in brain deceases or disorders such as Parkinson, Alzheimer’s, and so forth, where external deep brain stimulation is provided to modify the behavior of certain parts of the brain [
18,
19,
20]. On the other hand, for the evolving couplings, we depart from an analogy with biological systems, where evolution is not only a change along time, but the individuals define it according to their necessity. That is, the couplings in a network can change obeying specific evolving rules defined by the interacting agents. In this sense, the structure of the network is changing along time but defined by the collectivity [
21].
Both cases imply possible changes in the network structure. However, the nature of those changes is an important issue. In the first case (non autonomous couplings) the dynamics of the nodes do not determine or define the coupling values and can be considered as an exogenous imposition in the network structure. On the other hand, when the couplings are defined or influenced by the dynamics of the connected nodes, the changes in the network structure are defined collectively, and it is then considered as “evolution” since the behavior of the node produces the network changes. An important issue behind of these cases is the synchronization, on the one hand, the exponential stability of the synchronization is reached since the couplings are imposed and in such a case the stability is controlled. On the other hand, in the case of evolutionary couplings the stability is not controlled and depends absolutely on the coupling behavior and it cannot be controlled.
For dynamical non autonomous couplings in the network, some important assumptions are considered, the functions that describe the coupling strength are chaotic. We also assume that these chaotic functions are positive definite since changes from positive to negative would disconnect the node for an unknown period (when the coupling is 0), and that would completely change the network size.
For the case of networks with evolutionary couplings, a function of the nodes states determines the coupling. It is assumed that the function that describes the coupling is the same for any set of two connected nodes, and this function is positive definite.
Research on specific networks that fit this case has been of great interest [
22]. Works regarding analysis techniques for human epileptic brain networks [
23] and the study of the behavior of dynamical networks that evolve slowly as a function of their system properties [
24] are examples of the interest in this phenomena.
Our contribution addresses the issue from a more general approach, trying to characterize the generality of these evolving networks and their ability to achieve synchronization.
The paper is organized as follows, in the next section the model for a dynamical network is presented, in
Section 3, models for networks with chaotic couplings and evolving couplings are presented, alongside numerical examples to corroborate our results, and in
Section 4 some conclusions are provided.
2. Synchronization of Dynamical Networks
Consider a typical dynamical network consisting of N identical linear and diffusively coupled nodes, with each node being an n-dimensional dynamical system. The state equations for the
i-th node in the network is described by
In this model
are the state variables of node
i,
represents the coupling strength,
is a continuous vector field which describes the dynamical behavior of the node and
is the outer coupling matrix linking coupled variables. The coupling matrix
A is composed by
, such that, if there is a connection between node
i and node
j then
; otherwise
and
. Suppose that matrix A satisfies
L is the Laplacian matrix of the network that satisfies the dissipative coupling condition given by
, therefore,
is an eigenvalue of
L associated to the eigenvector
the eigenvalues of matrix
L have all real part less than or equal to 0 (
) and if the Laplacian
L is irreducible, then the eigenvalue
has multiplicity 1.
It is said that the nodes in a network synchronize if and since by the dissipative condition of the Laplacian matrix, the synchronization is achieved in an isolated node and the solution can be an equilibrium point, a periodic orbit or a chaotic attractor. Therefore, the stability of the network synchronization depends on the isolated node dynamics, the coupling matrix given by L, the coupling strength and the inner coupling matrix. The following Theorem is a well known result which provides conditions for the synchronization.
Theorem 1 ([
9,
10]).
Let us consider the network Equation (1) and let the eigenvalues of L areNow suppose that there exists an diagonal matrix and there are two constants and such thatFor all , with the Jacobean matrix at the isolated solution and moreover ifthen the synchronization state is exponentially stable. Mainly these are the conditions from where the synchronization of complex networks departs [
9,
10,
25]. Note that the Laplacian matrix has constant elements, this implies that the connection between nodes does not changes along time which means that there is not any possibility of changing the connection and the strength between nodes. At this point, two possible changes in the connection and the strength are identified. On the one hand, a non autonomous coupling, which only depends on time; this idea in the network represents an arbitrary change where the node dynamic does not interfere, this is, the node dynamics in the network does not determine the network structure, in this case, it is called non autonomous (time-varying) coupling. On the second hand, the couplings depend on the dynamics of the nodes, for this case, the network structure is defined by the node dynamics along the time, therefore it can be seen as an evolution of the network structure, since the nodes define their coupling according to the network requirements, in such a case it is called evolutionary couplings.
3. Synchronization of Non Autonomous (Time-Varying) and Autonomous Evolutionary Couplings Network
In this section, the two scenarios previously described are presented. The first one considers a connection that varies chaotically, this is, a linking between nodes which is non autonomous, however, this variation is dictated by an external subsystem, recall that in such a case the node dynamics in the network do not define the connection. In the second scenario, the connections are considered as functions of the states of the connected nodes, that is, an evolutionary coupling, in this case, the nodes define the strength between nodes and have the possibility of defines the network structure.
To represent the two couplings, from (
1)
c and
in the coupling matrix
A are given by
, such that the connection between any nodes vary along time as a function of an external subsystem or evolve according to the node dynamics. Actually, the couplings can be a function or depend on both, an external subsystem and as a function of the nodes in the network, therefore the model of the network is given by
where
are the state variables for the node
i;
is the vector field of the external subsystem. Let
be the matrix couplings in which
is defined as follows: if there is a connection between nodes
i and
j at time
t, then
; otherwise
for
. The diagonal elements in
are defined similar to the diagonal elements on matrix
L in (
1) given by
Let us consider that functions in matrix
are continuous, bounded, and positive definite since a zero-crossing change the network size (network growth or shrinking). Now the main objective is to achieve synchronization of the nodes in the presence of non autonomous or evolving couplings. Thus, from the model (
5) and following Theorem 1, we propose the following
Proposition 1. The synchronization manifold of the network (5) with varying/evolutionary couplings is exponentially stable if there exists a constant value such that . Proof. Let us first consider that
in (
4), this means that the coupling strength is now in
. Then now let us assume that there exists a constant value
such that
, from the properties of the Laplacian matrix
, which is dissipative for all
t, it has eigenvalues
and
for
and
thus, without loss generality
for all
t and from Theorem 1, we can consider
which satisfies (
3) and then the synchronization manifold is exponentially stable. □
Note that the value is negative for all and this value depends on the weight of the connection between nodes, this can be seen as if there is a particular coupling strength between two nodes and it has a particular weight, then the synchronization relay on the weight of the coupling and it moves the values of the bounds of .
3.1. Network Synchronization with Time Varying Couplings
We seek complete synchronization in the nodes of a network with non-linear nodes with chaotic couplings where the synchronization manifold is a result of the collective behavior and is not known a priory. We consider a simple scale-free like network with
. The considered systems for the nodes are the following Lorenz systems given by
where
,
and
, and
i indicates the node.
The topology of the network illustrated in
Figure 1 is given by the coupling matrix
and is defined as
To model the variation of the couplings a chaotic unpredictable behavior is considered. The main idea is to induce couplings that are not defined by the connected nodes, this is, couplings that are dictated by an external system or imposed externally. To this end, the Rössler chaotic system is proposed as an external coupling provider since the Rössler system possesses a well-known chaotic behavior. The diagonal elements of the coupling matrix (
8) are defined as in (
6), functions that comprise the rest of the matrix are given by external functions of the generalized Rössler system given by:
where
with
, one for each coupling, and the parameters
,
,
, and
. In order to have positive chaotic functions, each solution of the Rössler system is modified by adding a positive offset, this is, it is known that if the couplings alternates its sign along time, there is a disconnection for a short period of time and more over the behavior of the network could be unstable, therefore a modification of these coupling functions is needed, therefore, the couplings are proposed as
, with
,
and
,
is chosen in such a way that the couplings are positive for all time and as a control parameter in the stability of the synchronization error (which will discuss later). Thus in this way there are five different signals as shown in
Figure 2. The projections of the chaotic attractor which define the couplings are illustrated in
Figure 3. In
Figure 4 the error of the state
are illustrated, where the synchronization error is calculated using the arithmetic average as follows
where
is the mean value of the states.
The behavior of the systems is presented in
Figure 5, where the synchronization is achieved even if the couplings are chaotic. In this example, the couplings are non autonomous and unpredictable, thus for some periods the connection between nodes is strong whereas the same coupling debilitates a time after. It is easy to propose any non autonomous couplings, perhaps the minimum value is positive, in some cases could be zero but the nodes lose the connection. The offset is given to satisfy that the eigenvalues of the connection matrix are one zero and the rest of the same sign. In this case, the error between states in the nodes does not define the couplings. This is, the couplings are independent of the behavior in the nodes. For these external non autonomous couplings the behavior of the
is depicted in
Figure 6.
3.2. Network Synchronization with Evolutionary Coupling
As previously stated, one of the biggest problems when trying to model real-world networks is that most of them have inherent structural dynamics. These dynamics are usually modeled as an external system or algorithm that sets coupling values for the network. That is—in most cases—an oversimplification of the problem. In this case, we propose a model with evolving couplings: couplings that vary accordingly to the states of the nodes they connect (which may be seen as their environment).
Note that this definition of the coupling is now a function of the connected systems, whereas, in the case of the non autonomous couplings, they were functions that do not depend on the network’s systems behavior. In this sense, the network equation is given by
where
are the state variables for the node
i.
is the evolving coupling between node
i and
j. The evolving couplings are defined as follows: if there is connection between nodes
i and
j at time
t, then
, otherwise
for
. The diagonal elements in the matrix
are defined as follows
As in the case of the varying couplings, for this case the following assumptions are considered the functions that comprise must be continuous, bounded and positive definite. Consider that the coupling function between node i and j is given by .
We present the complete synchronization in the nodes of a network of Lorenz systems with evolving couplings where the synchronization manifold is a result of the collective behavior and is not known a priory.
The topology of the network is the same as in
Figure 1 and the coupling matrix is given by
defined as follows
where the diagonal elements are defined by (
12), and the functions that comprise the rest of the matrix are given by a function of the states on the nodes.
Figure 7 illustrates the couplings time evolution, note that this couplings are close to zero when the behavior between nodes are closed each other, but in such a case it seems that the nodes are weakly connected and synchronization could be lost, when synchronization is weak the coupling increases and the synchronization is recovered. Also note that in
Figure 7c), the magnitude of the coupling is greater than those of
Figure 7a,b,d) this is since
is the coupling that connects nodes 1, 2 and 3 to 4, 5 and 6 and requires a little more coupling force.
Now, with this configuration of the network the systems were initiated without connecting them, this is there is a free behavior until
s.
Figure 8 illustrates the synchronization of the network of Lorenz systems, as it can be seen the synchronization is sustained despite the value of the evolutionary couplings.
To verify the stability of the synchronization state, from the Proposition 1, the parameter
is obtained in order to reach exponential stability of the synchronization regime. First the Jacobian matrix of system (
7) at the virtual solution
is given by
Matrix
with
. Now from Equation (
3) the inequality is given by
Therefore, to verify the exponential stability of the synchronization state a value for parameter
satisfying
is required.
To determine the stability of the synchronization state, let us depart from the Sylvester’s criterion over the matrix (
15), from where if the leading principal minors satisfy
for
then matrix (
15) is negative definite. For this case
are given as follows
for these inequalities to hold, note that
,
and considering
and
satisfies the criterion and the matrix is negative definite, since
, and for this example
is illustrated in
Figure 9; note that the eigenvalue does not satisfies the criterion for exponential stability for all time, however the synchronization is still stable. This behavior in the eigenvalue is due to the coupling function, since the coupling evolves accordingly to the behavior of the connected nodes, note that when the coupling in
Figure 7 are close to zero, the value for
is also close to zero which means that the nodes are weakly connected and then the synchronization loses stability. This behavior on the stability of the synchronization is consequence of the evolutionary character of the couplings since these are not influenced or imposed to the network.
It is important to express that this result on the value of the parameter d also applies to the case of time variation couplings (chaotic couplings), since the analysis is over the Jacobian matrix and are independent to the connectivity matrix, therefore, for that case, the exponential stability is completely achieved, this is for all .
4. Discussion
The synchronization between nodes is achieved in both cases, one for imposed or external (exogenous), and the other defined by the connected nodes (endogenous). However, there is an important issue to be discussed, even if the synchronization is achieved, for the case of non autonomous exogenous couplings, the synchronization state is exponentially stable, in
Figure 6 the behavior of the second larges eigenvalue is illustrated and the condition
is satisfied for all time. On the other hand, the synchronization with evolutionary couplings is achieved, however, the synchronization is not exponentially stable for all time, this is the condition
is not satisfied for all
as it is shown in
Figure 9, where the red line dives the exponential stability from the stability of the synchronization state, this is due to the parameters which define the value
d. The synchronization is not lost since
and there are not isolated nodes.
Note that for non autonomous exogenous couplings, the synchronization is controlled by the selection of the couplings, in other words it is always possible to fix the value of the couplings in order to satisfy the condition for all whereas in case of evolutionary couplings, the synchronization is reached by the collectivity of the nodes and as the couplings are not controlled or imposed the synchronization can or cannot be exponentially. These scenarios are directly consequence of the couplings properties.
It is important to mention that the synchronization is achieved in presence of non autonomous couplings or couplings depending on the behavior of the network, the stability of the synchronization manifold is achieved depending on the magnitude of the coupling strength between nodes, however, the meaning is different, on the one hand, non autonomous couplings describes an imposition of the behavior of the coupling strength, this is, that the nodes in the network obey the strength and the nodes do not have any possibility of changing their coupling, on the other hand, when the couplings depend on their own dynamics they define the couplings strength along time, this is called the evolution of the coupling and it is an intrinsic property without any external influence. An example for the non autonomous coupling case which also can be seen as an imposition of the strength of the couplings in a network is found in treatments where deep brain stimulation is performed to treat Alzheimer’s, Parkinson’s disease, tremor, and so forth, this stimulation is imposed or dictated to the network [
18,
19,
20]. On the other hand, evolving couplings are considered when the nodes define their own coupling strength according to their requirements, this is, the relationships between agents are determined by them, for instance in social networks [
21], in problems of formation control of UAV [
26] and so forth.