# Barabási–Albert model

(Redirected from BA Model)
Display of three graphs generated with the Barabasi-Albert (BA) model. Each has 20 nodes and a parameter of attachment m as specified. The color of each node is dependent upon its degree (same scale for each graph).

The Barabási–Albert (BA) model is an algorithm for generating random scale-free networks using a preferential attachment mechanism. Several natural and human-made systems, including the Internet, the world wide web, citation networks, and some social networks are thought to be approximately scale-free and certainly contain few nodes (called hubs) with unusually high degree as compared to the other nodes of the network. The BA model tries to explain the existence of such nodes in real networks. The algorithm is named for its inventors Albert-László Barabási and Réka Albert and is a special case of a more general model called Price's model.[1]

## Concepts

Many observed networks (at least approximately) fall into the class of scale-free networks, meaning that they have power-law (or scale-free) degree distributions, while random graph models such as the Erdős–Rényi (ER) model and the Watts–Strogatz (WS) model do not exhibit power laws. The Barabási–Albert model is one of several proposed models that generate scale-free networks. It incorporates two important general concepts: growth and preferential attachment. Both growth and preferential attachment exist widely in real networks.

Growth means that the number of nodes in the network increases over time.

Preferential attachment means that the more connected a node is, the more likely it is to receive new links. Nodes with higher degree have stronger ability to grab links added to the network. Intuitively, the preferential attachment can be understood if we think in terms of social networks connecting people. Here a link from A to B means that person A "knows" or "is acquainted with" person B. Heavily linked nodes represent well-known people with lots of relations. When a newcomer enters the community, they are more likely to become acquainted with one of those more visible people rather than with a relative unknown. The BA model was proposed by assuming that in the World Wide Web, new pages link preferentially to hubs, i.e. very well known sites such as Google, rather than to pages that hardly anyone knows. If someone selects a new page to link to by randomly choosing an existing link, the probability of selecting a particular page would be proportional to its degree. The BA model claims that this explains the preferential attachment probability rule. However, in spite of being quite a useful model, empirical evidence suggests that the mechanism in its simplest form does not apply to the World Wide Web as shown in "Technical Comment to 'Emergence of Scaling in Random Networks'".

Later, the Bianconi–Barabási model works to address this issue by introducing a "fitness" parameter. Preferential attachment is an example of a positive feedback cycle where initially random variations (one node initially having more links or having started accumulating links earlier than another) are automatically reinforced, thus greatly magnifying differences. This is also sometimes called the Matthew effect, "the rich get richer". See also autocatalysis.

## Algorithm

The steps of the growth of the network according to the Barabasi–Albert model (${\displaystyle m_{0}=m=2}$)

The network begins with an initial connected network of ${\displaystyle m_{0}}$ nodes.

New nodes are added to the network one at a time. Each new node is connected to ${\displaystyle m\leq m_{0}}$ existing nodes with a probability that is proportional to the number of links that the existing nodes already have. Formally, the probability ${\displaystyle p_{i}}$ that the new node is connected to node ${\displaystyle i}$ is[2]

${\displaystyle p_{i}={\frac {k_{i}}{\sum _{j}k_{j}}},}$

where ${\displaystyle k_{i}}$ is the degree of node ${\displaystyle i}$ and the sum is made over all pre-existing nodes ${\displaystyle j}$ (i.e. the denominator results in twice the current number of edges in the network). Heavily linked nodes ("hubs") tend to quickly accumulate even more links, while nodes with only a few links are unlikely to be chosen as the destination for a new link. The new nodes have a "preference" to attach themselves to the already heavily linked nodes.

A network generated according to the Barabasi Albert model. The network is made of 50 vertices with initial degrees ${\displaystyle m_{0}=1}$.

## Properties

### Degree distribution

The degree distribution of the BA Model, which follows a power law. In loglog scale the power law function is a straight line.[3]

The degree distribution resulting from the BA model is scale free, in particular, it is a power law of the form

