Exact Solutions of a Generalized Weighted Scale Free Network

We investigate a class of generalized weighted scale-free networks, where the new vertex connects to m pairs of vertices selected preferentially. The key contribution of this paper is that, from the standpoint of random processes, we provide rigorous analytic solutions for the steady state distributions, including the vertex degree distribution, the vertex strength distribution and the edge weight distribution. Numerical simulations indicate that this network model yields three power law distributions for the vertex degrees, vertex strengths and edge weights, respectively.


Introduction
Recently, complex networks have become a hot topic due to the two most important discoveries of real networks, that is, the small world phenomenon raised by Watts and Strogatz [1], and the scale-free property presented by Barabási and Albert [2]. For most real-world complex systems, the proportion of vertices with degree obeys a power-law form, that is, ( ) ∝ − , where is a constant (see Albert et al. [3] for detail).
To describe such a property, Barabási et al. [4] proposed a model (BA model for short) with growth and preferential attachment. The mean field method was used in [4] to get approximate expressions of the degree distribution. Since then, Krapivsky et al. [5] expanded the preferential attachment of the BA model to be nonlinear. With the rate equation method, it is found that power-law degree distribution can only be obtained in the case of linear preferential attachment. Meanwhile, Dorogovtsev et al. [6] generalized the BA model to include initial attractiveness. Except for preferential selection of vertices, Dorogovtsev et al. [7] introduced a simple model based on selection of edges. This model starts with three vertices connected to each other; at each time step, an edge is randomly selected, and both its ends are connected to the new vertex. With the master equation method, the degree distribution of this model was shown to behave as a power-law form [7].
All the models of [4][5][6][7] share a common property that all links are equivalent. However, it is widely known that interaction strengths can vary widely, and such variations are essential to the network's ability to carry on its basic functions [8]. Barrat et al. [9] proposed a model for the growth of weighted networks that couples the establishment of new edges and vertices, and the weights' dynamical evolution. Tanaka and Aoyagi [10] generalized the model of [9] through the weight-driven preferential attachment of new vertices to existing vertices and the growth of the weights of existing links, generating scale-free networks with variable power-law exponents.
As one of the most important characteristics of networks, the degree distribution is always a big consideration. However, there is no systematic and rigorous method for solving degree distribution, the methods in [4][5][6][7] have to assume the existence of ( ) or assume the vertex number ( ) with degree be continuous. But Bollobás et al. [11] gave a rigorous method, but this only applies to networks with multiple edges and loops. Stem from this consideration, in this paper, we try to provide a systematic and rigorous method to solve the steady state distribution from a new perspective.
Firstly, we shall give a new weighted network model. The model starts with a complete graph, at each step, the new vertex is connected to ( ≥ 1) pairs of vertices which are selected preferentially; and meantime, the weights between the selected pairs of vertices are strengthened. Based on the Stolz theorem, rigorous analytic derivations to the steady state degree distribution of the model is provided. What's more, the method of solving the degree distribution can be applied to study the steady state of both the vertex strength distribution and the edge weight distribution. Both theoretical derivations and numerical simulations show that the model displays power-law behavior for the degree, the strength and the weight distribution. Lastly, we show that the clustering coefficient is larger than that of the BA model.
The rest of the paper is organized as follows. Section 2 is the description of the model, and Section 3 is devoted to discuss some properties of the network. The analytical results are demonstrated by numerical simulations in the last section.

The Model
The detailed model construction algorithm is described as follows: (i) Initialization. The initial network, denoted as 0 , is a complete graph with vertices, which are denoted by 01 , 02 , . . . , 0 . The weight of each link is assigned To construct +1 from , we add a new vertex +1 and then add 2 (2 ≤ − 1) edges between +1 and vertices of . We choose pairs of vertices 1 , 1 , 2 , 2 , . . . , , according to a preferential attachment rule. If ( ) is the weight between vertices and in , then the vertices and are selected with probability independently for each vertex pair , ∈ 1, 2, . . . , . Note that this allows the possibility that we choose the same vertex more than once, and hence the graphs may have multiple edges.
The weight of each new edge is 0 . Meanwhile, the weight between each selected pair of vertices is strengthened by adding 0 to it. If there are multiple edges added between two vertices, then the weight between them will increase multiple times. For example, if there are new edges added between vertices and , then the weight between them will increase 0 . In this case, the strength of vertex can be expressed as The number of vertices in is obviously = + . Because 2 edges are added at each stage, the total number of edges in is = ( − 1)/2 + 2 . Moreover, the total strength of vertices in is = ( − 1) 0 + 6 0 . Before analyzing the steady state distribution of the network, let us firstly introduce some concepts and symbols. Let the random variable ( ) be the number of vertices with degree in , and moreover let the network degree distribution be the average over all its vertices at time , namely ( ( ))/ . If lim → ∞ ( ( ( ))/ ) exists, we say that the steady state degree distribution of the network exists and it can be quantified by ( ) = lim → ∞ ( ( ( ))/ ).
Let ( ) represent the number of vertices with strength in Similarly, the steady state distribution of vertex strength is denoted as ( ) = lim → ∞ ( ( ( ))/ ).

