Ferroelectrics

Monte Carlo Simulation of First Order Phase Transitions

Chunlei Wang, J. C. Li, M. L. Zhao, Manuel I. Marqués, C. Aragó, Julio A. Gonzalo

Published Oct 29, 2010• 10 citations

Abstract

Monte Carlo simulation with Metropolis algorithm has been performed on a cubic Ising lattice with four-spin couplings. Temperature dependence of the spin average and average energy has been obtained. We found that the phase transition is of first order when the strength of four-spin coupling is stronger a critical value, and this critical value is of lattice size dependence. Binder cumulants still can be used to determine the second order phase transition temperature. For first order phase transition, the spin average disappears abruptly at the temperature of the average energy being a third of the ground state energy.

DOI: 10.1080/00150191003670291

Related papers

Computer simulation of phase transitions in low-dimensional Potts models

Akai Kurbanovich Murtazaev, Albert Babaevich Babaev, М. А. Магомедов

The Monte Carlo method is used to study phase transitions in the two-dimensional Potts model with the number of states spin q=3 on square and hexagonal lattices. Considered systems with linear dimensions Lx L=N, L=21/102. The obtained numerical data indicate that in the considered In the Potts model, a second-order phase transition is observed in accordance with the analytical theory. The Binder cumulants method of the fourth order determines the values of the critical points of the Potts model on different lattices. Keywords: Phase transitions, Potts model, Monte Carlo, thermodynamic parameters, critical temperature.

Physics of the Solid StateJan 1, 20220 citations

Verification of Ising phase transitions in the three-dimensional Ashkin-Teller model using Monte Carlo simulations

G. Szukowski, G. Kamieniarz, Grzegorz Musiał

Monte Carlo simulations in the three-dimensional (3D) Ashkin-Teller model on a cubic lattice are performed in the regions of the two-parameter space diagram where Ising-type phase transitions are expected. The scaling behavior of the Binder cumulant Q and the magnetic susceptibility in the critical region are exploited. In simulations the periodic boundary conditions and the Metropolis algorithm are used. Starting from Ising critical exponents and applying the finite-size-scaling analysis of the cumulant Q with nonlinear corrections, the accurate positions of the critical couplings on the continuous phase transition lines are calculated. For these couplings the critical exponent y(h) is calculated analyzing the magnetic susceptibility in the framework of the finite-size scaling with corrections. The values of y(h) agree with the 3D Ising model value along two lines determined by the order parameter [ssigma] .

Physical Review EMar 21, 200812 citations

Monte-Carlo simulation of phase transition in 2D and 3D ising model

Rajendra Neupane, Hari Ram Krishna Gauli, Krishna Bahadur Rai +1 more

In this work, Markov Chain-Monte Carlo technique was used to study the phase transition in two and three dimensional Ising Model (IM) in a square and cubic lattice. The study of temperature dependence of average magnetization and specific heat in different magnetic fields has been carried out in the 3x3 and 3x3x3 lattice with periodic boundary. Critical temperature point kBTc / J for 2D and 3D Ising Model has been observed at around 2.2 and 4.3 respectively at zero field. Our work satisfies Onsager’s critical value in 2D IM. The simulation suggests bifurcation in average magnetization below critical temperature Tc. Temperature plays the role of increasing randomness of spins. We found that Ising Model in small lattice size still retains interesting features like spontaneous magnetization and symmetry breaking below Tc at B = 0. At a non-zero field, the likelihood of spins to prefer certain alignment depends on the direction of the external field and magnitude of magnetization depends on magnitude of field ±B. Specific heat Cv, which gives us fluctuation of energy at particular temperature, has been found to attain maximum value at critical regions.

Scientific WorldAug 4, 20237 citations

Phase transitions in the antiferromagnetic layered Ising model on a cubic lattice

Рамазанов Магомедшейх Курбанович, Akai Kurbanovich Murtazaev

Phase transitions of the layered antiferromagnetic Ising model on a cubic lattice with allowance for intralayer next nearest neighbor interactions in the range of 0.0 ⩽ r < 1.0 have been studied on the basis of the replica algorithm by the Monte Carlo method. The character of phase transitions has been analyzed on the basis of the histogram method and the Binder cumulant method. The phase diagram of the dependence of the critical temperature on the next nearest neighbor interaction has been constructed. It has been shown that a secondorder phase transition is observed in the interval of 0.0 ⩽ r ⩽ 0.5 and for the value r = 0.9 and a first-order phase transition is observed in the interval of 0.6 ⩽ r < 0.8.

Journal of Experimental and Theoretical Physics LettersApr 1, 201627 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, 20260 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 PhysicsJan 1, 20260 citations