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.
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.
Disentangling transformations play a central role in the classical simulation of quantum many-body systems, yet their analytic structure and underlying mechanism remain largely unexplored. Here, we study the structure of the disentangler in the Haldane phase of spin-1 systems using generalized Clifford circuits. To this end, we extend the Clifford-circuit-augmented matrix product states (CAMPS)-based density-matrix renormalization group (DMRG) method to spin-1 systems. Within this framework, we find that the local disentanglers optimized for the Haldane phase implement the generalized Kramers--Wannier (KW) transformation, and we analytically verify its optimality for the Affleck--Kennedy--Lieb--Tasaki (AKLT) state. Beyond reducing entanglement, the KW transformation maps the Haldane phase to a phase with spontaneously broken Z2 symmetry. This mapping is distinct from the Kennedy--Tasaki transformation and provides a new unitary route from symmetry-protected topological order to symmetry breaking.
Clifford circuits can be utilized to disentangle quantum states with polynomial cost, thanks to the Gottesman-Knill theorem. Based on this idea, the Clifford Circuits Augmented Matrix Product States (CAMPS) method, which is a seamless integration of Clifford circuits within the density-matrix renormalization group algorithm, was proposed recently and was shown to be able to reduce entanglement in various quantum systems. In this work, we further explore the power of the CAMPS method in critical spin chains described by conformal field theories (CFTs) in the scaling limit. We find that the optimized disentanglers correspond to duality transformations, which significantly reduce the entanglement entropy in the ground state. For the critical quantum Ising spin chain governed by the Ising CFT with self-duality, the Clifford circuits found by CAMPS coincide with the duality transformation, i.e., the Kramer-Wannier self-duality in the critical Ising chain. It reduces the entanglement entropy by mapping the free conformal boundary condition to the fixed one. In the more general case of the XXZ chain, the CAMPS gives rise to a duality transformation mapping the model to the quantum Ashkin-Teller spin chain. Our results highlight the potential of the framework as a versatile tool for uncovering hidden dualities and simplifying the entanglement structure of critical quantum systems.
Determining the quantum-classical boundary between quantum circuits that can be efficiently simulated classically and those that cannot remains a fundamental question. One approach to classical simulation is to represent the output of a quantum circuit as a Clifford-augmented matrix product state (CAMPS), which, via a disentangling algorithm, decomposes the wave function into Clifford and MPS components, and from which Pauli expectation values can be computed in time polynomial in the MPS bond dimension. In this work, we develop an optimization-free disentangling (OFD) algorithm for Clifford circuits either doped with T gates or, equivalently, preceded by multiqubit gates of the form α I + β P . We give a simple and easily computed algebraic criterion that characterizes the individual quantum circuits for which OFD generates an efficient CAMPS—the bond dimension is exponential in the null space of a Galois binary field GF(2) matrix induced by a tableau of the twisted Pauli strings P . This significantly increases the number of circuits with rigorous polynomial-time classical simulations. We also give evidence that the typical N qubit random Clifford circuit doped with N uniformly distributed T gates of polylogarithmic depth or greater has a CAMPS with polynomial bond dimension. In addition, we compare OFD against disentangling by optimization. We further explore the representability of CAMPS for random Clifford circuits doped with more than N T gates. We also propose algorithms for sampling, probability, and amplitude estimation of bitstrings, and evaluation of entanglement Rényi entropy from CAMPS, which, though still having exponential complexity, are more efficient than standard MPS simulations. This work establishes a versatile framework for understanding classical simulatability of Clifford + T circuits and explores the interplay between quantum entanglement and quantum magic in quantum systems.
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.
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.
The main question we address is how to probe the fractionalized excitations of a quantum spin liquid (QSL), for example, in the Kitaev honeycomb model. By analyzing the energy spectrum and entanglement entropy, for antiferromagnetic couplings and a field along either [111] or [001], we find a gapless QSL phase sandwiched between the non-Abelian Kitaev QSL and polarized phases. Increasing the field strength towards the polarized limit destroys this intermediate QSL phase, resulting in a considerable reduction in the number of frequency modes and the emergence of a beating pattern in the local dynamical correlations, possibly observable in pump-probe experiments.