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.
Using the coupled cluster method for high orders of approximation and complementary exact diagonalization studies we investigate the ground state properties of the spin-1/2
J1−−J2 frustrated Heisenberg antiferromagnet on the square lattice. We have calculated the ground state energy, the magnetic order parameter, the spin stiffness, and several generalized susceptibilities to probe magnetically disordered quantum valence-bond phases. We have found that the quantum critical points for both the N\'eel and collinear orders are ${J}_{2}^{c1}\ensuremath{\approx}(0.44\ifmmode\pm\else\textpm\fi{}0.01){J}_{1}$ and ${J}_{2}^{c2}\ensuremath{\approx}(0.59\ifmmode\pm\else\textpm\fi{}0.01){J}_{1}$, respectively, which are in good agreement with the results obtained by other approximations. In contrast to the recent study by [Sirker et al. Phys. Rev. B 73, 184420 (2006)], our data do not provide evidence for the transition from the N\'eel to the valence-bond solid state to be first order. Moreover, our results are in favor of the deconfinement scenario for that phase transition. We also discuss the nature of the magnetically disordered quantum phase.
The zero-temperature quantum phase diagram of the spin-
21 J1--
J2--
J1⊥ model on an
AA-stacked bilayer honeycomb lattice is investigated using the coupled cluster method (CCM). The model comprises two monolayers in each of which the spins, residing on honeycomb-lattice sites, interact via both nearest-neighbor (NN) and frustrating next-nearest-neighbor isotropic antiferromagnetic (AFM) Heisenberg exchange iteractions, with respective strengths
J1>0 and
J2≡κJ1>0. The two layers are coupled via a comparable Heisenberg exchange interaction between NN interlayer pairs, with a strength
J1⊥≡δJ1. The complete phase boundaries of two quasiclassical collinear AFM phases, namely the Neel and Neel-II phases, are calculated in the
κδ half-plane with
κ>0. Whereas on each monolayer in the Neel state all NN pairs of spins are antiparallel, in the Neel-II state NN pairs of spins on zigzag chains along one of the three equivalent honeycomb-lattice directions are antiparallel, while NN interchain spins are parallel. We calculate directly in the thermodynamic (infinite-lattice) limit both the magnetic order parameter
M and the excitation energy
Δ from the
sTz=0 ground state to the lowest-lying
∣sTz∣=1 excited state (where
sTz is the total
z component of spin for the system as a whole, and where the collinear ordering lies along the
z direction) for both quasiclassical states used (separately) as the CCM model state, on top of which the multispin quantum correlations are then calculated to high orders (
n≤10) in a systematic series of approximations involving
n-spin clusters. The sole approximation made is then to extrapolate the sequences of
nth-order results for
M and
Δ to the exact limit,
n→∞.
We report on recent results for strongly frustrated quantum J1 - J2 antiferromagnets in dimensionality d = 1, 2, 3 obtained by the coupled cluster method (CCM). We demonstrate that the CCM in high orders of approximation allows us to investigate quantum phase transitions driven by frustration and to discuss novel quantum ground states. In detail we consider the ground-state properties of (i) the Heisenberg spin-1/2 antiferromagnet on the cubic lattice in d = 1, 2, 3, and use the results for the energy, the sublattice magnetization and the spin stiffness as a benchmark test for the precision of the method; (ii) coupled frustrated spin chains (the quasi-one-dimensional J1 - J2 model) and discuss the influence of the quantum fluctuations and the interchain coupling on the incommensurate spiral state present in the classical model; (iii) the Shastry-Sutherland antiferromagnet on the square lattice; and (iv) a stacked frustrated square-lattice Heisenberg antiferromagnet (the quasi-two-dimensional J1 - J2 model), and discuss the influence of the interlayer coupling on the quantum paramagnetic ground-state phase that is present for the strictly two-dimensional model.
We use a combination of analytical and numerical techniques to study the phase diagram of the frustrated Heisenberg model on the bilayer honeycomb lattice. Using the Schwinger-boson description of the spin operators followed by a mean-field decoupling, the magnetic phase diagram is studied as a function of the frustration coupling
J2 and the interlayer coupling ${J}_{\ensuremath{\perp}}$. The presence of both magnetically ordered and disordered phases is investigated by means of the evaluation of ground-state energy, spin gap, local magnetization, and spin-spin correlations. We observe a phase with a spin gap and short-range N\'eel correlations that survives for nonzero next-nearest-neighbor interaction and interlayer coupling. Furthermore, we detect signatures of a reentrant behavior in the melting of the N\'eel phase and symmetry restoring when the system undergoes a transition from an on-layer nematic valence-bond crystal phase to an interlayer valence-bond crystal phase. We complement our work with exact diagonalization on small clusters and dimer-series expansion calculations, together with a linear spin-wave approach to study the phase diagram as a function of the spin
S, the frustration, and the interlayer couplings.
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.