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:2410.14360 · Physical review. E

Spin- S Ising models with multispin interactions on the one-dimensional chain and two-dimensional square lattice

Kohei Suzuki

Published Feb 19, 2025• 2 citations
arXivHTMLPDF

Abstract

We study spin-S Ising models with p-spin interactions on the one-dimensional chain and the two-dimensional square lattice. Here, S denotes the magnitude of the spin and p represents the number of spins involved in each interaction. The analysis is performed for S=1/2,1,3/2,2, and p=3,4,5. For the one-dimensional model, we formulate transfer matrices and numerically diagonalize them to analyze the temperature dependence of the free energy and spin-spin correlations. In the case of S=1/2, the free energy does not depend on p, and the spin-spin correlations are uniformly enhanced across all temperature scales as p increases. In contrast, for S≥1, the free energy varies with p, and the spin-spin correlations are significantly enhanced at lower temperatures as p increases. For the two-dimensional model, by using multicanonical simulations, we analyze physical quantities such as an order parameter, internal energy, and specific heat. In addition, we define and examine an order parameter to distinguish ordered and disordered phases. It is found that a first-order phase transition occurs at finite temperatures for all S and p≥3, and increasing p strengthens its nature. We present S and p dependence of the transition temperature and latent heat, and discuss effects of higher-order interactions on the nature of phase transitions.

DOI
10.1103/physreve.111.024132

Related papers

Investigation of Finite-Size 2D Ising Model with a Noisy Matrix of Spin-Spin Interactions

Boris Kryzhanovsky, Magomed Malsagov, Iakov Karandashev

We analyze changes in the thermodynamic properties of a spin system when it passes from the classical two-dimensional Ising model to the spin glass model, where spin-spin interactions are random in their values and signs. Formally, the transition reduces to a gradual change in the amplitude of the multiplicative noise (distributed uniformly with a mean equal to one) superimposed over the initial Ising matrix of interacting spins. Considering the noise, we obtain analytical expressions that are valid for lattices of finite sizes. We compare our results with the results of computer simulations performed for square N = L × L lattices with linear dimensions L = 50 ÷ 1000. We find experimentally the dependencies of the critical values (the critical temperature, the internal energy, entropy and the specific heat) as well as the dependencies of the energy of the ground state and its magnetization on the amplitude of the noise. We show that when the variance of the noise reaches one, there is a jump of the ground state from the fully correlated state to an uncorrelated state and its magnetization jumps from 1 to 0. In the same time, a phase transition that is present at a lower level of the noise disappears.

arXiv:2002.00619•Entropy•Aug 7, 2018•7 citations

Crossover anomalies in frustrated two-dimensional Ising spin system with four-spin interaction

J. Szṅajd

As is well known, one-dimensional spin systems with short-range interactions cannot undergo a phase transition at finite temperature. However, in many 1D Ising spin models with frustration, at a certain crossover point, the thermodynamic quantities show anomalies that resemble such a transition W. Yin [Phys. Rev. Res. 6, 013331 (2024)2643-156410.1103/PhysRevResearch.6.013331]. Recently, Yin [Phys. Rev. Res. 6, 013331 (2024)2643-156410.1103/PhysRevResearch.6.013331] has shown, using decorated ladders as an example, that a one-dimensional Ising model can be arbitrarily close to the genuine phase transition. Using the transfer matrix method and real space renormalization group technique, we find that crossover resembling a phase transition accompanied by gigantic specific heat can also occur in two-dimensional systems. We consider coupled Ising spin chains with a zig-zag structure and a four-spin j_{4} interaction, which allows us to control the temperature at which the spins become unbound in the model with short-range interactions. The crossover point is defined by two parameters: the crossover temperature t=t^{*} and the crossover interaction j_{4}=j_{4}^{*}, and occurs when the temperature at which the effective interaction between the spin subsystems zeros out coincides with the temperature at which the spins become unbound. In systems in which a genuine phase transition occurs at the critical temperature t=t_{c}, there are two possibilities: if t_{c}>t^{*} the system undergoes a sequence of three phase transitions with a reentrant transition to disordered phase, while for t_{c}<t^{*} the phase transition is preceded by the considered crossover. We study the key role in the appearance of anomalies in the thermodynamic quantities of the four-spin interactions which are suspected to have a significant impact on the properties of magnetic two-dimensional lattices, but possible also on the construction of social systems modeling the phenomenon of social validation or peer pressure.

Physical review. E•Feb 13, 2025•4 citations

Dependence of Critical Parameters of 2D Ising Model on Lattice Size

B. V. Kryzhanovsky, M. Yu. Malsagov, Iakov Karandashev

