We use the coupled cluster method (CCM) to study a frustrated spin-12J1−J2−J1⊥ Heisenberg antiferromagnet on a bilayer honeycomb lattice with AA stacking. Both nearest-neighbor (NN) and frustrating next-nearest-neighbor antiferromagnetic (AFM) exchange interactions are present in each layer, with respective exchange coupling constants J1 > 0 and J2 ≡ κJ1 > 0. The two layers are coupled with NN AFM exchanges with coupling strength J1⊥≡δJ1>0. We calculate to high orders of approximation within the CCM the zero-field transverse magnetic susceptibility χ in the Néel phase. We thus obtain an accurate estimate of the full boundary of the Néel phase in the κδ plane for the zero-temperature quantum phase diagram. We demonstrate explicitly that the phase boundary derived from χ is fully consistent with that obtained from the vanishing of the Néel magnetic order parameter. We thus conclude that at all points along the Néel phase boundary quasiclassical magnetic order gives way to a nonclassical paramagnetic phase with a nonzero energy gap. The Néel phase boundary exhibits a marked reentrant behavior, which we discuss in detail.
We use the coupled cluster method (CCM) to study a frustrated spin-12J1−J2−J1⊥ Heisenberg antiferromagnet on a bilayer honeycomb lattice with AA stacking. Both nearest-neighbor (NN) and frustrating next-nearest-neighbor antiferromagnetic (AFM) exchange interactions are present in each layer, with respective exchange coupling constants J1 > 0 and J2 ≡ κJ1 > 0. The two layers are coupled with NN AFM exchanges with coupling strength J1⊥≡δJ1>0. We calculate to high orders of approximation within the CCM the zero-field transverse magnetic susceptibility χ in the Néel phase. We thus obtain an accurate estimate of the full boundary of the Néel phase in the κδ plane for the zero-temperature quantum phase diagram. We demonstrate explicitly that the phase boundary derived from χ is fully consistent with that obtained from the vanishing of the Néel magnetic order parameter. We thus conclude that at all points along the Néel phase boundary quasiclassical magnetic order gives way to a nonclassical paramagnetic phase with a nonzero energy gap. The Néel phase boundary exhibits a marked reentrant behavior, which we discuss in detail.
Quantum phase diagram of a frustrated spin-21J1--J2--J1⊥ Heisenberg model on an AA-stacked bilayer honeycomb lattice
The zero-temperature quantum phase diagram of the spin-21J1--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 study a frustrated spin-$\frac{1}{2}\phantom{\rule{0.16em}{0ex}}{J}_{1}\text{\ensuremath{-}}{J}_{2}\text{\ensuremath{-}}{J}_{3}\text{\ensuremath{-}}{J}_{1}^{\ensuremath{\perp}}$ Heisenberg antiferromagnet on an AA-stacked bilayer honeycomb lattice. In each layer we consider nearest-neighbor (NN), next-nearest-neighbor, and next-next-nearest-neighbor antiferromagnetic (AFM) exchange couplings J1,J2, and J3, respectively. The two layers are coupled with an AFM NN exchange coupling ${J}_{1}^{\ensuremath{\perp}}\ensuremath{\equiv}\ensuremath{\delta}{J}_{1}$. The model is studied for arbitrary values of $\ensuremath{\delta}$ along the line ${J}_{3}={J}_{2}\ensuremath{\equiv}\ensuremath{\alpha}{J}_{1}$ that includes the most highly frustrated point at $\ensuremath{\alpha}=\frac{1}{2}$, where the classical ground state is macroscopically degenerate. The coupled cluster method is used at high orders of approximation to calculate the magnetic order parameter and the triplet spin gap. We are thereby able to give an accurate description of the quantum phase diagram of the model in the $\ensuremath{\alpha}\ensuremath{\delta}$ plane in the window $0\ensuremath{\le}\ensuremath{\alpha}\ensuremath{\le}1,\phantom{\rule{0.16em}{0ex}}0\ensuremath{\le}\ensuremath{\delta}\ensuremath{\le}1$. This includes two AFM phases with N\'eel and striped order, and an intermediate gapped paramagnetic phase that exhibits various forms of valence-bond crystalline order. We obtain accurate estimations of the two phase boundaries, $\ensuremath{\delta}={\ensuremath{\delta}}_{{c}_{i}}(\ensuremath{\alpha})$, or equivalently, $\ensuremath{\alpha}={\ensuremath{\alpha}}_{{c}_{i}}(\ensuremath{\delta})$, with i=1 (N\'eel) and 2 (striped). The two boundaries exhibit an ``avoided crossing'' behavior with both curves being re-entrant. Thus, in this $\ensuremath{\alpha}\ensuremath{\delta}$ window, N\'eel order exists only for values of $\ensuremath{\delta}$ in the range ${\ensuremath{\delta}}_{{c}_{1}}^{<}(\ensuremath{\alpha})<\ensuremath{\delta}<{\ensuremath{\delta}}_{{c}_{1}}^{>}(\ensuremath{\alpha})$, with ${\ensuremath{\delta}}_{{c}_{1}}^{<}(\ensuremath{\alpha})=0$ for $\ensuremath{\alpha}<{\ensuremath{\alpha}}_{{c}_{1}}(0)\ensuremath{\approx}0.46(2)$ and ${\ensuremath{\delta}}_{{c}_{1}}^{<}(\ensuremath{\alpha})>0$ for ${\ensuremath{\alpha}}_{{c}_{1}}(0)<\ensuremath{\alpha}<{\ensuremath{\alpha}}_{1}^{>}\ensuremath{\approx}0.49(1)$, and striped order similarly exists only for values of $\ensuremath{\delta}$ in the range ${\ensuremath{\delta}}_{{c}_{2}}^{<}(\ensuremath{\alpha})<\ensuremath{\delta}<{\ensuremath{\delta}}_{{c}_{2}}^{>}(\ensuremath{\alpha})$, with ${\ensuremath{\delta}}_{{c}_{2}}^{<}(\ensuremath{\alpha})=0$ for $\ensuremath{\alpha}>{\ensuremath{\alpha}}_{{c}_{2}}(0)\ensuremath{\approx}0.600(5)$ and ${\ensuremath{\delta}}_{{c}_{2}}^{<}(\ensuremath{\alpha})>0$ for ${\ensuremath{\alpha}}_{{c}_{2}}(0)>\ensuremath{\alpha}>{\ensuremath{\alpha}}_{2}^{<}\ensuremath{\approx}0.56(1)$.
We use the coupled cluster method to study the zero-temperature quantum phase diagram of the spin-$\frac{1}{2}{J}_{1}\text{\ensuremath{-}}{J}_{2}\text{\ensuremath{-}}{J}_{1}^{\ensuremath{\perp}}$ model on the honeycomb bilayer lattice. In each layer, we include both nearest-neighbor and frustrating next-nearest-neighbor antiferromagnetic exchange couplings, of strength J1>0 and ${J}_{2}\ensuremath{\equiv}\ensuremath{\kappa}{J}_{1}>0$, respectively. The two layers are coupled by an interlayer nearest-neighbor exchange, with coupling constant ${J}_{1}^{\ensuremath{\perp}}\ensuremath{\equiv}\ensuremath{\delta}{J}_{1}>0$. We calculate directly in the infinite-lattice limit both the ground-state energy per spin and the N\'eel magnetic order parameter, as well as the triplet spin gap. By implementing the method to very high orders of approximation we obtain an accurate estimate for the full boundary of the N\'eel phase in the $\ensuremath{\kappa}\ensuremath{\delta}$ plane. For each value $\ensuremath{\delta}<{\ensuremath{\delta}}_{c}^{>}(0)\ensuremath{\approx}1.70(5)$, we find an upper critical value ${\ensuremath{\kappa}}_{c}(\ensuremath{\delta})$, such that N\'eel order is present for $\ensuremath{\kappa}<{\ensuremath{\kappa}}_{c}(\ensuremath{\delta})$. Conversely, for each value $\ensuremath{\kappa}<{\ensuremath{\kappa}}_{c}(0)\ensuremath{\approx}0.19(1)$, we find an upper critical value ${\ensuremath{\delta}}_{c}^{>}(\ensuremath{\kappa})$, such that N\'eel order persists for $0<\ensuremath{\delta}<{\ensuremath{\delta}}_{c}^{>}(\ensuremath{\kappa})$. Most interestingly, for values of $\ensuremath{\kappa}$ in the range ${\ensuremath{\kappa}}_{c}(0)<\ensuremath{\kappa}<{\ensuremath{\kappa}}^{>}\ensuremath{\approx}0.215(2)$, we find a reentrant behavior such that N\'eel order exists only in the range ${\ensuremath{\delta}}_{c}^{<}(\ensuremath{\kappa})<\ensuremath{\delta}<{\ensuremath{\delta}}_{c}^{>}(\ensuremath{\kappa})$, with ${\ensuremath{\delta}}_{c}^{<}(\ensuremath{\kappa})>0$. These latter upper and lower critical values coalesce when $\ensuremath{\kappa}={\ensuremath{\kappa}}^{>}$, such that ${\ensuremath{\delta}}_{c}^{<}({\ensuremath{\kappa}}^{>})={\ensuremath{\delta}}_{c}^{>}({\ensuremath{\kappa}}^{>})\ensuremath{\approx}0.25(5)$.
We study the zero-temperature quantum phase diagram of a spin-1 Heisenberg antiferromagnet on the honeycomb lattice with both nearest-neighbor exchange coupling J1>0 and frustrating next-nearest-neighbor coupling ${J}_{2}\ensuremath{\equiv}\ensuremath{\kappa}{J}_{1}>0$, using the coupled cluster method implemented to high orders of approximation, and based on model states with different forms of classical magnetic order. For each we calculate directly in the bulk thermodynamic limit both ground-state low-energy parameters (including the energy per spin, magnetic order parameter, spin stiffness coefficient, and zero-field uniform transverse magnetic susceptibility) and their generalized susceptibilities to various forms of valence-bond crystalline (VBC) order, as well as the energy gap to the lowest-lying spin-triplet excitation. In the range $0<\ensuremath{\kappa}<1$ we find evidence for four distinct phases. Two of these are quasiclassical phases with antiferromagnetic long-range order, one with two-sublattice N\'eel order for $\ensuremath{\kappa}<{\ensuremath{\kappa}}_{{c}_{1}}=0.250(5)$, and another with four-sublattice N\'eel-II order for $\ensuremath{\kappa}>{\ensuremath{\kappa}}_{{c}_{2}}=0.340(5)$. Two different paramagnetic phases are found to exist in the intermediate region. Over the range ${\ensuremath{\kappa}}_{{c}_{1}}<\ensuremath{\kappa}<{\ensuremath{\kappa}}_{c}^{i}=0.305(5)$ we find a gapless phase with no discernible magnetic order, which is a strong candidate for being a quantum spin liquid, while over the range ${\ensuremath{\kappa}}_{c}^{i}<\ensuremath{\kappa}<{\ensuremath{\kappa}}_{{c}_{2}}$ we find a gapped phase, which is most likely a lattice nematic with staggered dimer VBC order that breaks the lattice rotational symmetry.
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.
The ground-state (GS) phase diagram of the frustrated spin- s J 1 -J 2 -J 3 Heisenberg antiferromagnet on the honeycomb lattice is studied using the coupled cluster method implemented to high orders of approximation, for spin quantum numbers s = 1, 3/2, 2 , 5/2. The model has antiferromagnetic (AFM) nearest-neighbour, next-nearest-neighbour and next-next-nearest-neighbour exchange couplings (with strength J 1 > 0, J 2 > 0 and J 3 > 0, respectively). We specifically study the case J 3 = J 2 = κJ 1 , in the range 0 < κ < 1 of the frustration parameter, which includes the point of maximum classical (s → ∞ ) frustration, viz., the classical critical point at κ cl = 1/2 , which separates the Neel phase for κ < κ cl and the collinear striped AFM phase for κ > κ cl . Results are presented for the GS energy, magnetic order parameter and plaquette valence-bond crystal (PVBC) susceptibility. For all spins s > 3/2 we find a quantum phase diagram very similar to the classical one, with a direct first-order transition between the two collinear AFM states at a value κ c (s) which is slightly greater than κ c l [e.g., κ c (3/2) ≈ 0.53(1)] and which approaches it monotonically as s → ∞ . By contrast, for the case s = 1 the transition is split into two such that the stable GS phases are one with Néel AFM order for κ < κ c1 = 0.485(5) and one with striped AFM order for κ > κ c2 = 0.528(5), just as in the case s = 1/2 (for which κ c1 ≈ 0.47 and κ c2 ≈ 0.60). For both the s = 1/2 and s = 1 models the transition at κ c2 appears to be of first-order type, while that at κ c1 appears to be continuous. However, whereas in the s = 1/2 case the intermediate phase appears to have PVBC order over the entire range κ c1 < κ < κ c2 , in the s = 1 case PVBC ordering either exists only over a very small part of the region or, more likely, is absent everywhere.