arXiv:0805.3107 · Physical Review B
Effect of anisotropy on the ground-state magnetic ordering of the spin-half quantum J 1 X X Z – J 2 X X Z model on the square lattice
R. F. Bishop, P. H. Y. Li, R. Darradi, J. Schulenburg, Johannes Richter
Abstract
We study the zero-temperature phase diagram of the two-dimensional quantum ${J}_{1}^{XXZ}--{J}_{2}^{XXZ}$ spin-1/2 anisotropic Heisenberg model on the square lattice. In particular, the effects of the anisotropy $\ensuremath{\Delta}$ on the $z$-aligned N\'eel and (collinear) stripe states, as well as on the $xy$-planar-aligned N\'eel and collinear stripe states, are examined. All four of these quasiclassical states are chosen in turn as model states, on top of which we systematically include the quantum correlations using a coupled cluster method analysis carried out to very high orders. We find strong evidence for two quantum triple points (QTPs) at $({\ensuremath{\Delta}}^{c}=\ensuremath{-}0.10\ifmmode\pm\else\textpm\fi{}0.15,\text{ }{J}_{2}^{c}/{J}_{1}=0.505\ifmmode\pm\else\textpm\fi{}0.015)$ and $({\ensuremath{\Delta}}^{c}=2.05\ifmmode\pm\else\textpm\fi{}0.15,\text{ }{J}_{2}^{c}/{J}_{1}=0.530\ifmmode\pm\else\textpm\fi{}0.015)$, between which an intermediate magnetically disordered phase emerges to separate the quasiclassical N\'eel and stripe collinear phases. Above the upper QTP $(\ensuremath{\Delta}\ensuremath{\gtrsim}2.0)$ we find a direct first-order phase transition between the N\'eel and stripe phases, exactly as for the classical case. The $z$-aligned and $xy$-planar-aligned phases meet precisely at $\ensuremath{\Delta}=1$, also as for the classical case. For all values of the anisotropy parameter between those of the two QTPs there exists a narrow range of values of ${J}_{2}/{J}_{1}$, ${\ensuremath{\alpha}}^{{c}_{1}}(\ensuremath{\Delta})<{J}_{2}/{J}_{1}<{\ensuremath{\alpha}}^{{c}_{2}}(\ensuremath{\Delta})$, centered near the point of maximum classical frustration, ${J}_{2}/{J}_{1}=\frac{1}{2}$, for which the intermediate phase exists. This range is widest precisely at the isotropic point, $\ensuremath{\Delta}=1$, where ${\ensuremath{\alpha}}^{{c}_{1}}(1)=0.44\ifmmode\pm\else\textpm\fi{}0.01$ and ${\ensuremath{\alpha}}^{{c}_{2}}(1)=0.59\ifmmode\pm\else\textpm\fi{}0.01$. The two QTPs are characterized by values $\ensuremath{\Delta}={\ensuremath{\Delta}}^{c}$ at which ${\ensuremath{\alpha}}^{{c}_{1}}({\ensuremath{\Delta}}^{c})={\ensuremath{\alpha}}^{{c}_{2}}({\ensuremath{\Delta}}^{c})$.