Network Analysis.
To analytically obtain the statistical properties of the network generated by the above algorithm, we give strict proofs from the perspective of random process.
Proof. For use in this proof, as well as later proofs, we define the -algebra F = ( ; 0 ≤ ≤ ).
For 0 ≤ ≤ , the probability that an old vertex , with degree and strength , to get edges from the new vertex +1 is Clearly, the above equation can also be rewritten as Note that for vertices 01 , 02 , . . . , 0 , the initial strength is ( − 1) 0 . If the degree of 0 ( = 1, 2, . . . , ) at time is 0 ( ), then it means that there are 0 ( ) − ( − 1) new edges linked to vertex 0 in the process of network growth, thus the strength of 0 at time is Journal of Applied Mathematics 3 Similarly, for vertices 1 to , we have where ( ) is the degree of vertex in . The above two cases can be rewritten together as For variable , ( ) is a random process. Thus, given F , the expected number of vertices with degree in +1 is In particular, for = 2 , we have Taking expectations on both sides and combing the property of conditional probability, Computing by Lemma 1 iteratively, we get Divide by on both sides, let then , > 0, +1 − > 0 and { } is a monotone increasing sequence with → ∞, lim → ∞ (( +1 − )/ ( +1 − )) = 3/(2 + 3), by Lemma 2, we have lim → ∞ ( / ) = 3/(2 + 3). Thus, lim → ∞ ( ( 2 ( ))/ ) exists and For > 2 , Taking expectation on both sides and using the properties of conditional expectations, This completes the proof.

Remark 4.
For large , the degree distribution is proportional to −5/2 , that is, ( ) ∝ −5/2 . This implies that the degree distribution of our model obeys a power-law form like the BA model but the degree exponent is no longer 3.
In order to obtain the vertex strength distribution analytically, we derive it in the same way as the degree distribution Proof. Following the arguments of Theorem 3, given F , the expected number of vertices with strength in +1 is Use the same procedure of Theorem 3, we know that lim → ∞ ( ( 2 0 ( ))/ ) exists and For > 2 0 , suppose lim → ∞ ( ( − 2 0 ( ))/ ) exists and equals ( −2 0 ), then lim → ∞ ( ( ( ))/ ) also exists and This, together with (22), leads to Set 0 = 1, the exact expression of vertex strength distribution has the form of as required. This completes the proof.
In what follows, we shall discuss the edge weight distribution.

Theorem 6 (Weight Distribution). The steady state edge weight distribution of the model exists and is given by
Proof. Since the proof is similar to that of Theorems 3 and 5, here we omit it.
Remark 7. By Theorem 6, we can conclude that the edge weight distribution follows power-law behavior with degree exponent = 3 + 1, which is related to .

Numerical Simulations
In this section, we perform several numerical simulations in regard to the vertex degree distribution, the vertex strength   results predicted by (3). Moreover, the legends here are the same as those in Figures 2 and 3. Figure 2 shows the results of the simulations and the analysis of the vertex strength distribution ( ) in log-log scale. Because of the linear relation between the degree and the strength of the vertex, their distributions are in the same form, that is, both obey a power-law form with degree exponent = 5/2. Figure 3 clearly verifies the theoretical results obtained by (26) with = 1. However, when it comes to = 2 and = 3, there are too few simulation data to explain the theoretical formula. All of these are resulted from the attribute of the negative exponential form of (26). With this form, the percentage of the edges dramatically decreases with the increase of and/or , and the number of steps required for more accurate simulation results increases exponentially. In fact, the number is at least (3 +1) /2 without considering the random error. For example, the minimum number of required step for simulation is 174,763, when = 4 and = 3; however, we can obtain = 24 by the same number of simulation step if = 1. Despite the fact that the data of = 2 and = 3 are so few, the simulations can explain the edge weight distribution in some degree; they can at least illustrate the trend of the distribution. It is also easy to see that has great influence on the weight distribution.
In Figure 4, we make a comparison of the simulation time under different values of . It is obvious that the time required increases dramatically with the increase of the network size.

Conclusions
A weighted network model based on preferential selection of weights between vertices is studied in this paper. Based on the Stolz theorem, rigorous and analytic proofs for the steady state distributions of the model are provided. The approach developed here is quite general, applicable to many other scale free types of complex networks. Moreover, the theoretical derivations are verified by computer simulations (see Figures  1-3 for details). Both theoretical derivations and numerical simulations show that the model displays power-law behavior for the degree, the strength and the edge weight distribution. In addition, to illustrate the complexity of generating complex networks, we compare the number of steps and the simulation time under different values of .