Xu Wenhui, Wang Junfeng, Hu Hao, Deng Youjin
School of Physics and Materials Science, Anhui University, Hefei, Anhui 230601, China.
Hefei National Laboratory for Physical Sciences at the Microscale and Department of Modern Physics, University of Science and Technology of China, Hefei, Anhui 230026, China.
Phys Rev E. 2021 Feb;103(2-1):022127. doi: 10.1103/PhysRevE.103.022127.
Exact or precise thresholds have been intensively studied since the introduction of the percolation model. Recently, the critical polynomial P_{B}(p,L) was introduced for planar-lattice percolation models, where p is the occupation probability and L is the linear system size. The solution of P_{B}=0 can reproduce all known exact thresholds and leads to unprecedented estimates for thresholds of unsolved planar-lattice models. In two dimensions, assuming the universality of P_{B}, we use it to study a nonplanar lattice model, i.e., the equivalent-neighbor lattice bond percolation, and the continuum percolation of identical penetrable disks, by Monte Carlo simulations and finite-size scaling analysis. It is found that, in comparison with other quantities, P_{B} suffers much less from finite-size corrections. As a result, we obtain a series of high-precision thresholds p_{c}(z) as a function of coordination number z for equivalent-neighbor percolation with z up to O(10^{5}) and clearly confirm the asymptotic behavior zp_{c}-1∼1/sqrt[z] for z→∞. For the continuum percolation model, we surprisingly observe that the finite-size correction in P_{B} is unobservable within uncertainty O(10^{-5}) as long as L≥3. The estimated threshold number density of disks is ρ_{c}=1.43632505(10), slightly below the most recent result ρ_{c}=1.43632545(8) of Mertens and Moore obtained by other means. Our work suggests that the critical polynomial method can be a powerful tool for studying nonplanar and continuum systems in statistical mechanics.
自从渗流模型被引入以来,精确阈值就受到了深入研究。最近,针对平面晶格渗流模型引入了临界多项式(P_{B}(p,L)),其中(p)是占据概率,(L)是线性系统大小。(P_{B}=0)的解可以重现所有已知的精确阈值,并对未解决的平面晶格模型的阈值给出前所未有的估计。在二维情况下,假设(P_{B})具有通用性,我们通过蒙特卡罗模拟和有限尺寸标度分析,用它来研究一个非平面晶格模型,即等效邻域晶格键渗流,以及相同可穿透圆盘的连续渗流。结果发现,与其他量相比,(P_{B})受有限尺寸修正的影响要小得多。因此,我们得到了一系列高精度的阈值(p_{c}(z)),它是等效邻域渗流的配位数(z)的函数,(z)高达(O(10^{5})),并清楚地证实了(z→∞)时(zp_{c}-1∼1/\sqrt{z})的渐近行为。对于连续渗流模型,我们惊人地发现,只要(L≥3),在不确定性(O(10^{-5}))范围内,(P_{B})中的有限尺寸修正就不可观测。估计的圆盘阈值数密度为(\rho_{c}=1.43632505(10)),略低于Mertens和Moore通过其他方法得到的最新结果(\rho_{c}=1.43632545(8))。我们的工作表明,临界多项式方法可以成为研究统计力学中非平面和连续系统的有力工具。