We study the zero-temperature ground-state (GS) phase diagram of a spin-half
J1--
J2 XY model on the honeycomb lattice with nearest-neighbor exchange coupling
J1>0 and frustrating next-nearest-neighbor exchange coupling ${J}_{2}\ensuremath{\equiv}\ensuremath{\kappa}{J}_{1}>0$, where both bonds are of the isotropic
XY type, using the coupled cluster method. Results are presented for the GS energy per spin, magnetic order parameter, and staggered dimer valence-bond crystalline (SDVBC) susceptibility, for values of the frustration parameter in the range $0\ensuremath{\le}\ensuremath{\kappa}\ensuremath{\le}1$. In this range, we find phases exhibiting, respectively, N\'eel
xy planar [N(p)], N\'eel
z-aligned [N(
z)], SDVBC, and N\'eel-II
xy planar [N-II(p)] orderings. The N\'eel-II states, which break the lattice rotational symmetry, are ones in which the spins of nearest-neighbor pairs along one of the three equivalent honeycomb directions are parallel, while those in the other two directions are antiparallel. The N(p) state, which is stable for the classical version of the model in the range $0\ensuremath{\le}\ensuremath{\kappa}\ensuremath{\le}\frac{1}{6}$, is found to form the GS phase out to a first quantum critical point at ${\ensuremath{\kappa}}_{{c}_{1}}=0.216(5)$, beyond which the stable GS phase has N(
z) order over the range ${\ensuremath{\kappa}}_{{c}_{1}}<\ensuremath{\kappa}<{\ensuremath{\kappa}}_{{c}_{2}}=0.355(5)$. For values $\ensuremath{\kappa}>{\ensuremath{\kappa}}_{{c}_{2}}$, we find a strong competition to form the GS phase between states with N-II(p) and SDVBC forms of order. Our best estimate, however, is that the stable GS phase over the range ${\ensuremath{\kappa}}_{{c}_{2}}<\ensuremath{\kappa}<{\ensuremath{\kappa}}_{{c}_{3}}\ensuremath{\approx}0.52(3)$ is a mixed state with both SDVBC and N-II(p) forms of order; and for values $\ensuremath{\kappa}>{\ensuremath{\kappa}}_{{c}_{3}}$ is the N-II(p) state, which is stable at the classical level only at the highly degenerate point $\ensuremath{\kappa}=\frac{1}{2}$. Over the range $0\ensuremath{\le}\ensuremath{\kappa}\ensuremath{\le}1$, we find no evidence for any of the spiral phases that are present classically for all values $\ensuremath{\kappa}>\frac{1}{6}$, nor for any quantum spin-liquid state.
Using the density matrix renormalization group on wide cylinders, we study the phase diagram of the spin-1/2 XY model on the honeycomb lattice, with first-neighbor (J1=1) and frustrating second-neighbor (J2>0) interactions. For the intermediate frustration regime 0.22≲J2≲0.36, we find a surprising antiferromagnetic Ising phase, with ordered moments pointing along the z axis, despite the absence of any S(z)S(z) interactions in the Hamiltonian. Surrounding this phase as a function of J2 are antiferromagnetic phases with the moments pointing in the x-y plane for small J2 and a close competition between an x-y plane magnetic collinear phase and a dimer phase for large values of J2. We do not find any spin-liquid phases in this model.
Searching for spin-liquid states has long been attracting both experimentalists and theorists. In this paper, we review recent density matrix renormalization group studies of the spin-½ XY model on the honeycomb lattice, with first-neighbor (J1 = 1) and frustrating second-neighbor (J2 > 0) interactions. For the intermediate frustration regime 0.22 ≲ J2 ≲ 0.36, there exists a surprising antiferromagnetic Ising phase, with ordered moments pointing along the z-axis, despite the absence of any SzSz interactions in the Hamiltonian. Surrounding this phase as a function of J2 are antiferromagnetic phases with the moments pointing in the xy-plane for small J2 and a close competition between an xy-plane magnetic collinear phase and a dimer phase for large values of J2. No spin-liquid phases was found in the XY model even with the third-neighbor (J3 > 0) interactions.
We use the coupled cluster method (CCM) to investigate the ground-state (GS) properties of the frustrated spin-
21 J1-
J2-
J3 model on the honeycomb lattice, with nearest-neighbor exchange coupling
J1 plus next-nearest-neighbor (
J2) and next-next-nearest-neighbor (
J3) exchanges of equal strength. In particular we find a direct first-order phase transition between the N\'eel-ordered antiferromagnetic phase and the ferromagnetic phase at a value ${J}_{2}/{J}_{1}=\ensuremath{-}1.17\ifmmode\pm\else\textpm\fi{}0.01$ when
J1>0, compared to the corresponding classical value of $\ensuremath{-}1$. We find no evidence for any intermediate phase. From this and our previous CCM studies of the model we present its full zero-temperature GS phase diagram.
We use the coupled-cluster method in high orders of approximation to make a comprehensive study of the ground-state (GS) phase diagram of the spin-1/2
J1-
J2-
J3 model on a two-dimensional honeycomb lattice with antiferromagnetic (AFM) interactions up to third-nearest neighbors. Results are presented for the GS energy and the average local onsite magnetization. With the nearest-neighbor coupling strength ${J}_{1}\ensuremath{\equiv}1$, we find four magnetically ordered phases in the parameter window ${J}_{2},{J}_{3}\ensuremath{\in}[0,1]$, namely, the N\'eel, striped, and N\'eel-II collinear AFM phases, plus a spiral phase. The N\'eel-II phase appears as a stable GS phase in the classical version of the model only for values
J3<0. Each of these four ordered phases shares a boundary with a disordered quantum paramagnetic (QP) phase, and at several widely separated points on the phase boundaries the QP phase has an infinite susceptibility to plaquette valence-bond crystalline order. We identify all of the phase boundaries with good precision in the parameter window studied, and we find three tricritical quantum critical points therein at (a) $({J}_{2}^{{c}_{1}},{J}_{3}^{{c}_{1}})=(0.51\ifmmode\pm\else\textpm\fi{}0.01,0.69\ifmmode\pm\else\textpm\fi{}0.01)$ between the N\'eel, striped, and QP phases; (b) $({J}_{2}^{{c}_{2}},{J}_{3}^{{c}_{2}})=(0.65\ifmmode\pm\else\textpm\fi{}0.02,0.55\ifmmode\pm\else\textpm\fi{}0.01)$ between the striped, spiral, and QP phases; and (c) $({J}_{2}^{{c}_{3}},{J}_{3}^{{c}_{3}})=(0.69\ifmmode\pm\else\textpm\fi{}0.01,0.12\ifmmode\pm\else\textpm\fi{}0.01)$ between the spiral, N\'eel-II, and QP phases.
We study the ground-state (GS) properties of the frustrated spin-
21 J1--
J2--
J3 Heisenberg model on the two-dimensional honeycomb lattice with ferromagnetic nearest-neighbor (${J}_{1}=\ensuremath{-}1$) exchange and frustrating antiferromagnetic next-nearest-neighbor (
J2>0) and next-next-nearest-neighbor (
J3>0) exchanges, for the case
J3=J2. We use the coupled-cluster method implemented to high orders of approximation, complemented by the Lanczos exact diagonalization of a large finite lattice with 32 sites, in order to calculate the GS energy, magnetic order parameter, and spin-spin correlation functions. In one scenario we find a quantum phase transition point between regions characterized by ferromagnetic order and a form of antiferromagnetic (``striped'') collinear order at ${J}_{2}^{c}\ensuremath{\approx}0.1095\ifmmode\pm\else\textpm\fi{}0.0005$, which is below the corresponding hypothetical transition point at
J2cl=71 ($\ensuremath{\approx}0.143$) for the classical version of the model, in which we momentarily ignore the intervening noncollinear spiral phase in the region
101<J2<51. Hence we see that quantum fluctuations appear to stabilize somewhat the collinear antiferromagnetic order in preference to the ferromagnetic order in this model. We compare results for the present ferromagnetic case (with ${J}_{1}=\ensuremath{-}1$) to previous results for the corresponding antiferromagnetic case (with
J1=+1). The magnetic order parameter is found to behave similarly for the ferromagnetic and the antiferromagnetic models for large values of the frustration parameter
J2. However, there are considerable differences in the behavior of the order parameters for the two models for ${J}_{2}/|{J}_{1}|\ensuremath{\lesssim}0.6$. For example, the quasiclassical collinear magnetic long-range order for the antiferromagnetic model (with
J1=+1) breaks down at ${J}_{2}^{{c}_{2}}\ensuremath{\approx}0.60$, whereas the ``equivalent'' point for the ferromagnetic model (with ${J}_{1}=\ensuremath{-}1$) occurs at ${J}_{2}^{c}\ensuremath{\approx}0.11$. Unlike in the antiferromagnetic model (with
J1=+1), where a plaquette valence-bond crystal phase intrudes between the two corresponding quasiclassical antiferromagnetic phases (with N\'eel and striped order) for
J2c1<J2<J2c2, with ${J}_{2}^{{c}_{1}}\ensuremath{\approx}0.47$, we find no clear indications at all in the ferromagnetic model for an intermediate magnetically disordered phase between the corresponding phases exhibiting ferromagnetic and striped order. Instead the evidence for the ferromagnetic model (with ${J}_{1}=\ensuremath{-}1$) points to one of two scenarios: either there is a direct first-order transition between the two magnetically ordered phases, as mentioned above; or there exists an intervening phase between them in the very narrow range $0.10\ensuremath{\lesssim}{J}_{2}\ensuremath{\lesssim}0.12$, which is probably a remnant of the spiral phase that exists in the classical counterpart of the model over the larger range
101<J2<51.