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.
The
S=1/2 Kitaev model provides an exactly solvable setting in which spin fractionalization produces itinerant Majorana fermions coupled to an emergent
Z2 gauge field. In realistic Kitaev magnets, however, non-Kitaev interactions break integrability and make their finite-temperature dynamics difficult to calculate. We develop a finite-temperature parton mean-field theory based on the Abrikosov-fermion representation and incorporate the effects of exchange interactions through the random phase approximation to calculate the dynamical spin structure factor. For the pure Kitaev model, the method reproduces a strong, nearly dispersionless peak at a finite low energy together with a broad continuum at higher energies. It also reproduces the momentum dependence of the dynamical spin structure factor and its temperature evolution, including the downward shift of the low-energy spectral weight, in qualitative agreement with continuous-time quantum Monte Carlo calculations. Within the present theory, the contraction of the continuum is traced to the temperature-induced reduction of the bond amplitudes, which narrows the parton bands and compresses the two-parton phase space. For ferromagnetic Kitaev coupling, the energy of the low-energy peak is minimized at the
Γ point for a ferromagnetic Heisenberg interaction and maximized there for an antiferromagnetic interaction. In the ferromagnetic Kitaev-
Γ model, the static response develops a finite-momentum softening, indicating an instability toward the period-three magnetic order reported in numerical studies. For uniform-ferromagnetic states selected by a
[111] magnetic field, we find that a broad parton continuum coexists with magnonlike dispersive structures, producing spectral features consistent with a spinon-to-magnon crossover.