# Random walks in modular scale-free networks with multiple traps

###### Abstract

Extensive empirical investigation has shown that a plethora of real networks synchronously exhibit scale-free and modular structure, and it is thus of great importance to uncover the effects of these two striking properties on various dynamical processes occurring on such networks. In this paper, we examine two cases of random walks performed on a class of modular scale-free networks with multiple traps located at several given nodes. We first derive a formula of the mean first-passage time (MFPT) for a general network, which is the mean of the expected time to absorption originating from a specific node, averaged over all non-trap starting nodes. Although the computation is complex, the expression of the formula is exact; moreover, the computational approach and procedure are independent of the number and position of the traps. We then determine analytically the MFPT for the two random walks being considered. The obtained analytical results are in complete agreement with the numerical ones. Our results show that the number and location of traps play an important role in the behavior of the MFPT, since for both cases the MFPT grows as a power-law function of the number of nodes, but their exponents are quite different. We demonstrate that the root of the difference in the behavior of MFPT is attributed to the modular and scale-free topologies of the networks. This work can deepen the understanding of diffusion on networks with modular and scale-free architecture and motivate relevant studies for random walks running on complex random networks with multiple traps.

###### pacs:

05.40.Fb, 89.75.Hc, 05.60.Cd, 89.75.Da## I Introduction

In the past decade, with a huge amount of data and computational resources available, scientists have processed and analyzed data of a wide variety of real systems in different areas, leading to important advances in the understanding of complex systems AlBa02 (); DoMe02 (); Ne03 (); BoLaMoChHw06 (). A large volume of empirical studies showed that scale-free feature BaAl99 () and modular structure GiNe02 (); RaSoMoOlBa02 (); RaBa03 () are two prominent properties that seem to be common to real networks, especially biological and social networks. The former implies that the networks obey a power-law degree distribution as with , while the latter means that the networks can be divided into groups (modules), within which nodes are more tightly connected with each other than with nodes outside. These two remarkable natures constitute our fundamental understanding of the structure of complex networks, which are relevant to other topological features (i.e., average distance ChLu02 (); CoHa03 () and clustering coefficient RaBa03 ()), and have led to many popular topics of research in network science, including explaining the origin of the scale-free phenomenon AlBa02 (); DoMe02 (), identifying the modules DaDuDiAr05 (); PaDeFaVi05 (); Ne06 (); Fo10 (); ZhZhZhGu11 () and finding their concrete applications Pe01 (); WuGaHa04 ().

It is well known that one of the ultimate goals for research on complex networks is to make clear how the underlying structural characteristics affect the dynamical processes defined on networks Ne03 (); DoGoMe08 (). Among various dynamical processes, random walks have held continual interest within the scientific community HaBe87 (); MeKl00 (); MeKl04 (); BuCa05 (); SoMaBl97 (); PaAm04 (); NoRi04 (); BeCoMoSuVo05 (); SoRebe05 (); Bobe05 (); GaSoHaMa07 (); BaBeWi08 (); BaCaPa08 (); KiCaHaAr08 (); ZhZhZhYiGu09 (); HaRo09 (); ZhAlHoZhCh11 () because of their relevance to a wide range of different applications to many fields We94 (); Hu95 (). In particular, as an integral subject of random walks, the trapping problem is closely related to numerous aspects in a great many disciplines JaBl01 (); Sh05 (); WhGo99 (); FoPiReSa07 (). Over the past decades, scholars in a large interdisciplinary community have made a huge effort to address the trapping problem in diverse networks, including regular lattices Mo69 (), regular fractals KaBa02PRE (); KaBa02IJBC (); ZhWuZhZhGuWa10 (); Ag08 (); ZhLiZhWuGu09 (), small-world networks CaAb08 (), and scale-free networks ZhQiZhXiGu09 (); ZhZhXiChLiGu09 (); ZhGuXiQiZh09 (); AgBu09 (); TeBeVo09 (); ZhXiZhGaGu09 (); ZhXiZhLiGu09 (), among other graphs CoBeTeVoKl07 (); BeMeTeVo08 (); CoTeVoBeKl08 ().

Thus far, most previous works on random walks in complex networks have focused on the case with a single trap fixed at a given location, while work on the case with multiple traps is much less common. In particular, research on the multiple-trap problem in complex networks with modular organization and scale-free structure is still lacking, despite the multiple-trap issue having obvious applications to various aspects YuAc03 () (description of particle-cluster aggregation WiSa82 (); WiSa83 (), for instance) and being relevant in diffusion-limited reactions in chemical field Ri85 () and modular and scale-free topologies having vital influence on dynamical processes taking place on networks Hi07 (); Ma07 (); Di08 (); ZhLiGoZhGuLi09 (); BuSp09 (). ¡¡

