Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical Review E

Finite-size scaling of energylike quantities in percolation

Youjin Deng, Xianqing Yang

Published Jun 13, 2006• 4 citations
Publisher

Abstract

We study the bond-percolation model in two and three dimensions by Monte Carlo simulation, and investigate the finite-size scaling behavior of several quantities that account for fluctuations of the total numbers of clusters and occupied bonds, ${N}_{c}$ and ${N}_{b}$, respectively. These quantities include ${C}_{2c}=(⟨{N}_{c}^{2}⟩\ensuremath{-}{⟨{N}_{c}⟩}^{2})∕{L}^{d}$ and ${C}_{cb}=(⟨{N}_{c}{N}_{b}⟩\ensuremath{-}⟨{N}_{c}⟩⟨{N}_{b}⟩)∕{L}^{d}$, where $L$ is the linear system size and $d$ is the spatial dimensionality. In statistical models with thermal fluctuations, ${C}_{2c}$ and ${C}_{cb}$ are specific heatlike quantities. Despite the absence of thermal fluctuations in percolation, we find that the leading finite-size scaling of ${C}_{2c}$ and ${C}_{cb}$ is described by the thermal critical exponent ${y}_{t}\text{\ensuremath{-}}d$. We also measure quantity ${\ensuremath{\kappa}}_{b}=2⟨{N}_{b}{S}_{2}⟩∕⟨{S}_{2}⟩\ensuremath{-}⟨{N}_{b}{\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{S}}_{4}⟩∕⟨{\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{S}}_{4}⟩\ensuremath{-}⟨{N}_{b}⟩$ and an analogous quantity ${\ensuremath{\kappa}}_{c}$ for ${N}_{c}$, where ${S}_{2}$ and ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{S}}_{4}$ are quantities associated with the second and the fourth moments of cluster sizes, respectively. At criticality, we show that ${\ensuremath{\kappa}}_{b}$ and ${\ensuremath{\kappa}}_{c}$ diverge as ${L}^{{y}_{t}}$ for $L\ensuremath{\rightarrow}\ensuremath{\infty}$. The analysis of the data of ${\ensuremath{\kappa}}_{b}$ and ${\ensuremath{\kappa}}_{c}$ yields ${y}_{t}=1.145(2)$ for the three-dimensional percolation, in good agreement with existing results.

DOI
10.1103/physreve.73.066116

Related papers

Percolation in the canonical ensemble

Hao Hu, Henk W. J. Blöte, Youjin Deng

We study the bond percolation problem under the constraint that the total number of occupied bonds is fixed, so that the canonical ensemble applies. We show via an analytical approach that at criticality, the constraint can induce new finite-size corrections with exponent y can = 2 y t − d both in energy-like and magnetic quantities, where y t = 1/ν is the thermal renormalization exponent and d is the spatial dimension. Furthermore, we find that while most of the universal parameters remain unchanged, some universal amplitudes, like the excess cluster number, can be modified and become non-universal. We confirm these predictions by extensive Monte Carlo simulations of the two-dimensional percolation problem which has y can = −1/2. This article is part of 'Lattice models and integrability', a special issue of Journal of Physics A: Mathematical and Theoretical in honour of F Y Wu's 80th birthday.

arXiv:1210.3463•Journal of Physics A Mathematical and Theoretical•Nov 27, 2012•20 citations

Bond and site percolation in three dimensions

Junfeng Wang, Zongzheng Zhou, Wei Zhang +2 more

We simulate the bond and site percolation models on a simple-cubic lattice with linear sizes up to $L=512$, and estimate the percolation thresholds to be ${p}_{c}(\mathrm{bond})=0.248\phantom{\rule{0.16em}{0ex}}811\phantom{\rule{0.16em}{0ex}}82(10)$ and ${p}_{c}(\mathrm{site})=0.311\phantom{\rule{0.16em}{0ex}}607\phantom{\rule{0.16em}{0ex}}7(2)$. By performing extensive simulations at these estimated critical points, we then estimate the critical exponents $1/\ensuremath{\nu}=1.1410(15)$, $\ensuremath{\beta}/\ensuremath{\nu}=0.477\phantom{\rule{0.16em}{0ex}}05(15)$, the leading correction exponent ${y}_{i}=\ensuremath{-}1.2(2)$, and the shortest-path exponent ${d}_{\mathrm{min}}=1.3756(3)$. Various universal amplitudes are also obtained, including wrapping probabilities, ratios associated with the cluster-size distribution, and the excess cluster number. We observe that the leading finite-size corrections in certain wrapping probabilities are governed by an exponent $\ensuremath{\approx}$$\ensuremath{-}2$, rather than ${y}_{i}\ensuremath{\approx}\ensuremath{-}1.2$.

