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:2505.05194 · Physical Review Letters

Two-Dimensional J 1 − J 2 Clock Model: Enhanced Symmetries, Emergent Orders, and Landau-Incompatible Transitions

Pulloor Kuttanikkad Vishnu, Abhishodh Prakash, Rajesh Narayanan, Titas Chanda

Published Nov 26, 2025• 0 citations
arXivHTMLPDF

Abstract

We present a comprehensive study on the frustrated J_{1}-J_{2} classical q-state clock model with even q>4 on a two-dimensional square lattice, revealing a rich ensemble of phases driven by competing interactions. In the unfrustrated regime (J_{1}>2J_{2}), the model reproduces the standard clock model phenomenology: a low-temperature Z_{q}-broken ferromagnet, an intermediate XY-like critical quasi-long-range-ordered phase with emergent U(1) symmetry, and a high-temperature paramagnet. For J_{1}<2J_{2}, frustration stabilizes five distinct regimes: the disordered paramagnet, a stripe-ordered phase breaking Z_{q}×Z_{2} symmetry, two Z_{2}-broken nematic phases (one with and one without quasi-long-range order), and an exotic stripe phase with emergent discrete Z_{q} spin degrees of freedom prohibited in the microscopic Hamiltonian. Remarkably, this seemingly forbidden Z_{q} order emerges via a relevant operator in the infrared long-wavelength limit, rather than from an irrelevant perturbation, highlighting a nonstandard route to emergence. Using large-scale corner transfer matrix renormalization group calculations, complemented by classical Monte Carlo simulations, we map the complete phase diagram and identify Berezinskii-Kosterlitz-Thouless, Ising, first-order, and unconventional Landau-incompatible transitions between different phases. Finally, we propose an effective field-theoretic framework that encompasses these emergent orders and their interwoven transitions.

DOI
10.1103/6jpk-yj56

Related papers

Phase transitions in the q -state clock model

Arpita Goswami, Ravi Kumar, Monikana Gope +1 more

The q-state clock model, sometimes called the discrete XY model, is known to show a second-order (symmetry-breaking) phase transition in two dimensions for q≤4 (q=2 corresponds to the Ising model). On the other hand, the q→∞ limit of the model corresponds to the XY model, which shows the infinite order (non-symmetry-breaking) Berezinskii-Kosterlitz-Thouless (BKT) phase transition in two dimensions. Interestingly, the two-dimensional clock model with q≥5 is predicted to show three different phases and two associated phase transitions. There are varying opinions about the actual characters of phases and the associated transitions. In this work, we develop the basic and higher-order mean-field (MF) theories to study the q-state clock model systematically. Our MF calculations reaffirm that, for large q, there are three phases: (broken) Z_{q} symmetric ferromagnetic phase at the low temperature, emergent U(1) symmetric BKT phase at the intermediate temperature, and paramagnetic (disordered) phase at the high temperature. The phase transition at the higher temperature is found to be of the BKT type, and the other transition at the lower temperature is argued to be a large-order spontaneous symmetry-breaking type (the largeness of transition order yields the possibility of having some of the numerical characteristics of a BKT transition). The higher-order MF theory developed here better characterizes phases by estimating the spin-spin correlation between two neighbors.

arXiv:2407.17507•Physical review. E•May 19, 2025•5 citations

Phases and phase transitions of the disordered quantum clock model

Pulloor Kuttanikkad Vishnu, Gaurav Khairnar, Rajesh Narayanan +1 more

We study the phases and phase transitions of a disordered one-dimensional quantum $q$-state clock Hamiltonian using large-scale Monte Carlo simulations. Making contact with earlier results, we confirm that the clean, translational invariant version of the model, for $q=6$, hosts an intermediate emergent quasi-long-range ordered (QLRO) phase between the symmetry-broken true long-range ordered (TLRO) phase and the disordered (paramagnetic) phase. With increasing disorder strength, the quasi-long-range ordered phase shrinks and finally vanishes at a multicritical point, beyond which there is a direct transition from the TLRO phase to the paramagnetic phase. After establishing the phase diagram, we characterize the critical behaviors of the various quantum phase transitions in the model. We find that weak disorder is an irrelevant perturbation of the Berezinskii-Kosterlitz-Thouless transitions that separate the QLRO phase from the TLRO and paramagnetic phases. For stronger disorder, some of the critical exponents become disorder-dependent already before the system reaches the multicritical point. We also show that beyond the multicritical point, the direct transition from the TLRO phase to the paramagnetic phase is governed by an infinite-randomness critical point in line with strong-disorder renormalization group predictions. While our numerical results are for $q=6$, we expect the qualitative features of the behavior to hold for all $q>4$.

arXiv:2410.07640•Physical review. B./Physical review. B•Mar 28, 2025•3 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

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

Competing anisotropies and phase transitions in the $q$-state clock model with a $p$-fold crystalline field

Milan ˇ Zukovič

