Complex network zeta function

Different definitions have been given for the dimension of a complex network or graph. For example, metric dimension is defined in terms of the resolving set for a graph. Dimension has also been defined based on the box covering method applied to graphs.^[1] Here we describe the definition based on the complex network zeta function.^[2] This generalises the definition based on the scaling property of the volume with distance.^[3] The best definition depends on the application.

Definition

One usually thinks of dimension for a set which is dense, like the points on a line, for example. Dimension makes sense in a discrete setting, like for graphs, only in the large system limit, as the size tends to infinity. For example, in Statistical Mechanics, one considers discrete points which are located on regular lattices of different dimensions. Such studies have been extended to arbitrary networks, and it is interesting to consider how the definition of dimension can be extended to cover these cases. A very simple and obvious way to extend the definition of dimension to arbitrary large networks is to consider how the volume (number of nodes within a given distance from a specified node) scales as the distance (shortest path connecting two nodes in the graph) is increased. For many systems arising in physics, this is indeed a useful approach. This definition of dimension could be put on a strong mathematical foundation, similar to the definition of Hausdorff dimension for continuous systems. The mathematically robust definition uses the concept of a zeta function for a graph. The complex network zeta function and the graph surface function were introduced to characterize large graphs. They have also been applied to study patterns in Language Analysis. In this section we will briefly review the definition of the functions and discuss further some of their properties which follow from the definition.

We denote by $\textstyle r_{ij}$ the distance from node $\textstyle i$ to node $\textstyle j$, i.e., the length of the shortest path connecting the first node to the second node. $\textstyle r_{ij}$ is $\textstyle \infty$ if there is no path from node $\textstyle i$ to node $\textstyle j$. With this definition, the nodes of the complex network become points in a metric space.^[2] Simple generalisations of this definition can be studied, e.g., we could consider weighted edges. The graph surface function, $\textstyle S(r )$, is defined as the number of nodes which are exactly at a distance $\textstyle r$ from a given node, averaged over all nodes of the network. The complex network zeta function $\textstyle \zeta_G ( \alpha )$ is defined as

$\zeta_G ( \alpha ) := \frac{1}{N}\sum_i \sum_{j\neq i }r^{-\alpha}_{ij},$

where $\textstyle N$ is the graph size, measured by the number of nodes. When $\textstyle \alpha$ is zero all nodes contribute equally to the sum in the previous equation. This means that $\textstyle \zeta_{G}(0)$ is $\textstyle N-1$, and it diverges when $\textstyle N \rightarrow \infty$. When the exponent $\textstyle \alpha$ tends to infinity, the sum gets contributions only from the nearest neighbours of a node. The other terms tend to zero. Thus, $\textstyle \zeta_G ( \alpha )$ tends to the average degree $\textstyle <k>$ for the graph as $\textstyle \alpha \rightarrow \infty$.

$\langle k \rangle = \lim_{\alpha \rightarrow \infty} \zeta_G ( \alpha ).$

The need for taking an average over all nodes can be avoided by using the concept of supremum over nodes, which makes the concept much easier to apply for formally infinite graphs.^[4] The definition can be expressed as a weighted sum over the node distances. This gives the Dirichlet series relation

$\zeta_G ( \alpha ) = \sum_r S(r)/r^{\alpha}.$

This definition has been used in the shortcut model to study several processes and their dependence on dimension.

Properties

$\textstyle \zeta_G ( \alpha )$ is a decreasing function of $\textstyle \alpha$, $\textstyle \zeta_G ( \alpha_1 ) > \zeta_G ( \alpha_2 )$, if $\textstyle \alpha_1 < \alpha_2$. If the average degree of the nodes (the mean coordination number for the graph) is finite, then there is exactly one value of $\textstyle \alpha$, $\textstyle \alpha_{transition}$, at which the complex network zeta function transitions from being infinite to being finite. This has been defined as the dimension of the complex network. If we add more edges to an existing graph, the distances between nodes will decrease. This results in an increase in the value of the complex network zeta function, since $\textstyle S(r)$ will get pulled inward. If the new links connect remote parts of the system, i.e., if the distances change by amounts which do not remain finite as the graph size $\textstyle N \rightarrow \infty$, then the dimension tends to increase. For regular discrete d-dimensional lattices $\textstyle \mathbf Z^d$ with distance defined using the $\textstyle L^1$ norm

$\|\vec{n}\|_1=\|n_1\|+\cdots +\|n_d\|,$

the transition occurs at $\textstyle \alpha = d$. The definition of dimension using the complex network zeta function satisfies properties like monotonicity (a subset has a lower or the same dimension as its containing set), stability (a union of sets has the maximum dimension of the component sets forming the union) and Lipschitz invariance,^[5] provided the operations involved change the distances between nodes only by finite amounts as the graph size $\textstyle N$ goes to $\textstyle \infty$. Algorithms to calculate the complex network zeta function have been presented.^[6]

Values for discrete regular lattices

