Fault Tolerance of Random Graphs with respect to Connectivity ...

1 downloads 0 Views 336KB Size Report
Dec 21, 2017 - (2014). [54] B. Karrer, M. E. J. Newman, and L. Zdeborová, Phys. Rev. Lett. ... [67] G. Paul, T. Tanizawa, S. Havlin, and H. E. Stanley, Eur. Phys.
Fault Tolerance of Random Graphs with respect to Connectivity: Phase Transition in Logarithmic Average Degree Satoshi Takabe,∗ Takafumi Nakano, and Tadashi Wadayama

arXiv:1712.07807v1 [cond-mat.dis-nn] 21 Dec 2017

Department of Computer Science, Nagoya Institute of Technology, Gokiso-cho, Showa-ku, Nagoya, Aichi, 466-8555, Japan (Dated: December 22, 2017) The fault tolerance of random graphs with unbounded degrees with respect to connectivity is investigated. It is related to the reliability of wireless sensor networks with unreliable relay nodes. The model evaluates the network breakdown probability that a graph is disconnected after stochastic node removal. To establish a mean-field approximation for the model, the cavity method for finite systems is proposed. Then the asymptotic analysis is applied. As a result, the former enables us to obtain an approximation formula for any number of nodes and an arbitrary and degree distribution. In addition, the latter reveals that the phase transition occurs on random graphs with logarithmic average degrees. Those results, which are supported by numerical simulations, coincide with the mathematical results, indicating successful predictions by mean-field approximation for unbounded but not dense random graphs.

Introduction. Connectivity is a simple but fundamental property of graphs, as apparent from the fact that the first work on Erd¨ os–R´enyi (ER) random graphs examined the connectivity probability and its critical phenomenon [1]. Connectivity also has been investigated as a measure of the reliability of wireless communication networks first introduced by Moore and Shannon [2]. In the model, one considers the connectivity of networks after communication links, edges in a graph, are removed stochastically. Combined with the notion of random graphs [3], studies of reliability defined by connectivity [4] have provided theoretical insight into various technologies such as software fault tolerance [5] and the random key predistribution scheme [6]. After the emergence of network science dealing with complex networks [7, 8], a fraction of a giant component (GC) has also been studied as a measure of network reliability against accidental fault of (or vulnerable attacks to) nodes or edges [9–11]. This model has been applied to real networks such as power grids [12, 13]. Theoretically the mean-field approximation in statistical mechanics contributes to those studies, revealing that network topology strongly affects fault-tolerance properties, e.g. a percolation threshold. From the viewpoint of statistical mechanics, some mean-field approximations for random graphs depend on the order of average degrees. For dense random graphs whose average degree is proportional to the number of nodes n, the replica method, which is known as a meanfield approximation, is often used, contributing to successful establishment of the physical pictures of the spinglass theory [14]. Additionally, it has been applied to various fields such as neural networks [14–16], wireless communications [17], and inference problems including compressed sensing [18, 19]. Numerous interesting systems are mapped to spin systems on sparse graphs with bounded, i.e., O(1), degree: error-correcting codes [20, 21], optimization prob-

lems [22, 23], community detection [24, 25], and so on. As mean-field approximations for sparse random graphs, the cavity method [26–29], and the generating function method [30, 31] have been used along with the Monasson’s replica method [32]. These methods provide good evaluations of properties and physical pictures of those systems (e.g., [33–35]). Moreover, they offer conjectures which are sometimes generalized to mathematically rigorous results. For example, percolation thresholds in the node fault model were predicted using mean-field analysis [11], as proved generally later [36]. Considering the intermediate regime, where the average degree is scaled as O(ln n), the mean-field analysis is a challenging task. Because the average degree is unbounded but not dense, the mean-field approximations for sparse or dense random graphs are naively inapplicable. In addition, the absence of the analysis might result from the fact that almost all known models exhibit phase transitions related to either sparse or dense graphs. Recently, however, probabilistic analyses have shown that some phase transitions in compressed sensing [37–39] and community detection [40, 41] occur in the intermediate regime. Providing an example of mean-field analyses for disordered systems in the regime is expected to be useful not only to extend the scope of applications of the approximation but also to predict mathematical statements. The fault tolerance of random graphs with respect to connectivity has attracted a great deal of interest along with development of IoT technologies [42]. In the model, disconnectivity probability of random graphs after stochastic node removals is evaluated, which describes reliability of the wireless sensor network with unreliable relay nodes [43, 44]. Nodes of this kind simplify, e.g., the energy harvesting relay node [45, 46], which accidentally fails activation. As for critical phenomena, the threshold is not estimated, although it is expected to lie in the O(ln n) regime [47], which suggests that the model is a

