A Novel Neighbor-Preferential Growth Scale-Free Network Model and its Properties

In this paper, we propose a novel neighbor-preferential growth (NPG) network model. Theoretical analysis and numerical simulations indicate the new model can reproduce not only a scale-free degree distribution and its power exponent is related to the edge-adding number m, but also a small-world effect which has large clustering coefficient and small average path length. Interestingly, the clustering coefficient of the model is close to that of globally coupled network, and the average path length is close to that of star coupled network. Meanwhile, the synchronizability of the NPG model is much stronger than that of BA scale-free network, even stronger than that of synchronization-optimal growth network.


Introduction
In the past two decades, complex networks have been extensively studied and have gained rich research results.In particular, discoveries of the small world effect [1] and the scale-free feature [2] of complex networks have promoted the research of network structure [3]- [11].And it has been known that the topology of a network often plays an important role in the feature of dynamical network.
Since Barabási and Albert have addressed that two key mechanisms-growth and preferential attachment.BA scale-free network model [2] was proposed to promote the understanding of the real network.Subsequently, an increasing number of studies have been devoted to creating a new model [12]- [19].Bollobas and Riordan proposed another version of the scale-free model, which is called the linearized chord diagram (LCD) model [14], in the model multiple-edges and self-loops can be allowed.Xiang Li et al. proposed a local-world evolving network model [15], which represents a transition between power-law and exponential scaling.Y.J. Cao et al. proposed a neighbourhood evolving model [16], which also yields a transition between power-law and exponential.
In order to enhance the synchronizability of the growth network, X.F.Wang et al. have constructed the synchronization-optimal growing network [12], it is easier to achieve synchronization than BA scale-free networks, however, it is more particularly fragile with respect to deliberate attacks.Hence, for improving the synchronization robustness against the selective removal of vertices in complex networks, X.F.Wang et al. have also proposed the synchronous preference of complex dynamical networks [13], and its synchronizability is investigated and found to be robust against both random and specific removal of nodes.
There are many scale-free growth network models [17] [19], all of which are obtained by changing the preferential attachment mechanism.
In the real-world networks, there is another priority connection mechanism, "neighbor" effect.For example, when you enter a strange city to visit your friends, you will just make friends with the friends of your present friends.Based on the growth and neighbor preferential principles, we propose a dynamical networks model, which is called the NPG model in this paper.To further clarify the topology of the NPG network, we confirm the statistical properties of the new network in three aspects-the degree probability distribution, the clustering coefficient and the average path length.In fact, the three concepts play an important role in the development of complex network theory [20] [21].Theoretical analysis and numerical simulations indicate the new network is not only a scale-free network whose degree distribution obeys the power-law and its power exponent is related to the edge-adding number m, but also a small-world network which has large clustering coefficient and small average path distance.
Meanwhile, its synchronizability is much stronger than that of BA scale-free network [2], even stronger than that of synchronization-optimal network [12].
And the NPG network is robust with respect to random attacks and is fragile to specific removal of a small fraction of nodes.

The Construction of the NPG Model
Our model is defined as follows: 1) Initial network: Starting with a globally coupled network with 0 m nodes.
2) Growth: At every time step a new vertex is added and is connected to m ( 0 m m ≤ ) nodes that already-existing in the system.
3) Preferential attachment: a) The new node is connected to an existing node i only for one edge.The connection probability i Π depends on the degree of the existing node: b) Select the first 1 m − nodes successively in the order of descending degree (when the degree of the nodes is the same, first select the older node) from all the neighbor nodes of the node i which is connected to the new node.Add

Degree Distribution
The degree distribution as one of the most important statistical features of complex networks.We calculate the degree probability distribution ( ) P k , using the mean-field approach [22].By the connection mechanism of the algorithm of the model, we know that the model has 1 m − "hubs" which are connected to all nodes and a big node which is connected to most nodes, as shown in Figure 1(a) (in the later simulation we also see the phenomenon ).For convenience, the big node is also called the "hub".And there is no doubt that the m "hubs" are connected to each other.
Suppose the initial network has 0 m nodes and the degree of the node i is ( ) i k t at time t.When t is large enough, we can ignore the number of edges in the initial network and 0 m t t + ≈ .Divided the nodes into two categories: the m "hubs" and the common nodes.When a new node is added to the network, the probability of the new node is connected to the node i is ( ) i Π .If node i belongs to the "hubs", the ( ) 1 i Π = , but we will ignore the node i in the degree probability distribution, due to the number of the "hubs" is much smaller than the total number of nodes in the network ( m t  ).Then we only need to consider when the node i belongs to the common nodes, the probability of the new node is connected to the node i is i Π in the preferential attachment (i), if the new node is not connected to node i, the preferential attachment (ii) will work, but the new node is also not connected to node i, because it will be connected to the "hubs".So when the node i belongs to the common nodes, only the preferential attachment (i) works.In a word, the connection probability of node i is ( ) ( ) ( ) ( ) Assume the degree of the node i is continuous, and thus the connection probability can be interpreted as the rate of change of the ith node's degree, we can write ( ) ( ) The solution of this equation, with the initial condition that node i was added to the system at time i t with connectivity ( ) the probability that a node has a connectivity ( ) Suppose we add nodes at equal time intervals to the network, the probability By ( 5) and ( 6), then the probability density for ( ) That means, ( ) 2 1 ~m P k k − − , the degree probability distribution of the NPG model follows power-law distribution, and the power exponent is 2 1 m − − , not fixed.When m is larger, the smaller the power exponent, the steeper the corresponding graph.In Figure 2(a) we can see that when 3, 5, 7 m = , the lines will become more and more steep, supporting the analytical result.In particular, there is a limiting case-when 1 m = , the preferential attachment (ii) is not effec- tive in the NPG model.Namely, in this limiting case the NPG model is the same as BA model, and Eq. ( 7) is valid for the limiting case.In Figure 2(b), we can also see when 1 m = , the graphs of BA and the NPG are almost coincident, supporting the analytical result.

