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:2203.10216 · Physical review. B./Physical review. B

Nature of quantum spin liquids of the S = 1 2 Heisenberg antiferromagnet on the triangular lattice: A parallel DMRG study

Yi‐Fan Jiang, Hong‐Chen Jiang

Published Apr 27, 2023• 20 citations
arXivHTMLPDF

Abstract

We study the ground-state properties of the quantum spin liquid (QSL) phases of the spin-$1/2$ antiferromagnetic Heisenberg model on the triangular lattice with nearest- $({J}_{1})$, next-nearest- $({J}_{2})$, and third-neighbor $({J}_{3})$ interactions by using the density-matrix renormalization group (DMRG) method. By combining parallel DMRG with $\mathrm{SU}(2)$ spin rotational symmetry, we are able to obtain accurate results on large cylinders with length up to ${L}_{x}=48$ and circumference ${L}_{y}=6--12$. Our results suggest that the QSL phase of the ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model is gapped, which is characterized by the absence of a gapless mode and by short-range spin-spin and dimer-dimer correlations. In the presence of ${J}_{3}$ interaction, we find that a critical QSL with a single gapless mode emerges. While both spin-spin and scalar chiral-chiral correlations are short ranged, dimer-dimer correlations are quasi long ranged and decay as a power law at long distances.

DOI
10.1103/physrevb.107.l140411

Related papers

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

Competing spin-liquid states in the spin- 1 2 Heisenberg model on the triangular lattice

Wen-Jun Hu, Shou-Shu Gong, Wei Kang Zhu +1 more

We study the spin-$\frac{1}{2}$ Heisenberg model on the triangular lattice with antiferromagnetic first- $({J}_{1})$ and second- $({J}_{2})$ nearest-neighbor interactions using density matrix renormalization group. By studying the spin correlation function, we find a ${120}^{\ensuremath{\circ}}$ magnetic order phase for ${J}_{2}\ensuremath{\lesssim}0.07{J}_{1}$ and a stripe antiferromagnetic phase for ${J}_{2}\ensuremath{\gtrsim}0.15{J}_{1}$. Between these two phases, we identify a spin-liquid region characterized by exponential decaying spin and dimer correlations, as well as large spin singlet and triplet excitation gaps on finite-size systems. We find two near degenerating ground states with distinct properties in two sectors, which indicates more than one spin-liquid candidate in this region. While the sector with spinons is found to respect time reversal symmetry, the even sector without spinons breaks such a symmetry for finite-size systems. Furthermore, we detect the signature of the fractionalization by following the evolution of different ground states with inserting spin flux into the cylinder system. Moreover, by tuning the anisotropic bond coupling, we explore the nature of the spin-liquid phase and find the optimal parameter region for gapped ${Z}_{2}$ spin liquids.

arXiv:1504.00654•Physical Review B•Oct 2, 2015•214 citations

Global phase diagram and quantum spin liquids in a spin- 12 triangular antiferromagnet

Shou-Shu Gong, Wei Kang Zhu, Jian‐Xin Zhu +2 more

We study the spin-$1/2$ Heisenberg model on the triangular lattice with the nearest-neighbor ${J}_{1}>0$, the next-nearest-neighobr ${J}_{2}>0$ Heisenberg interactions, and the additional scalar chiral interaction ${J}_{\ensuremath{\chi}}({\stackrel{P\vec}{S}}_{i}\ifmmode\times\else\texttimes\fi{}{\stackrel{P\vec}{S}}_{j})\ifmmode\cdot\else\textperiodcentered\fi{}{\stackrel{P\vec}{S}}_{k}$ for the three spins in all the triangles using large-scale density matrix renormalization group calculation on cylinder geometry. With increasing ${J}_{2}$ (${J}_{2}/{J}_{1}\ensuremath{\le}0.3$) and ${J}_{\ensuremath{\chi}}$ (${J}_{\ensuremath{\chi}}/{J}_{1}\ensuremath{\le}1.0$) interactions, we establish a quantum phase diagram with the magnetically ordered ${120}^{\ensuremath{\circ}}$, stripe, and noncoplanar tetrahedral phase. In between these magnetic order phases, we find a chiral spin liquid (CSL) phase, which is identified as a $\ensuremath{\nu}=1/2$ bosonic fractional quantum Hall state with possible spontaneous rotational symmetry breaking. By switching on the chiral interaction, we find that the previously identified spin liquid in the ${J}_{1}\ensuremath{-}{J}_{2}$ triangular model ($0.08\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\lesssim}0.15$) shows a phase transition to the CSL phase at very small ${J}_{\ensuremath{\chi}}$. We also compute the spin triplet gap in both spin liquid phases, and our finite-size results suggest a large gap in the odd topological sector but a small or vanishing gap in the even sector. We discuss the implications of our results on the nature of the spin liquid phases.

