In this paper, we do a complete classification of valence-bond crystals (VBCs) on the kagomé lattice based on general arguments of symmetry only and thus identify many new VBCs for different unit cell sizes. For the spin-1/2 Heisenberg antiferromagnet, we study the relative energetics of competing gapless spin liquids (SLs) and VBC phases within the class of Gutzwiller-projected fermionic wave functions using variational Monte Carlo techniques, hence implementing exactly the constraint of one fermion per site. By using a state-of-the-art optimization method, we conclusively show that the U(1) Dirac SL is remarkably stable towards dimerizing into all 6-, 12- and 36-site unit cell VBCs. This stability is also preserved on addition of a next-nearest-neighbor super-exchange coupling of both antiferromagnetic and ferromagnetic (FM) type. However, we find that a 36-site unit cell VBC is stabilized on addition of a very small next-nearest-neighbor FM super-exchange coupling, i.e. | J 2 | ≈ 0.045, and this VBC is the same in terms of space-group symmetry as that obtained in an effective quantum dimer model study. It breaks reflection symmetry, has a nontrivial flux pattern and is a strong dimerization of the uniform RVB SL.
The low-energy singlet dynamics of the quantum Heisenberg antiferromagnet on the kagome lattice is described by a quantitative quantum dimer model. Using advanced numerical tools, the latter is shown to exhibit valence-bond crystal (VBC) order with a large 36-site unit cell and hidden degeneracy between even odd parities. Evidences are given that this ground state lies in the vicinity of a
Z2 dimer-liquid region separated by a quantum critical point (QCP). Implications regarding numerical analysis and experiments are discussed.
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.
The doping of quantum Heisenberg antiferromagnets on the kagome lattice by nonmagnetic impurities is investigated within the framework of a generalized quantum dimer model (QDM) describing (a) the valence-bond crystal (VBC), (b) the dimer liquid, and (c) the critical region on equal footing. Following the approach by Ralko et al. [Phys. Rev. Lett. 101, 117204 (2008)] for the square and triangular lattices, we introduce the (minimal) extension of the QDM on the kagome lattice to account for spontaneous creation of mobile
S=1/2 spinons at finite magnetic field. Modulations of the dimer density (at zero or finite magnetic field) and of the local field-induced magnetization in the vicinity of impurities are computed using Lanczos exact diagonalization techniques on small clusters (48 and 75 sites). The VBC is clearly revealed from its pinning by impurities, while, in the dimer liquid, crystallization around impurities involves only two neighboring dimers. We also find that a next-nearest-neighbor ferromagnetic coupling favors VBC order. Unexpectedly, a small-size spinon-impurity bound state appears in some region of the dimer-liquid phase. In contrast, in the VBC phase the spinon delocalizes within a large region around the impurity, revealing the weakness of the VBC confining potential. Lastly, we observe that an impurity concentration as small as 4% enhances dimerization substantially. These results are confronted to the valence-bond glass scenario [R. R. P. Singh, Phys. Rev. Lett. 104, 177203 (2010)] and implications to the interpretation of the nuclear magnetic resonance spectra of the Herbertsmithite compound are outlined.
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.
The highly frustrated spin-
21 quantum Heisenberg model with both nearest (
J1) and next-nearest (
J2) 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 small ferromagnetic next-nearest-neighbor exchange coupling
∣J2∣>0.09J1 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 nontrivial flux pattern and large dimerization.