Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1708.06162 · Physical review. B./Physical review. B

High-order study of the quantum critical behavior of a frustrated spin- 1 2 antiferromagnet on a stacked honeycomb bilayer

R. F. Bishop, P. H. Y. Li

Published Dec 12, 2017• 12 citations
arXivHTMLPDF

Abstract

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 ${J}_{1},\phantom{\rule{0.16em}{0ex}}{J}_{2}$, and ${J}_{3}$, 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)$.

DOI
10.1103/physrevb.96.224416

Related papers

Frustrated honeycomb-bilayer Heisenberg antiferromagnet: The spin- 12 J1−J2−J1⊥ model

R. F. Bishop, P. H. Y. Li

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 ${J}_{1}>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)$.

arXiv:1611.03287•Physical review. B./Physical review. B•Apr 11, 2017•22 citations

Quantum phase diagram of a frustrated spin-$\frac{1}{2}$ $J_{1}$--$J_{2}$--$J_{1}^{\perp}$ Heisenberg model on an $AA$-stacked bilayer honeycomb lattice

P. H. Y. Li, R. F. Bishop

The zero-temperature quantum phase diagram of the spin-$\frac{1}{2}$ $J_{1}$--$J_{2}$--$J_{1}^{\perp}$ 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 $J_{1} > 0$ and $J_{2} \equiv \kappa J_{1}>0$. The two layers are coupled via a comparable Heisenberg exchange interaction between NN interlayer pairs, with a strength $J_{1}^{\perp} \equiv \delta J_{1}$. The complete phase boundaries of two quasiclassical collinear AFM phases, namely the Neel and Neel-II phases, are calculated in the $\kappa \delta$ half-plane with $\kappa > 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 $\Delta$ from the $s^{z}_{T}=0$ ground state to the lowest-lying $|s^{z}_{T}|=1$ excited state (where $s^{z}_{T}$ 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 \leq 10$) in a systematic series of approximations involving $n$-spin clusters. The sole approximation made is then to extrapolate the sequences of $n$th-order results for $M$ and $\Delta$ to the exact limit, $n \to \infty$.

•May 3, 2018•0 citations

Phase diagram of a frustrated Heisenberg antiferromagnet on the honeycomb lattice: The J1 - J2 - J3 model

P. H. Y. Li, R. F. Bishop, DAMIAN J. J. FARNELL +1 more

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 ${J}_{1}$-${J}_{2}$-${J}_{3}$ 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 ${J}_{3}<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.

arXiv:1204.2390•Physical Review B•Oct 9, 2012•99 citations

Ground-state phases of the frustrated spin- 12 J1 – J2 – J3 Heisenberg ferromagnet ( J1<0 ) on the honeycomb lattice with J3=J2>0

P. H. Y. Li, R. F. Bishop, DAMIAN J. J. FARNELL +2 more

