The Kitaev honeycomb model realizes a quantum spin liquid where spin degrees of freedom are fractionalized into itinerant Majorana fermions and localized vortices. The entanglement among these particles holds promise for executing fault-tolerant topological quantum computation. While several materials have been proposed to host a quantum spin liquid phase, concrete evidence confirming its existence remains elusive. To aid in identifying this phase, we propose an electrical test using the charge profile, offering a signature for the Z2 vortices in Kitaev materials. This study explores the electric charge response in Kitaev materials, considering a generalized spin model including Heisenberg, off-diagonal exchange interactions, and magnetic Zeeman field. Majorana mean-field theory is employed to analyze spin correlations around vortices, revealing spatial anisotropy in the charge distribution, quantified by the electric quadrupole moment. This mean-field approach remains exact in the pure Kitaev limit and captures the fundamental physics across the Kitaev spin liquid phase. The quadrupole-quadrupole interaction between distant vortices is shown to be either repulsive or attractive, depending on the parameters. We predict that electrically biased scanning probe tips can help create vortices at predetermined locations. The findings of our study pave the way for electrically manipulating Ising anyons in Kitaev spin liquids, offering a promising path for technological applications of Kitaev materials.
Abstract We formulate a Majorana mean-field theory for the extended J-K-Gamma Kitaev model in a magnetic Zeeman field of arbitrary direction, and apply it for studying spatially inhomogeneous states harboring vortices. This mean-field theory is exact in the pure Kitaev limit and captures the essential physics throughout the Kitaev spin liquid phase. We determine the charge profile around vortices and the corresponding quadrupole tensor. The quadrupole-quadrupole interaction between distant vortices is shown to be either repulsive or attractive, depending on parameters. We predict that electrically biased scanning probe tips enable the creation of vortices at preselected positions. Our results open new perspectives for the electric manipulation of Ising anyons in Kitaev spin liquids.
We show that spin-spin correlations in a non-Abelian Kitaev spin liquid are associated with a characteristic inhomogeneous charge density distribution in the vicinity of Z_{2} vortices. This density profile and the corresponding local electric fields are observable, e.g., by means of surface probe techniques. Conversely, by applying bias voltages to several probe tips, one can stabilize Ising anyons (Z_{2} vortices harboring a Majorana zero mode) at designated positions, where we predict a clear Majorana signature in energy absorption spectroscopy.
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
Highly entangled excitations such as Majorana fermions of Kitaev quantum spin liquids have been proposed to be utilized for future quantum science and technology, and a deeper understanding of such excitations has been strongly desired. Here we demonstrate that Majorana fermion's mass and associated topological quantum phase transitions in the Kitaev quantum spin liquids may be manipulated by using electric fields in sharp contrast to the common belief that an insulator is inert under weak electric fields due to charge energy gaps. Using general symmetry analysis with perturbation and exact diagonalization, we uncover the universal phase diagrams with electric and magnetic fields. We also provide distinctive experimental signatures to identify Kitaev quantum spin liquids with electric fields, especially in connection with the candidate materials such as $\ensuremath{\alpha}\text{\ensuremath{-}}{\mathrm{RuCl}}_{3}$.
The Kitaev model on a honeycomb lattice with bond-dependent Ising interactions offers an exactly solvable model of a quantum spin liquid (QSL) with fractionalized excitations: gapped
Z2 fluxes and gapless linearly dispersing Majorana fermions in the isotropic limit (
Kx=Ky=Kz). We explore the phase diagram along two axes, an external magnetic field
h, applied out of plane of the honeycomb, and anisotropic interactions
Kz larger than the other two (${K}_{x}={K}_{y}\ensuremath{\equiv}K$). For ${K}_{z}/K\ensuremath{\gg}2$ and
h=0, the matter Majorana fermions have the largest gap, and the system is described by a gapped
Z2 toric code. One of the central questions we address is whether the fractionalized excitations in the different phases have sharp signatures that can be detected in experiments. We show that while the response to single-spin excitations is broad, the spectral function corresponding to two-spin excitations across a bond has sharp signatures that can be attributed to specific anyons. In the toric code regime, the $\ensuremath{\epsilon}=e\ifmmode\times\else\texttimes\fi{}m$ fermion, formed from the bosonic Ising electric (
e) and magnetic (
m) charges, disperses along a specific one-dimensional direction that provides a fingerprint of fractionalization. At lower
Kz in the center of the Abelian phase, in a regime we dub the primordial fractionalized regime, the field generates a hybridization between the $\ensuremath{\epsilon}$ fermion and the Majorana matter fermion, resulting in a $\ensuremath{\psi}$ fermion which too has a distinct quasi-one-dimensional dispersion. All the other phases in the field-anisotropy plane are naturally obtained from this primordial soup. These highly constrained fractonlike dispersions can be observable by inelastic light and neutron scattering, thereby providing ``smoking gun'' signatures of fractionalization in the QSL phase. Our analysis is based on calculations of susceptibilities, topological entanglement entropy, and excitation dynamics, obtained using exact diagonalization and density matrix renormalization group, and supported by perturbation theory.
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.