The properties of ground state of spin- \frac{1}{2} 1 2 kagome antiferromagnetic Heisenberg (KAFH) model have attracted considerable interest in the past few decades, and recent numerical simulations reported a spin liquid phase. The nature of the spin liquid phase remains unclear. For instance, the interplay between symmetries and Z_2 Z 2 topological order leads to different types of Z_2 Z 2 spin liquid phases. In this paper, we develop a numerical simulation method based on symmetric projected entangled-pair states (PEPS), which is generally applicable to strongly correlated model systems in two spatial dimensions. We then apply this method to study the nature of the ground state of the KAFH model. Our results are consistent with that the ground state is a U(1) U ( 1 ) Dirac spin liquid rather than a Z_2 Z 2 spin liquid.
Competing Spin Liquid Phases in the S=21 Heisenberg Model on the Kagome Lattice
The properties of ground state of spin- \frac{1}{2} 1 2 kagome antiferromagnetic Heisenberg (KAFH) model have attracted considerable interest in the past few decades, and recent numerical simulations reported a spin liquid phase. The nature of the spin liquid phase remains unclear. For instance, the interplay between symmetries and Z_2 Z 2 topological order leads to different types of Z_2 Z 2 spin liquid phases. In this paper, we develop a numerical simulation method based on symmetric projected entangled-pair states (PEPS), which is generally applicable to strongly correlated model systems in two spatial dimensions. We then apply this method to study the nature of the ground state of the KAFH model. Our results are consistent with that the ground state is a U(1) U ( 1 ) Dirac spin liquid rather than a Z_2 Z 2 spin liquid.
Z2 spin liquid ground state of the spin-21 Kagome antiferromagnetic Heisenberg model
We show that the best variational ground state of the spin-21 Kagome antiferromagnetic Heisenberg model with nearest-neighboring exchange coupling(NN-KAFH) is a Z2 spin liquid state, rather than the widely believed U(1) Dirac spin liquid state. The spinon excitation in the Z2 spin liquid state has a small gap of about 1/40 of the spinon band width. We find that while the Z2 and the U(1) spin liquid state have a large overlap on finite clusters and are thus very close in energy, they host totally different spinon excitation spectrum. The strength of the RVB theory becomes particularly clear in such a situation, as it provides not only a variational understanding of the ground state structure, but also a comprehensive picture for the excitation spectrum of the system. Our result indicates that the spin-21 NN-KAFH should be better understood as a nearly critical system, rather than a prototypical gapped Z2 spin liquid system.
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.
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.
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 provide strong numerical evidence, using improved variational wave functions, for a ground state with vanishing spin gap in the spin-1/2 quantum Heisenberg model on the kagome lattice. Starting from the algebraic U(1) Dirac spin liquid state proposed by Ran et al. [Phys. Rev. Lett. 98, 117205 (2007)] and iteratively applying a few Lanczos steps, we compute the lowest S=2 excitation constructed by exciting spinons close to the Dirac nodes. Our results are compatible with a vanishing spin gap in the thermodynamic limit and in consonance with a power-law decay of long distance spin-spin correlations in real space. The competition with a gapped (topological) spin liquid is discussed.
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 [Z subscript 2] 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 [Z subscript 2] SLs on the kagome lattice. Among them we found only one gapped Z2 [Z subscript 2] SL (which we call the Z2[0,π]β [Z subscript 2 [0,pi] 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 Z2[0,π]β [Z subscript 2 [0,pi] Beta] state to be the numerically discovered SL ground state of the HKLM.