Using an SU(2) invariant finite-temperature tensor network algorithm, we provide strong numerical evidence in favor of an Ising transition in the collinear phase of the spin-1/2 J_{1}-J_{2} Heisenberg model on the square lattice. In units of J_{2}, the critical temperature reaches a maximal value of T_{c}/J_{2}≃0.18 around J_{2}/J_{1}≃1.0. It is strongly suppressed upon approaching the zero-temperature boundary of the collinear phase J_{2}/J_{1}≃0.6, and it vanishes as 1/log(J_{2}/J_{1}) in the large J_{2}/J_{1} limit, as predicted by Chandra et al., [Phys. Rev. Lett. 64, 88 (1990)PRLTAO0031-900710.1103/PhysRevLett.64.88]. Enforcing the SU(2) symmetry is crucial to avoid the artifact of finite-temperature SU(2) symmetry breaking of U(1) algorithms, opening new perspectives in the investigation of the thermal properties of quantum Heisenberg antiferromagnets.
Using state of the art tensor network computations combined with spin wave theory, we compute the finite temperature phase diagram of the spin 1/2
J1−J2 Heisenberg model on the square lattice with ferromagnetic
J1<0 and antiferromagnetic
J2>0. We localize the first order line originating from the first order zero temperature point as well as the second order thermal transition at large
J2 associated with spontaneous lattice rotation symmetry breaking. This second order line terminates at a critical end point where it meets the first order line. We found no evidence for an intermediate phase.
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 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. 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
J2c1. Remarkably, all the observed critical exponents are consistent with the $J\text{\ensuremath{-}}Q$ model.
We apply the plaquette renormalization scheme of tensor network states [Phys. Rev. E 83, 056703 (2011)] to study the spin-1/2 frustrated Heisenberg
J1-
J2 model on an $L\ifmmode\times\else\texttimes\fi{}L$ square lattice with
L = 8, 16, and 32. By treating tensor elements as variational parameters, we obtain the ground states for different
J2/J1 values, and investigate staggered magnetizations, nearest-neighbor spin-spin correlations, and plaquette order parameters. In addition to the well-known N\'eel order and collinear order at low and high
J2/J1, we observe a plaquettelike order at ${J}_{2}/{J}_{1}\ensuremath{\approx}0.5$. A continuous transition between the N\'eel order and the plaquettelike order near ${J}_{2}^{{c}_{1}}\ensuremath{\approx}0.40{J}_{1}$ is observed. The collinear order emerges at ${J}_{2}^{{c}_{2}}\ensuremath{\approx}0.62{J}_{1}$ through a first-order phase transition.
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 Heisenberg model on the square lattice with the antiferromagnetic nearest-neighbor
J1 and the next-nearest-neighbor
J2 couplings by using the infinite projected entangled-pair state (iPEPS) ansatz and the density matrix renormalization group (DMRG) calculation. The iPEPS simulation, which studies the model directly in the thermodynamic limit, finds a crossing of the ground state from the N\'eel magnetic state to the stripe magnetic state at ${J}_{2}/{J}_{1}\ensuremath{\simeq}0.549$, showing a direct phase transition. In the finite-size DMRG calculation on the cylinder geometry up to the cylinder width
Ly=10, we find that around the same critical point the N\'eel and the stripe orders are strongly suppressed, which implies the absence of an intermediate phase. Both calculations identify that the stripe order comes with a first-order transition at ${J}_{2}/{J}_{1}\ensuremath{\simeq}0.549$. Our results indicate that unlike the spin-1/2 ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ square model, quantum fluctuations in the spin-1 model may not be strong enough to stabilize an intermediate nonmagnetic phase.
The classical Heisenberg model in two spatial dimensions constitutes one of the most paradigmatic spin models, taking an important role in statistical and condensed matter physics to understand magnetism. Still, despite its paradigmatic character and the widely accepted ban of a (continuous) spontaneous symmetry breaking, controversies remain whether the model exhibits a phase transition at finite temperature. Importantly, the model can be interpreted as a lattice discretization of the
O(3) non-linear sigma model in
1+1 dimensions, one of the simplest quantum field theories encompassing crucial features of celebrated higher-dimensional ones (like quantum chromodynamics in
3+1 dimensions), namely the phenomenon of asymptotic freedom. This should also exclude finite-temperature transitions, but lattice effects might play a significant role in correcting the mainstream picture. In this work, we make use of state-of-the-art tensor network approaches, representing the classical partition function in the thermodynamic limit over a large range of temperatures, to comprehensively explore the correlation structure for Gibbs states. By implementing an
SU(2) symmetry in our two-dimensional tensor network contraction scheme, we are able to handle very large effective bond dimensions of the environment up to
χEeff∼1500, a feature that is crucial in detecting phase transitions. With decreasing temperatures, we find a rapidly diverging correlation length, whose behaviour is apparently compatible with the two main contradictory hypotheses known in the literature, namely a finite-
T transition and asymptotic freedom, though with a slight preference for the second.