${\displaystyle P(k)\sim k^{-3}\,}$

### Hirsch index distribution

The h-index or Hirsch index distribution was shown to also be scale free and was proposed as the lobby index, to be used as a centrality measure[4]

${\displaystyle H(k)\sim k^{-6}\,}$

Furthermore, an analytic result for the density of nodes with h-index 1 can be obtained in the case where ${\displaystyle m_{0}=1}$

${\displaystyle H(1){\Big |}_{m_{0}=1}=4-\pi \,}$

### Average path length

The average path length of the BA model increases approximately logarithmically with the size of the network. The actual form has a double logarithmic correction and goes as[5]

${\displaystyle \ell \sim {\frac {\ln N}{\ln \ln N}}.}$

The BA model has a systematically shorter average path length than a random graph.

### Node degree correlations

Correlations between the degrees of connected nodes develop spontaneously in the BA model because of the way the network evolves. The probability, ${\displaystyle n_{k\ell }}$, of finding a link that connects a node of degree ${\displaystyle k}$ to an ancestor node of degree ${\displaystyle \ell }$ in the BA model for the special case of ${\displaystyle m=1}$ is given by

${\displaystyle n_{k\ell }={\frac {4\left(\ell -1\right)}{k\left(k+1\right)\left(k+\ell \right)\left(k+\ell +1\right)\left(k+\ell +2\right)}}+{\frac {12\left(\ell -1\right)}{k\left(k+\ell -1\right)\left(k+\ell \right)\left(k+\ell +1\right)\left(k+\ell +2\right)}}.}$

This is certainly not the result expected if the distributions were uncorrelated, ${\displaystyle n_{k\ell }=k^{-3}\ell ^{-3}}$.[2]

For general ${\displaystyle m}$, the fraction of links who connect a node of degree ${\displaystyle k}$ to a node of degree ${\displaystyle \ell }$ is[6]

${\displaystyle p(k,\ell )={\frac {2m(m+1)}{k(k+1)\ell (\ell +1)}}\left[1-{\frac {{\binom {2m+2}{m+1}}{\binom {k+\ell -2m}{\ell -m}}}{\binom {k+\ell +2}{\ell +1}}}\right].}$

Also, the nearest-neighbor degree distribution ${\displaystyle p(\ell \mid k)}$, that is, the degree distribution of the neighbors of a node with degree ${\displaystyle k}$, is given by[6]

${\displaystyle p(\ell \mid k)={\frac {m(k+2)}{k\ell (\ell +1)}}\left[1-{\frac {{\binom {2m+2}{m+1}}{\binom {k+\ell -2m}{\ell -m}}}{\binom {k+\ell +2}{\ell +1}}}\right].}$

In other words, if we select a node with degree ${\displaystyle k}$, and then select one of its neighbors randomly, the probability that this randomly selected neighbor will have degree ${\displaystyle \ell }$ is given by the expression ${\displaystyle p(\ell |k)}$ above.

### Clustering coefficient

While there is no analytical result for the clustering coefficient of the BA model, the empirically determined clustering coefficients are generally significantly higher for the BA model than for random networks. The clustering coefficient also scales with network size following approximately a power law

${\displaystyle C\sim N^{-0.75}.\,}$

Edit: Analytical result for the clustering coefficient of the BA model was obtained by Klemm and Eguíluz[7] and proven by Bollobás.[8][9] A mean-field approach to study the clustering coefficient was applied by Fronczak, Fronczak and Holyst.[10]

This behavior is still distinct from the behavior of small-world networks where clustering is independent of system size. In the case of hierarchical networks, clustering as a function of node degree also follows a power-law,

${\displaystyle C(k)=k^{-1}.\,}$

This result was obtained analytically by Dorogovtsev, Goltsev and Mendes.[11]

### Spectral properties

The spectral density of BA model has a different shape from the semicircular spectral density of random graph. It has a triangle-like shape with the top lying well above the semicircle and edges decaying as a power law. [12]

