Ferroelectrics

Phase Transitions in Three Dimensional Ising Systems under Free and Periodic Boundary Conditions

Ma Felisa Martínez Ruiz, Francisca Pajuelo, Jorge Garcı́a, C. Aragó, Julio A. Gonzalo

Published Jan 1, 2004• 0 citations

Abstract

The Ising model (and also the Ising model in a transverse field) has been used extensively to investigate ferroelectric phase transitions. Numerical Monte Carlo simulations of phase transitions in three-dimensional Ising systems with variable size (2 ≤ L ≤ 100) have been performed. The susceptibility maximum and peak temperature Tc 3D(L,Lc) for free and periodic boundary conditions have been determinated and have been observed to be well fitted by simple expressions as a function of L. For L → ∞, as expected [Tc 3D]fbc = [Tc 3D]pbc being fbc = free boundary conditions and pbc = periodic boundary conditions. Scaling plots MLβ/v vs |ϵ|Ll/v for larger L with both boundary conditions were obteined and compared. The scaling function for our data up to L ≈ 100 was examined closely and an excellent fit at T → Tc was obtenied.

DOI: 10.1080/00150190490455197

Related papers

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

Effects of Shape and Boundary Conditions on Finite-Size Scaling Functions for Anisotropic Three-Dimensional Ising Systems

Kazuhisa Kaneda, Yutaka Okabe, Macoto Kikuchi

By use of Monte Carlo simulations, we study the finite-size scaling (FSS) functions for anisotropic three-dimensional (3D) Ising models of size L1 ×L1 ×aL1 (a: anisotropy parameter). Paying attention to the a dependence, we have shown that the FSS functions for the magnetization distribution function p(m) at the critical temperature change from a two-peak structure to a single-peak one by increasing or decreasing a from 1. We also study the effects of boundary conditions on the FSS functions for the anisotropic 3D Ising models.

Progress of Theoretical Physics SupplementJan 1, 20001 citation

Critical behavior of a cubic-lattice 3D Ising model for systems with quenched disorder

Akai Kurbanovich Murtazaev, Ibragimkhan K. Kamilov, Albert Babaevich Babaev

A Monte Carlo method is applied to simulate the static critical behavior of a cubic-lattice 3D Ising model for systems with quenched disorder. Numerical results are presented for the spin concentrations of p = 1.0, 0.95, 0.9, 0.8, 0.6 on L × L × L lattices with L = 20–60 under periodic boundary conditions. The critical temperature is determined by the Binder cumulant method. A finite-size scaling technique is used to calculate the static critical exponents α, β, γ, and ν (for specific heat, susceptibility, magnetization, and correlation length, respectively) in the range of p under study. Universality classes of critical behavior are discussed for three-dimensional diluted systems.

Journal of Experimental and Theoretical PhysicsDec 1, 200435 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

Thermodynamic characteristics of three-dimensional Ising-like systems as functions of microscopic parameters. The ρ6 model approximation

Ihor V. Pylyuk

The microscopic approach to calculating and investigating the three-dimensional (3D) Ising-like system thermodynamics above and below the phase transition temperature Tc is developed using the collective variables method.Mathematical description is performed within the framework of the higher non-Gaussian approximation ( 6 model) taking into account the corrections to scaling.A nonuniversal factor depending on microscopic parameters of the system is singled out in the expressions for leading critical amplitudes and correction-to-scaling amplitudes.Numerical estimates of the critical region size, phase transition temperature as well as the graphs of temperature dependences of entropy and other thermodynamic characteristics near Tc are given for various values of eective radius of the exponentially decreasing interaction potential.The variation of the free energy of the system at the phase transition point, average spin moment, and speci c heat with increasing ratio of the potential eective radius to the simple cubic lattice constant is traced.The results of calculations and their comparison with the other authors' data show that the 6 model provides a more adequate quantitative description of the critical behaviour of a 3D Ising-like system than the 4 model.

Journal of Physical StudiesJan 1, 20010 citations

On ergodic relaxation time in the three-dimensional Ising model

Robertas Grigalaits, Saulius Lapinskas, J. Banys +1 more

We have studied the dynamical decay of the autocorrelation function of the 3D Ising model for different sizes L = 20–52 of spin cluster-cubes. The behaviour of the longest, ergodic relaxation time, τe, of a finite domain below the phase transition temperature Tc was mostly considered for two types of phase transition dynamics. A study of the scaling properties of τe demonstrates a negligible difference between the types of dynamics used, but a considerable difference for different boundary conditions. In contrast to the known result for periodic boundary conditions (τe ~ Lz exp [const(Lєν)2], where z and ν are the dynamical and correlation length exponents, respectively, and є = 1 – T/Tc), the ergodic relaxation time for open boundary conditions is proportional to Lz exp [const(Lєν)2k] with coeffcient k for lattices explored in this work slightly decreasing with L in between 1.65 and 1.58. This result implies that only the lattices of sizes close to or exceeding L = 300 with open boundary conditions might have ergodic relaxation times similar to those with perodic boundary conditions.

Lithuanian Journal of PhysicsJan 1, 20131 citation