In this paper, we study the classic random-walk problem for a category of modular scale-free networks RaSoMoOlBa02 (); RaBa03 () with several given nodes being occupied by immobile traps, which absorb all particles visiting it. The basic quantity we are interested in is the mean first-passage time (MFPT) Re01 () characterizing the trapping process, which is defined as the average of expected time for a particle starting off from a particular node until first visiting one of the traps, averaged over all nontrap source nodes. The networks we study are of a deterministic family, which has proven to be an important tool in the field of complex networks and has recently attracted much interest DoGoMe02 (); JuKiKa02 (); CoOzPe00 (); ZhRoGo06 (); HiBe06 (); RoHaAv07 (); ZhZhChGu07 (); ZhZhFaGuZh07 (); BoGoGu08 (); CoMi09 (); MoRoDa09 (); ZhGuDiChZh09 (); MiCoChZh10 ().

We first study the MFPT in a generic network and reduce the problem of computing the MFPT to finding the sum of elements of a matrix associated with the trapping issue. Then, based on the deterministic recursive construction of the modular scale-free networks being considered, we investigate analytically the key of quantity MFPT for two cases of particular arrangements of traps. In the first case, traps are placed on peripheral nodes; in the other case, traps are fixed on those nodes farthest from the hub. For both cases, we derive exactly the dominant scalings for the MFPT and show that they produce a power-law function of the network size with their exponents smaller than 1 but being different, which is confirmed by the numerical results obtained via inverting related matrices. The obtained results indicate that both the number and the location of traps have a significant impact on the behavior of the trapping. We demonstrate that the high efficiency of both trapping processes and the distinction between the behavior of the MFPT for the two random walks are rested with the scale-free property and the modular structure of the networks under consideration.

## Ii Modular scale-free networks

We first introduce the model for the modular scale-free networks, which are built in an iterative way RaSoMoOlBa02 (); RaBa03 (). Let stand for the network model after () iterations (i.e., number of generations). Initially (), the model is composed of () nodes linked by edges forming a complete graph, among which a node (e.g., the central node in Fig. 1) is called hub (or root) node, and the other nodes are named peripheral nodes. At the second generation (), replicas of are created with the peripheral nodes of each copy being connected to the root of the original . In this way, we obtain , the hub and peripheral nodes of which are the hub of the original and the peripheral nodes in the duplicates of , respectively. Supposing one has , the next generation network can be obtained by adding copies of to the primal , with all peripheral nodes of the replicas being linked to the hub of the original unit. The hub of the original and the peripheral nodes of the copies of form the hub node and peripheral nodes of , respectively. Repeating indefinitely the two steps of replication and connection, one obtains the modular scale-free networks. Figure 1 illustrates a network for the particular case of .

Many interesting quantities and properties of the model can be determined explicitly RaBa03 (); No03 (). In , the network size (number of nodes), denoted by , is . All these nodes can be classified into four distinct sets No03 (); NoRi04a (): the peripheral node set , the locally peripheral node set (), the set containing only the hub node of , and the local hub set (). The number of nodes in each of these four sets is

(1) |

(2) |

(3) |

and

(4) |

respectively. For , all nodes in a set have the same degree. It has been obtained exactly that the degree for a node in sets , , , and is, respectively,

(5) |

(6) |

(7) |

and

(8) |

In addition, it is easy to obtain that the average degree of all nodes is approximately equal to a constant in the limit of infinite , showing that the networks are sparse.

The model under consideration is in fact an extension of the one proposed in Ref. BaRaVi01 () and studied in much detail in Refs. IgYa05 (); ZhLiGaZhGu09 (); AgBuMa10 (). It presents some typical features observed in a variety of real-world systems RaBa03 (); No03 (). Its degree distribution follows a power-law scaling with a general exponent belonging to the interval . Its average clustering coefficient tends to a large constant dependent on ; and its average distance grows logarithmically with the network order ZhLi10 (), both of which show that the model is small world WaSt98 (). In addition, the betweenness distribution of nodes also obeys a power-law behavior with the exponent independent of the parameter . In particular, the whole class of networks shows a remarkable modular structure. These peculiar structural properties make the networks unique within the category of complex networks. It is thus interesting to address dynamical processes happening on them. The main purpose of this work is to study random walks on this network family with multiple traps located on some special nodes.

## Iii Formulation of random walks with multiple traps on a network

In this section, we formulate the problem of random walks on the network with multiple traps, which is a discrete-time random walk of a particle in the presence of several perfect traps placed on certain nodes. At each time step, the particle jumps with equal probability from its current location to one of its nearest neighbors. If the particle meets one of the traps, then it is absorbed. At last, the particle will be inevitably absorbed by the traps, regardless of its starting position KeSn76 (); AlFi99 ().

