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 singlet dynamics, which plays a major role in the physics of the spin-1/2 quantum Heisenberg antiferromagnet (QHAF) on the kagome lattice, can be approximately described by projecting onto the nearest-neighbor valence bond (NNVB) singlet subspace. We revisit here the effective quantum dimer model, which originates from the latter NNVB-projected Heisenberg model via a nonperturbative Rokhsar-Kivelson-like scheme. By using Lanczos exact diagonalization on a 108-site cluster supplemented by a careful symmetry analysis, it is shown that a previously found 36-site valence bond crystal (VBC) in fact competes with a new type of 12-site ``resonating-columnar'' VBC. Interestingly, these two VBCs ``emerge'' in different topological sectors. The exceptionally large degeneracy of the ground-state multiplets (144 on our 108-site cluster) and the proximity of a
Z2 dimer liquid have implications for the interpretation of numerical results on the QHAF, which are outlined. The possibility of a chiral VBC (i.e., spontaneously breaking time-reversal symmetry) is also discussed.
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.
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.
We examine the perennial quantum spin liquid candidate
S=21 Heisenberg antiferromagnet on the kagome lattice. Our paper is based on achieving Lanczos diagonalization of the Hamiltonian on a 48 site cluster in sectors with dimensions as large as $5\ifmmode\times\else\texttimes\fi{}{10}^{11}$. The results reveal intricate structures in the low-lying energy spectrum. These structures by no means unambiguously support a
Z2 spin liquid ground state, but instead appear compatible with several scenarios, including fourfold topological degeneracy, inversion symmetry breaking, and a combination thereof. We discuss finite-size effects, such as the apparent absence of the eigenstate thermalization hypothesis, and note that, while considerably reduced, some are still present for the largest cluster. Finally, we observe that an XXZ model in the Ising limit reproduces remarkably well the most striking features of finite-size spectra.
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 order with a large 36-site unit cell and hidden degeneracy between even and odd parities. Evidences are given that this groundstate lies in the vicinity of a
Z2 dimer liquid region separated by a Quantum Critical Point. Implications regarding numerical analysis and experiments are discussed.