Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Node-node distance distribution for growing networks

Node-node distance distribution for growing networks We present the simulation of the time evolution of the distance matrix. The result is the node-node distance distribution for various kinds of net- works. For the exponential trees, analytical formulas are derived for the moments of the distance distribution. 1 Introduction A graph is defined as a set of nodes (vertices) and a set of links between nodes (edges) [1, 2, 3, 4]. By graph evolution or growth we mean subsequent attach- ing of new nodes with m edges to previously existing nodes [5]. Such growing graphs may reflect some features of real evolving networks, e.g. a network of collaborators, a network of citations of scientific papers, some biological net- works (food chains or sexual relations) or Internet and world-wide-web pages with links between them [5, 6, 7, 8, 9]. The distance between nodes is the shortest number of edges which leads from one node to the other. The node-node distance (NND) distribution depends on how subsequent nodes are attached. If each node is connected with only one of preexisting nodes (m = 1) a tree appears. When a new node is attached to several different nodes with m > 1 edges, the growing structure is called a simple graph. We may choose nodes to which new nodes are attached in preferential way or randomly. In the latter case we deal with exponential networks. If the probability of choosing a node is proportional to its degree (e.g. to the number of its nearest neighbors) the growing structure is called scale-free or Albert– Baraba ´si networks [10]. In this paper, the numerical algorithm for the network growth — basing on distance matrix evolution — is presented both, for exponential and scale-free networks (m = 1, 2) [11, 12]. The NND distribution and its characteristics are calculated. For the exponential trees the iterative formulas for n-th ordinary moments are derived. 2 Computer simulations A graph with edges of unit length may be fully characterized by its distance matrix S, an element s of which is equal to the shortest path between nodes ij presented at the 37-th Polish Physicists’ Meeting, Gdan ´sk, Poland, 15-19 Sep. 2003 arXiv:cond-mat/0309255v2 [cond-mat.stat-mech] 24 Feb 2004 i and j. This matrix representation is also particularly useful when computer simulations for graph evolution are applied. Attaching subsequent node with one edge (m = 1) to previously existing network of N nodes corresponds to adding a new (N + 1)-th row and a new column to N × N large distance matrix S. The distance from newly added (N + 1)-th node to all others via selected — labeled as p — node is larger by one than the distance between p-th and all others. Thus, new (N + 1)-th row/column is a simple copy of p-th row/column but with all of its elements incremented [11]: ∀ 1 ≤ i ≤ N : s (N + 1, i) = s (i, N + 1) = s (p, i) + 1. (1) N+1 N+1 N Similarly, when new node is attached to the network with two edges (m = 2) to two different nodes — which are labeled as p and q — but the distance from all other nodes i to the newly added (N + 1)-th is one plus the smaller distance between p–i or q–i nodes pairs [12]: ∀ 1 ≤ i ≤ N : s (N +1, i) = s (i, N +1) = min s (p, i), s (q, i) +1. (2) N+1 N+1 N N In the case mentioned above of the growth of the simple graphs also the reevaluation of distances between nodes i and j must be done to check if adding a new node provides the shortcut [12]: ∀ 1 ≤ i, j ≤ N : s (i, j) = N+1 (3) min s (i, j), s (i, p) + 2 + s (q, j), s (i, q) + 2 + s (p, j) . N N N N N In both cases diagonal elements of new row/column are zero [11, 12]: s (N + 1, N + 1) = 0. (4) N+1 Selecting rows/columns (nodes to which we attempt to add a new node) may be random or preferential. In the latter case an additional evolving vector is in- troduced, which contains the node labels. These labels occur as vector elements with a probability proportional to the degree of the node. Random selection of elements of such a vector correspond to Albert–Baraba ´si construction rule. The procedure is known as the Kert´ esz algorithm [13]. 3 Analytical calculations Let us define n-th moments of the NND distribution for all distances N N XX n n n ℓ ≡ [{s (i, j)}] = [s (i, j)], (5) i=1 j=1 and only for non-zero distances N N XX n n n d ≡ [hs (i, j)i] = [s (i, j)], (6) N (N − 1) i=1 j=1 j6=i where {· · · }, h· · · i and [· · · ] denote an average over N matrix elements, an average over N (N − 1) non-diagonal matrix elements, and an average over N run 2 Table 1: The mean distance d(N ) = a ln N + b for different evolving networks. exponential exponential scale-free scale-free m 1 2 1 2 a 2.00 0.672 1.00 0.462 b −2.84 0.296 −0.08 0.889 Table 2: The dispersion σ (N ) = c ln N +d for the exponential and the scale-free trees (m = 1). exponential scale-free c 2.00 1.00 d −1.44 −1.64 independent realizations of the evolution process (matrices), respectively. Mo- ments (5) and (6) for n = 1 are sometimes called the network diameter. Both double sums in r.h.s. of Eqs. (5) and (6) are equal, due to obvious fact, that s(i, i) = 0. That allows to derive simple dependence between averages {· · · } and h· · · i: n n Nℓ = (N − 1)d . (7) N N For the exponential trees — basing on s(i, i) = 0 and distance matrix symmetry s(i, j) = s(j, i) — we are able to construct iterative equations for ℓ as N+1 dependent on ℓ (k = 1, · · · , n) N+1 N+1 N X X X 2 n n 2 n (N + 1) ℓ = [s (i, j)] = N ℓ + 2 1 + [s(i, q)] , (8) N+1 N i=1 j=1 i=1 where q is the number of the randomly selected row/column of the distance matrix S. Combining Eq. (8) with Eq. (7) gives the desired iterative formula: n−1 (N + 2)(N − 1) 2 2(N − 1) n n n k d = d + + d . (9) N+1 N N (N + 1)N N + 1 (N + 1)N k k=1 4 Results and conclusions 1 2 2 1 2 For the trees the mean of the NND d and its dispersion σ ≡ d − (d ) N N N grow logarithmically with N (see Tabs. 1, 2) [11]. For the graphs only the first cumulant (the average of the NND d ) grows logarithmically (see Tab. 2) [12]. Such a slow increase of d with number of network nodes is known as the small-world effect [14]. The histogram of NND is presented in Fig. 1. As we expected the NND for the graphs are more condensed than the NND for the trees, as well as, the scale-free graphs (trees) are more condensed than the exponential graphs (trees). Knowing the moments d — the averages of the n-th powers of the non- diagonal distance matrix elements (6) — allows to build all statistical parame- ters which characterize the NND distribution, e.g. the average distance d, the 3 5 4 N=1000, N =10 (m=1), N =10 (m=2) run run 0.45 exponential tree scale-free tree 0.4 Poisson distribution, <s(i,j)>=10.985 Poisson distribution, <s(i,j)>=6.788 exponential graphs, <s(i,j)>=4.928 0.35 scale-free graphs, <s(i,j)>=4.067 0.3 0.25 0.2 0.15 0.1 0.05 0 5 10 15 20 25 node-node distance s(i,j) Figure 1: The NND distribution for different types trees and graphs. exponential trees, N =10 run 1e+07 [<s(i,j)>] [<s (i,j)>] [<s (i,j)>] 1e+06 [<s (i,j)>] iterative equations 1 10 100 1000 10000 100000 1e+06 1e+07 1e+08 1e+09 1e+10 Figure 2: The main moments d (k = 1, · · · , 4) for the exponential trees given by Eq. (9) (lines) and from the direct simulations (symbols). The latter are averaged over N = 10 independent evolution process realizations. run P(s(i,j)) node-node distance distribution characteristics for exponential tress mean variation skewness kurtosis 0.1 1 10 100 1000 10000 100000 1e+06 1e+07 1e+08 1e+09 1e+10 Figure 3: The NND distribution characteristics for the exponential trees as derived from iterative Eq. (9). distance dispersion σ , its skewness 3 2 1 1 3 d − 3d d + 2(d ) N N N N ν ≡ , and kurtosis 4 3 1 2 1 2 1 4 d − 4d d + 6d (d ) − 3(d ) N N N N N N κ ≡ . The values of such characteristics of the NND for the exponential trees obtained via Eq. (9) are presented in Figs. 2 and 3. For trees, the distributions are similar to the Poisson distribution (see Fig. 1). However, even for large N the skewness and kurtosis do not vanish as one may expect for the normal distribution [15]. Acknowledgments K.M.’s participation in the 37-th Polish Physicists’ Meeting in Gdan ´sk was financed by Krak´ ow Branch of the Polish Physical Society (PTF). The numerical calculations were carried out in ACK–CYFRONET–AGH. The machine time on SGI 2800 is financed by the State Committee for Scientific Research (KBN) in Poland with grant No. KBN/SGI2800/AGH/018/2003. References [1] J.Clark, D.A.Holton, A First Look at Graph Theory, World Scientific, 1991. [2] R.J.Wilson, Introduction to Graph Theory, Longman Scientific and Tech- nical, New York, 1987. 5 [3] C.Berge, Graphs, North Holland, 1985. [4] F.Harary, Graph Theory, Addison-Wesley, 1969. [5] S.N.Dorogovtsev, J.F.F.Mendes, Adv. Phys. 51 (2002) 1079. [6] S.N.Dorogovtsev, J.F.F.Mendes, in From the Genome to the Internet, eds. S.Bornholdt and H.G.Schuster, Viley-VCH, Berlin, 2002. [7] S.N.Dorogovtsev, J.F.F.Mendes, A.N.Samukhin, Nucl. Phys. B653 (2003) [8] M.E.J.Newman, SIAM Review 45 (2003) 167. [9] R.Albert, A.-L.Barab´ asi, Rev. Mod. Phys. 286 (2002) 47. [10] A.-L.Barab´ asi, R.Albert, Science 286 (1999) 509. [11] K.Malarz, J.Czaplicki, B.Kawecka-Magiera, K.Kula kowski, Int. J. Mod. Phys. C14 (2003) 1201. [12] K.Malarz, K.Kula kowski, cond-mat/0304693. [13] D.Stauffer, private communication. [14] S.Milgram, Psychol. Today 2 (1967) 60. [15] G.Szab´ o, M.Alava, J.Kert´ esz, Phys. Rev. E66 (2002) 026101. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Condensed Matter arXiv (Cornell University)

