We present an extensive quantum Monte Carlo study of a nearest-neighbor, singlet-projection model on the kagome lattice that exhibits SO(
N) symmetry and is sign-problem-free. We find that in contrast to the previously studied SU(
N) variations of this model, the non-bipartite lattices appear to harbor spin-liquid phases for intermediate values of
N, a result also seen on the triangular lattice with this same model. Unlike the triangular lattice, however, the kagome lattice appears to admit no valence bond solid (VBS) phase for large
N, instead remaining a spin-liquid. We also observe that the spin-ordered phase survives to a relatively large value of
N, at least
N=8, and that it is gone for
N=10; the fate of
N=9 remains unclear.
The S=1/2 kagome-lattice antiferromagnet is one of interesting frustrated quantum spin systems. The system is supposed to exhibit the quantum spin liquid in the ground state, which was proposed as a possible origin of the high-temprature superconductivity. The spin gap is an important physical quantity to characterize the spin liquid behavior. Whether the S=1/2 kagome-lattice antiferromagnet is gapless or has a finite spin gap, is still an unsolved issue. Because any recently developped numerical calculation methods are not enough to determine it in the thermodynamic limit. Our large-scale numerical diazonalization up to 42-spin clusters and a finite-size scaling analysis indicated that the S=1/2 kagome-lattice antiferromagnet is gapless in the thremodynamic limit[1]. It is consistent with the U(1) Dirac spin liquid theory of the kagome-lattice antiferromagnet[2,3]. On the other hand, some density matrix renormalization group (DMRG) calculations supported the gapped Z2 topological spin liquid theory[4,5]. Our recent numerical diagonalization analysis on the magnetization process of a distorted kagome-lattice antiferromagnet indicated that the perfect kagome-lattice system is just on a quantum critical point[6]. It would be a possible reason why it is difficult to determine whether the perfect kagome-lattice antiferromagnet is gapless or gapped. In the present study, a field derivative analysis is applied for the spin gap issue of the kagome-lattice antiferromanget, using the numerical diagonalization data up to 42-spin clusters. It is consistent with our previous conclusion that it is gapless[7]. [1]H. Nakano and T. Sakai: J. Phys. Soc. Jpn. 80 (2011) 053704. [2]Y. Ran, M. Hermele, P. A. Lee and X. -G. Wen: Phys. Rev. Lett. 98 (2007) 117205. [3]Y. Iqbal, F. Becca, S. Sorella and D. Poilblanc: Phys. Rev. B 87 (2013) 060405(R). [4]S. Yan, D. A. Huse and S. R. White: Science 332 (2011) 1173. [5]S. Nishimoto, N. Shibata and C. Hotta: Nat. Commun. 4 (2013) 2287. [6]H. Nakano and T. Sakai: J. Phys. Soc. Jpn. 83 (2014) 104710. [7]H. Nakano and T. Sakai: in preparation.
We investigate the spin-
21 antiferromagnetic Heisenberg model with a Dzyaloshinskii-Moriya interaction on kagome lattice, making use of the variational Monte Carlo technique. An exotic quantum spin state is found to arise from a melting of the $\mathbit{Q}=0$ long-range magnetic order by a topological transition, when a small anisotropic third-nearest-neighbor antiferromagnetic Heisenberg interaction is turned on. This state is a gapped quantum spin liquid, characterized by a topological order with ground-state degeneracy
ng=4 and topological entanglement entropy $\ensuremath{\gamma}=ln2$, suggesting it is an Abelian topological phase. Furthermore, the Chern numbers of the spin-up (-down) spinon occupied bands of this state are ${C}_{\ensuremath{\uparrow}\ensuremath{\downarrow}}=\ifmmode\pm\else\textpm\fi{}1$, respectively. From this perspective, this state is also a time-reversal symmetric (total Chern number
Ctotal=0) topological insulator with spinons as the chiral edge states, which carry opposite spin and move in the opposite direction. It is analogous to the quantum spin Hall state but the spin current is carried by deconfined spinons in a quantum spin liquid, and so is dubbed as the spinon quantum spin Hall state.
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.
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.
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.