Chaos Synchronization in Time-Dependent Duplex Networks

Liu, Qian; Han, Wenchen; Lei, Lixing; Dai, Qionglin; Yang, Junzhong

doi:https://doi.org/10.1155/2019/6159365

Complexity

On this page

Abstract Introduction Model Conclusion Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2019 | Article ID 6159365 | https://doi.org/10.1155/2019/6159365

Chaos Synchronization in Time-Dependent Duplex Networks

Qian Liu,¹Wenchen Han,²Lixing Lei,¹Qionglin Dai,¹and Junzhong Yang¹

Academic Editor: Cornelio Posadas-Castillo

Received05 Jan 2019

Accepted20 Mar 2019

Published04 Apr 2019

Abstract

Though the complete chaos synchronization on single-layer networks has been well understood, it is still a challenge on multiplex networks. In this work, we study the complete chaos synchronization on time-dependent duplex networks in which interaction pattern among oscillators alternates periodically between two single-layer networks. The alternations between two layers are characterized by the offset strength and the switching frequency . We find that there are two dynamical regimes depending on . For high , the critical for the stable complete synchronization is independent of and the fast-switching approximation suggests that the time-dependent duplex networks can be approximated by the time-independent duplex networks with effective coupling strength. For low , the critical depends on nonmonotonically. At extremely low , the estimation of the critical can be obtained by a single-mode approximation taking one dominant transversal network mode to complete synchronization into considerations.

1. Introduction

Collective behaviors are ubiquitous in natural systems and human society. Among these, synchronization, one of the most interesting collective behaviors, has been investigated since the early days in natural science and has drawn great attention in the last decades [1–4]. There are various forms of synchronization, such as complete synchronization (CS) [5], cluster synchronization [6], remote synchronization [7], partial synchronization [8], phase synchronization [9], generalized synchronization [10], and lag synchronization [11]. Among these various types of synchronization, CS is the simplest one and its stability has been well analyzed by using the master stability function method [12, 13]. When complex networks are involved, it has been found that both the coupling function and the network structure can affect the stability of CS.

Previous studies on network science focused on single-layer networks. However, there are many complex systems ranging from social systems to technological systems which have to be modeled as multilayer networks other than single-layer ones [14–17]. For example, in social networks people may interact with each other through different channels such as WeChat, Twitter, and Facebook [18]. Each channel consists of an independent layer on which people interact with the coupling function of their own layer and all layers supported by different channels are integrated into a multilayer network. Among multilayer networks, multiplex networks refer to the ones in which different layers share the same group of nodes and each layer of the networks owns its unique interaction pattern [19–21]. Recently, CS has been investigated in multilayer networks [14, 22–24]. Different from single-layer networks on which the onset of the instability of CS is induced by the most unstable transversal network modes, it is still a challenge to figure out the mechanism for CS on multiplex networks to lose its stability. Except for some peculiar multiplex networks, there are no rigorous results for the stability of CS though some approximation methods have been proposed [25].

Interaction patterns may be time-dependent in many situations such as power transmission system [26], consensus problem [27], and person to person communications [28]. In modeling, the time-varying interaction patterns may be realized in several ways, for example, by rewiring connections in networks in the course of time [29], by assuming oscillators to move in space [30], and by allowing the switch of network topologies among several given interaction patterns [31]. Synchronization on time-varying networks has been investigated when oscillators may move in space [32, 33]. It has been shown that synchronization time may depend nonmonotonically on the mobility of oscillators and the mechanism driving synchronization is different for different dynamical regimes [34–37]. However, in these literatures concerning with synchronization on time-varying networks, the coupling function is always same. In this paper, we construct a time-varying duplex networks by switching the interactions pattern and coupling function among oscillators between two single-layer networks with different coupling functions and study CS among chaotic Lorenz oscillators on them.

