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:cond-mat/9609007 · Physical review. B, Condensed matter

Zero-temperature phase transitions of an antiferromagnetic Ising model of general spin on a triangular lattice

Chen Zeng, Christopher L. Henley

Published Jun 1, 1997• 61 citations
arXivHTMLPDF

Abstract

We map the ground-state ensemble of antiferromagnetic Ising model of spin-S on a triangular lattice to an interface model whose entropic fluctuations are proposed to be described by an effective Gaussian free energy, which enables us to calculate the critical exponents of various operators in terms of the stiffness constant of the interface. Monte Carlo simulations for the ground-state ensemble utilizing this interfacial representation are performed to study both the dynamical and the static properties of the model. This method yields more accurate numerical results for the critical exponents. By varying the spin magnitude in the model, we find that the model exhibits three phases with a Kosterlitz-Thouless phase transition at ${\mathrm{S}}_{\mathrm{KT}}$2 and a locking phase transition at ${\mathrm{S}}_{\mathrm{L}}$\ensuremath{\leqslant}3. The phase diagram at finite temperatures is also discussed.

DOI
10.1103/physrevb.55.14935

Related papers

CRITICAL COARSENING DYNAMICS OF SPIN S =1/2, 3/2 ANTIFERROMAGNETIC ISING MODEL ON TRIANGULAR LATTICE

Kwanghoon Chung, Mookyung Cheon, Iksoo Chang

The critical coarsening dynamics of the spin S =1/2, 3/2 antiferromagnetic Ising model on a triangular lattice is studied by the dynamic Monte Carlo simulation using a heat bath algorithm. The triangular antiferromagnetic Ising (TAI) model possesses an intrinsic geometrical frustration and a large degeneracy of ground state which may affect the equilibrium and non-equilibrium critical behaviors. The S =1/2 TAI has no phase transition at a finite temperature, but it was suggested that the S =3/2 TAI has the Kosterlitz–Thouless (KT)-type phase transition at a finite temperature. We employ a finite size scaling approach for the correlation function from the short-time dynamics of the S =1/2, 3/2 TAI to calculate the values of the static critical exponent η and the dynamic exponent z. For the S =1/2 TAI, our dynamic scaling analysis provides η =0.498±0.006 and z =2.278±0.020 at T =0, agreeing with the previous results. For the S =3/2 TAI, after identifying a KT-transition temperature TKT =0.51±0.01, we find the temperature-dependent η ranging from 0.301±0.008 at T =0.51 to 0.224±0.016 at T =0 along the KT-line whereas the value of z =2.20±0.06 is constant for T≤TKT. It is shown that the spin S =3/2 TAI model and the two-dimensional XY model, sharing the KT-type phase transition, exhibit similar static critical and coarsening dynamics behavior. For both the S =1/2, 3/2 TAI, the value of z (η) is compatible with (larger than) that of the Ising model at Tc and the XY model for T≤TKT in two-dimension. Our results imply that although the quasi-long-range order disappears with ηXY =0 in the two-dimensional XY model at T =0, the S =3/2 TAI still maintains it with η TAI =0.224 due to the effect of a frustration and a high degeneracy of ground state.

International Journal of Modern Physics C•Apr 1, 2006•1 citation

Determination of the Critical Points of Antiferromagnetic Ising Model with Next Nearest Neighbour Interactions on the Triangular Lattice

Seiji Miyashita, Hidetsugu Kitatani, Yasumasa Kanada

Positions of critical temperatures of antiferromagnetic Ising model with ferromagnetic next nearest neighbour interactions on the triangular lattice are investigated by making use of their properties of the Kosterlitz-Thouless type of phase transition. In particular, finite-size scaling method and Monte Carlo Renormalization Group methods (MCRG) with flow diagrams of the anomalous dimension η show the fixed line structure where the correlation length diverges. Another MCRG with flow diagram of the stiffness constant also shows a critical property clearly. The temperature dependence of the long range order is also discussed in order to locate the lower critical point.

Journal of the Physical Society of Japan•May 15, 1991•25 citations

Antiferromagnetic triangular Ising model: Critical behavior of the ground state

Henk W. J. Blöte, M. P. Nightingale

The critical ground state of the antiferromagnetic triangular Ising model is studied with finite-size scaling, with the use of exact techniques and numerical transfer-matrix methods. Exact solution of the model in zero field yields a conformal anomaly equal to 1, and values of critical exponents. For nonzero fields, dominant eigenvalues of the transfer matrix are computed for systems with linear dimensions up to 27 with the use of an efficient coding of spin states as a system of nonintersecting strings. The results are in perfect agreement with the Gaussian-model--Coulomb-gas scenario proposed by Nienhuis et al. for this model. Quantitative agreement is found with the predictions for exponents of spin-wave and vortex operators. Furthermore, we locate the field-induced Kosterlitz-Thouless transition to a long-range-ordered state at a field ${\mathit{H}}_{\mathit{c}}$=0.266\ifmmode\pm\else\textpm\fi{}0.010.

