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

Physical review. B, Condensed matter

Spin and dimer correlations forS=1 bilinear-biquadratic antiferromagnetic Heisenberg chains

J. J. Deisz

Published Aug 1, 1992• 4 citations
Publisher

Abstract

Using the world-line quantum Monte Carlo method, we calculate the spin 〈${\mathit{S}}_{\mathit{R}}^{\mathit{z}}$${\mathit{S}}_{0}^{\mathit{z}}$〉 and dimer 〈(${\mathbf{S}}_{0}$\ensuremath{\cdot}${\mathbf{S}}_{1}$)(${\mathbf{S}}_{\mathit{R}}$\ensuremath{\cdot}${\mathbf{S}}_{\mathit{R}+1}$)〉 correlation functions for the S=1 biquadratic Heisenberg antiferromagnetic chain. These results show the relevance of the Affleck-Kennedy-Lieb-Tasaki state for describing the nonclassical behavior of the S=1 antiferromagnetic Heisenberg model. We also include results for S=3/2 and 2, which show the inhibition of staggered dimerization as the classical limit (S\ensuremath{\rightarrow}\ensuremath{\infty}) is approached, where N\'eel order is favored.

DOI
10.1103/physrevb.46.2885

Related papers

Haldane phases and phase diagrams of the S=32 and S=1 bilinear-biquadratic Heisenberg model on the orthogonal dimer chain

Ke Ren, Muwei Wu, Shou-Shu Gong +2 more

We systematically study the effects of higher-order interactions on the $S=\frac{3}{2},1$ orthogonal dimer chains using exact diagonalization and density matrix renormalization group. Due to frustration and higher spin, there are rich quantum phases, including three Haldane phases, two gapless phases, and several magnetically ordered phases. To characterize these phases and their phase transitions, we study various physical quantities such as energy gap, energy-level crossing, fidelity susceptibility, spin correlation, entanglement spectrum, and central charge. According to our calculations, the biquadratic term can enhance the Haldane phase regions. In particular, we numerically identify that a Haldane phase in $S=\frac{3}{2}$ case is adiabatically connected to the exact AKLT (Affleck, Kennedy, Lieb, and Tasaki) point when adding a bicubic term. Our study on the orthogonal dimer model, which is a 1D version of the Shastry-Sutherland model, provides insights into understanding the possible $S=\frac{3}{2},1$ Haldane phases in quasi-1D and 2D frustrated magnetic materials.

arXiv:2307.06371•Physical review. B./Physical review. B•Dec 1, 2023•3 citations

Monte Carlo study of the spin-1 antiferromagnetic chain with both bilinear and biquadratic exchange interactions

M. Uchinami