The paper is organized as follows. In Section 2, we present the model. The time-varying duplex networks are characterized by the parameters, the switching frequency and the offset strength . In Section 3, we first consider two single-layer networks and study CS on each of them. Then we consider the time-varying duplex networks constructed from these two single-layer ones and investigate the dependence of CS on the switch frequency and the offset strength. The numerical results and some discussions are presented. In Section 4, a summary is made.

2. Model

We consider a network composed of chaotic oscillators whose motion equations are governed by where is the 3-dimensional state variable of oscillator . We consider chaotic Lorenz dynamics which is described as . and with the step function if and otherwise. and are the adjacency matrices in the two layers, respectively. () when there exists a connection between oscillators and in the layer () and, otherwise, (). The parameter is the offset strength. The time-dependent duplex network reduces to the single layer when and to the single layer for . At any given time, oscillators interact with each other either through layer or through layer for . Increasing the offset strength prolongs the layer to be the acting network. The parameter is the switching frequency accounting for the rate of the network topology alternating between the layers and . In the layers and , oscillators interact with each other through the coupling function and . Unless specified, we set and . and account for the coupling strengths in the layers and , respectively.

3. Numerical Simulations

Before investigating CS on time-dependent duplex networks, we take a look at CS on single-layer networks. For simplicity, we let . We generate two random networks and which are illustrated in the insets in Figures 1(a) and 1(b). Using the master stability function method, the stability of CS on or can be analyzed by decoupled linear differential equations with . are the eigenvalues of the Laplacian matrix constructed from and we suppose . Except for the eigenvalue accounting for the synchronous network mode, other nonzero represent the transversal network modes to CS. and are the Jacobian matrices of the corresponding functions evaluated at CS. The stability of CS requires that the largest transversal Lyapunov exponent, , for each nonzero to be negative. Figures 1(a) and 1(b) show for transversal network modes for the single-layer networks and , respectively. Clearly, the network supports stable CS for while CS is always unstable on . For an arbitrary perturbation to CS, it can always be expanded as a linear combination of all transversal network modes. With time elapsing, it is the network mode with the largest transversal Lyapunov exponent that takes responsibility for determining the stability of CS. According to Figures 1(a) and 1(b), we may find that, for the layer , the network mode is the dominant one which always has the largest . On the other hand, the dominant network mode depends on the coupling strength for layer . For example, the dominant mode is for while it changes to for . Furthermore, we calculate the Lyapunov exponent spectrum for the model dynamics equation (1) on the networks and . The first four largest Lyapunov exponents against the coupling strength and are presented in Figures 1(c) and 1(d), respectively. Beyond the critical coupling strength , CS is established for the network since all Lyapunov exponents are negative except for the largest Lyapunov exponent taking the value close to that at and the second largest Lyapunov exponent to be zero. On the other hand, CS is impossible for the network and the model dynamics may be either chaotic or time-independent for sufficiently strong .

(a)

(b)

(c)

(d)

Then, we consider the time-dependent duplex networks constructed from the networks and in Figure 1 and focus on the effects of the parameters and on the model dynamics in (1). To measure the synchronization in the model, we consider the synchronization error , which is given by refers to the Euclidean norm between oscillators and and stands for the time average. For each parameter combination, the synchronization error is averaged over time units after a transient time around time units.

We first set where CS is an attractor on the layer and where chaotic desynchronous state is realized on the layer . Increasing means that oscillators spend more time on the layer , which tends to disfavor CS on the time-dependent duplex network. Compatible with the intuition, Figure 2(a) shows that always increases with . There seems to exist a critical for each above which becomes nonzero and, to be stressed, increases rapidly at low and seems to saturate to a value at high . Figure 2(b) presents against for different . For close to 1, keeps nonzero and CS is impossible. On the other hand, is always possible for high when is not close to 1. Next, we set where equilibrium is the attractor on the layer . The dependence of against and are presented in Figures 2(c) and 2(d) and we observe the similar behaviors to those at . Thus, CS on the time-dependent duplex network is insensitive to chaotic or regular dynamics on the isolated single-layer networks. To be mentioned, Figure 1(d) suggests a quenched state at on the network where the maximum Lyapunov exponent is negative. Figure 2(b) further suggests that the quenched state is actually an inhomogeneous one according to the nonzero at .