We study the ground-state (GS) properties of the frustrated spin-$\frac{1}{2}$ ${J}_{1}$--${J}_{2}$--${J}_{3}$ Heisenberg model on the two-dimensional honeycomb lattice with ferromagnetic nearest-neighbor (${J}_{1}=\ensuremath{-}1$) exchange and frustrating antiferromagnetic next-nearest-neighbor (${J}_{2}>0$) and next-next-nearest-neighbor (${J}_{3}>0$) exchanges, for the case ${J}_{3}={J}_{2}$. 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 ${J}_{2}^{\mathrm{cl}}=\frac{1}{7}$ ($\ensuremath{\approx}0.143$) for the classical version of the model, in which we momentarily ignore the intervening noncollinear spiral phase in the region $\frac{1}{10}<{J}_{2}<\frac{1}{5}$. 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 ${J}_{1}=+1$). The magnetic order parameter is found to behave similarly for the ferromagnetic and the antiferromagnetic models for large values of the frustration parameter ${J}_{2}$. 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 ${J}_{1}=+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 ${J}_{1}=+1$), where a plaquette valence-bond crystal phase intrudes between the two corresponding quasiclassical antiferromagnetic phases (with N\'eel and striped order) for ${J}_{2}^{{c}_{1}}<{J}_{2}<{J}_{2}^{{c}_{2}}$, 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 $\frac{1}{10}<{J}_{2}<\frac{1}{5}$.

arXiv:1109.6229•Physical Review B•Feb 21, 2012•31 citations

Transverse magnetic susceptibility of a frustrated spin-12J1−J2−J1⊥ Heisenberg antiferromagnet on a bilayer honeycomb lattice

P. H. Y. Li, R. F. Bishop

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.

arXiv:1709.08205•AIP conference proceedings•Jan 1, 2017•2 citations

Ground-state phases of the spin- 1 J 1 − J 2 Heisenberg antiferromagnet on the honeycomb lattice

P. H. Y. Li, R. F. Bishop

We study the zero-temperature quantum phase diagram of a spin-1 Heisenberg antiferromagnet on the honeycomb lattice with both nearest-neighbor exchange coupling ${J}_{1}>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.

arXiv:1602.08915•Physical review. B./Physical review. B•Jun 28, 2016•21 citations
We study a frustrated spin-12J1\ensuremath-J2\ensuremath-J3\ensuremath-J1\ensuremath⊥\frac{1}{2}\phantom{\rule{0.16em}{0ex}}{J}_{1}\text{\ensuremath{-}}{J}_{2}\text{\ensuremath{-}}{J}_{3}\text{\ensuremath{-}}{J}_{1}^{\ensuremath{\perp}}21​J1​\ensuremath-J2​\ensuremath-J3​\ensuremath-J1\ensuremath⊥​ Heisenberg antiferromagnet on an AAAAAA-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{J}_{1},\phantom{\rule{0.16em}{0ex}}{J}_{2}J1​,J2​, and J3{J}_{3}J3​, respectively. The two layers are coupled with an AFM NN exchange coupling J1\ensuremath⊥\ensuremath≡\ensuremathδJ1{J}_{1}^{\ensuremath{\perp}}\ensuremath{\equiv}\ensuremath{\delta}{J}_{1}J1\ensuremath⊥​\ensuremath≡\ensuremathδJ1​. The model is studied for arbitrary values of \ensuremathδ\ensuremath{\delta}\ensuremathδ along the line J3=J2\ensuremath≡\ensuremathαJ1{J}_{3}={J}_{2}\ensuremath{\equiv}\ensuremath{\alpha}{J}_{1}J3​=J2​\ensuremath≡\ensuremathαJ1​ that includes the most highly frustrated point at \ensuremathα=12\ensuremath{\alpha}=\frac{1}{2}\ensuremathα=21​, 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α\ensuremathδ\ensuremath{\alpha}\ensuremath{\delta}\ensuremathα\ensuremathδ plane in the window 0\ensuremath≤\ensuremathα\ensuremath≤1,0\ensuremath≤\ensuremathδ\ensuremath≤10\ensuremath{\le}\ensuremath{\alpha}\ensuremath{\le}1,\phantom{\rule{0.16em}{0ex}}0\ensuremath{\le}\ensuremath{\delta}\ensuremath{\le}10\ensuremath≤\ensuremathα\ensuremath≤1,0\ensuremath≤\ensuremathδ\ensuremath≤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δ=\ensuremathδci(\ensuremathα)\ensuremath{\delta}={\ensuremath{\delta}}_{{c}_{i}}(\ensuremath{\alpha})\ensuremathδ=\ensuremathδci​​(\ensuremathα), or equivalently, \ensuremathα=\ensuremathαci(\ensuremathδ)\ensuremath{\alpha}={\ensuremath{\alpha}}_{{c}_{i}}(\ensuremath{\delta})\ensuremathα=\ensuremathαci​​(\ensuremathδ), with i=1i=1i=1 (N'eel) and 2 (striped). The two boundaries exhibit an ``avoided crossing'' behavior with both curves being re-entrant. Thus, in this \ensuremathα\ensuremathδ\ensuremath{\alpha}\ensuremath{\delta}\ensuremathα\ensuremathδ window, N'eel order exists only for values of \ensuremathδ\ensuremath{\delta}\ensuremathδ in the range \ensuremathδc1<(\ensuremathα)<\ensuremathδ<\ensuremathδc1>(\ensuremathα){\ensuremath{\delta}}_{{c}_{1}}^{<}(\ensuremath{\alpha})<\ensuremath{\delta}<{\ensuremath{\delta}}_{{c}_{1}}^{>}(\ensuremath{\alpha})\ensuremathδc1​<​(\ensuremathα)<\ensuremathδ<\ensuremathδc1​>​(\ensuremathα), with \ensuremathδc1<(\ensuremathα)=0{\ensuremath{\delta}}_{{c}_{1}}^{<}(\ensuremath{\alpha})=0\ensuremathδc1​<​(\ensuremathα)=0 for \ensuremathα<\ensuremathαc1(0)\ensuremath≈0.46(2)\ensuremath{\alpha}<{\ensuremath{\alpha}}_{{c}_{1}}(0)\ensuremath{\approx}0.46(2)\ensuremathα<\ensuremathαc1​​(0)\ensuremath≈0.46(2) and \ensuremathδc1<(\ensuremathα)>0{\ensuremath{\delta}}_{{c}_{1}}^{<}(\ensuremath{\alpha})>0\ensuremathδc1​<​(\ensuremathα)>0 for \ensuremathαc1(0)<\ensuremathα<\ensuremathα1>\ensuremath≈0.49(1){\ensuremath{\alpha}}_{{c}_{1}}(0)<\ensuremath{\alpha}<{\ensuremath{\alpha}}_{1}^{>}\ensuremath{\approx}0.49(1)\ensuremathαc1​​(0)<\ensuremathα<\ensuremathα1>​\ensuremath≈0.49(1), and striped order similarly exists only for values of \ensuremathδ\ensuremath{\delta}\ensuremathδ in the range \ensuremathδc2<(\ensuremathα)<\ensuremathδ<\ensuremathδc2>(\ensuremathα){\ensuremath{\delta}}_{{c}_{2}}^{<}(\ensuremath{\alpha})<\ensuremath{\delta}<{\ensuremath{\delta}}_{{c}_{2}}^{>}(\ensuremath{\alpha})\ensuremathδc2​<​(\ensuremathα)<\ensuremathδ<\ensuremathδc2​>​(\ensuremathα), with \ensuremathδc2<(\ensuremathα)=0{\ensuremath{\delta}}_{{c}_{2}}^{<}(\ensuremath{\alpha})=0\ensuremathδc2​<​(\ensuremathα)=0 for \ensuremathα>\ensuremathαc2(0)\ensuremath≈0.600(5)\ensuremath{\alpha}>{\ensuremath{\alpha}}_{{c}_{2}}(0)\ensuremath{\approx}0.600(5)\ensuremathα>\ensuremathαc2​​(0)\ensuremath≈0.600(5) and \ensuremathδc2<(\ensuremathα)>0{\ensuremath{\delta}}_{{c}_{2}}^{<}(\ensuremath{\alpha})>0\ensuremathδc2​<​(\ensuremathα)>0 for \ensuremathαc2(0)>\ensuremathα>\ensuremathα2<\ensuremath≈0.56(1){\ensuremath{\alpha}}_{{c}_{2}}(0)>\ensuremath{\alpha}>{\ensuremath{\alpha}}_{2}^{<}\ensuremath{\approx}0.56(1)\ensuremathαc2​​(0)>\ensuremathα>\ensuremathα2<​\ensuremath≈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{J}_{1}>0J1​>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)$.
Quantum phase diagram of a frustrated spin-12\frac{1}{2}21​ J1J_{1}J1​--J2J_{2}J2​--J1⊥J_{1}^{\perp}J1⊥​ Heisenberg model on an AAAAAA-stacked bilayer honeycomb lattice
The zero-temperature quantum phase diagram of the spin-12\frac{1}{2}21​ J1J_{1}J1​--J2J_{2}J2​--J1⊥J_{1}^{\perp}J1⊥​ model on an AAAAAA-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>0J_{1} > 0J1​>0 and J2≡κJ1>0J_{2} \equiv \kappa J_{1}>0J2​≡κJ1​>0. The two layers are coupled via a comparable Heisenberg exchange interaction between NN interlayer pairs, with a strength J1⊥≡δJ1J_{1}^{\perp} \equiv \delta J_{1}J1⊥​≡δJ1​. The complete phase boundaries of two quasiclassical collinear AFM phases, namely the Neel and Neel-II phases, are calculated in the κδ\kappa \deltaκδ half-plane with κ>0\kappa > 0κ>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 MMM and the excitation energy Δ\DeltaΔ from the sTz=0s^{z}_{T}=0sTz​=0 ground state to the lowest-lying ∣sTz∣=1|s^{z}_{T}|=1∣sTz​∣=1 excited state (where sTzs^{z}_{T}sTz​ is the total zzz component of spin for the system as a whole, and where the collinear ordering lies along the zzz 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≤10n \leq 10n≤10) in a systematic series of approximations involving nnn-spin clusters. The sole approximation made is then to extrapolate the sequences of nnnth-order results for MMM and Δ\DeltaΔ to the exact limit, n→∞n \to \inftyn→∞.
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{J}_{1}J1​-J2{J}_{2}J2​-J3{J}_{3}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{J}_{3}<0J3​<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-12\frac{1}{2}21​ J1{J}_{1}J1​--J2{J}_{2}J2​--J3{J}_{3}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{J}_{2}>0J2​>0) and next-next-nearest-neighbor (J3>0{J}_{3}>0J3​>0) exchanges, for the case J3=J2{J}_{3}={J}_{2}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=17{J}_{2}^{\mathrm{cl}}=\frac{1}{7}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 110<J2<15\frac{1}{10}<{J}_{2}<\frac{1}{5}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{J}_{1}=+1J1​=+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{J}_{2}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{J}_{1}=+1J1​=+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{J}_{1}=+1J1​=+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{J}_{2}^{{c}_{1}}<{J}_{2}<{J}_{2}^{{c}_{2}}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 110<J2<15\frac{1}{10}<{J}_{2}<\frac{1}{5}101​<J2​<51​.
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 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{J}_{1}>0J1​>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.