We use the coupled cluster method (CCM) to study the zero-temperature phase diagram of a two-dimensional frustrated spin-half antiferromagnet, the so-called Union Jack model. It is defined on a square lattice such that all nearest-neighbor pairs are connected by bonds with a strength J 1>0, but only half the next-nearest-neighbor pairs are connected by bonds with a strength J2 ≡κJ1>0. The bonds are arranged such that on the 2×2 unit cell they form the pattern of the Union Jack flag. Alternating sites on the square lattice are thus four-connected and eight-connected. We find strong evidence for a first phase transition between a Néel antiferromagnetic phase and a canted ferrimagnetic phase at a critical coupling κc1 =0.66±0.02. The transition is an interesting one, at which the energy and its first derivative seem continuous, thus providing a typical scenario of a second-order transition (just as in the classical case for the model), although a weakly first-order transition cannot be excluded. By contrast, the average on-site magnetization approaches a nonzero value Mc1 =0.195±0.005 on both sides of the transition, which is more typical of a first-order transition. The slope, dM/dκ, of the order parameter curve as a function of the coupling strength κ, also appears to be continuous, or very nearly so, at the critical point κc1, thereby providing further evidence of the subtle nature of the transition between the Néel and canted phases. Our CCM calculations provide strong evidence that the canted ferrimagnetic phase becomes unstable at large values of κ, and hence we have also used the CCM with a model collinear semistripe-ordered ferrimagnetic state in which alternating rows (and columns) are ferromagnetically and antiferromagnetically ordered, and in which the spins connected by J2 bonds are antiparallel to one another. We find tentative evidence, based on the relative energies of the two states, for a second zero-temperature phase transition between the canted and semistripe-ordered ferrimagnetic states at a large value of the coupling parameter around κc2 ≈ 125±5. This prediction, however, is based on an extrapolation of the CCM results for the canted state into regimes where the solutions have already become unstable and the CCM equations based on the canted state at any level of approximation beyond the lowest have no solutions. Our prediction for κ c2 is hence less reliable than that for κ c1. Nevertheless, if this second transition at κ c2 does exist, our results clearly indicate it to be of first-order type.
|Original language||English (US)|
|Journal||Physical Review B - Condensed Matter and Materials Physics|
|State||Published - Jul 15 2010|