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

Journal of Applied Physics

m=3 Ashkin–Teller-like cubic model on an FCC lattice

Ronald Fisch

Published Apr 15, 1996• 2 citations
Publisher

Abstract

Monte Carlo simulations have been used to study a model consisting of three Ising models with nearest neighbor exchange J on a face-centered-cubic lattice, which are coupled together by a constraint. The constraint requires that the vector sum of the four spins located at the corners of any elementary tetrahedron of the lattice not add to zero. Finite-size scaling analysis of the Monte Carlo results for L×L×L lattices with L=16, 32, and 64 has been used to find that Tc /J=10.84±0.01. The finite-size scaling functions for the magnetization and magnetic susceptibility near Tc are consistent with standard m=3 Heisenberg critical exponents, as predicted by renormalization group theory. Strengthening the four-spin constraint increases the value of the effective negative Ashkin–Teller four-spin coupling, and drives Tc to infinity.

DOI
10.1063/1.361310

Related papers

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 E•Mar 21, 2008•12 citations

Monte Carlo renormalization-group calculations of critical behavior in the simple-cubic Ising model

Godfrey Stuart Pawley, Robert H. Swendsen, D. J. Wallace +1 more

The Monte Carlo renormalization group is applied to the three-dimensional Ising model on simple cubic lattices with ${8}^{3}$, ${16}^{3}$, ${32}^{3}$, and ${64}^{3}$ sites. The comparison of block-spin correlation functions from the largest lattices yields the nearest-neighbor critical coupling ${K}_{1}^{c}=0.221654(6)$. After allowing for (i) interpolation to this best estimate for ${K}_{1}^{c}$, (ii) an apparent finite-size effect in the renormalization-group transformation due to the measurement of correlation functions of too few (seven) operators, and (iii) the extrapolation for the effect of a slow transient towards the fixed point, the values $\ensuremath{\nu}=0.629(4)$ and $\ensuremath{\eta}=0.031(5)$ are obtained for the thermal and magnetic exponents. The correction-to-scaling exponent $\ensuremath{\omega}$ is estimated to be around 1; to obtain an accuracy competitive with other methods requires measurements with more than seven operators. We briefly review the problem of redundant operators and indicate the future prospects for this kind of calculation.

Physical review. B, Condensed matter•Apr 1, 1984•318 citations

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

Grzegorz Musiał, L. Dębski, G. Kamieniarz

The results of the large-scale Monte Carlo simulations of the continuous phase transitions in the three-dimensional Ashkin-Teller model on a cubic lattice are presented in the two-parameter space. The invariance of the ratio of the square of the second moment of the order parameter to its fourth moment in the critical region and the Metropolis algorithm are exploited. Using the universality hypothesis and the finite-size-scaling analysis, the accurate positions of the continuous phase-transition lines are calculated and their Ising character is confirmed. Approximate locations of the tricritical points are also estimated.

Physical review. B, Condensed matter•Jun 28, 2002•27 citations

Worm-type Monte Carlo simulation of the Ashkin-Teller model on the triangular lattice

Jian-Ping Lv, Youjin Deng, Qing‐Hu Chen

We investigate the symmetric Ashkin-Teller (AT) model on the triangular lattice in the antiferromagnetic two-spin coupling region (J<0). In the J→-∞ limit, we map the AT model onto a fully packed loop-dimer model on the honeycomb lattice. On the basis of this exact transformation and the low-temperature expansion, we formulate a variant of worm-type algorithms for the AT model, which significantly suppress the critical slowing down. We analyze the Monte Carlo data by finite-size scaling, and locate a line of critical points of the Ising universality class in the region J<0 and K>0, with K the four-spin interaction. Further, we find that, in the J→-∞ limit, the critical line terminates at the decoupled point K=0. From the numerical results and the exact mapping, we conjecture that this "tricritical" point (J→-∞,K=0) is Berezinsky-Kosterlitz-Thouless-like and the logarithmic correction is absent. The dynamic critical exponent of the worm algorithm is estimated as z=0.28(1) near (J→-∞,K=0).