### Dynamic scaling

Generalized degree distribution ${\displaystyle F(q,t)}$ of the BA model for ${\displaystyle m=1}$.
The same data is plotted in the self-similar coordinates ${\displaystyle t^{1/2}F(q,N)}$ and ${\displaystyle q/t^{1/2}}$ and it gives an excellent collapsed revealing that ${\displaystyle F(q,t)}$ exhibit dynamic scaling.

By definition, the BA model describes a time developing phenomenon and hence, besides its scale-free property, one could also look for its dynamic scaling property. In the BA network nodes can also be characterized by generalized degree ${\displaystyle q}$, the product of the square root of the birth time of each node and their corresponding degree ${\displaystyle k}$, instead of the degree ${\displaystyle k}$ alone since the time of birth matters in the BA network. We find that the generalized degree distribution ${\displaystyle F(q,t)}$ has some non-trivial features and exhibits dynamic scaling

${\displaystyle F(q,t)\sim t^{-1/2}\phi (q/t^{1/2}).}$

It implies that the distinct plots of ${\displaystyle F(q,t)}$ vs ${\displaystyle q}$ would collapse into a universal curve if we plots ${\displaystyle F(q,t)t^{1/2}}$ vs ${\displaystyle q/t^{1/2}}$.[13]

## Limiting cases

### Model A

Model A retains growth but does not include preferential attachment. The probability of a new node connecting to any pre-existing node is equal. The resulting degree distribution in this limit is geometric,[14] indicating that growth alone is not sufficient to produce a scale-free structure.

### Model B

Model B retains preferential attachment but eliminates growth. The model begins with a fixed number of disconnected nodes and adds links, preferentially choosing high degree nodes as link destinations. Though the degree distribution early in the simulation looks scale-free, the distribution is not stable, and it eventually becomes nearly Gaussian as the network nears saturation. So preferential attachment alone is not sufficient to produce a scale-free structure.

The failure of models A and B to lead to a scale-free distribution indicates that growth and preferential attachment are needed simultaneously to reproduce the stationary power-law distribution observed in real networks.[2]

## History

The first use of a preferential attachment mechanism to explain power-law distributions appears to have been by Yule in 1925.[15] The modern master equation method, which yields a more transparent derivation, was applied to the problem by Herbert A. Simon in 1955[16] in the course of studies of the sizes of cities and other phenomena. It was first applied to the growth of networks by Derek de Solla Price in 1976[17] who was interested in the networks of citation between scientific papers. The name "preferential attachment" and the present popularity of scale-free network models is due to the work of Albert-László Barabási and Réka Albert, who rediscovered the process independently in 1999 and applied it to degree distributions on the web.[3]

## References

