Gapless Spin Liquid Research Articles

Possible gapless spin liquid in the rare-earth kagome lattice magnet Tm3Sb3Zn2O14

We report the thermodynamic and muon spin relaxation ($\mu$SR) evidences for a possible gapless spin liquid in Tm$_{3}$Sb$_{3}$Zn$_{2}$O$_{14}$, with the rare-earth ions Tm$^{3+}$ forming a two-dimensional kagom\'{e} lattice. We extract the magnetic specific heat of Tm$_{3}$Sb$_{3}$Zn$_{2}$O$_{14}$ by subtracting the phonon contribution of the non-magnetic isostructural material La$_{3}$Sb$_{3}$Zn$_{2}$O$_{14}$. We obtain a clear linear-$T$ temperature dependence of magnetic specific heat at low temperatures, with the heat-capacity coefficient $\gamma$ = 26.6(1) mJ mol-Tm$^{-1}$ K$^{-2}$ in the zero temperature limit. No long-range magnetic order was observed down to 0.35 K in the heat capacity measurement. A broad hump around 9 K was observed in the zero field magnetic specific heat and is gradually suppressed in the magnetic fields that also create a spin gap in the specific heat. The absence of magnetic order is further confirmed by the $\mu$SR measurement down to 20 mK. We find that the spin-lattice relaxation time remains constant down to the lowest temperatures. We point out that the physics in Tm$_{3}$Sb$_{3}$Zn$_{2}$O$_{14}$ is fundamentally different from the Cu-based herbertsmithite and propose spin liquid ground states with non-Kramers doublets on the kagom\'{e} lattice to account for the experimental results.

Spectral signatures of fractionalization in the frustrated Heisenberg model on the square lattice

We employ a variational Monte Carlo approach to efficiently obtain the dynamical structure factor for the spin-1/2 $J_1-J_2$ Heisenberg model on the square lattice. Upon increasing the frustrating ratio $J_2/J_1$, the ground state undergoes a continuous transition from a N\'eel antiferromagnet to a $\mathbb{Z}_{2}$ gapless spin liquid. We identify the characteristic spectral features in both phases and highlight the existence of a broad continuum of excitations in the proximity of the spin-liquid phase. The magnon branch, which dominates the spectrum of the unfrustrated Heisenberg model, gradually loses its spectral weight, thus releasing nearly-deconfined spinons, whose signatures are visible even in the magnetically ordered state. Our results show how free spinons emerge across a quantum critical point, providing evidence for the fractionalization of magnons into deconfined spinons.

Critical Level Crossings and Gapless Spin Liquid in the Square-Lattice Spin-1/2 J_{1}-J_{2} Heisenberg Antiferromagnet.

We use the density matrix renormalization group method to calculate several energy eigenvalues of the frustrated S=1/2 square-lattice J_{1}-J_{2} Heisenberg model on 2L×L cylinders with L≤10. We identify excited-level crossings versus the coupling ratio g=J_{2}/J_{1} and study their drifts with the system size L. The lowest singlet-triplet and singlet-quintuplet crossings converge rapidly (with corrections ∝L^{-2}) to different g values, and we argue that these correspond to ground-state transitions between the Néel antiferromagnet and a gapless spin liquid, at g_{c1}≈0.46, and between the spin liquid and a valence-bond solid at g_{c2}≈0.52. Previous studies of order parameters were not able to positively discriminate between an extended spin liquid phase and a critical point. We expect level-crossing analysis to be a generically powerful tool in density matrix renormalization group studies of quantum phase transitions.

Spinon Fermi Surface in a Cluster Mott Insulator Model on a Triangular Lattice and Possible Application to 1T-TaS_{2}.

1T-TaS_{2} is a cluster Mott insulator on the triangular lattice with 13Ta atoms forming a star of David cluster as the unit cell. We derive a two-dimensional XXZ spin-1/2 model with a four-spin ring exchange term to describe the effective low energy physics of a monolayer 1T-TaS_{2}, where the effective spin-1/2 degrees of freedom arises from the Kramers degenerate spin-orbital states on each star of David. A large scale density matrix renormalization group simulation is further performed on this effective model and we find a gapless spin liquid phase with a spinon Fermi surface at a moderate to large strength region of the four-spin ring exchange term. All peaks in the static spin structure factor are found to be located on the "2k_{F}" surface of a half-filled spinon on the triangular lattice. Experiments to detect the spinon Fermi surface phase in 1T-TaS_{2} are discussed.

