Schwinger boson theory for S = 1 Kitaev quantum spin liquids

The Kitaev model is an exactly solvable model that exhibits a quantum spin liquid as its ground state. While this model was originally proposed as an S = 1 / 2 spin model on a honeycomb lattice, extensions to higher-spin systems have recently attracted attention. In contrast to the S = 1 / 2 case, such higher- S models are not exactly solvable and remain poorly understood, particularly for spin excitations at finite temperatures. Here, we focus on the S = 1 Kitaev model, which has been proposed to host bosonic quasiparticles. We investigate this model using Schwinger boson mean-field theory, where bosonic spinons are introduced as fractional quasiparticles by extending bond operators to address anisotropic spin interactions. We determine the mean-field parameters that realize a quantum spin liquid in both ferromagnetic and antiferromagnetic Kitaev models. Based on this mean-field ansatz, we calculate the dynamical and equal-time spin structure factors. We find that, when one uses the conventional decoupling scheme based on Wick decoupling with respect to spinons to calculate spin correlations, the resultant spin structure factors exhibit a momentum dependence that is not consistent with the sign structure expected from the exchange interaction: They possess a strong spectral weight indicating ferromagnetic (antiferromagnetic) correlations in the antiferromagnetic (ferromagnetic) Kitaev model. To resolve this issue, we propose an alternative scheme for evaluating the spin correlations, which is based on decoupling with respect to the bond operators. We demonstrate that, in our scheme, this discrepancy is removed, and the momentum dependence of the spin structure factors is consistent with the sign of the exchange constant. We also calculate the temperature evolution of the dynamical spin structure factor and find that the continuum observed at zero temperature splits into two distinct structures as the temperature increases, which can be understood in terms of the bandwidth narrowing of spinons. Finally, we clarify the origin of why the two distinct decoupling schemes result in different momentum dependences of the spin structure factors and discuss their relationship to results obtained in previous studies.