It is well known that an arbitrary network can be completely represented by its adjacency matrix. For , its adjacency matrix A is a matrix consisting of entries 0 or 1, with an order ( is the number of nodes in ). The element of A is defined as follows: if and are neighbors and otherwise. Then the degree of node is given by , the diagonal degree matrix Z associated with is , and the corresponding normalized Laplacian matrix of is defined to be , where I is the identity matrix with order .

We use to denote the set of traps and to represent the number of traps. We are concerned with the expected time the particle spends, starting from a source node, before it falls on one of the traps for the first time. Let be the expected time, frequently called first-passage time (FPT) or trapping time, for a particle first arriving at any one of the traps, given that it starts from node . It is clear that for any node , we have . The set of this important quantity satisfies the relation

(9) |

where and is transition probability for the particle of going from node to node . According to the definition of the random-walk problem, it is not difficult to know that , which is exactly the element of the matrix .

In order to facilitate the description, we distinguish all nodes in by assigning each of them a unique number. We label consecutively all nodes, excluding those in , from 1 to and trap nodes are numbered from to . Then Eq. (9) can be rewritten in matrix form as

(10) |

where (the superscript of the vector represents transpose) is an -dimensional vector, e is the -dimensional unit vector , and is the transition matrix corresponding to the trapping problem. Equation (10) can be further recast as

(11) |

where

(12) |

with being the identity matrix.

It should be mentioned that the considered random walk is in fact a Markov process, and Eq. (12) is the fundamental matrix of the Markov chain representing such an unbiased random walk. We also note that the matrix on the right-hand side of Eq. (12) is actually a submatrix of the normalized discrete Laplacian matrix L of , which is obtained from L by suppressing the last rows and columns that correspond to the trap nodes.

Equation (11) shows that trapping time can be expressed in terms of the entries of inverse matrix of (i.e., a submatrix of L). Concretely, is provided by

(13) |

which accounts for the Markov chain representing the random walk: The entry of the fundamental matrix for the Markov process represents the expected number of times the particle visits node in the case that it starts off from node (see Ref. BaKl98 () for a single trap). Then the MFPT , which is defined as the average of over all initial nodes distributed uniformly over nodes in including the traps note (), is given by

(14) | |||||

Thus, on the basis of the definition of the unbiased random walks, we have derived the numerical yet exact solution to the MFPT for random walks on any network with multiple traps, independently of the number and location of the traps. We note that our derivation is a reformulation of the backward equation satisfied by the MFPT and that Eq. (14) can also be found in the literature in several equivalent forms KeSn76 (); AlFi99 ().

Equation (14) is very important, since it reduces the problem of calculating the MFPT to computing the sum of the elements of the matrix and can be used to check the results for derived by other methods, at least for networks with a small number of nodes. However, it is notable that although the above computational method, process, and result are applicable to the trapping issue on all networks, the derivation of Eq. (14) requires inverting the matrix with an order . Since the computation of inverting the matrix puts heavy demands on time and memory for large networks, by using Eq. (14) we can directly compute only for networks with small size. In particular, by applying the approach of inverting the matrix, it appears very difficult, even impossible, to get the exact dominating scaling of characterizing the efficiency of the trapping problem. Therefore, it is of significant practical importance to seek an alternative method of computing even for specific networks, which is able to reduce the computational effort of the method of inverting the matrix.

## Iv Scalings of the MFPT for random walks on modular networks with multiple traps

Here we study two particular trapping problems defined in the modular scale-free networks . We first address the case that traps are located at all peripheral nodes; then we consider the case that traps are fixed on those nodes farthest from the main hub. We will show that the special recursive construction of the modular scale-free networks and the particular selections made for the trap locations allow for an analytical treatment of the MFPT to the traps.

### iv.1 Determination of intermediate variables

Prior to studying the MFPT to the traps, we first define some intermediate variables and determine their values. We denote by the FPT for a walker starting from an arbitrary peripheral node of to visit the hub for the first time and by the FPT spent by a particle initially located at the hub to first visit any peripheral node. In Appendix A, we give detailed derivations for and , which read

(15) |

and

(16) |

respectively.

Equations (15) and (16) are very useful for the following derivation of the exact formula for the MFPT to the targets. We note that Eqs. (15) and (16) have also been derived in Ref. ZhLiGoZhGuLi09 () by using the technique of generating functions Wi94 (), but the approach used here is different from and relatively easier than the previous one.

### iv.2 Exact solution to the MFPT for random walks with traps located at peripheral nodes

After obtaining the intermediate quantities, we are now in a position to consider random walks on networks with all the peripheral nodes being occupied by traps. Our goal in this case is to determine the MFPT denoted by , which is the average of the FPT for a walker originating from a node in to first visit any target over all starting points including the traps. In order to find , we introduce another quantity defined as the FPTs of all nodes to the hub. From the structure of the networks, we can easily establish the following recursive relations for and :

(17) |

and

(18) |

Equations (17) and (18) can be rewritten as

(19) |

and

