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

arXiv:0910.2218 · Physical Review E

Finite-size analysis of a two-dimensional Ising model within a nonextensive approach

Nuno Crokidakis, Diogo O. Soares-Pinto, Mario S. Reis, Alexandre M. Souza, Roberto S. Sarthour, Ivan S. Oliveira

Published Nov 2, 2009• 14 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physreve.80.051101

Related papers

Critical Behavior Simulation of the Two-Dimensional Ising Model—A Computational Physics Study Based on the Monte Carlo Method

佳钢 李

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.

Modern Physics•Jan 1, 2026•0 citations

Finite-size-scaling estimation of critical exponents in the two-dimensional axial next-nearest-neighbor Ising model

Marcelo D. Grynberg, Horacio Ceva

The magnetic critical exponent of the two-dimensional axial next-nearest-neighbor Ising model is estimated by means of the finite-size scaling of the order parameter. For a finite system without ordering field, this magnitude is computed from the two leading eigenvectors of the transfer matrix. The ferromagnetic-paramagnetic transition temperatures estimated with this method exhibit a small lattice-size dependence. Moreover, they are in good agreement with those obtained using phenomenological renormalization, and also with analytical approximations. The correlation-length critical exponent is calculated by means of phenomenological renormalization. The results found with these numerical techniques indicate that the ferromagnetic-paramagnetic transition of the model seems to be continuous, in contradiction with a previous result obtained using the persistence-length criterion, which indicated that the transition is always first order. We found that the phase transition is first order only when the ratio of the next-nearest-neighbor and nearest-neighbor interactions equals one-half.

Physical review. B, Condensed matter•Nov 1, 1988•6 citations

Phase transitions in the antiferromagnetic ising model on a square lattice with next-nearest-neighbor interactions

Akai Kurbanovich Murtazaev, Рамазанов Магомедшейх Курбанович, Felix A. Kassan-Ogly +1 more

The phase transitions in the two-dimensional Ising model on a square lattice are studied using a replica algorithm, the Monte Carlo method, and histogram analysis with allowance for the next-nearest-neighbor interactions in the range 0.1 ≤ r < 1.0. A phase diagram is constructed for the dependence of the critical temperature on the next-nearest-neighbor interaction. A second-order phase transition is detected in this range and the model under study.

Journal of Experimental and Theoretical Physics•Dec 1, 2013•28 citations

Monte Carlo studies of finite-size scaling in the Gibbs ensemble at criticality for the two-dimensional Ising model

K. K. Mon

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.

Physical review. B, Condensed matter•Mar 1, 1993•8 citations

Monte Carlo Studies of the Square Lattice Ising System with the Two-Dimensional Dipolar Interaction

Seizo Watarai, Keizô Yamamoto, Akira Nakanishi

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.

Journal of the Physical Society of Japan•May 15, 1992•0 citations

Critical scaling of the mutual information in two-dimensional disordered Ising models

Sriluckshmy, P., Mandal, I.

Renyi mutual information, computed from second Renyi entropies, can identify classical phase transitions from their finite-size scaling at critical points. We apply this technique to examine the presence or absence of finite temperature phase transitions in various two-dimensional models on a square lattice, which are extensions of the conventional Ising model by adding a quenched disorder. When the quenched disorder causes the nearest neighbor bonds to be both ferromagnetic and antiferromagnetic, (a) a spin glass phase exists only at zero temperature, and (b) a ferromagnetic phase exists at a finite temperature when the antiferromagnetic bond distributions are sufficiently dilute. Furthermore, finite temperature paramagnetic-ferromagnetic transitions can also occur when the disordered bonds involve only ferromagnetic couplings of random strengths. In our numerical simulations, the 'zero temperature only' phase transitions are identified when there is no consistent finite-size scaling of the Renyi mutual information curves, while for finite temperature critical points, the curves can identify the critical temperature T-c by their crossings at T-c and 2 T-c.

•Apr 3, 2018•8 citations
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≤0.5q\ensuremath{\le}0.5q\ensuremath≤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≤1.00.5<q\ensuremath{\le}1.00.5<q\ensuremath≤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α\ensuremath{\alpha}\ensuremathα, \ensuremathβ\ensuremath{\beta}\ensuremathβ, and \ensuremathγ\ensuremath{\gamma}\ensuremathγ depend on the entropic index qqq in the range 0.5<q\ensuremath≤1.00.5<q\ensuremath{\le}1.00.5<q\ensuremath≤1.0 in the form \ensuremathα(q)=(10q2\ensuremath−33q+23)/20\ensuremath{\alpha}(q)=(10{q}^{2}\ensuremath{-}33q+23)/20\ensuremathα(q)=(10q2\ensuremath−33q+23)/20, \ensuremathβ(q)=(2q\ensuremath−1)/8\ensuremath{\beta}(q)=(2q\ensuremath{-}1)/8\ensuremathβ(q)=(2q\ensuremath−1)/8, and \ensuremathγ(q)=(q2\ensuremath−q+7)/4\ensuremath{\gamma}(q)=({q}^{2}\ensuremath{-}q+7)/4\ensuremathγ(q)=(q2\ensuremath−q+7)/4. On the other hand, the critical exponent \ensuremathν\ensuremath{\nu}\ensuremathν does not depend on qqq. This suggests a violation of the scaling relations 2\ensuremathβ+\ensuremathγ=d\ensuremathν2\ensuremath{\beta}+\ensuremath{\gamma}=d\ensuremath{\nu}2\ensuremathβ+\ensuremathγ=d\ensuremathν and \ensuremathα+2\ensuremathβ+\ensuremathγ=2\ensuremath{\alpha}+2\ensuremath{\beta}+\ensuremath{\gamma}=2\ensuremathα+2\ensuremathβ+\ensuremathγ=2 and a nonuniversality of the critical exponents along the ferro-paramagnetic frontier.
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.
The magnetic critical exponent of the two-dimensional axial next-nearest-neighbor Ising model is estimated by means of the finite-size scaling of the order parameter. For a finite system without ordering field, this magnitude is computed from the two leading eigenvectors of the transfer matrix. The ferromagnetic-paramagnetic transition temperatures estimated with this method exhibit a small lattice-size dependence. Moreover, they are in good agreement with those obtained using phenomenological renormalization, and also with analytical approximations. The correlation-length critical exponent is calculated by means of phenomenological renormalization. The results found with these numerical techniques indicate that the ferromagnetic-paramagnetic transition of the model seems to be continuous, in contradiction with a previous result obtained using the persistence-length criterion, which indicated that the transition is always first order. We found that the phase transition is first order only when the ratio of the next-nearest-neighbor and nearest-neighbor interactions equals one-half.
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.
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.
Renyi mutual information, computed from second Renyi entropies, can identify classical phase transitions from their finite-size scaling at critical points. We apply this technique to examine the presence or absence of finite temperature phase transitions in various two-dimensional models on a square lattice, which are extensions of the conventional Ising model by adding a quenched disorder. When the quenched disorder causes the nearest neighbor bonds to be both ferromagnetic and antiferromagnetic, (a) a spin glass phase exists only at zero temperature, and (b) a ferromagnetic phase exists at a finite temperature when the antiferromagnetic bond distributions are sufficiently dilute. Furthermore, finite temperature paramagnetic-ferromagnetic transitions can also occur when the disordered bonds involve only ferromagnetic couplings of random strengths. In our numerical simulations, the 'zero temperature only' phase transitions are identified when there is no consistent finite-size scaling of the Renyi mutual information curves, while for finite temperature critical points, the curves can identify the critical temperature T-c by their crossings at T-c and 2 T-c.