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.
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 perform highly accurate density matrix renormalization group (DMRG) simulations to investigate the ground-state properties of the spin-
1/2 antiferromagnetic square lattice Heisenberg
J1-
J2 model. Based on studies of numerous long cylinders with circumferences of up to 14 lattice spacings, we obtain strong evidence for a topological quantum spin liquid state in the region $0.41\ensuremath{\le}{J}_{2}/{J}_{1}\ensuremath{\le}0.62$, separating conventional N\'eel and striped antiferromagnetic states for smaller and larger
J2/J1, respectively. The quantum spin liquid is characterized numerically by the absence of magnetic or valence bond solid order, and nonzero singlet and triplet energy gaps. Furthermore, we positively identify its topological nature by measuring a nonzero topological entanglement entropy $\ensuremath{\gamma}=0.70\ifmmode\pm\else\textpm\fi{}0.02$, extremely close to $\ensuremath{\gamma}=\mathrm{ln}(2)\ensuremath{\approx}0.69$ (expected for a
Z2 quantum spin liquid) and a nontrivial finite size dimerization effect depending upon the parity of the circumference of the cylinder. We also point out that a valence bond solid, and indeed any discrete symmetry breaking state, would be expected to show a constant correction to the entanglement entropy of opposite sign to the topological entanglement entropy.
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 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 study the ground-state phase diagram of the quantum ${J}_{1}\ensuremath{-}{J}_{2}$ model on a square lattice by means of an entangled-plaquette variational Ansatz. In the range $0\ensuremath{\leqslant}{J}_{2}/{J}_{1}\ensuremath{\leqslant}1$, we find classical magnetic order of N\'eel and collinear type, for ${J}_{2}/{J}_{1}\ensuremath{\lesssim}0.5$ and ${J}_{2}/{J}_{1}\ensuremath{\gtrsim}0.6$, respectively. For intermediate values of
J2/J1 the ground state is a spin liquid (i.e., paramagnetic with no valence-bond crystalline order). Our estimates of the entanglement entropy show that such a spin liquid is topological.
The nature of the quantum spin liquid (QSL) phase in the frustrated
J1-
J2 Heisenberg model on the triangular lattice remains an open and actively debated problem. In this work, we employ the infinite projected entangled-pair state (PEPS) to systematically investigate the model under different symmetry constraints. Our simulations reveal a direct transition from the
120∘ Néel state to a putative QSL at
J2/J1≈0.08, signaled by the collapse of magnetic order. We further show that, through either an appropriate unitary rotation or spontaneous spin long-range order, the stripe antiferromagnetic phase can also be accurately captured within the infinite PEPS framework. A central focus of our study is the role played by the PEPS symmetry in approximating the QSL ground-state sandwiched between the two magnetic phases. We first found that a fully-symmetric topological
Z2 Resonating Valence Bond state, which can be written as a simple PEPS with bond dimension
D=3, exhibits a reasonably good variational energy. Motivated by this finding, we have further constructed generic
Z2-symmetric PEPS of larger bond dimension (up to
D=7). We found that, under wavefunction optimization, spinons condense and, simultaneously, topological vison excitations get confined, hence precluding
Z2 topological order. This strongly indicates the gapless (or critical) nature of the QSL phase, which is most naturally consistent with a U(1) Dirac spin liquid scenario.