We investigate the fate of the Kitaev spin liquid (KSL) under the influence of an external magnetic field h in the [001] direction and upon tuning bond anisotropy of the Kitaev coupling Kz keeping Kx=Ky=K . Guided by density matrix renormalization group, exact diagonalization, and with insights from parton mean-field theory, we uncover a field-induced gapless-to-gapless Lifshitz transition from the nodal KSL to an intermediate gapless phase. The intermediate phase sandwiched between hc1 and hc2 , which persists for a wide range of anisotropy Kz/K>0 , is composed of weakly coupled one-dimensional quantum critical chains. This intermediate phase is a dimensional crossover, which asymptotically leads to the one-dimensional quantum Ising criticality characterized by the (1+1)D conformal field theory as the field reaches the phase transition at hc2 . Beyond hc2 the system enters a partially polarized phase describable as effectively decoupled bosonic chains in which spin waves propagate along the one-dimensional zigzag direction. Our findings provide a comprehensive phase diagram and offer insights into the unusual physics of dimensional reduction generated by a uniform magnetic field in an otherwise two-dimensional quantum spin liquid. Published by the American Physical Society 2024
We investigate the fate of the Kitaev spin liquid (KSL) under the influence of an external magnetic field h in the [001] direction and upon tuning bond anisotropy of the Kitaev coupling Kz keeping Kx=Ky=K . Guided by density matrix renormalization group, exact diagonalization, and with insights from parton mean-field theory, we uncover a field-induced gapless-to-gapless Lifshitz transition from the nodal KSL to an intermediate gapless phase. The intermediate phase sandwiched between hc1 and hc2 , which persists for a wide range of anisotropy Kz/K>0 , is composed of weakly coupled one-dimensional quantum critical chains. This intermediate phase is a dimensional crossover, which asymptotically leads to the one-dimensional quantum Ising criticality characterized by the (1+1)D conformal field theory as the field reaches the phase transition at hc2 . Beyond hc2 the system enters a partially polarized phase describable as effectively decoupled bosonic chains in which spin waves propagate along the one-dimensional zigzag direction. Our findings provide a comprehensive phase diagram and offer insights into the unusual physics of dimensional reduction generated by a uniform magnetic field in an otherwise two-dimensional quantum spin liquid. Published by the American Physical Society 2024
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
Kitaev magnets have emerged as pivotal systems for investigating frustrated magnetism, providing a unique platform to explore quantum phases governed by the interplay between bond-dependent anisotropy and external magnetic fields. However, the quantum phase diagrams, particularly near the dominant antiferromagnetic Kitaev regime, remain puzzling despite extensive studies. In this work, we perform unbiased exact diagonalization calculations of the Kitaev-Γ model in a [111] magnetic field on a C6-symmetric 24-site cluster. By calculating the Z2 flux density and the topological entanglement entropy, we reveal multiple phase transitions and identify signatures of both scalar and vector chiral orders in the intermediate-field regime between the Kitaev spin liquid and the polarized phase. As the negative Γ interaction increases, we discover a proximate quantum spin liquid featured by a three-peak specific heat and a spin-flop phase at a moderate magnetic field. Our findings provide insight into the field-induced intermediate phases in the antiferromagnetic Kitaev model and pave the way for the hunt for emergent phases in real materials.
The Kitaev model with an applied magnetic field in the [Formula: see text] direction shows two transitions: from a nonabelian gapped quantum spin liquid (QSL) to a gapless QSL at [Formula: see text] and a second transition at a higher field [Formula: see text] to a gapped partially polarized phase, where K is the strength of the Kitaev exchange interaction. We identify the intermediate phase to be a gapless U(1) QSL and determine the spin structure function [Formula: see text] and the Fermi surface [Formula: see text] of the gapless spinons using the density matrix renormalization group (DMRG) method for large honeycomb clusters. Further calculations of static spin-spin correlations, magnetization, spin susceptibility, and finite temperature-specific heat and entropy corroborate the gapped and gapless nature of the different field-dependent phases. In the intermediate phase, the spin-spin correlations decay as a power law with distance, indicative of a gapless phase.
Topological spin liquids in two spatial dimensions are stable phases in the presence of a small magnetic field, but may give way to field-induced phenomena at intermediate field strengths. Sandwiched between the low-field spin liquid physics and the high-field spin-polarized phase, the exploration of magnetic phenomena in this intermediate regime, however, often remains elusive to controlled analytical approaches. Here we numerically study such intermediate-field magnetic phenomena for two representative Kitaev models (on the square-octagon and decorated honeycomb lattice) that exhibit either Abelian or non-Abelian topological order in the low-field limit. Using a combination of exact diagonalization and density matrix renormalization group techniques, as well as linear spin-wave theory, we establish the generic features of Kitaev spin liquids in an external magnetic field. While ferromagnetic models typically exhibit a direct transition to the polarized state at a relatively low field strength, antiferromagnetic couplings not only substantially stabilizes the topological spin liquid phase, but generically lead to the emergence of a distinct field-induced intermediate regime, separated by a crossover from the high-field polarized regime. Our results suggest that, for most lattice geometries, this regime generically exhibits significant spin canting, antiferromagnetic spin-spin correlations, and an extended proximate spin liquid regime at finite temperatures. Notably, we identify a symmetry obstruction in the original honeycomb Kitaev model that prevents, at least for certain field directions, the formation of such canted magnetism without breaking symmetries---consistent with the recent numerical observation of an extended gapless spin liquid in this case.
Motivated by intense research on two-dimensional spin-21 Kitaev materials, Kitaev spin chains and ladders, though geometrically limited, have been studied for their numerical simplicity and insights into extended Kitaev models. The phase diagrams under the magnetic field were also explored for these quasi-one-dimensional models. For an isotropic Kitaev chain, it was found that a magnetic field polarizes the ground state except along the symmetric field angle, where the chain is found to remain gapless up to a critical field strength where it enters an intriguing soliton phase before reaching the polarized state at higher field strengths. Here we study an anisotropic Kitaev chain under a magnetic field using the density matrix renormalization group technique, where the ground state has a macroscopic degeneracy with a finite gap in the absence of the magnetic field. When the field is mainly aligned parallel to the strong bond, four-site and large unit-cell ordered phases arise. In a certain angle of the field, another ordered phase characterized by a uniform chirality with six-site periodicity emerges. We employ a perturbation theory to understand such field-induced ordered phases. The effective model uncovers the presence of transverse Ising and Dzyaloshinskii-Moriya interactions between unit cells, as well as further-neighbor Ising interaction induced by the magnetic field, which collectively explain the mechanisms behind these ordered states. Open questions and challenges are also discussed.
Recent numerical studies indicate that the antiferromagnetic Kitaev honeycomb lattice model undergoes a magnetic-field-induced quantum phase transition into a new spin-liquid phase. This intermediate-field phase has been previously characterized as a gapless spin liquid. By implementing a recently developed variational approach based on the exact fractionalized excitations of the zero-field model, we demonstrate that the field-induced spin liquid is gapped and belongs to Kitaev's 16-fold way. Specifically, the low-field non-Abelian liquid with Chern number C = ±1 transitions into an Abelian liquid with C = ±4. The critical field and the field-dependent behaviors of key physical quantities are in good quantitative agreement with published numerical results. Furthermore, we derive an effective field theory for the field-induced critical point which readily explains the ostensibly gapless nature of the intermediate-field spin liquid.