Kitaev Model Research Articles

Magnetic-field induced melting of long-range magnetic order akin to Kitaev insulators in the metallic compound Tb5Si3

There have been constant efforts to find ‘exotic’ quantum spin-liquid (QSL) materials. Some of the transition metal insulators dominated by the direction-dependent anisotropic exchange interaction (‘Kitaev model’ for honeycomb network of magnetic ions) are considered to be promising cases for the same. In such Kitaev insulators, QSL is achieved from the zero-field antiferromagnetic state by the application of magnetic-field, suppressing other exchange interactions responsible for magnetic order. Here, we show that the features attributable to long-range magnetic ordering of the intermetallic compound, Tb5Si3, (T N = 69 K), containing honey-comb network of Tb ions, are completely suppressed by a critical applied field, H cr, in heat-capacity and magnetization data, mimicking the behavior of Kitaev physics candidates. The neutron diffraction patterns as a function of H reveal that it is an incommensurate magnetic structure that gets suppressed, showing peaks arising from multiple wave vectors beyond H cr. Increasing magnetic entropy as a function of H with a peak in the magnetically ordered state is in support of some kind of magnetic disorder in a narrow field range after H cr. Such a high-field behavior for a metallic heavy rare-earth system to our knowledge has not been reported in the past and therefore is intriguing.

Z_{2} Spin Liquids in the Higher Spin-S Kitaev Honeycomb Model: An Exact Deconfined Z_{2} Gauge Structure in a Nonintegrable Model.

The higher spin Kitaev model prominently features the extensive locally conserved quantities the same as the spin-1/2 Kitaev honeycomb model, although it is not exactly solvable. It remains an open question regarding the physical meaning of these conserved quantities in the higher spin model. In this Letter, by introducing a Majorana parton construction for a general spin-S we uncover that these conserved quantities are exactly the Z_{2} gauge fluxes in the general spin-S model, including the case of spin-1/2. Particularly, we find an even-odd effect that the Z_{2} gauge charges are fermions in the half integer spin model, but are bosons in the integer spin model. We further prove that the fermionic Z_{2} gauge charges are always deconfined; hence, the half integer spin Kitaev model would have nontrivial spin liquid ground states regardless of interaction strengths in the Hamiltonian. The bosonic Z_{2} gauge charges of the integer spin model, on the other hand, could condense, leading to a trivial product state, and this is indeed the case at the anisotropic limit of the model.

Variational quantum simulation of the quantum critical regime

The quantum critical regime marks a zone in the phase diagram where quantum fluctuation around the critical point plays a significant role at finite temperatures. While it is of great physical interest, simulation of the quantum critical regime can be difficult on a classical computer due to its intrinsic complexity. Herein, we propose a variational approach, which minimizes the variational free energy, to simulate and locate the quantum critical regime on a quantum computer. The variational quantum algorithm adopts an ansatz by performing an unitary operator on a product of a single-qubit mixed state, in which the entropy can be analytically obtained from the initial state, and thus the free energy can be accessed conveniently. With numeral simulation, using the one-dimensional Kitaev model as a demonstration we show that the quantum critical regime can be identified by accurately evaluating the temperature crossover line. Moreover, the dependencies of both the correlation length and the phase coherence time with temperature are evaluated for the thermal states. Our work suggests a practical way as well as a first step for investigating quantum critical systems at finite temperatures on quantum devices with few qubits.

Nonlinear response of the Kitaev honeycomb lattice model in a weak magnetic field

We investigate the nonlinear response of the Kitaev honeycomb lattice model in a weak magnetic field using the theory of two-dimensional (2D) coherent spectroscopy. We observe that at the isotropic point in the non-Abelian phase of this model, the nonlinear spectrum in the 2D frequency domain consists of sharp signals that originate from the flux excitations and Majorana bound states. Signatures of different flux excitations can be clearly observed in this spectrum, such that one can observe evidences of flux states with 4-adjacent, 2-nonadjacent, and 4-far-separated fluxes, which are not visible in linear response spectroscopy such as neutron scattering experiments. Moreover, in the Abelian phase we perceive that the spectrum in the frequency domain is composed of streak signals. These signals, as in the nonlinear response of the pure Kitaev model, represent a distinct signature of itinerant Majorana fermions. However, deep in the Abelian phase whenever a Kitaev exchange coupling is much stronger than the others, the streak signals are weakened and only single sharp spots are seen in the response, which resembles the dispersionless response of the conventional toric code.

