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

Tensor network simulation for the frustrated J1−J2 Ising model on the square lattice

Hong Li, Li-Ping Yang

Published Aug 16, 2021• 33 citations
arXivHTMLPDF

Abstract

By using extensive tensor network calculations, we map out the phase diagram of the frustrated J_{1}-J_{2} Ising model on the square lattice. In particular, we focus on the cases with controversy in the phase diagram, especially the stripe transition in the regime g=|J_{2}/J_{1}|>1/2 (J_{2}>0,J_{1}<0). While recent studies claimed that the phase transition is of first order when 1/2<g<g^{*} (with the smallest g^{*} being 0.67), our simulations suggest that if there is such a first-order region, it is smaller than those found in earlier studies by other methods. Combining with the analysis of critical properties, we provide evidence that the classical J_{1}-J_{2} model evolves continuously from two decoupled Ising models (g→∞ with central charge c=1) to a point belonging to the tricritical Ising universality class (with c=0.7) as g decreases to g^{*}≃0.54.

DOI
10.1103/physreve.104.024118

Related papers

Phase transitions in the frustrated Ising model on the square lattice

Songbo Jin, Arnab Sen, Wenan Guo +1 more

We consider the thermal phase transition from a paramagnetic to a stripe-antiferromagnetic phase in the frustrated two-dimensional square-lattice Ising model with competing interactions ${J}_{1}<0$ (nearest neighbor, ferromagnetic) and ${J}_{2}>0$ (second neighbor, antiferromagnetic). The striped phase breaks a ${Z}_{4}$ symmetry and is stabilized at low temperatures for $g={J}_{2}/|{J}_{1}|>1/2$. Despite the simplicity of the model, it has proved difficult to precisely determine the order and the universality class of the phase transitions. This was done convincingly only recently by Jin et al. [Phys. Rev. Lett. 108, 045702 (2012)]. Here, we further elucidate the nature of these transitions and their anomalies by employing a combination of cluster mean-field theory, Monte Carlo simulations, and transfer-matrix calculations. The ${J}_{1}$-${J}_{2}$ model has a line of very weak first-order phase transitions in the whole region $1/2<g<{g}^{*}$, where ${g}^{*}=0.67\ifmmode\pm\else\textpm\fi{}0.01$. Thereafter, the transitions from $g={g}^{*}$ to $g\ensuremath{\rightarrow}\ensuremath{\infty}$ are continuous and can be fully mapped, using universality arguments, to the critical line of the well-known Ashkin-Teller model from its four-state Potts point to the decoupled Ising limit. We also comment on the pseudo-first-order behavior at the Potts point and its neighborhood in the Ashkin-Teller model on finite lattices, which in turn leads to the appearance of similar effects in the vicinity of the multicritical point ${g}^{*}$ in the ${J}_{1}$-${J}_{2}$ model. The continuous transitions near ${g}^{*}$ can therefore be mistaken for first-order transitions, and this realization was the key to understanding the paramagnetic-striped transition for the full range of $g>1/2$. Most of our conclusions are based on Monte Carlo calculations, while the cluster mean-field and transfer-matrix results provide useful methodological benchmarks for weakly first-order behaviors and Ashkin-Teller criticality.

arXiv:1212.5339•Physical Review B•Apr 9, 2013•116 citations

Emergent intermediate phase in the J 1 − J 2 XY model from tensor network approaches

Feng-Feng Song, Hanggai Nuomin, Naoki Kawashima

We investigate the finite-temperature phase diagram of the classical ${J}_{1}\ensuremath{-}{J}_{2}$ XY model on a square lattice using a tensor network approach designed for frustrated spin systems. This model, characterized by competing nearest-neighbor and next-to-nearest-neighbor interactions, exhibits a complex interplay between $U(1)$ and ${Z}_{2}$ symmetries. Our study reveals an emergent intermediate phase around ${J}_{2}/{J}_{1}\ensuremath{\sim}0.505$, which is characterized by a ${Z}_{2}$ long-range stripe order without phase coherence in the XY spins. The intermediate phase features two well-separated phase transitions: a higher-temperature Ising transition and a lower-temperature Berezinskii-Kosterlitz-Thouless (BKT) transition. The relative separation between these transitions is significantly larger than previously reported, enabling a clearer investigation of their distinct thermodynamic properties. For $0.5<{J}_{2}/{J}_{1}<0.501$, two transitions merge into a single first-order phase transition, a phenomenon that cannot be explained solely by mapping to the Ising-XY model. As ${J}_{2}/{J}_{1}\ensuremath{\rightarrow}\ensuremath{\infty}$, the transition evolves continuously into the BKT universality class. These findings advance the understanding of the mechanisms driving phase transitions in frustrated spin systems, and they suggest potential experimental realizations in platforms such as ultracold atoms, Josephson junction arrays, and optical lattices.