arXiv:1302.0421•Physical Review E•May 7, 2013•170 citations

Critical behavior of the geometrical spin clusters and interfaces in the two-dimensional thermalized bond Ising model

S Davatolhagh, M Moshfeghian, A A Saberi

The fractal dimensions and the percolation exponents of geometrical spin clusters of like sign at criticality are obtained numerically for an Ising model with temperature dependent annealed bond dilution, also known as the thermalized bond Ising model (TBIM), in two dimensions. For this purpose, a modified Wolff single-cluster Monte Carlo simulation is used to generate equilibrium spin configurations on square lattices in the critical region. A tie-breaking rule is employed to identify nonintersecting spin cluster boundaries along the edges of the dual lattice. The values obtained for the fractal dimensions of the spanning geometrical clusters D c , and their interfaces D I , are in perfect agreement with those reported for the standard two-dimensional ferromagnetic Ising model. Furthermore, the variance of the winding angles results in a diffusivity κ = 3 for the two-dimensional thermalized bond Ising model, thus placing it in the universality class of the regular Ising model. A finite-size scaling analysis of the largest geometrical clusters results in a reliable estimation of the critical percolation exponents for the geometrical clusters in the limit of an infinite lattice size. The percolation exponents thus obtained are also found to be consistent with those reported for the regular Ising model. These consistencies are explained in terms of the Fisher renormalization relations, which express the thermodynamic critical exponents of systems with annealed bond dilution in terms of those of the regular model system.

arXiv:1202.6683•Journal of Statistical Mechanics Theory and Experiment•Feb 21, 2012•2 citations

Critical droplets and phase transitions in two dimensions

Santo Fortunato

In two space dimensions, the percolation point of the pure-site clusters of the Ising model coincides with the critical point ${T}_{c}$ of the thermal transition and the percolation exponents belong to a special universality class. By introducing a bond probability ${p}_{B}<1,$ the corresponding site-bond clusters keep on percolating at ${T}_{c}$ and the exponents do not change, until ${p}_{B}{=p}_{\mathrm{CK}}=1\ensuremath{-}\mathrm{exp}(\ensuremath{-}2J/kT):$ for this special expression of the bond weight the critical percolation exponents switch to the two-dimensional Ising universality class. We show here that the result is valid for a wide class of bidimensional models with a continuous magnetization transition: there is a critical bond probability ${p}_{c}$ such that, for any ${p}_{B}>~{p}_{c},$ the onset of percolation of the site-bond clusters coincides with the critical point of the thermal transition. The percolation exponents are the same for ${p}_{c}<{p}_{B}<~1$ but, for ${p}_{B}{=p}_{c},$ they suddenly change to the thermal exponents, so that the corresponding clusters are critical droplets of the phase transition. Our result is based on Monte Carlo simulations of various systems near criticality.

arXiv:cond-mat/0205211•Physical review. B, Condensed matter•Jan 15, 2003•16 citations

Surface critical phenomena in three-dimensional percolation

Youjin Deng, Henk W. J. Blöte

Using Monte Carlo methods and finite-size scaling, we investigate surface critical phenomena in the bond-percolation model on the simple-cubic lattice with two open surfaces in one direction. We decompose the whole lattice into percolation clusters and sample the surface and bulk dimensionless ratios Q1 and Qb, defined on the basis of the moments of the cluster-size distribution. These ratios are used to determine critical points. At the bulk percolation threshold pbc, we determine the surface bond-occupation probability at the special transition as p(s)1c = 0.418 17(2), and further obtain the corresponding surface thermal and magnetic exponents as y(s)t1 = 0.5387(2) and y(s)h1 = 1.8014(6), respectively. Next, from the pair correlation function on the surfaces, we determine y(o)h1 = 1.0246(4) and y(e)h1 = 1.25(6) for the ordinary and the extraordinary transition, respectively, of which the former is consistent with the existing result y(o)h1 = 1.024(4). We also numerically derive the line of surface phase transitions occurring at pb < pbc, and determine the pertinent asymptotic values of the universal ratios Q1 and Qb.

