Four-Spin Terms and the Origin of the Chiral Spin Liquid in Mott Insulators on the Triangular Lattice

Quantum spin liquid signatures in monolayer 1T-NbSe2

Quantum spin liquids (QSLs) are in a quantum disordered state that is highly entangled and has fractional excitations. As a highly sought-after state of matter, QSLs were predicted to host spinon excitations and to arise in frustrated spin systems with large quantum fluctuations. Here we report on the experimental observation and theoretical modeling of QSL signatures in monolayer 1T-NbSe2, which is a newly emerging two-dimensional material that exhibits both charge-density-wave (CDW) and correlated insulating behaviors. By using scanning tunneling microscopy and spectroscopy (STM/STS), we confirm the presence of spin fluctuations in monolayer 1T-NbSe2 by observing the Kondo resonance as monolayer 1T-NbSe2 interacts with metallic monolayer 1H-NbSe2. Subsequent STM/STS imaging of monolayer 1T-NbSe2 at the Hubbard band energy further reveals a long-wavelength charge modulation, in agreement with the spinon modulation expected for QSLs. By depositing manganese-phthalocyanine (MnPc) molecules with spin S = 3/2 onto monolayer 1T-NbSe2, new STS resonance peaks emerge at the Hubbard band edges of monolayer 1T-NbSe2. This observation is consistent with the spinon Kondo effect induced by a S = 3/2 magnetic impurity embedded in a QSL. Taken together, these experimental observations indicate that monolayer 1T-NbSe2 is a new promising QSL material.

Ising Ladder with Four-Spin Plaquette Interaction in a Transverse Magnetic Field

The spin-1/2 quantum transverse Ising model, defined on a ladder structure, with nearest-neighbor and four-spin interaction on a plaquette, was studied by using exact diagonalization on finite ladders together with finite-size-scaling procedures. The quantum phase transition between the ferromagnetic and paramagnetic phases has then been obtained by extrapolating the data to the thermodynamic limit. The critical transverse field decreases as the antiferromagnetic four-spin interaction increases and reaches a multicritical point. However, the exact diagonalization approach was not able to capture the essence of the dimer phase beyond the multicritical transition.

Chasing the spin gap through the phase diagram of a frustrated Mott insulator

The quest for entangled spin excitations has stimulated intense research on frustrated magnetic systems. For almost two decades, the triangular-lattice Mott insulator κ-(BEDT-TTF)₂Cu₂(CN)₃ has been one of the hottest candidates for a gapless quantum spin liquid with itinerant spinons. Very recently, however, this scenario was overturned as electron-spin-resonance (ESR) studies unveiled a spin gap, calling for reevaluation of the magnetic ground state. Here we achieve a precise mapping of this spin-gapped phase through the Mott transition by ultrahigh-resolution strain tuning. Our transport experiments reveal a reentrance of charge localization below T^⋆ = 6 K associated with a gap size of 30-50 K. The negative slope of the insulator-metal boundary, dT^⋆/dp < 0, evidences the low-entropy nature of the spin-singlet ground state. By tuning the enigmatic '6K anomaly' through the phase diagram of κ-(BEDT-TTF)₂Cu₂(CN)₃, we identify it as the transition to a valence-bond-solid phase, in agreement with previous thermal expansion and magnetic resonance studies. This spin-gapped insulating state persists at T → 0 until unconventional superconductivity and metallic transport proliferate.

Simulating Chiral Spin Liquids with Projected Entangled-Pair States

Doubts have been raised on the representation of chiral spin liquids exhibiting topological order in terms of projected entangled pair states (PEPSs). Here, starting from a simple spin-1/2 chiral frustrated Heisenberg model, we show that a faithful representation of the chiral spin liquid phase is in fact possible in terms of a generic PEPS upon variational optimization. We find a perfectly chiral gapless edge mode and a rapid decay of correlation functions at short distances consistent with a bulk gap, concomitant with a gossamer long-range tail originating from a PEPS bulk-edge correspondence. For increasing bond dimension, (i) the rapid decrease of spurious features-SU(2) symmetry breaking and long-range tails in correlations-together with (ii) a faster convergence of the ground state energy as compared to state-of-the-art cylinder matrix-product state simulations involving far more variational parameters, prove the fundamental relevance of the PEPS ansatz for simulating systems with chiral topological order.

Nonlocal Interactions in Moiré Hubbard Systems

Moiré materials formed in two-dimensional semiconductor heterobilayers are quantum simulators of Hubbard-like physics with unprecedented electron density and interaction strength tunability. Compared to atomic scale Hubbard-like systems, electrons or holes in moiré materials are less strongly attracted to their effective lattice sites because these are defined by finite-depth potential extrema. As a consequence, nonlocal interaction terms like interaction-assisted hopping and intersite exchange are more relevant. We theoretically demonstrate the possibility of tuning the strength of these coupling constants to favor unusual states of matter, including spin liquids, insulating ferromagnets, and superconductors.

Quantum Phases of Transition Metal Dichalcogenide Moiré Systems

Moiré systems provide a rich platform for studies of strong correlation physics. Recent experiments on heterobilayer transition metal dichalcogenide Moiré systems are exciting in that they manifest a relatively simple model system of an extended Hubbard model on a triangular lattice. Inspired by the prospect of the hetero-transition metal dichalcogenide Moiré system's potential as a solid-state-based quantum simulator, we explore the extended Hubbard model on the triangular lattice using the density matrix renormalization group. Specifically, we explore the two-dimensional phase space spanned by the key tuning parameters in the extended Hubbard model, namely, the kinetic energy strength and the further-range Coulomb interaction strengths. We find competition between Fermi fluid, chiral spin liquid, spin density wave, and charge order. In particular, our finding of the optimal further-range interaction for the chiral correlation presents a tantalizing possibility.

Spin Vortex Crystal Order in Organic Triangular Lattice Compound

Organic salts represent an ideal experimental playground for studying the interplay between magnetic and charge degrees of freedom, which has culminated in the discovery of several spin-liquid candidates such as κ-(ET)_{2}Cu_{2}(CN)_{3} (κ-Cu). Recent theoretical studies indicate the possibility of chiral spin liquids stabilized by ring exchange, but the parent states with chiral magnetic order have not been observed in this material family. In this Letter, we discuss the properties of the recently synthesized κ-(BETS)_{2}Mn[N(CN)_{2}]_{3} (κ-Mn). Based on analysis of specific heat, magnetic torque, and NMR measurements combined with ab initio calculations, we identify a spin-vortex crystal order. These observations definitively confirm the importance of ring exchange in these materials and support the proposed chiral spin-liquid scenario for triangular lattice organics.

Thermal Properties and Instability of a U(1) Spin Liquid on the Triangular Lattice

We study the effect of Dzyaloshinskii-Moriya (DM) interaction on the triangular lattice U(1) quantum spin liquid (QSL) which is stabilized by ring-exchange interactions. A weak DM interaction introduces a staggered flux to the U(1) QSL state and changes the density of states at the spinon Fermi surface. If the DM vector contains in-plane components, then the spinons gain nonzero Berry phase. The resultant thermal conductances κ_{xx} and κ_{xy} qualitatively agree with the experimental results on the material EtMe_{3}Sb[Pd(dmit)_{2}]_{2}. Furthermore, owing to perfect nesting of the Fermi surface, a spin density wave state is triggered by larger DM interactions. On the other hand, when the ring-exchange interaction decreases, another antiferromagnetic (AFM) phase with 120° order shows up which is proximate to a U(1) Dirac QSL. We discuss the difference of the two AFM phases from their static structure factors and excitation spectra.