Percolation on multifractal, scale-free weighted planar stochastic porous lattice

This paper introduces the Weighted Planar Stochastic Porous Lattice (WPSPL), a multifractal, scale-free porous substrate, and demonstrates through analytical and numerical methods that bond percolation on this lattice exhibits a continuous family of distinct universality classes with critical exponents that vary with porosity while satisfying the Rushbrooke inequality.

The Gist

Imagine you are a city planner trying to understand how traffic flows through a city. Usually, you'd study a perfect grid of streets, like Manhattan, where every block is the same size and every intersection looks identical. But real cities aren't perfect. They have parks, dead ends, huge skyscrapers, and tiny alleyways. They are messy, irregular, and full of "holes."

This paper introduces a new way to model that messy reality and studies how "connectivity" spreads through it. Here is the breakdown in simple terms:

1. The "Magic Square" Game (The Lattice)

The authors created a digital model called the Weighted Planar Stochastic Porous Lattice (WPSPL). Think of it as a game played with a giant square piece of paper.

• The Rules: You pick a piece of the paper (a "block") to cut. But you don't pick randomly; you are more likely to pick a big piece than a small one.
• The Cut: You slice that big piece into four smaller pieces.
• The Twist (Porosity): Here is the special part. After cutting, you flip a coin.
  • If it's heads (controlled by a number called $q$), you keep all four pieces.
  • If it's tails, you throw one of the four pieces away, leaving a hole (a void) in your city.
• The Result: As you repeat this thousands of times, you get a complex, jigsaw-puzzle-like map full of different-sized blocks and holes. It looks chaotic, but it follows strict mathematical rules.

2. Why This Map is Special (Multifractals & Scale-Free)

Most maps are simple. If you zoom in on a regular grid, it looks the same. But this "Magic Square" map is multifractal.

• The Analogy: Imagine a coastline. If you look at it from a satellite, it's jagged. If you zoom in with a drone, it's still jagged. If you zoom in with a microscope, it's still jagged. This map is like that, but even more complex. It has "self-similarity," meaning the pattern of big blocks and small blocks repeats itself at every level of zoom.
• The "Rich Get Richer" Effect: In this map, big blocks are more likely to get cut again, which means their neighbors get more connections. This creates a "scale-free" network, similar to how a few famous people on social media have millions of followers while most people have very few.

3. The Big Question: When Does the City Get "Connected"?

The authors wanted to know: At what point does this messy city become fully connected?

Imagine you are building a bridge between every pair of neighboring blocks. You start with no bridges. You slowly add bridges one by one.

• Early on: You have tiny islands of connected blocks.
• The Tipping Point: Suddenly, a "Giant Component" appears. One massive cluster of blocks connects from one side of the map to the other. This is called Percolation.

The authors asked: Does the amount of "holes" (porosity) in the city change when this tipping point happens?

4. The Surprising Discovery

In normal, perfect grids (like a chessboard), the tipping point happens at a specific, unchangeable number. But in this messy, hole-filled world, the tipping point changes depending on how many holes you have.

• More Holes ($q$ is low): It takes more bridges to connect the city because the holes break the paths.
• Fewer Holes ($q$ is high): It's easier to connect.

The "Universe" of Rules:
Usually, scientists believe that all 2D grids belong to the same "Universality Class"—meaning they all behave the same way near the tipping point, regardless of small details.

• The Breakthrough: This paper proves that for this specific type of messy, porous lattice, every different amount of porosity creates a completely new "Universe Class."
• The rules of how the city connects change continuously as you change the number of holes. It's like saying that every different type of terrain (sand, mud, rock) has its own unique physics for how water flows through it.

5. The "Thermodynamics" of a Non-Thermal System

To study this, the authors used a clever trick. They treated the probability of adding a bridge ($p$) like Temperature in a physical system.

• Low $p$ (Cold): Everything is frozen and disconnected.
• High $p$ (Hot): Everything is fluid and connected.
• The "Heat Capacity": They calculated how "disordered" the system gets as you add bridges. They found that even though this is just a geometric puzzle, it behaves exactly like a physical substance undergoing a phase change (like ice melting into water).

The Bottom Line

This paper shows that geometry matters more than we thought.

If you have a perfectly ordered world, the rules of connection are simple and predictable. But if you introduce disorder, holes, and irregular shapes (like in real-world porous rocks, social networks, or biological tissues), the rules change completely. The "critical point" where things connect becomes a sliding scale, not a fixed number.

In everyday language:
If you are trying to spread a rumor, a virus, or a fire, you can't just look at the average number of connections. You have to look at the shape of the network. If the network is full of holes and irregular sizes, the way the "infection" spreads is fundamentally different from a neat, orderly grid. This research gives us the math to predict exactly how that happens.

Technical Summary

Here is a detailed technical summary of the paper "Percolation on multifractal, scale-free weighted planar stochastic porous lattice" by Proshanto Kumar and Md. Kamrul Hassan.

1. Problem Statement

Percolation theory traditionally relies on regular spatial lattices (e.g., square or triangular grids) to model connectivity transitions. However, many real-world systems—such as porous media, epidemic spreading, and complex networks—exhibit intrinsic geometric disorder, heterogeneity, and dynamic evolution that regular lattices fail to capture. While previous studies introduced the Weighted Planar Stochastic Lattice (WPSL) to model scale-free, self-similar structures, it lacked the ability to model porosity (voids) and the specific geometric disorder found in porous materials.

The authors address the gap by investigating bond percolation on a new substrate: the Weighted Planar Stochastic Porous Lattice (WPSPL). The core problem is to determine how the introduction of tunable porosity (controlled by a parameter $q$) and the resulting multifractal, scale-free geometry affect critical phenomena, percolation thresholds, and universality classes.

