Faster calculation of the percolation correlation length on spatial networks

Michael M. Danziger, Bnaya Gross, and Sergey V. Buldyrev

Phys. Rev. E 101, 013306 – Published 13 January 2020

Abstract

The divergence of the correlation length $ξ$ at criticality is an important phenomenon of percolation in two-dimensional systems. Substantial speed-ups to the calculation of the percolation threshold and component distribution have been achieved by utilizing disjoint sets, but existing algorithms of this sort cannot measure the correlation length. Here we utilize the parallel axis theorem to track the correlation length as nodes are added to the system, allowing us to utilize disjoint sets to measure $ξ$ for the entire percolation process with arbitrary precision in a single sweep. This algorithm enables direct measurement of the correlation length in lattices as well as spatial network topologies and provides an important tool for understanding critical phenomena in spatial systems.

Received 26 March 2019

DOI:https://doi.org/10.1103/PhysRevE.101.013306

Physics Subject Headings (PhySH)

Continuous percolation transition Network phase transitions Network structure Percolation Percolation phase transition

Networks & random structures

Networks Analysis Tools Percolation theory

NetworksGeneral PhysicsStatistical Physics & ThermodynamicsCondensed Matter, Materials & Applied Physics

Authors & Affiliations

Michael M. Danziger ¹, Bnaya Gross ², and Sergey V. Buldyrev³

¹Network Science Institute, Northeastern University, Boston, Massachusetts 02115, USA
²Department of Physics, Bar Ilan University, Ramat Gan 5290002, Israel
³Department of Physics, Yeshiva University, New York, New York 10033, USA

Article Text (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand

Issue

Vol. 101, Iss. 1 — January 2020

Reuse & Permissions

Access Options

Author publication services for translation and copyediting assistance advertisement

Images

Figure 1
Tracking inertia and radius of gyration as links are added. In the left panel, there are two clusters, with known inertia, and centers of mass marked with an “x”. In the right panel, a new link is added which joins the two clusters. Instead of calculating the inertia from scratch, it can be updated using the parallel axis theorem as shown. In this way, as links (or nodes) are added to the network we can continually track the inertia, the radius of gyration, and, by averaging the radius of gyration across clusters, the correlation length $ξ$ .
Reuse & Permissions
Figure 2
Node renumbering in light of periodic boundary conditions. Here a $4 \times 4$ lattice with periodic boundary conditions is shown, along with its replicas on all sides. The figure shows a system with cluster structure defined by the solid links, as the dashed links (a), (b), and (c) are added. Each node has two coordinates: the upper “lattice” coordinate which is based on the lattice location and not affected by the cluster structure, and a lower “relative” coordinate which is updated to maintain consistent distance measurements within clusters. Link (a) will require no renumbering because neither previous cluster (the single nodes) crosses the boundary. In this case, correction $C = [(1, 3) - (1, 2)] - [(1, 3) - (1, 2)] = (0, 0)$ . Link (b) will cause the smaller cluster (in blue, on the left) to be renumbered to match the larger cluster (in green, on the right) which crosses a boundary. In this case the correction $C = [(- 2, 3) - (2, 2)] - [(2, 3) - (2, 2)] = - (4, 0)$ . For link (c) correction $C = (- 4, 0)$ , the same as for link (b), but because link (c) connects the cluster to itself, the renumbering is not possible, which means that now this cluster spans the system from top to bottom and, hence, must be excluded from the calculation of the correlation length.
Reuse & Permissions
Figure 3
The percolation correlation length $ξ$ as a function of $p_{c} - p$ below criticality ( $p < p_{c}$ ). It is seen that both the burning and disjoint sets algorithms give similar results showing the expected scaling of $ξ \sim {(p_{c} - p)}^{- ν}$ with $ν = 4 / 3$ . The inset shows $ξ$ as function of $p$ above and below the scaling region. Simulations shown for $N = 10^{6}$ on a next-nearest lattice ( $z = 8$ ) averaged over 1000 realizations.
Reuse & Permissions
Figure 4
Running time as function of the resolution. The continuous lines represent the computational running time of the burning algorithm for a given system size, and the dashed lines represent the running time of our disjoint sets algorithm for the same system size. One can see that for resolution higher than 10, our disjoint sets algorithm is faster, and the running time remains constant while the burning algorithm running time increases polynomially.
Reuse & Permissions
Figure 5
Convolved results $ξ_{p}$ compared to raw results $ξ_{n}$ . As expected by the theory, the slope of 4/3 is more closely matched by our numerical measurements when they are convolved with the binomial distribution. Shown here is an average of 100 runs for $L = 1000$ and $z = 4$ , with sampling step size of 2.
Reuse & Permissions
Figure 6
Measurement of the correlation length in a spatial network with exponential link lengths. In contrast to previous optimizations, the algorithm presented here can efficiently measure the correlation length in networks with spatial disorder. The network model was introduced in Ref. [15] and further analyzed in Refs. [16, 17, 19].
Reuse & Permissions

Physical Review E

covering statistical, nonlinear, biological, and soft matter physics