Mean Hitting Time for Random Walks on a Class of Sparse Networks.

For random walks on a complex network, the configuration of a network that provides optimal or suboptimal navigation efficiency is meaningful research. It has been proven that a complete graph has the exact minimal mean hitting time, which grows linearly with the network order. In this paper, we present a class of sparse networks G(t) in view of a graphic operation, which have a similar dynamic process with the complete graph; however, their topological properties are different. We capture that G(t) has a remarkable scale-free nature that exists in most real networks and give the recursive relations of several related matrices for the studied network. According to the connections between random walks and electrical networks, three types of graph invariants are calculated, including regular Kirchhoff index, M-Kirchhoff index and A-Kirchhoff index. We derive the closed-form solutions for the mean hitting time of G(t), and our results show that the dominant scaling of which exhibits the same behavior as that of a complete graph. The result could be considered when designing networks with high navigation efficiency.

1. Introduction

A complex network is recognized as a powerful tool for revealing the mysteries of complex systems [1]. It is widely used in metabolic networks [2], software engineering [3], ecosystems [4] and so on. In addition to some topological parameters, such as power-law degree distribution, average path length and clustering coefficient of complex network, the random walks also received widespread attention because the research of random walk theory can disclose dynamic processes on complex networks. As a key quantity of random walks, the hitting time is related to the mixing rate of an irreducible Markov chain, and it is also considered when calculating the expected time of mixing the Markov chain [5]. The hitting time can be used to measure the navigation efficiency of the network [6,7], and it has a core position in different disciplines, including mathematics, computer, biology, physics, control science and engineering [8,9,10,11,12].

Most of the previous research on random walks about complex networks focuses on two aspects: one is that the nodes of the studied network have identical walking rules [13,14], and the other is to study random walks on heterogeneous networks that set a trap on the node with the largest degree and have scale-free characteristics [15,16,17]. Since many real networks have a scale-free nature, every node in the network can be a trap. Thus, we construct a deterministic network that satisfies the above restrictions to approximate the real network, and it is more beneficial for us to evaluate the dynamic process of the network.

It has been proven that the complete graph has the minimum mean hitting time among all undirected networks, which shows that its propagation is quite efficient [18]. For the purpose of constructing a highly efficient network and controlling its trapping process, it is necessary to explore and design some networks with a small mean hitting time. Since most real networks are sparse networks, the average degree of these networks is much less than that of a complete graph. In this paper, we design and analyze a class of sparse networks with scale-free properties; their topological properties are different from the complete graph. We prove that the dominant scaling of the mean hitting time exhibits the same behavior as that of a complete graph, and they can also have high navigation efficiency.

The main contents in the other sections of this paper are as follows. In Section 2, we propose a graphic operation and design a class of sparse networks, show their differences in several topological parameters, including average degree, degree distribution, clustering coefficient and diameter. In Section 3, we present some lemmas about electrical networks and random walks. In Section 4, we analytically obtain the closed-form solution of the mean hitting time according to the connections between the Kirchhoff index and the mean hitting time. In Section 5, we conclude our work with a concise narrative.

2. Topological Characteristics of the Network

Before proceeding, we propose a graphic operation called the rhombus operation and construct a network by iterating the rhombus operation. Then, we compare the topological features of and the complete graph with the same network order.

Rhombus operation: For a given edge with two endnodes i and j, add two new nodes to both sides of this edge, denoted by u and v, and then connect edge , respectively. Figure 1a shows the operation process of a rhombus operation.

Figure 1

(a) The illustration of rhombus operation; (b) the network at time step ; (c) the network at time step .

With the preparation of a graphic operation, we show the construction rule of networks as follows. Initial state, , is only an edge. For , can be born from by performing a rhombus operation on every edge in . Figure 1b,c illustrate the topological structure of and .

The iterative construction allows us to precisely analyze relevant topological properties of the network. Let and be the node set and edge set of ; in more detail, the new node set and new edge set at time step t are denoted as and , which means , and the node belonging to is called the old node. Then, the number of nodes and edges of the network are denoted as and , respectively. The recursive relation is obviously established according to the rhombus operation, and we can obtain due to . Additionally, we have , so it is easy to verify that . Let be the degree of node i in that was generated at iteration , which satisfies that the degree of node i at time step t is 3 times the degree of the previous time step ; that is, .

For the network

2.1. Average Degree

When , the condition is clearly established, so our network model is a sparse network according to literature [19]. However, the complete graph is not sparse. For a complete graph with the same number of nodes, the average degree of is due to the degree of each node is .

2.2. Cumulative Degree Distribution

