We study a generalization of the well-known classical two-dimensional square lattice compass model of XY spins (sometimes referred to as the ${90}^{\ensuremath{\circ}}$ compass model), which interpolates between the XY model and the compass model. Our model possesses the combined
C4 lattice and spin rotation symmetry of the compass model but is free of its fine-tuned subsystem symmetries. Using both field theoretic arguments and Monte Carlo simulations, we find that our model possesses a line of critical points with continuously varying exponents of the Ashkin-Teller type terminating at the four-state Potts point. Further, our Monte Carlo study uncovers that beyond the four-state Potts point, the line of phase transition is connected to the lattice-nematic Ising phase transition in the square lattice compass model through a region of first-order transitions.
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 (nearest neighbor, ferromagnetic) and
J2>0 (second neighbor, antiferromagnetic). The striped phase breaks a
Z4 symmetry and is stabilized at low temperatures for
g=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-
J2 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
J1-
J2 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.
We study the thermal phase transition of the fourfold degenerate phases (the plaquette and single-stripe states) in the two-dimensional frustrated Ising model on the Shastry-Sutherland lattice using Monte Carlo simulations. The critical Ashkin-Teller-like behavior is identified both in the plaquette phase region and the single-stripe phase region. The four-state Potts critical end points differentiating the continuous transitions from the first-order ones are estimated based on finite-size-scaling analyses. Furthermore, a similar behavior of the transition to the fourfold single-stripe phase is also observed in the anisotropic triangular Ising model. Thus, this work clearly demonstrates that the transitions to the fourfold degenerate states of two-dimensional Ising antiferromagnets exhibit similar transition behavior.
We investigate thermal phase transitions to a valence-bond solid phase in
SU(N) Heisenberg models with four- or six-body interactions on a square or honeycomb lattice, respectively. In both cases, a thermal phase transition occurs that is accompanied by rotational symmetry breaking of the lattice. We perform quantum Monte Carlo calculations in order to clarify the critical properties of the models. The estimated critical exponents indicate that the universality classes of the square- and honeycomb-lattice cases are identical to those of the classical
XY model with a
Z4 symmetry-breaking field and the three-state Potts model, respectively. In the square-lattice case, the thermal exponent, $\ensuremath{\nu}$, monotonically increases as the system approaches the quantum critical point, while the values of the critical exponents, $\ensuremath{\eta}$ and $\ensuremath{\gamma}/\ensuremath{\nu}$, remain constant. From a finite-size scaling analysis, we find that the system exhibits weak universality, because the
Z4 symmetry-breaking field is always marginal. In contrast, $\ensuremath{\nu}$ in the honeycomb-lattice case exhibits a constant value, even in the vicinity of the quantum critical point, because the
Z3 field remains relevant in the SU(3) and SU(4) cases.
We investigate the classical chiral Ashkin-Teller model on a square lattice with the corner transfer matrix renormalization group algorithm. We show that the melting of the period-4 phase in the presence of a chiral perturbation takes different forms depending on the coefficient of the four-spin term in the Ashkin-Teller model. Close to the clock limit of two decoupled Ising models, the system undergoes a two-step commensurate-incommensurate transition as soon as the chirality is introduced, with an intermediate critical floating phase bounded by a Kosterlitz-Thouless transition at high temperature and a Pokrovsky-Talapov transition at low temperature. By contrast, close to the four-states Potts model, we argue for the existence of a unique commensurate-incommensurate transition that belongs to the chiral universality class, and for the presence of a Lifshitz point where the ordered, disordered, and floating phases meet. Finally, we map the whole phase diagram, which turns out to be in qualitative agreement with the 40-year-old prediction by Huse and Fisher.
A comprehensive study of the two-dimensional (2D) compass model on the square lattice is performed for classical and quantum spin degrees of freedom using Monte Carlo and quantum Monte Carlo methods. We employ state-of-the-art implementations using Metropolis, stochastic series expansion, and parallel tempering techniques to obtain the critical ordering temperatures and critical exponents. In a preinvestigation we reconsider the classical compass model where we study and contrast the finite-size scaling behavior of ordinary periodic boundary conditions against annealed boundary conditions. It is shown that periodic boundary conditions suffer from extreme finite-size effects which might be caused by closed-loop excitations on the torus. These excitations also appear to have severe effects on the Binder parameter. On this footing we report on a systematic Monte Carlo study of the quantum compass model. Our numerical results are at odds with recent literature on the subject which we trace back to neglecting the strong finite-size effects on periodic lattices. The critical temperatures are obtained as
Tc=0.1464(2)J and
Tc=0.055(1)J for the classical and quantum versions, respectively, and our data support a transition in the 2D Ising universality class for both cases.