(20) |

From Eq. (19), we further have

(21) |

which, together with Eqs. (19) and (20), yields

(22) | |||||

that is,

(23) |

Substituting Eqs. (15) and (16) into Eq. (23) and considering the initial condition , we can solve inductively Eq. (23) to obtain the following rigorous expression:

(24) | |||||

Plugging this result for into Eq. (18) and using the initial condition , Eq. (18) is solved to yield

(25) | |||||

To confirm our analytic formulas, we have compared them with the numerical values from the method of inverting the matrix provided by Eq. (14); see Fig. 2. For various values of and , the results for and obtained separately from Eqs. (24) and (25) are in complete agreement with those from Eq. (14). This agreement serves as a mutual test of our numerical solution and analytical formulas, providing an important evidence of the validity of Eqs. (14), (24) and (25).

We proceed to represent and as functions of network size to obtain their dependence on . From we have and , which enables us to recast Eqs. (24) and (25) in terms of as

(26) | |||||

and

Equations (26) and (IV.2) imply that in the limit of large network size (i.e., ), both and grow asymptotically as power-law functions of network size with the same exponent :

(28) |

and

(29) |

Obviously, the exponent is smaller than 1, showing that both and scale sublinearly with the network size.

We note that the scaling in Eq. (28) has been previously derived in Ref. ZhLiGoZhGuLi09 () by using the theory of generating functions, the computation process of which is a little complex. Equation (28) shows that when the hub node is considered an immobile trap, the trapping efficiency is high (even the highest among all networks TeBeVo09 (); ZhAlHoZhCh11 ()), which can be elaborated as follows. In the average distance between the hub and other nodes is only half of the average distance between all pairs on nodes ZhLi10 (), suggesting that the hub node is spatially closer than any other node. In contrast, the degree of the hub node is the highest, which is why the MFPT to the hub is very low.

### iv.3 Behavior of the MFPT for random walks with traps fixed at farthest nodes

We now focus on the trapping problem with traps being placed on the nodes farthest from the main hub, which are expected be more difficult to visit compared with the peripheral nodes TeBeVo11 ().

#### iv.3.1 Related definitions and quantities

In the maximum value of the distance from the main hub to other nodes is . We let denote the set of those nodes in at a distance from the main hub of , hereafter called the farthest nodes of , and denote the cardinality (number of elements in a set) of . By construction, is composed of a primal and copies of , denoted separately by (). For , its farthest nodes are exactly its peripheral nodes; for , its farthest nodes correspond to the hub nodes of all . Proceeding analogously, for the farthest nodes of must belong to all subgraphs , and the farthest nodes of the primal central subgraphs (i.e., ) forming constitute . Thus, we have

(30) |

Considering and , the recursive relation can be solved to obtain

(31) |

Next we concentrate on the MFPT from the hub to the farthest nodes in , which will be denoted by henceforth since, as will be shown, it has the same scaling as that of the average of MFPTs to the farthest nodes , taken over all starting points. For a convenient description of the computation for the MFPT to the farthest nodes, we introduce more variables. For , let and express the sets of the main hub and peripheral nodes, respectively. In addition, for those nodes of that belong to , we can further classify them in the following way. Let () denote the set of those local hubs that are directly connected to classes of local peripheral nodes in and () stand for the set of the local peripheral nodes whose neighbors are different local hubs belonging to . It is easy to verify that the respective degrees of nodes in and are and , respectively.

#### iv.3.2 Exact solution to the MFPT from the hub to farthest nodes

According to the structure of , for a walker starting from the main hub, in order to reach the farthest nodes, it should follow the path . Then it is natural to define the following quantities. Let and represent, respectively, the FPT from a node in to any of its neighboring nodes in and the FPT from a node in to any of its neighbors belonging to . In Appendix B, we report the derivation for and , the exact expressions for which are

(32) | |||||

and

(33) | |||||

respectively.

Using the obtained intermediate quantities, we can derive an exact formula for . We distinguish two cases: (i) is odd and (ii) is even. For odd we have

(34) | |||||

By plugging Eqs. (32) and (33) into Eq. (34) and doing some algebra, we find a closed-form solution to given by

(35) | |||||

When is even, it is not difficult to reach the following expression:

(36) | |||||

Inserting Eqs. (32) and (33) into Eq. (36), after some algebra, the explicit expression for is obtained, which reads

(37) | |||||

To check the validity of Eqs. (35) and (37), we also compute numerically by using the approach of inverting the related matrix; see Eq. (13). The results obtained by analytical and numerical methods completely agree with each other. The comparison is shown in Fig. 3. Equation (35), together with Eq. (37), indicates that for large networks, i.e., ,

(38) |

with the exponent smaller than 1.