Theory of a quantum spin liquid in the hydrogen-intercalated honeycomb iridate H3LiIr2O6

We propose a theoretical model for a gapless spin liquid phase that may have been observed in a recent experiment on $\mathrm{H_3Li Ir_2 O_6}$. Despite the insulating and non-magnetic nature of the material, the specific heat coefficient $C/T \sim 1/\sqrt{T}$ in zero magnetic field and $C/T \sim T/ B^{3/2}$ with finite magnetic field $B$ have been observed. In addition, the NMR relaxation rate shows $1/(T_1T) \sim (C/T)^2$. Motivated by the fact that the interlayer/in-plane lattice parameters are reduced/elongated by the hydrogen-intercalation of the parent compound $\mathrm{Li_2 Ir O_3}$, we consider four layers of the Kitaev honeycomb lattice model with additional interlayer exchange interactions. It is shown that the resulting spin liquid excitations reside mostly in the top and bottom layers of such a layered structure and possess a quartic dispersion. In an applied magnetic field, each quartic mode is split into four Majorana cones with the velocity $v \sim B^{3/4}$. We suggest that the spin liquid phase in these "defect" layers, placed between different stacking patterns of the honeycomb layers, can explain the major phenomenology of the experiment, which can be taken as evidence that the Kitaev interaction plays the primary role in the formation of a quantum spin liquid in this material.

Gapless chiral spin liquid from coupled chains on the kagome lattice

Using a perturbative renormalization group approach, we show that the extended (J_1-J_2-J_dJ1−J2−Jd) Heisenberg model on the kagome lattice with a staggered chiral interaction (J_\chiJχ) can exhibit a gapless chiral quantum spin liquid phase. Within a coupled-chain construction, this phase can be understood as a chiral sliding Luttinger liquid with algebraic decay of spin correlations along the chain directions. We calculate the low-energy properties of this gapless chiral spin liquid using the effective field theory and show that they are compatible with the predictions from parton mean-field theories with symmetry-protected line Fermi surfaces. These results may be relevant to the state observed in the kapellasite material.

