We give an exact form of the density matrix of the spin-1/2 Kitaev spin liquid represented in terms of spin operators and study the entanglement structures of the Kitaev honeycomb model within the spin framework. We show that the density matrix is naturally organized by equivalence classes of string operators associated with the underlying gauge structure of the model. With the explicit form of the density matrix, plus the exact Gauss law of the emergent gauge theory and the exact 1-form Wilson symmetry in the Kitaev model, we demonstrate the existence of the underlying symmetry-resolved block-diagonal structure of the reduced density matrix, which gives rise to the extensive degeneracy in the entanglement spectrum. The block-diagonal structure is then proven to be responsible for the separability of the entanglement entropy into the gauge and matter parts. Furthermore, we extend the formalism to subsystems with an odd number of lattice sites, revealing a relation between the entanglement spectrum and the fermion parity that is seldom mentioned in the literature.
While most of the features of the Kitaev honeycomb model are already clarified from the combined
Z2 plaquette and Majorana fermion picture, we show that retrieving the exact quantum many-body form gives a clue to fix the remaining issues on its spin liquid property. We provide an exact two-dimensional tensor network description of a full set of solutions from a random product state, by sequential quantum-number projection and Gutzwiller projection, both expressed as local matrix-product operators of dimension
D=2. The ground and excited states share the same bond dimension yet exhibit distinct Schmidt spectra, distinguishing their area-law and volume-law entanglement. We demonstrate that the topological entanglement entropy emerges {\it for all eigenstates only after} the Gutzwiller projection onto the physical spin Hilbert space, where the gauge constraint removes exactly half of the Schmidt states. This establishes that excited states are likewise topologically ordered. We further show that there exists a projection-induced collapse of the wave function into a single gauge sector in the Parton mean-field states that is known to give a good description of the spin liquids. The demonstrations on the
S=1/2 kagome and triangular Heisenberg models suggest that the present perspective is applied to a wider class of models hosting
Z2 spin liquids. Our finding visualizes a bulk-boundary correspondence with the Wilson-loop characterization of topological order.
A rigorous analysis is presented for the entanglement spectrum of quantum many-body states possessing a higher-form group-representation symmetry generated by topological Wilson loops, which is generally non-invertible. A general framework based on elementary algebraic topology and category theory is developed to determine the block structure of reduced density matrices for arbitrary bipartite manifolds on which the states are defined. Within this framework, we scrutinize the impact of topology on the entanglement structure for low-dimensional manifolds, including especially the torus, the Klein bottle, and lens spaces. By further incorporating gauge invariance, we refine our framework to determine the entanglement structure for topological gauge theories in arbitrary dimensions. In particular, in two dimensions, it is shown for the Kitaev quantum double model that not only the topological entanglement entropy can be reproduced, but also the Li-Haldane conjecture concerning the full entanglement spectrum holds exactly.
A spin fractionalizes into matter and gauge fermions in the Kitaev spin liquid on the honeycomb lattice, which is topologically ordered for certain parameters. To what extent does only one constituent of a fractionalized topological system signal the existence of topological order? To answer this, the authors focus on the reduced density matrix of the gauge sector on a cylinder after tracing out the matter degrees of freedom. The ensuing gauge entanglement Hamiltonian contains infinitely long-range correlations along the symmetry axis of the cylinder. Rather small correctly chosen gauge partitions can still account for the topological entanglement entropy in spite of long-range correlations in the gauge entanglement Hamiltonian.
Abstract In a quantum many-body system that possesses an additive conserved quantity, the entanglement entropy of a subsystem can be resolved into a sum of contributions from different sectors of the subsystem’s reduced density matrix, each sector corresponding to a possible value of the conserved quantity. Recent studies have discussed the basic properties of these symmetry-resolved contributions, and calculated them using conformal field theory and numerical methods. In this work we employ the generalized Fisher–Hartwig conjecture to obtain exact results for the characteristic function of the symmetry-resolved entanglement (‘flux-resolved entanglement’) for certain 1D spin chains, or, equivalently, the 1D fermionic tight binding and the Kitaev chain models. These results are true up to corrections of order where L is the subsystem size. We confirm that this calculation is in good agreement with numerical results. For the gapless tight binding chain we report an intriguing periodic structure of the characteristic functions, which nicely extends the structure predicted by conformal field theory. For the Kitaev chain in the topological phase we demonstrate the degeneracy between the even and odd fermion parity sectors of the entanglement spectrum due to virtual Majoranas at the entanglement cut. We also employ the Widom conjecture to obtain the leading behavior of the symmetry-resolved entanglement entropy in higher dimensions for an ungapped free Fermi gas in its ground state.
The Kitaev honeycomb model, which is exactly solvable by virtue of an extensive number of conserved quantities, supports a gapless quantum spin liquid phase as well as gapped descendants relevant for fault-tolerant quantum computation. We show that the anomalous edge modes of one-dimensional (1D) cluster-state-like symmetry-protected topological (SPT) phases provide natural building blocks for a variant of the Kitaev model that enjoys only a subextensive number of conserved quantities. The symmetry of our variant allows a single additional nearest-neighbor perturbation, corresponding to an anisotropic version of the $\mathrm{\ensuremath{\Gamma}}$ term studied in the context of Kitaev materials. We determine the phase diagram of the model using exact diagonalization. Additionally, we use the density matrix renormalization group to show that the underlying 1D SPT building blocks can emerge from a ladder Hamiltonian exhibiting only two-spin interactions supplemented by a Zeeman field. Our approach may provide a pathway toward realizing Kitaev honeycomb spin liquids in spin-orbit-coupled Mott insulators.
These lecture notes focus on the Kitaev honeycomb model, a two-dimensional quantumsystem with three distinct Ising-type interactions, depending on the direction of eachbond. This model has proven particularly important for describing quantum spin liquids,as well as for applications in topological quantum computing, due to its exotictopological properties. It constitutes the first exactly solvable theoretical model thatfully captures the behavior of a quantum spin liquid. We begin by presenting the geometryof the model, the Hamiltonian that describes the system, and Alexei Kitaev's effort toobtain an exact analytical solution, based on the decomposition of the Hilbert spaceinto flux sectors. We then present Ettore Majorana's observation that the Dirac equationalso admits real solutions, leading to the theoretical prediction of Majorana fermions-- particles that are identical to their antiparticles. We examine how the spatialseparation of Majorana operators leads to the rise of unbound Majorana fermions, andshowcase the exotic topological properties that manifest in physical systems such assuperconducting wires. Next, we represent the Hamiltonian in the basis of Majoranafermions, by expressing spin operators using four Majorana operators. The transitionfrom the extended space
L~ to the physical subspace
Lis achieved through a gauge transformation of the
Z2 group. We thenreformulate the Hamiltonian in the basis of complex bond and matter fermions, transitionto momentum space, and assume that the system lies in the zero-flux sector. Through theBogoliubov--de Gennes transformation, we arrive at the final diagonalized Hamiltonianin the quasiparticle basis. Finally, we represent the possible phases -- the gaplessphase (phase
B) and the gapped phases (phase
Aαi) -- compute thecorresponding energy spectrum, and present the complete phase diagram. By employing thecorrelation function, we show that all phases are characterized by strong quantumentanglement and spin fractionalization, classifying the system as a quantum spin liquid.