Phase transitions in the frustrated Ising model on the square lattice

TL;DR

This paper investigates the nature of phase transitions in the frustrated J1-J2 Ising model on a square lattice, revealing a line of weak first-order transitions and a continuous transition line mapped to the Ashkin-Teller model.

Contribution

It provides a detailed analysis of the phase transition order and universality class in the J1-J2 Ising model, clarifying the transition from weak first-order to continuous behavior.

Findings

Line of weak first-order transitions for 1/2<g<g*

Continuous transitions mapped to Ashkin-Teller model

Pseudo-first-order effects near the multicritical point g*

Abstract

We consider the thermal phase transition from a paramagnetic to 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. [PRL 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(1). Thereafter, the transitions from g above g* are continuous and can be fully mapped, using universality arguments, to the critical line of the well known Ashkin-Teller model from its 4-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 to be 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 results are based on Monte Carlo calculations, while the cluster mean-field and transfer-matrix results provide useful methodological bench-marks for weakly first-order behaviors and Ashkin-Teller criticality.