Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1505.06276 · Physical Review B

Variational Monte Carlo study of a gapless spin liquid in the spin- 12 XXZ antiferromagnetic model on the kagome lattice

Wen-Jun Hu, Shou-Shu Gong, Federico Becca, D. N. Sheng

Published Nov 11, 2015• 6 citations
arXivHTMLPDF

Abstract

By using the variational Monte Carlo technique, we study the spin-$\frac{1}{2}$ XXZ antiferromagnetic model (with easy-plane anisotropy) on the kagome lattice. A class of Gutzwiller projected fermionic states with a spin Jastrow factor is considered to describe either spin liquids [with $U(1)$ or ${Z}_{2}$ symmetry] or magnetically ordered phases [with $\mathbf{q}=(0,0)$ or $\mathbf{q}=(4\ensuremath{\pi}/3,0)$]. We find that the magnetic states are not stable in the thermodynamic limit. Moreover, there is no energy gain to break the gauge symmetry from $U(1)$ to ${Z}_{2}$ within the spin-liquid states, as previously found in the Heisenberg model. The best variational wave function is therefore the $U(1)$ Dirac state, supplemented by the spin Jastrow factor. Furthermore, a vanishing $S=2$ spin gap is obtained at the variational level, in the whole regime from the $XY$ to the Heisenberg model.

DOI
10.1103/physrevb.92.201105

Related papers

Gapless spin-liquid phase in the kagome spin- 12 Heisenberg antiferromagnet

Yasir Iqbal, Federico Becca, Sandro Sorella +1 more

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) ${\mathbb{Z}}_{2}$ spin liquids. By performing a finite-size extrapolation of the ground-state energy, we obtain an energy per site $E/J=\ensuremath{-}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.

arXiv:1209.1858•Physical Review B•Feb 14, 2013•331 citations

Persistence of the gapless spin liquid in the breathing kagome Heisenberg antiferromagnet

Yasir Iqbal, Didier Poilblanc, Ronny Thomale +1 more

The nature of the ground state of the spin $S=1/2$ Heisenberg antiferromagnet on the kagome lattice with breathing anisotropy (i.e., with different superexchange couplings ${J}_{▵}$ and ${J}_{▿}$ within elementary up- and down-pointing triangles) is investigated within the framework of Gutzwiller projected fermionic wave functions and Monte Carlo methods. We analyze the stability of the $\text{U}(1)$ Dirac spin liquid with respect to the presence of fermionic pairing that leads to a gapped ${\mathbb{Z}}_{2}$ spin liquid. For several values of the ratio ${J}_{▿}/{J}_{▵}$, the size scaling of the energy gain due to the pairing fields and the variational parameters are reported. Our results show that the energy gain of the gapped spin liquid with respect to the gapless state either vanishes for large enough system size or scales to zero in the thermodynamic limit. Similarly, the optimized pairing amplitudes (responsible for opening the spin gap) are shown to vanish in the thermodynamic limit. Our outcome is corroborated by the application of one and two Lanczos steps to the gapless and gapped wave functions, for which no energy gain of the gapped state is detected when improving the quality of the variational states. Finally, we discuss the competition with the ``simplex'' ${\mathbb{Z}}_{2}$ resonating-valence-bond spin liquid, valence-bond crystal, and nematic states in the strongly anisotropic regime, i.e., ${J}_{▿}\ensuremath{\ll}{J}_{▵}$.

arXiv:1712.04579•Physical review. B./Physical review. B•Mar 13, 2018•40 citations

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

Yasir Iqbal, Federico Becca, Didier Poilblanc

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 ${\mathbb{Z}}_{2}$ 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 ${\mathbb{Z}}_{2}$ 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 ${\mathbb{Z}}_{2}$ 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 ${\mathbb{Z}}_{2}$ 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.

arXiv:1105.0341•Physical Review B•Jul 20, 2011•172 citations

Variational Monte Carlo study of a chiral spin liquid in the extended Heisenberg model on the kagome lattice

Wen-Jun Hu, Wei Kang Zhu, Yi Zhang +3 more

We investigate the extended Heisenberg model on the kagome lattice by using Gutzwiller projected fermionic states and the variational Monte Carlo technique. In particular, when both second- and third-neighbor superexchanges are considered, we find that a gapped spin liquid described by nontrivial magnetic fluxes and long-range chiral-chiral correlations is energetically favored compared to the gapless U(1) Dirac state. Furthermore, the topological Chern number, obtained by integrating the Berry curvature, and the degeneracy of the ground state, by constructing linearly independent states, lead us to identify this flux state as the chiral spin liquid with a $C=1/2$ fractionalized Chern number.

arXiv:1411.1327•Physical Review B•Jan 30, 2015•78 citations

Comment on "$Z_{2}$ spin liquid phase on the kagome lattice: a new saddle point", by Tao Li [arXiv:1601.02165 (2016)]