For the 2D Ising model, we analyzed dependences of thermodynamic characteristics on number of spins by means of computer simulations. We compared experimental data obtained using the Fisher-Kasteleyn algorithm on a square lattice with N = l × l spins and the asymptotic Onsager solution ( N → ∞). We derived empirical expressions for critical parameters as functions of N and generalized the Onsager solution on the case of a finite-size lattice. Our analytical expressions for the free energy and its derivatives (the internal energy, the energy dispersion and the heat capacity) describe accurately the results of computer simulations. We showed that when N increased the heat capacity in the critical point increased as ln N . We specified restrictions on the accuracy of the critical temperature due to finite size of our system. Also in the finite-dimensional case, we obtained expressions describing temperature dependences of the magnetization and the correlation length. They are in a good qualitative agreement with the results of computer simulations by means of the dynamic Metropolis Monte Carlo method.

arXiv:2002.00615•Optical Memory and Neural Networks•Jan 1, 2018•7 citations

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

Monte Carlo Investigation of Magnetic Phase Transitions in the Two-Dimensional Ising Model

Satyendra Lohar, Neha Srivas, Akhilesh Dwivedi

Many problems in statistical physics resist exact solution once the number of interacting particles becomes large, and the two-dimensional Ising model is a classic case where numerical simulation fills that gap. This paper reports a Monte Carlo study of the two-dimensional Ising model using the Metropolis algorithm, undertaken to examine how magnetization, internal energy, specific heat, and magnetic susceptibility change with temperature. Spin configurations were generated on a square lattice under periodic boundary conditions, with each proposed spin flip accepted or rejected according to the Boltzmann weight of the resulting energy change. The simulation reproduces the expected transition from an ordered ferromagnetic state at low temperature to a disordered paramagnetic state at high temperature, with magnetization falling sharply and both specific heat and susceptibility peaking near a critical temperature of Tc ≈ 2.3–2.4 (in units of J/kB), consistent with the exact value obtained analytically for this system. These results confirm that the Metropolis algorithm, despite its simplicity, remains a dependable tool for locating second-order phase transitions and extracting thermodynamic quantities in lattice spin systems, and the same computational framework can be adapted to more elaborate models used in condensed matter and materials research.

•Aug 14, 2026•0 citations

Spin correlations in low-dimensional spin-1/2 antiferromagnets

Biao Hao, Chang-De Gong

We study the long-range behavior of spin correlations in low-dimensional spin-1/2 antiferromagnets with both Ising anisotropy and spatial anisotropy. We represent the spin operators in terms of spin-density operators and spin-phase operators. Then for spin-1/2 antiferromagnets, within a simple harmonic approximation, we obtain expressions for the spin correlations in the ground state. For the one-dimensional linear chain, we find that the two leading terms in spin correlation are (-1${)}^{\mathit{r}}$c/${\mathit{r}}^{\mathrm{\ensuremath{\eta}}}$-${\mathit{c}}_{1}$/${\mathit{r}}^{2}$ for large distance r. The power \ensuremath{\eta} depends on the Ising anisotropy, and \ensuremath{\eta}=0.405 at XY point and \ensuremath{\eta}=1.062 at Heisenberg point. For the two-dimensional square lattice, there is long-range staggered magnetization, and the dependence of magnetization on Ising anisotropy is given. For the crossover from one-dimensional linear chain to two-dimensional square lattice, we find that the disorder-order transition occurs as long as the interchain coupling is introduced. We believe that our method gives a unified approach to the long-range behavior of spin correlations in low-dimensional spin-1/2 antiferromagnets.

