Floquet quantum error-correcting codes provide an operationally economical route to fault tolerance by dynamically generating stabilizer structures using only two-body Pauli measurements. But while it is well established that stabilizer codes in higher spatial dimensions gain additional levels of intrinsic robustness, higher-dimensional Floquet codes have hitherto been explored only in limited scope. Here we introduce a 3d generalization of a Floquet code whose instantaneous stabilizer group realizes a 3d fermionic toric code, while crucially preserving all three logical qubits throughout the entire measurement sequence. One central ingredient is the identification of a 3d lattice geometry that generalizes the features of the Kekulé lattice underlying the 2d Hastings-Haah code - specifically, a structure where deleting any one edge color yields a two-color subgraph that decomposes into short, closed loops rather than homologically nontrivial chains. This loop property avoids the collapse of logical information that plagues naive sequential two-color measurement schedules on many 3d lattices. Although, for our lattice geometry, a simple 3-round cycle that sequentially measures the three types of parity checks does not expose the full error syndrome set, we show that one can append a measurement sequence to extract the missing syndromes without disturbing the logical subspace. Beyond code design, 3d tricoordinated lattice geometries define a family of 3d monitored Kitaev models, in which random measurements of the non-commuting parity checks give rise to dynamically created entangled phases with nontrivial topology. In discussing the general structure of their underlying phase diagrams and, in particular, the existence of certain quantum critical points, we again make a connection to the general preservation of logical information in time-ordered Floquet protocols.
A prominent feature of quantum spin liquids is fractionalization of the spin degree of freedom. Fractionalized excitations have their own dynamics in different energy scales, and hence, affect finite-temperature (
T) properties in a peculiar manner even in the paramagnetic state harboring the quantum spin liquid state. We here present a comprehensive theoretical study of the spin dynamics in a wide
T range for the Kitaev model on a honeycomb lattice, whose ground state is a quantum spin liquid. In this model, the fractionalization occurs to break up quantum spins into itinerant matter fermions and localized
Z2 fluxes, which results in two crossovers at very different
T scales. Extending the previous study for the isotropic coupling case [J. Yoshitake, J. Nasu, and Y. Motome, Phys. Rev. Lett. 117, 157203 (2016)], we calculate the dynamical spin structure factor $S(\mathbf{q},\ensuremath{\omega})$, the NMR relaxation rate
1/T1, and the magnetic susceptibility $\ensuremath{\chi}$ while changing the anisotropy in the exchange coupling constants, by using the dynamical mean-field theory based on a Majorana fermion representation. We describe the details of the methodology including the continuous-time quantum Monte Carlo method for computing dynamical spin correlations and the maximum entropy method for analytic continuation. We confirm that the combined method provides accurate results in a wide
T range including the region where the spins are fractionalized. We find that also in the anisotropic cases the system exhibits peculiar behaviors below the high-
T crossover whose temperature is comparable to the average of the exchange constants: $S(\mathbf{q},\ensuremath{\omega})$ shows an inelastic response at the energy scale of the averaged exchange constant,
1/T1 continues to grow even though the equal-time spin correlations are saturated and almost
T independent, and $\ensuremath{\chi}$ deviates from the Curie-Weiss behavior. In particular, when the exchange interaction in one direction is stronger than the other two, the dynamical quantities exhibit qualitatively different
T dependences from the isotropic case at low
T, reflecting the opposite parity between the flux-free ground state and the flux-excited state, and a larger energy cost for flipping a spin in the strong interaction direction. On the other hand, when the exchange anisotropy is in the opposite way, the results are qualitatively similar to those in the isotropic case. All these behaviors manifest the spin fractionalization in the paramagnetic region. Among them, the dichotomy between the static and dynamical spin correlations is unusual behavior hardly seen in conventional magnets. We discuss the relation between the dichotomy and the spatial configuration of the
Z2 gauge fields. Our results could stimulate further experimental and theoretical analyses of candidate materials for the Kitaev quantum spin liquids.
Finite-temperature (
T) properties of a Kitaev model defined on a honeycomb lattice are investigated by a quantum Monte Carlo simulation, from the viewpoint of fractionalization of quantum
S=1/2 spins into two types of Majorana fermions, itinerant and localized. In this system, the entropy is released successively at two well-separated
T scales, as a clear indication of the thermal fractionalization. We show that the high-
T crossover, which is driven by itinerant Majorana fermions, is closely related with the development of nearest-neighbor spin correlations. On the other hand, the low-
T crossover originates in thermal fluctuations of fluxes composed of localized Majorana fermions, by which the spectrum of itinerant Majorana fermions is significantly disturbed. As a consequence, in the intermediate-
T range, the system exhibits coherent transport of Majorana fermions and
T-linear behavior in the specific heat, despite the Dirac semimetallic spectrum in the low-
T limit. We also show that the flux fluctuations affect the gapless-gapped phase boundary. Our results indicate that the fractionalization is experimentally observable in spin correlations, specific heat, and transport properties.
Arrays of Rydberg atoms constitute a highly tunable, strongly interacting venue for the pursuit of exotic states of matter. We develop a strategy for accessing a family of fractionalized phases known as quantum spin liquids in two-dimensional Rydberg arrays. We specifically use effective field theory methods to study arrays assembled from Rydberg chains tuned to an Ising phase transition that famously hosts emergent fermions propagating within each chain. This highly entangled starting point allows us to naturally access spin liquids familiar from Kitaev's honeycomb model---albeit from an entirely different framework. In particular, we argue that finite-range repulsive Rydberg interactions, which frustrate nearby symmetry-breaking orders, can enable coherent propagation of emergent fermions between the chains in which they were born. Delocalization of emergent fermions across the full two-dimensional Rydberg array yields a gapless
Z2 spin liquid with a single massless Dirac cone. Here, the Rydberg occupation numbers exhibit universal power-law correlations that provide a straightforward experimental diagnostic of this phase. We further show that explicitly breaking symmetries perturbs the gapless spin liquid into gapped, topologically ordered descendants: Breaking lattice symmetries generates toric-code topological order, whereas introducing Floquet-mediated chirality generates non-Abelian Ising topological order. In the toric-code phase, we analytically construct microscopic incarnations of non-Abelian defects, which can be created and transported by dynamically controlling the atom positions in the array. Our work suggests that appropriately tuned Rydberg arrays provide a cold-atoms counterpart of solid-state ``Kitaev materials'' and, more generally, it spotlights a different angle for pursuing experimental platforms for Abelian and non-Abelian fractionalization.
The formation of quantum spin liquids in frustrated magnets represents an exciting possibility due to the rather exotic features they harbor, including fractionalized excitations and emergent gauge fields. \nUnfortunately, they are notoriously difficult to study as there are often no good analytical methods available and quantum Monte Carlo simulations are hindered by the negative sign problem. \nThe Kitaev honeycomb model is a notable exception of a frustrated quantum model which is exactly solvable and which hosts a number of distinct quantum spin liquid ground states. \nAs such, it allows for a rare opportunity to study the physics of spin liquids with full analytical control. \n \nIn this thesis, we study the fractionalization of spin-1/2 moments into Majorana fermions and an emergent Z2 gauge field in a generalization of the Kitaev honeycomb model to a number of three-dimensional lattices. \nWhile the excitations of the gauge field are always gapped, the fermionic quasiparticles may exhibit a gapless dispersion, forming fully two-dimensional Fermi surfaces, symmetry protected nodal lines, or topological Weyl nodes. \nWe show that one can deduce rather general constraints on the possible gapless excitations by making use of an object called the projective symmetry group. \nIn doing so we provide a scheme for classifying the various gapless Kitaev spin liquids. \nA thorough analysis is carried out for a number of these spin liquids, primarily investigating the stability of the gapless modes and the novel features resulting from their sometimes non-trivial topology, as well as their effects on certain equal-time correlation functions.
The Kitaev model on a honeycomb lattice with bond-dependent Ising interactions offers an exactly solvable model of a quantum spin liquid (QSL) with fractionalized excitations: gapped
Z2 fluxes and gapless linearly dispersing Majorana fermions in the isotropic limit (
Kx=Ky=Kz). We explore the phase diagram along two axes, an external magnetic field
h, applied out of plane of the honeycomb, and anisotropic interactions
Kz larger than the other two (${K}_{x}={K}_{y}\ensuremath{\equiv}K$). For ${K}_{z}/K\ensuremath{\gg}2$ and
h=0, the matter Majorana fermions have the largest gap, and the system is described by a gapped
Z2 toric code. One of the central questions we address is whether the fractionalized excitations in the different phases have sharp signatures that can be detected in experiments. We show that while the response to single-spin excitations is broad, the spectral function corresponding to two-spin excitations across a bond has sharp signatures that can be attributed to specific anyons. In the toric code regime, the $\ensuremath{\epsilon}=e\ifmmode\times\else\texttimes\fi{}m$ fermion, formed from the bosonic Ising electric (
e) and magnetic (
m) charges, disperses along a specific one-dimensional direction that provides a fingerprint of fractionalization. At lower
Kz in the center of the Abelian phase, in a regime we dub the primordial fractionalized regime, the field generates a hybridization between the $\ensuremath{\epsilon}$ fermion and the Majorana matter fermion, resulting in a $\ensuremath{\psi}$ fermion which too has a distinct quasi-one-dimensional dispersion. All the other phases in the field-anisotropy plane are naturally obtained from this primordial soup. These highly constrained fractonlike dispersions can be observable by inelastic light and neutron scattering, thereby providing ``smoking gun'' signatures of fractionalization in the QSL phase. Our analysis is based on calculations of susceptibilities, topological entanglement entropy, and excitation dynamics, obtained using exact diagonalization and density matrix renormalization group, and supported by perturbation theory.