Disorder effects in the Kitaev-Heisenberg model

We study the interplay of disorder and Heisenberg interactions in the Kitaev model on a honeycomb lattice. The effect of disorder on the transition between Kitaev spin liquid and magnetic ordered states as well as the stability of magnetic ordering is investigated. Using Lanczos exact diagonalization we discuss the consequences of two types of disorder: (i) random-coupling disorder and (ii) singular-coupling disorder. They exhibit qualitatively similar effects in the pure Kitaev-Heisenberg model without long-range interactions. The range of spin-liquid phases is reduced and the transition to magnetic ordered phases becomes more crossoverlike. Furthermore, the long-range zigzag and stripy orderings in the clean system are replaced by their three domains with different ordering direction. Especially in the crossover range the coexistence of magnetically ordered and Kitaev spin-liquid domains is possible. With increasing the disorder strength the area of domains becomes smaller and the system goes into a spin-glass state. However, the disorder effect is different in magnetically ordered phases caused by long-range interactions. The stability of such magnetic ordering is diminished by singular-coupling disorder and, accordingly, the range of the spin-liquid regime is extended. This mechanism may be relevant to materials like $\ensuremath{\alpha}\text{\ensuremath{-}}{\mathrm{RuCl}}_{3}$ and ${\mathrm{H}}_{3}{\mathrm{LiIr}}_{2}{\mathrm{O}}_{6}$ where the zigzag ground state is stabilized by weak long-range interactions. We also find that the flux gap closes at a critical disorder strength and vortices appears in the flux arrangement. Interestingly, the vortices tend to form kinds of commensurate ordering.

Free energy subadditivity for symmetric random Hamiltonians

We consider a random Hamiltonian H:Σ→R defined on a compact space Σ that admits a transitive action by a compact group G. When the law of H is G-invariant, we show its expected free energy relative to the unique G-invariant probability measure on Σ, which obeys a subadditivity property in the law of H itself. The bound is often tight for weak disorder and relates free energies at different temperatures when H is a Gaussian process. Many examples are discussed, including branching random walks, several spin glasses, random constraint satisfaction problems, and the random field Ising model. We also provide a generalization to quantum Hamiltonians with applications to the quantum Sherrington–Kirkpatrick and Sachdev–Ye–Kitaev models.

A Z2 invariant for chiral and particle–hole symmetric topological chains

We define a Z2-valued topological and gauge invariant associated with any one-dimensional, translation-invariant topological insulator that satisfies either particle–hole symmetry or chiral symmetry. The invariant can be computed from the Berry phase associated with a suitable basis of Bloch functions that is compatible with the symmetries. We compute the invariant in the Su–Schrieffer–Heeger model for chiral symmetric insulators and in the Kitaev model for particle–hole symmetric insulators. We show that in both cases, the Z2 invariant predicts the existence of zero-energy boundary states for the corresponding truncated models.

Generalized black hole entropy in two dimensions

The Bekenstein–Hawking entropy of a black hole is proportional to its horizon area, hence in [Formula: see text] spacetime dimensions it is constant because the horizon degenerates into two points. This fact is consistent with Einstein’s gravity becoming topological in two dimensions. In [Formula: see text] gravity, which is non-trivial even in [Formula: see text], we find that the entropy is constant, as for Bekenstein–Hawking. As shown in Europhys. Lett. 139(6) (2022) 69001, arXiv: 2208.10146, two-dimensional [Formula: see text] gravity is equivalent to Jackiw–Teitelboim gravity, in turn, equivalent to the Sachdev–Ye–Kitaev model where the entropy becomes constant in the large [Formula: see text] limit. Several recently proposed entropies are functions of the Bekenstein–Hawking entropy and become constant in [Formula: see text], but in two-dimensional dilaton gravity entropies are not always constant. We study general dilaton gravity and obtain arbitrary static black hole solutions for which the non-constant entropies depend on the mass, horizon radius, or Hawking temperature, and constitute new proposals for a generalized entropy.