Physical review. B, Condensed matter•Jun 1, 1993•64 citations

Phase diagram of the quantum Ising model on a triangular lattice under external field

Yuan Da Liao, Han Li, Zheng Yan +4 more

Quantum Ising model on a triangular lattice hosts a finite temperature Berezinskii-Kosterlitz-Thouless (BKT) phase with emergent U(1) symmetry, and it will transit into an up-up-down (UUD) phase with ${C}_{3}$ symmetry breaking upon an infinitesimal external field along the longitudinal direction, but the overall phase diagram spanned by the axes of external field and temperature remains opaque due to the lack of systematic investigations with controlled methodologies. By means of quantum Monte Carlo at finite temperature and ground state density matrix renormalization group simulations, we map out the phase diagram of triangular quantum Ising model. Starting from the upper BKT temperature at zero field, we obtain the phase boundary between the UUD and paramagnetic phases with its 2D $q=3$ Potts universality at weak field and weakly first order transition at strong field. Originated from the lower BKT temperature at zero field, we analyze the low temperature phase boundary between the clock phase and the UUD phase with Ising symmetry breaking at weak fields and the quantum phase transition between the UUD and fully polarized phases at strong fields. The accurate many-body numerical results are consistent with our field theoretical analysis. The experimental relevance towards the BKT magnet ${\mathrm{TmMgGaO}}_{4}$ and programmable quantum simulators are also discussed.

arXiv:2101.11610•Physical review. B./Physical review. B•Mar 11, 2021•26 citations

Monte Carlo Simulation of the Antiferromagnetic Ising Model on the Triangular Lattice with the First and Second Neighbour Interactions

Sumiyoshi Fujiki, Keisei Shutoh, Yoshihiko Abe +1 more

Monte Carlo simulations of the Ising model on the triangular lattice with the antiferromagnetic nn interaction J and the ferromagnetic nnn interaction J ' are carried out. Several physical quantities (magnetizations, energies, susceptibilities, specific heats of sublattices and of the total system) are calculated. A phase which resembles the partially disordered phase proposed by Mekata is observed in the finite temperature region. But it is shown in this region that the sublattice magnetizations vanish at the thermodynamic limit and the Kosterlitz-Thouless phase transition is predicted by the finite size analysis.

Journal of the Physical Society of Japan•May 1, 1983•61 citations

Interface Approach to Phase Transitions and Ordering by Monte Carlo Simulations and Its Applications to Three-Dimensional Antiferromagnetic Potts Models

Yohtaro Ueno, Gang Sun, Ikuo Ono

We develop an interface approach to bulk phase transitions and ordered states by using Monte Garlo simulations, and apply it to antiferromagnetic q -state Potts models with q =3∼6 on the simple cubic lattice. A stiffness exponent a defined by Δ F ∼ L a , where Δ F is the interface free energy for a system of size L , is introduced as a measure of the stiffness of the ordered phase against an external stress. Applying finite-size scaling to Δ F and to the squared order parameter enable one to determine each of T c , ν, β and γ in order where T c is the critical temperature and the others are critical exponents. This approach also provides a means to study properties of the ordered state through a and interface profiles. In the q =3 and 4 models the estimated exponents ν, β, γ and a indicate new universality classes. The q =5 model shows a phase transition, while the q =6 model does not. The q =3∼5 models have non-integer values of a at low temperatures. The phase transitions in these models are confirmed to be due to entropy gains.

