Motivated by recent density-matrix renormalization group (DMRG) calculations [Yan, Huse, and White, Science 332, 1173 (2011)], which claimed that the ground state of the nearest-neighbor spin-
1/2 Heisenberg antiferromagnet on the kagome lattice geometry is a fully gapped spin liquid with numerical signatures of
Z2 gauge structure, and a further theoretical work [Lu, Ran, and Lee, Phys. Rev. B 83, 224413 (2011)], which gave a classification of all Schwinger-fermion mean-field fully symmetric
Z2 spin liquids on the kagome lattice, we have thoroughly studied Gutzwiller-projected fermionic wave functions by using quantum variational Monte Carlo techniques, hence implementing exactly the constraint of one fermion per site. In particular, we investigated the energetics of all
Z2 candidates (gapped and gapless) that lie in the neighborhood of the energetically competitive U(1) gapless spin liquids. By using a state-of-the-art optimization method, we were able to conclusively show that the U(1) Dirac state is remarkably stable with respect to all
Z2 spin liquids in its neighborhood, and in particular for opening a gap toward the so-called ${\mathbb{Z}}_{2}[0,\ensuremath{\pi}]\ensuremath{\beta}$ state, which was conjectured to describe the ground state obtained by the DMRG method. Finally, we also considered the addition of a small second nearest-neighbor exchange coupling of both antiferromagnetic and ferromagnetic type, and obtained similar results, namely, a U(1) Dirac spin-liquid ground state.
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.
We present a general review of the projective symmetry group classification of fermionic quantum spin liquids for lattice models of spin
S=1/2. We then introduce a systematic generalization of the approach for symmetric
Z2 quantum spin liquids to the one of chiral phases (i.e., singlet states that break time reversal and lattice reflection, but conserve their product). We apply this framework to classify and discuss possible chiral spin liquids on triangular and kagome lattices. We give a detailed prescription on how to construct quadratic spinon Hamiltonians and microscopic wave functions for each representation class on these lattices. Among the chiral
Z2 states, we study the subset of U(1) phases variationally in the antiferromagnetic ${J}_{1}\text{\ensuremath{-}}{J}_{2}\text{\ensuremath{-}}{J}_{d}$ Heisenberg model on the kagome lattice. We discuss static spin structure factors and symmetry constraints on the bulk spectra of these phases.