Node-node distance distribution for growing networks

Malarz, K.; Karpinska, J.; Kardas, A.; Kulakowski, K.
Condensed Matter , Volume 2021 (cond-mat/0309255) – Sep 10, 2003
Download PDF
6 pages

Loading next page...
 
/lp/arxiv-cornell-university/node-node-distance-distribution-for-growing-networks-fiXz9894nt

References

References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.

ISSN
0096-3003
eISSN
ARCH-3331
DOI
10.1016/j.amc.2012.03.091
Publisher site
See Article on Publisher Site

Abstract

We present the simulation of the time evolution of the distance matrix. The result is the node-node distance distribution for various kinds of net- works. For the exponential trees, analytical formulas are derived for the moments of the distance distribution. 1 Introduction A graph is defined as a set of nodes (vertices) and a set of links between nodes (edges) [1, 2, 3, 4]. By graph evolution or growth we mean subsequent attach- ing of new nodes with m edges to previously existing nodes [5]. Such growing graphs may reflect some features of real evolving networks, e.g. a network of collaborators, a network of citations of scientific papers, some biological net- works (food chains or sexual relations) or Internet and world-wide-web pages with links between them [5, 6, 7, 8, 9]. The distance between nodes is the shortest number of edges which leads from one node to the other. The node-node distance (NND) distribution depends on how subsequent nodes are attached. If each node is connected with only one of preexisting nodes (m = 1) a tree appears. When a new node is attached to several different nodes with m > 1 edges, the growing structure is called a simple graph. We may choose nodes to which new nodes are attached in preferential way or randomly. In the latter case we deal with exponential networks. If the probability of choosing a node is proportional to its degree (e.g. to the number of its nearest neighbors) the growing structure is called scale-free or Albert– Baraba ´si networks [10]. In this paper, the numerical algorithm for the network growth — basing on distance matrix evolution — is presented both, for exponential and scale-free networks (m = 1, 2) [11, 12]. The NND distribution and its characteristics are calculated. For the exponential trees the iterative formulas for n-th ordinary moments are derived. 2 Computer simulations A graph with edges of unit length may be fully characterized by its distance matrix S, an element s of which is equal to the shortest path between nodes ij presented at the 37-th Polish Physicists’ Meeting, Gdan ´sk, Poland, 15-19 Sep. 2003 arXiv:cond-mat/0309255v2 [cond-mat.stat-mech] 24 Feb 2004 i and j. This matrix representation is also particularly useful when computer simulations for graph evolution are applied. Attaching subsequent node with one edge (m = 1) to previously existing network of N nodes corresponds to adding a new (N + 1)-th row and a new column to N × N large distance matrix S. The distance from newly added (N + 1)-th node to all others via selected — labeled as p — node is larger by one than the distance between p-th and all others. Thus, new (N + 1)-th row/column is a simple copy of p-th row/column but with all of its elements incremented [11]: ∀ 1 ≤ i ≤ N : s (N + 1, i) = s (i, N + 1) = s (p, i) + 1. (1) N+1 N+1 N Similarly, when new node is attached to the network with two edges (m = 2) to two different nodes — which are labeled as p and q — but the distance from all other nodes i to the newly added (N + 1)-th is one plus the smaller distance between p–i or q–i nodes pairs [12]: ∀ 1 ≤ i ≤ N : s (N +1, i) = s (i, N +1) = min s (p, i), s (q, i) +1. (2) N+1 N+1 N N In the case mentioned above of the growth of the simple graphs also the reevaluation of distances between nodes i and j must be done to check if adding a new node provides the shortcut [12]: ∀ 1 ≤ i, j ≤ N : s (i, j) = N+1 (3) min s (i, j), s (i, p) + 2 + s (q, j), s (i, q) + 2 + s (p, j) . N N N N N In both cases diagonal elements of new row/column are zero [11, 12]: s (N + 1, N + 1) = 0. (4) N+1 Selecting rows/columns (nodes to which we attempt to add a new node) may be random or preferential. In the latter case an additional evolving vector is in- troduced, which contains the node labels. These labels occur as vector elements with a probability proportional to the degree of the node. Random selection of elements of such a vector correspond to Albert–Baraba ´si construction rule. The procedure is known as the Kert´ esz algorithm [13]. 3 Analytical calculations Let us define n-th moments of the NND distribution for all distances N N XX n n n ℓ ≡ [{s (i, j)}] = [s (i, j)], (5) i=1 j=1 and only for non-zero distances N N XX n n n d ≡ [hs (i, j)i] = [s (i, j)], (6) N (N − 1) i=1 j=1 j6=i where {· · · }, h· · · i and [· · · ] denote an average over N matrix elements, an average over N (N − 1) non-diagonal matrix elements, and an average over N run 2 Table 1: The mean distance d(N ) = a ln N + b for different evolving networks. exponential exponential scale-free scale-free m 1 2 1 2 a 2.00 0.672 1.00 0.462 b −2.84 0.296 −0.08 0.889 Table 2: The dispersion σ (N ) = c ln N +d for the exponential and the scale-free trees (m = 1). exponential scale-free c 2.00 1.00 d −1.44 −1.64 independent realizations of the evolution process (matrices), respectively. Mo- ments (5) and (6) for n = 1 are sometimes called the network diameter. Both double sums in r.h.s. of Eqs. (5) and (6) are equal, due to obvious fact, that s(i, i) = 0. That allows to derive simple dependence between averages {· · · } and h· · · i: n n Nℓ = (N − 1)d . (7) N N For the exponential trees — basing on s(i, i) = 0 and distance matrix symmetry s(i, j) = s(j, i) — we are able to construct iterative equations for ℓ as N+1 dependent on ℓ (k = 1, · · · , n) N+1 N+1 N X X X 2 n n 2 n (N + 1) ℓ = [s (i, j)] = N ℓ + 2 1 + [s(i, q)] , (8) N+1 N i=1 j=1 i=1 where q is the number of the randomly selected row/column of the distance matrix S. Combining Eq. (8) with Eq. (7) gives the desired iterative formula: n−1 (N + 2)(N − 1) 2 2(N − 1) n n n k d = d + + d . (9) N+1 N N (N + 1)N N + 1 (N + 1)N k k=1 4 Results and conclusions 1 2 2 1 2 For the trees the mean of the NND d and its dispersion σ ≡ d − (d ) N N N grow logarithmically with N (see Tabs. 1, 2) [11]. For the graphs only the first cumulant (the average of the NND d ) grows logarithmically (see Tab. 2) [12]. Such a slow increase of d with number of network nodes is known as the small-world effect [14]. The histogram of NND is presented in Fig. 1. As we expected the NND for the graphs are more condensed than the NND for the trees, as well as, the scale-free graphs (trees) are more condensed than the exponential graphs (trees). Knowing the moments d — the averages of the n-th powers of the non- diagonal distance matrix elements (6) — allows to build all statistical parame- ters which characterize the NND distribution, e.g. the average distance d, the 3 5 4 N=1000, N =10 (m=1), N =10 (m=2) run run 0.45 exponential tree scale-free tree 0.4 Poisson distribution, <s(i,j)>=10.985 Poisson distribution, <s(i,j)>=6.788 exponential graphs, <s(i,j)>=4.928 0.35 scale-free graphs, <s(i,j)>=4.067 0.3 0.25 0.2 0.15 0.1 0.05 0 5 10 15 20 25 node-node distance s(i,j) Figure 1: The NND distribution for different types trees and graphs. exponential trees, N =10 run 1e+07 [<s(i,j)>] [<s (i,j)>] [<s (i,j)>] 1e+06 [<s (i,j)>] iterative equations 1 10 100 1000 10000 100000 1e+06 1e+07 1e+08 1e+09 1e+10 Figure 2: The main moments d (k = 1, · · · , 4) for the exponential trees given by Eq. (9) (lines) and from the direct simulations (symbols). The latter are averaged over N = 10 independent evolution process realizations. run P(s(i,j)) node-node distance distribution characteristics for exponential tress mean variation skewness kurtosis 0.1 1 10 100 1000 10000 100000 1e+06 1e+07 1e+08 1e+09 1e+10 Figure 3: The NND distribution characteristics for the exponential trees as derived from iterative Eq. (9). distance dispersion σ , its skewness 3 2 1 1 3 d − 3d d + 2(d ) N N N N ν ≡ , and kurtosis 4 3 1 2 1 2 1 4 d − 4d d + 6d (d ) − 3(d ) N N N N N N κ ≡ . The values of such characteristics of the NND for the exponential trees obtained via Eq. (9) are presented in Figs. 2 and 3. For trees, the distributions are similar to the Poisson distribution (see Fig. 1). However, even for large N the skewness and kurtosis do not vanish as one may expect for the normal distribution [15]. Acknowledgments K.M.’s participation in the 37-th Polish Physicists’ Meeting in Gdan ´sk was financed by Krak´ ow Branch of the Polish Physical Society (PTF). The numerical calculations were carried out in ACK–CYFRONET–AGH. The machine time on SGI 2800 is financed by the State Committee for Scientific Research (KBN) in Poland with grant No. KBN/SGI2800/AGH/018/2003. References [1] J.Clark, D.A.Holton, A First Look at Graph Theory, World Scientific, 1991. [2] R.J.Wilson, Introduction to Graph Theory, Longman Scientific and Tech- nical, New York, 1987. 5 [3] C.Berge, Graphs, North Holland, 1985. [4] F.Harary, Graph Theory, Addison-Wesley, 1969. [5] S.N.Dorogovtsev, J.F.F.Mendes, Adv. Phys. 51 (2002) 1079. [6] S.N.Dorogovtsev, J.F.F.Mendes, in From the Genome to the Internet, eds. S.Bornholdt and H.G.Schuster, Viley-VCH, Berlin, 2002. [7] S.N.Dorogovtsev, J.F.F.Mendes, A.N.Samukhin, Nucl. Phys. B653 (2003) [8] M.E.J.Newman, SIAM Review 45 (2003) 167. [9] R.Albert, A.-L.Barab´ asi, Rev. Mod. Phys. 286 (2002) 47. [10] A.-L.Barab´ asi, R.Albert, Science 286 (1999) 509. [11] K.Malarz, J.Czaplicki, B.Kawecka-Magiera, K.Kula kowski, Int. J. Mod. Phys. C14 (2003) 1201. [12] K.Malarz, K.Kula kowski, cond-mat/0304693. [13] D.Stauffer, private communication. [14] S.Milgram, Psychol. Today 2 (1967) 60. [15] G.Szab´ o, M.Alava, J.Kert´ esz, Phys. Rev. E66 (2002) 026101.

Journal

Condensed Matter – arXiv (Cornell University)

Published: Sep 10, 2003

References

  • Complex networks: structure and dynamics
    Boccaletti
  • Analyzing and modeling real-world phenomena with complex networks: a survey of applications
    da Fontoura Costa