Abstract Quantum spin liquids (QSLs) give rise to exotic emergent particles by weaving intricate entanglement patterns in the underlying electrons. Bipartite measures between subregions can detect the presence of anyons, but little is known about the full entanglement structure of QSLs. Here, we study the multiparty entanglement of QSLs via entanglement microscopy. We find that in contrast to conventional matter, the genuine multiparty entanglement (GME) between spins is absent in the smallest subregions, a phenomenon we call ``entanglement frustration.'' Instead, GME is more collective, and arises solely in loops. By exploiting exact results and large-scale numerics, we confirm these properties in various gapped and gapless QSLs realized in physically motivated Hamiltonians, as well as with string-net wavefunctions hosting abelian or non-abelian anyons. Our results shed new light on the phase diagram of Kitaev’s honeycomb model in a Zeeman field, and the Kagome Heisenberg model under various perturbations. Going beyond QSLs, we provide evidence that entanglement loops are a universal property of quantum gauge theories. This leads to a new understanding of fractionalization, and the means by which gauge bosons encode quantum information.
Abstract Quantum spin liquids (QSLs) give rise to exotic emergent particles by weaving intricate entanglement patterns in the underlying electrons. Bipartite measures between subregions can detect the presence of anyons, but little is known about the full entanglement structure of QSLs. Here, we study the multiparty entanglement of QSLs via entanglement microscopy. We find that in contrast to conventional matter, the genuine multiparty entanglement (GME) between spins is absent in the smallest subregions, a phenomenon we call ``entanglement frustration.'' Instead, GME is more collective, and arises solely in loops. By exploiting exact results and large-scale numerics, we confirm these properties in various gapped and gapless QSLs realized in physically motivated Hamiltonians, as well as with string-net wavefunctions hosting abelian or non-abelian anyons. Our results shed new light on the phase diagram of Kitaev’s honeycomb model in a Zeeman field, and the Kagome Heisenberg model under various perturbations. Going beyond QSLs, we provide evidence that entanglement loops are a universal property of quantum gauge theories. This leads to a new understanding of fractionalization, and the means by which gauge bosons encode quantum information.
Abstract Many-body entanglement properties of quantum spin liquids (QSLs), persisting at arbitrarily long distances, have been intensely explored over the past two decades, but mostly for QSLs viewed as closed quantum systems. However, in experiments and potential quantum computing applications, candidate materials for this exotic phase of quantum matter will always interact with a dissipative environment, such as the one generated by bosonic quasiparticles in solids at finite temperature. Here we investigate both the stability and spatial distribution of entanglement for the Kitaev model of QSL, which is made open by its sudden coupling to an infinite bosonic bath of Caldeira–Leggett type and then time-evolved in both Markovian and non-Markovian regimes. From the time-dependent density matrix of QSL subregions, we extract genuine multipartite negativity (GMN), quantum Fisher information, spin–spin correlators, and the expectation value (EV) of the Wilson loop operator. In particular, time dependence of GMN offers the most penetrating insights: ( i ) in the Markovian regime, it remains nonzero only in hexagonal loopy subregions of QSL (as also discovered very recently for closed QSLs), eventually vanishing on the same timescale on which the EV of the Wilson loop operator vanishes; ( ii ) in the non-Markovian regime with pronounced memory effects, surprisingly, GMN remains nonzero up to much higher temperatures. In addition, the non-Markovian dynamics generates emergent interactions between spins, thereby opening avenues for tailoring properties of QSL via engineering of dissipation.
Quantum spin liquids (QSLs) are long-range entangled phases of matter and a natural setting for exploring how quantum correlations generate quantum complexity. Motivated by the emergence of magic, or nonstabilizerness, as a diagnostic of many-body complexity beyond entanglement, we study magic of QSLs by computing the stabilizer Rényi entropy (SRE) of the Kitaev honeycomb model. We derive a correspondence between Pauli strings and products of itinerant Majorana operators in the model's free-fermion description; this enables the sampling of SRE using an optimized algorithm for Gaussian states in systems with thousands of spins. Our results show that magic is largest in the gapless phase, where subleading volume-law corrections indicate nonlocal contributions, while in the gapped phases, SRE decreases in quantitative agreement with a perturbative expansion that we develop in the anisotropic limit. At the topological phase transition, SRE exhibits universal scaling dictated by critical exponents.
Quantum spin liquids (QSLs) are exotic phases of matter that host fractionalized excitations. It is difficult for local probes to characterize QSL, whereas quantum entanglement can serve as a powerful diagnostic tool due to its nonlocality. The kagome antiferromagnetic Heisenberg model is one of the most studied and experimentally relevant models for QSL, but its solution remains under debate. Here, we perform a numerical Aharonov-Bohm experiment on this model and uncover universal features of the entanglement entropy. By means of the density matrix renormalization group, we reveal the entanglement signatures of emergent Dirac spinons, which are the fractionalized excitations of the QSL. This scheme provides qualitative insights into the nature of kagome QSL and can be used to study other quantum states of matter. As a concrete example, we also benchmark our methods on an interacting quantum critical point between a Dirac semimetal and a charge-ordered phase.
Quantum spin liquids (QSL) are phases of matter which are distinguished not by the symmetries they break, but rather by the patterns of entanglement within them. Although these entanglement properties have been widely discussed for ground states, the way in which QSL form at finite temperature remains an open question. Here we introduce a method of characterizing both the depth and spatial structure of entanglement, and use this to explore how patterns of entanglement form as temperature is reduced in two widely studied models of QSL, the Kitaev honeycomb model, and the spin-1/2 Heisenberg antiferromagnet on a Kagome lattice. These results enable us to evaluate both the temperature at which spins within the high-temperature paramagnet first become entangled, and the temperature at which the system first develops the structured, multipartite entanglement characteristic of its QSL ground state.
Frustrated magnetism provides a rich platform for realizing exotic quantum phases in condensed matter systems. Among these, quantum spin liquids (QSLs) are a particularly intriguing class of states that exhibit no magnetic ordering and preserve all symmetries down to zero temperature. Just as conventional magnets can be thought of as spin crystals,QSLs represent their liquid analogs. This thesis explores two broad types of QSLs in two-dimensional magnets, with a focus on their dynamic properties and the nature of transitions from conventional magnetic states into QSL phases. QSLs are topological phases that host fractionalized excitations known as spinons. In these states, mobile spinons are coupled to an emergent gauge field, forming a highly entangled phase of matter. The first part of this thesis focuses on Dirac spin liquids (DSLs), a class of QSLs where spinons exhibit a gapless Dirac dispersion. To study the dynamicsof these spin liquids beyond leading-order mean-field theory, we develop a novel large-N fermionic limit of gapless spinons. In this limit, all gauge excitations are suppressed, and the leading-order corrections instead arise from order parameter fluctuations, which may condense and drive transitions into symmetry-broken phases from the parent DSL. Incorporating random phase approximation (RPA) corrections to the spin susceptibility, we find that sharp, spin-wave–like modes are recovered in the excitation spectrum. This methodology successfully captures characteristic features of the dynamics near phase transitions and in weakly symmetry-broken phases proximate to a DSL. We apply this approach to microscopic models, including the J1-J2 XXZ model on the triangular lattice and the J1-J3 XXZ model on the honeycomb lattice. These models closely represent a wide range of frustrated materials, such as triangular and honeycomb cobaltates and triangular ytterbium compounds. Our theoretical predictions show excellent agreement with experimentally measured inelastic neutron scattering spectra from these materials. Gapless spinons can also sometimes become gapped through chiral symmetry breaking, leading to the formation of a chiral spin liquid (CSL). The CSL is an archetypal bosonic realization of a fractional quantum Hall (FQH) state at half-filling, and has been proposed as a time-reversal–breaking ground state in a variety of frustrated spin models. The second part of this thesis addresses the intriguing possibility of realizing a CSL through the "melting" of a chiral spin crystal—a conventional magnet with non-coplanar ordering. Such a transition offers a physically appealing analogue to the melting of a Wigner crystal into a fermionic FQH liquid. Focusing on the kagome lattice, where chiral interactions are knownto stabilize a CSL ground state, we study the effects of additional long-range Heisenberg interactions and demonstrate a phase transition into a non-coplanar octahedral magnetic order. To describe the universal class of transition between a CSL and a non-coplanar anti-ferromagnetic order, we then develop a novel Chern-Simons–matter field theory framework, constrained by symmetry principles and symmetry fractionalization patterns. This thesis thus provides new insights into the dynamical signatures and phase transitions of quantum spin liquids, offering a bridge between theoretical models and experimental observations in frustrated magnetic systems.
In this work, we investigate the Kitaev honeycomb model employing the recently developed Clifford Circuits Augmented Matrix Product States (CAMPS) method. While the model in the gapped phase is known to reduce to the toric code model - whose ground state is entirely constructible from Clifford circuits - we demonstrate that the very different gapless quantum spin liquid (QSL) phase can also be significantly disentangled with Clifford circuits. Specifically, CAMPS simulations reveal that approximately two-thirds of the entanglement entropy in the isotropic point arises from Clifford-circuit contributions, enabling dramatically more efficient computations compared to conventional matrix product state (MPS) methods. Crucially, this finding implies that the Kitaev QSL state retains significant Clifford-simulatable structure, even in the gapless phase with non-abelian anyon excitations when time reversal symmetry is broken. This property not only enhances classical simulation efficiency significantly but also suggests substantial resource reduction for preparing such states on quantum devices. As an application, we leverage CAMPS to study the Kitaev-Heisenberg model and determine the most accurate phase boundary between the anti-ferromagnetic phase and the Kitaev QSL phase in the model. Our results highlight how Clifford circuits can effectively disentangle the intricate entanglement of Kitaev QSLs, opening avenues for efficiently simulating related and similar strongly correlated models.