Thus far we have found the rigorous formula for the MFPT to farthest nodes in and its dependence on network size . We stress that the analytical computation for the MFPT to the farthest nodes that average all starting points in is rather lengthy and awkward. However, it is easy to infer that when is large enough, the dominant term of also increases as a power-law function of network size with an exponent identical to that of , which can be understood from the following heuristic explanation. Note that consists of subgraphs, which are copies of . For those nodes in the central subgraph, their MFPT to the farthest nodes is equal to , the dominant term of which is ; for nodes in each of the fringe subgraphs (), their MFPT to the farthest nodes is identical but smaller than . Hence, for all nodes in , the dominating term of the MFPT is proportional to but its prefactor may be different from that of .

### iv.4 Result analysis

Equations (28) and (38) show that when traps are positioned at several particular nodes, the MFPTs to the target node are very small, which scale sublinearly with the network order. When either peripheral nodes or farthest nodes are occupied by traps, the characteristic exponent or is a decreasing function of : When the parameter increases from 3 to , both and drop and are close to zero. Therefore, the efficiency of the random-walk process is reliant on : The larger the parameter , the more efficient the random-walk process. The fact that both trapping processes are very efficient demonstrates that the modular scale-free networks being studied exhibit an efficient configuration for random walks with traps positioned at certain given nodes.

In contrast, for each given parameter , is smaller than , which implies that when traps are located at peripheral nodes, the trapping efficiency is higher than that of the case when traps are placed on farthest nodes. Thus, the two trapping processes defined on the networks under consideration display rich behavior in the context of MFPTs to the traps. The difference between and shows that the number and location of traps sensitively affect the behavior of random walks on the modular scale-free networks.

Actually, the intrinsic structure of the modular scale-free networks is responsible for the high efficiency of random walks performing on them with certain nodes being occupied by traps. In these networks, there are many small highly integrated clusters, which group into a few larger but less compact modules linked by local hub nodes; see Fig. 1. These relatively large modules combine to form even larger and fewer groups, which are further joined to shape a fine modular and scale-free architecture, a topology that accounts for the fast diffusion phenomenon in .

In the case that traps are placed on peripheral nodes, when a particle originates from a node in a duplicate (an element of ), it will either be directly trapped by one of the traps or jump to local hub nodes in a few steps. These local hubs play a bridge role linking different small modules together. After arriving at local hub nodes, the particle can be easily trapped in a short time. In contrast, if the particle starts off from a node in the original (the central part of ), it will easily visit local hub nodes or the hub first, through which it can find the way to one of the traps quickly. Thus, the particle can drop into the traps very fast wherever it starts to jump, which can be understood from the above heuristic argument based on the inherent structure of the considered networks.

When the traps are fixed on the farthest nodes, to find a garget, the walker must first visit the local hubs and local peripheral nodes of a larger and sparser cluster, starting from which it continues to arrive at the local hubs and local peripheral nodes of smaller and denser groups. From Eqs. (32) and (33) we know that the expected time between local hub nodes and local peripheral nodes in inner subgraphs rely on their size or deepness (i.e., ): The smaller the value of , the smaller the size of inner subgraphs, and the higher the expected time. This can account for the main reason the farthest nodes are more difficult to reach than the peripheral nodes.

## V Conclusions

We have studied the random-walk dynamics on a family of modular scale-free networks with multiple traps, which exhibit remarkable characteristics observed for various real-life networks, such as social and biological networks. We first deduced a general formula for the MFPT to the traps in a generic network, which is expressed in terms of several elements of a matrix associated with the trapping problem. Then we studied the MFPT for two trapping issues on the studied networks with two different arrangements of targets. In the first case, peripheral nodes are treated as traps; in the second case, farthest nodes work as traps.

For the two trapping problems, we studied both numerically and analytically the MFPT to traps, the results of which are compatible with each other. Our results show that in both cases, the MFPT varies as a power-law function of network size with the exponent depending on the parameter , which is lower than 1 in the full range of . Thus, the studied networks display an efficient architecture in favor of diffusion. Moreover, we demonstrated that, compared with the second case, the diffusion is faster in the first case, which indicates that the transport efficiency relies on the number and location of the absorbing nodes. We also showed that the modular topology, together with the scale-free behavior, is responsible for the quick diffusion processes, as well as the scaling difference of the MFPT for the two trappings running on the networks addressed. We expect that our work can provide insight into designing networks with a structure in favor of diffusion. Finally, it should be mentioned that the method developed here applies only to very specific sets of traps and is hard to generalize to other sets of traps.

### Acknowledgment

This research was supported by the National Natural Science Foundation of China under Grant No. 61074119.

## Appendix A Derivation of and

According to the particular structure of the networks, for any , the two quantities and obey the following recursion relations:

(39) | |||||

and

(40) | |||||