arXiv:2412.18892•Physical review. B./Physical review. B•Mar 20, 2025•1 citation

Analysis of the phase transition for the Ising model on the frustrated square lattice

Ansgar Kalz, A. Honecker, Marion Moliner

We analyze the phase transition of the frustrated ${J}_{1}$-${J}_{2}$ Ising model with antiferromagnetic nearest- and strong next-nearest-neighbor interactions on the square lattice. Using extensive Monte Carlo simulations we show that the nature of the phase transition for $1/2<{J}_{2}/{J}_{1}\ensuremath{\lesssim}1$ is not of the weakly universal type---as commonly believed---but we conclude from the clearly doubly peaked structure of the energy histograms that the transition is of weak first order. Motivated by these results, we analyze the phase transitions via field-theoretic methods; i.e., we calculate the central charge of the underlying field theory via transfer-matrix techniques and present, furthermore, a field-theoretic discussion on the phase-transition behavior of the model. Starting from the conformally invariant fixed point of two decoupled critical Ising models (${J}_{1}=0$), we calculate the effect of the nearest-neighbor coupling term perturbatively using operator product expansions. As an effective action we obtain the Ashkin-Teller model.

arXiv:1105.4836•Physical Review B•Nov 9, 2011•55 citations

First-Order Transition and Marginal Critical Behavior in a 2D Frustrated Ising Model with Efficient Tensor-Network Implementation

Christophe Chatelain

The phase diagram of a novel two-dimensional frustrated Ising model with both anti-ferromagnetic and ferromagnetic couplings is studied using Tensor-Network Renormalization-Group techniques. This model can be seen as two anti-ferromagnetic Ising replicas coupled by non-local spin-spin interactions, designed in such a way that the continuum limit matches that of the still debated J1 -J2 model and induces a marginal critical behavior. Our model has the advantage of having more symmetries than the J1 -J2 model and of allowing a more straightforward implementation of Tensor-Network Renormalization-Group algorithms We demonstrate the existence of two transition lines, featuring both first and second-order regimes. In the latter, the central charge and the critical exponents are shown to be compatible with the Ashkin-Teller universality class. This picture is consistent with that given by Monte Carlo simulations of the J1 -J2 model but not with recent studies with Tensor-Network techniques.

•Oct 2, 2024•0 citations

Higher-order tensor renormalization group study of the J1−J2 Ising model on a square lattice

Kota Yoshiyama, Koji Hukushima

Phase transitions of the J_{1}-J_{2} Ising model on a square lattice are studied using the higher-order tensor renormalization group (HOTRG) method. This system involves a competition between the ferromagnetic interaction J_{1} and antiferromagnetic interaction J_{2}, and in previous studies, weak first-order and second-order transitions were observed near the ratio g=J_{2}/|J_{1}|=1/2. It has also been suggested that the universality class of the second-order phase transition connected to the first-order transition line for g>1/2 belongs to the Ashkin-Teller class, which is characterized by a continuously varying critical exponent with g, as predicted by field-theoretical and other studies. Our results, based on the HOTRG calculations for significantly larger sizes, indicate that the region of the first-order transition is marginally narrower than that in previous studies. Furthermore, it is suggested that the region where the critical exponent changes does not necessarily coincide with the Ashkin-Teller region.

arXiv:2303.07733•Physical review. E•Nov 15, 2023•22 citations

Tensor network simulation of the phase diagram of the frustrated J1 - J2 Heisenberg model on a checkerboard lattice

Yang‐Hao Chan, Yong‐Jian Han, L.-M. Duan