Yasir Iqbal, Didier Poilblanc, Federico Becca

In a recent paper [arXiv:1601.02165], Tao Li claimed that a gapped $\mathbb{Z}_2$ spin liquid, obtained using projected Gutzwiller fermionic wave functions, can be stabilized in the Heisenberg model on the kagome lattice. We perform very accurate variational calculations and confirm with unprecedented accuracy the fact, that the $\mathbb{Z}_{2}$ spin liquid is a local energy minimum that goes away with system size, thus reaffirming the scenario of a gapless $U(1)$ Dirac spin liquid ground state.

arXiv:1606.02255•Jun 5, 2016•6 citations

Variational Monte Carlo study of chiral spin liquid in quantum antiferromagnet on the triangular lattice

Wen-Jun Hu, Shou-Shu Gong, D. N. Sheng

By using Gutzwiller projected fermionic wave functions and variational Monte Carlo technique, we study the spin-$1/2$ Heisenberg model with the first-neighbor (${J}_{1}$), second-neighbor (${J}_{2}$), and additional scalar chiral interaction ${J}_{\ensuremath{\chi}}{\mathbf{S}}_{i}\ifmmode\cdot\else\textperiodcentered\fi{}({\mathbf{S}}_{j}\ifmmode\times\else\texttimes\fi{}{\mathbf{S}}_{k})$ on the triangular lattice. In the nonmagnetic phase of the ${J}_{1}\ensuremath{-}{J}_{2}$ triangular model with $0.08\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\lesssim}0.16$, recent density-matrix renormalization group (DMRG) studies [Zhu and White, Phys. Rev. B 92, 041105(R) (2015) and Hu, Gong, Zhu, and Sheng, Phys. Rev. B 92, 140403(R) (2015)] find a possible gapped spin liquid with the signal of a competition between a chiral and a ${Z}_{2}$ spin liquid. Motivated by the DMRG results, we consider the chiral interaction ${J}_{\ensuremath{\chi}}{\mathbf{S}}_{i}\ifmmode\cdot\else\textperiodcentered\fi{}({\mathbf{S}}_{j}\ifmmode\times\else\texttimes\fi{}{\mathbf{S}}_{k})$ as a perturbation for this nonmagnetic phase. We find that with growing ${J}_{\ensuremath{\chi}}$, the gapless U(1) Dirac spin liquid, which has the best variational energy for ${J}_{\ensuremath{\chi}}=0$, exhibits the energy instability towards a gapped spin liquid with nontrivial magnetic fluxes and nonzero chiral order. We calculate topological Chern number and ground-state degeneracy, both of which identify this flux state as the chiral spin liquid with fractionalized Chern number $C=1/2$ and twofold topological degeneracy. Our results indicate a positive direction to stabilize a chiral spin liquid near the nonmagnetic phase of the ${J}_{1}\ensuremath{-}{J}_{2}$ triangular model.