We study the two-dimensional $q$-state clock model in the presence of an additional $p$-fold symmetry-breaking crystalline field using Monte Carlo simulations. While the pure clock model exhibits Berezinskii--Kosterlitz--Thouless (BKT) transitions for sufficiently large $q$, the effect of competing discrete anisotropies on this topological phase remains nontrivial. We show that even weak crystalline fields qualitatively modify the phase diagram by suppressing the BKT phase and inducing transitions to states with true long-range order. The resulting behavior depends sensitively on the interplay between the intrinsic $\mathbb{Z}_q$ symmetry and the imposed $\mathbb{Z}_p$ anisotropy. In particular, in the six-state clock model for $p=2$ we observe qualitatively different scenarios depending on the sign of the field: a single transition for $h_2>0$ and a two-step ordering process for $h_2<0$ with an intermediate ordered phase. For $p=3$, the system exhibits a direct transition consistent with three-state Potts criticality. These results demonstrate that the phase structure cannot be inferred from symmetry considerations alone, but is governed by the competition between distinct locking mechanisms. Our findings provide a discrete counterpart to the multi-frequency sine-Gordon description of generalized $XY$ models and illustrate how additional anisotropies reshape topological phase transitions in two dimensions.

arXiv:2605.20990•May 20, 2026•0 citations

Phase diagrams of two-dimensional frustrated spin systems

Ansgar Kalz

In this work the ordering processes in classical frustrated Ising models and the ground-state behavior of frustrated quantum spin models are investigated. In the first part the phase transition from a high-temperature paramagnetic phase to a antiferromagnetic ground state is analyzed for the square-lattice Ising model with competing antiferromagnetic interactions between nearest- and next-nearest-neighbor spins. In the case of spatially isotropic interactions the phase transition to a collinear state is of special interest. Using Monte-Carlo simulations on finite lattices as well as field-theoretical methods for the continuous model it is shown that depending on the frustration three different sectors can be distinguished that show (i) Ising-universal behavior (ii) Ashkin-Teller-like non-universal behavior (varying critical exponents) and (iii) non-critical first-order behavior. For spatially anisotropic interactions the emergence of an intermediate incommensurate phase is verified by means of Monte-Carlo simulations which has been predicted by earlier mean-field calculations. The nature of this phase is studied analyzing spin-spin correlation functions and its ordering wave-vector reveals a strong temperature and size dependence. In the second part the influence of quantum fluctuations on classical antiferromagnetic ground states are investigated in anisotropic quantum spin-1/2 Heisenberg models with competing interactions. Using quantum Monte-Carlo simulations, exact diagonalization and perturbation theory the ground-state phase diagrams are calculated for a square lattice with nearest- and next-nearest-neighbor interactions and a honeycomb lattice with up to third-nearest-neighbor interactions whereby the fluctuations are chosen to be non-frustrating. For both cases the stability regions of classical magnetic states are identified and an intermediate phase is detected that shows no sign of long-range magnetic order. Further calculations also rule out valence-bond-crystal or topological order for this state, thus a quantum disordered ground state is stabilized in these models.