The three terms on the right-hand side (rhs) of Eq. (39) can be explained as follows. The first term is based on the fact that the walker takes only one time step to first reach the hub. The second term describes the process by which the particle first jumps to one of its neighbors belonging to in one time step and then takes more steps to first get to the target node. The last term accounts for the fact that the walker first makes a jump to a local hub node belonging to , then takes time steps, starting off from the local hub, to reach any node in , and continues to jump more steps to reach the target node for the first time.

Analogously, the two terms on the rhs of Eq. (40) are based on the following two processes. The first term describes the fact that the walker, starting from the hub, requires only one time step to hit a peripheral node. The second term explains such a process that the walker, starting off from the hub, first jumps to a local peripheral node belonging to in one time step, then makes jumps to the hub, and proceeds to any node in , taking more time steps.

After merging similar items, Eqs. (39) and (40) can be rewritten as

(41) |

and

(42) |

respectively. Equations (41) and (42) lead to

(43) |

and

(44) |

According to Eq. (44), we obtain

(45) | |||||

where the relation provided in Eq. (43) was used. Applying the initial condition , we solve Eq. (45) to obtain

(46) |

Inserting the result for into Eq. (43), we arrive the exact formula for given by

(47) |

## Appendix B Derivation of and

For the two quantities and , the following relations hold:

(48) | |||||

and

(49) | |||||

Equation (48) can be elaborated as follows. Originating from a node in , the particle can jump to one of the neighboring nodes belonging to , from which it continues to jump steps to first visit a target; this is accounted for by the first term on the rhs. Alternatively, the walker can go to a local hub belonging to , then takes time to reach a neighbor in , and proceeds to bounce steps to hit a target for the first time, this process is explained by the second term. The third term describes the process by which the walker goes directly to a target node. The fourth term represents the process that the walker first jumps to a neighbor belonging to , makes a move returning to a node in , and then walks continuously in time to arrive at a destination node. Finally, the last sum term explains the fact that the particle goes to a local hub in () from which it takes an average time to return to one of its neighbors in , and then moves on average steps to get to a target. Analogously, we can explain Eq. (49).

## References

