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

Magnetization of the spin- 12 Heisenberg antiferromagnet on the triangular lattice

Qian Li, Li Hong, Jize Zhao, Hong‐Gang Luo, Z. Y. Xie

Published May 16, 2022• 23 citations
arXivHTMLPDF

Abstract

After decades of debate, now there is a rough consensus that at zero temperature the spin-$1/2$ Heisenberg antiferromagnet on the triangular lattice is three-sublattice ${120}^{\ensuremath{\circ}}$ magnetically ordered, in contrast to a quantum spin liquid as originally proposed. However, there remains considerable discrepancy in the magnetization reported among various methods. To resolve this issue, in this work we revisit this model by the tensor-network state algorithm. The ground-state energy per bond ${E}_{b}$ and magnetization per spin ${M}_{0}$ in the thermodynamic limit are obtained with high precision. The former is estimated to be ${E}_{b}=\ensuremath{-}0.18334(10)$. This value agrees well with that from the series expansion. The three-sublattice magnetic order is firmly confirmed and the magnetization is determined as ${M}_{0}=0.161(5)$. It is about $32%$ of its classical value and slightly below the lower bound from the series expansion. In comparison with the best estimated value by Monte Carlo and density-matrix renormalization group, our result is about $20%$ smaller. This magnetic order is consistent with further analysis of the three-body correlation. Our work thus provides benchmark results for this prototypical model.

DOI
10.1103/physrevb.105.184418

Related papers

On the magnetization of the 120° order of the spin-1/2 triangular lattice Heisenberg model: a DMRG revisited

Jiale Huang, Xiangjian Qian, Mingpu Qin

We revisit the issue about the magnetization of the 120° order in the spin-1/2 triangular lattice Heisenberg model with density matrix renormalization group (DMRG). The accurate determination of the magnetization of this model is challenging for numerical methods and its value exhibits substantial disparities across various methods. We perform a large-scale DMRG calculation of this model by employing bond dimension as large asD=24000and by studying the system with width as large asLy=12. With careful extrapolation with truncation error and suitable finite size scaling, we give a conservative estimation of the magnetization asM0=0.208(8). The ground state energy per site we obtain isEg=-0.5503(8). Our results provide valuable benchmark values for the development of new methods in the future.

arXiv:2310.11774•Journal of Physics Condensed Matter•Jan 23, 2024•9 citations

Triangular J 1 − J 2 Heisenberg antiferromagnet in a magnetic field

Thomas Bader, Shi Feng, Sasank Budaraju +3 more

The behavior of the paradigmatic J 1 − J 2 triangular lattice Heisenberg antiferromagnet in a magnetic field remains unsettled despite decades of study. We map out the phase diagram using three complementary approaches, including self-consistent nonlinear spin-wave theory, density-matrix renormalization group, and variational Monte Carlo. This combined analysis resolves the competition among different field-induced magnetic orders and magnetization plateaus across the classically frustrated parameter range. In particular, there is a finite range in the parameter regime around J 2 / J 1 = 1 8 in which (i) upon the application of the external field, the gapless quantum spin liquid acquires a finite density of monopoles, and (ii) by further increasing the field, two plateaus are clearly obtained at m = 1 3 and m = 1 2 . We discuss the experimental importance of the consecutive magnetization plateau transitions as a signature of an underlying quantum spin-liquid phase.

arXiv:2512.08768•Physical review. B./Physical review. B•May 22, 2026•3 citations

Quasiclassical magnetic order and its loss in a spin- 1 2 Heisenberg antiferromagnet on a triangular lattice with competing bonds

P. H. Y. Li, R. F. Bishop, C. E. Campbell

