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

1 The spin-1/2 J1-J2 Heisenberg antiferromagnet on the square lattice: Exact diagonalization for N = 40 spins

Johannes Richter, J. Schulenburg

Published Jan 1, 2009• 91 citations
OpenAlex

Abstract

We present numerical exact results for the ground state and the low-lying excitations for the spin-1/2 J1-J2 Heisenberg antiferromagnet on finite square lattices of up to N = 40 sites. Using finite-size extrapolation we determine the ground-state energy, the magnetic order parameters, the spin gap, the uniform susceptibility, as well as the spin-wave velocity and the spin stiffness as functions of the frustration parameter J2/J1. In agreement with the generally excepted scenario we find semiclassical magnetically ordered phases for J2 < J c1 2 and J2> J c2 2 separated by a gapful quantum paramagnetic phase. We estimate Jc1 2 ≈ 0.35J1 and J c2 2 ≈ 0.66J1. PACS codes: 75.10.Jm Quantized spin models 75.45.+j Macroscopic quantum phenomena in magnetic systems 1

Related papers

Ground state phases of the spin-1/2 J1–J2 Heisenberg antiferromagnet on the square lattice: A high-order coupled cluster treatment

R. Darradi, Oleg Derzhko, Ronald Zinke +3 more

Using the coupled cluster method for high orders of approximation and complementary exact diagonalization studies we investigate the ground state properties of the spin-1/2 ${J}_{1}--{J}_{2}$ frustrated Heisenberg antiferromagnet on the square lattice. We have calculated the ground state energy, the magnetic order parameter, the spin stiffness, and several generalized susceptibilities to probe magnetically disordered quantum valence-bond phases. We have found that the quantum critical points for both the N\'eel and collinear orders are ${J}_{2}^{c1}\ensuremath{\approx}(0.44\ifmmode\pm\else\textpm\fi{}0.01){J}_{1}$ and ${J}_{2}^{c2}\ensuremath{\approx}(0.59\ifmmode\pm\else\textpm\fi{}0.01){J}_{1}$, respectively, which are in good agreement with the results obtained by other approximations. In contrast to the recent study by [Sirker et al. Phys. Rev. B 73, 184420 (2006)], our data do not provide evidence for the transition from the N\'eel to the valence-bond solid state to be first order. Moreover, our results are in favor of the deconfinement scenario for that phase transition. We also discuss the nature of the magnetically disordered quantum phase.

arXiv:0806.3825•Physical Review B•Dec 11, 2008•183 citations

Direct calculation of the spin stiffness in the J1 - J2 Heisenberg antiferromagnet

Torbjörn Einarsson, Heinz J. Schulz

We calculate the spin stiffness ${\mathrm{\ensuremath{\rho}}}_{\mathit{s}}$ for the frustrated spin-1/2 Heisenberg antiferromagnet on a square lattice by exact diagonalizations on finite clusters of up to 36 sites followed by extrapolations to the thermodynamic limit. For the nonfrustrated case, we find that ${\mathrm{\ensuremath{\rho}}}_{\mathit{s}}$=(0.183\ifmmode\pm\else\textpm\fi{}0.003)${\mathit{J}}_{1}$, in excellent agreement with the best results obtained by other means. Turning on frustration, the extrapolated stiffness vanishes for 0.4\ensuremath{\lesssim}${\mathit{J}}_{2}$/${\mathit{J}}_{1}$\ensuremath{\lesssim}0.6. In this intermediate region, the finite-size scaling works poorly---an additional sign that there is neither N\'eel nor collinear magnetic order. Using a hydrodynamic relation, and previous results for the transverse susceptibility, we also estimate the spin-wave velocity in the N\'eel-ordered region.

arXiv:cond-mat/9410090•Physical review. B, Condensed matter•Mar 1, 1995•81 citations

Tensor-product state approach to spin- 12 square J1−J2 antiferromagnetic Heisenberg model: Evidence for deconfined quantum criticality

Ling Wang, Zheng‐Cheng Gu, Frank Verstraete +1 more