(a)

(b)

(c)

(d)

To get an overview on CS in the model (1), we present in the plane of and . Figures 3(a) and 3(b) display the results for the parameter combinations (, ) and (, ), respectively. It is clear that CS may be realized in a large range of domain in the parameter plane, for example, high and not close to 1. Actually, there exist two dynamical regimes depending on . For high , the critical separating CS from desynchronous states is roughly independent of , which is around 0.8 for (, ) and 0.6 for (, ). In contrast, for low , the parameter regimes supporting CS are narrow. The enlargement of the parameter plane in Figures 3(c) and 3(d) shows irregular boundary between CS and desynchronous states at low .

Figure 3

The contour plot of the synchronization error in the plane of and in (a) for and and in (b) for and . (c) and (d) The enlargement view of (a) and (b) in the range of , respectively. (e) The phase diagram in the plane of and for the time-independent duplex network with two layers represented by and . The shaded area corresponds to the CS region. The dot lines denote the effective and acquired by the fast-switching approximation. (f) The critical are plotted against with in black, in blue, and acquired from the synchronization error in red. In addition, against in the time-dependent duplex network with replaced by is plotted in green. for the solid symbols and for the open symbols. The inset in (f) shows the single-layer network .

The results at high can be understood as follows. For high , the interaction pattern among oscillators switches between the layers and rapidly. In this case, the fast-switching approximation [34, 38], which averages out the effect of alternation between and , may be applied. Consequently, the interaction pattern among oscillators may be approximated by a time-independent duplex network with and being its two layers and with the effective coupling strengths and where and is the time duration that oscillators interact with each other in the time-dependent network through . Just like social networks where each person may simultaneously communicate with others via Facebook network and Twitter network, in this static duplex network as an approximation, each oscillator may simultaneously interact with others through two different networks assisted by the coupling functions and , respectively. Since is determined by , it is unchanged with for a given , which means that and are determined by uniquely. We present the phase diagram in the plane of and for the time-independent duplex network with and in Figure 3(e) to show the CS region. Then we present the effective coupling strengths and acquired from the coupling strengths and on the time-dependent duplex network by changing from to with the increment . The green dots in Figure 3(e) are for the set (, ) while the red dots for the set (, ). According to the plot, we find the critical is around for and for , which are in agreement with Figures 3(a) and 3(b).

On the other hand, the dynamics is quite complicated at low . However, for the extremely low , we could present some plausible explanations. For the extremely low , the system has to spend sufficiently long time on one layer before switching to the other one. Consequently, the transversal network modes on the layer and on the layer will have different effects on model dynamics. To sort out these effects from different network modes, we firstly take the transversal network modes with the largest transversal Lyapunov exponent on each layer into considerations, for example, on the layer and, on the layer , for and for . Supposing the transversal Lyapunov exponent of the network mode in the layer to be and in the layer to be , any perturbation away from CS would evolve in the form . For the perturbation to decay in one period , should be larger than . That is, it requires with obeying against is presented in Figure 3(f) for and . Clearly, at is higher than that at since at is closer to 0. Then, we numerically calculate against for different and the critical for the onset of desynchronization acquired in this way is also presented in Figure 3(f). There is strong discrepancy between and in the range of where may reach a much higher value than . The discrepancy lies at the inappropriate approximation that only the network mode with the largest in each layer is considered. Actually, all transversal network modes collectively determine how the perturbation away from synchronous chaos evolves. Considering the fact that the inner products, , of network modes on different layers determine the redistribution of the perturbation on over at the time when the two layers are switching, we think that, in the approximation of one network mode, which network mode in one layer should be taken into considerations depends on its inner products with the network modes on the other layer and its transversal Lyapunov exponents. Along this line, the network mode is the one in the layer crucial to . On the layer , though has the largest Lyapunov exponent for , the product is much smaller than and . Therefore, the dominant network mode on the layer is the mode for and otherwise according to Figure 1(d). Substituting with the transversal Lyapunov exponent with the corresponding dominant mode in (4), we obtain a new critical (denoted as ). As Figure 3(f) shows, may reach a much higher value in the range of .