Physical Review E•Jan 12, 2005•11 citations

Logarithmic finite-size scaling of the four-dimensional Ising model

Zhiyi Li, Tianning Xiao, Zongzheng Zhou +2 more

Field-theoretical calculations predict that, at the upper critical dimension ${d}_{c}=4$, the finite-size scaling (FSS) behaviors of the Ising model would be modified by multiplicative logarithmic corrections with thermal and magnetic correction exponents $({\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{y}}_{t},{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{y}}_{h})=(1/6,1/4)$. Using high-efficient cluster algorithms and the lifted worm algorithm, we present a systematic study to the FSS of the four-dimensional Ising model at criticality in the Fortuin-Kasteleyn (FK) bond and loop representations. In the FK representation, the size of the largest cluster is observed to scale as ${C}_{1}\ensuremath{\sim}{L}^{3}{(lnL)}^{{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{y}}_{h}}$, while the size of the second-largest cluster scales as ${C}_{2}\ensuremath{\sim}{L}^{3}{(lnL)}^{{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{y}}_{h2}}$ with ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{y}}_{h2}=\ensuremath{-}1/4$ a new correction exponent not yet predicted from field theory. In the loop representation, we observe that the size of the largest loop cluster scales as ${F}_{1}\ensuremath{\sim}{L}^{2}{(lnL)}^{{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{y}}_{t}}$, and the specific heat scales as ${c}_{\mathrm{E}}\ensuremath{\sim}{(lnL)}^{2{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{y}}_{t}}$. This clarifies the long-standing open question that whether the specific heat for the critical Ising model at ${d}_{c}=4$ diverges logarithmically.

arXiv:2408.15230•Physical review. E•Dec 20, 2024•2 citations
We study the bond-percolation model in two and three dimensions by Monte Carlo simulation, and investigate the finite-size scaling behavior of several quantities that account for fluctuations of the total numbers of clusters and occupied bonds, Nc{N}_{c}Nc​ and Nb{N}_{b}Nb​, respectively. These quantities include C2c=(⟨Nc2⟩\ensuremath−⟨Nc⟩2)∕Ld{C}_{2c}=(⟨{N}_{c}^{2}⟩\ensuremath{-}{⟨{N}_{c}⟩}^{2})∕{L}^{d}C2c​=(⟨Nc2​⟩\ensuremath−⟨Nc​⟩2)∕Ld and Ccb=(⟨NcNb⟩\ensuremath−⟨Nc⟩⟨Nb⟩)∕Ld{C}_{cb}=(⟨{N}_{c}{N}_{b}⟩\ensuremath{-}⟨{N}_{c}⟩⟨{N}_{b}⟩)∕{L}^{d}Ccb​=(⟨Nc​Nb​⟩\ensuremath−⟨Nc​⟩⟨Nb​⟩)∕Ld, where LLL is the linear system size and ddd is the spatial dimensionality. In statistical models with thermal fluctuations, C2c{C}_{2c}C2c​ and Ccb{C}_{cb}Ccb​ are specific heatlike quantities. Despite the absence of thermal fluctuations in percolation, we find that the leading finite-size scaling of C2c{C}_{2c}C2c​ and Ccb{C}_{cb}Ccb​ is described by the thermal critical exponent yt\ensuremath-d{y}_{t}\text{\ensuremath{-}}dyt​\ensuremath-d. We also measure quantity \ensuremathκb=2⟨NbS2⟩∕⟨S2⟩\ensuremath−⟨NbS\ifmmode~\else˜\fi4⟩∕⟨S\ifmmode~\else˜\fi4⟩\ensuremath−⟨Nb⟩{\ensuremath{\kappa}}_{b}=2⟨{N}_{b}{S}_{2}⟩∕⟨{S}_{2}⟩\ensuremath{-}⟨{N}_{b}{\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{S}}_{4}⟩∕⟨{\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{S}}_{4}⟩\ensuremath{-}⟨{N}_{b}⟩\ensuremathκb​=2⟨Nb​S2​⟩∕⟨S2​⟩\ensuremath−⟨Nb​S\ifmmode~\else˜\fi​4​⟩∕⟨S\ifmmode~\else˜\fi​4​⟩\ensuremath−⟨Nb​⟩ and an analogous quantity \ensuremathκc{\ensuremath{\kappa}}_{c}\ensuremathκc​ for Nc{N}_{c}Nc​, where S2{S}_{2}S2​ and S\ifmmode~\else˜\fi4{\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{S}}_{4}S\ifmmode~\else˜\fi​4​ are quantities associated with the second and the fourth moments of cluster sizes, respectively. At criticality, we show that \ensuremathκb{\ensuremath{\kappa}}_{b}\ensuremathκb​ and \ensuremathκc{\ensuremath{\kappa}}_{c}\ensuremathκc​ diverge as Lyt{L}^{{y}_{t}}Lyt​ for L\ensuremath→\ensuremath∞L\ensuremath{\rightarrow}\ensuremath{\infty}L\ensuremath→\ensuremath∞. The analysis of the data of \ensuremathκb{\ensuremath{\kappa}}_{b}\ensuremathκb​ and \ensuremathκc{\ensuremath{\kappa}}_{c}\ensuremathκc​ yields yt=1.145(2){y}_{t}=1.145(2)yt​=1.145(2) for the three-dimensional percolation, in good agreement with existing results.
We study the bond percolation problem under the constraint that the total number of occupied bonds is fixed, so that the canonical ensemble applies. We show via an analytical approach that at criticality, the constraint can induce new finite-size corrections with exponent y can = 2 y t − d both in energy-like and magnetic quantities, where y t = 1/ν is the thermal renormalization exponent and d is the spatial dimension. Furthermore, we find that while most of the universal parameters remain unchanged, some universal amplitudes, like the excess cluster number, can be modified and become non-universal. We confirm these predictions by extensive Monte Carlo simulations of the two-dimensional percolation problem which has y can = −1/2. This article is part of 'Lattice models and integrability', a special issue of Journal of Physics A: Mathematical and Theoretical in honour of F Y Wu's 80th birthday.
We simulate the bond and site percolation models on a simple-cubic lattice with linear sizes up to L=512L=512L=512, and estimate the percolation thresholds to be pc(bond)=0.24881182(10){p}_{c}(\mathrm{bond})=0.248\phantom{\rule{0.16em}{0ex}}811\phantom{\rule{0.16em}{0ex}}82(10)pc​(bond)=0.24881182(10) and pc(site)=0.3116077(2){p}_{c}(\mathrm{site})=0.311\phantom{\rule{0.16em}{0ex}}607\phantom{\rule{0.16em}{0ex}}7(2)pc​(site)=0.3116077(2). By performing extensive simulations at these estimated critical points, we then estimate the critical exponents $1/\ensuremath{\nu}=1.1410(15)$, $\ensuremath{\beta}/\ensuremath{\nu}=0.477\phantom{\rule{0.16em}{0ex}}05(15)$, the leading correction exponent ${y}_{i}=\ensuremath{-}1.2(2)$, and the shortest-path exponent dmin=1.3756(3){d}_{\mathrm{min}}=1.3756(3)dmin​=1.3756(3). Various universal amplitudes are also obtained, including wrapping probabilities, ratios associated with the cluster-size distribution, and the excess cluster number. We observe that the leading finite-size corrections in certain wrapping probabilities are governed by an exponent $\ensuremath{\approx}$$\ensuremath{-}2$, rather than ${y}_{i}\ensuremath{\approx}\ensuremath{-}1.2$.
The fractal dimensions and the percolation exponents of geometrical spin clusters of like sign at criticality are obtained numerically for an Ising model with temperature dependent annealed bond dilution, also known as the thermalized bond Ising model (TBIM), in two dimensions. For this purpose, a modified Wolff single-cluster Monte Carlo simulation is used to generate equilibrium spin configurations on square lattices in the critical region. A tie-breaking rule is employed to identify nonintersecting spin cluster boundaries along the edges of the dual lattice. The values obtained for the fractal dimensions of the spanning geometrical clusters D c , and their interfaces D I , are in perfect agreement with those reported for the standard two-dimensional ferromagnetic Ising model. Furthermore, the variance of the winding angles results in a diffusivity κ = 3 for the two-dimensional thermalized bond Ising model, thus placing it in the universality class of the regular Ising model. A finite-size scaling analysis of the largest geometrical clusters results in a reliable estimation of the critical percolation exponents for the geometrical clusters in the limit of an infinite lattice size. The percolation exponents thus obtained are also found to be consistent with those reported for the regular Ising model. These consistencies are explained in terms of the Fisher renormalization relations, which express the thermodynamic critical exponents of systems with annealed bond dilution in terms of those of the regular model system.
In two space dimensions, the percolation point of the pure-site clusters of the Ising model coincides with the critical point Tc{T}_{c}Tc​ of the thermal transition and the percolation exponents belong to a special universality class. By introducing a bond probability pB<1,{p}_{B}<1,pB​<1, the corresponding site-bond clusters keep on percolating at Tc{T}_{c}Tc​ and the exponents do not change, until ${p}_{B}{=p}_{\mathrm{CK}}=1\ensuremath{-}\mathrm{exp}(\ensuremath{-}2J/kT):$ for this special expression of the bond weight the critical percolation exponents switch to the two-dimensional Ising universality class. We show here that the result is valid for a wide class of bidimensional models with a continuous magnetization transition: there is a critical bond probability pc{p}_{c}pc​ such that, for any pB> pc,{p}_{B}>~{p}_{c},pB​> pc​, the onset of percolation of the site-bond clusters coincides with the critical point of the thermal transition. The percolation exponents are the same for pc<pB< 1{p}_{c}<{p}_{B}<~1pc​<pB​< 1 but, for pB=pc,{p}_{B}{=p}_{c},pB​=pc​, they suddenly change to the thermal exponents, so that the corresponding clusters are critical droplets of the phase transition. Our result is based on Monte Carlo simulations of various systems near criticality.
Using Monte Carlo methods and finite-size scaling, we investigate surface critical phenomena in the bond-percolation model on the simple-cubic lattice with two open surfaces in one direction. We decompose the whole lattice into percolation clusters and sample the surface and bulk dimensionless ratios Q1 and Qb, defined on the basis of the moments of the cluster-size distribution. These ratios are used to determine critical points. At the bulk percolation threshold pbc, we determine the surface bond-occupation probability at the special transition as p(s)1c = 0.418 17(2), and further obtain the corresponding surface thermal and magnetic exponents as y(s)t1 = 0.5387(2) and y(s)h1 = 1.8014(6), respectively. Next, from the pair correlation function on the surfaces, we determine y(o)h1 = 1.0246(4) and y(e)h1 = 1.25(6) for the ordinary and the extraordinary transition, respectively, of which the former is consistent with the existing result y(o)h1 = 1.024(4). We also numerically derive the line of surface phase transitions occurring at pb < pbc, and determine the pertinent asymptotic values of the universal ratios Q1 and Qb.
Field-theoretical calculations predict that, at the upper critical dimension dc=4{d}_{c}=4dc​=4, the finite-size scaling (FSS) behaviors of the Ising model would be modified by multiplicative logarithmic corrections with thermal and magnetic correction exponents $({\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{y}}_{t},{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{y}}_{h})=(1/6,1/4)$. Using high-efficient cluster algorithms and the lifted worm algorithm, we present a systematic study to the FSS of the four-dimensional Ising model at criticality in the Fortuin-Kasteleyn (FK) bond and loop representations. In the FK representation, the size of the largest cluster is observed to scale as ${C}_{1}\ensuremath{\sim}{L}^{3}{(lnL)}^{{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{y}}_{h}}$, while the size of the second-largest cluster scales as ${C}_{2}\ensuremath{\sim}{L}^{3}{(lnL)}^{{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{y}}_{h2}}$ with ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{y}}_{h2}=\ensuremath{-}1/4$ a new correction exponent not yet predicted from field theory. In the loop representation, we observe that the size of the largest loop cluster scales as ${F}_{1}\ensuremath{\sim}{L}^{2}{(lnL)}^{{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{y}}_{t}}$, and the specific heat scales as ${c}_{\mathrm{E}}\ensuremath{\sim}{(lnL)}^{2{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{y}}_{t}}$. This clarifies the long-standing open question that whether the specific heat for the critical Ising model at dc=4{d}_{c}=4dc​=4 diverges logarithmically.