2 good example to examine the phase transition in the intermediate regime. The motivation of this Letter is to establish a meanfield analysis for the node fault model of random graphs with logarithmic average degrees and investigate its phase transition. In the analysis, we first obtain a selfconsistent cavity equation for disordered finite systems. It enables us to obtain the approximation formula with finite n, which coincides asymptotically with the numerical results. The asymptotic analysis then reveals that the model exhibits a discontinuous phase transition with a threshold lying in the O(ln n) regime. Finally, we discuss the validity of the mean-field analysis and some future avenues of research. Model Definition Let G = (V, E) be an undirected simple graph, i.e., a graph without self-loops and multiple edges, where V = {1, . . . , N } and E ⊂ V 2 respectively represent a vertex set and an edge set. Connectivity is defined as a graph property by which there exists a path between any pair of nodes in the graph. We consider the node fault model [47], also known as the random node breakdown [11], in which each node is removed independently with constant probability , called the node breakdown probability. The resultant graph is named a survival graph. The network breakdown probability P () (G) is defined as the probability that a survival graph is disconnected. Our goal is to evaluate the network breakdown probability averaged over random graphs with unbounded degrees and to investigate its phase transition. A random graph ensemble Ωn is defined as a probn abilistic space (G n , 2G , Q) where G n is a set of simple n graphs with n nodes, 2G represents is the power set of n G n , and Q is a probability measure of 2G . As a representative model, we restrict ourselves to the case in which Ωn is characterized solely by the number of nodes n and degree distribution pn (k) (0 ≤ k ≤ n − 1) [30], which indicates that graphs sampled from Ωn have no correlation such as degree–degree correlation [50], on average. The network breakdown probability PΩn () is then defined by the average of probability P () (G) over a random graph ensemble Ωn . Similar to the network breakdown probability, let us denote the average fraction of GC by ρΩn . It is noteworthy that the phase transition of ρΩn in the node fault model is known as network resilience [11]. Mean-field Analysis Because connectivity is realized when all nodes in survival graphs belong to the GC, we first evaluate the fraction of the GC ρΩn . The mean-field analysis of ρΩn has been investigated in [51] using the generating function method, in [52–54] using the cavity method, and in [55] using nonbacktracking expansions. Herein, we briefly describe mean-field analysis based on the cavity method. We introduce a classical spin model on a given graph G = (V, E) as in [53]. We define a binary variable si for the i-th node by si = 1 if the node is active and si = 0 otherwise. The other binary variable σi is introduced to

describe whether the i-th node is contained to the GC: σi = 1 if the node does not belong to the GC, and σi = 0 otherwise. After some nodes have broken down, i.e., the sequence {si } is fixed, a constraint for the i-th node to belong to the GC is given as Y σi = (1 − sj + sj σj ), (1) j∈∂i

where ∂i represents a set of neighbors of the i-th node. To find the giant component satisfying Eq. (1), we use the sum-product algorithm. Given that G is a tree, σi can be obtained exactly as Y σi = (1 − sj + sj mj→i ), (2) j∈∂i

where mj→i is a message from the j-th node to the i-th node, which takes a value of one if σj = 1 holds on the cavity graph G\i, and zero otherwise. These messages satisfy cavity equations for the GC [53], expressed as Y mi→j = (1 − sk + sk mk→i ). (3) k∈∂i\j

Although the sum-product algorithm is no longer exact for graphs with cycles, it is regarded as a kind of Bethe– Peierls approximation [56, 57]. To calculate the random graph average, the replica symmetric (RS) cavity method is applied to the system. Although earlier works directly take the large-n limit, as with the case of the conventional RS cavity method, the method for random graphs with finite n is proposed here to consider random graphs with unbounded average degrees. Considering that the messages {mi→j } are updated re(t) cursively by the cavity equations (3), let In be the probability that the message mi→j takes one for randomly chosen graph G and its edge (i, j) from Ωn at the t-th iteration step. If the correlation between any pair of spins (t) is negligible, then In satisfies the following relations. In(t+1) =

