Quantum spin liquids (QSLs) are long-range entangled phases of frustrated magnets exhibiting fractionalised spin excitations. In two dimensions, there is limited analytical understanding of their excitation spectra beyond parton mean-field theories, which fail to capture many features of the finite frequency dynamical response. We use a self-consistent mean-field theory combined with the random phase approximation (RPA) for the
J1-
J2 Heisenberg model on the triangular lattice to describe the strong spinon-spinon interactions of the U(1) Dirac QSL. We obtain quantitative results for the dynamical spin structure factor and phase diagram compatible with comprehensive numerical efforts. We show the continuum response at the Brillouin zone corners is described by a spinon-exciton hybridised with monopole gauge excitations. We further show that the transition from the QSL to 120-degree coplanar order is in the QED
3 O(3)-GN universality class. We extend the method to chiral QSLs and XXZ anisotropy, and discuss its broad range of applicability to other models and for describing inelastic neutron scattering experiments in KYbSe
2.
Understanding the nature of the excitation spectrum in quantum spin liquids is of fundamental importance, in particular for the experimental detection of candidate materials.However, current theoretical and numerical techniques have limited capabilities, especially in obtaining the dynamical structure factor, which gives a crucial characterization of the ultimate nature of the quantum state and may be directly assessed by inelastic neutron scattering.In this work, we investigate the low-energy properties of the S ¼ 1=2 Heisenberg model on the triangular lattice, including both nearest-neighbor J 1 and next-nearestneighbor J 2 superexchanges, by a dynamical variational Monte Carlo approach that allows accurate results on spin models.For J 2 ¼ 0, our calculations are compatible with the existence of a well-defined magnon in the whole Brillouin zone, with gapless excitations at K points (i.e., at the corners of the Brillouin zone).The strong renormalization of the magnon branch (also including rotonlike minima around the M points, i.e., midpoints of the border zone) is described by our Gutzwiller-projected state, where Abrikosov fermions are subject to a nontrivial magnetic π flux threading half of the triangular plaquettes.When increasing the frustrating ratio J 2 =J 1 , we detect a progressive softening of the magnon branch at M, which eventually becomes gapless within the spin-liquid phase.This feature is captured by the band structure of the unprojected wave function (with two Dirac points for each spin component).In addition, we observe an intense signal at low energies around the K points, which cannot be understood within the unprojected picture and emerges only when the Gutzwiller projection is considered, suggesting the relevance of gauge fields for the low-energy physics of spin liquids.
We study the spin-
21 Heisenberg model on the triangular lattice with antiferromagnetic first-
(J1) and second-
(J2) nearest-neighbor interactions using density matrix renormalization group. By studying the spin correlation function, we find a ${120}^{\ensuremath{\circ}}$ magnetic order phase for ${J}_{2}\ensuremath{\lesssim}0.07{J}_{1}$ and a stripe antiferromagnetic phase for ${J}_{2}\ensuremath{\gtrsim}0.15{J}_{1}$. Between these two phases, we identify a spin-liquid region characterized by exponential decaying spin and dimer correlations, as well as large spin singlet and triplet excitation gaps on finite-size systems. We find two near degenerating ground states with distinct properties in two sectors, which indicates more than one spin-liquid candidate in this region. While the sector with spinons is found to respect time reversal symmetry, the even sector without spinons breaks such a symmetry for finite-size systems. Furthermore, we detect the signature of the fractionalization by following the evolution of different ground states with inserting spin flux into the cylinder system. Moreover, by tuning the anisotropic bond coupling, we explore the nature of the spin-liquid phase and find the optimal parameter region for gapped
Z2 spin liquids.
The
U(1) Dirac spin liquid might realize an exotic phase of matter whose low-energy properties are described by quantum electrodynamics in
2+1 dimensions, where gapless modes exist but spinons and gauge fields are strongly coupled. Its existence was proposed in frustrated Heisenberg models in the presence of frustrating superexchange interactions by the (Abrikosov) fermionic representation of the spin operators [Wen, Phys. Rev. B 65, 165113 (2002)], supplemented by the Gutzwiller projection. Here, we construct charge-
Q monopole excitations in the Heisenberg model on the triangular lattice with nearest-neighbor (
J1) and next-nearest-neighbor (
J2) couplings. In the highly frustrated regime, singlet and triplet monopoles with
Q=1 become gapless in the thermodynamic limit; in addition, the energies for generic
Q agree with field-theoretical predictions, obtained for a large number of gapless fermion modes. Finally, we consider localized gauge excitations, in which magnetic $\ensuremath{\pi}$ fluxes are concentrated in the triangular plaquettes (in analogy with
Z2 visons), showing that these kinds of states do not play a relevant role at low energies. All our findings lend support to a stable
U(1) Dirac spin liquid in the ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on the triangular lattice.
We use a projective symmetry group analysis to determine all symmetric spin liquids on the stuffed honeycomb lattice Heisenberg model. This lattice interpolates between honeycomb, triangular, and dice lattices, always preserving hexagonal symmetry, and it already has one spin-liquid candidate,
TbInO3, albeit with strong spin-orbit coupling not considered here. In addition to the stuffed honeycomb lattice itself, we gain valuable insight into potential spin liquids on the honeycomb and triangular lattices, as well as how they might be connected. For example, the sublattice pairing state proposed on the honeycomb lattice connects to the uniform spinon Fermi surface that may be relevant for the triangular lattice with ring exchange, while there are no spin liquids competitive on both the ${J}_{1}\ensuremath{-}{J}_{2}$ honeycomb and triangular lattice limits. In particular, we find three stuffed honeycomb descendants of the U(1) Dirac spin liquid widely believed to be found on the ${J}_{1}\ensuremath{-}{J}_{2}$ triangular lattice. We also discuss how spin liquids near the honeycomb limit can potentially explain the physics of
LiZn2Mo3O8.
By using Gutzwiller projected fermionic wave functions and variational Monte Carlo technique, we study the spin-
1/2 Heisenberg model with the first-neighbor (
J1), second-neighbor (
J2), and additional scalar chiral interaction ${J}_{\ensuremath{\chi}}{\mathbf{S}}_{i}\ifmmode\cdot\else\textperiodcentered\fi{}({\mathbf{S}}_{j}\ifmmode\times\else\texttimes\fi{}{\mathbf{S}}_{k})$ on the triangular lattice. In the nonmagnetic phase of the ${J}_{1}\ensuremath{-}{J}_{2}$ triangular model with $0.08\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\lesssim}0.16$, recent density-matrix renormalization group (DMRG) studies [Zhu and White, Phys. Rev. B 92, 041105(R) (2015) and Hu, Gong, Zhu, and Sheng, Phys. Rev. B 92, 140403(R) (2015)] find a possible gapped spin liquid with the signal of a competition between a chiral and a
Z2 spin liquid. Motivated by the DMRG results, we consider the chiral interaction ${J}_{\ensuremath{\chi}}{\mathbf{S}}_{i}\ifmmode\cdot\else\textperiodcentered\fi{}({\mathbf{S}}_{j}\ifmmode\times\else\texttimes\fi{}{\mathbf{S}}_{k})$ as a perturbation for this nonmagnetic phase. We find that with growing ${J}_{\ensuremath{\chi}}$, the gapless U(1) Dirac spin liquid, which has the best variational energy for ${J}_{\ensuremath{\chi}}=0$, exhibits the energy instability towards a gapped spin liquid with nontrivial magnetic fluxes and nonzero chiral order. We calculate topological Chern number and ground-state degeneracy, both of which identify this flux state as the chiral spin liquid with fractionalized Chern number
C=1/2 and twofold topological degeneracy. Our results indicate a positive direction to stabilize a chiral spin liquid near the nonmagnetic phase of the ${J}_{1}\ensuremath{-}{J}_{2}$ triangular model.