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:1411.0789 · Physical Review B

Absence of classical long-range order in an S=12 Heisenberg antiferromagnet on a triangular lattice

Nobuo Suzuki, Fumitaka Matsubara, Sumiyoshi Fujiki, Takayuki Shirakura

Published Nov 14, 2014• 12 citations
arXivHTMLPDF

Abstract

We study the quantum phase transition of an $S=\frac{1}{2}$ anisotropic $\ensuremath{\alpha}$ $(\ensuremath{\equiv}{J}_{z}/{J}_{xy})$ Heisenberg antiferromagnet on a triangular lattice. We calculate the sublattice magnetization and the long-range helical order parameter and their Binder ratios on finite systems with $N\ensuremath{\le}36$ sites. The $N$ dependence of the Binder ratios reveals that the classical ${120}^{\ensuremath{\circ}}$ N\'eel state occurs for $\ensuremath{\alpha}\ensuremath{\lesssim}0.55$, whereas a critical collinear state occurs for $1/\ensuremath{\alpha}\ensuremath{\lesssim}0.6$. This result is at odds with a widely held belief that the ground state of a Heisenberg antiferromagnet is the ${120}^{\ensuremath{\circ}}$ N\'eel state, but it also provides a possible mechanism explaining experimentally observed spin liquids.

DOI
10.1103/physrevb.90.184414

Related papers

Phase transition of theS=1/2 quantum anisotropic Heisenberg antiferromagnet on the triangular lattice

Naoki Suzuki, Fumitaka Matsubara

We study the S=1/2 anisotropic Heisenberg antiferromagnet on finite triangular lattices with N\ensuremath{\le}24 sites: H=2J${\mathcal{J}}_{\mathrm{〈}\mathit{i},}$${\mathit{j}}_{\mathrm{〉}}$(${\mathit{S}}_{\mathit{i}}^{\mathit{x}}$${\mathit{S}}_{\mathit{j}}^{\mathit{x}}$ +${\mathit{S}}_{\mathit{i}}^{\mathit{y}}$${\mathit{S}}_{\mathit{j}}^{\mathit{y}}$ +\ensuremath{\Delta}${\mathit{S}}_{\mathit{i}}^{\mathit{z}}$${\mathit{S}}_{\mathit{j}}^{\mathit{z}}$) with 0\ensuremath{\le}\ensuremath{\Delta}\ensuremath{\le}1. The specific heat C and the chiral-order parameter 〈${\mathrm{\ensuremath{\chi}}}^{2}$〉 are calculated using a quantum transfer Monte Carlo method. Remarkable differences in the size dependences of C and 〈${\mathrm{\ensuremath{\chi}}}^{2}$〉 are found between the cases of \ensuremath{\Delta}\ensuremath{\le}0.4 and of \ensuremath{\Delta}\ensuremath{\ge}0.6. For \ensuremath{\Delta}\ensuremath{\le}0.4, the peak height of C increases with increasing N and an extrapolation of 〈${\mathrm{\ensuremath{\chi}}}^{2}$〉 to the thermodynamic limit gives a finite, nonzero value at low temperatures. In contrast with these, for \ensuremath{\Delta}\ensuremath{\ge}0.6, the peak height does not increase and the extrapolation of 〈${\mathrm{\ensuremath{\chi}}}^{2}$〉 gives a smaller value even at very low temperatures indicating the absence of a long-range chiral order. From the results, we suggest that the chiral-ordered phase transition occurs at a finite, nonzero temperature when \ensuremath{\Delta}\ensuremath{\le}${\mathrm{\ensuremath{\Delta}}}_{\mathit{c}}$ with ${\mathrm{\ensuremath{\Delta}}}_{\mathit{c}}$\ensuremath{\gtrsim}0.4. The phase diagram of the model is predicted.

Physical review. B, Condensed matter•Mar 1, 1995•8 citations

MAGNETIC ORDERING OF ANTIFERROMAGNETS ON A SPATIALLY ANISOTROPIC TRIANGULAR LATTICE

R. F. Bishop, P. H. Y. Li, DAMIAN J. J. FARNELL +1 more

We study the spin-1/2 and spin-1 [Formula: see text] Heisenberg antiferromagnets (HAFs) on an infinite, anisotropic, two-dimensional triangular lattice, using the coupled cluster method. With respect to an underlying square-lattice geometry the model contains antiferromagnetic (J1 > 0) bonds between nearest neighbours and competing [Formula: see text] bonds between next-nearest-neighbours across only one of the diagonals of each square plaquette, the same diagonal in each square. In a topologically equivalent triangular-lattice geometry the model has two sorts of nearest-neighbour bonds, with [Formula: see text] bonds along parallel chains and with J1 bonds providing an interchain coupling. The model thus interpolates between an isotropic HAF on the square lattice at one extreme (κ = 0) and a set of decoupled chains at the other (κ → ∞), with the isotropic HAF on the triangular lattice in between at κ = 1. For the spin-1/2 [Formula: see text] model, we find a weakly first-order (or possibly second-order) quantum phase transition from a Néel-ordered state to a helical state at a first critical point at κc1 = 0.80 ± 0.01, and a second critical point at κc2 = 1.8 ± 0.4 where a first-order transition occurs between the helical state and a collinear stripe-ordered state. For the corresponding spin-1 model we find an analogous transition of the second-order type at κc1 = 0.62 ± 0.01 between states with Néel and helical ordering, but we find no evidence of a further transition in this case to a stripe-ordered phase.

International Journal of Modern Physics B•Oct 20, 2010•4 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

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

Ground state properties of the S=1/2 Heisenberg antiferromagnet on the triangular lattice

