We use the density matrix renormalization group to perform accurate calculations of the ground state of the nearest-neighbor quantum spin S = 1/2 Heisenberg antiferromagnet on the kagome lattice. We study this model on numerous long cylinders with circumferences up to 12 lattice spacings. Through a combination of very-low-energy and small finite-size effects, our results provide strong evidence that, for the infinite two-dimensional system, the ground state of this model is a fully gapped spin liquid.
We perform a density-matrix renormalization group (DMRG) study of the S=1/2 Heisenberg antiferromagnet on the kagome lattice to identify the conjectured spin liquid ground state. Exploiting SU(2) spin symmetry, which allows us to keep up to 16,000 DMRG states, we consider cylinders with circumferences up to 17 lattice spacings and find a spin liquid ground state with an estimated per site energy of -0.4386(5), a spin gap of 0.13(1), very short-range decay in spin, dimer and chiral correlation functions and finite topological entanglement gamma consistent with gamma=ld 2, ruling out gapless, chiral or non-topological spin liquids. At the same time, DMRG results provide strong evidence for a gapped topological Z2 spin liquid.
The spin-
21 Heisenberg model on the kagome lattice, which is closely realized in layered Mott insulators such as
ZnCu3(OH)6Cl2, is one of the oldest and most enigmatic spin-
21 lattice models. While the numerical evidence has accumulated in favor of a quantum spin liquid, the debate is still open as to whether it is a
Z2 spin liquid with very short-range correlations (some kind of resonating valence bond spin liquid), or an algebraic spin liquid with power-law correlations. To address this issue, we have pushed the program started by Rokhsar and Kivelson in their derivation of the effective quantum dimer model description of Heisenberg models to unprecedented accuracy for the spin-
21 kagome, by including all the most important virtual singlet contributions on top of the orthogonalization of the nearest-neighbor valence bond singlet basis. Quite remarkably, the resulting picture is a competition between a
Z2 spin liquid and a diamond valence bond crystal with a 12-site unit cell, as in the density-matrix renormalization group simulations of Yan et al. Furthermore, we found that, on cylinders of finite diameter
d, there is a transition between the
Z2 spin liquid at small
d and the diamond valence bond crystal at large
d, the prediction of the present microscopic description for the two-dimensional lattice. These results show that, if the ground state of the spin-
21 kagome antiferromagnet can be described by nearest-neighbor singlet dimers, it is a diamond valence bond crystal, and, a contrario, that, if the system is a quantum spin liquid, it has to involve long-range singlets, consistent with the algebraic spin liquid scenario.
We revisit the description of the low-energy singlet sector of the spin-
21 Heisenberg antiferromagnet on kagome in terms of an effective quantum dimer model. With the help of exact diagonalizations of appropriate finite-size clusters, we show that the embedding of a given process in its kagome environment leads to dramatic modifications of the amplitudes of the elementary loop processes, an effect not accessible to the standard approach based on the truncation of the Hamiltonian to the nearest-neighbor valence-bond basis. The resulting parameters are consistent with a
Z2 spin liquid rather than with a valence-bond crystal, in agreement with the last density matrix renormalization group results.
Gapped
Z2 spin liquids have been proposed as candidates for the ground state of the
S=1/2 quantum antiferromagnet on the kagome lattice. We extend the use of projected entangled pair states to construct (on the cylinder) resonating valence bond (RVB) states including both nearest-neighbor and next-nearest-neighbor singlet bonds. Our ansatz---dubbed ``simplex spin liquid''---allows for an asymmetry between the two types of triangles (of order
2 in the energy density after optimization) leading to the breaking of inversion symmetry. We show that the topological
Z2 structure is still preserved and, by considering the presence or the absence of spinon and vison lines along an infinite cylinder, we explicitly construct four orthogonal RVB minimally entangled states. The spinon and vison coherence lengths are extracted from a finite size scaling with regard to the cylinder perimeter of the energy splittings of the four sectors and are found to be of the order of the lattice spacing. The entanglement spectrum of a partitioned (infinite) cylinder is found to be gapless, suggesting the occurrence, on a cylinder with real open boundaries, of gapless edge modes formally similar to Luttinger liquid (nonchiral) spin and charge modes. When inversion symmetry is spontaneously broken, the RVB spin liquid exhibits an extra Ising degeneracy, which might have been observed in recent exact diagonalization studies.
Due to strong geometric frustration and quantum fluctuation, the
S=1/2 quantum Heisenberg antiferromagnet on the
kagome lattice has long been considered as an ideal platform to realize a spin liquid (SL), a phase exhibiting fractionalized excitations without any symmetry breaking. A recent numerical study (Yan et al., e-print arXiv:1011.6114) of the Heisenberg
S=1/2,kagome lattice model (HKLM) shows, in contrast to earlier results, that the ground state is a singlet-gapped SL with signatures of
Z2 topological order. Motivated by this numerical discovery, we use the projective symmetry group to classify all 20 possible Schwinger fermion mean-field states of
Z2 SLs on the
kagome lattice. Among them we found only one gapped
Z2 SL (which we call the ${\mathbb{Z}}_{2}[0,\ensuremath{\pi}]\ensuremath{\beta}$ state) in the neighborhood of the U(1) Dirac SL state. Since its parent state, i.e., the U(1) Dirac SL, was found [Ran et al., Phys. Rev. Lett. 98, 117205 (2007)] to be the lowest among many other candidate U(1) SLs, including the uniform resonating-valence-bond states, we propose this ${\mathbb{Z}}_{2}[0,\ensuremath{\pi}]\ensuremath{\beta}$ state to be the numerically discovered SL ground state of the HKLM.