We use the coupled cluster method (CCM) to study the zero-temperature ground-state (GS) properties of a spin-$\frac{1}{2}$ ${J}_{1}\ensuremath{-}{J}_{2}$ Heisenberg antiferromagnet on a triangular lattice with competing nearest-neighbor and next-nearest-neighbor exchange couplings ${J}_{1}>0$ and ${J}_{2}\ensuremath{\equiv}\ensuremath{\kappa}{J}_{1}>0$, respectively, in the window $0\ensuremath{\le}\ensuremath{\kappa}<1$. The classical version of the model has a single GS phase transition at ${\ensuremath{\kappa}}^{\mathrm{cl}}=\frac{1}{8}$ in this window from a phase with 3-sublattice antiferromagnetic (AFM) ${120}^{\ensuremath{\circ}}$ N\'eel order for $\ensuremath{\kappa}<{\ensuremath{\kappa}}^{\mathrm{cl}}$ to an infinitely degenerate family of 4-sublattice AFM N\'eel phases for $\ensuremath{\kappa}>{\ensuremath{\kappa}}^{\mathrm{cl}}$. This classical accidental degeneracy is lifted by quantum fluctuations, which favor a 2-sublattice AFM striped phase. For the quantum model we work directly in the thermodynamic limit of an infinite number of spins, with no consequent need for any finite-size scaling analysis of our results. We perform high-order CCM calculations within a well-controlled hierarchy of approximations, which we show how to extrapolate to the exact limit. In this way we find results for the case $\ensuremath{\kappa}=0$ of the spin-$\frac{1}{2}$ model for the GS energy per spin, $E/N=\ensuremath{-}0.5521(2){J}_{1}$, and the GS magnetic order parameter, $M=0.198(5)$ (in units where the classical value is ${M}^{\mathrm{cl}}=\frac{1}{2}$), which are among the best available. For the spin-$\frac{1}{2}$ ${J}_{1}\ensuremath{-}{J}_{2}$ model we find that the classical transition at $\ensuremath{\kappa}={\ensuremath{\kappa}}^{\mathrm{cl}}$ is split into two quantum phase transitions at ${\ensuremath{\kappa}}_{1}^{c}=0.060(10)$ and ${\ensuremath{\kappa}}_{2}^{c}=0.165(5)$. The two quasiclassical AFM states (viz., the ${120}^{\ensuremath{\circ}}$ N\'eel state and the striped state) are found to be the stable GS phases in the regime $\ensuremath{\kappa}<{\ensuremath{\kappa}}_{1}^{c}$ and $\ensuremath{\kappa}>{\ensuremath{\kappa}}_{2}^{c}$, respectively, while in the intermediate regimes ${\ensuremath{\kappa}}_{1}^{c}<\ensuremath{\kappa}<{\ensuremath{\kappa}}_{2}^{c}$ the stable GS phase has no evident long-range magnetic order.

arXiv:1410.6003•Physical Review B•Jan 22, 2015•86 citations

Ground-state phases of the spin- 1 J 1 − J 2 Heisenberg antiferromagnet on the honeycomb lattice

P. H. Y. Li, R. F. Bishop

We study the zero-temperature quantum phase diagram of a spin-1 Heisenberg antiferromagnet on the honeycomb lattice with both nearest-neighbor exchange coupling ${J}_{1}>0$ and frustrating next-nearest-neighbor coupling ${J}_{2}\ensuremath{\equiv}\ensuremath{\kappa}{J}_{1}>0$, using the coupled cluster method implemented to high orders of approximation, and based on model states with different forms of classical magnetic order. For each we calculate directly in the bulk thermodynamic limit both ground-state low-energy parameters (including the energy per spin, magnetic order parameter, spin stiffness coefficient, and zero-field uniform transverse magnetic susceptibility) and their generalized susceptibilities to various forms of valence-bond crystalline (VBC) order, as well as the energy gap to the lowest-lying spin-triplet excitation. In the range $0<\ensuremath{\kappa}<1$ we find evidence for four distinct phases. Two of these are quasiclassical phases with antiferromagnetic long-range order, one with two-sublattice N\'eel order for $\ensuremath{\kappa}<{\ensuremath{\kappa}}_{{c}_{1}}=0.250(5)$, and another with four-sublattice N\'eel-II order for $\ensuremath{\kappa}>{\ensuremath{\kappa}}_{{c}_{2}}=0.340(5)$. Two different paramagnetic phases are found to exist in the intermediate region. Over the range ${\ensuremath{\kappa}}_{{c}_{1}}<\ensuremath{\kappa}<{\ensuremath{\kappa}}_{c}^{i}=0.305(5)$ we find a gapless phase with no discernible magnetic order, which is a strong candidate for being a quantum spin liquid, while over the range ${\ensuremath{\kappa}}_{c}^{i}<\ensuremath{\kappa}<{\ensuremath{\kappa}}_{{c}_{2}}$ we find a gapped phase, which is most likely a lattice nematic with staggered dimer VBC order that breaks the lattice rotational symmetry.