2. Methodology

A. Construction of the WPSPL

The authors define the WPSPL through an iterative stochastic growth algorithm starting from a unit square:

1. Selection: At each step, a block is selected with a probability proportional to its area.
2. Subdivision: The selected block is divided into four smaller sub-blocks by two mutually perpendicular cuts.
3. Porosity Rule: One of the four new sub-blocks is subjected to a retention test. It is retained with probability $q$ or removed (creating a void) with probability $1-q$. The other three are always retained.
4. Iteration: The process repeats, generating a lattice with tunable porosity. When $q=1$, the lattice is non-porous (WPSL); when $0 < q < 1$, it is porous.

B. Analytical Framework

• Multifractal Analysis: The authors employ the kinetics of planar fragmentation and Mellin transforms to derive the evolution of the block-size distribution. They identify an infinite set of conserved quantities, proving the lattice is multifractal.
• Dynamic Scaling: They analyze the block-area distribution $C(a, t)$ to establish temporal self-similarity, showing that the distribution collapses onto a universal scaling function.
• Network Dual: The lattice is mapped to a dual network where block centers are nodes and shared boundaries are links. The authors analyze the degree distribution of this network.

C. Percolation Simulation

• Algorithm: Bond percolation is simulated on the dual network using the Newman-Ziff (N-Z) algorithm, which is efficient for determining critical thresholds.
• Thermodynamic Analogy: To define critical exponents, the authors map percolation variables to thermodynamics:
  • Occupation probability $p$ $\leftrightarrow$ External ordering field.
  • $1-p$$\leftrightarrow$ Effective temperature.
  • Shannon entropy of cluster configurations $\leftrightarrow$ Entropy (used to define specific heat).
  • Derivative of the order parameter w.r.t. $p$ $\leftrightarrow$ Susceptibility.
• Finite-Size Scaling (FSS): Critical exponents ($\alpha, \beta, \gamma$) are extracted by analyzing the scaling behavior of observables (spanning probability, order parameter, specific heat, susceptibility) across different system sizes and time steps.

3. Key Contributions

1. Introduction of WPSPL: The paper formalizes a new geometric model that combines stochastic fragmentation, scale-free coordination numbers, and tunable porosity.
2. Proof of Multifractality: Analytically demonstrates that the WPSPL possesses infinitely many nontrivial conserved quantities, leading to a family of intertwined multifractal spectra rather than a single fractal dimension.
3. Scale-Free Topology: Confirms that the dual network of the WPSPL exhibits a power-law degree distribution $P(k) \sim k^{-\gamma}$, where the exponent $\gamma$ depends on the porosity parameter $q$.
4. Universality Class Discovery: Establishes that percolation on the WPSPL does not belong to the standard universality class of 2D regular lattices. Instead, it defines a continuous family of distinct universality classes, where critical exponents vary continuously with $q$.
5. Validation of Scaling Relations: Despite the non-universality, the authors demonstrate that the Rushbrooke inequality ($\alpha + 2\beta + \gamma \geq 2$) is satisfied with near-equality, confirming the validity of static scaling hypotheses even in this highly disordered, multifractal system.

4. Key Results

• Percolation Threshold ($p_c$): The threshold increases with porosity ($q$). For example, $p_c \approx 0.389$ for $q=0.90$ and $p_c \approx 0.416$ for $q=0.85$. Higher porosity reduces connectivity, requiring a higher bond occupation probability to percolate.
• Critical Exponents: The exponents vary continuously with $q$:
  • Order Parameter ($\beta$): Decreases as $q$ decreases (e.g., $\beta \approx 0.170$ for $q=0.90$ vs. $0.194$for$q=1.0$).
  • Susceptibility ($\gamma$): Increases as $q$ decreases (e.g., $\gamma \approx 0.782$ for $q=0.90$ vs. $0.757$for$q=1.0$).
  • Specific Heat ($\alpha$): Increases as $q$ decreases (e.g., $\alpha \approx 0.931$ for $q=0.90$ vs. $0.861$for$q=1.0$).
• Dimensionality: The global (Hausdorff) dimension of the lattice skeleton is $d_f = \sqrt{3+q} - 1$. Even for the non-porous case ($q=1$), where $d_f=2$, the system lies outside the standard 2D universality class due to its scale-free coordination disorder.
• Data Collapse: The study provides robust numerical evidence (via data collapse plots) for the dynamic scaling of block sizes and the finite-size scaling of percolation observables.

5. Significance

This work significantly advances the understanding of critical phenomena in disordered and complex media.

• Beyond Regular Lattices: It challenges the assumption that 2D systems must belong to a single universality class, showing that geometric disorder and multifractality can tune critical exponents continuously.
• Modeling Real Porous Media: The WPSPL offers a more realistic geometric substrate for modeling fluid flow, transport, and connectivity in porous rocks, biological tissues, and fractured materials than regular grids.
• Theoretical Robustness: The finding that the Rushbrooke inequality holds despite the continuous variation of exponents reinforces the universality of scaling laws, even when the underlying geometry is non-Euclidean and multifractal.
• Methodological Impact: The successful application of thermodynamic analogies (entropy as specific heat, derivative of order parameter as susceptibility) to a non-thermal, geometric system provides a robust framework for analyzing criticality in other complex networks.

In summary, the paper demonstrates that geometric disorder, multifractality, and porosity are not just perturbations but fundamental drivers that can generate a continuous spectrum of universality classes, fundamentally altering the critical behavior of percolation transitions.