Actually, the above analysis can be further confirmed by replacing with another network, the layer [shown in the inset of Figure 3(f)], which keeps the overall interaction pattern unchanged but with rearranged oscillators on the layer. The Laplacian of has the same eigenvalues with and, therefore, the stability of CS on is exactly the same as that on . Different from , the product between and is around 0.66, which assures that the network mode is the dominant one in determining . Consequently, , acquired according to , in the range of is greatly lowered and the discrepancy between and is greatly reduced, as shown in Figure 3(f).

We have found that CS is possibly established on a time-dependent duplex network when one of the layers does not support CS on its own. Actually, CS can be realized on time-dependent duplex networks even when two layers do not allow for the stable CS on their own. To illustrate it, we consider the coupling functions and where CS is unstable on either or . We present the synchronization error in the plane of and in Figure 4(a) with and Figure 4(b) with and . Except for the parameter regimes close to and , CS can always be produced at high where the critical for CS is insensitive to . Figure 4(c) shows the stability diagram in the plane of and for the time-independent duplex network. The critical in Figures 4(a) and 4(b) can be acquired by plotting the effective coupling strengths and at different in Figure 4(c).

(a)

(b)

(c)

4. Conclusion

In this work, we constructed time-dependent duplex networks in which the network topology alternates periodically between two single-layer networks. The temporal properties of the networks are characterized by the offset strength and the switching frequency . We studied the complete chaos synchronization on this type of time-dependent duplex networks. By monitoring the synchronization error, we find that there are two dynamical regimes depending on the switching frequency . For high where the fast-switching approximation is valid, the time-dependent duplex networks can be approximated by time-independent duplex network with effective coupling strength. The critical for the onset of stable CS is independent of . For low , the critical strongly depends on and the coupling strengths on the two layers. At extremely low , we find that the dependence of the critical on the coupling strengths can be roughly explained based on the picture in which only one dominant network mode in each layer is taken into considerations. We found that which network mode is the dominant one is determined by its transversal Lyapunov exponent and its inner product with the dominant one on the other layer.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. 11575036 and Grant no. 11805021).

References