Luca Capriotti, A. E. Trumper, Sandro Sorella

We investigate the ground state properties of the S=1/2 Heisenberg antiferromagnet on the triangular lattice using mainly exact diagonalization (for lattices up to N=36 sites) and the recently developed Green Function Monte Carlo with stochastic reconfiguration (up to N=144). A systematic size scaling of the order parameter and of the spin gap gives strong evidences for the presence of long-range Néel order in the ground state of this highly frustrated model.

Journal of Applied Physics•May 1, 2000•2 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
We study the quantum phase transition of an S=12S=\frac{1}{2}S=21​ anisotropic \ensuremathα\ensuremath{\alpha}\ensuremathα (\ensuremath≡Jz/Jxy)(\ensuremath{\equiv}{J}_{z}/{J}_{xy})(\ensuremath≡Jz​/Jxy​) Heisenberg antiferromagnet on a triangular lattice. We calculate the sublattice magnetization and the long-range helical order parameter and their Binder ratios on finite systems with N\ensuremath≤36N\ensuremath{\le}36N\ensuremath≤36 sites. The NNN dependence of the Binder ratios reveals that the classical 120\ensuremath∘{120}^{\ensuremath{\circ}}120\ensuremath∘ N'eel state occurs for \ensuremathα\ensuremath≲0.55\ensuremath{\alpha}\ensuremath{\lesssim}0.55\ensuremathα\ensuremath≲0.55, whereas a critical collinear state occurs for 1/\ensuremathα\ensuremath≲0.61/\ensuremath{\alpha}\ensuremath{\lesssim}0.61/\ensuremathα\ensuremath≲0.6. This result is at odds with a widely held belief that the ground state of a Heisenberg antiferromagnet is the 120\ensuremath∘{120}^{\ensuremath{\circ}}120\ensuremath∘ N'eel state, but it also provides a possible mechanism explaining experimentally observed spin liquids.
We study the S=1/2 anisotropic Heisenberg antiferromagnet on finite triangular lattices with N\ensuremath{\le}24 sites: H=2JJ〈i,{\mathcal{J}}_{\mathrm{〈}\mathit{i},}J〈i,​j〉{\mathit{j}}_{\mathrm{〉}}j〉​(Six{\mathit{S}}_{\mathit{i}}^{\mathit{x}}Six​Sjx{\mathit{S}}_{\mathit{j}}^{\mathit{x}}Sjx​ +Siy{\mathit{S}}_{\mathit{i}}^{\mathit{y}}Siy​Sjy{\mathit{S}}_{\mathit{j}}^{\mathit{y}}Sjy​ +\ensuremath{\Delta}Siz{\mathit{S}}_{\mathit{i}}^{\mathit{z}}Siz​Sjz{\mathit{S}}_{\mathit{j}}^{\mathit{z}}Sjz​) with 0\ensuremath{\le}\ensuremath{\Delta}\ensuremath{\le}1. The specific heat C and the chiral-order parameter 〈${\mathrm{\ensuremath{\chi}}}^{2}$〉 are calculated using a quantum transfer Monte Carlo method. Remarkable differences in the size dependences of C and 〈${\mathrm{\ensuremath{\chi}}}^{2}$〉 are found between the cases of \ensuremath{\Delta}\ensuremath{\le}0.4 and of \ensuremath{\Delta}\ensuremath{\ge}0.6. For \ensuremath{\Delta}\ensuremath{\le}0.4, the peak height of C increases with increasing N and an extrapolation of 〈${\mathrm{\ensuremath{\chi}}}^{2}$〉 to the thermodynamic limit gives a finite, nonzero value at low temperatures. In contrast with these, for \ensuremath{\Delta}\ensuremath{\ge}0.6, the peak height does not increase and the extrapolation of 〈${\mathrm{\ensuremath{\chi}}}^{2}$〉 gives a smaller value even at very low temperatures indicating the absence of a long-range chiral order. From the results, we suggest that the chiral-ordered phase transition occurs at a finite, nonzero temperature when \ensuremath{\Delta}\ensuremath{\le}${\mathrm{\ensuremath{\Delta}}}_{\mathit{c}}$ with ${\mathrm{\ensuremath{\Delta}}}_{\mathit{c}}$\ensuremath{\gtrsim}0.4. The phase diagram of the model is predicted.
We study the spin-1/2 and spin-1 [Formula: see text] Heisenberg antiferromagnets (HAFs) on an infinite, anisotropic, two-dimensional triangular lattice, using the coupled cluster method. With respect to an underlying square-lattice geometry the model contains antiferromagnetic (J1 > 0) bonds between nearest neighbours and competing [Formula: see text] bonds between next-nearest-neighbours across only one of the diagonals of each square plaquette, the same diagonal in each square. In a topologically equivalent triangular-lattice geometry the model has two sorts of nearest-neighbour bonds, with [Formula: see text] bonds along parallel chains and with J1 bonds providing an interchain coupling. The model thus interpolates between an isotropic HAF on the square lattice at one extreme (κ = 0) and a set of decoupled chains at the other (κ → ∞), with the isotropic HAF on the triangular lattice in between at κ = 1. For the spin-1/2 [Formula: see text] model, we find a weakly first-order (or possibly second-order) quantum phase transition from a Néel-ordered state to a helical state at a first critical point at κc1 = 0.80 ± 0.01, and a second critical point at κc2 = 1.8 ± 0.4 where a first-order transition occurs between the helical state and a collinear stripe-ordered state. For the corresponding spin-1 model we find an analogous transition of the second-order type at κc1 = 0.62 ± 0.01 between states with Néel and helical ordering, but we find no evidence of a further transition in this case to a stripe-ordered phase.
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 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 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.