No Arabic abstract
We use the state-of-the-art tensor network state method, specifically, the finite projected entangled pair state (PEPS) algorithm, to simulate the global phase diagram of spin-$1/2$ $J_1$-$J_2$ Heisenberg model on square lattices up to $24times 24$. We provide very solid evidences to show that the nature of the intermediate nonmagnetic phase is a gapless quantum spin liquid (QSL), whose spin-spin and dimer-dimer correlations both decay with a power law behavior. There also exists a valence-bond solid (VBS) phase in a very narrow region $0.56lesssim J_2/J_1leq0.61$ before the system enters the well known collinear antiferromagnetic phase. We stress that our work gives rise to the first solid PEPS results beyond the well established density matrix renormalization group (DMRG) through one-to-one direct benchmark for small system sizes. Thus our numerical evidences explicitly demonstrate the huge power of PEPS for solving long-standing 2D quantum many-body problems. The physical nature of the discovered gapless QSL and potential experimental implications are also addressed.
The spin-1/2 $J_1$-$J_2$ Heisenberg model on square lattices are investigated via the finite projected entangled pair states (PEPS) method. Using the recently developed gradient optimization method combining with Monte Carlo sampling techniques, we are able to obtain the ground states energies that are competitive to the best results. The calculations show that there is no Neel order, dimer order and plaquette order in the region of 0.42 $lesssim J_2/J_1lesssim$ 0.6, suggesting a single spin liquid phase in the intermediate region. Furthermore, the calculated staggered spin, dimer and plaquette correlation functions all have power law decay behaviours, which provide strong evidences that the intermediate nonmagnetic phase is a single gapless spin liquid state.
Liu et al. [Phys.Rev.B 98, 241109 (2018)] used Monte Carlo sampling of the physical degrees of freedom of a Projected Entangled Pair State (PEPS) type wave function for the $S=1/2$ frustrated $J_1$-$J_2$ Heisenberg model on the square lattice and found a non-magnetic state argued to be a gapless spin liquid when the coupling ratio $g=J_2/J_1$ is in the range $g in [0.42,0.6]$. Here we show that their definition of the order parameter for another candidate ground state within this coupling window---a spontaneously dimerized state---is problematic. The order parameter as defined will not detect dimer order when lattice symmeties are broken due to open boundaries or asymmetries originating from the calculation itself. Thus, a dimerized phase for some range of $g$ cannot be excluded (and is likely based on several other recent works).
Strongly correlated systems with geometric frustrations can host the emergent phases of matter with unconventional properties. Here, we study the spin $S = 1$ Heisenberg model on the honeycomb lattice with the antiferromagnetic first- ($J_1$) and second-neighbor ($J_2$) interactions ($0.0 leq J_2/J_1 leq 0.5$) by means of density matrix renormalization group (DMRG). In the parameter regime $J_2/J_1 lesssim 0.27$, the system sustains a N{e}el antiferromagnetic phase. At the large $J_2$ side $J_2/J_1 gtrsim 0.32$, a stripe antiferromagnetic phase is found. Between the two magnetic ordered phases $0.27 lesssim J_2/J_1 lesssim 0.32$, we find a textit{non-magnetic} intermediate region with a plaquette valence-bond order. Although our calculations are limited within $6$ unit-cell width on cylinder, we present evidence that this plaquette state could be a strong candidate for this non-magnetic region in the thermodynamic limit. We also briefly discuss the nature of the quantum phase transitions in the system. We gain further insight of the non-magnetic phases in the spin-$1$ system by comparing its phase diagram with the spin-$1/2$ system.
We study the phase diagram of the frustrated Heisenberg model on the triangular lattice with nearest and next-nearest neighbor spin exchange coupling, on 3-leg ladders. Using the density-matrix renormalization-group method, we obtain the complete phase diagram of the model, which includes quasi-long-range $120^circ$ and columnar order, and a Majumdar-Ghosh phase with short-ranged correlations. All these phases are non-chiral and planar. We also identify the nature of phase transitions.
We numerically study the Heisenberg models on triangular lattices by extending it from the simplest equilateral lattice with only the nearest-neighbor exchange interaction. We show that, by including an additional weak next-nearest-neighbor interaction, a quantum spin-liquid phase is stabilized against the antiferromagnetic order. The spin gap (triplet excitation gap) and spin correlation at long distances decay algebraically with increasing system size at the critical point between the antiferromagnetic phase and the spin-liquid phase. This algebraic behavior continues in the spin-liquid phase as well, indicating the presence of an unconventional critical (algebraic spin-liquid) phase characterized by the dynamical and anomalous critical exponents $z+etasim1$. Unusually small triplet and singlet excitation energies found in extended points of the Brillouin zone impose constraints on this algebraic spin liquid.