arXiv:1602.08915•Physical review. B./Physical review. B•Jun 28, 2016•21 citations

Spin liquid phase of the S = 1 2 J 1 − J 2 Heisenberg model on the triangular lattice

Zhenyue Zhu, Steven R. White

We study the $S=1/2$ Heisenberg model on the triangular lattice with nearest- and next-nearest-neighbor interactions ${J}_{1}$ and ${J}_{2}$ with the density matrix renormalization group, on long open cylinders with widths up to nine lattice spacings. In an intermediate ${J}_{2}$ region $0.06\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\lesssim}0.17$, we find evidence for a spin liquid (SL) state with short range spin-spin, bond-bond, and chiral correlation lengths, bordered by a classical ${120}^{\ensuremath{\circ}}$ N\'eel ordered state at small ${J}_{2}$ and by a two sublattice collinear magnetically ordered state at larger ${J}_{2}$. Focusing on ${J}_{2}/{J}_{1}=0.1$, we find a number of signatures of a gapped SL phase.

arXiv:1502.04831•Physical Review B•Jul 7, 2015•275 citations

Nontrivial three-sublattice magnetization in the easy-axis spin- 1 2 XXZ antiferromagnet on the triangular lattice

Masahiro Kadosawa, Masaaki Nakamura, Y. Ohta +1 more

We investigate the ground-state magnetic structure of the spin- 1 / 2 XXZ antiferromagnet on the triangular lattice in the easy-axis regime using the density-matrix renormalization group. By applying spiral boundary conditions, we exactly map finite L × L clusters onto one-dimensional chains while avoiding the spatial anisotropy inherent in cylindrical geometries. From symmetry-broken local magnetization profiles, we extract the three-sublattice moments and track their evolution with anisotropy. At the isotropic point, we obtain a positive sublattice moment of 0.217(3), consistent with previous numerical estimates. In the easy-axis regime ( Δ = J z / J ⊥ > 1 ), the ordered moments remain close to a Y-like zero-magnetization three-sublattice state, whose z -component pattern is of the form ( 2 m , − m , − m ) , over a broad range of Δ . Extrapolation in 1 / Δ shows that the positive sublattice moment stays well below the classical saturation value 1 / 2 , approaching 0.419(7) as Δ → ∞ , while the magnitude of the negative sublattice moment approaches 0.209(4). We further compare the energies of the Y state and the up-down-down state and find that the Y state is favored at zero field. Independent thermodynamic-limit energy calculations, performed without assuming any particular ordered pattern, yield an energy consistent with the Y-state solution. These results show that the easy-axis ground state does not simply cross over to a trivially saturated collinear Ising state, but instead remains a nontrivial three-sublattice ordered state selected from the macroscopically degenerate Ising manifold by quantum fluctuations.

