We study a spatial network model with exponentially distributed link lengths on an underlying grid of points, undergoing a structural crossover from a random, Erdős-Rényi graph, to a d-dimensional lattice at the characteristic interaction range ζ. We find that, whilst far from the percolation threshold the random part of the giant component scales linearly with ζ, close to criticality it extends in space until the universal length scale ζ6/(6-d), for d<6, before crossing over to the spatial one. We demonstrate the universal behavior of the spatiotemporal scales characterizing this critical stretching phenomenon of mean-field regimes in percolation and in dynamical processes on d=2 networks, and we discuss its general implications to real-world phenomena, such as neural activation, traffic flows or epidemic spreading.
Israel Science Foundation10.13039/501100003977Office of Naval Research10.13039/100000006Defense Threat Reduction Agency10.13039/100000774HDTRA-1-10-1-0014United States-Israel Binational Science Foundation10.13039/501100001742National Science Foundation10.13039/1000000012015781Army Research Office10.13039/100000183Ministry of Science, Technology and Space10.13039/501100001738Japan Science FoundationMinistero degli Affari Esteri e della Cooperazione Internazionale10.13039/501100006601Bar-Ilan University10.13039/501100002744
Two decades after its inception, network theory has grown to become an exhaustive framework for taming the complexity of a multitude of systems at different scales [1–3]. Despite the many advances made in the study of critical phenomena [4] on random networks, their analysis in the presence of a geometric embedding has remained an important theoretical challenge, so far grounded on very little analytical results [5–9]. This represents today a fundamental barrier in the understanding of complex systems of great interest, like networks of neurons [10,11], communication [12] and transportation systems [13], or power-grids [14], whose functioning strongly intertwines their spatial constraints.
This Letter aims to explore the effects that a tunable characteristic link length ζ has on the universality of critical phenomena in a realistic model of a spatial network. We consider a variant on a lattice of the so-called Waxman model [15], where N=Ld nodes are distributed on an hypercube of linear size L with spacing a≡1, and edges are drawn with probabilities 𝒫ij∝κe-dij/ζ, for every pair of nodes i, j placed at the Euclidean distance dij. Like in Erdős-Rényi (ER) networks, the parameter κ≡⟨k⟩/(N-1) controls the density of links, with ⟨k⟩ the average connectivity, while ζ∈(0,L] is a characteristic link length. Although simplistic, the exponential wiring cost function of this model correctly fits the spatial distribution of link lengths observed in communication systems [16], like the Internet [17,18] or transport networks [19,20], interactions based on similarity in social networks [21], and recent experimental results [22–26] support the existence of scale-invariant organizational principles in the connectome of mammals based on it.
This disordered lattice is a natural benchmark for investigating crossover phenomena [27,28] in spatial networks. Depending on the values of ζ, in fact, our model crosses over between two limiting regimes: for ζ∼O(L), all link lengths are equally probable (i.e., 𝒫ij∼κ) and we recover a random, ER-like structure, whilst for ζ∼O(1/ρ) (where ρ is the average density of nodes), long links are prohibited, leading instead to a latticelike topology. By performing site percolation [29–31], we disclose an intriguing scenario: while far from criticality the random part of the giant component crosses over to a spatial regime at the characteristic scale ζ, close to the percolation threshold it surprisingly extends in space until the intrinsic scale ξ*(ζ)=ζ6/(6-d), for d<6. We call this structural expansion of mean-field regimes the critical stretching, and we explore its universal features by means of scaling arguments, and simulations in d=2. In particular, because the structure affects the collective phenomena occurring on it, the critical stretching should be reflected in processes evolving on this spatial network. We verify this claim on an SIR (susceptible-infective-removed) process [32], where the mean-field (diffusive) regime of the outbreak front’s velocity undergoes a temporal critical stretching, making the epidemic invade the system at “slow” rates for significantly longer times.
Crossover and critical stretching.—To investigate the structural crossover, we perform site percolation by randomly removing a fraction 1-p of the network’s nodes, for given average degree ⟨k⟩=4 and varying [33] characteristic link-lengths. For increasing ζ, a rapid monotonic decrease of the percolation threshold pc from the lattice value to the ER one pcMF=1/⟨k⟩ is observed [Fig. 1(b), and Ref. [20] ]. To analyze the response the scaling of percolation observables have to this variation, we consider the correlation length ξ2(p)≔∫r2g(r,p)dr/∫g(r,p)dr, where g(r,p) is the pair connectivity function for a given occupation probability p. By definition, ξ measures the mean distance between sites belonging to the same cluster, and it diverges as p approaches pc like ξ∼|p-pc|-νd. Numerical results [Fig. 1(a)] for networks in d=2, display excellent agreement with the critical exponents expected in the two limiting cases, i.e., ν2D=4/3 for ζ∼O(1) and νMF=1/2 for ζ∼O(L), further showing that a crossover from the former to the latter regime occurs closer and closer to pc for increasing ζ.
This critical crossover can be understood by invoking the Levanyuk-Ginzburg (LG) criterion [34]. Since the spectrum of our spatial network satisfies the classical LG assumptions [6], we expect fluctuations to be negligible as long as u3-d/2ζd≫1 for d<6, where u≡|p/pc-1| is the displacement of p from its critical value. To obtain this condition, we noticed that the scale up to which the mean-field regime extends in the giant cluster is ξMF∝ζξℓ1/2, where ξℓ(u)∼u-1 is the (dimensionless) correlation length in chemical space, and the factor ζ takes into account the characteristic length of links [35]. The critical crossover should then take place at u*(ζ)≡ζ-2d/(6-d) and, in light of the above, at the intrinsic scale ξ*(ζ)≡ζ6/(6-d). For such variation to occur, the amplitude modulating the scaling of ξ with u must itself become anomalous in ζ as soon as fluctuations set in, so that ξd∝ζΨ(d)u-νd. The exponent Ψ(d) can be found by invoking a smooth variation at u*, i.e., ξMF|u*=ξd|u*, yielding Ψ(d)=[(6-2dνd)/(6-d)]. Factoring out ξ*, as well as u/u*, we obtain the generalized scaling ξ(u,ζ)ξ*(ζ)=(uu*(ζ))-½F(uu*(ζ)),where F(x)∝x-νd+1/2 for x≪1, and constant otherwise. After specializing these results to the case d=2, the data presented in Fig. 1(a) have been accordingly rescaled. The data collapse depicted in Figs. 1(c) and 1(d) supports our scaling arguments, and further highlights an interesting prediction of the LG criterion. Although one would have expected the random-to-spatial crossover to occur at the input characteristic scale ζ—as it does far from criticality—here it stretches until ζ6/(6-d). This scale can be regarded as the network analogue at the crossover of what Binder calls “thermodynamic length” in spin systems [27]. However, in contrast to thermal phase transitions, ξ* identifies here a “structural length” at which geometric fluctuations set in. This means that ξ* does not only rule the critical scaling of the percolation observables, but also the network’s geometry at the onset of the giant component, stretching nonlinearly with ζ its mean-field regimes. For this reason, we call this phenomenon critical stretching and in what follows we characterize its geometric and dynamic properties.
Geometric stretching at criticality.—From a geometric standpoint, the crossover scaling in Eq. (1) implies a nontrivial structural variation of the large-scale organization of the network close to its percolation threshold. To put this on firmer bases, consider the scaling of the pair connectivity function g(r,u,ζ)∝r-d+2-ηe-r/ξ(u,ζ), where η is the anomalous dimension. The correlation length ξ(u,ζ) is itself the representative of a crossover in the giant component, separating a self-similar organization where g(r)∼r-d+2-η, from a “noncritical” one marked instead by a decreasing exponential [31]. Then again, we know from Eq. (1) that ξ(u,ζ) crosses over between the mean-field and the lattice singular behaviors, depending on how close the network is to pc. Let us then fix the occupation probability to a value u≪u*, so that geometric fluctuations are free to set in above ξ*. In turn, this means that ξ|u≪u*≫ξ* and so the cutoff e-r/ξ(u,ζ) becomes significant only if r∼O(ξ*). Therefore, reading gu(r,ζ) as a function of r for fixed u≪u* and variable values of ζ, we advance the crossover relation gu(r,ζ)g*(ζ)=(rξ*(ζ))-4G(rξ*(ζ),ξ|u≪u*ξ*(ζ)),where g*(ζ)≡ζ-24/(6-d), and G(x,y)∝xΦ(d)e-x/y for x≫1 and constant otherwise, with y≫1 and Φ(d)=6-d-η. To obtain Eq. (2), we have rescaled the argument of the exponential cutoff to make explicit the ratio r/ξ*, and factorized out the mean-field regime gMF∼r-4.
Equation (2) provides a geometric formulation of the LG criterion at criticality, upon which density fluctuations of the spanning cluster are negligible as far as r≪ζ6/(6-d). This result enables testing the critical stretching of the high-dimensional (mean-field) regime of the networks’ structure by comparing its geometric features at and away from criticality. With this aim, we focus on two exponents: (i) the shortest-path fractal dimension dmin, and (ii) the mass fractal dimension df of the incipient cluster [30]. Although the analysis could be carried out in any dimension d<6, we restrict here to the case d=2 for simplicity of computation. Let us start from the behavior at p=pc. (i) By definition, dmin regulates the scaling of the mean Euclidean distance ⟨r⟩ between sites in the giant cluster separated by the chemical distance ℓ, so that ℓ∼⟨r⟩dmin for ⟨r⟩ large enough. For future reference, we consider the reciprocal scaling ⟨r⟩∼ℓ1/dmin, for ℓ large enough. In chemical space, ξ* translates [35] into the quantity ξℓ*≡⌊u*-1⌋=⌊ζ⌋ (where the brackets denote the floor function) which, consistently with the scaling ansatz in Eq. (2), identifies the intrinsic chemical length at which the geometric crossover occurs. In defining dmin, we must thus consider two regimes: ℓ≪ξℓ* where dminMF=2 (i.e., random walk), and ℓ≫ξℓ* where instead dmin2D≃1.13077(2)[36]. Merging the two behaviors into a single equation, we propose the scaling ⟨r⟩(ℓ,ζ)ξ*(ζ)=(ℓξℓ*(ζ))½R(ℓξℓ*(ζ)),where R(x)∝x1/dmin2D-1/2 for x≫1 and constant otherwise. We tested Eq. (3) by using a breadth-first search algorithm whose results, displayed in Fig. 2(a), support our scaling arguments. (ii) To investigate the mass scaling M∼rdf at pc, we first need to determine the dependence on ζ of the giant cluster’s amplitude. This quantity is found by calculating the ζ correction to the size S∞ of the giant component in the spatial regime which, according to the LG criterion, scales as S∞2D∝ℬ(ζ)u5/36 with ℬ(ζ)≡ζ-31/36. In light of this, we obtain the scaling M2D∝ℬ(ζ)r91/48 for r≫ξ* which, combined with what expected in the high-dimensional (mean-field) part, i.e., MMF∼r4[30], yields the ansatz M(r,ζ)M*(ζ)=(rξ*(ζ))4M(rξ*(ζ)),where M(x)∝x-101/48 for x≫1, constant otherwise, and M*(ζ)≡ζ571/288 is the incipient cluster’s mass at ξ*. Equation (4) states in more practical terms the prediction of Eq. (2), namely that the fractal geometry of the critical cluster in our network crosses over from the one of an ER graph [37] (i.e., dfMF=4) to that of a d=2 lattice (for which df2D=91/48), where the former stretches up until the intrinsic length ζ3/2. As shown in Fig. 2(b), Eq. (4) exquisitely matches the numerical data.
110.1103/PhysRevLett.123.088301.f1
Crossover of the correlation length. (a) Below criticality ξ∼|p-pc|-νd, where the lattice exponent ν2D=4/3 (dashed line) appears at large distances (i.e., p-pc≪1), while the ER one νMF=1/2 (thick line) is found otherwise. For ζ large enough, the scaling of ξ crosses over between these extremes. (b) Shift of pc for increasing ζ; notice that pc≃0.25=pcMF already at ζ≃50. (c) Collapse of the data in (a) according to Eq. (1) with d=2. (d) Fitted slope: u*=ζ-ω with ω≃1.01, confirming the LG predictions. The curves are obtained for N=108 by averaging over 200 realizations.
210.1103/PhysRevLett.123.088301.f2
Geometric crossover and stretching. (a) Scaled data of ⟨r⟩(ℓ) at pc according to Eq. (3); linear fit (inset) yields ξℓ*=ζψ, ψ≃1.03. (b) Mass scaling at pc: scaled data are consistent with Eq. (4) and were obtained by adopting the spatial cluster growing algorithm introduced in Ref. [8]. (Inset) The fitted slope of ξ*=ζϕ yields ϕ≃1.48. (c) Mass scaling at p=1: for r≳ζ the network has the lattice’s geometry, while for r≲ζ it embeds onto a Euclidean space whose dimension increases with ζ. (Inset) Inverse maximal slope of the mass scaling in the small-world regime vs 1/log(ζ): the network’s dimension grows unbounded for ζ≫1. (d) Shortest-path scaling at p=1, and crossover length ξℓ×∝logζ (inset) in semilog scale. Simulation settings as in Fig. 1.
To complete the picture, we must verify that far from the percolation threshold, ξ×≡ζ is indeed the characteristic length at which the geometric crossover occurs. Let us start from the mass scaling (ii) at p=1. In this case, the Euclidean dimension d=2 is expected to regulate the scaling of M for distances r∼O(ξ×), while an infinite-dimensional (small-world) behavior should be observed for r∼o(ξ×). The numerical results in Fig. 2(c) corroborate this intuition, further showing [Fig. 2(c), inset] that below ξ× the random structure embeds onto an Euclidean space whose dimension grows unbounded with ζ. (i) The shortest-path’s scaling at p=1 supports this large-to small-world crossover: while for r≲ξ×any path shall behave like a random walk [30], i.e., ⟨r⟩∝ζℓ1/2, for ⟨r⟩≳ξ× (lattice regime) the lengths ℓ and ⟨r⟩ must share the same metric, i.e., ⟨r⟩∼ℓ. We can then advance the scaling ⟨r⟩=ζℓ1/2f(ℓ) where f(ℓ)∝ℓ1/2 for ℓ≳ξℓ×(ζ) and constant otherwise. The quantity ξℓ× can be regarded as the chemical length up to which the small-world extends, so ξℓ×∝logNζ∼logζ since Nζ∼ζ2. The data collapse in Fig. 2(d) is consistent with our scaling, and further supports [Fig. 2(d), inset] the approximation ξℓ×∝logζ. Altogether, these results confirm the prediction that far from criticality the network’s geometry is governed by the characteristic scale ζ, corroborating the surprising effect of the critical stretching.
Stretching effects in spreading dynamics.—Since collective phenomena are strongly affected by the geometry of the underlying interaction structure, the critical stretching will accordingly influence the dynamics of processes acting on the network. As an illustrative example, let us consider an SIR process with infection rate β and recovery time γ≡1, spreading on our network model in d=2. In light of the mapping to percolation [32], the clusters of nodes removed by the outbreak are geometrically equivalent to those originated from a static deletion, whereas the shortest-path ℓ identifies now the number of time steps t separating infected nodes in the ℓth shell from an initial spreader. Setting the relative rate β′≡1-β/βc to a value β′≪ζ-1, we thus expect the spreading of the disease to undergo a temporal crossover at the intrinsic time scale τ*≡⌊ζ⌋ from a mean-field behavior to a spatial one. In particular, interpreting ⟨r⟩ as the average distance of the infective nodes from an initial spreader, we can deduce from Eq. (3) that the spreading velocity v≔d⟨r⟩/dt∼t1/dmin-1 of the outbreak front will accordingly satisfy the temporal scaling v(t,ζ)v*(ζ)=(tτ*(ζ))-½V(tτ*(ζ)),where V(x)∝x1/dmin2D-1/2 for x≫1, constant otherwise, and v*(ζ)≡ζ1/2 is the crossover critical velocity. The data collapse in Fig. 3 validates our premises, confirming that, if critical, the epidemic front undergoes a spatiotemporal crossover at the distances (ξ*,τ*) from the infective spreader, varying from diffusive to superdiffusive propagations. Since the crossover time scale stretches from a logarithmic—τ×∝logζ far from criticality—to a linear dependence with ζ, we disclose an intriguing scenario: if critical, an outbreak will invade the system at “slow” (i.e., diffusive) rates for longer times, before turning superdiffusive. Besides of theoretical interest, these results are relevant for the analysis of spreading processes in real-world systems, e.g., in virus propagation on wireless networks [16], in the spreading of diseases through transportation systems [19,20] or of rumors in social networks [21], whose spatial constraints are well captured by the exponential wiring cost function here adopted.
Discussion.—In this Letter, we have investigated the effects that a tunable characteristic link length ζ has on the critical phenomena in spatial networks. We studied the spatiotemporal location where a crossover from random (mean-field) to spatial regimes takes place when considering the network away from (p≫pc) or close to (p→pc) its percolation threshold. We found that while this crossover occurs at ζ when p=1, it surprisingly extends until ζ6/(6-d) close to pc, signaling that the mean-field (ER) structure expands beyond the input characteristic link length ζ near criticality. We named this phenomenon critical stretching, and we demonstrated its geometric and dynamic features by analyzing the Euclidean correlation length ξ(p,ζ) and its chemical counterpart ξℓ(p,ζ), whose crossover scales are summarized in Table I.
310.1103/PhysRevLett.123.088301.f3
Dynamical stretching in SIR-epidemics. The scaling of the infected mass I (calculated as in Ref. [8]), undergoes the same geometric crossover and critical stretching of M as in Eq. (4). (Inset) Spreading velocity of the epidemic front wave, scaled according to Eq. (5). (Subinset) Linear fit of τ*=ζϖ yielding ϖ≈0.99. Results are obtained for N=2.5×107, and averaged over 104 realizations.
I10.1103/PhysRevLett.123.088301.t1
Universal length scales in spatial networks. Dependence of the crossover correlation length ξ and of its chemical counterpart ξℓ on the characteristic link length ζ and the dimension of the embedding space, d<6, close to (right) and away from (left) the network’s percolation threshold.
p≫pc
p→pc
ξcross
ζ
ζ6/(6-d)
ξℓcross
logζ
ζ2d/(6-d)
Some relevant comments are in order. (i) The critical stretching length scales ξ* and ξℓ* are univocally defined by the critical exponents of percolation and the dimension of the embedding space. As such, they are not influenced by the system’s finite size [35] and can be regarded as a set of universal spatiotemporal scales for any spatial network featuring a characteristic scale of link lengths and a homogeneous degree distribution. (ii) In light of the expressions in Table I, the stretching effects become stronger at higher embedding dimensions. In d=3 networks, e.g., the intrinsic scales ξ*=ξℓ*=ζ2 at pc are significantly larger if compared to the crossover scales ξ×=ζ and ξℓ×∝logζ away from criticality. This supports the possibility that the stretching effects we found can be experimentally observed in real-world systems. In fact, besides accurately modeling communication [16–18,21] and transport networks [19,20], recent evidences [23–26] support the existence of scale-invariant organization principles in mammal brains based on the exponential wiring cost function examined here. Our results offer then a fundamental keystone for understanding, e.g., the traffic flow [13] or the neural functioning—e.g., in reference to the observed activity before epileptic seizures [38]—and they raise the perspective of developing criteria (say, measuring the amount of stretching in the propagation of the action potentials for neural networks) for testing how close these systems are to their tipping points.
I. B. and B. G. contributed equally to this work. S. H. acknowledges financial support from the ISF, ONR, DTRA: HDTRA-1-10-1-0014, BSF-NSF: 2015781, ARO, the Israeli Ministry of Science, Technology and Space (MOST) in joint collaboration with the Japan Science Foundation (JSF), and the Italian Ministry of Foreign Affairs and International Cooperation (MAECI), and the Bar-Ilan University Center for Research in Applied Cryptography and Cyber Security. I. B. thanks S. V. Buldyrev, G. Sicuro, and M. C. Strinati for valuable discussions.
A.-L. Barabási and M. Pósfai, (Cambridge University Press, Cambridge, 2016); V. Latora, V. Nicosia, and G. Russo, (Cambridge University Press, Cambridge, England, 2017); M. Newman, (Oxford University Press, Oxford, 2018); G. Bianconi, (Oxford University Press, Oxford, 2018).S. Fortunato,
486, 75 (2010).PRPLCM0370-157310.1016/j.physrep.2009.11.002A.-L. Barabási, N. Gulbahce, and J. Loscalzo,
12, 56 (2011).1471-005610.1038/nrg2918S. N. Dorogovtsev, A. V. Goltsev, and J. F. F. Mendes,
80, 1275 (2008).RMPHAT0034-686110.1103/RevModPhys.80.1275K. Kosmidis, S. Havlin, and A. Bunde,
82, 48005 (2008); EULEEJ0295-507510.1209/0295-5075/82/48005T. Emmerich, A. Bunde, S. Havlin, G. Li, and D. Li,
87, 032802 (2013).PRESCM1539-375510.1103/PhysRevE.87.032802S. Bradde, F. Caccioli, L. Dall’Asta, and G. Bianconi,
104, 218701 (2010).PRLTAO0031-900710.1103/PhysRevLett.104.218701M. Barthélemy,
499, 1 (2011).PRPLCM0370-157310.1016/j.physrep.2010.11.002L. Daqing, K. Kosmidis, A. Bunde, and S. Havlin,
7, 481 (2011).NPAHAX1745-247310.1038/nphys1932M. Barthélemy,
19, 205 (2018).CRPOBN1631-070510.1016/j.crhy.2018.10.006A. Fornito, A. Zalesky, and E. Bullmore, (Academic Press, New York, NY, 2016).M. Breakspear,
20, 340 (2017).10.1038/nn.4497R. Lambiotte, V. D. Blondel, C. de Kerchove, E. Huens, C. Prieur, Z. Smoreda, and P. Van Dooren,
387A, 5317 (2008).PHYADX0378-437110.1016/j.physa.2008.05.014G. Zeng, D. Li, S. Guo, L. Gao, Z. Gao, H. Eugene Stanley, and S. Havlin,
116, 23 (2019).PNASA60027-842410.1073/pnas.1801545116Y. Yang, T. Nishikawa, and A. E. Motter,
358, eaan3184 (2017).SCIEAS0036-807510.1126/science.aan3184B. M. Waxman,
6, 1617 (1988).ISACEM0733-871610.1109/49.12889M. Haenggi,
27, 1029 (2009).10.1109/JSAC.2009.090902E. W. Zegura, K. L. Calvert, and S. Bhattacharjee,
2, 594 (1996).10.1109/INFCOM.1996.493353S.-H. Yook, H. Jeong, and A.-L. Barabási,
99, 13382 (2002).PNASA60027-842410.1073/pnas.172501399A. Halu, S. Mukherjee, and G. Bianconi,
89, 012806 (2014).PRESCM1539-375510.1103/PhysRevE.89.012806M. Danziger, L. M. Shekhtman, Y. Berezin, and S. Havlin,
115, 36002 (2016).EULEEJ0295-507510.1209/0295-5075/115/36002D. J. Watts, P. S. Dodds, and M. E. Newman,
296, 1302 (2002).SCIEAS0036-807510.1126/science.1070120E. Bullmore and O. Sporns,
13, 336 (2012).NRNAAN1471-003X10.1038/nrn3214N. T. Markov,
24, 17 (2012).10.1093/cercor/bhs270M. Ercsey-Ravasz, N. T. Markov, C. Lamy, D. C. Van Essen, K. Knoblauch, Z. Toroczkai, and H. Kennedy,
80, 184 (2013).NERNET0896-627310.1016/j.neuron.2013.07.036S. Horvát,
14, e1002512 (2016).10.1371/journal.pbio.1002512R. Gămănuţ, H. Kennedy, Z. Toroczkai, M. Ercsey-Ravasz, D. C. Van Essen, K. Knoblauch, and A. Burkhalter,
97, 698 (2018).NERNET0896-627310.1016/j.neuron.2017.12.037K. Binder and E. Luijten,
344, 179 (2001).PRPLCM0370-157310.1016/S0370-1573(00)00127-7A. Pelissetto and E. Vicari,
368, 549 (2002).PRPLCM0370-157310.1016/S0370-1573(02)00219-3A. Coniglio,
15, 3829 (1982).JPHAC50305-447010.1088/0305-4470/15/12/032A. Bunde and S. Havlin, (Springer Sc. & Business Media, Berlin, 1996).D. Stauffer and A. Aharony, (Taylor & Francis, London, 1994).R. Pastor-Satorras, C. Castellano, P. Van Mieghem, and A. Vespignani,
87, 925 (2015).RMPHAT0034-686110.1103/RevModPhys.87.925
We consider ζ∈[1,L/2), where ζmax=L/2 reflects the choice of periodic boundary conditions (and avoids links’ re-drawing), while the lower bound ζmin=1 describes the lattice-dominated regime, since ρ∼1/ad for L≫1, where a≡1 is the lattice spacing.
A. P. Levanyuk,
36, 810 (1959) ZETFA70044-4510[A. P. Levanyuk9, 571 (1959)]; V. L. Ginzburg,
2, 2031 (1960) [V. L. Ginzburg2, 1824 (1961)].See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevLett.123.088301 for details on the percolation LG criterion, and for the finite-size analysis of the critical stretching length scales.Z. Zhou, J. Yang, Y. Deng, and R. M. Ziff,
86, 061101 (2012).PRESCM1539-375510.1103/PhysRevE.86.061101R. Cohen and S. Havlin,
336A, 6 (2004).PHYADX0378-437110.1016/j.physa.2004.01.005N. Vladimirov, R. D. Traub, and Y. Tu,
6, e20536 (2011).POLNCL1932-620310.1371/journal.pone.0020536