The ground state phase of a spin-$\frac{1}{2} {J}_{1}\text{\ensuremath{-}}{J}_{2}$ antiferromagnetic Heisenberg model on a square lattice around the maximally frustrated regime $({J}_{2}\ensuremath{\sim}0.5{J}_{1})$ has been debated for decades. Here we study this model using the cluster update algorithm for tensor-product states (TPSs). The ground state energies at finite sizes and in the thermodynamic limit (with finite size scaling) are in good agreement with exact diagonalization study. Through finite size scaling of the spin correlation function, we find the critical point ${J}_{2}^{{c}_{1}}=0.572(5){J}_{1}$ and critical exponents $\ensuremath{\nu}=0.50(8),\phantom{\rule{0.16em}{0ex}}{\ensuremath{\eta}}_{s}=0.28(6)$. In the range of $0.572<{J}_{2}/{J}_{1}\ensuremath{\le}0.6$ we find a paramagnetic ground state with an exponentially decaying spin-spin correlation. Up to a $24\ifmmode\times\else\texttimes\fi{}24$ system size, we observe power law decaying dimer-dimer and plaquette-plaquette correlations with an anomalous plaquette scaling exponent ${\ensuremath{\eta}}_{p}=0.24(1)$ and an anomalous columnar scaling exponent ${\ensuremath{\eta}}_{c}=0.28(1)$ at ${J}_{2}/{J}_{1}=0.6$. These results are consistent with a potential gapless $U(1)$ spin-liquid phase. However, since the $U(1)$ spin liquid is unstable due to the instanton effect, a valence bond solid order with very small amplitude might develop in the thermodynamic limit. Thus, our numerical results strongly indicate a deconfined quantum critical point at ${J}_{2}^{{c}_{1}}$. Remarkably, all the observed critical exponents are consistent with the $J\text{\ensuremath{-}}Q$ model.

arXiv:1112.3331•Physical review. B./Physical review. B•Aug 19, 2016•102 citations

Spin-liquid Phase in Spin-1/2 square J_1-J_2 Heisenberg Model: A Tensor Product State Approach

Ling Wang

The ground state phase of spin-1/2 J_1-J_2 antiferromagnetic Heisenberg model on square lattice in the maximally frustrated regime (J_2 ~ 0.5J_1) has been debated for decades. Here we study this model by using a recently proposed novel numerical method - the cluster update algorithm for tensor product states (TPSs). The ground state energies at finite sizes and in the thermodynamic limit (with finite size scaling) are in good agreement with the state of art exact diagonalization study, and the energy differences between these two studies are of the order of 10^{-3}J_1 per site. At the largest bond dimension available D (D=9), we find a paramagnetic ground state without any valence bond solid order in the thermodynamic limit in the range of 0.5 < J_2/J_1 < 0.6, which implies the emergence of a spin-liquid phase. Furthermore, we investigate the topologically ordered nature of such a spin-liquid phase by measuring a nonzero topological entanglement entropy.

•Apr 25, 2012•7 citations

Possible spin liquid state in the spin 1/2 J1-J2 antiferromagnetic Heisenberg model on square lattice: A tensor product state approach

Ling Wang, Zheng‐Cheng Gu, Frank Verstraete +1 more

We study the spin 1/2 J1-J2 antiferromagnetic Heisenberg model on square lattice by using a recently proposed cluster update method for tensor product states (TPSs). The ground state energies in the thermodynamic limit are in good agreement with the state of art exact diagonalization study, and the energy differences between these two studies are of the order of 10J1 per site. For small bond dimension D (D < 5), we find a paramagnetic ground state with mixed columnar and staggered dimer orders of a magnitude 10−2 in the range of 0.47 < J2/J1 < 0.6, which agrees with previous results. For large bond dimension D (D ≥ 5), we observe all these orders vanish, implying the emergence of a spin liquid phase.

•Dec 15, 2011•2 citations

J1−J2 frustrated two-dimensional Heisenberg model: Random phase approximation and functional renormalization group

Johannes Reuther, P. Wölfle