arXiv:2604.14767•Physical review. B./Physical review. B•Jun 12, 2026•1 citation
After decades of debate, now there is a rough consensus that at zero temperature the spin-1/21/21/2 Heisenberg antiferromagnet on the triangular lattice is three-sublattice 120\ensuremath∘{120}^{\ensuremath{\circ}}120\ensuremath∘ magnetically ordered, in contrast to a quantum spin liquid as originally proposed. However, there remains considerable discrepancy in the magnetization reported among various methods. To resolve this issue, in this work we revisit this model by the tensor-network state algorithm. The ground-state energy per bond Eb{E}_{b}Eb​ and magnetization per spin M0{M}_{0}M0​ in the thermodynamic limit are obtained with high precision. The former is estimated to be Eb=\ensuremath−0.18334(10){E}_{b}=\ensuremath{-}0.18334(10)Eb​=\ensuremath−0.18334(10). This value agrees well with that from the series expansion. The three-sublattice magnetic order is firmly confirmed and the magnetization is determined as M0=0.161(5){M}_{0}=0.161(5)M0​=0.161(5). It is about 3232%32 of its classical value and slightly below the lower bound from the series expansion. In comparison with the best estimated value by Monte Carlo and density-matrix renormalization group, our result is about 2020%20 smaller. This magnetic order is consistent with further analysis of the three-body correlation. Our work thus provides benchmark results for this prototypical model.
We revisit the issue about the magnetization of the 120° order in the spin-1/2 triangular lattice Heisenberg model with density matrix renormalization group (DMRG). The accurate determination of the magnetization of this model is challenging for numerical methods and its value exhibits substantial disparities across various methods. We perform a large-scale DMRG calculation of this model by employing bond dimension as large asD=24000and by studying the system with width as large asLy=12. With careful extrapolation with truncation error and suitable finite size scaling, we give a conservative estimation of the magnetization asM0=0.208(8). The ground state energy per site we obtain isEg=-0.5503(8). Our results provide valuable benchmark values for the development of new methods in the future.
The behavior of the paradigmatic J 1 − J 2 triangular lattice Heisenberg antiferromagnet in a magnetic field remains unsettled despite decades of study. We map out the phase diagram using three complementary approaches, including self-consistent nonlinear spin-wave theory, density-matrix renormalization group, and variational Monte Carlo. This combined analysis resolves the competition among different field-induced magnetic orders and magnetization plateaus across the classically frustrated parameter range. In particular, there is a finite range in the parameter regime around J 2 / J 1 = 1 8 in which (i) upon the application of the external field, the gapless quantum spin liquid acquires a finite density of monopoles, and (ii) by further increasing the field, two plateaus are clearly obtained at m = 1 3 and m = 1 2 . We discuss the experimental importance of the consecutive magnetization plateau transitions as a signature of an underlying quantum spin-liquid phase.
We use the coupled cluster method (CCM) to study the zero-temperature ground-state (GS) properties of a spin-12\frac{1}{2}21​ ${J}_{1}\ensuremath{-}{J}_{2}$ Heisenberg antiferromagnet on a triangular lattice with competing nearest-neighbor and next-nearest-neighbor exchange couplings J1>0{J}_{1}>0J1​>0 and ${J}_{2}\ensuremath{\equiv}\ensuremath{\kappa}{J}_{1}>0$, respectively, in the window $0\ensuremath{\le}\ensuremath{\kappa}<1$. The classical version of the model has a single GS phase transition at ${\ensuremath{\kappa}}^{\mathrm{cl}}=\frac{1}{8}$ in this window from a phase with 3-sublattice antiferromagnetic (AFM) ${120}^{\ensuremath{\circ}}$ N\'eel order for $\ensuremath{\kappa}<{\ensuremath{\kappa}}^{\mathrm{cl}}$ to an infinitely degenerate family of 4-sublattice AFM N\'eel phases for $\ensuremath{\kappa}>{\ensuremath{\kappa}}^{\mathrm{cl}}$. This classical accidental degeneracy is lifted by quantum fluctuations, which favor a 2-sublattice AFM striped phase. For the quantum model we work directly in the thermodynamic limit of an infinite number of spins, with no consequent need for any finite-size scaling analysis of our results. We perform high-order CCM calculations within a well-controlled hierarchy of approximations, which we show how to extrapolate to the exact limit. In this way we find results for the case $\ensuremath{\kappa}=0$ of the spin-12\frac{1}{2}21​ model for the GS energy per spin, $E/N=\ensuremath{-}0.5521(2){J}_{1}$, and the GS magnetic order parameter, M=0.198(5)M=0.198(5)M=0.198(5) (in units where the classical value is Mcl=12{M}^{\mathrm{cl}}=\frac{1}{2}Mcl=21​), which are among the best available. For the spin-12\frac{1}{2}21​ ${J}_{1}\ensuremath{-}{J}_{2}$ model we find that the classical transition at $\ensuremath{\kappa}={\ensuremath{\kappa}}^{\mathrm{cl}}$ is split into two quantum phase transitions at ${\ensuremath{\kappa}}_{1}^{c}=0.060(10)$ and ${\ensuremath{\kappa}}_{2}^{c}=0.165(5)$. The two quasiclassical AFM states (viz., the ${120}^{\ensuremath{\circ}}$ N\'eel state and the striped state) are found to be the stable GS phases in the regime $\ensuremath{\kappa}<{\ensuremath{\kappa}}_{1}^{c}$ and $\ensuremath{\kappa}>{\ensuremath{\kappa}}_{2}^{c}$, respectively, while in the intermediate regimes ${\ensuremath{\kappa}}_{1}^{c}<\ensuremath{\kappa}<{\ensuremath{\kappa}}_{2}^{c}$ the stable GS phase has no evident long-range magnetic order.
We study the zero-temperature quantum phase diagram of a spin-1 Heisenberg antiferromagnet on the honeycomb lattice with both nearest-neighbor exchange coupling J1>0{J}_{1}>0J1​>0 and frustrating next-nearest-neighbor coupling ${J}_{2}\ensuremath{\equiv}\ensuremath{\kappa}{J}_{1}>0$, using the coupled cluster method implemented to high orders of approximation, and based on model states with different forms of classical magnetic order. For each we calculate directly in the bulk thermodynamic limit both ground-state low-energy parameters (including the energy per spin, magnetic order parameter, spin stiffness coefficient, and zero-field uniform transverse magnetic susceptibility) and their generalized susceptibilities to various forms of valence-bond crystalline (VBC) order, as well as the energy gap to the lowest-lying spin-triplet excitation. In the range $0<\ensuremath{\kappa}<1$ we find evidence for four distinct phases. Two of these are quasiclassical phases with antiferromagnetic long-range order, one with two-sublattice N\'eel order for $\ensuremath{\kappa}<{\ensuremath{\kappa}}_{{c}_{1}}=0.250(5)$, and another with four-sublattice N\'eel-II order for $\ensuremath{\kappa}>{\ensuremath{\kappa}}_{{c}_{2}}=0.340(5)$. Two different paramagnetic phases are found to exist in the intermediate region. Over the range ${\ensuremath{\kappa}}_{{c}_{1}}<\ensuremath{\kappa}<{\ensuremath{\kappa}}_{c}^{i}=0.305(5)$ we find a gapless phase with no discernible magnetic order, which is a strong candidate for being a quantum spin liquid, while over the range ${\ensuremath{\kappa}}_{c}^{i}<\ensuremath{\kappa}<{\ensuremath{\kappa}}_{{c}_{2}}$ we find a gapped phase, which is most likely a lattice nematic with staggered dimer VBC order that breaks the lattice rotational symmetry.
We study the S=1/2S=1/2S=1/2 Heisenberg model on the triangular lattice with nearest- and next-nearest-neighbor interactions J1{J}_{1}J1​ and J2{J}_{2}J2​ with the density matrix renormalization group, on long open cylinders with widths up to nine lattice spacings. In an intermediate J2{J}_{2}J2​ region $0.06\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\lesssim}0.17$, we find evidence for a spin liquid (SL) state with short range spin-spin, bond-bond, and chiral correlation lengths, bordered by a classical ${120}^{\ensuremath{\circ}}$ N\'eel ordered state at small J2{J}_{2}J2​ and by a two sublattice collinear magnetically ordered state at larger J2{J}_{2}J2​. Focusing on J2/J1=0.1{J}_{2}/{J}_{1}=0.1J2​/J1​=0.1, we find a number of signatures of a gapped SL phase.
We investigate the ground-state magnetic structure of the spin- 1 / 2 XXZ antiferromagnet on the triangular lattice in the easy-axis regime using the density-matrix renormalization group. By applying spiral boundary conditions, we exactly map finite L × L clusters onto one-dimensional chains while avoiding the spatial anisotropy inherent in cylindrical geometries. From symmetry-broken local magnetization profiles, we extract the three-sublattice moments and track their evolution with anisotropy. At the isotropic point, we obtain a positive sublattice moment of 0.217(3), consistent with previous numerical estimates. In the easy-axis regime ( Δ = J z / J ⊥ > 1 ), the ordered moments remain close to a Y-like zero-magnetization three-sublattice state, whose z -component pattern is of the form ( 2 m , − m , − m ) , over a broad range of Δ . Extrapolation in 1 / Δ shows that the positive sublattice moment stays well below the classical saturation value 1 / 2 , approaching 0.419(7) as Δ → ∞ , while the magnitude of the negative sublattice moment approaches 0.209(4). We further compare the energies of the Y state and the up-down-down state and find that the Y state is favored at zero field. Independent thermodynamic-limit energy calculations, performed without assuming any particular ordered pattern, yield an energy consistent with the Y-state solution. These results show that the easy-axis ground state does not simply cross over to a trivially saturated collinear Ising state, but instead remains a nontrivial three-sublattice ordered state selected from the macroscopically degenerate Ising manifold by quantum fluctuations.