Valence-bond crystal in the extended kagome spin-1/2 quantum Heisenberg antiferromagnet: A variational Monte Carlo approach

The highly-frustrated spin-1/2 quantum Heisenberg model with both nearest ($J_1$) and next-nearest ($J_2$) neighbor exchange interactions is revisited by using an extended variational space of projected wave functions that are optimized with state-of-the-art methods. Competition between modulated valence-bond crystals (VBCs) proposed in the literature and the Dirac spin liquid (DSL) is investigated. We find that the addition of a {\it small} ferromagnetic next-nearest-neighbor exchange coupling $|J_2|>0.09 J_1$ leads to stabilization of a 36-site unit cell VBC, although the DSL remains a local minimum of the variational parameter landscape. This implies that the VBC is not trivially connected to the DSL: instead it possesses a non-trivial flux pattern and large dimerization.Comment: 5 pages, 4 figure

Projected wave function study of Z2 spin liquids on the kagome lattice for the spin-1/2 quantum Heisenberg antiferromagnet

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 Z2[0,{\pi}]{\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.Comment: 5 pages + supplementary material (2 pages), 3 figures, 1 Table: Final published version, selected as an Editor's suggestio

Spin-12\frac{1}{2} Heisenberg J1J_1-J2J_2 antiferromagnet on the kagome lattice

We report variational Monte Carlo calculations for the spin-$\frac{1}{2}$ Heisenberg model on the kagome lattice in the presence of both nearest-neighbor $J_1$ and next-nearest-neighbor $J_2$ antiferromagnetic superexchange couplings. Our approach is based upon Gutzwiller projected fermionic states that represent a flexible tool to describe quantum spin liquids with different properties (e.g., gapless and gapped). We show that, on finite clusters, a gapped $\mathbb{Z}_{2}$ spin liquid can be stabilized in the presence of a finite $J_2$ superexchange, with a substantial energy gain with respect to the gapless $U(1)$ Dirac spin liquid. However, this energy gain vanishes in the thermodynamic limit, implying that, at least within this approach, the $U(1)$ Dirac spin liquid remains stable in a relatively large region of the phase diagram. For $J_2/J_1 \gtrsim 0.3$, we find that a magnetically ordered state with ${\bf q}={\bf 0}$ overcomes the magnetically disordered wave functions, suggesting the end of the putative gapless spin-liquid phase.Comment: 6 pages, 4 figures. Published versio

Vanishing spin gap in a competing spin-liquid phase in the kagome Heisenberg antiferromagnet

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 Y. 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.Comment: 5 pages, 3 figures, 2 tables. Published versio

Gapless spin-liquid phase in the kagome spin-1/2 Heisenberg antiferromagnet

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=-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.Comment: Final published version. Main Paper (6 pages, 4 figures, 1 table) + Supplementary Material (4 pages, 8 figures