Average Path Length
The average path length is also an important statistical characteristics for studying complex networks, especially when describing the synchronizability of networks.The shorter the average distance length, the stronger the synchronizability [23] [24] [25], Generally.For the NPG network, the nodes are divided into three sets.Set one: the 1 m − "hubs" which are connected to all nodes; Set two: a big node which is connected to most nodes; Set three: the common nodes.As shown in Figure 1(a).The path length can be expressed mathematically as follows: , respectively.
where ij S is the sum of the shortest path lengths of all nodes in set i to all nodes in set j, where x k is the degree of the big node.So the average path length of the NPG network is S When N is large enough, the average length is two.From Figure 3(b), we can get that the numerical results are consistent with the theoretical predictions.In Figure 3(b), we also show that the comparison results of the average path length for the NPG model and BA model.We find the average path length of the NPG network is significantly shorter than that of the BA network.

The Clustering Coefficient
The small-world characteristic consists of two properties: large clustering coefficient and small average path length.We have demonstrated that the NPG network has a small average path length.The clustering coefficient is C where i C is the local clustering coefficient for node i. ( ) where ( ) E i is the number of edges in the neighbor set of the node i, and i k is the degree of node i.First of all, we simply analyze the structure of the NPG network (see Figure 1): 1) When the degree of node i is equal to the edge-adding number m ( i k m = ), the clustering coefficient of node i is one.
2) Get rid of the "hubs" from the neighbor set of any node i which exist in Set 2 and Set 3 to form a new neighbor set in which edges is not exists between nodes.
Simply prove ii), if there is a connection between node i and j in this new neighbor set, then node i or node j with 1 m + edges into the network system, which is impossible.Namely, conclusion ii) is established.The local clustering coefficients of the nodes in Set 2 and Set 3 are i C The local clustering coefficients of the nodes in Set 1 are i C′ , according to the degree of the nodes in the network, it is easy to prove the following equation The degree i k of the "hubs" satisfies ′ ≈ .Finally, for the NPG network, the clustering coefficient of node i is 2 1 by (21),it can be seen when the degree of the node becomes larger, the local clustering scales as ( ) , where

( )
C k denotes the average clustering values of nodes with degree k.In Figure 4, this phenomenon is also seen.When calculating the clustering coefficient C, we only consider the common nodes.We have ( ) ( ) l is a constant and l N  ，since the degree distribution is ( ) , the average clustering coefficient C can be rewritten as The numerical results are demonstrated in Figure 3(a), we can get that a little fluctuations between the numerical results and analysis results, which may  attribute to the fluctuations of the power-law exponent of degree distribution.
But it does not affect the network with a large clustering coefficient.In Figure 3(a), we show the comparison result of the clustering coefficient for BA model and the NPG model.We find that the clustering coefficient of the proposed model is larger than that of BA model obviously.And the clustering coefficient of the new model will increase as m increases.
The above simulation results show that the new network is a scale-free network and its power exponent is related to the edge-adding number m.What's more, it has a large clustering coefficient and a short average path length, which represent a small-world effect.So the model has the characterization of both scale-free and small-world networks.

