We investigate bond- and site-percolation models on several two-dimensional lattices numerically, by means of transfer-matrix calculations and Monte Carlo simulations. The lattices include the square, triangular, honeycomb kagome, and diced lattices with nearest-neighbor bonds, and the square lattice with nearest- and next-nearest-neighbor bonds. Results are presented for the bond-percolation thresholds of the kagome and diced lattices, and the site-percolation thresholds of the square, honeycomb, and diced lattices. We also include the bond- and site-percolation thresholds for the square lattice with nearest- and next-nearest-neighbor bonds. We find that corrections to scaling behave according to the second temperature dimension X_{t2}=4 predicted by the Coulomb gas theory and the theory of conformal invariance. In several cases there is evidence for an additional term with the same exponent, but modified by a logarithmic factor. Only for the site-percolation problem on the triangular lattice does such a logarithmic term appear to be small or absent. The amplitude of the power-law correction associated with X_{t2}=4 is found to be dependent on the orientation of the lattice with respect to the cylindrical geometry of the finite systems.
We investigate bond- and site-percolation models on several two-dimensional lattices numerically, by means of transfer-matrix calculations and Monte Carlo simulations. The lattices include the square, triangular, honeycomb kagome, and diced lattices with nearest-neighbor bonds, and the square lattice with nearest- and next-nearest-neighbor bonds. Results are presented for the bond-percolation thresholds of the kagome and diced lattices, and the site-percolation thresholds of the square, honeycomb, and diced lattices. We also include the bond- and site-percolation thresholds for the square lattice with nearest- and next-nearest-neighbor bonds. We find that corrections to scaling behave according to the second temperature dimension X_{t2}=4 predicted by the Coulomb gas theory and the theory of conformal invariance. In several cases there is evidence for an additional term with the same exponent, but modified by a logarithmic factor. Only for the site-percolation problem on the triangular lattice does such a logarithmic term appear to be small or absent. The amplitude of the power-law correction associated with X_{t2}=4 is found to be dependent on the orientation of the lattice with respect to the cylindrical geometry of the finite systems.
We study the site-bond-correlated percolation model in two and three dimensions. In this model, two first-neighbor sites are considered to be joined by an active bond if both are present and if their first neighbors, along the line joining them, are also present. Using a Monte Carlo renormalization-group approach, we calculate percolation thresholds and critical exponents for both square and simple cubic lattices. Our results indicate that the correlated model is in the same universality class as usual percolation, in two and three dimensions; the critical probabilities for the square and simple cubic lattices are found to be 0.7402\ifmmode\pm\else\textpm\fi{}0.0002 and 0.5660\ifmmode\pm\else\textpm\fi{}0.0004, respectively.
We investigate the site-percolation problem on the square and simple-cubic lattices by means of a Monte Carlo algorithm that in fact simulates systems with size ${L}^{d\ensuremath{-}1}\ifmmode\times\else\texttimes\fi{}\ensuremath{\infty}$, where L specifies the linear system size. This algorithm can be regarded either as an extension of the Hoshen-Kopelman method or as a special case of the transfer-matrix Monte Carlo technique. Various quantities, such as the magnetic correlation function, are sampled in the finite directions of the above geometry. Simulations are arranged such that both bulk and surface quantities can be sampled. On the square lattice, we locate the percolation threshold at pc=0.5927465(4), and determine two universal quantities as Qgbc=0.93034(1) and Qgsc=0.79372(3), which are associated with bulk and surface correlations, respectively. These values agree well with the exact values ${2}^{\ensuremath{-}5∕48}$ and ${2}^{\ensuremath{-}1∕3}$, respectively, which follow from conformal invariance. On the simple-cubic lattice, we locate the percolation threshold at pc=0.3116077(4). We further determine the bulk thermal and magnetic exponents as yt=1.1437(6) and yh=2.5219(2), respectively, and the surface magnetic exponent at the ordinary phase transition as yhs(o)=1.0248(3).
In the preceding paper, one of us (F. Y. Wu) considered the Potts model and bond and site percolation on two general classes of two-dimensional lattices, the triangular-type and kagome-type lattices, and obtained closed-form expressions for the critical frontier with applications to various lattice models. For the triangular-type lattices Wu's result is exact, and for the kagome-type lattices Wu's expression is under a homogeneity assumption. The purpose of the present paper is twofold: First, an essential step in Wu's analysis is the derivation of lattice-dependent constants A,B,C for various lattice models, a process which can be tedious. We present here a derivation of these constants for subnet networks using a computer algorithm. Second, by means of a finite-size scaling analysis based on numerical transfer matrix calculations, we deduce critical properties and critical thresholds of various models and assess the accuracy of the homogeneity assumption. Specifically, we analyze the q -state Potts model and the bond percolation on the 3-12 and kagome-type subnet lattices (n×n):(n×n) , n≤4 , for which the exact solution is not known. Our numerical determination of critical properties such as conformal anomaly and magnetic correlation length verifies that the universality principle holds. To calibrate the accuracy of the finite-size procedure, we apply the same numerical analysis to models for which the exact critical frontiers are known. The comparison of numerical and exact results shows that our numerical values are correct within errors of our finite-size analysis, which correspond to 7 or 8 significant digits. This in turn infers that the homogeneity assumption determines critical frontiers with an accuracy of 5 decimal places or higher. Finally, we also obtained the exact percolation thresholds for site percolation on kagome-type subnet lattices (1×1):(n×n) for 1≤n≤6 .
We investigate the influence of the range of interactions in the two-dimensional bond percolation model by means of Monte Carlo simulations. We locate the phase transitions for several interaction ranges, as expressed by the number z of equivalent neighbors. We also consider the $z\ensuremath{\rightarrow}\ensuremath{\infty}$ limit, i.e., the complete graph case, where percolation bonds are allowed between each pair of sites, and the model becomes mean-field-like. All investigated models with finite z are found to belong to the short-range universality class. There is no evidence of a tricritical point separating the short-range and long-range behavior, such as is known to occur for q=3 and q=4 Potts models. We determine the renormalization exponent describing a finite-range perturbation at the mean-field limit as ${y}_{r}\ensuremath{\approx}2/3$. Its relevance confirms the continuous crossover from mean-field percolation universality to short-range percolation universality. For finite interaction ranges, we find approximate relations between the coordination numbers and the amplitudes of the leading correction terms as found in the finite-size scaling analysis.
We theoretically investigate the quantum percolation problem on Lieb lattices in two and three dimensions. We study the statistics of the energy levels through random matrix theory and determine the level spacing distributions, which, with the aid of finite-size scaling theory, allows us to obtain accurate estimates for site- and bond-percolation thresholds and critical exponents. Our numerical investigation supports a localized-delocalized transition at finite threshold, which decreases as the average coordination number increases. The precise determination of the localization length exponent enables us to claim that quantum site- and bond-percolation problems on Lieb lattices belong to the same universality class, with ν decreasing with lattice dimensionality d, similarly to the classical percolation problem. In addition, we verify that, in three dimensions, quantum percolation on Lieb lattices belongs to the same universality class as the Anderson impurity model.
Using suitable Monte Carlo methods and finite-size scaling, we investigate the Blume-Capel model on the square lattice. We construct percolation clusters by placing nearest-neighbor bonds between vacancies with a variable bond probability p(b) . At the tricritical point, we locate the percolation threshold of these vacancy clusters at p(bc) =0.706 33 (6) . At this point, we determine the fractal dimension of the vacancy clusters as Xf =0.1308 (5) approximately equal to 21/160, and the exponent governing the renormalization flow in the p(b) direction as y(p) =0.426 (2) approximately equal to 17/40 . For bond probability p(b) > p(bc) , the vacancy clusters maintain strong critical correlations; the fractal dimension is Xf =0.0750 (2) approximately equal to 3/40 and the leading correction exponent is y(p) =-0.45 (2) approximately equal to -19/40 . The above values fit well in the Kac table for the tricritical Ising model. These vacancy clusters have much analogy with those consisting of Ising spins of the same sign, although the associated quantities rho and magnetization m are energylike and magnetic quantities, respectively. However, along the critical line of the Blume-Capel model, the vacancies are more or less uniformly distributed over the whole lattice. In this case, no critical percolation correlations are observed in the vacancy clusters, at least in the physical region p(b) < or = 1 .