1. ^ Albert, Réka; Barabási, Albert-László (2002). "Statistical mechanics of complex networks". Reviews of Modern Physics. 74 (1): 47–97. arXiv:cond-mat/0106096. Bibcode:2002RvMP...74...47A. CiteSeerX 10.1.1.242.4753. doi:10.1103/RevModPhys.74.47. ISSN 0034-6861.
2. ^ a b c Albert, Réka; Barabási, Albert-László (2002). "Statistical mechanics of complex networks" (PDF). Reviews of Modern Physics. 74 (47): 47–97. arXiv:cond-mat/0106096. Bibcode:2002RvMP...74...47A. CiteSeerX 10.1.1.242.4753. doi:10.1103/RevModPhys.74.47. Archived from the original (PDF) on 2015-08-24.
3. ^ a b Barabási, Albert-László; Albert, Réka (October 1999). "Emergence of scaling in random networks" (PDF). Science. 286 (5439): 509&ndash, 512. arXiv:cond-mat/9910332. Bibcode:1999Sci...286..509B. doi:10.1126/science.286.5439.509. PMID 10521342. Archived from the original (PDF) on 2012-04-17.
4. ^ Korn, A.; Schubert, A.; Telcs, A. (2009). "Lobby index in networks". Physica A. 388 (11): 2221–2226. arXiv:0809.0514. Bibcode:2009PhyA..388.2221K. doi:10.1016/j.physa.2009.02.013.
5. ^ Cohen, Reuven; Havlin, Shlomo (2003). "Scale-Free Networks Are Ultrasmall". Physical Review Letters. 90 (5): 058701. arXiv:cond-mat/0205476. Bibcode:2003PhRvL..90e8701C. doi:10.1103/PhysRevLett.90.058701. ISSN 0031-9007. PMID 12633404.
6. ^ a b Fotouhi, Babak; Rabbat, Michael (2013). "Degree correlation in scale-free graphs". The European Physical Journal B. 86 (12): 510. arXiv:1308.5169. Bibcode:2013EPJB...86..510F. doi:10.1140/epjb/e2013-40920-6.
7. ^ Klemm, K.; Eguíluz, V. C. (2002). "Growing scale-free networks with small-world behavior". Physical Review E. 65 (5): 057102. arXiv:cond-mat/0107607. Bibcode:2002PhRvE..65e7102K. doi:10.1103/PhysRevE.65.057102. hdl:10261/15314. PMID 12059755.
8. ^ Bollobás, B. (2003). "Mathematical results on scale-free random graphs". Handbook of Graphs and Networks. pp. 1–37. CiteSeerX 10.1.1.176.6988.
9. ^ http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.176.6988
10. ^ Albert, Reka; Barabasi, Albert-Laszlo; Hołyst, Janusz A (2003). "Mean-field theory for clustering coefficients in Barabasi-Albert networks". Phys. Rev. E. 68 (4): 046126. arXiv:cond-mat/0306255. doi:10.1103/PhysRevE.68.046126. PMID 14683021.
11. ^ Dorogovtsev, S.N.; Goltsev, A.V.; Mendes, J.F.F. (25 June 2002). "Pseudofractal scale-free web". Physical Review E. 65 (6): 066122. arXiv:cond-mat/0112143. Bibcode:2002PhRvE..65f6122D. doi:10.1103/PhysRevE.65.066122. PMID 12188798.
12. ^ Farkas, I.J.; Derényi, I.; Barabási, A.-L.; Vicsek, T. (20 July 2001) [19 February 2001]. "Spectra of "real-world" graphs: Beyond the semicircle law". Physical Review E. 64 (2): 1–12. arXiv:cond-mat/0102335. Bibcode:2001PhRvE..64b6704F. doi:10.1103/PhysRevE.64.026704. hdl:2047/d20000692. PMID 11497741.
13. ^ M. K. Hassan, M. Z. Hassan and N. I. Pavel, “Dynamic scaling, data-collapseand Self-similarity in Barabasi-Albert networks” J. Phys. A: Math. Theor. 44 175101 (2011) https://dx.doi.org/10.1088/1751-8113/44/17/175101
14. ^ Pekoz, Erol; Rollin, A.; Ross, N. (2012). "Total variation and local limit error bounds for geometric approximation". Bernoulli.
15. ^ Yule, G. Udny; Yule, G. Udny (1925). "A Mathematical Theory of Evolution Based on the Conclusions of Dr. J. C. Willis, F.R.S". Journal of the Royal Statistical Society. 88 (3): 433–436. doi:10.2307/2341419. JSTOR 2341419.
16. ^ Simon, Herbert A. (December 1955). "On a Class of Skew Distribution Functions". Biometrika. 42 (3–4): 425&ndash, 440. doi:10.1093/biomet/42.3-4.425.
17. ^ Price, D.J. de Solla (September 1976). "A general theory of bibliometric and other cumulative advantage processes". Journal of the American Society for Information Science. 27 (5): 292–306. CiteSeerX 10.1.1.161.114. doi:10.1002/asi.4630270505.