In real life, there are few fully connected networks like a complete graph. Most real networks exhibit scale-free nature, and their nodes with a large degree are fewer, but nodes with a small degree are the majority of the network. A network is said to be scale-free when its cumulative degree distribution obeys , where , and the cumulative degree distribution represents the probability that the degree of a node is equal to or greater than k, where is a probability of a randomly selected node with k neighbors in the network .

The cumulative degree distribution of the sparse network

The degree of node i will increase by a factor 3, that is, , which shows that the degree spectrum of is discrete. In Table 1, we enumerate the degree k and the number of nodes with degree k, then the cumulative degree distribution of is calculated by where has been substituted into the above formula, and for large t, the cumulative degree distribution follows a power law with exponent . Therefore, we have proven that the network is a scale-free network. On the other hand, the degree of all nodes in the complete graph is the same, so is not a scale-free network. It can be seen that our network is more suitable for simulating scale-free real networks. □

Table 1

The degree and clustering coefficient of node in .

t	k	n(k)	c(k)
0	3^t	2	2/3^t
1	2×3^t−1	2×5^0	1/3^t−1
2	2×3^t−2	2×5^1	1/3^t−2
⋯	⋯	⋯	⋯
t_i	2×3^t−t_i	2×5^t_i−1	1/3^t−t_i
⋯	⋯	⋯	⋯
t−1	2×3^1	2×5^t−2	1/3^1
t	2×3^0	2×5^t−1	1/3^0

2.3. Clustering Coefficient

The clustering coefficient is used to describe the tightness of clumps between nodes in the graph. Specifically, it is a measure of the dense or sparse connection between neighbors of each node. The clustering coefficient of each node i is the ratio between the number of edges that actually exist in all nearest neighbors and the number of all possible edges between them, and it is expressed as . The whole clustering coefficient of the network is the average value of over all nodes in the network, which can be written as .

For the network

Since the clustering coefficient of each node with the same degree in is also the same, let be the number of nodes with degree k, and represents the clustering coefficient of each node with degree k. For each degree k, we can calculate the clustering coefficient and the corresponding number of nodes with degree k, as shown in Table 1, so the whole clustering coefficient of network can be calculated as follows: we have when , and is not a highly clustered network. However, the clustering coefficient of the complete graph is equal to 1, and there is a significant difference between and on this topological parameter. □

2.4. Diameter

The diameter is defined as the maximum of the shortest distances between all pairs of nodes in network , denoted as , and it is often used to characterize the longest communication delay in complex network.

For

When is a small network, we can easily enumerate its diameter, such as , . At time step , the distance between two new nodes is the longest, and the path through them must contain an old node, so . For , we can find that the diameter refers to the distance between a pair of new nodes. For simplicity of description, we denote the newly generated node at time as ; therefore, we only need to consider the maximum value of the shortest distance between two nodes t. According to the structure of , when , the shortest path between new nodes must be the path , which contains nodes generated at time step . For , the shortest path between two new nodes must be a path . Hence, is true for , which shows that the diameter scales logarithmically with the network order. For a complete graph with the same number of nodes, it is well known that its diameter is equal to 1, which means that the diameter of is larger than that of . □

3. Random Walks and Electrical Networks

In this section, we aim to show the closed-form solution for the mean hitting time of our networks. Firstly, we introduce several notions and lemmas about electrical networks and random walks, then we provide the relationships between the mean hitting time and Kirchhoff index. The electrical network corresponding to a graph can be constructed by replacing each edge in with a unit resistor, but we still denote the resulting electrical network as . The effective resistance between any two distinct nodes is defined as the potential difference between them when a unit current from i to j is kept; when , we set .

([20]). For an electrical network with  nodes, the sum of effective resistances between all pairs of adjacent nodes can be written as

([21]). For any pair of distinct node i and j in an electrical network

According to [22], the regular Kirchhoff index of a network is defined as the sum of the effective resistances of all pairs of disordered nodes in :

Taking into account the influence of degree on the Kirchhoff index, the M-Kirchhoff index and the A-Kirchhoff index have been proposed in [23,24], respectively, and they are interpreted by the formula as and

The unbiased discrete time random walks means that the particle starting from the current location jumps to each of its neighboring nodes with equal probability at every time step [25]. The hitting time of network is a key quantity pertaining to random walk, and it is defined as the expected time taken by a particle jumping to the ending node j from the starting node i for the first time. The mean hitting time is the average of hitting times over all node pairs [26], and it can be solved by the and the network order and size of .

For a network

The Kirchhoff index can be represented in terms of the non-zero eigenvalues of the Laplacian matrix for network as [24]. The mean hitting time of network is the average of hitting times over all node pairs, and it is expressed as , where is the hitting time from node i to another node j. In addition, can be expressed in terms of the non-zero eigenvalues of the Laplacian matrix [18], that is , then we can obtain the equation by combining the above equations. □