Synchronization Stability Criterion of Complex Dynamical Networks
Consider a complex dynamical network consisting of N identical linearly and diffusively coupled vertices, with each vertex being an n-dimensional dynamical oscillator.The state equation of the network can be written as [25] [26] [27]: ( ) where , , , are the state variables of vertex i and the constant coupling strength c is assumed positive.Furthermore, , , , there is a connection between vertex i and vertex j ( ) where the degree i k of vertex i is defined to be the number of its outreaching connections.Defined that dynamical network ( 24) is (asymptotically) synchro- where ( ) s t is a solution of an isolated vertex [28].Suppose that the network is connected in the sense that there are no isolated clusters.Then, a symmetric and irreducible matrix, whose eigenvalues are . It has been proved that under certain assumptions the synchronized state ( 26) is exponentially stable if where d is a constant, which was further specified to be max h , the maximal Lya- punov exponent of an individual n-dimensional chaotic dynamical system [27].
Given the dynamics of an isolated vertex, the synchronizability of network (24) with respect to a specific coupling configuration A is said to be strong if the network can synchronize with a small coupling strength c.Inequality ( 27) implies that the synchronizability of network ( 24) can be characterized by the second-largest eigenvalue of its coupling matrix, i.e., the smaller the secondlargest eigenvalue, the stronger the synchronizability of a network.respectively [12].This indicates that the synchronizability of the NPG network is stronger than that of the synchronization-optimal growing network.What's more, the time complexity of the algorithm of the NPG network model is much lower than that of synchronization-optimal growing network model.

Synchronizability of the NPG Network
We find the second-largest eigenvalue 2np ) of the multi-center growth model [12].Because the topology of the NPG network is so similar to that of the multi-center network, both of them have several "hubs".
The second-largest eigenvalues of dynamical network with BA scale-free topology, synchronization-optimal growth topology, NPG topology and multicenter topology are shown in the Table 1.
Table 1.A comparision of the second-largest eigenvalues of different model.The new second-largest eigenvalue

Robustness and Fragility
Clearly, when we remove some nodes in a network will change its coupling matrix.If the second-largest eigenvalue of the coupling matrix remains unchanged, namely, the synchronizability of the network will remain unchanged after the removal of a small fraction of nodes, which indicates the robustness of the synchronizability of the network.If the second-largest eigenvalue of the coupling matrix changes greatly after the removal of some its nodes, which implies the fragility of the synchronizability of the network.Now we have two ways to remove a small fraction f ( 0 1 f <  ) of nodes in the work: random or specific.

We let
be the coupling matrices of the original network with N nodes and the new network after removal of fN     nodes, respectively.The second-largest eigenvalues of A and A  are denoted by 2 λ and 2 λ  , respectively.
In the simulation, we take 1000 N = and 0 3 m m = = . It has been shown that even when as many as 5% of randomly chosen nodes are removed from a BA scale-free network, the new second-largest eigenvalue 2ba λ  will not be changed too much (Figure 6(a)).And the new second-largest eigenvalue 2np λ  of the NPG network remains unchanged when 5% of vertices are randomly removed.This indicates that both the BA scale-free model and the NPG model are robust against random attacks, the robustness of the NPG model is stronger, especially.However, the NPG model is more fragile to specific removal of a small fraction of nodes than a BA network (Figure 6(b)), due to the topology of the NPG network is very similar to that of the multi-center network, and the whole network is broken into isolate when several "biggest" vertices are removed.

Conclusion
In this paper, we have proposed the neighbor preferential mechanism and construct the neighbor-preferential growth network model.Theoretical analysis and numerical simulations show the NPG model is a scale-free small-world network model.What's more, its power exponent will increase drastically as the edgeadding number m increases.And the NPG model can reproduce star coupled network feature and globally coupled network property, for its average path length is close to two and its clustering coefficient of is close to one.Moreover, the synchronizability of the NPG network model is much stronger than that of BA scale-free network [2], even stronger than that of synchronization-optimal network [12].Due to the topology of the NPG network is so similar to that of the multi-center network, the NPG model is robust against random attacks and is more fragile to specific failures.Perhaps this model helps to future understand the real-world networks.However, there are still many problems need to be further studied, for example, how to adjust the number of the "hubs"?How does the clustering coefficient change with the number of the "hubs"?How many the "hubs" in the model are more in line with the real-world networks?And so on.
There questions are waiting for us to investigate.

Figure 1 .
Figure 1.The simplified diagram of the NPG model with 3 m = .

Figure 2 .
Figure 2. The degree probability distribution of the network ( 5000 N = ) in the log-log scale.(a) the NPG network with 3,5,7 m = , respectively.(b) the NPG and BA network with 1,3 m =, respectively.

Figure 3 .
Figure 3. (a) the comparison of clustering coefficient for the NPG and BA model with 3,5,7 m = , respectively.(b) the comparison of the average path length for the NPG model and BA model with 3,5,7 m =, respectively.

Figure 4 .
Figure 4.The average clustering coefficient of nodes with degree k when 3,5,7 m =, respectively.The data are averaged over 20 independent runs of network size of 5000 N = .

ForFigure 5 .3
Figure 5.The second-largest eigenvalue of the coupling matrix of the BA (green) and the NPG(red) network for 0 3,5,7 m m = = , respectively.

Figure 6 .
Figure 6.The new second-largest eigenvalue of the BA scale-free networks (red solid line with circle), the NPG networks(blue solid line with stars)when a fraction f of vertices are removed.(a) randomly selected vertices.(b) the most connected vertices.