We use the recently developed tensor network algorithm based on infinite projected entangled pair states (iPEPS) to study the phase diagram of frustrated antiferromagnetic ${J}_{1}$-${J}_{2}$ Heisenberg model on a checkerboard lattice. The simulation indicates a N\'eel ordered phase when ${J}_{2}<0.88{J}_{1}$, a plaquette valence bond solid state when $0.88<{J}_{2}/{J}_{1}<1.11$, and a stripe phase when ${J}_{2}>1.11{J}_{1}$, with two first-order transitions across the phase boundaries. The calculation shows that the cross dimer state proposed before is unlikely to be the ground state of the model, although such a state indeed arises as a metastable state in some parameter region.

arXiv:1108.1595•Physical Review B•Dec 12, 2011•16 citations
By using extensive tensor network calculations, we map out the phase diagram of the frustrated J_{1}-J_{2} Ising model on the square lattice. In particular, we focus on the cases with controversy in the phase diagram, especially the stripe transition in the regime g=|J_{2}/J_{1}|>1/2 (J_{2}>0,J_{1}<0). While recent studies claimed that the phase transition is of first order when 1/2<g<g^{} (with the smallest g^{} being 0.67), our simulations suggest that if there is such a first-order region, it is smaller than those found in earlier studies by other methods. Combining with the analysis of critical properties, we provide evidence that the classical J_{1}-J_{2} model evolves continuously from two decoupled Ising models (g→∞ with central charge c=1) to a point belonging to the tricritical Ising universality class (with c=0.7) as g decreases to g^{*}≃0.54.
We consider the thermal phase transition from a paramagnetic to a stripe-antiferromagnetic phase in the frustrated two-dimensional square-lattice Ising model with competing interactions J1<0{J}_{1}<0J1​<0 (nearest neighbor, ferromagnetic) and J2>0{J}_{2}>0J2​>0 (second neighbor, antiferromagnetic). The striped phase breaks a Z4{Z}_{4}Z4​ symmetry and is stabilized at low temperatures for g=J2/∣J1∣>1/2g={J}_{2}/|{J}_{1}|>1/2g=J2​/∣J1​∣>1/2. Despite the simplicity of the model, it has proved difficult to precisely determine the order and the universality class of the phase transitions. This was done convincingly only recently by Jin et al. [Phys. Rev. Lett. 108, 045702 (2012)]. Here, we further elucidate the nature of these transitions and their anomalies by employing a combination of cluster mean-field theory, Monte Carlo simulations, and transfer-matrix calculations. The J1{J}_{1}J1​-J2{J}_{2}J2​ model has a line of very weak first-order phase transitions in the whole region 1/2<g<g∗1/2<g<{g}^{*}1/2<g<g∗, where ${g}^{*}=0.67\ifmmode\pm\else\textpm\fi{}0.01$. Thereafter, the transitions from g=g∗g={g}^{*}g=g∗ to $g\ensuremath{\rightarrow}\ensuremath{\infty}$ are continuous and can be fully mapped, using universality arguments, to the critical line of the well-known Ashkin-Teller model from its four-state Potts point to the decoupled Ising limit. We also comment on the pseudo-first-order behavior at the Potts point and its neighborhood in the Ashkin-Teller model on finite lattices, which in turn leads to the appearance of similar effects in the vicinity of the multicritical point g∗{g}^{*}g∗ in the J1{J}_{1}J1​-J2{J}_{2}J2​ model. The continuous transitions near g∗{g}^{*}g∗ can therefore be mistaken for first-order transitions, and this realization was the key to understanding the paramagnetic-striped transition for the full range of g>1/2g>1/2g>1/2. Most of our conclusions are based on Monte Carlo calculations, while the cluster mean-field and transfer-matrix results provide useful methodological benchmarks for weakly first-order behaviors and Ashkin-Teller criticality.
We investigate the finite-temperature phase diagram of the classical ${J}_{1}\ensuremath{-}{J}_{2}$ XY model on a square lattice using a tensor network approach designed for frustrated spin systems. This model, characterized by competing nearest-neighbor and next-to-nearest-neighbor interactions, exhibits a complex interplay between U(1)U(1)U(1) and Z2{Z}_{2}Z2​ symmetries. Our study reveals an emergent intermediate phase around ${J}_{2}/{J}_{1}\ensuremath{\sim}0.505$, which is characterized by a Z2{Z}_{2}Z2​ long-range stripe order without phase coherence in the XY spins. The intermediate phase features two well-separated phase transitions: a higher-temperature Ising transition and a lower-temperature Berezinskii-Kosterlitz-Thouless (BKT) transition. The relative separation between these transitions is significantly larger than previously reported, enabling a clearer investigation of their distinct thermodynamic properties. For 0.5<J2/J1<0.5010.5<{J}_{2}/{J}_{1}<0.5010.5<J2​/J1​<0.501, two transitions merge into a single first-order phase transition, a phenomenon that cannot be explained solely by mapping to the Ising-XY model. As ${J}_{2}/{J}_{1}\ensuremath{\rightarrow}\ensuremath{\infty}$, the transition evolves continuously into the BKT universality class. These findings advance the understanding of the mechanisms driving phase transitions in frustrated spin systems, and they suggest potential experimental realizations in platforms such as ultracold atoms, Josephson junction arrays, and optical lattices.
We analyze the phase transition of the frustrated J1{J}_{1}J1​-J2{J}_{2}J2​ Ising model with antiferromagnetic nearest- and strong next-nearest-neighbor interactions on the square lattice. Using extensive Monte Carlo simulations we show that the nature of the phase transition for $1/2<{J}_{2}/{J}_{1}\ensuremath{\lesssim}1$ is not of the weakly universal type---as commonly believed---but we conclude from the clearly doubly peaked structure of the energy histograms that the transition is of weak first order. Motivated by these results, we analyze the phase transitions via field-theoretic methods; i.e., we calculate the central charge of the underlying field theory via transfer-matrix techniques and present, furthermore, a field-theoretic discussion on the phase-transition behavior of the model. Starting from the conformally invariant fixed point of two decoupled critical Ising models (J1=0{J}_{1}=0J1​=0), we calculate the effect of the nearest-neighbor coupling term perturbatively using operator product expansions. As an effective action we obtain the Ashkin-Teller model.
The phase diagram of a novel two-dimensional frustrated Ising model with both anti-ferromagnetic and ferromagnetic couplings is studied using Tensor-Network Renormalization-Group techniques. This model can be seen as two anti-ferromagnetic Ising replicas coupled by non-local spin-spin interactions, designed in such a way that the continuum limit matches that of the still debated J1 -J2 model and induces a marginal critical behavior. Our model has the advantage of having more symmetries than the J1 -J2 model and of allowing a more straightforward implementation of Tensor-Network Renormalization-Group algorithms We demonstrate the existence of two transition lines, featuring both first and second-order regimes. In the latter, the central charge and the critical exponents are shown to be compatible with the Ashkin-Teller universality class. This picture is consistent with that given by Monte Carlo simulations of the J1 -J2 model but not with recent studies with Tensor-Network techniques.
Phase transitions of the J_{1}-J_{2} Ising model on a square lattice are studied using the higher-order tensor renormalization group (HOTRG) method. This system involves a competition between the ferromagnetic interaction J_{1} and antiferromagnetic interaction J_{2}, and in previous studies, weak first-order and second-order transitions were observed near the ratio g=J_{2}/|J_{1}|=1/2. It has also been suggested that the universality class of the second-order phase transition connected to the first-order transition line for g>1/2 belongs to the Ashkin-Teller class, which is characterized by a continuously varying critical exponent with g, as predicted by field-theoretical and other studies. Our results, based on the HOTRG calculations for significantly larger sizes, indicate that the region of the first-order transition is marginally narrower than that in previous studies. Furthermore, it is suggested that the region where the critical exponent changes does not necessarily coincide with the Ashkin-Teller region.
We use the recently developed tensor network algorithm based on infinite projected entangled pair states (iPEPS) to study the phase diagram of frustrated antiferromagnetic J1{J}_{1}J1​-J2{J}_{2}J2​ Heisenberg model on a checkerboard lattice. The simulation indicates a N\'eel ordered phase when J2<0.88J1{J}_{2}<0.88{J}_{1}J2​<0.88J1​, a plaquette valence bond solid state when 0.88<J2/J1<1.110.88<{J}_{2}/{J}_{1}<1.110.88<J2​/J1​<1.11, and a stripe phase when J2>1.11J1{J}_{2}>1.11{J}_{1}J2​>1.11J1​, with two first-order transitions across the phase boundaries. The calculation shows that the cross dimer state proposed before is unlikely to be the ground state of the model, although such a state indeed arises as a metastable state in some parameter region.