In this paper we study finite-size scaling at criticality in the Gibbs ensemble with Monte Carlo simulations. A recent lattice-gas version of Gibbs-ensemble sampling for the two-dimensional Ising-model magnet with nearest-neighbor ferromagnetic interactions on two subsystems of L\ifmmode\times\else\texttimes\fi{}L lattices is considered. Both square and triangular lattices with a range of lattice sizes are simulated at bulk criticality. We are interested in the question of universality for the fourth-order cumulant of the order-parameter distribution in the Gibbs ensemble. Estimates of -1.954\ifmmode\pm\else\textpm\fi{}0.005 and -1.9543\ifmmode\pm\else\textpm\fi{}0.005 are obtained for the square and triangular lattices, respectively. These values confirm universality within the Gibbs ensemble, but are statistically different from the canonical-ensemble values of \ensuremath{\sim}-1.83. We point out that one must interpret this apparent ensemble dependence with care. Unlike in the canonical ensemble, the order parameter in the Gibbs ensemble is the subsystem magnetization. Its fourth-order cumulant is not simply related to the usual fourth magnetic-field derivative of the total free energy. Standard finite-size scaling analysis of data for the total specific heat, subsystem magnetization, and susceptibility are also considered and found to be in good agreement with finite-size scaling.
In this paper we study finite-size scaling at criticality in the Gibbs ensemble with Monte Carlo simulations. A recent lattice-gas version of Gibbs-ensemble sampling for the two-dimensional Ising-model magnet with nearest-neighbor ferromagnetic interactions on two subsystems of L\ifmmode\times\else\texttimes\fi{}L lattices is considered. Both square and triangular lattices with a range of lattice sizes are simulated at bulk criticality. We are interested in the question of universality for the fourth-order cumulant of the order-parameter distribution in the Gibbs ensemble. Estimates of -1.954\ifmmode\pm\else\textpm\fi{}0.005 and -1.9543\ifmmode\pm\else\textpm\fi{}0.005 are obtained for the square and triangular lattices, respectively. These values confirm universality within the Gibbs ensemble, but are statistically different from the canonical-ensemble values of \ensuremath{\sim}-1.83. We point out that one must interpret this apparent ensemble dependence with care. Unlike in the canonical ensemble, the order parameter in the Gibbs ensemble is the subsystem magnetization. Its fourth-order cumulant is not simply related to the usual fourth magnetic-field derivative of the total free energy. Standard finite-size scaling analysis of data for the total specific heat, subsystem magnetization, and susceptibility are also considered and found to be in good agreement with finite-size scaling.
We study the higher-order cumulant universality of the order-parameter distribution at criticality for the two-dimensional nearest-neighbor ferromagnetic Ising model with Monte Carlo methods. These cumulants are interesting because they are a quantitative measure of the shape of the critical order parameter distribution. This shape has been predicted to be universal. Up to now, only the fourth-order cumulant has been studied in any detail. The set of higher-order cumulants as a whole is more sensitive to the details of the distribution and hence is a stronger test. It cannot be sampled efficiently with standard Monte Carlo sampling. To the best of our knowledge, no definite numerical results have been reported in the literature. Using the umbrella-sampling technique, we are able to obtain accurate estimates for the fourth-, sixth-, eighth-, and tenth-order cumulants for both the square and triangular lattices of different lattice sizes. These results support universality.
To demonstrate the application of computational physics in modern statistical physics, this paper investigates the critical behavior of a finite-size ferromagnetic system based on the two-dimensional square-lattice Ising model.A nearest-neighbor Ising Hamiltonian is constructed, and periodic boundary conditions are adopted to reduce finite-boundary effects.The Metropolis Monte Carlo algorithm is used to sample spin configurations at different temperatures, while energy, magnetization, specific heat, magnetic susceptibility and Binder cumulant are calculated to describe the transition from an ordered state to a disordered state.The simulation results show that the system exhibits clear long-range order at low temperatures, where the order parameter remains high.As temperature increases, thermal fluctuations become stronger, magnetization decreases rapidly, and the system gradually evolves into a disordered state.The specific heat curve shows a pronounced peak at approximately (T = 2.269), while the susceptibility peak shifts slightly toward the high-temperature side due to finite-size effects and the discrete temperature grid.Overall, the simulated critical region is consistent with the theoretical critical temperature of the two-dimensional square-lattice Ising model.Spin-configuration analysis further reveals enhanced fluctuations and domain reconstruction near the critical region.In addition, a thermalization criterion and an extended sampling test for L = 48 near the critical region are added to examine the influence of finite sampling and equilibration.These results indicate that Monte Carlo simulation of the two-dimensional Ising model can effectively reproduce continuous phase transition, critical fluctuation and finite-size effects, providing a basic computational case for understanding critical phenomena in complex systems and data-driven physical analysis.
On Critical Temperature and Finite Size Scaling of Continuous Spin 2d Ising Model
In this paper, we have studied the critical temperature Tc of continuous spin 2d square-lattice Ising model using Monte-Carlo simulation. We have considered spins s in a bounded interval, where s∈[−1,+1] in square-lattice configuration with periodic boundary condition. We have observed that the critical temperature Tc is approximately 0.925, showing a clear second order phase transition. Considering finite size scaling, we have also obtained the critical exponents associated with susceptibility, specific heat, magnetization and we find that these values are in good agreement with the corresponding values obtained for the standard 2d Ising universality class.
The corrections to finite-size scaling in the critical two-point correlation function G(r) of 2D Ising model on a square lattice have been studied numerically by means of exact transfer-matrix algorithms. The systems have been considered, including up to 800 spins. The calculation of G(r) at a distance r equal to the half of the system size L shows the existence of an amplitude correction proportional to 1/L^2. A nontrivial correction ~1/L^0.25 of a very small magnitude also has been detected, as it can be expected from our recently developed GFD (grouping of Feynman diagrams) theory. Monte Carlo simulations of the squared magnetization of 3D Ising model have been performed by Wolff's algorithm in the range of the reduced temperatures t =< 0.000086 and system sizes L =< 410. The effective critical exponent beta_eff(t) tends to increase above the currently accepted numerical values. The critical coupling K_c=0.22165386(51) has been extracted from the Binder cumulant data within 96 =< L =< 384. The critical exponent 1/nu, estimated from the finite-size scaling of the derivatives of the Binder cumulant, tends to decrease slightly below the RG value 1.587 for the largest system sizes. The finite-size scaling of accurately simulated maximal values of the specific heat C_v in 3D Ising model confirms a logarithmic rather than power-like critical singularity of C_v.
In this work we present a thorough analysis of the phase transitions that occur in a ferromagnetic two-dimensional Ising model, with only nearest-neighbors interactions, in the framework of the Tsallis nonextensive statistics. We performed Monte Carlo simulations on square lattices with linear sizes L ranging from 32 up to 512. The statistical weight of the Metropolis algorithm was changed according to the nonextensive statistics. Discontinuities in the m(T) curve are observed for $q\ensuremath{\le}0.5$. However, we have verified only one peak on the energy histograms at the critical temperatures, indicating the occurrence of continuous phase transitions. For the $0.5<q\ensuremath{\le}1.0$ regime, we have found continuous phase transitions between the ordered and the disordered phases, and determined the critical exponents via finite-size scaling. We verified that the critical exponents $\ensuremath{\alpha}$, $\ensuremath{\beta}$, and $\ensuremath{\gamma}$ depend on the entropic index q in the range $0.5<q\ensuremath{\le}1.0$ in the form $\ensuremath{\alpha}(q)=(10{q}^{2}\ensuremath{-}33q+23)/20$, $\ensuremath{\beta}(q)=(2q\ensuremath{-}1)/8$, and $\ensuremath{\gamma}(q)=({q}^{2}\ensuremath{-}q+7)/4$. On the other hand, the critical exponent $\ensuremath{\nu}$ does not depend on q. This suggests a violation of the scaling relations $2\ensuremath{\beta}+\ensuremath{\gamma}=d\ensuremath{\nu}$ and $\ensuremath{\alpha}+2\ensuremath{\beta}+\ensuremath{\gamma}=2$ and a nonuniversality of the critical exponents along the ferro-paramagnetic frontier.
Monte Carlo simulations are performed on the square lattice Ising system with both the nearest neighbor and the two-dimensional dipolar interactions in order to investigate the critical property of ferromagnetic phase transitions in this 2-dimensional system. Applying the finite size scaling theory to Monte Carlo data just near the transition point, the critical indices of this dipolar system are calculated and α=0.31±0.04, β =0.20±0.02, γ=1.08±0.20 and ν=0.72±0.10 are obtained. However, in the temperature region considerably far from the transition point, different values of α, β and γ seem to be adequate, which are rather close to those of the rigorous result of the 2-dimensional Ising model.