For a one-dimensional regular lattice the graph surface function $\textstyle S_{1}(r)$ is exactly two for all values of $\textstyle r$ (there are two nearest neighbours, two next-nearest neighbours, and so on). Thus, the complex network zeta function $\textstyle \zeta_G ( \alpha )$ is equal to $\textstyle 2\zeta(\alpha)$, where $\textstyle \zeta(\alpha)$ is the usual Riemann zeta function. By choosing a given axis of the lattice and summing over cross-sections for the allowed range of distances along the chosen axis the recursion relation below can be derived

$S_{d+1}(r) = 2+S_d (r)+2\sum^{r-1}_{i=1}S_d(i).$

From combinatorics the surface function for a regular lattice can be written^[7] as

$S_d(r) = \sum^{d-1}_{i=0}(-1)^{i}2^{d-i}{d \choose i} { d+r-i-1 \choose d-i-1 }.$

The following expression for the sum of positive integers raised to a given power $\textstyle k$ will be useful to calculate the surface function for higher values of $\textstyle d$:

$\sum^{r}_{i=1}i^k = \frac{r^{k+1}}{(k+1)}+ \frac{r^k}{2}+\sum^{\lfloor (k+1)/2 \rfloor}_{j=1}\frac{(-1)^{j+1}2\zeta(2j)k!r^{k+1-2j}}{(2\pi)^{2j}(k+1-2j)!}.$

Another formula for the sum of positive integers raised to a given power $\textstyle k$ is

$\sum^{n}_{k=1}\bigl(\begin{smallmatrix} n+1\ k \end{smallmatrix}\bigr)\sum^{r}_{i=1}i^{k} = (r+1)((r+1)^{n}-1).$

$\textstyle S_{d}(r) \rightarrow O(2^{d}r^{d-1}/\Gamma(d))$ as $\textstyle r\rightarrow \infty$.

The Complex network zeta function for some lattices is given below.

$\textstyle d=1$  : $\textstyle \zeta_G(\alpha)=2\zeta(\alpha)$

$\textstyle d=2$  : $\textstyle \zeta_G(\alpha)=4\zeta(\alpha-1)$

$\textstyle d=3$  : $\textstyle \zeta_G(\alpha)=4\zeta(\alpha-2)+2\zeta(\alpha$)

$\textstyle d=4$  : $\textstyle \zeta_G(\alpha)=\frac{8}{3}\zeta(\alpha-3)+\frac{16}{3}\zeta(\alpha-1)$

$\textstyle r\rightarrow \infty$ : $\textstyle \zeta_{G}(\alpha)=2^{d}\zeta(\alpha-d+1)/\Gamma(d)$ (for $\alpha$ near the transition point.)

Random graph zeta function

Random graphs are networks having some number $\textstyle N$ of vertices, in which each pair is connected with probability $\textstyle p$, or else the pair is disconnected. Random graphs have a diameter of two with probability approaching one, in the infinite limit ($\textstyle N \rightarrow \infty$). To see this, consider two nodes $\textstyle A$ and $\textstyle B$. For any node $\textstyle C$ different from $\textstyle A$ or $\textstyle B$, the probability that $\textstyle C$ is not simultaneously connected to both $\textstyle A$ and $\textstyle B$ is $\textstyle (1-p^2)$. Thus, the probability that none of the $\textstyle N-2$ nodes provides a path of length $\textstyle 2$ between nodes $\textstyle A$ and $\textstyle B$ is $\textstyle (1-p^2)^{N-2}$. This goes to zero as the system size goes to infinity, and hence most random graphs have their nodes connected by paths of length at most $\textstyle 2$. Also, the mean vertex degree will be $\textstyle p(N-1)$. For large random graphs almost all nodes are at a distance of one or two from any given node, $\textstyle S(1)$ is $\textstyle p(N-1)$, $\textstyle S(2)$ is $\textstyle (N-1)(1-p)$, and the graph zeta function is

$\zeta_G ( \alpha ) = p(N-1) + (N-1)(1-p)2^{-\alpha}.$

References

1. K.-I. Goh, G. Salvi, B. Kahng and D. Kim, Phys. Rev. Lett. 96, 018701 (2006).
2. O. Shanker (2007). "Graph Zeta Function and Dimension of Complex Network". Modern Physics Letters B 21 (11): 639–644. Bibcode:2007MPLB...21..639S. doi:10.1142/S0217984907013146.
3. O. Shanker (2007). "Defining Dimension of a Complex Network". Modern Physics Letters B 21 (6): 321–326. Bibcode:2007MPLB...21..321S. doi:10.1142/S0217984907012773.
4. O. Shanker (2010). "Complex Network Dimension and Path Counts". Theoretical Computer Science 411 (26–28): 2454–2458. doi:10.1016/j.tcs.2010.02.013.
5. K. Falconer, Fractal Geometry: Mathematical Foundations and Applications, Wiley, second edition, 2003
6. O. Shanker, (2008). "Algorithms for Fractal Dimension Calculation". Modern Physics Letters B 22 (7): 459–466. Bibcode:2008MPLB...22..459S. doi:10.1142/S0217984908015048.
7. O. Shanker (2008). "Sharp dimension transition in a shortcut model". J. Phys. A: Math. Theor. 41 (28): 285001. Bibcode:2008JPhA...41B5001S. doi:10.1088/1751-8113/41/28/285001.