n X kpn (k)

hki

k=1

( + (1 − )In(t) )k−1 ,

(4)

Pn where hki ≡ k=1 kpn (k) represents the average degree of random graphs. A fixed-point solution In is obtained in the large-t limit, which is given as In =

n X kpn (k) k=1

hki

( + (1 − )In )k−1 .

(5)

The average fraction of the GC under the mean-field approximation is then obtained as ρMF Ωn =

n X k=1

pn (k)( + (1 − )In )k .

(6)

3

k=0

1 Network breakdown prob. Pλ,n (ϵ)

The breakdown probability PΩn () is then approximated using the evaluation presented above. As a graph is connected when each active node belongs to the GC of the survival graph, the approximation formula of the network breakdown probability is given as n PΩMF () = 1 −  + ρMF Ωn n !n n X k = 1 − 1 − (1 − ) pn (k)[ + (1 − )In ] .

0.8 Analysis: n = 28 n = 29

0.6

n = 211 n = 213 0.4

n = 220 Simulation: n = 28 n = 29

0.2

n = 211

(7)

n = 213 0

Equation (4) has a trivial solution In = 1 correspondent to a “broken phase” in which ρMF Ωn = 0 holds. It ˜ also possibly has another solution In in [0, 1). Because  is assumed to be constant, PΩMF () converges to one if I˜n n remains positive in the large-n limit. We therefore examine the case in which the average degree and resultant I˜n depend on n using the RS cavity method for finite-size systems. Asymptotic Analysis. To investigate the phase transition of connectivity, here we execute an asymptotic analysis of the approximation formula. Unfortunately, it is difficult to obtain an explicit form of the critical threshold in general. Alternatively, we describe the analytical results for some well-known random graph ensembles. (i) ER random graphs. The ER random graph with n nodes is characterized by the binomial degree distribu n−k−1 k tion, pn (k) = n−1 p , where p ≡ c/(n − 1) k (1 − p) and c is the average degree. Then Eq. (5) reads as n−2

In = [1 − p(1 − )(1 − In )]

.

(8)

MF () The approximate network breakdown probability Pn,c is given as h in n−1 MF Pn,c () = 1 − 1 − (1 − ) {1 − p(1 − )(1 − In )} .

(9) If we set c = d ln n, then the asymptotic analysis of Eq. (8) shows that the leading term of I˜n is given as n−d(1−) . This asymptotic form is then substituted to In in Eq. (9). Using the fact that, in the large-n limit, (1 − an−b )n converges to zero when b < 1, to e−a when b = 1 and to one otherwise, the threshold ∗ of node breaking probability is obtained by ∗ = 1−1/d. In other MF () jumps from zero to one at the threshold words, Pn,c ∗  . This result suggests that ER random graphs loses connectivity even when  = 0 if c < ln n holds, which is proved in [1]. (ii) Regular random graphs. The regular random (RR) graph is characterized by pn (k) = δk,λ where δn,m represents Kronecker’s delta. Although RR graphs themselves have connectivity with high probability if λ ≥ 3 [58], results show that connectivity is lost by node breakdowns for any  > 0. As in the case of ER random graphs, we scale the degree as λ = d ln n. The asymptotic analyses

0.2

0.3

0.4 0.5 0.6 Node breakdown prob. ϵ

0.7

0.8

FIG. 1. Network breakdown probability Pλ,n () as a function of node breakdown probability  for regular random graphs with n nodes of which the degree is λ = log2 n. Solid lines represent evaluations by the mean-field approximation. Symbols show numerical results averaged over 106 randomly generated survival graphs. The vertical dashed line represents the phase transition threshold ∗ = 1/2 by asymptotic analysis.