Journal of the Physical Society of Japan•Apr 1, 1989•47 citations
We map the ground-state ensemble of antiferromagnetic Ising model of spin-S on a triangular lattice to an interface model whose entropic fluctuations are proposed to be described by an effective Gaussian free energy, which enables us to calculate the critical exponents of various operators in terms of the stiffness constant of the interface. Monte Carlo simulations for the ground-state ensemble utilizing this interfacial representation are performed to study both the dynamical and the static properties of the model. This method yields more accurate numerical results for the critical exponents. By varying the spin magnitude in the model, we find that the model exhibits three phases with a Kosterlitz-Thouless phase transition at SKT{\mathrm{S}}_{\mathrm{KT}}SKT​2 and a locking phase transition at SL{\mathrm{S}}_{\mathrm{L}}SL​\ensuremath{\leqslant}3. The phase diagram at finite temperatures is also discussed.
The critical coarsening dynamics of the spin S =1/2, 3/2 antiferromagnetic Ising model on a triangular lattice is studied by the dynamic Monte Carlo simulation using a heat bath algorithm. The triangular antiferromagnetic Ising (TAI) model possesses an intrinsic geometrical frustration and a large degeneracy of ground state which may affect the equilibrium and non-equilibrium critical behaviors. The S =1/2 TAI has no phase transition at a finite temperature, but it was suggested that the S =3/2 TAI has the Kosterlitz–Thouless (KT)-type phase transition at a finite temperature. We employ a finite size scaling approach for the correlation function from the short-time dynamics of the S =1/2, 3/2 TAI to calculate the values of the static critical exponent η and the dynamic exponent z. For the S =1/2 TAI, our dynamic scaling analysis provides η =0.498±0.006 and z =2.278±0.020 at T =0, agreeing with the previous results. For the S =3/2 TAI, after identifying a KT-transition temperature TKT =0.51±0.01, we find the temperature-dependent η ranging from 0.301±0.008 at T =0.51 to 0.224±0.016 at T =0 along the KT-line whereas the value of z =2.20±0.06 is constant for T≤TKT. It is shown that the spin S =3/2 TAI model and the two-dimensional XY model, sharing the KT-type phase transition, exhibit similar static critical and coarsening dynamics behavior. For both the S =1/2, 3/2 TAI, the value of z (η) is compatible with (larger than) that of the Ising model at Tc and the XY model for T≤TKT in two-dimension. Our results imply that although the quasi-long-range order disappears with ηXY =0 in the two-dimensional XY model at T =0, the S =3/2 TAI still maintains it with η TAI =0.224 due to the effect of a frustration and a high degeneracy of ground state.
Positions of critical temperatures of antiferromagnetic Ising model with ferromagnetic next nearest neighbour interactions on the triangular lattice are investigated by making use of their properties of the Kosterlitz-Thouless type of phase transition. In particular, finite-size scaling method and Monte Carlo Renormalization Group methods (MCRG) with flow diagrams of the anomalous dimension η show the fixed line structure where the correlation length diverges. Another MCRG with flow diagram of the stiffness constant also shows a critical property clearly. The temperature dependence of the long range order is also discussed in order to locate the lower critical point.
The critical ground state of the antiferromagnetic triangular Ising model is studied with finite-size scaling, with the use of exact techniques and numerical transfer-matrix methods. Exact solution of the model in zero field yields a conformal anomaly equal to 1, and values of critical exponents. For nonzero fields, dominant eigenvalues of the transfer matrix are computed for systems with linear dimensions up to 27 with the use of an efficient coding of spin states as a system of nonintersecting strings. The results are in perfect agreement with the Gaussian-model--Coulomb-gas scenario proposed by Nienhuis et al. for this model. Quantitative agreement is found with the predictions for exponents of spin-wave and vortex operators. Furthermore, we locate the field-induced Kosterlitz-Thouless transition to a long-range-ordered state at a field Hc{\mathit{H}}_{\mathit{c}}Hc​=0.266\ifmmode\pm\else\textpm\fi{}0.010.
Quantum Ising model on a triangular lattice hosts a finite temperature Berezinskii-Kosterlitz-Thouless (BKT) phase with emergent U(1) symmetry, and it will transit into an up-up-down (UUD) phase with C3{C}_{3}C3​ symmetry breaking upon an infinitesimal external field along the longitudinal direction, but the overall phase diagram spanned by the axes of external field and temperature remains opaque due to the lack of systematic investigations with controlled methodologies. By means of quantum Monte Carlo at finite temperature and ground state density matrix renormalization group simulations, we map out the phase diagram of triangular quantum Ising model. Starting from the upper BKT temperature at zero field, we obtain the phase boundary between the UUD and paramagnetic phases with its 2D q=3q=3q=3 Potts universality at weak field and weakly first order transition at strong field. Originated from the lower BKT temperature at zero field, we analyze the low temperature phase boundary between the clock phase and the UUD phase with Ising symmetry breaking at weak fields and the quantum phase transition between the UUD and fully polarized phases at strong fields. The accurate many-body numerical results are consistent with our field theoretical analysis. The experimental relevance towards the BKT magnet TmMgGaO4{\mathrm{TmMgGaO}}_{4}TmMgGaO4​ and programmable quantum simulators are also discussed.
Monte Carlo simulations of the Ising model on the triangular lattice with the antiferromagnetic nn interaction J and the ferromagnetic nnn interaction J ' are carried out. Several physical quantities (magnetizations, energies, susceptibilities, specific heats of sublattices and of the total system) are calculated. A phase which resembles the partially disordered phase proposed by Mekata is observed in the finite temperature region. But it is shown in this region that the sublattice magnetizations vanish at the thermodynamic limit and the Kosterlitz-Thouless phase transition is predicted by the finite size analysis.
We develop an interface approach to bulk phase transitions and ordered states by using Monte Garlo simulations, and apply it to antiferromagnetic q -state Potts models with q =3∼6 on the simple cubic lattice. A stiffness exponent a defined by Δ F ∼ L a , where Δ F is the interface free energy for a system of size L , is introduced as a measure of the stiffness of the ordered phase against an external stress. Applying finite-size scaling to Δ F and to the squared order parameter enable one to determine each of T c , ν, β and γ in order where T c is the critical temperature and the others are critical exponents. This approach also provides a means to study properties of the ordered state through a and interface profiles. In the q =3 and 4 models the estimated exponents ν, β, γ and a indicate new universality classes. The q =5 model shows a phase transition, while the q =6 model does not. The q =3∼5 models have non-integer values of a at low temperatures. The phase transitions in these models are confirmed to be due to entropy gains.