These lecture notes focus on the Kitaev honeycomb model, a two-dimensional quantumsystem with three distinct Ising-type interactions, depending on the direction of eachbond. This model has proven particularly important for describing quantum spin liquids,as well as for applications in topological quantum computing, due to its exotictopological properties. It constitutes the first exactly solvable theoretical model thatfully captures the behavior of a quantum spin liquid. We begin by presenting the geometryof the model, the Hamiltonian that describes the system, and Alexei Kitaev's effort toobtain an exact analytical solution, based on the decomposition of the Hilbert spaceinto flux sectors. We then present Ettore Majorana's observation that the Dirac equationalso admits real solutions, leading to the theoretical prediction of Majorana fermions-- particles that are identical to their antiparticles. We examine how the spatialseparation of Majorana operators leads to the rise of unbound Majorana fermions, andshowcase the exotic topological properties that manifest in physical systems such assuperconducting wires. Next, we represent the Hamiltonian in the basis of Majoranafermions, by expressing spin operators using four Majorana operators. The transitionfrom the extended space
L~ to the physical subspace
Lis achieved through a gauge transformation of the
Z2 group. We thenreformulate the Hamiltonian in the basis of complex bond and matter fermions, transitionto momentum space, and assume that the system lies in the zero-flux sector. Through theBogoliubov--de Gennes transformation, we arrive at the final diagonalized Hamiltonianin the quasiparticle basis. Finally, we represent the possible phases -- the gaplessphase (phase
B) and the gapped phases (phase
Aαi) -- compute thecorresponding energy spectrum, and present the complete phase diagram. By employing thecorrelation function, we show that all phases are characterized by strong quantumentanglement and spin fractionalization, classifying the system as a quantum spin liquid.
We investigate an S=3/2 quantum spin model on a two-dimensional honeycomb lattice that continuously interpolates between two paradigmatic quantum disordered states with distinct entanglement structures: the Kitaev quantum spin liquid and the Affleck-Kennedy-Lieb-Tasaki (AKLT) valence bond solid. Combining classical, semi-classical, and exact diagonalization approaches, we map out the ground-state phase diagram and elucidate the role of quantum fluctuations across the entire parameter range. While classical and semi-classical frameworks predict noncoplanar orders competing with a collinear Néel state, we find these phases to be fragile: once full quantum fluctuations are included, they melt into a quantum-entangled state characterized by suppressed spin correlations and enhanced entanglement entropy. Our findings highlight how competition between qualitatively different quantum disordered phases provides a fertile playground for unconventional phases emerging from their interplay and quantum fluctuations.
We study the ground-state phase diagram of the spin-
1/2 Kitaev-Heisenberg model on the bilayer honeycomb lattice with large-scale tensor network calculations based on the infinite projected entangled pair state technique as well as high-order series expansions. We find that beyond various magnetically ordered phases, including ferromagnetic, zigzag, antiferromagnetic (AFM) and stripy states, two extended quantum spin liquid phases arise in the proximity of the Kitaev limit. While these ordered phases also appear in the monolayer Kitaev-Heisenberg model, our results further show that a valence bond solid state emerges in a relatively narrow range of parameter space between the AFM and stripy phases, which can be adiabatically connected to isolated Heisenberg dimers. Our results highlight the importance of considering interlayer interactions on the emergence of novel quantum phases in the bilayer Kitaev materials.
We explore the phase diagram of the Kitaev-Heisenberg model with nearest neighbor interactions on the honeycomb lattice using the exact diagonalization of finite systems combined with the cluster mean field approximation, and supplemented by the insights from analytic approaches: the linear spin-wave and second-order perturbation theories. This study confirms that by varying the balance between the Heisenberg and Kitaev term, frustrated exchange interactions stabilize in this model either one of four phases with magnetic long range order: N\'eel phase, ferromagnetic phase, and two other phases with coexisting antiferromagnetic and ferromagnetic bonds, zigzag and stripy phase, or one of two distinct spin-liquid phases. Out of these latter disordered phases, the one with ferromagnetic Kitaev interactions has a substantially broader range of stability as the neighboring competing ordered phases, ferromagnetic and stripy, have very weak quantum fluctuations. Focusing on the quantum spin-liquid phases, we study spatial spin correlations and dynamic spin structure factor of the model by the exact diagonalization technique, and discuss the evolution of gapped low-energy spin response across the quantum phase transitions between the disordered spin liquid and phases with long range magnetic order.
We investigate the properties of the Kitaev honeycomb model with both site dilution and exchange randomness. Embarking on this work, we review disorder in some strongly correlated electron systems, including spin-½ and spin-1 Heisenberg antiferromagnetic chains, two dimensional Heisenberg antiferromagnets, the cuprates and graphene. We outline some aspects of resonating valence bond phases, valence bond solids, spin liquids and quantum computation that are pertinent to an understanding of the Kitaev model. The properties of the Kitaev model without disorder are discussed and it is found to be a critical spin liquid, with algebraic correlations in two spin operators sigma^{alpha}_{i}sigma^{alpha}_{j}, where i and j are either end of a link of type alpha = x, y or z on the honeycomb lattice. The Kitaev model is exactly solvable and we show that this remains so in the presence of site dilution and exchange randomness. We find that vacancies bind a flux. In the gapped phase, a vacancy forms an effective paramagnetic moment. With two or more vacancies we describe the interaction of their effective moments and show that a finite density of vacancies leads to a divergent macroscopic susceptibility at small fields. In the gapless phase the effective moment has a susceptibility that is, to leading order at small fields, chi(h)~log(1/h). Interaction between the moments from two vacancies on opposite sublattices cuts off this divergence in susceptibility at a large but finite constant. Two vacancies on the same sublattice behave quite differently and we find the combined susceptibility is parametrically larger than that of an isolated vacancy, chi(h)sim [h(log(1/h))^{3/2}]^{-1}. We also investigate the effects of slowly varying, quenched disorder in exchange coupling. We demonstrate that this does not qualitatively affect the susceptibility but show that the heat capacity C ~ T^{2/z}, where z is a measure of the disorder and increases from one with increasing disorder strength.
The Hubbard model encapsulates the physics of strongly correlated quantum systems in its most basic form. It has been studied intensively in the context of the high-temperature superconductivity. A number of novel phases were recently proposed for Hubbard-like models on the honeycomb lattice, the structure of graphene. We analyzed the Hubbard model of spin-
21 fermions on the honeycomb lattice at half-filling using large-scale quantum Monte Carlo simulations. We find that the weak coupling semimetal and the antiferromagnetic Mott insulator at strong interaction are separated by an extended gapped phase in an intermediate coupling regime. Exploring excitation gaps, various correlation functions as well as probing for flux quantization, we conclude that a quantum spin liquid, lacking any conventional order, emerges with local charge and spin correlations, best described by a resonating valence bonds state.