Spectrum of a Gross-Neveu Yukawa model with flavor disorder in three dimensions

We show that a variant of the Gross-Neveu Yukawa model with disorder provides a real, nonsupersymmetric generalization of the Sachdev--Ye Kitaev (SYK) model to three dimensions. The model contains $M$ real scalar fields and $N$ Dirac (or Majorana) fermions, interacting via a Yukawa interaction with a local Gaussian random coupling in three dimensions. In the limit where $M$ and $N$ are both large, and the ratio $M/N$ is held fixed, the model defines a line of infrared fixed points parametrized by $M/N$, reducing to the Gross-Neveu vector model when $M/N=0$. When $M/N$ is nonzero, the model is dominated by melonic diagrams and gives rise to SYK-like physics. We compute the spectrum of single-trace operators in the theory, and find that it is real for all values of $M/N$.

Emergence of bond-dependent highly anisotropic magnetic interactions in Sr4RhO6 : A theoretical study

The quantum spin liquid states as a natural ground state of the Kitaev model has led to a quest for new materials candidates hosting Kitaev physics. Yet, there are very few material candidates in this category. Using a combination of $ab$ $initio$ and model Hamiltonian methods, we propose that Ruddlesden-Popper compound Sr$_4$RhO$_6$ belongs to this category. With a tight-binding model and exact diagonalization approach, we show that despite substantial trigonal-like distortion, the electronic and magnetic properties of Sr$_4$RhO$_6$ can be well described in terms of pseudo-spin = 1/2 states. Magnetic interactions among pseudo-spins, estimated using the second-order perturbation method are highly bond-dependent anisotropic in nature with two particularly noticeable features, antiferromagnetic Kitaev and Dzyaloshinskii-Moriya interactions. The gaped spin-wave spectra of Sr$_4$RhO$_6$ obtained with linear spin-wave theory is consistent with the underlying magnetic frustration. Additional analysis of the role of individual or a particular combination of magnetic interactions reveals that the spin-wave spectra of Sr$_4$RhO$_6$ is a combined effect of the highly anisotropic interactions and a relatively simpler minimal model may not be plausible in the current case. The crucial insights about coupling between the local structural features and magnetic properties of Sr$_4$RhO$_6$ obtained in this study may be helpful for future studies belonging to this class.

Topological phase in a Kitaev chain with spatially separated pairing processes

The dynamic balance between pair creation and annihilation processes takes a crucial role to the topological superconductivity in the Kitaev model. Here we study the effect of spatial separation of creation and annihilation terms, i.e., sources and drains of pair are arranged alternatively. In this regard, a non-Hermitian Hamiltonian is naturally considered, which may possess complex energy branches. However, when the Bardeen-Cooper-Schrieffer pair excitation is only considered, the spectrum in such a subspace is fully real. In particular, the Zak phases extracted from the pair wave function are quantized and therefore able to characterize the different phases regardless of the breaking of time-reversal symmetry. For open chain system, the corresponding Majorana lattice is investigated. We find that although there are complex modes, all the edge modes have zero energy and obey the bulk-boundary correspondence. This results in the Kramer-type degeneracy of both the real and complex levels as a signature of the topologically nontrivial phase.

Variational optimization of tensor-network states with the honeycomb-lattice corner transfer matrix

We develop a method of variational optimization of the infinite projected entangled pair states on the honeycomb lattice. The method is based on the automatic differentiation of the honeycomb-lattice corner transfer matrix renormalization group. We apply the approach to the antiferromagnetic Heisenberg spin-1/2 and ferromagnetic Kitaev models on the honeycomb lattice. The developed formalism gives quantitatively accurate results for the main physical observables and has a necessary potential for further extensions.

Honeycomb oxide heterostructure as a candidate host for a Kitaev quantum spin liquid