A Monte Carlo simulation study based on the heat-bath algorithm on a checkerboard arrangement is applied to the spin-1 antiferromagnetic chain model with both bilinear and biquadratic exchange interactions. Changing the parameter \ensuremath{\lambda} which denotes strength of the biquadratic interaction against the bilinear one from -1.0 to 2.0, we obtain the following results according to the measurements of energy and two-point correlation function. At \ensuremath{\lambda}=1 and -(1/3 where some physical quantities are derived rigorously, our Monte Carlo simulation results agree well with the rigorous ones. For ${\ensuremath{\lambda}}_{c}$\ensuremath{\le}\ensuremath{\lambda}\ensuremath{\le}${\ensuremath{\lambda}}_{1}$ with taking 0.25\ensuremath{\le}${\ensuremath{\lambda}}_{c}$\ensuremath{\le}0.5 and ${\ensuremath{\lambda}}_{1}$\ensuremath{\ge}2.0, its model has the gapless ground state accompanied by a power-law decaying correlation function. This implies the existence of a critical line. For -1.0${\ensuremath{\lambda}}_{c}$, on the other hand, its model has the ground state with energy gap accompanied with exponentially decaying correlation function. This indicates that the valence-bond-solid (VBS) state asserted by Affleck et al. is stabilized there. As the point of \ensuremath{\lambda}=0 comes to the standard Heisenberg antiferromagnetic model and is in the region of -1.0${\ensuremath{\lambda}}_{c}$, the Haldane's conjecture is supported that this ground state is described by the VBS state for the case of integer spin, e.g., S=1.

Physical review. B, Condensed matter•Mar 1, 1989•7 citations

Gutzwiller projected wave functions in the fermionic theory of S=1 spin chains

Zheng-Xin Liu, Yi Zhou, Hong-Hao Tu +2 more

We study in this paper a series of Gutzwiller projected wave functions for $S=1$ spin chains obtained from a fermionic mean-field theory for general $S>1/2$ spin systems [Liu, Zhou, and Ng, Phys. Rev. B 81, 224417 (2010)] applied to the bilinear-biquadratic ($J$-$K$) model. The free-fermion mean-field states before the projection are 1D paring states. By comparing the energies and correlation functions of the projected pairing states with those obtained from known results, we show that the optimized Gutzwiller projected wave functions are very good trial ground-state wave functions for the antiferromagnetic bilinear-biquadratic model in the regime $K<J$ ($\ensuremath{-}3\ensuremath{\pi}/4<\ensuremath{\theta}<\ensuremath{\pi}/4$). We find that different topological phases of the free-fermion paring states correspond to different spin phases: the weak pairing (topologically nontrivial) state gives rise to the Haldane phase, whereas the strong pairing (topologically trivial) state gives rise to the dimer phase. In particular, the mapping between the Haldane phase and Gutwziller wave function is exact at the Affleck-Kennedy-Lieb-Tasaki (AKLT) point $K/J=1/3$ [$\ensuremath{\theta}={\mathrm{tan}}^{\ensuremath{-}1}(1/3)$]. The transition point between the two phases determined by the optimized Gutzwiller projected wave function is in good agreement with the known result. The effect of ${Z}_{2}$ gauge fluctuations above the mean-field theory is analyzed.

arXiv:1202.5723•Physical Review B•May 29, 2012•26 citations

Variational and density-matrix renormalization-group studies of the frustrated antiferromagnetic Heisenberg S=1 quantum spin chain

Alexei K. Kolezhuk, R. Roth, Ulrich Schollwöck

In this paper we study a frustrated antiferromagnetic isotropic Heisenberg S=1 quantum spin chain H=${\ensuremath{\sum}}_{\mathrm{i}}$${\mathbf{S}}_{\mathrm{i}}$${\mathbf{S}}_{\mathrm{i}+1}$+\ensuremath{\alpha}${\ensuremath{\sum}}_{\mathrm{i}}$${\mathbf{S}}_{\mathrm{i}}$${\mathbf{S}}_{\mathrm{i}+2}$, using a variational ansatz starting from valence bond states and the density-matrix renormalization group. We find both methods to give results in very good qualitative and good quantitative agreement, which clarify the phase diagram as follows: At ${\mathrm{\ensuremath{\alpha}}}_{\mathrm{D}}$=0.284(1), there is a disorder point of the first kind, marking the onset of incommensurate spin-spin correlations in the chain. At ${\mathrm{\ensuremath{\alpha}}}_{\mathrm{L}}$=0.3725(25) there is a Lifshitz point, at which the excitation spectrum is found to develop a particular doubly degenerate structure. These points are the quantum remnants of the transition from antiferromagnetic to spiral order in the classical frustrated chain. At ${\mathrm{\ensuremath{\alpha}}}_{\mathrm{T}}$=0.7444(6) there is a first-order phase transition from an Affleck-Kennedy-Lieb-Tasaki (AKLT) phase characterized by nonvanishing string order to a phase which can be understood as a next-nearest-neighbor generalization of the AKLT model. At the transition, the string order parameter shows (to numerical precision) a discontinuous jump of 0.085 to zero; the correlation length and the gap are both finite at the transition. The problem of edge states in open frustrated chains is discussed at length.

arXiv:cond-mat/9609064•Physical review. B, Condensed matter•Apr 1, 1997•51 citations

Dimer-Quadrupolar Quantum Phase Transition in the Quasi-One-Dimensional Heisenberg Model with Biquadratic Interaction

Kenji Harada, Naoki Kawashima, Matthias Troyer

The quasi-one-dimensional S=1 Heisenberg antiferromagnet with a biquadratic term is investigated at zero temperature by quantum Monte Carlo simulation. As the magnitude of the inter-chain coupling is increased, the system undergoes a phase transition from a spontaneously dimerized phase to a Néel ordered or spin nematic phase. The numerical results suggest the possibility of an unconventional second-order transition in which the symmetry group characterizing one phase is not a subgroup of the other.

arXiv:cond-mat/0608446•Journal of the Physical Society of Japan•Dec 25, 2006•40 citations

Exact asymptotic correlation functions of bilinear spin operators of the Heisenberg antiferromagnetic spin- 12 chain

T. Vekua, Gaoyong Sun

Exact asymptotic expressions of the uniform parts of the two-point correlation functions of bilinear spin operators in the Heisenberg antiferromagnetic spin-$\frac{1}{2}$ chain are obtained. Apart from the algebraic decay, the logarithmic contribution is identified, and the numerical prefactor is determined. We also confirm numerically the multiplicative logarithmic correction of the staggered part of the bilinear spin operators $\ensuremath{\langle}\ensuremath{\langle}{S}_{0}^{a}{S}_{1}^{a}{S}_{r}^{b}{S}_{r+1}^{b}\ensuremath{\rangle}\ensuremath{\rangle}={(\ensuremath{-}1)}^{r}d/(r{ln}^{\frac{3}{2}}r)+(3{\ensuremath{\delta}}_{a,b}\ensuremath{-}1){ln}^{2}r/(12{\ensuremath{\pi}}^{4}{r}^{4})$, and estimate the numerical prefactor as $d\ensuremath{\simeq}0.067$. The relevance of our results for ground-state fidelity susceptibility at the Berezinskii-Kosterlitz-Thouless quantum phase transition points in one-dimensional systems is discussed at the end of our work.

arXiv:1604.00939•Physical review. B./Physical review. B•Jul 14, 2016•7 citations
Using the world-line quantum Monte Carlo method, we calculate the spin 〈{\mathit{S}}_{\mathit{R}}^{\mathit{z}}$${\mathit{S}}_{0}^{\mathit{z}}〉 and dimer 〈(S0{\mathbf{S}}_{0}S0​\ensuremath{\cdot}S1{\mathbf{S}}_{1}S1​)(SR{\mathbf{S}}_{\mathit{R}}SR​\ensuremath{\cdot}SR+1{\mathbf{S}}_{\mathit{R}+1}SR+1​)〉 correlation functions for the S=1 biquadratic Heisenberg antiferromagnetic chain. These results show the relevance of the Affleck-Kennedy-Lieb-Tasaki state for describing the nonclassical behavior of the S=1 antiferromagnetic Heisenberg model. We also include results for S=3/2 and 2, which show the inhibition of staggered dimerization as the classical limit (S\ensuremath{\rightarrow}\ensuremath{\infty}) is approached, where N'eel order is favored.
We systematically study the effects of higher-order interactions on the S=32,1S=\frac{3}{2},1S=23​,1 orthogonal dimer chains using exact diagonalization and density matrix renormalization group. Due to frustration and higher spin, there are rich quantum phases, including three Haldane phases, two gapless phases, and several magnetically ordered phases. To characterize these phases and their phase transitions, we study various physical quantities such as energy gap, energy-level crossing, fidelity susceptibility, spin correlation, entanglement spectrum, and central charge. According to our calculations, the biquadratic term can enhance the Haldane phase regions. In particular, we numerically identify that a Haldane phase in S=32S=\frac{3}{2}S=23​ case is adiabatically connected to the exact AKLT (Affleck, Kennedy, Lieb, and Tasaki) point when adding a bicubic term. Our study on the orthogonal dimer model, which is a 1D version of the Shastry-Sutherland model, provides insights into understanding the possible S=32,1S=\frac{3}{2},1S=23​,1 Haldane phases in quasi-1D and 2D frustrated magnetic materials.
A Monte Carlo simulation study based on the heat-bath algorithm on a checkerboard arrangement is applied to the spin-1 antiferromagnetic chain model with both bilinear and biquadratic exchange interactions. Changing the parameter \ensuremath{\lambda} which denotes strength of the biquadratic interaction against the bilinear one from -1.0 to 2.0, we obtain the following results according to the measurements of energy and two-point correlation function. At \ensuremath{\lambda}=1 and -(1/3 where some physical quantities are derived rigorously, our Monte Carlo simulation results agree well with the rigorous ones. For ${\ensuremath{\lambda}}_{c}$\ensuremath{\le}\ensuremath{\lambda}\ensuremath{\le}${\ensuremath{\lambda}}_{1}$ with taking 0.25\ensuremath{\le}${\ensuremath{\lambda}}_{c}$\ensuremath{\le}0.5 and ${\ensuremath{\lambda}}_{1}$\ensuremath{\ge}2.0, its model has the gapless ground state accompanied by a power-law decaying correlation function. This implies the existence of a critical line. For -1.0${\ensuremath{\lambda}}_{c}$, on the other hand, its model has the ground state with energy gap accompanied with exponentially decaying correlation function. This indicates that the valence-bond-solid (VBS) state asserted by Affleck et al. is stabilized there. As the point of \ensuremath{\lambda}=0 comes to the standard Heisenberg antiferromagnetic model and is in the region of -1.0${\ensuremath{\lambda}}_{c}$, the Haldane's conjecture is supported that this ground state is described by the VBS state for the case of integer spin, e.g., S=1.
We study in this paper a series of Gutzwiller projected wave functions for S=1S=1S=1 spin chains obtained from a fermionic mean-field theory for general S>1/2S>1/2S>1/2 spin systems [Liu, Zhou, and Ng, Phys. Rev. B 81, 224417 (2010)] applied to the bilinear-biquadratic (JJJ-KKK) model. The free-fermion mean-field states before the projection are 1D paring states. By comparing the energies and correlation functions of the projected pairing states with those obtained from known results, we show that the optimized Gutzwiller projected wave functions are very good trial ground-state wave functions for the antiferromagnetic bilinear-biquadratic model in the regime K<JK<JK<J ($\ensuremath{-}3\ensuremath{\pi}/4<\ensuremath{\theta}<\ensuremath{\pi}/4$). We find that different topological phases of the free-fermion paring states correspond to different spin phases: the weak pairing (topologically nontrivial) state gives rise to the Haldane phase, whereas the strong pairing (topologically trivial) state gives rise to the dimer phase. In particular, the mapping between the Haldane phase and Gutwziller wave function is exact at the Affleck-Kennedy-Lieb-Tasaki (AKLT) point K/J=1/3K/J=1/3K/J=1/3 [$\ensuremath{\theta}={\mathrm{tan}}^{\ensuremath{-}1}(1/3)$]. The transition point between the two phases determined by the optimized Gutzwiller projected wave function is in good agreement with the known result. The effect of Z2{Z}_{2}Z2​ gauge fluctuations above the mean-field theory is analyzed.
In this paper we study a frustrated antiferromagnetic isotropic Heisenberg S=1 quantum spin chain H=${\ensuremath{\sum}}_{\mathrm{i}}$Si{\mathbf{S}}_{\mathrm{i}}Si​Si+1{\mathbf{S}}_{\mathrm{i}+1}Si+1​+\ensuremath{\alpha}${\ensuremath{\sum}}_{\mathrm{i}}$Si{\mathbf{S}}_{\mathrm{i}}Si​Si+2{\mathbf{S}}_{\mathrm{i}+2}Si+2​, using a variational ansatz starting from valence bond states and the density-matrix renormalization group. We find both methods to give results in very good qualitative and good quantitative agreement, which clarify the phase diagram as follows: At ${\mathrm{\ensuremath{\alpha}}}_{\mathrm{D}}$=0.284(1), there is a disorder point of the first kind, marking the onset of incommensurate spin-spin correlations in the chain. At ${\mathrm{\ensuremath{\alpha}}}_{\mathrm{L}}$=0.3725(25) there is a Lifshitz point, at which the excitation spectrum is found to develop a particular doubly degenerate structure. These points are the quantum remnants of the transition from antiferromagnetic to spiral order in the classical frustrated chain. At ${\mathrm{\ensuremath{\alpha}}}_{\mathrm{T}}$=0.7444(6) there is a first-order phase transition from an Affleck-Kennedy-Lieb-Tasaki (AKLT) phase characterized by nonvanishing string order to a phase which can be understood as a next-nearest-neighbor generalization of the AKLT model. At the transition, the string order parameter shows (to numerical precision) a discontinuous jump of 0.085 to zero; the correlation length and the gap are both finite at the transition. The problem of edge states in open frustrated chains is discussed at length.
The quasi-one-dimensional S=1 Heisenberg antiferromagnet with a biquadratic term is investigated at zero temperature by quantum Monte Carlo simulation. As the magnitude of the inter-chain coupling is increased, the system undergoes a phase transition from a spontaneously dimerized phase to a Néel ordered or spin nematic phase. The numerical results suggest the possibility of an unconventional second-order transition in which the symmetry group characterizing one phase is not a subgroup of the other.
Exact asymptotic expressions of the uniform parts of the two-point correlation functions of bilinear spin operators in the Heisenberg antiferromagnetic spin-12\frac{1}{2}21​ chain are obtained. Apart from the algebraic decay, the logarithmic contribution is identified, and the numerical prefactor is determined. We also confirm numerically the multiplicative logarithmic correction of the staggered part of the bilinear spin operators $\ensuremath{\langle}\ensuremath{\langle}{S}_{0}^{a}{S}_{1}^{a}{S}_{r}^{b}{S}_{r+1}^{b}\ensuremath{\rangle}\ensuremath{\rangle}={(\ensuremath{-}1)}^{r}d/(r{ln}^{\frac{3}{2}}r)+(3{\ensuremath{\delta}}_{a,b}\ensuremath{-}1){ln}^{2}r/(12{\ensuremath{\pi}}^{4}{r}^{4})$, and estimate the numerical prefactor as $d\ensuremath{\simeq}0.067$. The relevance of our results for ground-state fidelity susceptibility at the Berezinskii-Kosterlitz-Thouless quantum phase transition points in one-dimensional systems is discussed at the end of our work.