Spin systems with frustrated anisotropic interactions are of significant interest due to possible exotic ground states. We have explored their phase diagram on a nearest-neighbor triangular lattice using the density-matrix renormalization group and mapped out the topography of the region that can harbor a spin liquid. We find that this spin-liquid phase is continuously connected to a previously discovered spin-liquid phase of the isotropic J_{1}-J_{2} model. The two limits show nearly identical spin correlations, making the case that their respective spin liquids are isomorphic to each other.
Spin systems with frustrated anisotropic interactions are of significant interest due to possible exotic ground states. We have explored their phase diagram on a nearest-neighbor triangular lattice using the density-matrix renormalization group and mapped out the topography of the region that can harbor a spin liquid. We find that this spin-liquid phase is continuously connected to a previously discovered spin-liquid phase of the isotropic J_{1}-J_{2} model. The two limits show nearly identical spin correlations, making the case that their respective spin liquids are isomorphic to each other.
We study the spin-21 Heisenberg model on the triangular lattice with antiferromagnetic first- (J1) and second- (J2) nearest-neighbor interactions using density matrix renormalization group. By studying the spin correlation function, we find a ${120}^{\ensuremath{\circ}}$ magnetic order phase for ${J}_{2}\ensuremath{\lesssim}0.07{J}_{1}$ and a stripe antiferromagnetic phase for ${J}_{2}\ensuremath{\gtrsim}0.15{J}_{1}$. Between these two phases, we identify a spin-liquid region characterized by exponential decaying spin and dimer correlations, as well as large spin singlet and triplet excitation gaps on finite-size systems. We find two near degenerating ground states with distinct properties in two sectors, which indicates more than one spin-liquid candidate in this region. While the sector with spinons is found to respect time reversal symmetry, the even sector without spinons breaks such a symmetry for finite-size systems. Furthermore, we detect the signature of the fractionalization by following the evolution of different ground states with inserting spin flux into the cylinder system. Moreover, by tuning the anisotropic bond coupling, we explore the nature of the spin-liquid phase and find the optimal parameter region for gapped Z2 spin liquids.
The antiferromagnetic Heisenberg model on the triangular lattice is perhaps the best known example of frustrated magnets, but it orders at low temperatures. Recent density matrix renormalization group (DMRG) calculations find that the next nearest neighbor interaction J_{2} enhances the frustration, and it leads to a spin liquid for J_{2}/J_{1}∈(0.08,0.15). In addition, a DMRG study of a dipolar Heisenberg model with longer range interactions gives evidence for a spin liquid at a small dipole tilting angle θ∈[0,10°). In both cases, the putative spin liquid region appears to be small. Here, we show that for the triangular lattice dipolar Heisenberg model, a robust quantum paramagnetic phase exists in a surprisingly wide region, θ∈[0,54°), for dipoles tilted along the lattice diagonal direction. We obtain the phase diagram of the model by functional renormalization group (RG), which treats all magnetic instabilities on equal footing. The quantum paramagnetic phase is characterized by a smooth continuous flow of vertex functions and spin susceptibility down to the lowest RG scale, in contrast to the apparent breakdown of RG flow in phases with stripe or spiral order. Our finding points to a promising direction to search for quantum spin liquids in ultracold dipolar molecules.
In a recent study [Phys. Rev. X 10, 021042 (2020)], we showed using large-scale density matrix renormalization group (DMRG) simulations on infinite cylinders that the triangular lattice Hubbard model has a chiral spin liquid phase. In this work, we introduce hopping anisotropy in the model, making one of the three distinct bonds on the lattice stronger or weaker compared with the other two. We implement the anisotropy in two inequivalent ways, one which respects the mirror symmetry of the cylinder and one which breaks this symmetry. In the full range of anisotropy, from the square lattice to weakly coupled one-dimensional chains, we find a variety of phases. Near the isotropic limit we find the three phases identified in our previous work: metal, chiral spin liquid, and ${120}^{\ensuremath{\circ}}$ spiral order; we note that a recent paper suggests the apparently metallic phase may actually be a Luther-Emery liquid, which would also be in agreement with our results. When one bond is weakened by a relatively small amount, the ground state quickly becomes the square lattice N\'eel order. When one bond is strengthened, the story is much less clear, with the phases that we find depending on the orientation of the anisotropy and on the cylinder circumference. While our work is to our knowledge the first DMRG study of the anisotropic triangular lattice Hubbard model, the overall phase diagram we find is broadly consistent with that found previously using other methods, such as variational Monte Carlo and dynamical mean field theory.
We investigate a highly frustrated spin-1 model on the triangular lattice, with nearest- and next-nearest-neighbor antiferromagnetic $S\ifmmode\cdot\else\textperiodcentered\fi{}S$ interactions and nearest-neighbor ${(S\ifmmode\cdot\else\textperiodcentered\fi{}S)}^{2}$ interactions. Using the density matrix renormalization group technique, we find three magnetically ordered phases, namely ${120}^{\ensuremath{\circ}}$ spiral order, stripe order, and tetrahedral order, as well as two spin nematic phases: ferroquadrupolar and antiferroquadrupolar. While our data could be consistent with a spin liquid phase between the ${120}^{\ensuremath{\circ}}$ spiral and antiferroquadrupolar orders, the more likely scenario is a direct continuous transition between these two orders.
We study the S=1/2 Heisenberg model on the triangular lattice with nearest- and next-nearest-neighbor interactions J1 and J2 with the density matrix renormalization group, on long open cylinders with widths up to nine lattice spacings. In an intermediate J2 region $0.06\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\lesssim}0.17$, we find evidence for a spin liquid (SL) state with short range spin-spin, bond-bond, and chiral correlation lengths, bordered by a classical ${120}^{\ensuremath{\circ}}$ N\'eel ordered state at small J2 and by a two sublattice collinear magnetically ordered state at larger J2. Focusing on J2/J1=0.1, we find a number of signatures of a gapped SL phase.
Topological Spin Liquid Phase in the Spin-1/2 Triangular Heisenberg Model
We study the spin-1/2 Heisenberg model on the triangular lattice with the antiferromagnetic first (J1) and second nearest-neighbor (J2) interactions using density matrix renormalization group. By studying the spin correlation function, we find a 120◦ magnetic order phase for J2 . 0.08J1 and a stripe antiferromagnetic phase for J2 & 0.16J2. Between these two phases, we identify a gapped spin liquid phase characterized by the exponential decaying spin and dimer correlations, as well as finite spin singlet and triplet excitation gaps. We also find two topological degenerating ground states in two topological sectors. Furthermore, we detect the signature of the fractionalization by following the evolution of different ground states with inserting spin flux into the cylinder system. While the sector with spinon is found to respect the time reveral symmetry, the even sector without a spinon may break such a symmetry. We discuss the implications of numerical results for the nature of the topological spin liquid phase.