Kitaev quantum spin liquid, massively quantum entangled states, is so scarce in nature that searching for new candidate systems remains a great challenge. Honeycomb heterostructure could be a promising route to realize and utilize such an exotic quantum phase by providing additional controllability of Hamiltonian and device compatibility, respectively. Here, we provide epitaxial honeycomb oxide thin film Na3Co2SbO6, a candidate of Kitaev quantum spin liquid proposed recently. We found a spin glass and antiferromagnetic ground states depending on Na stoichiometry, signifying not only the importance of Na vacancy control but also strong frustration in Na3Co2SbO6. Despite its classical ground state, the field-dependent magnetic susceptibility shows remarkable scaling collapse with a single critical exponent, which can be interpreted as evidence of quantum criticality. Its electronic ground state and derived spin Hamiltonian from spectroscopies are consistent with the predicted Kitaev model. Our work provides a unique route to the realization and utilization of Kitaev quantum spin liquid.

Field-induced Kitaev multipolar liquid in spin-orbit coupled d2 honeycomb Mott insulators

The Kitaev model, characterized by bond-dependent Ising spin interactions among spin-orbit entangled dipole moments in Mott insulators, offered a new approach to quantum spin liquids. Motivated by another type of bond-dependent interaction among quadrupole moments in $5{d}^{2}$ Mott insulators, we provide a microscopic route to uncover the Kitaev multipolar liquid, featuring fractionalized excitations out of non-Kramers doublets carrying multipole moments. The key ingredient is the magnetic field that allows for bond-anisotropic quadrupole-octupole interactions via mixing with the excited triplet states. The conditions to realize signatures of this phase in real materials are also discussed.

Connection between the Kitaev chain and the Gutzwiller-projected BCS model

We show that the topological Kitaev model can be placed in the context of phenomena associated with strongly correlated electron systems. This can be achieved implementing the Gutzwiller projection of the conventional s-wave BCS superconductor. Within that approach, neither spin–orbit coupling nor time-reversal symmetry breaking are needed to produce such an effect.

Geometrical frustration versus Kitaev interactions in BaCo2(AsO4)2.

Recently, Co-based honeycomb magnets have been proposed as promising candidate materials to host the Kitaev spin liquid (KSL) state. One of the front-runners is BaCo2(AsO4)2 (BCAO), where it was suggested that the exchange processes between Co2+ ions via the surrounding edge-sharing oxygen octahedra could give rise to bond-dependent Kitaev interactions. In this work, we present and analyze a comprehensive inelastic neutron scattering (INS) study of BCAO with fields in the honeycomb plane. Combining the constraints from the magnon excitations in the high-field polarized state and the inelastic spin structure factor measured in zero magnetic field, we examine two leading theoretical models: the Kitaev-type [Formula: see text] model and the XXZ[Formula: see text]model. We show that the existing experimental data can be consistently accounted for by the XXZ[Formula: see text]model but not by the [Formula: see text] model, and we discuss the implications of these results for the realization of a spin liquid phase in BCAO and more generally for the realization of the Kitaev model in cobaltates.

Thermal Hall conductivity with sign change in the Heisenberg–Kitaev kagome magnet

The Heisenberg–Kitaev (HK) model on various lattices has attracted a lot of attention because it may lead to exotic states such as quantum spin liquid and topological orders. The rare-earth-based kagome lattice (KL) compounds and have q = 0, 120° order and canted ferromagnetic (CFM) order, respectively. Interestingly, the HK model on the KL has the same ground state long-range orders. In the theoretical phase diagram, the CFM phase resides in a continuous parameter region and there is no phase change across special parameter points, such as the Kitaev ferromagnetic (KFM) point, the ferromagnetic (FM) point and its dual FM point. However, a ground state property cannot distinguish a system with or without topological nontrivial excitations and related phase transitions. Here, we study the topological magnon excitations and related thermal Hall conductivity in the HK model on the KL with CFM order. The CFM phase can be divided into two regions related by the Klein duality, with the self dual KFM point as their boundary. We find that the scalar spin chirality, which is intrinsic in the CFM order, changes sign across the KFM point. This leads to the opposite Chern numbers of corresponding magnon bands in the two regions, and also the sign change of the magnon thermal Hall conductivity.