arXiv:1603.03365•Physical review. B./Physical review. B•Aug 15, 2016•64 citations
By using the variational Monte Carlo technique, we study the spin-12\frac{1}{2}21​ XXZ antiferromagnetic model (with easy-plane anisotropy) on the kagome lattice. A class of Gutzwiller projected fermionic states with a spin Jastrow factor is considered to describe either spin liquids [with U(1)U(1)U(1) or Z2{Z}_{2}Z2​ symmetry] or magnetically ordered phases [with q=(0,0)\mathbf{q}=(0,0)q=(0,0) or q=(4\ensuremathπ/3,0)\mathbf{q}=(4\ensuremath{\pi}/3,0)q=(4\ensuremathπ/3,0)]. We find that the magnetic states are not stable in the thermodynamic limit. Moreover, there is no energy gain to break the gauge symmetry from U(1)U(1)U(1) to Z2{Z}_{2}Z2​ within the spin-liquid states, as previously found in the Heisenberg model. The best variational wave function is therefore the U(1)U(1)U(1) Dirac state, supplemented by the spin Jastrow factor. Furthermore, a vanishing S=2S=2S=2 spin gap is obtained at the variational level, in the whole regime from the XYXYXY to the Heisenberg model.
We study the energy and the static spin structure factor of the ground state of the spin-1/21/21/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)U(1)U(1) Dirac spin liquid is competitive with previously proposed gapped (topological) Z2{\mathbb{Z}}_{2}Z2​ spin liquids. By performing a finite-size extrapolation of the ground-state energy, we obtain an energy per site $E/J=\ensuremath{-}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.
The nature of the ground state of the spin S=1/2S=1/2S=1/2 Heisenberg antiferromagnet on the kagome lattice with breathing anisotropy (i.e., with different superexchange couplings J▵{J}_{▵}J▵​ and J▿{J}_{▿}J▿​ within elementary up- and down-pointing triangles) is investigated within the framework of Gutzwiller projected fermionic wave functions and Monte Carlo methods. We analyze the stability of the U(1)\text{U}(1)U(1) Dirac spin liquid with respect to the presence of fermionic pairing that leads to a gapped Z2{\mathbb{Z}}_{2}Z2​ spin liquid. For several values of the ratio J▿/J▵{J}_{▿}/{J}_{▵}J▿​/J▵​, the size scaling of the energy gain due to the pairing fields and the variational parameters are reported. Our results show that the energy gain of the gapped spin liquid with respect to the gapless state either vanishes for large enough system size or scales to zero in the thermodynamic limit. Similarly, the optimized pairing amplitudes (responsible for opening the spin gap) are shown to vanish in the thermodynamic limit. Our outcome is corroborated by the application of one and two Lanczos steps to the gapless and gapped wave functions, for which no energy gain of the gapped state is detected when improving the quality of the variational states. Finally, we discuss the competition with the ``simplex'' Z2{\mathbb{Z}}_{2}Z2​ resonating-valence-bond spin liquid, valence-bond crystal, and nematic states in the strongly anisotropic regime, i.e., ${J}_{▿}\ensuremath{\ll}{J}_{▵}$.
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/21/21/2 Heisenberg antiferromagnet on the kagome lattice geometry is a fully gapped spin liquid with numerical signatures of Z2{\mathbb{Z}}_{2}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{\mathbb{Z}}_{2}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{\mathbb{Z}}_{2}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{\mathbb{Z}}_{2}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 investigate the extended Heisenberg model on the kagome lattice by using Gutzwiller projected fermionic states and the variational Monte Carlo technique. In particular, when both second- and third-neighbor superexchanges are considered, we find that a gapped spin liquid described by nontrivial magnetic fluxes and long-range chiral-chiral correlations is energetically favored compared to the gapless U(1) Dirac state. Furthermore, the topological Chern number, obtained by integrating the Berry curvature, and the degeneracy of the ground state, by constructing linearly independent states, lead us to identify this flux state as the chiral spin liquid with a C=1/2C=1/2C=1/2 fractionalized Chern number.
Comment on "Z2Z_{2}Z2​ spin liquid phase on the kagome lattice: a new saddle point", by Tao Li [arXiv:1601.02165 (2016)]
In a recent paper [arXiv:1601.02165], Tao Li claimed that a gapped Z2\mathbb{Z}_2Z2​ spin liquid, obtained using projected Gutzwiller fermionic wave functions, can be stabilized in the Heisenberg model on the kagome lattice. We perform very accurate variational calculations and confirm with unprecedented accuracy the fact, that the Z2\mathbb{Z}_{2}Z2​ spin liquid is a local energy minimum that goes away with system size, thus reaffirming the scenario of a gapless U(1)U(1)U(1) Dirac spin liquid ground state.
By using Gutzwiller projected fermionic wave functions and variational Monte Carlo technique, we study the spin-1/21/21/2 Heisenberg model with the first-neighbor (J1{J}_{1}J1​), second-neighbor (J2{J}_{2}J2​), and additional scalar chiral interaction ${J}_{\ensuremath{\chi}}{\mathbf{S}}_{i}\ifmmode\cdot\else\textperiodcentered\fi{}({\mathbf{S}}_{j}\ifmmode\times\else\texttimes\fi{}{\mathbf{S}}_{k})$ on the triangular lattice. In the nonmagnetic phase of the ${J}_{1}\ensuremath{-}{J}_{2}$ triangular model with $0.08\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\lesssim}0.16$, recent density-matrix renormalization group (DMRG) studies [Zhu and White, Phys. Rev. B 92, 041105(R) (2015) and Hu, Gong, Zhu, and Sheng, Phys. Rev. B 92, 140403(R) (2015)] find a possible gapped spin liquid with the signal of a competition between a chiral and a Z2{Z}_{2}Z2​ spin liquid. Motivated by the DMRG results, we consider the chiral interaction ${J}_{\ensuremath{\chi}}{\mathbf{S}}_{i}\ifmmode\cdot\else\textperiodcentered\fi{}({\mathbf{S}}_{j}\ifmmode\times\else\texttimes\fi{}{\mathbf{S}}_{k})$ as a perturbation for this nonmagnetic phase. We find that with growing ${J}_{\ensuremath{\chi}}$, the gapless U(1) Dirac spin liquid, which has the best variational energy for ${J}_{\ensuremath{\chi}}=0$, exhibits the energy instability towards a gapped spin liquid with nontrivial magnetic fluxes and nonzero chiral order. We calculate topological Chern number and ground-state degeneracy, both of which identify this flux state as the chiral spin liquid with fractionalized Chern number C=1/2C=1/2C=1/2 and twofold topological degeneracy. Our results indicate a positive direction to stabilize a chiral spin liquid near the nonmagnetic phase of the ${J}_{1}\ensuremath{-}{J}_{2}$ triangular model.