arXiv:1705.00510•Physical review. B./Physical review. B•Aug 9, 2017•107 citations

Competing States in the S = 1 / 2 Triangular-Lattice J 1 − J 2 Heisenberg Model: A Dynamical Density-Matrix Renormalization Group Study

Shengtao Jiang, Steven R. White, Steven Allan Kivelson +1 more

Previous studies of the S=1/2 triangular-lattice J_{1}-J_{2} Heisenberg antiferromagnet have inferred the existence of a nonmagnetic ground-state phase for an intermediate range of J_{2}, but disagree concerning whether it is a gapped Z_{2} quantum spin liquid (QSL), a gapless (Dirac) QSL, or a weakly symmetry-broken phase. Using an improved dynamical density-matrix renormalization group method, we investigate the relevant intermediate J_{2} regime for cylinders with circumferences from 6 to 9. Depending on the initial state and boundary conditions, we find two distinct variational states. The higher energy state is consistent with a Dirac QSL. In the lower-energy state, both the static and dynamical properties are qualitatively similar to the magnetically ordered state at J_{2}=0, suggestive of either a weakly magnetically ordered non-QSL or a gapped QSL proximate to a continuous transition to such an ordered state.

arXiv:2602.14892•Physical Review Letters•Jun 16, 2026•1 citation

Investigation of the $J_1$-$J_2$ Heisenberg model on the triangular lattice: A study with projected entangled-pair states

L R Ma, Wei-Lin Tu, Didier Poilblanc +1 more

The nature of the quantum spin liquid (QSL) phase in the frustrated $J_1$-$J_2$ Heisenberg model on the triangular lattice remains an open and actively debated problem. In this work, we employ the infinite projected entangled-pair state (PEPS) to systematically investigate the model under different symmetry constraints. Our simulations reveal a direct transition from the $120^\circ$ Néel state to a putative QSL at $J_2/J_1\approx 0.08$, signaled by the collapse of magnetic order. We further show that, through either an appropriate unitary rotation or spontaneous spin long-range order, the stripe antiferromagnetic phase can also be accurately captured within the infinite PEPS framework. A central focus of our study is the role played by the PEPS symmetry in approximating the QSL ground-state sandwiched between the two magnetic phases. We first found that a fully-symmetric topological $\mathbb{Z}_2$ Resonating Valence Bond state, which can be written as a simple PEPS with bond dimension $D=3$, exhibits a reasonably good variational energy. Motivated by this finding, we have further constructed generic $\mathbb{Z}_2$-symmetric PEPS of larger bond dimension (up to $D=7$). We found that, under wavefunction optimization, spinons condense and, simultaneously, topological vison excitations get confined, hence precluding $\mathbb{Z}_2$ topological order. This strongly indicates the gapless (or critical) nature of the QSL phase, which is most naturally consistent with a U(1) Dirac spin liquid scenario.

arXiv:2606.31021•Jun 30, 2026•0 citations

Gapless spin liquid ground state of the spin- 12 J1−J2 Heisenberg model on square lattices

Wen-Yuan Liu, Shaojun Dong, Chao Wang +4 more