yield that the leading term of I˜n is given as n−d ln(1/) . The threshold is given as ∗ = e−1/d . Figure 1 presents numerical results and mean-field evaluations of the network breakdown probability for RR graphs with degree λ = log2 n. The numerical simulation is performed by generating 103 random graphs and by checking the connectivity of 103 survival graphs per random graph using depth-first search. Results show that the mean-field evaluations by Eqs. (5) and (7) are asymptotically close to the numerical results, indicating that the analytical result ∗ = 1/2 is correct. (iii) Scale-free networks. We also consider scale-free (SF) networks, with the degree distribution given as a power-law distribution pn (k) ∝ k −γ . To parameterize the ensemble with γ and average degree hki, we consider ˜ The random SF networks with the minimum degree k. degree distribution is therefore written as ˜ − ζ(γ, n), pn (k) = Cn−1 k −γ (k˜ ≤ k < n), Cn = ζ(γ, k) (10) P∞ −s where ζ(s, a) = (k + a) is the Hurwitz zeta k=0 function. Then, the average degree is given as hki = ˜ − ζ(γ − 1, n)], which is well-defined when Cn−1 [ζ(γ − 1, k) γ > 2. As the above cases illustrate, it is apparent that the possible threshold is trivial, i.e., ∗ = 0 when hki = O(1). If the average degree is again scaled as d ln n, then the minimum degree also depends on n because hki = (γ − 1)/(γ − 2)k˜ + O(1) holds. The selfconsistent equation of In , Eq. (5), is expressed as   1 ˜ [ + (1 − )I]k−1 Φ  + (1 − )In , γ − 1, k˜ , Cn hki (11) P∞ where Φ(z, s, a) = k=0 z k /(k + a)s is called the Lerch In =

4 1 RR ER 0.8

SF (γ = 6) SF (γ = 3) SF (γ = 2.1)

ϵ∗

0.6

0.4

0.2

0

0

1

2

d=

⟨k⟩ ln n

3

4

5

FIG. 2. Phase transition threshold ∗ as a function of a factor of the average degree hki given as d = hki/ ln n.

transcendent. In this case, the leading term of I˜n is found ˜ to be n−d ln(1/) by an asymptotic form of Φ(z, s, a) [59]. Finally, the threshold is given as ∗ = e−(γ−1)/[(γ−2)d] . The asymptotic analysis enables us to compare the reliability of each random graph ensemble. Figure 2 presents a summary of the thresholds ∗ of random graphs of three types as a function of d ≡ hki/ ln n. Although the threshold for ER random graphs exists only when d ≥ 1, other ensembles have the threshold for any d. With factor d fixed, the threshold of SF networks increases monotonously as the power γ increases. Particularly, ∗ for RR graphs coincides with that for SF network in the large-γ limit, which is always larger than that for ER random graphs. This fact is intuitive because the minimum degree, which is maximized by eliminating the variance of degree distribution, influences connectivity. We can thereby conclude that RR graphs is the most reliable ensemble among them with respect to connectivity in the node fault model when average degrees are fixed. Summary and Discussion. In this Letter, we presented consideration of the node fault model defined by graph connectivity. Our analysis comprises two stages: a mean-field analysis for finite systems based on the cavity method and its asymptotic analysis to obtain the phase transition threshold. The mean-field analysis yields the approximation formula available for any number of nodes n and an arbitrary degree distribution. It is also revealed that the model exhibits a discontinuous phase transition in the intermediate regime where the average degree is given as O(ln n). Below the threshold ∗ , the network breakdown probability vanishes in the large-system limit, although it converges to one when  > ∗ . The accuracy of our approximation formula and thresholds are confirmed by numerical simulations. The mathematical rigor of our analytical results is worth noting. For ER random graphs, in fact, we can prove the correctness of our prediction [60] by extending

the result by Erd¨os and R´enyi [1]. In contrast, because the proof is specific to the ensemble, the rigor of our predictions on other ensembles is left as an open problem. We emphasize, however, that the correctness of the specific case is strong evidence of the validity of our analyses. Next we can discuss the validity of the approximation itself. In the view of spin-glass theory, the cavity method is recognized as the Bethe–Peierls approximation. In the case of sparse random graphs, the use of approximation is justified by the locally tree-like structure of graphs that the length of cycles are typically scaled logarithmically [35]. In this manner, our analysis is justified because the length of cycles is evaluated as O(ln n/ ln ln n) [7]. From a mathematical perspective, however, the justification remains an open problem because the known rigorous results [61] only treat random graphs with bounded degrees. In other words, our results provide not only conjectures related to analytical results, but also a challenging open problem in the mean-field theory. The node fault model presented in this Letter has some possible extensions. The first direction is consideration of intentional attacks to nodes and correlations in graphs. These are analogous to the original node fault model with respect to percolation [62, 63]. Studies of these extensions reveal the network topology dependence on reliability. It is another direction to examine the node fault model on random graphs in metric space. The representative model is known as random geometric graphs [64], which is often used as an abstract model of wireless ad hoc networks [65, 66]. These extensions enable us to optimize networks while conserving reliability, as in the case of percolation [67]. Our results suggest that the mean-field approximation in this Letter will be applicable to other systems with critical phenomena, possibly in the intermediate regime. We hope that this work will be helpful in elucidating interesting critical phenomena related to complex networks with logarithmic average degrees and resultant extensions of statistical mechanics for disordered systems. The authors warmly thank T. Nozaki, T. Hasegawa, K. Hukushima, and Y. Kabashima for fruitful discussions. ST also thank T. Kawamoto for carefully reading the manuscript. This research is partially supported by JSPS KAKENHI Grant Nos. 16H02878, 16K14267 (TW), and 17H06758 (ST).



s [email protected] [1] P. Erd¨ os, and A. R´enyi, Publ. Math. Debrecen 6, 290 (1959). [2] E. F. Moore and C. E. Shannon, J. Franklin Inst. 262, 191 (1956). [3] B. Bollob´ as, Random Graphs (Academic Press, London, 1985).

5 [4] M. O. Ball, C. J. Colbourn, and J. S. Provan, in Handbook of Operations Research and Management Science, 7 673 (1995). [5] A. Avizienis, IEEE Trans. on Soft. Eng. SE-11, 1491 (1985). [6] H. Chan, A. Perrig, and D. Song, in Proceedings of the 2003 IEEE Symposium on Security and Privacy, 2003 (Washington, DC, 2003) p.197. [7] D. J. Watts and S. H. Strogatz, Nature 393, 440 (1998). [8] A. L. Barab´ asi and R. Albert, Science 286, 509 (1999). [9] R. Albert, H. Jeong, and A. L. Barab´ asi, Nature 406, 378 (2000). [10] D. S. Callaway, M. E. J. Newman, S. H. Strogatz, and D. J. Watts, Phys. Rev. Lett. 85, 5468 (2000). [11] R. Cohen, K. Erez, D. Ben-Avraham, and S. Havlin, Phys. Rev. Lett. 85, 4626 (2000). [12] R. Albert, I. Albert, and G. L. Nakarado, Phys. Rev. E 69, 1 (2004). [13] M. Rosas-Casals, S. Valverde, and R. V. Sol´e, Int. J. Bifurcation Chaos 17, 2465 (2007). ´ Virasoro Spin Glass [14] M. M´ezard, G. Parisi, and M. A. Theory and Beyond (World Scientific, Singapore, 1987). [15] J. J. Hopfield, Proc. Natl. Acad. Sci. USA 79 2554 (1982). [16] E. Gardner, J. Phys. A 21, 257 (1988). [17] T. Tanaka, IEEE Trans. Inf. Theory 48 2888 (2002). [18] D. L. Donoho, A. Maleki, and A. Montanari, Proc. Natl. Acad. Sci. USA 106, 18914 (2009). [19] Y. Kabashima, T. Wadayama, and T. Tanaka, J. Stat. Mech. 2009, L09003 (2009). [20] Y. Kabashima and D. Saad, Eur. Phys. Lett. 45, 97 (1999). [21] T. J. Richardson and R. L. Urbanke, IEEE Trans. Inf. Theory 47, 599 (2001). [22] R. Monasson, R. Zecchina, S. Kirkpatrick, B. Selman, and L. Troyanskyk, Nature 400, 133 (1999). [23] M. Weigt and A. K. Hartmann, Phys. Rev. Lett. 84, 6118 (2000). [24] M. E. J. Newman, Phys. Rev. E 74, 036104 (2006). [25] A. Decelle, F. Krzakala, C. Moore, and L. Zdeborov´ a, Phys. Rev. Lett. 107, 065701 (2011). [26] S. Katsura, S. Fujiki, and S. Inawashiro, J. Phys. C 12, 2839 (1979). [27] M. M´ezard, G. Parisi, N. Sourlas, G. Toulouse, and M. Virasoro, in Spin Glass Theory and Beyond (World Scientific, Singapore, 1987), p. 199. [28] J. Pearl, Probabilistic Reasoning in Intelligent Systems: Networks of Plausible Inference (Morgan Kaufmann Publishers, San Francisco, CA, 1991). [29] M. M´ezard and G. Parisi, Eur. Phys. J. B 20, 217 (2001). [30] M. E. J. Newman, S. H. Strogatz, and D. J. Watts, Phys. Rev. E 64, 26118 (2001). [31] M. E. J. Newman, Phys. Rev. Lett. Rev. 89, 278071 (2002). [32] R. Monasson, J. Phys. A. 31, 513 (1999). [33] R. Albert and A. L. Barab´ asi, Rev. Mod. Phys. 74, 47 (2002). [34] S. N. Dorogovtsev, A. V. Goltsev, and J. F. F. Mendes, Rev. Mod. Phys. 80, 1275 (2008).

[35] M. M´ezard and A. Montanari, Information, Physics, and Computation (Oxford University Press, Oxford, 2009). [36] N. Fountoulakis, Internet Math. 4, 329 (2007). [37] E. L. Cand`es, Comp. Rendus Math. 346, 589 (2008). [38] R. Baraniuk, M. Davenport, R. DeVore, and M. Wakin, Constructive Approximation 28, 253 (2008). [39] D. Donoho and J. Tanner, J. Amer. Math. Soc. 22, 1 (2009). [40] E. Abbe, A. S. Bandeira, and G. Hall, IEEE Trans. Inf. Theory 62, 471 (2016). [41] E. Abbe, arXiv:1703.10146 (2017). [42] J. Yick, M. Biswanath, and G. Dipak, Comp. Networks 52, 2292 (2008). [43] O. R. Theologou and J. G. Carlier, IEEE Trans. Reliability 40, 210 (1991). [44] S. Y. Kuo, F. M. Yeh, and H. Y. Lin, IEEE Trans. Reliability 56, 288 (2007). [45] S. Sudevalayam and K. Purushottam, IEEE Comm. Surveys and Tutorials 13 443 (2011). [46] H. Kawabata, K. Ishibashi, S. Vuppala, and G. T. F. de Abreu, IEEE Internet of Things J. 4, 384 (2017). [47] T. Nozaki, T. Nakano, and T. Wadayama, in Proceedings of 2017 IEEE International Symposium on Information Theory (ISIT), 2017 (Aachen, Germany, 2017) p. 481. [48] R. Albert and A. L. Barabasi, Rev. Mod. Phys. 74, 47 (2002). [49] P. Holme, B. J. Kim, C. Yoon, and S. K. Han, Phys. Rev. E 65, 56109 (2002). [50] M. E. J. Newman, Phys. Rev. Lett. 89, 208701 (2002). [51] S. V. Buldyrev, N. W. Shere, and G. A. Cwilich, Phys. Rev. E 83, 016112 (2011). [52] Y. Shiraki and Y. Kabashima, Phys. Rev. E 82, 036101 (2010). [53] S. Watanabe and Y. Kabashima, Phys. Rev. E 89, 12808 (2014). [54] B. Karrer, M. E. J. Newman, and L. Zdeborov´ a, Phys. Rev. Lett. 113, 208702 (2014). [55] G. Tim´ ar, R. A. da Costa, S. N. Dorogovtsev, and J. F. F. Mendes, Phys. Rev. E, 95, 042322 (2017). [56] H. A. Bethe, Proc. Roy. Soc. London A 150, 552 (1935). [57] Y. Kabashima and D. Saad, Eur. Phys. Lett. 44, 668 (1998). [58] N. C. Wormald, J. Comb. Th., Series B 31, 156 (1981). [59] C. Ferreira and J. L. Lop´ez, J. Math. Anal. and Appl. 298, 210 (2004). [60] S. Takabe, T. Wadayama (in preparation). [61] A. Dembo and A. Montanari. Ann. Appl. Prob. 20, 562 (2010). [62] R. Cohen, K. Erez, D. ben-Avraham, and S. Havlin, Phys. Rev. Lett. 86, 3682 (2001). [63] A. V´ azquez and Y. Moreno, Phys. Rev. E 67, 015101 (2003). [64] E. N. Gilbert, J. Soc. Indust. Appl. Math. 9, 533 (1961). [65] P. Santi, ACM Comput. Surv. 37, 164 (2005). [66] C. P. Dettmann and O. Georgiou, Phys. Rev. E 93, 1 (2016). [67] G. Paul, T. Tanizawa, S. Havlin, and H. E. Stanley, Eur. Phys. J. B 38, 187 (2004).