Motivated by experimental studies that have found signatures of a quantum spin liquid phase in organic crystals whose structure is well described by the two-dimensional triangular lattice, we study the Hubbard model on this lattice at half filling using the infinite-system density matrix renormalization group (iDMRG) method. On infinite cylinders with finite circumference, we identify an intermediate phase between observed metallic behavior at low interaction strength and Mott insulating spin-ordered behavior at strong interactions. Chiral ordering from spontaneous breaking of time-reversal symmetry, a fractionally quantized spin Hall response, and characteristic level statistics in the entanglement spectrum in the intermediate phase provide strong evidence for the existence of a chiral spin liquid in the full two-dimensional limit of the model.
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.
Two-dimensional density-matrix renormalization group method is employed to examine the ground-state phase diagram of the Hubbard model on the triangular lattice at half-filling. The calculation reveals two discontinuities in the double occupancy with increasing the repulsive Hubbard interaction
U at ${U}_{\mathrm{c}1}\ensuremath{\sim}7.8t$ and ${U}_{\mathrm{c}2}\ensuremath{\sim}9.9t$ (
t being the hopping integral), indicating that there are three phases separated by first-order transitions. The absence of any singularity in physical quantities for $0\ensuremath{\le}U<{U}_{\mathrm{c}1}$ implies a metallic phase in this regime. For
U>Uc2, the local spin density induced by an applied pinning magnetic field exhibits a three sublattice feature, which is compatible with the ${120}^{\ensuremath{\circ}}$ N\'eel-ordered state realized in the limit of $U\ensuremath{\rightarrow}\ensuremath{\infty}$. For
Uc1<U<Uc2, a response to the applied pinning magnetic field is comparable to that in the metallic phase with a relatively large spin correlation length, but showing neither valence bond nor chiral magnetic order, which therefore resembles gapless spin liquid. However, the spin structure factor for the intermediate phase exhibits the maximum at the
K and ${\mathrm{K}}^{\ensuremath{'}}$ points in the momentum space, which is not compatible to spin liquid with a large spinon Fermi surface. The calculation also finds that the pairing correlation function monotonically decreases with increasing
U and thus the superconductivity is unlikely in the intermediate phase.
Quantum spin liquids emerging near the metal-insulator transition have been a longstanding issue in physics. Here, the author identifies a chiral spin liquid in the half-filled triangular-lattice Hubbard model using state-of-the-art non-Abelian density matrix renormalization. By computing the spin chirality order parameter with unprecedented system sizes and high accuracy, the findings extrapolate toward the two-dimensional limit, impacting theoretical research and experimental realizations in ultracold-atom optical lattices and transition-metal dichalcogenide moir\'e materials.
A chiral spin liquid (CSL) phase has been recently reported in the Hubbard model on a triangular lattice. It emerges in an intermediate coupling regime at half-filling, which is sandwiched between a ${120}^{\ensuremath{\circ}}$ antiferromagnetic (AFM) phase and a metallic phase as a function of on-site repulsion
U. In this work, we examine the mechanism of the CSL and complex phase diagram via analytic analysis and numerical density matrix renormalization group (DMRG) method. First, we identify an exact Berry-phase-like sign structure in the partition function of the model at arbitrary
U, which is originated from the Fermi sign structure at
U=0.The spin and charge degrees of freedom are generally entangled via a singular phase string in such a many-body sign structure. In the large-
U limit, the suppression of the charge fluctuation at half-filling can render the phase string ineffective, resulting in the AFM order. However, if one precisely switches off such a phase string in the DMRG calculation, the ${120}^{\ensuremath{\circ}}$ AFM is shown to survive all the way down to a much weaker
U without the emergence of the CSL and metallic phases. It indicates that the phase-string sign structure plays the key role to mediate the mutual interaction between the charge and spin fluctuations, which results in the CSL and metallic phases at finite
U. General implications for the Mott physics in the Hubbard model will also be discussed.