The spin-$\frac{1}{2}$ ${J}_{1}\ensuremath{-}{J}_{2}$ Heisenberg model on square lattices is investigated via the finite projected entangled pair states (PEPS) method. Using the recently developed gradient optimization method combined with Monte Carlo sampling techniques, we are able to obtain the ground state energies that are competitive with the best results. The calculations show that there is no N\'eel order, dimer order, or plaquette order in the region of $0.42\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\lesssim}0.6$, suggesting a single spin liquid phase in the intermediate region. Furthermore, the calculated staggered spin, dimer, and plaquette correlation functions all have power law decay behaviors, which provide strong evidence that the intermediate nonmagnetic phase is a single gapless spin liquid state.

arXiv:1806.07031•Physical review. B./Physical review. B•Dec 17, 2018•71 citations
We study the ground-state properties of the quantum spin liquid (QSL) phases of the spin-1/21/21/2 antiferromagnetic Heisenberg model on the triangular lattice with nearest- (J1)({J}_{1})(J1​), next-nearest- (J2)({J}_{2})(J2​), and third-neighbor (J3)({J}_{3})(J3​) interactions by using the density-matrix renormalization group (DMRG) method. By combining parallel DMRG with SU(2)\mathrm{SU}(2)SU(2) spin rotational symmetry, we are able to obtain accurate results on large cylinders with length up to Lx=48{L}_{x}=48Lx​=48 and circumference Ly=6−−12{L}_{y}=6--12Ly​=6−−12. Our results suggest that the QSL phase of the J1\ensuremath-J2{J}_{1}\text{\ensuremath{-}}{J}_{2}J1​\ensuremath-J2​ Heisenberg model is gapped, which is characterized by the absence of a gapless mode and by short-range spin-spin and dimer-dimer correlations. In the presence of J3{J}_{3}J3​ interaction, we find that a critical QSL with a single gapless mode emerges. While both spin-spin and scalar chiral-chiral correlations are short ranged, dimer-dimer correlations are quasi long ranged and decay as a power law at long distances.
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.
We study the spin-12\frac{1}{2}21​ Heisenberg model on the triangular lattice with antiferromagnetic first- (J1)({J}_{1})(J1​) and second- (J2)({J}_{2})(J2​) nearest-neighbor interactions using density matrix renormalization group. By studying the spin correlation function, we find a ${120}^{\ensuremath{\circ}}$ magnetic order phase for ${J}_{2}\ensuremath{\lesssim}0.07{J}_{1}$ and a stripe antiferromagnetic phase for ${J}_{2}\ensuremath{\gtrsim}0.15{J}_{1}$. Between these two phases, we identify a spin-liquid region characterized by exponential decaying spin and dimer correlations, as well as large spin singlet and triplet excitation gaps on finite-size systems. We find two near degenerating ground states with distinct properties in two sectors, which indicates more than one spin-liquid candidate in this region. While the sector with spinons is found to respect time reversal symmetry, the even sector without spinons breaks such a symmetry for finite-size systems. Furthermore, we detect the signature of the fractionalization by following the evolution of different ground states with inserting spin flux into the cylinder system. Moreover, by tuning the anisotropic bond coupling, we explore the nature of the spin-liquid phase and find the optimal parameter region for gapped Z2{Z}_{2}Z2​ spin liquids.
We study the spin-1/21/21/2 Heisenberg model on the triangular lattice with the nearest-neighbor J1>0{J}_{1}>0J1​>0, the next-nearest-neighobr J2>0{J}_{2}>0J2​>0 Heisenberg interactions, and the additional scalar chiral interaction ${J}_{\ensuremath{\chi}}({\stackrel{P\vec}{S}}_{i}\ifmmode\times\else\texttimes\fi{}{\stackrel{P\vec}{S}}_{j})\ifmmode\cdot\else\textperiodcentered\fi{}{\stackrel{P\vec}{S}}_{k}$ for the three spins in all the triangles using large-scale density matrix renormalization group calculation on cylinder geometry. With increasing J2{J}_{2}J2​ (${J}_{2}/{J}_{1}\ensuremath{\le}0.3$) and ${J}_{\ensuremath{\chi}}$ (${J}_{\ensuremath{\chi}}/{J}_{1}\ensuremath{\le}1.0$) interactions, we establish a quantum phase diagram with the magnetically ordered ${120}^{\ensuremath{\circ}}$, stripe, and noncoplanar tetrahedral phase. In between these magnetic order phases, we find a chiral spin liquid (CSL) phase, which is identified as a $\ensuremath{\nu}=1/2$ bosonic fractional quantum Hall state with possible spontaneous rotational symmetry breaking. By switching on the chiral interaction, we find that the previously identified spin liquid in the ${J}_{1}\ensuremath{-}{J}_{2}$ triangular model ($0.08\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\lesssim}0.15$) shows a phase transition to the CSL phase at very small ${J}_{\ensuremath{\chi}}$. We also compute the spin triplet gap in both spin liquid phases, and our finite-size results suggest a large gap in the odd topological sector but a small or vanishing gap in the even sector. We discuss the implications of our results on the nature of the spin liquid phases.
Previous studies of the S=1/2 triangular-lattice J_{1}-J_{2} Heisenberg antiferromagnet have inferred the existence of a nonmagnetic ground-state phase for an intermediate range of J_{2}, but disagree concerning whether it is a gapped Z_{2} quantum spin liquid (QSL), a gapless (Dirac) QSL, or a weakly symmetry-broken phase. Using an improved dynamical density-matrix renormalization group method, we investigate the relevant intermediate J_{2} regime for cylinders with circumferences from 6 to 9. Depending on the initial state and boundary conditions, we find two distinct variational states. The higher energy state is consistent with a Dirac QSL. In the lower-energy state, both the static and dynamical properties are qualitatively similar to the magnetically ordered state at J_{2}=0, suggestive of either a weakly magnetically ordered non-QSL or a gapped QSL proximate to a continuous transition to such an ordered state.
Investigation of the J1J_1J1​-J2J_2J2​ Heisenberg model on the triangular lattice: A study with projected entangled-pair states
The nature of the quantum spin liquid (QSL) phase in the frustrated J1J_1J1​-J2J_2J2​ Heisenberg model on the triangular lattice remains an open and actively debated problem. In this work, we employ the infinite projected entangled-pair state (PEPS) to systematically investigate the model under different symmetry constraints. Our simulations reveal a direct transition from the 120∘120^\circ120∘ Néel state to a putative QSL at J2/J1≈0.08J_2/J_1\approx 0.08J2​/J1​≈0.08, signaled by the collapse of magnetic order. We further show that, through either an appropriate unitary rotation or spontaneous spin long-range order, the stripe antiferromagnetic phase can also be accurately captured within the infinite PEPS framework. A central focus of our study is the role played by the PEPS symmetry in approximating the QSL ground-state sandwiched between the two magnetic phases. We first found that a fully-symmetric topological Z2\mathbb{Z}_2Z2​ Resonating Valence Bond state, which can be written as a simple PEPS with bond dimension D=3D=3D=3, exhibits a reasonably good variational energy. Motivated by this finding, we have further constructed generic Z2\mathbb{Z}_2Z2​-symmetric PEPS of larger bond dimension (up to D=7D=7D=7). We found that, under wavefunction optimization, spinons condense and, simultaneously, topological vison excitations get confined, hence precluding Z2\mathbb{Z}_2Z2​ topological order. This strongly indicates the gapless (or critical) nature of the QSL phase, which is most naturally consistent with a U(1) Dirac spin liquid scenario.
The spin-12\frac{1}{2}21​ ${J}_{1}\ensuremath{-}{J}_{2}$ Heisenberg model on square lattices is investigated via the finite projected entangled pair states (PEPS) method. Using the recently developed gradient optimization method combined with Monte Carlo sampling techniques, we are able to obtain the ground state energies that are competitive with the best results. The calculations show that there is no N\'eel order, dimer order, or plaquette order in the region of $0.42\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\lesssim}0.6$, suggesting a single spin liquid phase in the intermediate region. Furthermore, the calculated staggered spin, dimer, and plaquette correlation functions all have power law decay behaviors, which provide strong evidence that the intermediate nonmagnetic phase is a single gapless spin liquid state.