The Kitaev model, renowned for its exact solvability and potential to host non-Abelian anyons, remains a focal point in the study of quantum spin liquids and topological phases. While much of the existing literature has employed Majorana fermion techniques to analyze the model, particularly at zero temperature, its finite-temperature behavior has been less thoroughly explored via alternative methods. In this paper, we investigate the finite-temperature properties of the Kitaev model using the spin Green's function formalism. This approach provides a unified framework for computing key physical quantities, such as spin correlations, magnetic susceptibility, and the dynamical spin structure factor, offering valuable insights into the system's thermal dynamics. To solve the equation of motion for the spin Green's function, we truncate the hierarchy of multispin Green's functions using the Tyablikov decoupling approximation, which is particularly accurate at high temperatures. Our results show several similarities with Majorana-based numerical simulations, though notable differences emerge. Specifically, both static and dynamical spin-spin correlation functions capture not only
Z2 flux excitations but also simple spin-flip excitations, with the latter dominating the response. Additionally, without explicitly assuming fractionalization, our results for the spin susceptibility and spin relaxation rate suggest the presence of fermionic degrees of freedom at low temperatures. Unlike Majorana-based approaches, which rely on exact solvability and are inherently limited to the pure Kitaev model, the spin Green's function formalism is well-suited for studying more realistic systems. It can naturally incorporate the effects of magnetic fields, including linear-order terms, and non-Kitaev interactions, making it applicable to a broader range of materials and experimental conditions. By providing an alternative method for analyzing the finite-temperature properties of the Kitaev model, this study complements existing approaches and lays the groundwork for future investigations into real materials and extended models.
Temperature evolution of the spin correlation and excitation spectrum is investigated for the Kitaev model defined on a decorated honeycomb lattice by using the quantum Monte Carlo simulation in the Majorana fermion representation. The ground state of this quantum spin model is given by two kinds of chiral spin liquids: one is topologically trivial with Abelian anyon excitations, and the other is topologically nontrivial with non-Abelian anyon excitations. While lowering temperature, the model exhibits several crossovers in the paramagnetic state, which originate from the fractionalization of quantum spins into Majorana fermions, in addition to a phase transition associated with time reversal symmetry breaking. We show that the spin correlation develops around the crossover temperatures, whereas it shows a slight change at the critical temperature, as in other Kitaev-type models. We also calculate the excitation spectrum in terms of Majorana fermions, and find that the excitation gap in the non-Abelian phase is fragile against thermal fluctuations of the Z 2 fluxes, while that in the Abelian phase remains open.
Temperature evolution of the spin correlation and excitation spectrum is investigated for a Kitaev model defined on a decorated honeycomb lattice by using the quantum Monte Carlo simulation in the Majorana fermion representation. The ground state of this quantum spin model is given by two kinds of chiral spin liquids: one is topologically trivial with Abelian anyon excitations, and the other is topologically nontrivial with non-Abelian anyon excitations. While lowering temperature, the model exhibits several crossovers in the paramagnetic state, which originate from the fractionalization of quantum spins into Majorana fermions, in addition to a phase transition associated with time reversal symmetry breaking. We show that the spin correlation develops around the crossover temperatures, whereas it shows a slight change at the critical temperature, as in other Kitaev-type models. We also calculate the excitation spectrum in terms of Majorana fermions, and find that the excitation gap in the non-Abelian phase is fragile against thermal fluctuations of the
Z2 fluxes, while that in the Abelian phase remains open.
We study the simulation of the Kitaev spin model on quantum computers. In particular we focus on the models defined on the honeycomb, and square-octagon lattices. Using a fermionic language to describe these models reveals a region of the parameter space that is exactly solvable. We explore an ansatz that is capable of expressing the ground state in the exactly solvable region of the parameter space and extend it outside this region with good accuracy. Doing the calculation using fermions, while requiring the introduction of a non-local map from the fermionic Hilbert space to that of qubits, offers the potentially interesting application of realizing non-abelian anyons on quantum computers, and can also lead to a reduction in the number of qubits required by half.
The Kitaev honeycomb model is a paradigm of exactly solvable models, showing nontrivial physical properties such as topological quantum order, Abelian and non-Abelian anyons, and chirality. Its solution is one of the most beautiful examples of the interplay of different mathematical techniques in condensed matter physics. In this paper, we show how to derive a tensor network (TN) description of the eigenstates of this spin-1/2 model in the thermodynamic limit, and in particular for its ground state. In our setting, eigenstates are naturally encoded by an exact 3d TN structure made of fermionic unitary operators, corresponding to the unitary quantum circuit building up the many-body quantum state. In our derivation we review how the different ``solution ingredients'' of the Kitaev honeycomb model can be accounted for in the TN language, namely, Jordan-Wigner transformation, braidings of Majorana modes, fermionic Fourier transformation, and Bogoliubov transformation. The TN built in this way allows for a clear understanding of several properties of the model. In particular, we show how the fidelity diagram is straightforward both at zero temperature and at finite temperature in the vortex-free sector. We also show how the properties of two-point correlation functions follow easily. Finally, we also discuss the pros and cons of contracting of our 3d TN down to a 2d projected entangled pair state (PEPS) with finite bond dimension. The results in this paper can be extended to generalizations of the Kitaev model, e.g., to other lattices, spins, and dimensions.
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.