- (1) R. Albert and A.-L. Barabási, Rev. Mod. Phys. 74, 47 (2002).
- (2) S. N. Dorogovtsev and J. F. F. Mendes, Adv. Phys. 51, 1079 (2002).
- (3) M. E. J. Newman, SIAM Rev. 45, 167 (2003).
- (4) S. Boccaletti, V. Latora, Y. Moreno, M. Chavez and D.-U. Hwanga, Phys. Rep. 424, 175 (2006).
- (5) A.-L. Barabási and R. Albert, Science 286, 509 (1999).
- (6) M. Girvan and M. E. J. Newman, Proc. Natl. Acad. Sci. U.S.A. 99, 7821 (2002).
- (7) E. Ravasz, A. L. Somera, D. A. Mongru. Z. N. Oltvai, and A.-L. Barabási, Science 297, 1551 (2002).
- (8) E. Ravasz and A.-L. Barabási, Phys. Rev. E 67, 026112 (2003).
- (9) F. Chung and L. Lu, Proc. Natl. Acad. Sci. U.S.A. 99, 15879 (2002).
- (10) R. Cohen and S. Havlin, Phys. Rev. Lett. 90, 058701 (2003).
- (11) L. Danon, J. Duch, A. Diaz-Guilera, and A. Arenas, J. Stat. Mech.: Theory Exp. (2005) P09008.
- (12) G. Palla, I. Derényi, I. Farkas, and T. Vicsek, Nature (London) 435, 814 (2005).
- (13) M. E. J. Newman, Proc. Natl. Acad. Sci. U.S.A. 103, 8577 (2006).
- (14) S. Fortunato, Phys. Rep. 486, 75 (2010).
- (15) W. H. Zhan, Z. Z. Zhang, J. H. Guan, and S. G. Zhou, Phys. Rev. E 83, 066120 (2011).
- (16) C. E. Perkins, Ad Hoc Networking, 1st ed. (Addison-Wesley, Reading, MA, 2001).
- (17) A. Y. Wu, M. Garland, and J. Han, in KDD ’04: Proceedings of the Tenth ACM SIGKDD International Conference on Knowledge Discovery and Data Mining (ACM, New York, 2004), pp. 719-724.
- (18) S. N. Dorogovtsev, A. V. Goltsev and J. F. F. Mendes, Rev. Mod. Phys. 80, 1275 (2008).
- (19) S. Havlin and D. ben-Avraham, Adv. Phys. 36, 695 (1987).
- (20) R. Metzler and J. Klafter, Phys. Rep. 339, 1 (2000).
- (21) R. Metzler and J. Klafter, J. Phys. A 37, R161 (2004).
- (22) R Burioni and D Cassi, J. Phys. A 38, R45 (2005).
- (23) I. M. Sokolov, J. Mai, and A. Blumen, Phys. Rev. Lett. 79, 857 (1997).
- (24) S. A. Pandit and R. E. Amritkar, Phys. Rev. E 63, 041104 (2001).
- (25) J. D. Noh and H. Rieger, Phys. Rev. Lett. 92, 118701 (2004).
- (26) O. Bénichou, M. Coppey, M. Moreau, P.-H. Suet, and R. Voituriez, Phys. Rev. Lett. 94, 198101 (2005).
- (27) V. Sood, S. Redner, and D. ben-Avraham, J. Phys. A 38, 109 (2005).
- (28) E. M. Bollt and D. ben-Avraham, New J. Phys. 7, 26 (2005).
- (29) L. K. Gallos, C. Song, S. Havlin, and H. A. Makse, Proc. Natl. Acad. Sci. U.S.A. 104, 7746 (2007).
- (30) P. Barthelemy, J. Bertolotti, and D. S. Wiersma, Nature (London) 453, 495 (2008).
- (31) A. Baronchelli, M. Catanzaro, and R. Pastor-Satorras, Phys. Rev. E 78, 011114 (2008).
- (32) A. Kittas, S. Carmi, S. Havlin, and P. Argyrakis, EPL 84, 40008 (2008).
- (33) Z. Z. Zhang, Y. C. Zhang, S. G. Zhou, M. Yin, and J. H. Guan, J. Math. Phys. 50, 033514 (2009).
- (34) C. P. Haynes and A. P. Roberts, Phys. Rev. Lett. 103, 020601 (2009).
- (35) Z. Z. Zhang, J. L. T. Alafate, B. Y. Hou, H. J. Zhang, and G. R. Chen, Eur. Phys. J. B 84, 691 (2011).
- (36) G. H. Weiss, Aspects and Applications of the Random Walk (North Holland, Amsterdam, 1994).
- (37) B. H. Hughes, Random Walks and Random Environments (Clarendon Press, Oxford, 1995), Vols. 1 and 2.
- (38) F. Jasch and A. Blumen, Phys. Rev. E 63, 041108 (2001).
- (39) M. F. Shlesinger, Nature (London) 443, 281 (2006).
- (40) J. Whitmarsh and J. A. Govindjee, in Concepts in Photobiology: Photosynthesis and Photomorphogenesis, edited by G. Singhal, G. Renger, S. Sopory, K.-D. Irrgang, and J. A. Govindjee (Narosa, New Delhi, 1999), pp. 11-51.
- (41) F. Fouss, A. Pirotte, J. M. Renders, and M. Saerens, IEEE Trans. Knowl. Data Eng. 19, 355 (2007).
- (42) E. W. Montroll, J. Math. Phys. 10, 753 (1969).
- (43) J. J. Kozak and V. Balakrishnan, Phys. Rev. E 65, 021105 (2002).
- (44) J. J. Kozak and V. Balakrishnan, Int. J. Bifurcation Chaos Appl. Sci. Eng. 12, 2379 (2002).
- (45) Z. Z. Zhang, B. Wu, H. J. Zhang, S. G. Zhou, J. H. Guan, and Z. G. Wang, Phys. Rev. E 81, 031118 (2010).
- (46) E. Agliari, Phys. Rev. E 77, 011128 (2008).
- (47) Z. Z. Zhang, Y. Lin, S. G. Zhou, B. Wu, and J. H. Guan, New J. Phys. 11, 103043 (2009).
- (48) A. Garcia Cantú and E. Abad, Phys. Rev. E 77, 031121 (2008).
- (49) Z. Z. Zhang, Y. Qi, S. G. Zhou, W. L. Xie, and J. H. Guan, Phys. Rev. E 79, 021127 (2009).
- (50) Z. Z. Zhang, S. G. Zhou, W. L. Xie, L. C. Chen, Y. Lin, and J. H. Guan, Phys. Rev. E 79, 061113 (2009).
- (51) Z. Z. Zhang, J. H. Guan, W. L. Xie, Y. Qi, and S. G. Zhou, EPL, 86, 10006 (2009).
- (52) E. Agliari and R. Burioni, Phys. Rev. E 80, 031125 (2009).
- (53) V. Tejedor, O. Bénichou, and R. Voituriez, Phys. Rev. E 80, 065104(R) (2009).
- (54) Z. Z. Zhang, W. L. Xie, S. G. Zhou, S. Y. Gao, and J. H. Guan, EPL 88, 10001 (2009).
- (55) Z. Z. Zhang, W. L. Xie, S. G. Zhou, M. Li, and J. H. Guan, Phys. Rev. E 80, 061111 (2009).
- (56) S. Condamin, O. Bénichou, V. Tejedor, R. Voituriez, and J. Klafter, Nature (London) 450, 77 (2007).
- (57) O. Bénichou, B. Meyer, V. Tejedor, and R. Voituriez, Phys. Rev. Lett. 101, 130601 (2008).
- (58) S. Condamin, V. Tejedor, R. Voituriez, O. Bénichou and J. Klafter, Proc. Natl. Acad. Sci. U.S.A. 105, 5675 (2008).
- (59) S. B. Yuste and L. Acedo, Phys. Rev. E 68, 036134 (2003).
- (60) T. A. Witten and L. M. Sander, Phys. Rev. Lett. 47, 1400 (1981).
- (61) T. A. Witten and L. M. Sander, Phys. Rev. B 27, 5686 (1983).
- (62) S. A. Rice, Diffusion-Limited Reactions, edited by C. H. Bamford, C. F. H. Tipper, and R. G. Compton, Comprehensive Chemical Kinetics Vol. 25 (Elsevier, Amsterdam, 1985).
- (63) M. Hinczewski, Phys. Rev. E 75, 061104 (2007).
- (64) S. Maslov, Nature Phys. 3, 18 (2007).
- (65) A. Diaz-Guilera, J. Phys. A 41, 224007 (2008).
- (66) Z. Z. Zhang, Y. Lin, S. Y. Gao, S. G. Zhou, J. H. Guan, and M. Li, Phys. Rev. E 80, 051120 (2009).
- (67) E. Bullmore and O. Sporns, Nature Rev. Neurosci. 10, 186 (2009).
- (68) S. Redner, A Guide to First-Passage Processes (Cambridge University Press, Cambridge, 2001).
- (69) S. N. Dorogovtsev, A. V. Goltsev, and J. F. F. Mendes, Phys. Rev. E 65, 066122 (2002).
- (70) S. Jung, S. Kim, and B. Kahng, Phys. Rev. E 65, 056101 (2002).
- (71) F. Comellas, J. Ozón, and J.G. Peters, Inf. Process. Lett. 76, 83 (2000)
- (72) Z. Z. Zhang, L. L Rong, and C. H. Guo, Physica A 363, 567 (2006).
- (73) M. Hinczewski and A. N. Berker, Phys. Rev. E 73, 066126 (2006).
- (74) H. D. Rozenfeld, S. Havlin, and D. ben-Avraham, New J. Phys. 9, 175 (2007).
- (75) Z. Z. Zhang, S. G. Zhou, T. Zou, L. C. Chen, and J. H. Guan, Eur. Phys. J. B 60, 259 (2007).
- (76) Z. Z. Zhang, S. G. Zhou, L. J. Fang, J. H. Guan, and Y. C. Zhang, EPL 79, 38007 (2007).
- (77) S. Boettcher, B. Gonçalves, and H. Guclu, J. Phys. A 41, 252001 (2008).
- (78) F. Comellas and A. Miralles, Physica A. 388, 2227 (2009).
- (79) K. L. Morrow, T. Rowland, and C. M. Danforth, Phys. Rev. E 80, 016103 (2009).
- (80) Z. Z. Zhang, J. H. Guan, B. L. Ding, L. C. Chen, and S. G. Zhou, New J. Phys. 11, 083007 (2009).
- (81) A. Miralles, F. Comellas, L. C. Chen, and Z. Z. Zhang, Physica A. 389, 1955 (2010).
- (82) J. D. Noh, Phys. Rev. E 67, 045103(R) (2003).
- (83) J. D. Noh and H. Rieger, Phys. Rev. E 69, 036111 (2004).
- (84) A.-L. Barabási, E. Ravasz, and T. Vicsek, Physica A 299, 559 (2001).
- (85) K. Iguchi and H. Yamada, Phys. Rev. E 71, 036144 (2005).
- (86) Z. Z. Zhang, Y. Lin, S. Y. Gao, S. G. Zhou, and J. H. Guan, J. Stat. Mech. (2009) P10022.
- (87) E. Agliari, R. Burioni, and A. Manzotti, Phys. Rev. E 82, 011118 (2010).
- (88) Z. Z. Zhang and Y. Lin, J. Stat. Mech. (2010) P12017.
- (89) D. J. Watts and H. Strogatz, Nature (London) 393, 440 (1998).
- (90) J. G. Kemeny and J. L. Snell, Finite Markov Chains (Springer, New York, 1976).
- (91) D. Aldous and J. Fill, http://www.stat.berkeley.edu/ aldous/RWG/Chap2.pdf
- (92) A. Bar-Haim and J. Klafter, J. Chem. Phys. 109, 5187 (1998).
- (93) Note that the definition of MFPT slightly differs from that used in some previous work where the average of is taken over all nodes other than the traps. However, when the concentration of traps (the ratio of the trap number to network size ) is very small, both definitions lead to the same behavior of MFPT .
- (94) H. S. Wilf, Generatingfunctionology, 2nd ed. (Academic, London, 1994).
- (95) V. Tejedor, O. Bénichou, and R. Voituriez, Phys. Rev. E 83, 066102 (2011).