Ground-state phase diagram of the triangular lattice Hubbard model by the density-matrix renormalization group method

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 Uc1 = 7.55 t and Uc2 = 9.65 t (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 < U < Uc1 implies that this phase corresponds to a metallic phase. The local spin density induced by an applied pinning magnetic field for U > Uc2 exhibits a three sublattice feature, which is compatible with the Neel ordered state realized in the strong coupling limit. For Uc1 < U < Uc2, a response to the applied pinning magnetic field is comparable to that in the metallic phase but a relatively large spin correlation length is found with neither valence bond nor chiral magnetic order, suggesting a paramagnetic nature which resembles gapless spin liquid. The calculation also finds that the pair- ing correlation function monotonically decreases with increasing U and thus the superconductivity is unlikely in the intermediate phase.

Field-Induced Instability of a Gapless Spin Liquid with a Spinon Fermi Surface.

The ground state of the quantum kagome antiferromagnet Zn-brochantite, ZnCu_{3}(OH)_{6}SO_{4}, which is one of only a few known spin-liquid (SL) realizations in two or three dimensions, has been described as a gapless SL with a spinon Fermi surface. Employing nuclear magnetic resonance in a broad magnetic-field range down to millikelvin temperatures, we show that in applied magnetic fields this enigmatic state is intrinsically unstable against a SL with a full or a partial gap. A similar instability of the gapless Fermi-surface SL was previously encountered in an organic triangular-lattice antiferromagnet, suggesting a common destabilization mechanism that most likely arises from spinon pairing. A salient property of this instability is that an infinitesimal field suffices to induce it, as predicted theoretically for some other types of gapless SLs.

Generalized electromagnetism of subdimensional particles: A spin liquid story

It has recently been shown that there exists a class of stable gapless spin liquids in 3+1 dimensions described by higher rank tensor U(1) gauge fields, giving rise to an emergent tensor electromagnetism. The tensor gauge field of these theories couples naturally to subdimensional particles (such as fractons), which are restricted by gauge invariance to move only along lower-dimensional subspaces of the system. We here work out some of the basic generalized electromagnetic properties of subdimensional particles coupled to tensor electromagnetism, such as generalized electrostatic fields, potential formulations, Lorentz forces, Maxwell equations, and Biot-Savart laws. Some concepts from conventional electromagnetism will carry over directly, while others require significant modification.

Thermodynamics of the quantum spin liquid state of the single-component dimer Mott system κ−H3(Cat−EDT−TTF)2

Thermodynamic properties of the proton-mediated single component dimer-Mott insulator of $\ensuremath{\kappa}\text{\ensuremath{-}}{\mathrm{H}}_{3}{(\mathrm{Cat}\text{\ensuremath{-}}\mathrm{EDT}\text{\ensuremath{-}}\mathrm{TTF})}_{2}$, which has a two-dimensional triangle lattice structure of $S=1/2$ spins, are reported. The extraordinary large electronic heat capacity coefficient $\ensuremath{\gamma}\phantom{\rule{0.16em}{0ex}}=\phantom{\rule{0.16em}{0ex}}58.8\phantom{\rule{0.28em}{0ex}}\mathrm{mJ}\phantom{\rule{0.16em}{0ex}}{\mathrm{K}}^{\ensuremath{-}2}\phantom{\rule{0.16em}{0ex}}\mathrm{mo}{\mathrm{l}}^{\ensuremath{-}1}$ observed by the low-temperature heat capacity measurements up to 6 T suggests the formation of the gapless spin liquid ground state. Although the magnetic interaction $J/{k}_{\mathrm{B}}$ is quite different from those of other dimer-Mott spin liquids, the thermodynamic feature scales well with the Wilson ratio of 1.4\char21{}1.6. The heat capacity measurements also detected that the deuteration of the proton-linkage changes the ground state to the nonmagnetic one with almost vanishing \ensuremath{\gamma}. Using the data of the deuterated compound, the accurate temperature dependence of the magnetic heat capacity reflecting on the low-energy excitations from the gapless spin liquids ground state is discussed.

Gapless Spin-Liquid Ground State in the S=1/2 Kagome Antiferromagnet.

The defining problem in frustrated quantum magnetism, the ground state of the nearest-neighbor S=1/2 antiferromagnetic Heisenberg model on the kagome lattice, has defied all theoretical and numerical methods employed to date. We apply the formalism of tensor-network states, specifically the method of projected entangled simplex states, which combines infinite system size with a correct accounting for multipartite entanglement. By studying the ground-state energy, the finite magnetic order appearing at finite tensor bond dimensions, and the effects of a next-nearest-neighbor coupling, we demonstrate that the ground state is a gapless spin liquid. We discuss the comparison with other numerical studies and the physical interpretation of this result.

Subdimensional particle structure of higher rank U(1) spin liquids

Spin liquids are conventionally described by gauge theories with a vector gauge field. However, there exists a wider class of spin liquids with higher rank tensors as the gauge variable. In this work, we focus on (3+1)-dimensional spin liquids described by U(1) symmetric tensor gauge theories, which have recently been shown to be stable gapless spin liquids. We investigate the particle structure of these tensor gauge theories and find that they have deep connections with the "fracton" models recently discovered by Vijay, Haah, and Fu. Tensor gauge theories have more conservation laws than the simple charge conservation law of rank 1 theories. These conservation laws place severe restrictions on the motion of particles. Particles in some models are fully immobile (fractons), while other models have particles restricted to motion along lower-dimensional subspaces.

Interplay of nonsymmorphic symmetry and spin-orbit coupling in hyperkagome spin liquids: Applications to Na4Ir3O8

Na$_4$Ir$_3$O$_8$ provides a material platform to study three-dimensional quantum spin liquids in the geometrically frustrated hyperkagome lattice of Ir$^{4+}$ ions. In this work, we consider quantum spin liquids on hyperkagome lattice for generic spin models, focusing on the effects of anisotropic spin interactions. In particular, we classify possible $\mathbb{Z}_2$ and $U(1)$ spin liquid states, following the projective symmetry group analysis in the slave-fermion representation. There are only three distinct $\mathbb{Z}_2$ spin liquids, together with 2 different $U(1)$ spin liquids. The non-symmorphic space group symmetry of hyperkagome lattice plays a vital role in simplifying the classification, forbidding "$\pi$-flux" or "staggered-flux" phases in contrast to symmorphic space groups. We further prove that both $U(1)$ states and one $Z_2$ state among all 3 are symmetry-protected gapless spin liquids, robust against any symmetry-preserving perturbations. Motivated by the "spin-freezing" behavior recently observed in Na$_4$Ir$_3$O$_8$ at low temperatures, we further investigate the nearest-neighbor spin model with dominant Heisenberg interaction subject to all possible anisotropic perturbations from spin-orbit couplings. We found a $U(1)$ spin liquid ground state with spinon fermi surfaces is energetically favored over $Z_2$ states. Among all spin-orbit coupling terms, we show that only Dzyaloshinskii-Moriya (DM) interaction can induce spin anisotropy in the ground state when perturbing from the isotropic Heisenberg limit. Our work paves the way for a systematic study of quantum spin liquids in various materials with a hyperkagome crystal structure.

Coherent hole propagation in an exactly solvable gapless spin liquid

We examine the dynamics of a single hole in the gapless phase of the Kitaev honeycomb model, focusing on the slow-hole regime where the bare hopping amplitude $t$ is much less than the Kitaev exchange energy $J$. In this regime, the hole does not generate gapped flux excitations and is dressed only by the gapless fermion excitations. Investigating the single-hole spectral function, we find that the hole propagates coherently with a quasiparticle weight that is finite but approaches zero as $t/J \to 0$. This conclusion follows from two approximate treatments, which capture the same physics in complementary ways. Both treatments use the stationary limit as an exactly solvable starting point to study the spectral function approximately (i) by employing a variational approach in terms of a trial state that interpolates between the limits of a stationary hole and an infinitely fast hole and (ii) by considering a special point in the gapless phase that corresponds to a simplified one-dimensional problem.

Fate of Spinons at the Mott Point.

Gapless spin liquids have recently been observed in several frustrated Mott insulators, with elementary spin excitations-"spinons"-reminiscent of degenerate Fermi systems. However, their precise role at the Mott point, where charge fluctuations begin to proliferate, remains controversial and ill understood. Here we present the simplest theoretical framework that treats the dynamics of emergent spin and charge excitations on the same footing, providing a new physical picture of the Mott metal-to-insulator transition at half filing. We identify a generic orthogonality mechanism leading to strong damping of spinons, arising as soon as the Mott gap closes. Our results indicate that spinons should not play a significant role within the high-temperature quantum critical regime above the Mott point-in striking agreement with all available experiments.

Emergent quasi-one-dimensionality in a kagome magnet: A simple route to complexity

We study the ground state phase diagram of the quantum spin-$1/2$ Heisenberg model on the kagom\'{e} lattice with first- ($J_1 < 0$), second- ($J_2 < 0$), and third-neighbor interactions ($J_d > 0$) by means of analytical low-energy field theory and numerical density-matrix renormalization group (DMRG) studies. The results offer a consistent picture of the $J_d$-dominant regime in terms of three sets of spin chains weakly coupled by the ferromagnetic inter-chain interactions $J_{1,2}$. When either $J_1$ or $J_2$ is dominant, the model is found to support one of two cuboctohedral phases, cuboc1 and cuboc2. These cuboc states host non-coplanar long-ranged magnetic order and possess finite scalar spin chirality. However, in the compensated regime $J_1 \simeq J_2$, a valence bond crystal phase emerges between the two cuboc phases. We find excellent agreement between an analytical theory based on coupled spin chains and unbiased DMRG calculations, including at a very detailed level of comparison of the structure of the valence bond crystal state. To our knowledge, this is the first such comprehensive understanding of a highly frustrated two-dimensional (2d) quantum antiferromagnet. We find no evidence of either the one-dimensional (1d) gapless spin liquid or the chiral spin liquids, which were previously suggested by parton mean field theories.

Classification of gaplessZ2spin liquids in three-dimensional Kitaev models

Frustrated quantum magnets can harbor unconventional spin liquid ground states in which the elementary magnetic moments fractionalize into new emergent degrees of freedom. While the fractionalization of quantum numbers is one of the recurring themes in modern condensed matter physics, it often remains a challenge to devise a controlled analytical framework tracking this phenomenon. A notable exception is the exactly solvable Kitaev model, in which spin degrees of freedom fractionalize into Majorana fermions and a Z2 gauge field. Here we discuss the physics of fractionalization in three-dimensional Kitaev models and demonstrate that the itinerant Majorana fermions generically form a (semi)metal which, depending on the underlying lattice structure, exhibits Majorana Fermi surfaces, nodal lines or topologically protected Weyl nodes. We show that the nature of these Majorana metals can be deduced from an elementary symmetry analysis of the projective time-reversal and inversion symmetries for a given lattice. This allows us to comprehensively classify the gapless spin liquids of Kitaev models for the most elementary tricoordinated lattices in three dimensions. We further expand this classification by addressing the effects of time-reversal symmetry breaking and additional interactions.

Variational Monte Carlo study of a gapless spin liquid in the spin-12XXZ antiferromagnetic model on the kagome lattice

By using the variational Monte Carlo technique, we study the spin-$1/2$ XXZ antiferromagnetic model (with easy-plane anisotropy) on the kagome lattice. A class of Gutzwiller projected fermionic states with a spin Jastrow factor is considered to describe either spin liquids (with $U(1)$ or $Z_2$ symmetry) or magnetically ordered phases (with ${\bf q}=(0,0)$ or ${\bf q}=(4\pi/3,0)$). We find that the magnetic states are not stable in the thermodynamic limit. Moreover, there is no energy gain to break the gauge symmetry from $U(1)$ to $Z_2$ within the spin-liquid states, as previously found in the Heisenberg model. The best variational wave function is therefore the $U(1)$ Dirac state, supplemented by the spin Jastrow factor. Furthermore, a vanishing $S=2$ spin gap is obtained at the variational level, in the whole regime from the $XY$ to the Heisenberg model.

Gapless quantum spin liquid ground state in the two-dimensional spin-1/2 triangular antiferromagnet YbMgGaO4.

Quantum spin liquid (QSL) is a novel state of matter which refuses the conventional spin freezing even at 0 K. Experimentally searching for the structurally perfect candidates is a big challenge in condensed matter physics. Here we report the successful synthesis of a new spin-1/2 triangular antiferromagnet YbMgGaO4 with symmetry. The compound with an ideal two-dimensional and spatial isotropic magnetic triangular-lattice has no site-mixing magnetic defects and no antisymmetric Dzyaloshinsky-Moriya (DM) interactions. No spin freezing down to 60 mK (despite θw ~ −4 K), the power-law temperature dependence of heat capacity and nonzero susceptibility at low temperatures suggest that YbMgGaO4 is a promising gapless (≤|θw|/100) QSL candidate. The residual spin entropy, which is accurately determined with a non-magnetic reference LuMgGaO4, approaches zero (<0.6%). This indicates that the possible QSL ground state (GS) of the frustrated spin system has been experimentally achieved at the lowest measurement temperatures.

Spin-Peierls Instability of Three-Dimensional Spin Liquids with Majorana Fermi Surfaces.

Three-dimensional (3D) variants of the Kitaev model can harbor gapless spin liquids with a Majorana Fermi surface on certain tricoordinated lattice structures such as the recently introduced hyperoctagon lattice. Here, we investigate Fermi surface instabilities arising from additional spin exchange terms (such as a Heisenberg coupling) which introduce interactions between the emergent Majorana fermion degrees of freedom. We show that independent of the sign and structure of the interactions, the Majorana surface is always unstable. Generically, the system spontaneously doubles its unit cell at exponentially small temperatures and forms a spin liquid with line nodes. Depending on the microscopics, further symmetries of the system can be broken at this transition. These spin-Peierls instabilities of a 3D spin liquid are closely related to BCS instabilities of fermions.