We perform a systematic and exact study of Majorana fermion dynamics in the Kitaev-Heisenberg-Γ model in a few finite-size clusters increasing in size up to twelve sites. We employ exact Jordan-Wigner transformations to evaluate certain measures of Majorana fermion correlation functions, which effectively capture matter and gauge Majorana fermion dynamics in different parameter regimes. An external magnetic field is shown to produce a profound effect on gauge fermion dynamics. Depending on certain non-zero choices of other non-Kitaev interactions, it can stabilise it to its non-interacting Kitaev limit. For all the parameter regimes, gauge fermions are seen to have slower dynamics, which could help build approximate decoupling schemes for appropriate mean-field theory. The probability of Majorana fermions returning to their original starting site shows that the Kitaev model in small clusters can be used as a test bed for the quantum speed limit.
We perform a systematic and exact study of Majorana fermion dynamics in the Kitaev-Heisenberg-Γ model in a few finite-size clusters increasing in size up to twelve sites. We employ exact Jordan-Wigner transformations to evaluate certain measures of Majorana fermion correlation functions, which effectively capture matter and gauge Majorana fermion dynamics in different parameter regimes. An external magnetic field is shown to produce a profound effect on gauge fermion dynamics. Depending on certain non-zero choices of other non-Kitaev interactions, it can stabilise it to its non-interacting Kitaev limit. For all the parameter regimes, gauge fermions are seen to have slower dynamics, which could help build approximate decoupling schemes for appropriate mean-field theory. The probability of Majorana fermions returning to their original starting site shows that the Kitaev model in small clusters can be used as a test bed for the quantum speed limit.
Recently, the Kitaev-Heisenberg-Γ system has been used to explore various aspects of Kitaev spin liquid physics. Here, we consider a few small clusters of up to twelve sites and study them in detail to unravel many interesting findings due to the competition between all possible signs and various magnitudes of these interactions under the influence of an external magnetic field. When Heisenberg interaction is taken anti-ferromagnetic, one obtains plateaus in correlation functions where, surprisingly, the exact groundstate reduces to the eigenstate of Heisenberg interaction as well. On the other hand, for ferromagnetic Heisenberg interaction, its competition with Kitaev interaction results in non-monotonicity in the correlation functions. We discuss, in detail, the competing effects on low energy spectrum, flux operator, magnetization, susceptibility, and specific heat. Finally, we discuss how our findings could be helpful to explain some of the recent experimental and theoretical findings in materials with Kitaev interactions.
Motivated by the appearance of Majorana fermions in a broad range of correlated and topological electronic systems, we develop a general method to compute the dynamical response of interacting Majorana fermions in the random-phase approximation (RPA). This can be applied self-consistently on top of Majorana mean-field theory backgrounds, thereby in particular providing a powerful tool to analyze behavior in the vicinity of (various heavily studied) exactly soluble models. Prime examples are quantum spin liquids (QSL) with emergent Majorana excitations, with the celebrated exact solution of Kitaev. We employ the RPA to study in considerable detail phase structure and dynamics of the extended Kitaev honeycomb K J Γ model, with and without an applied field. First, we benchmark our method with Kitaev's exactly soluble model, finding a remarkable agreement. The interactions between Majorana fermions even turn out to mimic the effect of local Z 2 flux excitations, which we explain analytically. Second, we show how small non-Kitaev couplings J and Γ induce Majorana bound states, resulting in sharp features in the dynamical structure factor in the presence of fractionalization: such “spinon excitons” naturally appear, and can coexist and interact with the broad Majorana continuum. Third, for increasing couplings or field, our theory predicts instabilities of the Kitaev QSL (KQSL) triggered by the condensation of the sharp modes. From the high-symmetry momenta of the condensation we can deduce which magnetically ordered phases surround the KQSL, in good agreement with previous finite-size numerics. We discuss implications for experiments and the broad range of applicability of our method to other QSL and Majorana systems.
We determine the phase diagrams of anisotropic Kitaev-Heisenberg models on the honeycomb lattice using parton mean-field theories based on different Majorana fermion representations of the S=1/2 spin operators. First, we use a two-dimensional Jordan-Wigner transformation (JWT) involving a semi-infinite snake string operator. To ensure that the fermionized Hamiltonian remains local, we consider the limit of extreme Ising exchange anisotropy in the Heisenberg sector. Second, we use the conventional Kitaev representation in terms of four Majorana fermions subject to local constraints, which we enforce through Lagrange multipliers. For both representations, we self-consistently decouple the interaction terms in the bond and magnetization channels and determine the phase diagrams as a function of the anisotropy of the Kitaev couplings and the relative strength of the Ising exchange. While both mean-field theories produce identical phase boundaries for the topological phase transition between the gapless and gapped Kitaev quantum spin liquids, the JWT fails to correctly describe the magnetic instability and finite-temperature behavior. Our results show that the magnetic phase transition is first order at low temperatures but becomes continuous above a certain temperature. At this energy scale, we also observe a finite-temperature crossover on the quantum-spin-liquid side, from a fractionalized paramagnet at low temperatures, in which gapped flux excitations are frozen out, to a conventional paramagnet at high temperatures. Published by the American Physical Society 2024
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.
We study the spin transport through the quantum spin liquid (QSL) by investigating the real-time and real-space dynamics of the Kitaev spin system with zigzag edges using the time-dependent Majorana mean-field theory. After the magnetic-field pulse is introduced to one of the edges, spin moments are excited in the opposite edge region although spin moments are never induced in the Kitaev QSL region. This unusual spin transport originates from the fact that the S=1/2 spins are fractionalized into the itinerant and localized Majorana fermions in the Kitaev system. Although both Majorana fermions are excited by the magnetic pulse, only the itinerant ones flow through the bulk regime without spin excitations, resulting in the spin transport in the Kitaev system despite the presence of a nonzero spin gap. We also demonstrate that this phenomenon can be observed in the system with small Heisenberg interactions using the exact diagonalization.
Emergent gauge theories take a prominent role in the description of quantum matter, supporting deconfined phases with topological order and fractionalized excitations. A common construction of \mathbb{Z}_2 ℤ 2 lattice gauge theories, first introduced by Wegner, involves Ising gauge spins placed on links and subject to a discrete \mathbb{Z}_2 ℤ 2 Gauss law constraint. As shown by Kitaev, \mathbb{Z}_2 ℤ 2 lattice gauge theories also emerge in the exact solution of certain spin systems with bond-dependent interactions. In this context, the \mathbb{Z}_2 ℤ 2 gauge field is constructed from Majorana fermions, with gauge constraints given by the parity of Majorana fermions on each site. In this work, we provide an explicit Jordan-Wigner transformation that maps between these two formulations on the square lattice, where the Kitaev-type gauge theory emerges as the exact solution of a spin-orbital (Kugel-Khomskii) Hamiltonian. We then apply our mapping to study local perturbations to the spin-orbital Hamiltonian, which correspond to anisotropic interactions between electric-field variables in the \mathbb{Z}_2 ℤ 2 gauge theory. These are shown to induce anisotropic confinement that is characterized by emergence of weakly-coupled one-dimensional spin chains. We study the nature of these phases and corresponding confinement transitions in both absence and presence of itinerant fermionic matter degrees of freedom. Finally, we discuss how our mapping can be applied to the Kitaev spin-1/2 model on the honeycomb lattice.