arXiv:1009.3172•Physical Review E•Aug 17, 2011•6 citations

Finite-size scaling of Monte Carlo simulations for the fcc Ising antiferromagnet: Effects of the low-temperature phase degeneracy

Ronja Stübel, Wolfhard Janke E

The Ising antiferromagnet on a face-centered cubic (fcc) lattice with nearest-neighbor interaction only is well known to exhibit a macroscopic (exponential in the system size $L$) ground-state degeneracy. With increasing temperature, this degeneracy is expected to be lifted and the model undergoes a first-order phase transition. For a model with an exponential degeneracy in the whole low-temperature phase, it was recently found that the finite-size scaling behavior is governed by leading correction terms $\ensuremath{\sim}{L}^{\ensuremath{-}2}$ instead of $\ensuremath{\sim}{L}^{\ensuremath{-}3}$ as usual. To test the conjecture that such a transmuted behavior may effectively persist also for the fcc antiferromagnet up to some crossover system size, we have performed parallel multicanonical Monte Carlo simulations for lattices of linear size $L\ensuremath{\le}18$ with periodic boundary conditions and determined various inverse pseudo phase transition temperatures, as well as the extremal values of the specific heat and the energetic Binder parameter. We indeed find that, for the simulated lattice sizes, the conjectured transmuted finite-size scaling ansatz fits the data better than the standard ansatz. On this basis, we extrapolate for the transition temperature an estimate of ${T}_{0}=1.735047(46)$.

arXiv:1810.05006•Physical review. B./Physical review. B•Nov 9, 2018•7 citations

Critical properties of the three-dimensional equivalent-neighbor model and crossover scaling in finite systems

Erik Luijten

Accurate numerical results are presented for the three-dimensional equivalent-neighbor model on a cubic lattice, for 12 different interaction ranges (coordination number between 18 and 250). These results allow the determination of the range dependences of the critical temperature and various critical amplitudes, which are compared to renormalization-group predictions. In addition, the analysis yields an estimate for the interaction range at which the leading corrections to scaling vanish for the spin-1 / 2 model, and confirms earlier conclusions that the leading Wegner correction must be negative for the three-dimensional (nearest-neighbor) Ising model. By complementing these results with Monte Carlo data for systems with coordination numbers as large as 52 514, the full finite-size crossover curves between classical and Ising-like behavior are obtained as a function of a generalized Ginzburg parameter. Also, the crossover function for the effective magnetic exponent is determined.

