arXiv:2004.00575v2 [cond-mat.stat-mech] 31 Aug 2020 oe ehv h rcinldffso nudrce net- undirected on diffusion fractional be- mentioned the walks have random we nonlocal fore, the to between addition In hops long-range 19]. define [18, use nodes to can distant information the strategies of this nonlocal architecture and all total play [15– into the neighbors comes cases, nearest network other its in of whereas one walker to 17], the way, node this one in from require and moves only node each them been of of have information some networks local decades, Several last on the examples. strategies in few walk introduced a random mention of to information types just useful [14], spread- extracting data for epidemic from algorithms [7–11], 13], human [12, settlements of ing urban understanding the in ac- [6], mobility and damage aging [3–5], of of web the cumulation field on topics searching challenging interdisciplinary and ranking with and like differ- connections important following to due an network research prob- the is that of diffusion strategies walker nodes random the ent a the particular, of visiting dynamics In hops among the to systems associated 2]. lem computer [1, and others in in many social applications impact important biology, significant with a physics, engineering have and networks science on place taking ∗ aperezr@fisica.unam.mx h td n nesadn fdnmclprocesses dynamical of understanding and study The olclbae admwlsadfatoa rnpr ndi on transport fractional and walks random biased Nonlocal ASnmes 97.c 54.b 25.r 05.60.Cd 02.50.-r, ar 05.40.Fb, connected 89.75.Hc, strongly numbers: more PACS structures. any these not on is network bias directed with rand a walker on when random dynamics directe the fractional and different a rings on of on walks transport biased random t the ergodic apply works, We nonlocal anal strategies. and We walk local random matrix. the eigenvector characterize Laplacian and to eigenvalues the times respective of t their powers with and process of matrices Markovian terms t discrete-time approach in a general expressed as a dynamics present the We ing networks. directed of matrix nti ae,w td olclrno aksrtge gene strategies walk random nonlocal study we paper, this In eatmnod ´sc,Fcla eCeca,Universida Ciencias, de F´ısica, Facultad de Departamento .INTRODUCTION I. NSUR79, lc use,722Prscdx0,France 05, cedex Paris 75252 Jussieu, place 7190,4 UMR CNRS nttt eFıia nvria ainlAtooad M´ Aut´onoma de Nacional F´ısica, Universidad de Instituto praoPsa 034 10 iddd M´exico, M´exico de Ciudad 01000 20-364, Postal Apartado iddUiestra iddd ´xc 41,M´exico M´exico 04510, de Ciudad Universitaria, Ciudad obneUieste ntttJa eRn d’Alembert, Rond le Jean Institut Universit´e, Sorbonne Dtd etme ,2020) 2, September (Dated: .Pizarro-Medina A. ..Michelitsch T.M. networks ..Riascos A.P. we n aro oe xssapt ffiielnt.A length. finite ( of nodes path of ‘ a called is exists graph nodes directed of pair any tween nudrce rp srfre oa ‘ as to referred is the graph of undirected nodes decay. An power-law all asymptotic between characteristic connections with network with network a to sgnrtdb rcinlpwr fteLpainma- Laplacian the of a powers generates fractional trix) by that function generated matrix the is (i.e. matrix Laplacian fractional ti motn ontc hti a in that notice to network important of is case It the to formalism [31]. networks this weighted all directed of networks extension on the dynamics require Other nonlocal of [30]. structures applications data potential algo- of learning classification network the semi-supervised a for rithms in of implemented Laplacian re- been fractional networks; has of on concept processes the the Beyond dynamical cently, nonlocal [20]. that networks of but small-world study trees in offers average and evident that large is lattices also like result with nodes a networks between nodes, in distances the advantage visit significant improv- to a [21] capacity network gener- the the possible to on all ing mechanism nodes of two This information connecting matrix paths the Laplacian 28]. combines the 27, nonlocality of [20, ate powers network of the terms of nodes in between dy- defined L´evy probabilities like are a transition the with where associated namics process a [20–29], works nohrwrsaynd a erahdfo n other any from reached be can node any words other in ∗ ,tema rtpsaetmsadglobal and times passage first mean the s, i praht h td fparticular of study the to approach his n na in and ewrs eepoecruatnet- circulant explore we networks; d ffcso roiiyls hc occur which loss ergodicity of Effects asto rbblte ewe nodes between probabilities ransition nlzn hs taeisb defin- by strategies these analyzing o mntok.W td h efficiency the study We networks. om ainlAtooad M´exico, Aut´onoma de Nacional d ij lodiscussed. also e ae ihtefatoa Laplacian fractional the with rated hr r ietdpaths directed are there ) z h lmnso h transition the of elements the yze ul once topology connected fully togycnetddrce graph directed connected strongly exico, togyconnected strongly once undirected connected i → connected j ffraypair any for if ’ corresponding rected and j → fbe- if ’ i or , the 2 node by a finite number of steps (see Ref. [31] for defi- dom walker to visit all the nodes or reach a particular nitions and outline of properties). Connected undirected target on the network. and strongly connected directed graphs fulfill the condi- tion of aperiodic ergodicity. Conversely aperiodic ergodic graphs always are either connected undirected networks II. FRACTIONAL LAPLACIAN OF DIRECTED or strongly connected directed graphs. NETWORKS The fractional Laplacian contains the complete informa- In this section, we introduce a generalization of the tion on the topology of the network. For an outline how fractional Laplacian of undirected networks (see Refs. long-range interactions of asymptotic power-law decay in [20, 27, 28]) to a general class of directed weighted net- harmonic systems modify their spectral properties and works. In terms of this operator, we define transition universal features, we refer to [32]. The spectral dimen- probabilities of a Markovian random walker associated sion of general networks is analyzed in the seminal paper to the biased fractional transport on networks. [33] and applications of this approach in spin models on graphs is outlined in the article [34], the Laplacian spec- We consider directed weighted networks with N nodes trum of simplicial complexes using the renormalization i = 1,...,N. The topology of the network is described group is discussed in Ref. [35]. For a general analysis of by an adjacency matrix A with elements Aij = 1 if there the spectral dimension of (unbiased) L´evy flights in the is an edge between the nodes i and j and A = 0 oth- Rd ij we invite the reader to consult Chapter 8 in Ref. [28]. erwise. In addition to the network structure, we have a N N matrix of weights Ω with elements Ωij 0. In the present paper, we explore the dynamics of a ran- × ≥ dom walker with transition probabilities defined in terms The matrix Ω can include information of the structure of the elements of the fractional Laplacian matrix in di- of the network or incorporate additional data describing rected networks. In the first part, general definitions the flow capacity of each link [6, 36, 37]. In the simplest and properties of the fractional transport, ergodicity, and case, Ω coincides with the adjacency matrix A. emergence of nonlocality on strongly connected directed Since the matrix of weights in general is not symmetric, networks are discussed. The formalism generalizes dif- we define two types of degrees associated to each node. ferent results and techniques developed in the context of First, we have the in-degree given by the fractional Laplacian of an undirected graph [28]. In contrast to the unbiased transport defined through sym- N (in) metric weighted matrices [36], in the directed case the ki = Ωli. (1) eigenvalues of fractional transition matrices are complex l=1 X numbers. We explore different types of directed struc- tures, especially circulant directed networks such as di- This degree determines the total flow to reach the node rected rings, but also random directed networks of the i from all the nodes. In a similar way, we have the out- Erd˝os-R´enyi type. In the case of directed rings, we show degree analytically the connection of fractional dynamics with N nonlocal random walks similar to L´evy flights where the (out) k = Ωil, (2) property of strong connectivity generates aperiodic er- i l=1 godicity in the fractional walk. We demonstrate that in X directed networks which are not strongly connected the that quantifies the total flow from the node i to all fractional Laplacian has zero matrix elements, and con- the nodes in the network. Without loss in the gener- versely if the fractional Laplacian matrix has uniquely ality of the formalism, we assume also that Ωii = 0 for non-zero entries, the directed graph is strongly connected i = 1, 2,...,N. In the following, we consider connected where the resulting walk is aperiodic ergodic. (out) directed networks for which ki > 0 for all the nodes. We also analyze the efficiency of fractional random walk In terms of the matrix of weights, we define the Laplacian strategies that emerge in directed networks using mean- matrix L with elements i, j, given by first passage times and global times that characterize the transport. The implementation of these measures allows (out) Lij = ki δij Ωij (3) identifying cases where the combination of biased trans- − port and nonlocality is an inefficient strategy to explore where δij denotes the Kronecker’s delta. Equation (3) a network; this particular result is in contrast with the is a generalization of the Laplacian matrix for binary efficiency on undirected networks for which the fractional undirected networks [38–40], to include the possibility dynamics always improve the speed of the exploration in of weights in the connections and asymmetry in the flow comparison with a local random walker [28]. The gen- on some lines, for these particular connections Ωij =Ωji. eral approach introduced reveals several cases where the It is noted that dynamical processes in directed networks6 combination of nonlocal displacements and the bias gen- have a greater variety than in the undirected case. For erated by the directions of lines produce a global effect example, for the matrix L, the diagonal could be defined that can either reduce or improve the efficiency of a ran- in terms of the in-degree in Eq. (1) producing a different 3
FIG. 1. (Color online) Fractional Laplacian of directed networks with N = 10 nodes. Graphs formed with two directed rings and (a) a directed edge from 5 to 6 and (b) a symmetric link connecting these two nodes. For each network we present the Laplacian matrix (local limit with γ = 1) and the fractional Laplacian Lγ with γ = 0.5, the values of entries i, j are codified in the colorbar. For the fractional Laplacian γ = 0.5, we evaluate numerically Eq. (5). process. Our choice in Eq. (3) is motivated by diffusive Jordan canonical form we refer to the article [31]. transport and the effect of nonlocality, similar nonlocal In Fig. 1, we show the fractional Laplacian Lγ for di- effects have been found in human mobility in different rected networks with N = 10. In this case, we choose Ω types of transport described by directed networks [9–11]. as the adjacency matrix of the network. These graphs are On the other hand, in the context of the fractional dif- formed by two directed rings (the first ring with nodes fusion on networks is introduced the fractional Lapla- 1 to 5 and the second one with nodes 6 to 10) and a cian matrix Lγ , where γ is a real number (0 <γ< 1). directed line connecting the two rings (see Fig. 1(a)) or The resulting operator models the fractional dynamics a symmetric link between nodes 5 and 6 (as depicted in on general networks [20, 27, 28]. Using Dirac’s nota- Fig. 1(b)). The network in Fig. 1(a) is not (strongly) tion for the eigenvectors, we have a set of right eigen- connected, see for example that there are no paths start- N vectors Ψj that satisfy the eigenvalue equation ing from one of the nodes 6 to 10 ending in nodes 1 to {| i}j=1 L Ψj = µj Ψj for j =1,...,N. With this information, 5, in these cases the distances between nodes are infinite. | i | i we define the matrix Q with elements Qij = i Ψj and In this directed graph, we observe that the respective el- h | i the diagonal matrix Λ = diag(µ ,µ ,...,µN ). These ements of the fractional Laplacian are also null (i.e. the 1 2 γ matrices satisfy LQ = QΛ, therefore block (L )rl = 0 for r =6,..., 10 and l =1,..., 5), this γ is a consequence of the fact that (L )ij incorporates in- 1 L = QΛQ− , (4) formation of all the possible paths connecting i with j, γ then when distance dij , (L )ij = 0 for 0 < γ 1 1 1 → ∞ ≤ where Q− is the inverse of Q. Using the matrix Q− , (see Ref. [28]). In contrast, the structure in Fig. 1(b) is ¯ N we define the set of left eigenvectors Ψi i=1 with com- strongly connected, hence, there is a path of finite length 1 { } ponents Ψ¯ i j = (Q− )ij . Therefore connecting any pair of nodes of the network, in this case, | γ we see that all the elements of (L )ij are non-null, drawn N for γ =0.5, as a result of the full connectivity of that net- γ γ 1 γ L = QΛ Q− = µ Ψm Ψ¯ m , (5) work. m | i m=1 X Back to the general case with strongly connected net- where Λγ = diag(µγ ,µγ ,...,µγ ) for 0 < γ 1. It is works, due to the asymmetry of Ω, the eigenvalues µl 1 2 N ≤ can take complex values. However, as a consequence of sufficient here for our aims to consider uniquely the diag- N onalizable cases. For an outline where the Laplacian has l=1 Lil = 0, the definition of the out-degree in Eq. P 4
(γ) (2) and the conditions Lii > 0, Lij 0 for i = j, the matrix wi j > 0 ( i = j) of the fractional walk are γ ≤ 6 fractional Laplacian L of a directed weighted network → ∀ 6 (γ) strictly positive (with wi i = 0 per construction) and has real entries and fulfills the following properties for → hence W(γ) fulfills the condition of aperiodic ergodicity. 0 <γ 1: ≤ It is worth noticing the role of the nonlocality generated (i) For the fractional out-degree, we have by Lγ . If the local random walker (γ = 1) can reach any (γ) γ γ node in the network in a finite number of steps starting k (L )ii = (L )im. (6) γ i ≡ − from any node (ergodic condition), L combines the in- m=i X6 formation of all these trajectories in the directed network to define a new nonlocal process that maintains ergodic- Condition (i) reflects the property that the zero eigen- ity. However, this is not the case when the process with value of the Laplacian matrix is conserved by the frac- γ = 1 is not ergodic as we saw in the example in Fig. 1 tional Laplacian (where the corresponding eigenvector (a). For this reason, in the next part, we maintain our has constant components). discussion only for strongly connected weighted networks (ii) The diagonal elements of Lγ are positive real values; for which W(γ) defines ergodic processes (see Appendix (γ) VII for a proof of aperiodic ergodicity in connected undi- in this way ki > 0 for i =1, 2,...,N. rected and strongly connected directed graphs). (iii) The non-diagonal elements of Lγ are real values sat- γ isfying (L )ij 0 for i = j. See Refs. [27, 28, 31] for Here we consider Markovian memoryless walks on di- a detailed discussion≤ on6 these properties for undirected rected graphs where at each time instant t = 0, 1, 2,... and directed networks, respectively. the fractional random walker makes a jump from one to another node on the network in a process without mem- Considering the following integral representation of the ory. The probability P (t; γ) to start at time t = 0 on fractional Laplacian matrix [27, 28] ij node i and to reach the node j at time t satisfies the master equation [16, 17, 28] γ 1 ∞ γ 1 Lt L = t− − (I e− ) dt, 0 <γ< 1 (7) −Γ( γ) − N − Z0 (γ) Pij (t + 1; γ)= Pim(t; γ)wm j . (9) Γ(1 γ) → ( Γ( γ) = − > 0) shows that the properties (i)- m=1 − − γ X (iii) of L are conserved in the interval of convergence In the following part, we analyze the consequences of γ (0, 1) of Eq. (7) (see Appendix VII for a brief demon- the fractional dynamics defined by Eqs. (8)-(9) on dif- ∈ stration). In this relation I = (δij ) indicates the N N ferent directed weighted networks that include circulant identity matrix and Γ(..) stands for the Gamma-function.× directed networks, biased transport on rings and random We observe that in the fractional interval γ (0, 1) all networks. We explore the characteristics of the eigenval- matrix elements of the fractional Laplacian in∈ Eq. (7) ues of the transition matrix, the probabilities of transi- are strictly non-zero if and only if the directed network tion in Eq. (8), mean first passage times and global times is strongly connected which is true for the graph in Fig. that describe the efficiency of the random walker to reach 1(b), however it is not true for the graph in Fig. 1(a) any node on the network. which is not strongly connected where some elements of the fractional Laplacian matrix are zero. The fractional Laplacian of a strongly connected structure for γ (0, 1) III. FRACTIONAL RANDOM WALKS ON γ γ ∈ DIRECTED CIRCULANT NETWORKS fulfills (L )ii > 0 and (L )ij < 0 (for i = j) where all entries are strictly non-zero. 6 In this section, we explore the fractional transport with The characteristics of the fractional Laplacian matrix al- transition probabilities given by Eq. (8). We analyze low to define the fractional diffusion on directed weighted directed networks defined by matrices of weights Ω with networks as a discrete-time Markovian process deter- a circulant matrix structure. (γ) (γ) mined by a transition matrix W with elements wi j representing the probability to hop from i to j given by→ A. Circulant matrices γ (γ) (L )ij wi j = δij (γ) 0 <γ 1. (8) → − ≤ A circulant matrix C is an n n matrix with the form ki [41, 42] × Above properties (i)-(iii) of the fractional Laplacian ma- c0 cn 1 cn 2 ... c1 trix indeed guarantee stochasticity of the fractional tran- − − c1 c0 cn 1 ... c2 sition matrix (8) in the interval 0 < γ 1. Further − ≤ C = c2 c1 c0 ... c3 , (10) for a strongly connected directed network it follows from . . . . ...... the above consideration that in the fractional interval . . . . . 0 <γ < 1 all off-diagonal elements of the transition cn 1 cn 2 cn 3 ... c0 − − − 5 with elements Cij . Thus, each column has real elements c0,c1,...,cn 1 ordered in such a way that c0 describes − the diagonal elements and Cij = c(i j)mod n. In addi- tion to C, the elementary circulating matrix− E is defined, which has all its null elements except c1 = 1. From E, the integer powers El for l =0, 1, 2,...,n 1 are also circu- − lant matrices with null elements except cl = 1. Therefore, Eq. (10) can be expressed as [42]
2 n 1 C = c0I + c1E + c2E + ... + cn 1E − − n 1 − m = cmE , (11) m=0 X where I = E0 is the n n identity matrix. Furthermore, the relation En = I requires× that the eigenvalues ν of E n satisfy ν = 1; therefore, those eigenvalues are given by FIG. 2. (Color online) Circulant directed networks with [42] N = 10 nodes. We define cyclic structures using circulant adjacency matrices A with particular non-null elements in − i 2π(l 1) Eq. (10). For each network we indicate the coefficients c νl = e n for l =1, . . . , n, (12) i (i = 0, 1,...,N − 1) taking the value 1. Arrows represent n the direction of each edge, and connections including both with i √ 1. The respective eigenvectors Ψm m=1 ≡ − 1 i 2π (l 1)({|m 1)i} directions are represented with a line. have the components l Ψm = e− n − − (see h | i √n Ref. [42] for details). Now, using Eq. (11), the eigenvec- tors Ψl satisfy C Ψl = ηl Ψl , where the eigenvalues ηl | i | i | i two nodes are connected with links in both directions. are given by [42] We represent this particular type of connection with lines between nodes as shown in Figs. 2(d)-(f). n 1 − i 2π (l 1) m ηl = cme n − (13) The advantage in the study of the fractional dynamics m=0 on circulant networks is that we know all the eigenval- X ues and eigenvectors of the Laplacian matrix L, since for l = 1, 2,...,n. This result defines the eigenvalues of this is also a circulant matrix. In addition, circulant net- C in terms of the coefficients c0,c1,...,cn 1. − works are regular with the same fractional out-degree (γ) 1 N γ k = N m=1 µm for all the nodes. Therefore, the γ eigenvalues of the transition matrix W(γ) = I L in B. Directed circulant networks P k(γ) Eq. (8) for a fractional random walker are −
The circulant matrix C defined in Eq. (10) allows us to µγ λ(γ) =1 i i =1, 2,...,N, (14) explore different directed structures with N nodes and an i − k(γ) adjacency matrix A with specific values c1,c2,...,cN 1 − equal to 0 or 1, the result is a network with a periodic where, applying the result in Eq. (13) for a network de- structure. When non-null elements appear in pairs ci fined by A with a set of lines with a particular sequence N 1 and cN i (i = 1, 2,...,N 1) in the adjacency matrix, of values 0 and 1 for the coefficients cm − , the eigen- − − m=1 the structure is an undirected network with symmetric values of the Laplacian matrix L are{ } A. Fractional dynamics on undirected circulant networks has been studied in detail for continuous-time random N 1 − i 2π (l 1) m walks [21, 28], quantum transport [22] and diffusion on µl = c0 cme N − , (15) − m=1 multilayer networks [29]. X In the following, we explore cases when the matrix of N 1 where the diagonal element c0 m−=1 cm is the out- weights Ω that defines the Laplacian in Eq. (3) is Ω = A, degree of each node. ≡ this adjacency matrix is not symmetric with a net direc- P tion in some lines. For example, the particular network The result in Eq. (15) shows that in circulant directed with non-null c1 = 1 (or cN 1 = 1) produces a directed networks the eigenvalues of the Laplacian matrix are ring. In Fig. 2, we illustrate− several directed circulant complex numbers and, as a consequence, the eigenval- γ structures with N = 10 nodes. In Fig. 2(a) we have a ues of W( ) in Eq. (14) are also complex numbers. In directed ring, whereas in Figs. 2(b)-(f) other networks Fig. 3, we illustrate the effect of the fractional parameter (γ) are generated by adding new sets of lines defined with γ on the eigenvalues λi in Eq. (14) for different struc- non-null elements ci (i =1, 2,...,N 1). In some cases, tures with topologies similar to the networks in Fig. 2. − 6
FIG. 3. (Color online) Eigenvalues of the transition matrix W(γ) for circulant directed networks with N = 100 represented in (γ) the complex plane. Each eigenvalue λi , i = 1, 2,...,N, is determined by Eq. (14) with eigenvalues µl given by Eq. (15). In (a)-(f) we specify the non-null coefficients cm = 1 that define the adjacency matrix of the network. For each γ, we have N = 100 points and, the effect of the fractional parameter γ modifies the set of eigenvalues represented with different colors codified in the colorbar in the interval 0 <γ ≤ 1. The studied networks have a similar topology to those presented in Fig. 2.
(γ) We show the real and imaginary parts of λi for directed mitian conjugate. Therefore, the elements of the frac- (γ) I Lγ networks with N = 100 nodes and 0 <γ 1. tional transition matrix W = γ are ≤ − k( ) N γ 2π i N (l 1)(j i) (γ) l=1 µl e − − C. Long-range dynamics and L´evy flights wi j = δij N γ . (17) → − P l=1 µl We have analyzed the spectral properties of circulant P γ With the spectral characteristics of the eigenvalues µl matrices associated to the Laplacian L and fractional (γ) and the respective λ illustrated in Fig. 3, the sums transition probabilities W(γ). In the following part, we l in Eq. (17) in the interval 0 < γ 1 produce well de- (γ) ≤ calculate the probabilities wi j and the relation with the fined probabilities of transition between nodes. In Fig. 4 → distance dij that gives the shortest-path length between (γ) we present the transition probabilities wi j for circulant nodes i and j. Using Eq. (15) and the respective eigen- 4 → vectors of circulant matrices, we have for the fractional networks with N = 10 nodes. Probabilities are pre- Laplacian sented as a function of the distance dij connecting the nodes i and j (dij is the length of the shortest path on N the directed structure). In the case of the directed ring γ γ (L )ij = µl i Ψl Ψl j (γ) 1 γ h | ih | i [Fig. 4(a)], we see the relation wi j d−ij − , a power- l=1 → ∝ X law decay also observed in the large-world network ex- N 1 γ i 2π (l 1)(j i) plored in Fig. 4(b) and exemplified in Fig. 2(b). This = µ e N − − . (16) N l particular relation between transition probabilities and l X=1 distances show an emergent dynamics generated through Here, we use the fact that left eigenvectors in circulant the fractional Laplacian related with a L´evy-like dynam- matrices satisfy Ψ¯ l = ( Ψl )†, where denotes the Her- ics in directed structures. L´evy flights in undirected | i †
7
(γ) FIG. 4. (Color online) Transition probabilities wi→j between two nodes as a function of the distance dij in directed circulant 4 networks with N = 10 nodes. (a) A directed ring defined with cN−1 = 1, (b) a network with cN−1 = cN−2 = 1. The (γ) probabilities wi→j are deduced from the analytical relation in Eq. (17) defined in terms of the eigenvalues in Eq. (15) for (γ) −1−γ γ = 0.25, 0.5, 0.75 and γ = 0.9. Dashed lines represent the inverse power-law relation wi→j ∝ dij .
γ γ networks have been explored in a series of works [18– In particular, (L )ij = 0 for jtriangular matrix where all capacity to explore a network by inducing dynamically entries below the diagonal are null. The infinite ring of the small-world property [18, 20]. Similar long-range our example hence is (unlike the finite ring) not any more strategies have been identified in human mobility [9], strongly connected and hence not any more ergodic, and in the movement of cyclists between stations in bike- the corresponding fractional walk on the infinite ring al- sharing systems in Chicago and New York [10], in taxi lows only jumps into the positive integer-direction (hav- trips in New York City [11] and in the infection spreading ing a triangular transition matrix given in Eq. (44)). through the United States’ highly-connected air travel Indeed this is a strictly increasing walk into the positive network [46]. integer direction and can be identified with the so-called Sibuya walk (See Appendix VIII for some derivations and for a profound analysis of properties consult [47]). D. Infinite directed ring On the other hand, in the limit j i, using the re- lation Γ(n + α) Γ(n)nα for n large,≫ we have for Eq. ≈ Now, we explore the fractional Laplacian matrix and (19) transition probabilities in a directed ring. This network 1 1 is illustrated in Fig. 2(a), in the limit N . In the (L)γ for j i (20) → ∞ ij ∼ Γ( γ) (j i)1+γ ≫ particular case of a directed ring with N nodes, Eq. (16) − − takes the form (we use the value cN 1 = 1) where in this relation γ (0, 1). Mention worthy in − ∈ this relation is the local limit γ 1− which becomes N γ γ 1 i 2π (l 1) i 2π (l 1)(j i) → (L )ij = 1 e− N − e N − − . (18) (j i = x) N − − l=1 γ γ γ 1 X lim (L )ij = lim x− − γ 1− γ 1− −Γ(1 γ) However, in the limit of N large, we can define a contin- → → 2π 2π − uous variable ϕ = N (l 1) and dϕ = N . Therefore, the (21) − d x γ d elements of the fractional Laplacian matrix for an infinite = lim − = δ(x)=0, x large. directed ring are given by γ 1− dx Γ(1 γ) dx → − 2π γ 1 iϕ γ iϕ(j i) The asymptotic result in Eq. (20) shows analytically that (L )ij = (1 e− ) e − dϕ, γ (0, 1] 2π 0 − ∈ the fractional dynamics in the infinite directed ring pro- Z (γ) 1 γ duces transition probabilities w d− − for j i j i+1 sin[π(j i γ)]Γ(j i γ) i j ij = ( 1) − csc(πγ) − − − − → ∝ ≫ − Γ( γ)Γ(j i + 1) with distances dij = j i, this relation is also valid − − for networks with N large− but not necessarily infinite Γ(j i γ) = − − , j i. (19) (see Appendix VIII for a detailed discussion). The be- Γ( γ)Γ(j i + 1) ≥ − − havior observed differs from the undirected ring, where 8
4 FIG. 5. (Color online) Mean first passage time hTij (γ)i as a function of the distance dij in circulant networks with N = 10 nodes. (a) A directed ring defined with cN−1 = 1, (b) an undirected ring with cN−1 = c1 = 1. Numerical values are obtained from Eq. (24) defined in terms of the eigenvalues in Eq. (14). We codified in the colorbar the different values 0 <γ ≤ 1. a similar analysis reveals a L´evy-like dynamics with This is the Kemeny’s constant that quantifies the ca- (γ) 1 2γ pacity of the process to explore the network in regular wi j d−ij − (see Refs. [20, 21, 27]), a result that in→ the∝ general case of undirected n dimensional lattices structures and only depends on the eigenvalues of the (γ) n 2γ − transition matrix W. is wi j dij− − [24, 26, 28]. → ∝ In the fractional transport on circulant networks, the eigenvectors of L, the fractional Laplacian Lγ and the E. Efficiency of the fractional transport in (γ) circulant structures transition probability matrix W coincide, since all of them are circulant matrices. Hence, i φl φ¯l j = 1 i 2π (l 1)(j i) h | i | i Ψl Ψl j = e N − − . Then, for the fractional In connected circulant networks, each node has the h | ih | i N same fractional degree k(γ). In addition, if the struc- transport in circulant networks, Eq. (22) takes the form (i = j) ture is strongly connected we have only one eigenvalue 6 (γ) λ1 = 1 and; in this particular case, the ergodic con- N 1 i 2π (l 1)(j i) dition is fulfilled since, for sufficiently large time, the Tij (γ) = 1 e N − − , h i (γ) − random walker can reach any node of the network in- l=2 1 λl dependently of the initial node. Therefore, well-known X − h i results for the mean first passage time T , which gives (γ) ij and Tii(γ) = N, λl is given by Eq. (14). Therefore the average number of steps of a discrete-timeh i random h i walker to start in node i and reach for the first time j, N N i 2π (l 1)(j i) 1 γ 1 e N − − still apply for circulant directed structures [2]. Tij (γ) = µ − γ , (24) h i N m µ m=2 ! l=2 l In terms of the left and right eigenvectors ( φ¯l and φl , X X | i respectively) of the transition matrix W of a Markovian where we use the result µ1 = 0. In a similar way, we random walker and the associated eigenvalues λl (we de- obtain for the Kemeny’s constant in Eq. (23) note λ1 = 1), we have for i = j [6, 18, 27, 28] 6 N N N 1 γ 1 1 j φl φ¯l j i φl φ¯l j (γ)= µ . (25) T = h | i | − h | i | , (22) T N m µγ ij m=2 ! l h i 1 λl j φ φ¯ j l=2 l=2 1 1 X X X − h | i | 1 Our findings in Eqs. (24)-(25) are valid for directed and and the mean first return time Tii = ( i φ1 φ¯1 i )− . In addition, in structures with theh samei fractionalh | i | degree undirected connected circulant networks and allow to cal- we have the global time [6, 18, 27, 28] culate analytically the mean first passage time (MFPT) and the Kemeny’s constant through the specification of N 1 the coefficients c0,c1,...,cN 1 in Eq. (15). = . (23) − T 1 λl l In Fig. 5, we depict the MFPT for a directed ring and an X=2 − 9
different directed structures with topologies and prop- erties discussed in Figs. 2-3. In Fig. 6, the curve (a) describes a fractional random walker in a directed ring with cN 1 = 1. We can see that the best strategy to explore the− network is defined by γ = 1. This is a deter- ministic limit where the walker visits a new node at each step. With the introduction of L´evy flights through the fractional dynamics with 0 <γ< 1, the stochastic trans- port increases the time (γ) and in the limit γ = 0 we 2T have complete = (N 1) /N also obtained for a complete networkT [28]. In curve− (b) we have a structure defined with the elements cN 1 = cN 2 = 1, we observe a similar behavior with an optimal− value− in the local-limit γ 1. In directed structures with more lines like in curves→ (c)- (f) we see a different behavior where a particular value of γ⋆ < 1 produces a maximum in the Kemeny’s constant; however, with the reduction of γ in the limit γ 0 all FIG. 6. (Color online) Kemeny’s constant T (γ) as a func- → 4 the cases approach to the value . The results tion of γ for different directed networks with N = 10 Tcomplete nodes. The networks have a topology similar to those show particular cases where the combination of nonlocal in Fig. 2 and are defined with the particular set of co- displacements and the bias generated by the directions efficients equal to 1: (a) {cN−1}, (b) {cN−1, cN−2}, (c) of the lines produce a global effect that reduces the effi- {cN−1, c2}, (d) {cN−1, c1, cN−2}, (e) {cN−1, cN−2, c2}, (f) ciency to visit all the nodes of the network. {cN−1, c1, cN−2, c2, cN−3}. The results were obtained with Eq. (25) and the Laplacian eigenvalues in Eq. (15). The 2 dashed horizontal line represents Tcomplete = (N − 1) /N for IV. BIASED TRANSPORT ON RINGS a random walker in a complete (fully connected) network.
In the previous section, we studied random walks on
4 circulant directed networks for which we considered the undirected ring, both with N = 10 nodes. The results weights Ω = A. We are now interested in the effects illustrate a completely different behavior in the fractional of the fractional transport when the matrix Ω describes dynamics in directed and undirected rings. First of all, as some type of bias determined by weights in the links. we can see in Fig. 5(a), for the directed ring defined with We explore the transport on a ring with transition prob- an adjacency matrix cN 1 = 1, the results show that for − abilities different from the directed and undirected rings γ = 1 and i = j, Tij (γ = 1) = dij , where dij = j i for studied before. j i. In the6 directedh ring, thei value γ = 1 produces− a de- terministic≥ dynamics where at each step the walker visits Let us now consider a probability 0 p 1 and a ring ≤ ≤ a new adjacent node with a cover time N. On the other with N nodes which are connected only to their first hand, in the interval 0 <γ< 1, the temporal evolution neighbors. In addition, the coefficients Ωij are defined is stochastic with a biased L´evy like dynamics increasing by a circulant matrix with non-null elements c1 = p and the MFPT but maintaining these times below or equal cN 1 = 1 p. In this way, we have the probability p − − to the value N. The results for the biased transport that to hop in one direction, and 1 p to the opposite direc- − emerge in the fractional dynamics on the directed ring tion. Using Eq. (15) for the eigenvalues of the Laplacian agree with previous studies showing that L´evy flights do matrix, and c0 = c1 + cN 1 = 1, we have − not always optimize the search problem in the presence iϕl iϕl µl =1 pe (1 p)e− (26) of an external drift [48]. In Fig. 5(b), we present the − − − results for times Tij (γ) in a symmetric ring. In this h i 2π case, the fractional dynamics with L´evy flights reduce where ϕl N (l 1) and 1 l N. ≡ − ≤ ≤ (γ) the MFPT found in the local limit γ = 1 for which the In the general case, the eigenvalues λl of the frac- random walker at each step moves with probability 1/2 tional transition matrix W(γ) are complex values. By from a node to one of its two neighbors. applying Eq. (14), we obtain Finally, it is important to have a global time that char- (γ) 1 γ acterizes the capacity of the random walker to explore iϕl iϕl λl =1 (γ) 1 pe (1 p)e− (27) the network. In structures such as circulant networks, − k − − − the Kemeny’s constant (γ) defined in Eq. (25) gives with a fractional degree a global value quantifyingT the efficiency of the fractional random walker to reach all the nodes. In Fig. 6, we N (γ) 1 iϕ iϕ γ present the values of (γ) as a function of γ (0 <γ 1) k = 1 pe l (1 p)e− l . (28) T 4 ≤ N − − − in circulant networks with N = 10 nodes. We explore l=1 X 10
FIG. 8. (Color online) Kemeny’s constant T (γ) as a function of γ for biased fractional transport on a ring with N = 104 nodes. The times T (γ) are calculated with the analytical result in Eq. (29). We explore different values of the bias pa- rameter p in the interval 0.5 ≤ p ≤ 1 codified in the colorbar. FIG. 7. (Color online) Eigenvalues of the transition matrix The case p = 1/2 recovers a symmetric random walker with W(γ) for biased fractional transport on a ring with N = 100 L´evy flights whereas the limit p = 1 describes the transport nodes. For γ = 1, the transition probabilities to the two on a directed ring. These two limits are explored in Figs. 4-5. nearest neighbors are p to move in one direction and 1 − p In the inset we present the value γ⋆, as a function of p, that in the opposite direction. We explore the effect of γ and p in maximizes T (γ) and is obtained numerically with Eq. (31). the eigenvalues of W(γ) for p = 0.9, 0.8, 0.7, 0.6. Each color represents a different set of eigenvalues, obtained with a given γ codified in the colorbar. The eigenvalues are calculated with for 0 < p 1 and 0 < γ 1. In particular, in the limit Eq. (27). The limit p = 1 represents a directed ring with p =1/2 we≤ recover the result≤ for the symmetric ring eigenvalues analyzed in Fig. 3(a).
N γ (γ) 1 S(γ)= k . (30) T 1 cos ϕm In particular, for p = 1/2 and γ = 1, we recover the m=2 − (γ=1) X eigenvalues λl = cos(ϕl) for the local random walk (γ) 1 N γ in a symmetric ring. with the fractional degree k = (1 cos ϕl) . N l=2 − (γ) The global time S(γ) is analyzed in detail in Ref. [28] to In Fig. 7 we show the eigenvalues λl of the transition characterize theT fractional transport onP undirected rings. matrix W(γ) for the biased fractional transport in a ring with N = 100 nodes. The particular limit p = 1 recov- In Fig. 8 we present the results obtained with Eq. (29) ers the transport on the directed ring presented in Fig. for the Kemeny’s constant describing the fractional trans- 3(a). We see how the bias modeled with the parame- port with bias in rings for different values of p. We ob- ter 0.6 p 0.9 reduces the imaginary component of serve in the behaviour of the Kemeny’s constant that, the eigenvalues≤ ≤ when p 1/2. In the limit p = 1/2 all for p in the interval 0.53 p 0.96, (γ) presents a → ≤ ≤ T the eigenvalues are real. We obtain similar results for maximum for a particular value γ denoted as γ∗. Cal- d 0 p 0.5 since, in this interval, the random walker has culating (γ) = 0, we obtain that the value γ∗ ≤ ≤ dγ T γ=γ⋆ the same dynamics but with a change in all the directions maximizing the Kemeny’s constant satisfies of the walker. N γ⋆ N γ⋆ Once deduced the eigenvalues of the Laplacian matrix, zl zl log(zl) = log(zm) (31) we can quantify globally the capacity of the fractional zm zm l,m l,m random walker to explore the network. Since the frac- X=2 X=2 tional degrees are the same for all the nodes, we use the iϕ iϕ with zl 1 pe l (1 p)e− l . In the inset in Fig. relation for the Kemeny’s constant in Eq. (25). There- ≡ − − − fore 8, we show the values γ∗ as a function of p, calculated numerically with Eq. (31) . This result determines which
N N iϕ iϕ γ L´evy flight strategy is the least efficient to explore the 1 1 pe l (1 p)e− l (γ)= − − − (29) network, i.e. for which value of γ the Kemeny’s constant T N 1 peiϕm (1 p)e iϕm l m=2 − (γ) has a maximum in the interval 0 <γ< 1. X=2 X − − − T 11
FIG. 9. (Color online) Fractional transport on directed networks with N = 50 nodes generated randomly with a probability p = 0.04 that define the probability to create a connection between two nodes. In (a)-(c) we present three connected networks with the respective eigenvalues of the transition matrix W(γ) in the complex plane. For each γ codified in the colorbar in the interval 0 < γ ≤ 1, we obtain N = 50 points representing the eigenvalues of W(γ); the eigenvalues for γ = 1 are presented with black dots. In (d) we explore the numerical values of the global time hT(γ)i in Eq. (32) for the three networks. The inset shows the results for the Kemeny’s constant T (γ).
V. DIRECTED RANDOM NETWORKS connectivity of the network, in the cases explored, the (γ) eigenvalue λ1 = 1 is unique for all γ. The remaining eigenvalues are represented as black dots in the complex In this section, we explore the fractional transport on plane for γ = 1, and with different colors we show how directed networks generated stochastically with an algo- the reduction of γ concentrates the eigenvalues around rithm similar to the Erd˝os-R´enyi (ER) model [49]. For (γ) N nodes we have an adjacency matrix A and, in each the origin with λl = 1/(N 1) for l = 2, 3,...,N in the limit γ 0. − − non-diagonal entry Aij , we decide to include the value 1 → or 0 randomly with probabilities p and 1 p, respectively. In addition to the eigenvalues, it is useful to quantify the − However, the difference with the traditional ER model is capacity of the fractional random walker to explore the that the choice of Aij is independent of Aji and, as a whole network. In this case, we use the global MFPT result, A is the non-symmetric matrical representation T(γ) defined as of a directed network where pN(N 1) is the average h i − N N number of links in this structure. 1 T(γ) = Tij (γ) (32) In Fig. 9, we explore three random networks with size h i N 2 h i i=1 j=1 N = 50 generated with the value p = 0.04, following X X the same approach presented for the analysis of regu- that gives the average of the mean first passage time lar topologies. We use connected networks for all the Tij (γ) considering all the initial nodes i and all the nodes to analyze the modifications in the eigenvalues in- hnodes targeti j. In Fig. 9(d), we show the results ob- troduced by the fractional transport. The eigenvalues tained with the numerical values of the eigenvalues and of the transition probabilities and the respective net- left and right eigenvectors of the transition matrix W(γ) works are shown in Figs. 9(a)-(c). In this representation, for 0 < γ 1. In addition, in the inset in this fig- for each value γ, we calculate numerically the fractional ure, we include≤ the values of the Kemeny’s constant Laplacian matrix given by Eq. (5) and in this way, we 1 N (γ) − (γ) (γ) = 1 λ that only includes informa- have the elements that define the transition matrix W , T l=2 − l (γ) (γ) for which we calculate the eigenvalues λl . Due to the tion of theP eigenvalues of W . 12
The results in Fig. 9 illustrate the different cases that can VII. APPENDIX A: PROPERTIES OF THE appear in the fractional transport on directed structures FRACTIONAL LAPLACIAN for the same type of random network. In the informa- tion in Fig. 9(d), the global MFPT T(γ) shows that In this appendix, we demonstrate briefly that the es- h i in some cases, like in the network (a), the effect of γ sential good properties (i)-(iii) of the Laplacian in Eq. can improve the efficiency of the random walker to reach (3) are conserved by the fractional Laplacian matrix a target, in comparison with the normal random walk Lγ in the interval 0 < γ 1 in order to guarantee γ 1. In networks (b) and (c) we see two cases where that the fractional transition≤ matrix (8) is stochastic, i.e. → the long-range dynamics increases the values of T(γ) in (γ) N (γ) h i 0 wi j 1 with l=1 wi l = 1. By the same demon- comparison with the normal dynamics. In addition, the stration≤ → we≤ show that for the→ class of strongly connected eigenvalues in (a)-(c) reveal different spectral properties directed graphs (i.e.P between each pair of nodes i, j at with the change of γ. However, in contrast with the reg- least one incoming and one outgoing connecting path of ular cases explored with circulant matrices, in this case finite length exists i j [31]) described by a Laplacian the Kemeny’s constant is not a good descriptor of the matrix of the form (3)↔ the fractional walk with transi- global activity in the fractional transport, revealing that (γ) tion matrix (8) is aperiodic ergodic. To this end let us in the networks analyzed the eigenvectors of W con- (out) I introduce Λ > max(ki ) thus the matrix Λ L has tain important information for a global characterization (out) − of the dynamical process. uniquely non-negative entries (Λ ki )δij +Ωij 0. n − t(ΛI L≥) Then it follows that (ΛI L) and hence e − = tn n Λt −tL ∞ I n=0 n! (Λ L) = e e− (t > 0) also conserve this property (where− for a strongly connected structure there P n exists a n0 such that all entries [(ΛI L) ]ij > 0 (n n0) are strictly positive, for non-connected− structures the≥ in- VI. CONCLUSIONS n equality [(ΛI L) ]ij 0 remains true for all n. Then for a strongly connected− ≥ structure since eΛt > 0 the matrix tL exponential (e− )ij > 0. It follows then that the non- In this work, we presented a general approach to ex- L L diagonal elements (I e t ) = (e t ) < 0 (i = j) amining ergodic Markovian nonlocal random walks gen- − ij − ij are uniquely negative− (in a non-connected− structure6 non- erated through the fractional Laplacian of directed net- positive). Since the zero eigenvalue to the constant eigen- works. This formalism is explored for different types of vector i Ψ = 1 of L is conserved by the matrix func- strongly connected networks defined in terms of circulant 1 √N h | tiL tion (I e− ) Ψ = 0, it follows that matrices and also for random directed networks. We an- − | 1i alyzed the eigenvalues of transition matrices that define N the random walk and the effect of the fractional dynamics tL tL (I e− )ii = (I e− )ij > 0 (33) in the spectrum, showing how eigenvalues are modified in − − − j=i the complex plane. In addition, for regular networks, we X6 analyze the Kemeny’s constant, which gives a good de- tL in a strongly connected directed graph (and (I e− )ii scription of the global dynamics. With this global time, 0 if the network is not strongly connected). Applying− (7)≥ defined only in terms of the eigenvalues, we analyzed the on both sides of this relation yields Eq. (6) and con- effect of biased nonlocal transport, showing that the ex- (γ) serves the signs in Eq. (33), i.e. k = (Lγ ) > 0 ploration of the network can be more effective in some i ii together with (Lγ ) < 0 for i = j in strongly connected particular cases reducing Kemeny’s constant. We also ij directed graphs. If the graph is6 disconnected properties identified different configurations for which the capacity (i)-(iii) still hold; however the Laplacian matrix contains of the fractional random walker to reach any node is re- blocks of zero entries which are conserved by all inte- duced in comparison with a random walker with hops to ger powers of the Laplacian matrix and hence also by neighbor nodes. This is a fundamental difference with the L the matrix exponential e t leading by virtue of Eq. (7) results observed in undirected networks where the frac- − that these blocks of zero-entries are still conserved in the tional dynamics always improve the transport through fractional Laplacian matrix, and hence in the resulting long-range displacements. Finally, we explored the trans- fractional walk ergodicity is lost. Compare especially the port of random directed networks; in this case, the global Laplacian and fractional Laplacian, respectively, in Figs. activity is analyzed using the average of the mean first 1(a)-(b). In Fig. 1(a) the block of zero entries of the passage time to all the nodes considering all the initial non-strongly connected Laplacian matrix is conserved by conditions. This quantity depends on the eigenvalues and the fractional Laplacian. eigenvectors of the transition matrix that defines the pro- cess. Our findings and methods introduced are general In this way we have demonstrated that the good prop- and pave the way to further extensions for the exploration erties (i)-(iii) of L are indeed conserved by the frac- of fractional dynamical processes on directed structures tional Laplacian matrix Lγ in the interval of convergence with possible applications in the understanding of human 0 <γ < 1 of the integral representation in Eq. (7) mobility, data analysis, synchronization, among others. where for strongly connected directed graphs the frac- 13 tional Laplacian generates an aperiodic ergodic walk. In and combining Eqs. (35)-(39), we obtain the decomposi- this proof we did not make use of whether or not the tion of the elements of the circulant fractional Laplacian Laplacian matrix of a directed graph is diagonalizable. matrix (36), i.e. the fractional Laplacian matrix for the It includes therefore also the cases where the Laplacian finite ring in Eq. (18), for N finite but not necessarily matrix has Jordan canonical form. For a more detailed large analysis (however focused on undirected graphs) we refer N to Ref. [28]. γ 1 iϕℓ γ iϕℓ(q p) (L )pq = (1 e− ) e − (40) N − ℓ X=1 VIII. APPENDIX B: FRACTIONAL (γ) q p+Nt ∞ γ LAPLACIAN FOR DIRECTED RINGS = q p = ( 1) − (41) A − − q p + Nt t=0 X − In this Appendix, we analyze the elements of Lγ for a for q p =0,...,N 1. Hence we get directed ring of finite size N where N is not necessarily − − large. Let us consider the Laplacian L for a directed ring γ (γ) q p γ (L )pq = q p = ( 1) − defined with elements c0 = cN 1 = 1 and cm = 0 for A − − q p m =1, 2,N 2. Using Eq. (15)− we have the eigenvalues − ∞ − q p+Nt γ + ( 1) − , γ (0, 1] (42) N − q p + Nt ∈ 2π t=1 iϕℓm − µℓ = cm(1 e ), ϕℓ = (ℓ 1). (34) X − N − m=1 (circulant). For the finite directed ring all matrix ele- X ments of the fractional Laplacian are non-vanishing, i.e. Therefore, the Laplacian eigenvalue is given by this structure hence is strongly connected and hence ape-
iϕℓ(N 1) iϕℓ riodic ergodic. We see in this result for the finite ring that µℓ =1 e − =1 e− (35) − − the first term in the series (t = 0), namely
In addition, the elements of the fractional Laplacian ma- q p γ ( 1) − q p , q p 0 trix given by Eq. (18) for 0 < γ 1 are determined γ ( ) − − − ≥ ≤ (L ) ∞ = (43) by pq 0, q p< 0 N − γ γ γ is the matrix element of Eq. (19) for the directed infinite (L )pq = (L )q p = µℓ p Ψℓ Ψℓ q − h | ih | i γ ( ) γ ( ) ℓ=1 ring. The fractional degree (L ) ∞ = 1 and (L ) ∞ < 0 X qq pq N γ ( ) for q>p whereas (L )pq∞ = 0 is null for q
γ 1 γ 1 γ ∞ (Nt) ∞ τ N − − − − dτ − . (45) Γ( γ) ≈ Γ( γ) ∼−Γ(1 γ) t=1 N X − Z − −
[1] A. Barrat, M. Barth´elemy, and A. Vespignani, Dynami- [21] A. P. Riascos and J. L. Mateos, cal Processes on Complex Networks (Cambridge Univer- J. Stat. Mech. 2015, P07015 (2015). sity Press, Cambridge, 2008). [22] A. P. Riascos and J. L. Mateos, [2] N. Masuda, M. A. Porter, and R. Lambiotte, Phys. Rev. E 92, 052814 (2015). Phys. Rep. 716–717, 1 (2017). [23] T. M. Michelitsch, B. A. Collet, A. P. Riascos, [3] S. Brin and L. Page, A. F. Nowakowski, and F. C. G. A. Nicolleau, Comput. Netw. ISDN Syst. 30, 107 (1998). Chaos Solitons & Fractals 92, 43 (2016). [4] J. Leskovec, A. Rajaraman, and J. D. Ullman, Mining of [24] T. M. Michelitsch, B. A. Collet, A. P. Riascos, Massive Datasets, 2nd ed. (Cambridge University Press, A. F. Nowakowski, and F. C. G. A. Nicolleau, Cambridge, 2014). J. Phys. A: Math. Theor. 50, 055003 (2017). [5] L. Ermann, K. M. Frahm, and D. L. Shepelyansky, [25] S. de Nigris, T. Carletti, and R. Lambiotte, Rev. Mod. Phys. 87, 1261 (2015). Phys. Rev. E 95, 022113 (2017). [6] A. P. Riascos, J. Wang-Michelitsch, and T. M. Miche- [26] T. M. Michelitsch, B. A. Collet, A. P. Riascos, litsch, Phys. Rev. E 100, 022312 (2019). A. F. Nowakowski, and F. C. G. A. Nicolleau, [7] M. Lenormand, A. Bassolas, and J. J. Ramasco, J. Phys. A: Math. Theor. 50, 505004 (2017). J. Transp. Geogr. 51, 158 (2016). [27] A. P. Riascos, T. M. Michelitsch, B. A. Collet, [8] H. Barbosa, M. Barth´elemy, G. Ghoshal, C. R. James, A. F. Nowakowski, and F. C. G. A. Nicolleau, M. Lenormand, T. Louail, R. Menezes, J. J. Ramasco, J. Stat. Mech. 2018, 043404 (2018). F. Simini, and M. Tomasini, Phys. Rep. 734, 1 (2018). [28] T. M. Michelitsch, A. P. Riascos, B. A. Collet, A. F. [9] A. P. Riascos and J. L. Mateos, Nowakowski, and F. C. G. A. Nicolleau, Fractional Dy- PLOS ONE 12, e0184532 (2017). namics on Networks and Lattices (ISTE/Wiley, London, [10] D. Loaiza-Monsalve and A. P. Riascos, 2019). PLOS ONE 14, e0213106 (2019). [29] A. Allen-Perkins and R. F. S. Andrade, [11] A. P. Riascos and J. L. Mateos, J. Stat. Mech. 2019, 123302 (2019). Sci. Rep. 10, 4022 (2020). [30] E. Bautista, P. Abry, and P. Gonalves, [12] V. Belik, T. Geisel, and D. Brockmann, Appl. Netw. Sci. 4, 57 (2019). Phys. Rev. X 1, 011001 (2011). [31] M. Benzi, D. Bertaccini, F. Durastante, and I. Simunec, [13] L. D. Valdez, L. A. Braunstein, and S. Havlin, J. Compl. Net. 8, cnaa017 (2020). Phys. Rev. E 101, 032309 (2020). [32] R. Burioni and D. Cassi, [14] P. Blanchard and D. Volchenkov, Random Walks and Mod. Phys. Lett. B 11, 1095 (1997). Diffusions on Graphs and Databases: An Introduction, [33] K. Hatori, T. Hatori, and H. Watanabe, Springer Series in Synergetics (Springer, Berlin, 2011). Prog. Theor. Phys. Suppl. 92, 108143 (1999). [15] E. W. Montroll and G. H. Weiss, [34] R. Burioni, D. Cassi, and A. Vezzani, J. Math. Phys. 6, 167 (1965). Phys. Rev. E 60, 1500 (1999). [16] B. D. Hughes, Random Walks and Random Environ- [35] G. Bianconi and S. N. Dorogovstev, ments: Vol. 1: Random Walks (Oxford University Press, J. Stat. Mech. 2020, 014005 (2020). Oxford, 1996). [36] A. P. Riascos and J. L. Mateos, “Random [17] J. D. Noh and H. Rieger, walks on weighted networks: Exploring local Phys. Rev. Lett. 92, 118701 (2004). and non-local navigation strategies,” (2019), [18] A. P. Riascos and J. L. Mateos, arXiv:1901.05609 [cond-mat.stat-mech]. Phys. Rev. E 86, 056110 (2012). [37] R. Lambiotte, R. Sinatra, J.-C. Delvenne, [19] E. Estrada, J.-C. Delvenne, N. Hatano, J. L. Ma- T. S. Evans, M. Barahona, and V. Latora, teos, R. Metzler, A. P. Riascos, and M. T. Schaub, Phys. Rev. E 84, 017102 (2011). J. Compl. Net. 6, 382 (2018). [38] M. E. J. Newman, Networks: An Introduction (Oxford [20] A. P. Riascos and J. L. Mateos, University Press, Oxford, 2010). Phys. Rev. E 90, 032809 (2014). [39] B. Mohar, in Graph Theory, Combinatorics, and Appli- cations, Vol. 2, edited by Y. Alavi, G. Chartrand, O. R. 15
Oellermann, and W. A. J. Schwenk (New York, Wiley, [43] T. Weng, M. Small, J. Zhang, and P. Hui, 1991) pp. 871–898. Sci. Rep. 5, 17309 (2015). [40] B. Mohar, in Graph Symmetry: Algebraic Methods and [44] Q. Guo, E. Cozzo, Z. Zheng, and Y. Moreno, Applications, NATO ASI Series (Series C: Mathematical Sci. Rep. 6, 37641 (2016). and Physical Sciences), Vol. 497, edited by G. Hahn and [45] A. Allen-Perkins, A. B. Serrano, T. A. de Assis, and G. Sabidussi (Springer, Dordrecht, 1997) pp. 227–275. R. F. S. Andrade, Phys. Rev. E 100, 030101(R) (2019). [41] R. Aldrovandi, Special Matrices of Mathematical Physics [46] K. B. Gustafson, B. S. Bayati, and P. A. Eckhoff, (Stochastic, Circulantand Bell matrices) (World Scien- Front. Ecol. Evol. 5, 35 (2017). tific, Singapore, 2001). [47] A. Pachon, F. Polito, and C. Ricciuti, [42] P. Van Mieghem, Graph Spectra for Complex Networks Discrete Contin. Dyn. Syst. Ser. B (In press) (2020), 10.3934/dcdsb.2020170. (Cambridge University Press, New York, 2011). [48] V. V. Palyulin, A. V. Chechkin, and R. Metzler, Proc. Natl. Acad. Sci. U.S.A 111, 2931 (2014). [49] P. Erd¨os and A. R´enyi, Publ. Math. (Debrecen) 6, 290 (1959).