•Jan 1, 2012•0 citations
We present a comprehensive study on the frustrated J_{1}-J_{2} classical q-state clock model with even q>4 on a two-dimensional square lattice, revealing a rich ensemble of phases driven by competing interactions. In the unfrustrated regime (J_{1}>2J_{2}), the model reproduces the standard clock model phenomenology: a low-temperature Z_{q}-broken ferromagnet, an intermediate XY-like critical quasi-long-range-ordered phase with emergent U(1) symmetry, and a high-temperature paramagnet. For J_{1}<2J_{2}, frustration stabilizes five distinct regimes: the disordered paramagnet, a stripe-ordered phase breaking Z_{q}×Z_{2} symmetry, two Z_{2}-broken nematic phases (one with and one without quasi-long-range order), and an exotic stripe phase with emergent discrete Z_{q} spin degrees of freedom prohibited in the microscopic Hamiltonian. Remarkably, this seemingly forbidden Z_{q} order emerges via a relevant operator in the infrared long-wavelength limit, rather than from an irrelevant perturbation, highlighting a nonstandard route to emergence. Using large-scale corner transfer matrix renormalization group calculations, complemented by classical Monte Carlo simulations, we map the complete phase diagram and identify Berezinskii-Kosterlitz-Thouless, Ising, first-order, and unconventional Landau-incompatible transitions between different phases. Finally, we propose an effective field-theoretic framework that encompasses these emergent orders and their interwoven transitions.
The q-state clock model, sometimes called the discrete XY model, is known to show a second-order (symmetry-breaking) phase transition in two dimensions for q≤4 (q=2 corresponds to the Ising model). On the other hand, the q→∞ limit of the model corresponds to the XY model, which shows the infinite order (non-symmetry-breaking) Berezinskii-Kosterlitz-Thouless (BKT) phase transition in two dimensions. Interestingly, the two-dimensional clock model with q≥5 is predicted to show three different phases and two associated phase transitions. There are varying opinions about the actual characters of phases and the associated transitions. In this work, we develop the basic and higher-order mean-field (MF) theories to study the q-state clock model systematically. Our MF calculations reaffirm that, for large q, there are three phases: (broken) Z_{q} symmetric ferromagnetic phase at the low temperature, emergent U(1) symmetric BKT phase at the intermediate temperature, and paramagnetic (disordered) phase at the high temperature. The phase transition at the higher temperature is found to be of the BKT type, and the other transition at the lower temperature is argued to be a large-order spontaneous symmetry-breaking type (the largeness of transition order yields the possibility of having some of the numerical characteristics of a BKT transition). The higher-order MF theory developed here better characterizes phases by estimating the spin-spin correlation between two neighbors.
We study the phases and phase transitions of a disordered one-dimensional quantum qqq-state clock Hamiltonian using large-scale Monte Carlo simulations. Making contact with earlier results, we confirm that the clean, translational invariant version of the model, for q=6q=6q=6, hosts an intermediate emergent quasi-long-range ordered (QLRO) phase between the symmetry-broken true long-range ordered (TLRO) phase and the disordered (paramagnetic) phase. With increasing disorder strength, the quasi-long-range ordered phase shrinks and finally vanishes at a multicritical point, beyond which there is a direct transition from the TLRO phase to the paramagnetic phase. After establishing the phase diagram, we characterize the critical behaviors of the various quantum phase transitions in the model. We find that weak disorder is an irrelevant perturbation of the Berezinskii-Kosterlitz-Thouless transitions that separate the QLRO phase from the TLRO and paramagnetic phases. For stronger disorder, some of the critical exponents become disorder-dependent already before the system reaches the multicritical point. We also show that beyond the multicritical point, the direct transition from the TLRO phase to the paramagnetic phase is governed by an infinite-randomness critical point in line with strong-disorder renormalization group predictions. While our numerical results are for q=6q=6q=6, we expect the qualitative features of the behavior to hold for all q>4q>4q>4.
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 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.
Competing anisotropies and phase transitions in the qqq-state clock model with a ppp-fold crystalline field
We study the two-dimensional qqq-state clock model in the presence of an additional ppp-fold symmetry-breaking crystalline field using Monte Carlo simulations. While the pure clock model exhibits Berezinskii--Kosterlitz--Thouless (BKT) transitions for sufficiently large qqq, the effect of competing discrete anisotropies on this topological phase remains nontrivial. We show that even weak crystalline fields qualitatively modify the phase diagram by suppressing the BKT phase and inducing transitions to states with true long-range order. The resulting behavior depends sensitively on the interplay between the intrinsic Zq\mathbb{Z}_qZq​ symmetry and the imposed Zp\mathbb{Z}_pZp​ anisotropy. In particular, in the six-state clock model for p=2p=2p=2 we observe qualitatively different scenarios depending on the sign of the field: a single transition for h2>0h_2>0h2​>0 and a two-step ordering process for h2<0h_2<0h2​<0 with an intermediate ordered phase. For p=3p=3p=3, the system exhibits a direct transition consistent with three-state Potts criticality. These results demonstrate that the phase structure cannot be inferred from symmetry considerations alone, but is governed by the competition between distinct locking mechanisms. Our findings provide a discrete counterpart to the multi-frequency sine-Gordon description of generalized XYXYXY models and illustrate how additional anisotropies reshape topological phase transitions in two dimensions.
In this work the ordering processes in classical frustrated Ising models and the ground-state behavior of frustrated quantum spin models are investigated. In the first part the phase transition from a high-temperature paramagnetic phase to a antiferromagnetic ground state is analyzed for the square-lattice Ising model with competing antiferromagnetic interactions between nearest- and next-nearest-neighbor spins. In the case of spatially isotropic interactions the phase transition to a collinear state is of special interest. Using Monte-Carlo simulations on finite lattices as well as field-theoretical methods for the continuous model it is shown that depending on the frustration three different sectors can be distinguished that show (i) Ising-universal behavior (ii) Ashkin-Teller-like non-universal behavior (varying critical exponents) and (iii) non-critical first-order behavior. For spatially anisotropic interactions the emergence of an intermediate incommensurate phase is verified by means of Monte-Carlo simulations which has been predicted by earlier mean-field calculations. The nature of this phase is studied analyzing spin-spin correlation functions and its ordering wave-vector reveals a strong temperature and size dependence. In the second part the influence of quantum fluctuations on classical antiferromagnetic ground states are investigated in anisotropic quantum spin-1/2 Heisenberg models with competing interactions. Using quantum Monte-Carlo simulations, exact diagonalization and perturbation theory the ground-state phase diagrams are calculated for a square lattice with nearest- and next-nearest-neighbor interactions and a honeycomb lattice with up to third-nearest-neighbor interactions whereby the fluctuations are chosen to be non-frustrating. For both cases the stability regions of classical magnetic states are identified and an intermediate phase is detected that shows no sign of long-range magnetic order. Further calculations also rule out valence-bond-crystal or topological order for this state, thus a quantum disordered ground state is stabilized in these models.