Physical review. B, Condensed matter•Jul 1, 1995•5 citations
We study spin-S Ising models with p-spin interactions on the one-dimensional chain and the two-dimensional square lattice. Here, S denotes the magnitude of the spin and p represents the number of spins involved in each interaction. The analysis is performed for S=1/2,1,3/2,2, and p=3,4,5. For the one-dimensional model, we formulate transfer matrices and numerically diagonalize them to analyze the temperature dependence of the free energy and spin-spin correlations. In the case of S=1/2, the free energy does not depend on p, and the spin-spin correlations are uniformly enhanced across all temperature scales as p increases. In contrast, for S≥1, the free energy varies with p, and the spin-spin correlations are significantly enhanced at lower temperatures as p increases. For the two-dimensional model, by using multicanonical simulations, we analyze physical quantities such as an order parameter, internal energy, and specific heat. In addition, we define and examine an order parameter to distinguish ordered and disordered phases. It is found that a first-order phase transition occurs at finite temperatures for all S and p≥3, and increasing p strengthens its nature. We present S and p dependence of the transition temperature and latent heat, and discuss effects of higher-order interactions on the nature of phase transitions.
We analyze changes in the thermodynamic properties of a spin system when it passes from the classical two-dimensional Ising model to the spin glass model, where spin-spin interactions are random in their values and signs. Formally, the transition reduces to a gradual change in the amplitude of the multiplicative noise (distributed uniformly with a mean equal to one) superimposed over the initial Ising matrix of interacting spins. Considering the noise, we obtain analytical expressions that are valid for lattices of finite sizes. We compare our results with the results of computer simulations performed for square N = L × L lattices with linear dimensions L = 50 ÷ 1000. We find experimentally the dependencies of the critical values (the critical temperature, the internal energy, entropy and the specific heat) as well as the dependencies of the energy of the ground state and its magnetization on the amplitude of the noise. We show that when the variance of the noise reaches one, there is a jump of the ground state from the fully correlated state to an uncorrelated state and its magnetization jumps from 1 to 0. In the same time, a phase transition that is present at a lower level of the noise disappears.
As is well known, one-dimensional spin systems with short-range interactions cannot undergo a phase transition at finite temperature. However, in many 1D Ising spin models with frustration, at a certain crossover point, the thermodynamic quantities show anomalies that resemble such a transition W. Yin [Phys. Rev. Res. 6, 013331 (2024)2643-156410.1103/PhysRevResearch.6.013331]. Recently, Yin [Phys. Rev. Res. 6, 013331 (2024)2643-156410.1103/PhysRevResearch.6.013331] has shown, using decorated ladders as an example, that a one-dimensional Ising model can be arbitrarily close to the genuine phase transition. Using the transfer matrix method and real space renormalization group technique, we find that crossover resembling a phase transition accompanied by gigantic specific heat can also occur in two-dimensional systems. We consider coupled Ising spin chains with a zig-zag structure and a four-spin j_{4} interaction, which allows us to control the temperature at which the spins become unbound in the model with short-range interactions. The crossover point is defined by two parameters: the crossover temperature t=t^{*} and the crossover interaction j_{4}=j_{4}^{*}, and occurs when the temperature at which the effective interaction between the spin subsystems zeros out coincides with the temperature at which the spins become unbound. In systems in which a genuine phase transition occurs at the critical temperature t=t_{c}, there are two possibilities: if t_{c}>t^{*} the system undergoes a sequence of three phase transitions with a reentrant transition to disordered phase, while for t_{c}<t^{*} the phase transition is preceded by the considered crossover. We study the key role in the appearance of anomalies in the thermodynamic quantities of the four-spin interactions which are suspected to have a significant impact on the properties of magnetic two-dimensional lattices, but possible also on the construction of social systems modeling the phenomenon of social validation or peer pressure.
For the 2D Ising model, we analyzed dependences of thermodynamic characteristics on number of spins by means of computer simulations. We compared experimental data obtained using the Fisher-Kasteleyn algorithm on a square lattice with N = l × l spins and the asymptotic Onsager solution ( N → ∞). We derived empirical expressions for critical parameters as functions of N and generalized the Onsager solution on the case of a finite-size lattice. Our analytical expressions for the free energy and its derivatives (the internal energy, the energy dispersion and the heat capacity) describe accurately the results of computer simulations. We showed that when N increased the heat capacity in the critical point increased as ln N . We specified restrictions on the accuracy of the critical temperature due to finite size of our system. Also in the finite-dimensional case, we obtained expressions describing temperature dependences of the magnetization and the correlation length. They are in a good qualitative agreement with the results of computer simulations by means of the dynamic Metropolis 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.
Many problems in statistical physics resist exact solution once the number of interacting particles becomes large, and the two-dimensional Ising model is a classic case where numerical simulation fills that gap. This paper reports a Monte Carlo study of the two-dimensional Ising model using the Metropolis algorithm, undertaken to examine how magnetization, internal energy, specific heat, and magnetic susceptibility change with temperature. Spin configurations were generated on a square lattice under periodic boundary conditions, with each proposed spin flip accepted or rejected according to the Boltzmann weight of the resulting energy change. The simulation reproduces the expected transition from an ordered ferromagnetic state at low temperature to a disordered paramagnetic state at high temperature, with magnetization falling sharply and both specific heat and susceptibility peaking near a critical temperature of Tc ≈ 2.3–2.4 (in units of J/kB), consistent with the exact value obtained analytically for this system. These results confirm that the Metropolis algorithm, despite its simplicity, remains a dependable tool for locating second-order phase transitions and extracting thermodynamic quantities in lattice spin systems, and the same computational framework can be adapted to more elaborate models used in condensed matter and materials research.
We study the long-range behavior of spin correlations in low-dimensional spin-1/2 antiferromagnets with both Ising anisotropy and spatial anisotropy. We represent the spin operators in terms of spin-density operators and spin-phase operators. Then for spin-1/2 antiferromagnets, within a simple harmonic approximation, we obtain expressions for the spin correlations in the ground state. For the one-dimensional linear chain, we find that the two leading terms in spin correlation are (-1)r{)}^{\mathit{r}})rc/${\mathit{r}}^{\mathrm{\ensuremath{\eta}}}$-c1{\mathit{c}}_{1}c1​/r2{\mathit{r}}^{2}r2 for large distance r. The power \ensuremath{\eta} depends on the Ising anisotropy, and \ensuremath{\eta}=0.405 at XY point and \ensuremath{\eta}=1.062 at Heisenberg point. For the two-dimensional square lattice, there is long-range staggered magnetization, and the dependence of magnetization on Ising anisotropy is given. For the crossover from one-dimensional linear chain to two-dimensional square lattice, we find that the disorder-order transition occurs as long as the interchain coupling is introduced. We believe that our method gives a unified approach to the long-range behavior of spin correlations in low-dimensional spin-1/2 antiferromagnets.