Unifying Kitaev Magnets, Kagomé Dimer Models, and Ruby Rydberg Spin Liquids

The exploration of quantum spin liquids (QSLs) has been guided by different approaches including the resonating valence bond (RVB) picture, deconfined lattice gauge theories and the Kitaev model. More recently, a spin liquid ground state was numerically established on the ruby lattice, inspired by the Rydberg blockade mechanism. Here we unify these varied approaches in a single parent Hamiltonian, in which local fluctuations of anyons stabilize deconfinement. The parent Hamiltonian is defined on kagom\'e triangles -- each hosting four RVB-like states -- and includes only Ising interactions and single-site transverse fields. In the weak-field limit, the ruby spin liquid and exactly soluble kagom\'e dimer models are recovered, while the strong-field limit reduces to the Kitaev honeycomb model, thereby unifying three seemingly different approaches to QSLs. We similarly obtain the chiral Yao-Kivelson model, honeycomb toric code and a new spin-1 quadrupolar Kitaev model. The last is shown to be in a QSL phase by a non-local mapping to the kagom\'e Ising antiferromagnet. We demonstrate various applications of our framework, including (a) an adiabatic deformation of the ruby lattice model to the exactly soluble kagom\'e dimer model, conclusively establishing the QSL phase in the former; and (b) demystifying the dynamical protocol for measuring off-diagonal strings in the Rydberg implementation of the ruby lattice spin liquid. More generally, we find an intimate connection between Kitaev couplings and the repulsive interactions used for emergent dimer models. For instance, we show how a spin-1/2 XXZ model on the ruby lattice encodes a Kitaev honeycomb model, providing a new route toward realizing the latter in cold-atom or solid-state systems.

Majorana nanowires for topological quantum computation

Majorana bound states are quasiparticle excitations localized at the boundaries of a topologically nontrivial superconductor. They are zero-energy, charge-neutral, particle–hole symmetric, and spatially-separated end modes which are topologically protected by the particle–hole symmetry of the superconducting state. Due to their topological nature, they are robust against local perturbations and, in an ideal environment, free from decoherence. Furthermore, unlike ordinary fermions and bosons, the adiabatic exchange of Majorana modes is noncommutative, i.e., the outcome of exchanging two or more Majorana modes depends on the order in which exchanges are performed. These properties make them ideal candidates for the realization of topological quantum computers. In this tutorial, I will present a pedagogical review of 1D topological superconductors and Majorana modes in quantum nanowires. I will give an overview of the Kitaev model and the more realistic Oreg–Lutchyn model, discuss the experimental signatures of Majorana modes, and highlight their relevance in the field of topological quantum computation. This tutorial may serve as a pedagogical and relatively self-contained introduction for graduate students and researchers new to the field, as well as an overview of the current state-of-the-art of the field and a reference guide to specialists.

Ground-state phase diagram of spin- S Kitaev-Heisenberg models

The Kitaev model, whose ground state is a quantum spin liquid (QSL), was originally conceived for spin $S=1/2$ moments on a honeycomb lattice. In recent years, the model has been extended to higher $S$ from both theoretical and experimental interests, but the stability of the QSL ground state has not been systematically clarified for general $S$, especially in the presence of other additional interactions, which inevitably exist in candidate materials. Here we study the spin-$S$ Kitaev-Heisenberg models by using an extension of the pseudofermion functional renormalization group method to general $S$. We show that, similar to the $S=1/2$ case, the phase diagram for higher $S$ contains the QSL phases in the vicinities of the pristine ferromagnetic and antiferromagnetic Kitaev models, in addition to four magnetically ordered phases. We find, however, that the QSL phases shrink rapidly with increasing $S$, becoming vanishingly narrow for $S\geq 2$, whereas the phase boundaries between the ordered phases remain almost intact. Our results provide a reference for the search of higher-$S$ Kitaev materials.