3.1. Related Matrices

All the nodes of a given network are marked as , respectively, and the adjacency relations between all nodes and edges are implicit in an adjacency matrix , where if node i and j are connected by an edge in the network, and if there is no edge between i and j. Let be the diagonal degree matrix of . Its i-th diagonal entry is the degree of node i, and the remaining entries are zero. The Laplacian matrix of is defined as .

Use and to abbreviate the old node set and the new node set in network , and the number of new nodes is . The network is generated iteratively by , then we show the recursive relationship between two consecutive time steps of these matrices. The adjacency matrix can be written in block form as where is obvious according to the definition of adjacency matrix, and is the zero matrix with order . On the other hand, the diagonal matrix satisfies where the symbol I represents the identity matrix with order . Equation (13) is based on the fact that the nodes contained in set are 2 degree nodes, the degree of every node in set increases by a factor 3. Thus, the Laplacian matrix of network can be expressed as

For the sparse network

The left side and right side of equation are denoted by and , respectively; thus, the entries of are

□

Our main task is to verify that the entries of are equivalent to those of . Matrix can be partitioned into column vectors for , that is , and we have due to , and we can calculate the product of the two matrices as .

The entries of can be determined by distinguishing two cases. (a) When , the diagonal entry is , and we can obtain . (b) When , the non-diagonal entry of matrix is equal to

Before going on, we introduce a concept about -inverse of a matrix [27]. Matrix M is called a -inverse of X if holds, let be one of the -inverses of X. A lemma about the -inverse of a block matrix is shown below.

([18]). Let matrix

3.2. Effective Resistances

For a connected network , the effective resistance between any pair of nodes can be obtained from the elements of any -inverse of its Laplacian matrix, and we can refer to the following lemma.

([28]). For a given

Next, we show that the effective resistance between any two nodes in can be represented in terms of effective resistances of node pairs in . In the following calculation process, we divide the nodes into the old node and the new node to investigate the effective resistance between them. To achieve this goal, we introduce some variables and define for any two subsets X and Y of set on . Then, for a node in , we define , where is the neighbors set of node i and .

For the effective resistance between the node pairs in the network

(1) Let

(2) Let

(3) Let

(4) Let

The detailed proof of Lemma 6 is given in Appendix A.

4. Mean Hitting Time

Based on the above preparations, we determine the mean hitting time for network using the connection between the mean hitting time and the Kirchhoff index. Firstly, we calculate the exact solutions of three auxiliary variables, including , , and for two subsets X and Y of set in network , and give the relationships between them. The following lemmas can support our main result.

For network

(a)

(b)

(a) Since each edge of network can generate two new nodes of network , then we have (b) For any old node , there are new nodes in that are adjacent to j, so is summed times. □

The M-Kirchhoff index and the A-Kirchhoff index of our network and

The proof of Lemma 8 is given in Appendix B.

The solution of regular Kirchhoff index of our network

Through the two types of Kirchhoff indexes obtained by Lemma 8, we deduce the relationship between the regular Kirchhoff index and them. Divide all the nodes in the network into new node and old node, then is equal to

□

Considering , plugging Equations (21) and (22) into Equation (24) yields the solution of Kirchhoff index of network , as shown in Equation (23). Figure 2 shows a schematic diagram of the three types of Kirchhoff indexes of network . We are ready to show the result for the mean hitting time of .

Figure 2

The schematic diagram of the three types of Kirchhoff indexes.

For then

Since Lemma 3 and the total number of edges of whole network is , we have

It is easy to obtain the result in Equation (25) by substituting Equation (23) into Equation (26). We continue to express as a function of the network order , it can be observed that from the exact value of ; hence, the mean hitting time of can be expressed in terms of network order as Therefore, when , for a large network, we have which shows increases linearly with the total number of nodes in our network; the mean hitting time of the random walks shown is similar to that of the complete graph; and they all have high transmission efficiency. □

5. Conclusions

In this paper, we have presented a class of a sparse network and have pointed out the differences between and the complete graph with the same order in several topological characteristics. The main differences are has a scale-free property, while does not. The scale-free feature is a shock discovery in real complex systems. is not sparse, but is sparse; it is rare to achieve a tight connections like a complete graph in a real network. It has been proven that the mean hitting time of the complete graph is minimal and increases linearly with the network order. Based on the relationship between the mean hitting time and the Kirchhoff index, we have calculated a closed-form solution to the mean hitting time of our network, and the result shows that the dominant scaling of which exhibits the same behavior as that of a complete graph. We hope that our work will be instructive for the design and construction of complex networks with efficient navigation.