#### Synchronizability of Small-World Networks Generated from a Two-Dimensional Kleinberg Model

Synchronizability of Small-World Networks Generated from a Two-Dimensional Kleinberg Model
Yi Zhao,1 Jianwen Feng,1 and Jingyi Wang2
1College of Mathematics and Computational Science, Shenzhen University, Shenzhen 518060, China
2College of Information and Engineering, Shenzhen University, Shenzhen 518060, China
Received 5 June 2013; Accepted 26 July 2013
Academic Editor: Wenwu Yu
Copyright © 2013 Yi Zhao 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.
Abstract
This paper investigates the synchronizability of small-world networks generated from a two-dimensional Kleinberg model, which is more general than NW small-world network. The three parameters of the Kleinberg model, namely, the distance of neighbors, the number of edge-adding, and the edge-adding probability, are analyzed for their impacts on its synchronizability and average path length. It can be deduced that the synchronizability becomes stronger as the edge-adding probability increases, and the increasing edge-adding probability could make the average path length of the Kleinberg small-world network go smaller. Moreover, larger distance among neighbors and more edges to be added could play positive roles in enhancing the synchronizability of the Kleinberg model. The lorentz oscillators are employed to verify the conclusions numerically.
1. Introduction
A complex network is a large set of nodes (or vertices) connected by a set of links (or edges) such as coupled biological and chemical system, neural networks, social interacting species, the Internet, and the World Wide Web. Recently, increasing interest has been devoted to the study of collective behaviors in complex networks for its widely applications in real world. Among the studies on the complex network, synchronization phenomena attract the interests of most scientists and engineers. Loosely speaking, synchronization is the process in which two (or more) dynamical systems seek to adjust a certain prescribed property of their motion to a common behavior in the limit as time tends to infinity either by virtue of coupling or by forcing. Synchronization of complex networks is an important mathematical problem in both the physical and biological sciences since it has potential applications to diverse fields such as communications security, seismology, and parallel image processing [1–6].
Complex networks could be classified as many types. Among them, synchronization on small-world networks has attracted considerable attention since the pioneering work of Stanley Milgram in the 1960s [7]. A small-world network can be generated by either random edge-rewiring, which gives WS small-world network [8], or random edge-adding, which yields the NW small-world network [9]. The Kleinberg small-world network, in which the edge-adding probabilities are proportional to the length of the edge to be added, could be seen as a more general NW small-world network.
Various literatures have already been devoted to the studies on synchronizability of small-world networks. In the research articles [10–13], the synchronizability of a small-world network generated by randomly adding a fraction of long-range shortcuts to a ring network is investigated. It can be deduced from the theoretical analysis and numerical simulation that the synchronizability of the small-world network becomes stronger as the edge-adding probability grows larger. In [14], Tang et al. found that the synchronizability of the network as a function of the distance is fluctuant and there exist some distances that have almost no impact on the synchronizability when they investigated the impact of edge-adding number and edge-adding distance on both synchronizability and average path length of NW small-world networks generated from ring networks via random edge-adding. Moreover, the relationship between the synchronizability and the average path length of a small-world network is studied in [15–18]. The analysis and numerical simulations show that the synchronizability of the small-world network grows as increases and the average path length becomes smaller as goes larger. Therefore, it can be deduced that the decreasing in the average path length may result in the increasing synchronizability. These phenomena are interesting, and a natural question is that whether other small-world networks have similar properties, which motivates us to take a two-dimensional Kleinberg small-world network [19] as an example and investigate the impact factors of such network. It should be mentioned that the synchronizability of an undirected Kleinberg small-world network was investigated in [20]. However, the Kleinberg model is built as a directed network in [19]. Thus, the directions of the edge-adding in building the Kleinberg model are considered in this paper. Moreover, [20] only discussed the relationship between the edge-adding probability and the synchronizability of the small-world network, while in this paper the three parameters of the Kleinberg model, namely, the distance of neighbours, the number of edge-adding, and the edge-adding probability, would be analyzed for their impacts on its synchronizability and average path length. Actually, this paper improves the results in [20].
In this paper, we investigate the impacts of the distance of neighbors, the number of edge-adding and the edge-adding probability on the synchronizability of the Kleinberg small-world network. The Kleinberg small-world network is an two-dimensional one. We add edges on the nodes with certain probability . Then, we could get some conclusions about impact factor on the synchronizability and the average path length of the Kleinberg small-world network, which are complementary to the studies on the synchronizability of the small-world networks.
2. Preliminaries
First of all, we build a Kleinberg small-world network in the way introduced in [19]. Figure 1 also comes from [19], and we redraw it in our case to be studied. A Kleinberg small-world network is composed of the set of lattice points in an square, which are denoted as . The lattice distance between two nodes and is defined to be the number of “lattice steps,” which could be written as . Let and be positive integers. The node is connected with every other node within lattice distance , and we name it local contact. We also construct edges from to other nodes using independent random trials, which are called the long-range contacts. The probability of edge connected between and is proportional to , where is a given constant. Precisely speaking, this probability of the connections between and is denoted as , and Figure 1 shows basic structures of a Kleinberg small-world network. In Figure 1(b), there are two long-rang contacts from a node “” to a node “” and a node “”. Actually, two long-range contacts are added to every node in this network if .
Figure 1: (a) A two-dimensional Kleinberg network with nodes, the distance of local contacts , and it has no long-range contacts thus . (b) The contacts of a node “” with the distance of local contacts and the number of long-range contacts . “” and “” are the two long-range contacts.
Actually, this model could be interpreted in the point of “geography” in [19]. Individuals live on a grid and know their neighbors for some number of steps in all directions; they also have some number of acquaintances distributed more broadly across the grid. If we fixed and and let the value of the exponent vary, we would have a one-parameter family of network models. When , the uniform distribution over long-range contacts could be obtained, which means long-range contacts are chosen independently of their position on the grid. In this sense, the Kleinberg small-world network could be seen as a kind of NW small-world network. As increases, the long-range contacts of a node become more and more clustered in their vicinity on the grid. Thus, could be seen as a basic structural parameter measuring how widely “networked” the underlying society of nodes is. Considering that could reveal some basic properties of the Kleinberg network, in this paper, we investigate the effect of adding-edges probabilities, namely, the effect of the parameter on the synchronizability of the Kleinberg small-world network. Moreover, two other parameters of the Kleinberg small-world network, namely, the distance of neighbors and the number of edge-adding, are also considered for their influences of synchronizability.
Next, the synchronizability analysis of the complex dynamical system follows [21]. The complex dynamical system considered in this paper consists of coupled continuous-time nonlinear oscillators. Since chaotic behaviors are common since the intrinsic nonlinearity exists in each individual oscillator, chaotic synchronization is addressed by choosing the parameters of each oscillator such that it exhibits a chaotic attractor in order to be general. Then, the network of identical dynamics nodes considered in this paper can be written as Here, if the node is located at in the network, we denote the index , and thus, . is the state vector of the th node in all nodes. is a positive constant coupling strength. is a well-defined nonlinear function and is a coupling function. is a coupling matrix determined by the connection topology. That is, if the node and have connections, and otherwise. Moreover, it is required that for any . Note that is not necessarily symmetric since our considered network is directed. Moreover, there is only one zero eigenvalue of the matrix such that the eigenvalues can be sorted as Let , and is automatically a solution of (2). Then, the synchronous state is defined as . Let . For the system described by (2), the variational equations governing the time evolution of the set of infinitesimal vectors about the synchronous solution are where and are the Jacobian matrices of the corresponding vector functions evaluated at , respectively. Let . Then, (4) can be rewritten as By using Jordan transformation with respect to the coupling matrix , we have where is composed of the eigenvectors of . And in (6) are blocks of the form where is one of the eigenvalues of .
Letting and employing (6), (5) could also be written as Then, each block of the Jordan canonical form corresponds to a subject of these columns in , which obeys a subset of equations in (8). For instance, if block is , and suppose the corresponding columns of are denoted by , which could be seen as the modes of perturbations in the generalized eigenspace associated with eigenvalue , then the equations have the following form where . Each block of the previous decoupled equation (9) is structurally the same with only the factor of being different [21].
Define as a normalized coupling parameter that takes values in the set . Denote as the largest Lyapunov exponent, which is determined from (9). In order to achieve the synchronization of the network, is required to be negative. It can be explained that a small disturbance from the synchronization state will diminish exponentially so that the synchronous solution is stable. On the contrary, the synchronous solution is unstable and cannot be realized physically if is positive because small perturbations from the synchronous state will lead to trajectories that diverge from the state. For the reasons mentioned above, in order to achieve the synchronization of the coupled oscillator network (2), all normalized coupling parameters for should fall in a region in which should be negative. This region is called the synchronized region. In case that the synchronized region is bounded, namely, , then from (3), we have When the spread in the eigenvalue goes smaller, it becomes easier that all the numbers fall into the bounded region , which means that the synchronizability of the network is better. Thus, we need the ratio of eigenvalue satisfying to be smaller. In case that the synchronized region is unbounded, then for . Thus, the synchronizability of the network is better if the eigenvalue is smaller. Thus, and are used as the measure to evaluate the synchronizability of the network.
3. Influencing Factors the Synchronizability of a Network
An extensive numerical analysis is employed to investigate the influences of the distance of neighbors, the number of edge-adding, and the edge-adding probability on the network synchronizability. The Kleinberg small-world networks we considered have nodes and nodes, respectively. Identical dynamics are assumed for all the nodes in the Kleinberg small-world network.
Let , which means that each node in the Kleinberg small-world network is connected with its nearest neighbors in distance 1 or in distance 2. They form the local contacts. Then, long-range contacts are added with the probability defined in (1). Thus, different and would result in different corresponding Laplacian matrix. The eigenvalues of such Laplacian matrix could be calculated. Since the nodes with long-range contacts added are chosen randomly, the Laplacian matrix would be different on each trial. Thus, in the simulation, and 225 different realizations were performed and the results were averaged.
Let , respectively, for the network with nodes and nodes. The parameter in (1) on the probability of adding the long-range contacts is chosen from 0.1 to 10 with step size 0.1. Their corresponding Laplacian eigenvalues and as a function of are found, which was shown in Figures 2 and 3, respectively.
Figure 2: The relationship between and for the Kleinberg small-world network with (a) nodes and (b) nodes.
Figure 3: The relationship between and for the Kleinberg small-world network with (a) nodes and (b) nodes.
Figures 2 and 3 reveal that the values of and are continuously and monotonically increasing as increases. It means that the synchronizability of the Kleinberg small-world network becomes stronger as the edge-adding probability increases. Note that when , Kleinberg small-world network could be seen as the NW small-world network. This result corresponds with the observations in references [10–13]. Also, we can see from Figures 2 and 3 that the synchronizability is enhanced as the number of long-range contacts and the distance of neighbours increase. It could be explained that the average path length is reduced as more long-range contacts are built and more short-range contacts are constructed. We illustrate this relationship in Figure 4.
Figure 4: The average path length as a function of with different values of and for the Kleinberg small-world network with (a) nodes and (b) nodes.
It is well known that as the distance of the edges added increases, the synchronizability of the network becomes stronger since the average path length is shortened. The probability of edge-adding between and is proportional to , which is to say that the shorter distance of the edge-adding has larger probabilities than the longer distance of edge-adding for any given . Then, it can be concluded that if is fixed, the synchronizability of the Kleinberg small-world network would be better though the probability of edge adding will become smaller. Figures 2, 3, and 4 show that if is fixed, in other words, edges would be added into the network, the synchronizability becomes weaker as we take a larger value of . Then, it means that if distance of the edge-adding and the number of edges are fixed, the synchronizability is enhanced as the probability of the edge-adding increases. Meanwhile, in the Kleinberg small-world network, the probabilities of the long-range contacts decrease as the distance between two nodes increases, especially for large . As increases, the long-range contacts of a node become more and more clustered in their vicinity on the grid.
4. Numerical Simulation for Lorentz Oscillators
The linearly coupled Kleinberg small-world network containing identical Lorentz oscillators is used for numerical simulations. Such oscillators can be written as for . Here, , , and . Let In the simulation, we consider a network with ( nodes). Define the error term as If , the complex network achieves synchronization. We consider that the synchronizability would be better if goes to zero faster.
For the network (2), let the coupling strength , the distance of local contacts , the number of long-range contacts , and the parameter in the probability of edge-adding (1), respectively. Figure 5 shows how evolved in network with the same initial values chosen randomly in the interval . It can be seen from Figure 5 that the synchronizability of the network is enhanced as the number of long-range contacts increases, and the synchronizability of the network becomes better as goes smaller, which means that longer distance of the edge-adding could still enhance the synchronizability though the probabilities of longer distance of the edge-adding are small compared with the probabilities of shorter distance of edge-adding. Moreover, it can be seen that the synchronizability goes better as the number of edges added becomes larger.
Figure 5: The errors of synchronization among different Lorentz dynamics oscillators over time for , , and , respectively. The inset is an amplified figure on some time interval.
5. Conclusion
The impact factors of synchronizability of two-dimensional Kleinberg small-world network are investigated in this paper. Through mathematical analysis and numerical simulations, we show that the Kleinberg small-world network shares similar properties as NW small-world networks but Kleinberg small-world network is more general. Namely, we see that synchronizability of two-dimensional Kleinberg small-world network is enhanced as the edge-adding probability increases, and the average path length of the Kleinberg small-world network decreases with the increasing edge-adding probability. Moreover, larger distance among neighbors and more edges to be added could play positive roles in enhancing the synchronizability of the Kleinberg model. A network of Lorentz oscillators is taken to make numerical simulations in order to verify the observed phenomena.
Acknowledgments
The authors thank the referees and the editor for their valuable comments on this paper. This work is supported by the Natural Science Foundation of China (Grant no. 61273220), Guangdong Education University Industry Cooperation Projects (Grant no. 2009B090300355) and the Shenzhen Basic Research Project (JC201006010743A, JCYJ20120613105730482), and 2011 Foundation for Distinguished Young Talents in Higher Education of Guangdong (LYM11115).
References S. H. Strogatz and I. Stewart, “Coupled oscillators and biological synchronization,” Scientific American, vol. 269, no. 6, pp. 102–109, 1993. View at Google Scholar · View at ScopusC. M. Gray, “Synchronous oscillations in neuronal systems: mechanisms and functions,” Journal of Computational Neuroscience, vol. 1, no. 1-2, pp. 11–38, 1994. View at Publisher · View at Google Scholar · View at ScopusL. Glass, “Synchronization and rhythmic processes in physiology,” Nature, vol. 410, no. 6825, pp. 277–284, 2001. View at Publisher · View at Google Scholar · View at ScopusM. S. de Vieira, “Chaos and synchronized chaos in an earthquake model,” Physical Review Letters, vol. 82, no. 1, pp. 201–204, 1999. View at Google Scholar · View at ScopusS. H. Wang, J. Y. Kuang, J. H. Li, Y. L. Luo, H. P. Lu, and G. Hu, “Chaos-based secure communications in a large community,” Physical Review E, vol. 66, no. 6, Article ID 065202, 2002. View at Google ScholarD. Yu, M. Righero, and L. Kocarev, “Estimating topology of networks,” Physical Review Letters, vol. 97, no. 18, Article ID 188701, 2006. View at Publisher · View at Google Scholar · View at ScopusS. Milgram, “The small world problem,” in Psychology Today, vol. 2, pp. 60–67, Sussex Publishers, 1967. View at Google ScholarD. J. Watts and S. H. Strogatz, “Collective dynamics of small-world networks,” Nature, vol. 393, no. 6684, pp. 440–442, 1998. View at Google Scholar · View at ScopusM. E. J. Newman and D. J. Watts, “Renormalization group analysis of the small-world network model,” Physics Letters A, vol. 263, no. 4–6, pp. 341–346, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNetX. F. Wang and G. Chen, “Synchronization in small-world dynamical networks,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 12, no. 1, pp. 187–192, 2002. View at Publisher · View at Google Scholar · View at ScopusF. Qi, Z. Hou, and H. Xin, “Ordering chaos by random shortcuts,” Physical Review Letters, vol. 91, no. 6, Article ID 064102, 2003. View at Google Scholar · View at ScopusV. N. Belykh, I. V. Belykh, and M. Hasler, “Connection graph stability method for synchronized coupled chaotic systems,” Physica D, vol. 195, no. 1-2, pp. 159–187, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNetI. 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 · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNetL. Tang, J. Lu, and G. Chen, “Synchronizability of small-world networks generated from ring networks with equal-distance edge additions,” Chaos, vol. 22, no. 2, Article ID 023121. View at Publisher · View at Google ScholarT. Nishikawa, A. E. Motter, Y.-C. Lai, and F. C. Hoppensteadt, “Heterogeneity in oscillator networks: are smaller worlds easier to synchronize?” Physical Review Letters, vol. 91, no. 1, Article ID 014101, 2003. View at Google Scholar · View at ScopusM. Zhao, T. Zhou, B.-H. Wang, G. Yan, H.-J. Yang, and W.-J. Bai, “Relations between average distance, heterogeneity and network synchronizability,” Physica A, vol. 371, no. 2, pp. 773–780, 2006. View at Publisher · View at Google Scholar · View at ScopusX. F. Wang, X. Li, and G. Chen, Complex Network and Application, Tsinghua University Press, Beijing, China, 2006. X. F. Wang, X. Li, and G. Chen, Network Science: An Introduction, Higher Education Press, Beijing, China, 2012. J. Kleinberg, “The small-world phenomenon: an algorithm perspective,” in Proceedings of the 32nd Annual ACM Symposium on Theory of Computing, 2000. View at Publisher · View at Google Scholar · View at MathSciNetY. Zhao, J. W. Feng, and J. Y. Wang, “On synchronizability of kleinberg small world networks,” in Proceedings of 8th International Conference on Computational Intelligence and Security, pp. 204–208, 2012. T. Nishikawa and A. E. Motter, “Maximum performance at minimum cost in network synchronization,” Physica D, vol. 224, no. 1-2, pp. 77–89, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet