We investigate the dependence of the critical Binder cumulant of the magnetization and the largest Fortuin-Kasteleyn cluster on the boundary conditions and aspect ratio of the underlying square Ising lattices. By means of the Swendsen-Wang algorithm, we generate numerical data for large system sizes and we perform a detailed finite-size scaling analysis for several values of the aspect ratio r, for both periodic and free boundary conditions. We estimate the universal probability density functions of the largest Fortuin-Kasteleyn cluster and we compare it to those of the magnetization at criticality. It is shown that these probability density functions follow similar scaling laws, and it is found that the values of the critical Binder cumulant of the largest Fortuin-Kasteleyn cluster are upper bounds to the values of the respective order-parameter's cumulant, with a splitting behavior for large values of the aspect ratio. We also investigate the dependence of the amplitudes of the magnetization and the largest Fortuin-Kasteleyn cluster on the aspect ratio and boundary conditions. We find that the associated exponents, describing the aspect-ratio dependencies, are different for the magnetization and the largest Fortuin-Kasteleyn cluster, but in each case are independent of boundary conditions.
We present a Monte Carlo study of the geometric properties of Fortuin-Kasteleyn (FK) clusters of the Ising model on square [two-dimensional (2D)] and simple-cubic [three-dimensional (3D)] lattices. The wrapping probability, a dimensionless quantity characterizing the topology of the FK clusters on a torus, is found to suffer from smaller finite-size corrections than the well-known Binder ratio and yields a high-precision critical coupling as
Kc(3D)=0.221654631(8). We then study other geometric properties of FK clusters at criticality. It is demonstrated that the distribution of the critical largest-cluster size
C1 follows a single-variable function as $P({C}_{1},L)\phantom{\rule{0.16em}{0ex}}d{C}_{1}=\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{P}(x)\phantom{\rule{0.16em}{0ex}}dx$ with $x\ensuremath{\equiv}{C}_{1}/{L}^{{d}_{\mathrm{F}}}$ (
L is the linear size), where the fractal dimension
dF is identical to the magnetic exponent. An interesting bimodal feature is observed in distribution $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{P}(x)$ in three dimensions, and attributed to the different approaching behaviors for $K\ensuremath{\rightarrow}{K}_{c}+{0}^{\ifmmode\pm\else\textpm\fi{}}$. To characterize the compactness of the FK clusters, we measure their graph distances and determine the shortest-path exponents as
dmin(3D)=1.25936(12) and
dmin(2D)=1.0940(2). Further, by excluding all the bridges from the occupied bonds, we obtain bridge-free configurations and determine the backbone exponents as
dB(3D)=2.1673(15) and
dB(2D)=1.7321(4). The estimates of the universal wrapping probabilities for the 3D Ising model and of the geometric critical exponents
dmin and
dB either improve over the existing results or have not been reported yet. We believe that these numerical results would provide a testing ground in the development of further theoretical treatments of the 3D Ising model.
Large-deviation tails of critical probability distributions provide a sensitive probe of universality beyond standard finite-size scaling. We study these tails for critical percolation and Fortuin--Kasteleyn Ising models on two-dimensional lattices, three-dimensional lattices, and complete graphs. We consider two rescaled order parameters: the magnetization-like variable
xm=∣M∣/⟨∣M∣⟩, including a signed cluster-mass analogue for percolation, and the largest-cluster variable
xC=C1/⟨C1⟩. For
xm, we test the expected stretched-exponential large-deviation tail and show that the same form applies to the percolation analogue. For
xC, guided by the exact complete-graph result and scaling arguments, we propose universal scaling forms for both tails of the cumulative distribution and test them by extensive Monte Carlo simulations. In the complete-graph FK-Ising model, the left tail is governed by rare configurations with percolation-like scaling rather than by the typical Ising scaling. Our results show that the tails of order-parameter distributions reveal universal features of critical fluctuations that are not captured by averaged observables alone.
The Fortuin-Kasteleyn (FK) random-cluster model, which can be exactly mapped from the q-state Potts spin model, is a correlated bond percolation model. By extensive Monte Carlo simulations, we study the FK bond representation of the critical Ising model (q=2) on a finite complete graph, i.e., the mean-field Ising model. We provide strong numerical evidence that the configuration space for q=2 contains an asymptotically vanishing sector in which quantities exhibit the same finite-size scaling as in the critical uncorrelated bond percolation (q=1) on the complete graph. Moreover, we observe that, in the full configuration space, the power-law behavior of the cluster-size distribution for the FK Ising clusters except the largest one is governed by a Fisher exponent taking the value for q=1 instead of q=2. This demonstrates the percolation effects in the FK Ising model on the complete graph.
We present an extensive Markov chain Monte Carlo study of the finite-size scaling behavior of the Fortuin-Kasteleyn Ising model on five-dimensional hypercubic lattices with periodic boundary conditions. We observe that physical quantities, which include the contribution of the largest cluster, exhibit complete graph asymptotics. However, for quantities where the contribution of the largest cluster is removed, we observe that the scaling behavior is mainly controlled by the Gaussian fixed point. Our results therefore suggest that both scaling predictions, i.e., the complete graph and the Gaussian fixed point asymptotics, are needed to provide a complete description for the five-dimensional finite-size scaling behavior on the torus.
The five-dimensional Ising model with free boundary conditions has recently received a renewed interest in a debate concerning the finite-size scaling of the susceptibility near the critical temperature. We provide evidence in favour of the conventional scaling picture, where the susceptibility scales as L2 inside a critical scaling window of width 1/L2. Our results are based on Monte Carlo data gathered on system sizes up to L=79 (ca. three billion spins) for a wide range of temperatures near the critical point. We analyse the magnetisation distribution, the susceptibility and also the scaling and distribution of the size of the Fortuin–Kasteleyn cluster containing the origin. The probability of this cluster reaching the boundary determines the correlation length, and its behaviour agrees with the mean field critical exponent δ=3, that the scaling window has width 1/L2.
Using Monte Carlo techniques, the critical Binder cumulant U * of isotropic nearest neighbour Ising models on square and triangular lattices is studied. For rectangular shapes, employing periodic boundary conditions, U * is found to show the same dependence on the aspect ratio for both lattice types. Similarly, applying free boundary conditions for systems with square as well as circular shapes for both lattices, the simulational findings are also consistent with the suggestion that, for isotropic Ising models with short range interactions, U * depends on the shape and the boundary condition, but not on the lattice structure.