arXiv:cond-mat/9811332•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•May 1, 1999•26 citations
Monte Carlo simulations have been used to study a model consisting of three Ising models with nearest neighbor exchange J on a face-centered-cubic lattice, which are coupled together by a constraint. The constraint requires that the vector sum of the four spins located at the corners of any elementary tetrahedron of the lattice not add to zero. Finite-size scaling analysis of the Monte Carlo results for L×L×L lattices with L=16, 32, and 64 has been used to find that Tc /J=10.84±0.01. The finite-size scaling functions for the magnetization and magnetic susceptibility near Tc are consistent with standard m=3 Heisenberg critical exponents, as predicted by renormalization group theory. Strengthening the four-spin constraint increases the value of the effective negative Ashkin–Teller four-spin coupling, and drives Tc to infinity.
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] .
The Monte Carlo renormalization group is applied to the three-dimensional Ising model on simple cubic lattices with 83{8}^{3}83, 163{16}^{3}163, 323{32}^{3}323, and 643{64}^{3}643 sites. The comparison of block-spin correlation functions from the largest lattices yields the nearest-neighbor critical coupling K1c=0.221654(6){K}_{1}^{c}=0.221654(6)K1c​=0.221654(6). After allowing for (i) interpolation to this best estimate for K1c{K}_{1}^{c}K1c​, (ii) an apparent finite-size effect in the renormalization-group transformation due to the measurement of correlation functions of too few (seven) operators, and (iii) the extrapolation for the effect of a slow transient towards the fixed point, the values $\ensuremath{\nu}=0.629(4)$ and $\ensuremath{\eta}=0.031(5)$ are obtained for the thermal and magnetic exponents. The correction-to-scaling exponent $\ensuremath{\omega}$ is estimated to be around 1; to obtain an accuracy competitive with other methods requires measurements with more than seven operators. We briefly review the problem of redundant operators and indicate the future prospects for this kind of calculation.
The results of the large-scale Monte Carlo simulations of the continuous phase transitions in the three-dimensional Ashkin-Teller model on a cubic lattice are presented in the two-parameter space. The invariance of the ratio of the square of the second moment of the order parameter to its fourth moment in the critical region and the Metropolis algorithm are exploited. Using the universality hypothesis and the finite-size-scaling analysis, the accurate positions of the continuous phase-transition lines are calculated and their Ising character is confirmed. Approximate locations of the tricritical points are also estimated.
We investigate the symmetric Ashkin-Teller (AT) model on the triangular lattice in the antiferromagnetic two-spin coupling region (J<0). In the J→-∞ limit, we map the AT model onto a fully packed loop-dimer model on the honeycomb lattice. On the basis of this exact transformation and the low-temperature expansion, we formulate a variant of worm-type algorithms for the AT model, which significantly suppress the critical slowing down. We analyze the Monte Carlo data by finite-size scaling, and locate a line of critical points of the Ising universality class in the region J<0 and K>0, with K the four-spin interaction. Further, we find that, in the J→-∞ limit, the critical line terminates at the decoupled point K=0. From the numerical results and the exact mapping, we conjecture that this "tricritical" point (J→-∞,K=0) is Berezinsky-Kosterlitz-Thouless-like and the logarithmic correction is absent. The dynamic critical exponent of the worm algorithm is estimated as z=0.28(1) near (J→-∞,K=0).
The Ising antiferromagnet on a face-centered cubic (fcc) lattice with nearest-neighbor interaction only is well known to exhibit a macroscopic (exponential in the system size LLL) ground-state degeneracy. With increasing temperature, this degeneracy is expected to be lifted and the model undergoes a first-order phase transition. For a model with an exponential degeneracy in the whole low-temperature phase, it was recently found that the finite-size scaling behavior is governed by leading correction terms $\ensuremath{\sim}{L}^{\ensuremath{-}2}$ instead of $\ensuremath{\sim}{L}^{\ensuremath{-}3}$ as usual. To test the conjecture that such a transmuted behavior may effectively persist also for the fcc antiferromagnet up to some crossover system size, we have performed parallel multicanonical Monte Carlo simulations for lattices of linear size $L\ensuremath{\le}18$ with periodic boundary conditions and determined various inverse pseudo phase transition temperatures, as well as the extremal values of the specific heat and the energetic Binder parameter. We indeed find that, for the simulated lattice sizes, the conjectured transmuted finite-size scaling ansatz fits the data better than the standard ansatz. On this basis, we extrapolate for the transition temperature an estimate of T0=1.735047(46){T}_{0}=1.735047(46)T0​=1.735047(46).
Accurate numerical results are presented for the three-dimensional equivalent-neighbor model on a cubic lattice, for 12 different interaction ranges (coordination number between 18 and 250). These results allow the determination of the range dependences of the critical temperature and various critical amplitudes, which are compared to renormalization-group predictions. In addition, the analysis yields an estimate for the interaction range at which the leading corrections to scaling vanish for the spin-1 / 2 model, and confirms earlier conclusions that the leading Wegner correction must be negative for the three-dimensional (nearest-neighbor) Ising model. By complementing these results with Monte Carlo data for systems with coordination numbers as large as 52 514, the full finite-size crossover curves between classical and Ising-like behavior are obtained as a function of a generalized Ginzburg parameter. Also, the crossover function for the effective magnetic exponent is determined.