The Kitaev honeycomb model (KHM) consists of spin-$1/2$ particles on a honeycomb lattice with direction-dependent Ising-like interactions. It can alternatively be described in terms of non-interacting Majorana fermions, can be solved exactly, and has a quantum spin-liquid ground state. Open boundaries then host Majorana zero modes (MZMs) that are robust against some types of disorder. We analyze the fate of the MZMs when they couple to an environment via a Lindblad master equation. By computing the time evolution of the density matrix, we find that when decoherence occurs, the steady state is mostly the maximally mixed state. Among the few exceptions is a parameter regime that realizes the superconducting Kitaev chain model with periodic boundary conditions. We consistently observe a quantum Zeno effect in the density matrix as well as in the entropy and fidelity, while it is not found in the energy gap of some gapped spin liquids. We thus present a comprehensive overview over MZMs coupled to a spin bath that is relevant to proposals to detect MZMs of Kitaev layers on surfaces using scanning tunneling microscopy (STM).
The Kitaev honeycomb model (KHM) consists of spin-1/2 particles on a honeycomb lattice with direction-dependent Ising-like interactions. It can alternatively be described in terms of non-interacting Majorana fermions, can be solved exactly, and has a quantum spin-liquid ground state. Open boundaries then host Majorana zero modes (MZMs) that are robust against some types of disorder. We analyze the fate of the MZMs when they couple to an environment via a Lindblad master equation. By computing the time evolution of the density matrix, we find that when decoherence occurs, the steady state is mostly the maximally mixed state. Among the few exceptions is a parameter regime that realizes the superconducting Kitaev chain model with periodic boundary conditions. We consistently observe a quantum Zeno effect in the density matrix as well as in the entropy and fidelity, while it is not found in the energy gap of some gapped spin liquids. We thus present a comprehensive overview over MZMs coupled to a spin bath that is relevant to proposals to detect MZMs of Kitaev layers on surfaces using scanning tunneling microscopy (STM).
Abstract The S=3/2 Kitaev honeycomb model (KHM) is a quantum spin liquid (QSL) state coupled to a static Z2 gauge field. Employing an SO(6) Majorana representation of spin3/2’s, we find an exact representation of the conserved plaquette fluxes in terms of static Z2 gauge fields akin to the S=1/2 KHM which enables us to treat the remaining interacting matter fermion sector in a parton mean-field theory. We uncover a ground-state phase diagram consisting of gapped and gapless QSLs. Our parton description is in quantitative agreement with numerical simulations, and is furthermore corroborated by the addition of a [001] single ion anisotropy (SIA) which continuously connects the gapless Dirac QSL of our model with that of the S=1/2 KHM. In the presence of a weak [111] SIA, we discuss an emergent chiral QSL within a perturbation theory.
The S=23 Kitaev honeycomb model (KHM) is unique among the spin-S Kitaev models due to a massive ground-state quasidegeneracy that hampered previous numerical and analytical studies. In a recent work [Jin et al., Nat. Commun. 13, 3813 (2022)], we showed how an SO(6) Majorana parton mean-field theory of the S=23 isotropic KHM explains the anomalous features of this Kitaev spin liquid (KSL) in terms of an emergent low-energy Majorana flat band. Away from the isotropic limit, the S=23 KSL generally displays a quadrupolar order with gapped or gapless Majorana excitations, features that were quantitatively confirmed by density-matrix renormalization group simulations. In this paper, we explore the connection between the S=23 KHM with Kugel-Khomskii models and discover exactly soluble examples for the latter. We perform a symmetry analysis for the variational parton mean-field Ans\"atze in the spin and orbital basis for different quantum liquid phases of the S=23 KHM. Finally, we investigate a proposed time-reversal symmetry-breaking spin liquid induced by [111] single-ion anisotropy and elucidate its topological properties as well as experimental signatures, e.g., an unquantized thermal Hall response.
Here we revisit the thermodynamics of the Kitaev quantum spin liquid realized on the honeycomb lattice. We address two main questions: First, we investigate whether there are observable thermodynamic signatures of the topological Majorana boundary modes of the Kitaev honeycomb model. We argue that for the time-reversal invariant case the residual low-temperature entropy is the primary thermodynamic signature of these Majorana edge modes, and verify using large-scale Monte Carlo simulations that this residual entropy is present in the full Kitaev model. When time-reversal symmetry is broken, the Majorana edge modes are potentially observable in more direct thermodynamic measurements such as the specific heat, though only at temperatures well below the bulk gap. % Second, we study the energetics, and the corresponding thermodynamic signatures, of the flux excitations in the Kitaev model. Specifically, we study the flux interactions on both cylinder and torus geometries numerically, and quantify their impact on the thermodynamics of the Kitaev spin liquid by using a polynomial fit for the average flux energy as a function of flux density and extrapolating it to the thermodynamic limit. By comparing this model to Monte Carlo simulations, we find that flux interactions have a significant quantitative impact on the shape and the position of the low-temperature peak in the specific heat.
Kitaev's model of spins interacting on a honeycomb lattice describes a quantum spin liquid, where an emergent static Z2 gauge field is coupled to Majorana fermions. In the presence of an external magnetic field and for a range of interaction strengths, the system behaves as a gapped, non-Abelian quantum spin liquid. In this phase, the vortex excitations of the emergent Z2 gauge field have Majorana zero modes bound to them. Motivated by recent experimental progress in measuring and characterizing real materials that could exhibit spin-liquid behavior, we analytically calculate the dynamical spin structure factor in the non-Abelian phase of the Kitaev's honeycomb model. In particular, we treat the case of quenched disorder in the vortex configurations. Our calculations reveal a peak in the low-energy dynamical structure factor that is a signature of the spin-liquid behavior. We map the effective Hamiltonian to that of a chiral p-wave superconductor by using the Jordan-Wigner transformation. Subsequently, we analytically calculate the wave functions of the Majorana zero modes, the energy splitting for finite separation of the vortices and, finally, the dynamical structure factor in presence of quenched disorder.
We study the effects of the flux configurations on the emergent Majorana fermions in the S=21 Kitaev model on a honeycomb lattice, where quantum spins are fractionalized into itinerant Majorana fermions and localized fluxes. A quantum spin liquid appears as the ground state of the Kitaev model in the flux-free sector, which has intensively been investigated so far. In this flux sector, the Majorana fermion system has linear dispersions and shows power-law behavior in the Majorana correlations. On the other hand, periodically arranged flux configurations yield low-energy excitations in the Majorana fermion system, which are distinctly different from those in the flux-free state. We find that one of the periodically arranged flux states results in the gapped Majorana dispersion and the exponential decay in the Majorana correlations. The Kitaev system with another flux configuration exhibits a semi-Dirac-type dispersion, leading to the power-law decay with a smaller power than that in the flux-free sector along symmetry axes. We also examine the effect of the randomness in the flux configurations and clarify that the Majorana density of states is filled by increasing the flux density, and power-law decay in the Majorana correlations remains. The present results could be important to control the motion of Majorana fermions, which carries the spin excitations, in the Kitaev candidate materials.
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.