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.
A complicating factor in the realization and observation of quantum spin liquids in materials is the ubiquitous presence of other degrees of freedom, in particular lattice distortion modes (phonons). These provide additional routes for relieving magnetic frustration, thereby possibly destabilizing spin-liquid ground states. In this work, we focus on triangular-lattice Heisenberg antiferromagnets, where recent numerical evidence suggests the presence of an extended U(1) Dirac spin liquid phase which is described by compact quantum electrodynamics in 2+1 dimensions (QED
3), featuring gapless spinons and monopoles as gauge excitations, and believed to flow to a strongly-coupled fixed point with conformal symmetry. Using complementary perturbation theory and scaling arguments, we show that a symmetry-allowed coupling between (classical) finite-wavevector lattice distortions and monopole operators of the U(1) Dirac spin liquid generally induces a spin-Peierls instability towards a (confining) 12-site valence-bond solid state. We support our theoretical analysis with state-of-the-art density matrix renormalization group simulations. Away from the limit of static distortions, we demonstrate that the phonon energy gap establishes a parameter regime where the spin liquid is expected to be stable, and show that the monopole-lattice coupling leads to softening of the phonon in analogy to the Kohn anomaly. We discuss the applicability of our results to similar systems, in particular the Dirac spin liquid on the Kagome lattice.
We study the spin-1/2 Heisenberg model on the triangular lattice with the antiferromagnetic first (J1) and second nearest-neighbor (J2) interactions using density matrix renormalization group. By studying the spin correlation function, we find a 120◦ magnetic order phase for J2 . 0.08J1 and a stripe antiferromagnetic phase for J2 & 0.16J2. Between these two phases, we identify a gapped spin liquid phase characterized by the exponential decaying spin and dimer correlations, as well as finite spin singlet and triplet excitation gaps. We also find two topological degenerating ground states in two topological sectors. Furthermore, we detect the signature of the fractionalization by following the evolution of different ground states with inserting spin flux into the cylinder system. While the sector with spinon is found to respect the time reveral symmetry, the even sector without a spinon may break such a symmetry. We discuss the implications of numerical results for the nature of the topological spin liquid phase.
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.
We study the energy and the static spin structure factor of the ground state of the spin-
1/2 quantum Heisenberg antiferromagnetic model on the kagome lattice. By the iterative application of a few Lanczos steps on accurate projected fermionic wave functions and the Green's function Monte Carlo technique, we find that a gapless (algebraic)
U(1) Dirac spin liquid is competitive with previously proposed gapped (topological)
Z2 spin liquids. By performing a finite-size extrapolation of the ground-state energy, we obtain an energy per site $E/J=\ensuremath{-}0.4365(2)$, which is equal, within three error bars, to the estimates given by the density-matrix renormalization group (DMRG). Our estimate is obtained for a translationally invariant system, and, therefore, does not suffer from boundary effects, like in DMRG. Moreover, on finite toric clusters at the pure variational level, our energies are lower compared to those from DMRG calculations.
The dynamical dimer structure factor is an observable probing spin-singlet excitations of quantum magnets distinct from those commonly studied by the spin structure factor. We report the dimer response for the extended spin-
1/2 antiferromagnetic Heisenberg model on the triangular lattice using large-scale GPU-accelerated matrix-product-state simulations. We investigate the ordered phases with
120∘ coplanar, collinear stripe, and tetrahedral spin order, as well as candidate quantum spin-liquid (QSL) regimes, comprising an expected gapless
U(1) Dirac QSL and a chiral QSL at finite spin-scalar-chirality coupling. In the ordered phases, we find low-energy modes below the onset of the two-magnon continuum illustrating avoided quasiparticle decay. Within the candidate gapless QSL, we observe absolute dispersion minima at momenta of half the Brillouin zone corners,
X≡K/2, in agreement with field-theory predictions that singlet monopole excitations of the
U(1) Dirac spin liquid become gapless at these points. Thus, the high-resolution dynamical dimer response provides support for a
U(1) Dirac QSL with singlet monopole excitations.