Random site percolation on honeycomb lattices with complex neighborhoods

Krzysztof Malarz

DOI: https://doi.org/10.1063/5.0099066

2024-02-12

Abstract:We present a rough estimation -- up to four significant digits, based on the scaling hypothesis and the probability of belonging to the largest cluster vs. the occupation probability -- of the critical occupation probabilities for the random site percolation problem on a honeycomb lattice with complex neighborhoods containing sites up to the fifth coordination zone. There are 31 such neighborhoods with their radius ranging from one to three and containing from three to 24 sites. For two-dimensional regular lattices with compact extended-range neighborhoods, in the limit of the large number $z$ of sites in the neighborhoods, the site percolation thresholds $p_c$ follow the dependency $p_c\propto 1/z$, as recently shown by Xun, Hao and Ziff [Physical Review E 105, 024105 (2022)]. On the contrary, noncompact neighborhoods (with holes) destroy this dependence due to the degeneracy of the percolation threshold (several values of $p_c$ corresponding to the same number $z$ of sites in the neighborhoods). An example of a single-value index $\zeta=\sum_i z_i r_i$ -- where $z_i$ and $r_i$ are the number of sites and radius of the $i$-th coordination zone, respectively -- characterizing the neighborhood and allowing avoiding the above-mentioned degeneracy is presented. The percolation threshold obtained follows the inverse square root dependence $p_c\propto 1/\sqrt\zeta$. The functions boundaries() (written in C) for basic neighborhoods (for the unique coordination zone) for the Newman and Ziff algorithm [Physical Review E 64, 016706 (2001)] are also presented.

Statistical Mechanics,Computational Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the critical occupation probabilities in the random site percolation problem with complex neighborhood structures on the honeycomb lattice. Specifically, the author estimates the critical occupation probabilities on the honeycomb lattice in complex neighborhoods containing up to the 5th coordination zone by numerical methods. The main objectives of the study include:

1. **Estimating the critical occupation probability**: For 31 different complex neighborhoods containing up to the 5th coordination zone, estimate their critical occupation probability \( p_c \). The radii of these neighborhoods range from 1 to 3, and the number of sites they contain ranges from 3 to 24.

2. **Verifying the validity of the index \( \xi \)**: Check whether a degeneracy phenomenon of the critical occupation probability \( p_c \) with respect to the total number of sites \( z \) in the neighborhood can be eliminated by introducing a scalar index \( \xi \). The index \( \xi \) is defined as:
\[
\xi = \sum_{i} \frac{z_i r_i^2}{i}
\]
where \( z_i \) and \( r_i \) are the number of sites and the radius of the \( i \)-th coordination zone, respectively.

3. **Proposing a new index \( \zeta \)**: If the index \( \xi \) is not applicable to the honeycomb lattice, then propose a new index \( \zeta \), defined as:
\[
\zeta = \sum_{i} z_i r_i
\]
and verify whether \( p_c \) is proportional to the reciprocal square root of \( \zeta \), that is:
\[
p_c \propto \frac{1}{\sqrt{\zeta}}
\]

Through these objectives, the author hopes to provide an effective method to distinguish the critical occupation probabilities under different complex neighborhoods and provide new insights for the study of percolation theory.