We study the ground-state properties of the two-dimensional spin-1/2 ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on a square lattice, within diagrammatic approximations using an auxiliary-fermion formulation with exact projection. In a first approximation, we assume a phenomenological width of the pseudofermion spectral function to calculate the magnetization, susceptibilities, and the spin-correlation length within random-phase approximation, demonstrating the appearance of a paramagnetic phase between the N\'eel-ordered and Collinear-ordered phases, at sufficiently large pseudofermion damping. Second we use a functional renormalization-group formulation. We find that the conventional truncation scheme omitting three-particle and higher-order vertices is not sufficient. We therefore include self-energy renormalizations in the single-scale propagator as recently proposed by Katanin, to preserve Ward identities in a better way. We find N\'eel order at $g={J}_{2}/{J}_{1}\ensuremath{\lesssim}{g}_{\text{c}1}\ensuremath{\approx}0.4\dots{}0.45$ and Collinear order at $g\ensuremath{\gtrsim}{g}_{\text{c}2}\ensuremath{\approx}0.66\dots{}0.68$, which is in good agreement with results obtained by numerical studies. In the intervening quantum paramagnetic phase, we find enhanced columnar dimer and plaquette fluctuations of equal strength.

arXiv:0912.0860•Physical Review B•Apr 8, 2010•161 citations
We present numerical exact results for the ground state and the low-lying excitations for the spin-1/2 J1-J2 Heisenberg antiferromagnet on finite square lattices of up to N = 40 sites. Using finite-size extrapolation we determine the ground-state energy, the magnetic order parameters, the spin gap, the uniform susceptibility, as well as the spin-wave velocity and the spin stiffness as functions of the frustration parameter J2/J1. In agreement with the generally excepted scenario we find semiclassical magnetically ordered phases for J2 < J c1 2 and J2> J c2 2 separated by a gapful quantum paramagnetic phase. We estimate Jc1 2 ≈ 0.35J1 and J c2 2 ≈ 0.66J1. PACS codes: 75.10.Jm Quantized spin models 75.45.+j Macroscopic quantum phenomena in magnetic systems 1
Using the coupled cluster method for high orders of approximation and complementary exact diagonalization studies we investigate the ground state properties of the spin-1/2 J1−−J2{J}_{1}--{J}_{2}J1​−−J2​ frustrated Heisenberg antiferromagnet on the square lattice. We have calculated the ground state energy, the magnetic order parameter, the spin stiffness, and several generalized susceptibilities to probe magnetically disordered quantum valence-bond phases. We have found that the quantum critical points for both the N\'eel and collinear orders are ${J}_{2}^{c1}\ensuremath{\approx}(0.44\ifmmode\pm\else\textpm\fi{}0.01){J}_{1}$ and ${J}_{2}^{c2}\ensuremath{\approx}(0.59\ifmmode\pm\else\textpm\fi{}0.01){J}_{1}$, respectively, which are in good agreement with the results obtained by other approximations. In contrast to the recent study by [Sirker et al. Phys. Rev. B 73, 184420 (2006)], our data do not provide evidence for the transition from the N\'eel to the valence-bond solid state to be first order. Moreover, our results are in favor of the deconfinement scenario for that phase transition. We also discuss the nature of the magnetically disordered quantum phase.
We calculate the spin stiffness ${\mathrm{\ensuremath{\rho}}}_{\mathit{s}}$ for the frustrated spin-1/2 Heisenberg antiferromagnet on a square lattice by exact diagonalizations on finite clusters of up to 36 sites followed by extrapolations to the thermodynamic limit. For the nonfrustrated case, we find that ${\mathrm{\ensuremath{\rho}}}_{\mathit{s}}$=(0.183\ifmmode\pm\else\textpm\fi{}0.003)J1{\mathit{J}}_{1}J1​, in excellent agreement with the best results obtained by other means. Turning on frustration, the extrapolated stiffness vanishes for 0.4\ensuremath{\lesssim}J2{\mathit{J}}_{2}J2​/J1{\mathit{J}}_{1}J1​\ensuremath{\lesssim}0.6. In this intermediate region, the finite-size scaling works poorly---an additional sign that there is neither N\'eel nor collinear magnetic order. Using a hydrodynamic relation, and previous results for the transverse susceptibility, we also estimate the spin-wave velocity in the N\'eel-ordered region.
The ground state phase of a spin-$\frac{1}{2} {J}_{1}\text{\ensuremath{-}}{J}_{2}$ antiferromagnetic Heisenberg model on a square lattice around the maximally frustrated regime $({J}_{2}\ensuremath{\sim}0.5{J}_{1})$ has been debated for decades. Here we study this model using the cluster update algorithm for tensor-product states (TPSs). The ground state energies at finite sizes and in the thermodynamic limit (with finite size scaling) are in good agreement with exact diagonalization study. Through finite size scaling of the spin correlation function, we find the critical point J2c1=0.572(5)J1{J}_{2}^{{c}_{1}}=0.572(5){J}_{1}J2c1​​=0.572(5)J1​ and critical exponents $\ensuremath{\nu}=0.50(8),\phantom{\rule{0.16em}{0ex}}{\ensuremath{\eta}}_{s}=0.28(6)$. In the range of $0.572<{J}_{2}/{J}_{1}\ensuremath{\le}0.6$ we find a paramagnetic ground state with an exponentially decaying spin-spin correlation. Up to a $24\ifmmode\times\else\texttimes\fi{}24$ system size, we observe power law decaying dimer-dimer and plaquette-plaquette correlations with an anomalous plaquette scaling exponent ${\ensuremath{\eta}}_{p}=0.24(1)$ and an anomalous columnar scaling exponent ${\ensuremath{\eta}}_{c}=0.28(1)$ at J2/J1=0.6{J}_{2}/{J}_{1}=0.6J2​/J1​=0.6. These results are consistent with a potential gapless U(1)U(1)U(1) spin-liquid phase. However, since the U(1)U(1)U(1) spin liquid is unstable due to the instanton effect, a valence bond solid order with very small amplitude might develop in the thermodynamic limit. Thus, our numerical results strongly indicate a deconfined quantum critical point at J2c1{J}_{2}^{{c}_{1}}J2c1​​. Remarkably, all the observed critical exponents are consistent with the $J\text{\ensuremath{-}}Q$ model.
The ground state phase of spin-1/2 J_1-J_2 antiferromagnetic Heisenberg model on square lattice in the maximally frustrated regime (J_2 ~ 0.5J_1) has been debated for decades. Here we study this model by using a recently proposed novel numerical method - the cluster update algorithm for tensor product states (TPSs). The ground state energies at finite sizes and in the thermodynamic limit (with finite size scaling) are in good agreement with the state of art exact diagonalization study, and the energy differences between these two studies are of the order of 10^{-3}J_1 per site. At the largest bond dimension available D (D=9), we find a paramagnetic ground state without any valence bond solid order in the thermodynamic limit in the range of 0.5 < J_2/J_1 < 0.6, which implies the emergence of a spin-liquid phase. Furthermore, we investigate the topologically ordered nature of such a spin-liquid phase by measuring a nonzero topological entanglement entropy.
We study the spin 1/2 J1-J2 antiferromagnetic Heisenberg model on square lattice by using a recently proposed cluster update method for tensor product states (TPSs). The ground state energies in the thermodynamic limit are in good agreement with the state of art exact diagonalization study, and the energy differences between these two studies are of the order of 10J1 per site. For small bond dimension D (D < 5), we find a paramagnetic ground state with mixed columnar and staggered dimer orders of a magnitude 10−2 in the range of 0.47 < J2/J1 < 0.6, which agrees with previous results. For large bond dimension D (D ≥ 5), we observe all these orders vanish, implying the emergence of a spin liquid phase.
We study the ground-state properties of the two-dimensional spin-1/2 ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on a square lattice, within diagrammatic approximations using an auxiliary-fermion formulation with exact projection. In a first approximation, we assume a phenomenological width of the pseudofermion spectral function to calculate the magnetization, susceptibilities, and the spin-correlation length within random-phase approximation, demonstrating the appearance of a paramagnetic phase between the N\'eel-ordered and Collinear-ordered phases, at sufficiently large pseudofermion damping. Second we use a functional renormalization-group formulation. We find that the conventional truncation scheme omitting three-particle and higher-order vertices is not sufficient. We therefore include self-energy renormalizations in the single-scale propagator as recently proposed by Katanin, to preserve Ward identities in a better way. We find N\'eel order at $g={J}_{2}/{J}_{1}\ensuremath{\lesssim}{g}_{\text{c}1}\ensuremath{\approx}0.4\dots{}0.45$ and Collinear order at $g\ensuremath{\gtrsim}{g}_{\text{c}2}\ensuremath{\approx}0.66\dots{}0.68$, which is in good agreement with results obtained by numerical studies. In the intervening quantum paramagnetic phase, we find enhanced columnar dimer and plaquette fluctuations of equal strength.