A. Arenas, A. D. Guilera, J. Kurths, Y. Moreno, and C. Zhou, “Synchronization in complex networks,” Physics Reports, vol. 469, no. 3, pp. 93–153, 2008.
View at: Publisher Site | Google Scholar | MathSciNet
J. Zhou and T. Chen, “Synchronization in general complex delayed dynamical networks,” IEEE Transactions on Circuits and Systems I: Regular Papers, vol. 53, no. 3, pp. 733–744, 2006.
View at: Publisher Site | Google Scholar | MathSciNet
X. F. Wang and G. Chen, “Synchronization in scale-free dynamical networks: robustness and fragility,” IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications, vol. 49, no. 1, pp. 54–62, 2002.
View at: Publisher Site | Google Scholar | MathSciNet
J. Lu, X. Yu, and G. Chen, “Chaos synchronization of general complex dynamical networks,” Physica A: Statistical Mechanics and its Applications, vol. 334, no. 1-2, pp. 281–302, 2004.
View at: Publisher Site | Google Scholar | MathSciNet
W. Zou and M. Zhan, “Complete synchronization in coupled limit-cycle systems,” Europhysics Letters. EPL, vol. 81, no. 1, Art. 10006, 6 pages, 2008.
View at: Publisher Site | Google Scholar | MathSciNet
V. N. Belykh, I. Belykh, and E. Mosekilde, “Cluster synchronization modes in an ensemble of coupled chaotic oscillators,” Physical Review E, vol. 63, Article ID 036216, 2001.
View at: Google Scholar
L. V. Gambuzza, A. Cardillo, A. Fiasconaro, L. Fortuna, J. Gómez-Gardeñes, and M. Frasca, “Analysis of remote synchronization in complex networks,” Chaos, vol. 23, no. 4, Article ID 043103, 2014.
View at: Google Scholar
B. Ao and Z. Zheng, “Partial synchronization on complex networks,” EPL (Europhysics Letters), vol. 74, no. 2, pp. 229–235, 2006.
View at: Publisher Site | Google Scholar | MathSciNet
M. G. Rosenblum, A. S. Pikovsky, and J. Kurths, “Phase synchronization of chaotic oscillators,” Physical Review Letters, vol. 76, no. 11, pp. 1804–1807, 1996.
View at: Publisher Site | Google Scholar
N. F. Rulkov, M. M. Sushchik, L. S. Tsimring, and H. D. I. Abarbanel, “Generalized synchronization of chaos in directionally coupled chaotic systems,” Physical Review E, vol. 51, no. 2, pp. 980–994, 1995.
View at: Publisher Site | Google Scholar
M. G. Rosenblum, A. S. Pikovsky, and J. Kurths, “From phase to lag synchronization in coupled chaotic oscillators,” Physical Review Letters, vol. 78, no. 22, pp. 4193–4196, 1997.
View at: Publisher Site | Google Scholar
L. M. Pecora and T. L. Carroll, “Master stability functions for synchronized coupled systems,” Physical Review Letters, vol. 80, no. 10, pp. 2109–2112, 1998.
View at: Publisher Site | Google Scholar
J. Yang, G. Hu, and J. Xiao, “Chaos synchronization in coupled chaotic oscillators with multiple positive lyapunov exponents,” Physical Review Letters, vol. 80, no. 3, pp. 496–499, 1998.
View at: Publisher Site | Google Scholar
C. I. Del Genio, J. Gómez-Gardeñes, I. Bonamassa, and S. Boccaletti, “Synchronization in networks with multiple interaction layers,” Science Advances, vol. 2, no. 11, Article ID e1601679, 2016.
View at: Google Scholar
J. Gómez-Gardeñes, I. Reinares, A. Arenas, and L. M. Floría, “Evolution of cooperation in multiplex networks,” Scientific Reports, vol. 2, no. 1, Article ID 620, 2012.
View at: Publisher Site | Google Scholar
A. Cardillo, J. Gomez-Gardeñes, M. Zanin et al., “Emergence of network features from multiplexity,” Scientific Reports, vol. 3, Article ID 1344, 2013.
View at: Google Scholar
C. Hong, J. Zhang, X. Cao, and W. Du, “Structural properties of the chinese air transporation multilayer network,” Chaos Solitons and Fractals, vol. 86, pp. 28–34, 2016.
View at: Google Scholar
M. Salehi, R. Sharma, M. Marzolla, M. Magnani, P. Siyari, and D. Montesi, “Spreading processes in multilayer networks,” IEEE Transactions on Network Science and Engineering, vol. 2, no. 2, pp. 65–83, 2015.
View at: Publisher Site | Google Scholar
M. De Domenico, C. Granell, M. A. Porter, and A. Arenas, “The physics of spreading processes in multilayer networks,” Nature Physics, vol. 12, no. 10, pp. 901–906, 2016.
View at: Publisher Site | Google Scholar
Z. Wang, L. Wang, A. Szolnoki, and M. Perc, “Evoluationary games on multilayer networks: a colloquium,” The European Physical Journal B, vol. 88, Article ID 124, 2015.
View at: Google Scholar
S. Boccaletti, G. Bianconi, R. Criado et al., “The structure and dynamics of multilayer networks,” Physics Reports, vol. 544, no. 1, pp. 1–122, 2014.
View at: Publisher Site | Google Scholar
M. Liu, Y. Shao, and X. Fu, “Complete synchronization on multi-layer center dynamical networks,” Chaos, Solitons & Fractals, vol. 41, no. 5, pp. 2584–2591, 2009.
View at: Publisher Site | Google Scholar | MathSciNet
F. Sorrentino, “Synchronization of hypernetworks of coupled dynamical systems,” New Journal of Physics, vol. 14, no. 3, Article ID 033035, 2012.
View at: Publisher Site | Google Scholar
Y. Li, X. Wu, J.-A. Lu, and J. Lu, “Synchronizability of duplex networks,” IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 63, no. 2, pp. 206–210, 2016.
View at: Publisher Site | Google Scholar
W. Han and J. Yang, “Approximation methods for the stability analysis of complete synchronization on duplex networks,” Chaos, Solitons & Fractals, vol. 106, pp. 23–27, 2018.
View at: Publisher Site | Google Scholar | MathSciNet
M. L. Sachtjen, B. A. Carreras, and V. E. Lynch, “Disturbances in a power transmission system,” Physical Review E, vol. 61, no. 5, pp. 4877–4882, 2000.
View at: Google Scholar
R. Olfati-Saber, J. A. Fax, and R. M. Murray, “Consensus and cooperation in networked multi-agent systems,” Proceedings of the IEEE, vol. 95, no. 1, pp. 215–233, 2007.
View at: Publisher Site | Google Scholar
J. Onnela, J. Saramaki, J. Hyvonen et al., “Structure and tie strengths in mobile communication networks,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 104, no. 18, pp. 7332–7336, 2007.
View at: Publisher Site | Google Scholar
S. Boccaletti, D.-U. Hwang, M. Chavez, A. Amann, J. Kurths, and L. M. Pecora, “Synchronization in dynamical networks: evolution along commutative graphs,” Physical Review E, vol. 74, no. 1, Article ID 016102, 5 pages, 2006.
View at: Publisher Site | Google Scholar | MathSciNet
M. Frasca, A. Buscarino, A. Rizzo, L. Fortuna, and S. Boccaletti, “Synchronization of moving chaotic agents,” Physical Review Letters, vol. 100, no. 4, Article ID 044102, 2008.
View at: Google Scholar
D. B. Kingston and R. W. Beard, “Discrete-time average-consensus under switching network topologies,” in Proceedings of the 2006 American Control Conference, pp. 3551–3556, USA, June 2006.
View at: Google Scholar
J. Gómez-Gardeñes, V. Nicosia, R. Sinatra, and V. Latora, “Motion-induced synchronization in metapopulations of mobile agents,” Physical Review E, vol. 87, no. 3, Article ID 032814, 2013.
View at: Google Scholar
Y.-H. Eom, S. Boccaletti, and G. Caldarelli, “Concurrent enhancement of percolation and synchronization in adaptive networks,” Scientific Reports, vol. 6, Article ID 27111, 2016.
View at: Google Scholar
N. Fujiwara, J. Kurths, and A. Díaz-Guilera, “Synchronization in networks of mobile oscillators,” Physical Review E, vol. 83, no. 2, Article ID 025101, 2011.
View at: Publisher Site | Google Scholar
N. Fujiwara, J. Kurths, and A. Díaz-Guilera, “Synchronization of mobile chaotic oscillator networks,” Chaos, vol. 26, no. 9, Article ID 094824, 8 pages, 2016.
View at: Publisher Site | Google Scholar | MathSciNet
L. Prignano, O. Sagarra, and A. Díaz-Guilera, “Tuning synchronization of integrate-and-fire oscillators through mobility,” Physical Review Letters, vol. 110, no. 11, Article ID 114101, 2013.
View at: Google Scholar
A. Beardo, L. Prignano, O. Sagarra, and A. Díaz-Guilera, “Influence of topology in the mobility enhancement of pulse-coupled oscillator synchronization,” Physical Review E, vol. 96, no. 6, Article ID 062306, 2017.
View at: Publisher Site | Google Scholar
I. V. Belykh, V. N. Belykh, and M. Hasler, “Blinking model and synchronization in small-world networks with a time-varying coupling,” Physica D, vol. 195, no. 1-2, pp. 188–206, 2004.
View at: Publisher Site | Google Scholar | MathSciNet

Copyright

Copyright © 2019 Qian Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies