No Arabic abstract
We study the effects of doping the Kitaev model on the honeycomb lattice where the spins interact via the bond-directional interaction $J_K$, which is known to have a quantum spin liquid as its exact ground state. The effect of hole doping is studied within the $t$-$J_K$ model on a three-leg cylinder using density-matrix renormalization group. Upon light doping, we find that the ground state of the system has quasi-long-range charge-density-wave correlations but short-range single-particle correlations. The dominant pairing channel is the even-parity superconducting pair-pair correlations with $d$-wave-like symmetry, which oscillate in sign as a function of separation with a period equal to that of the spin-density wave and two times the charge-density wave. Although these correlations fall rapidly (possibly exponentially) at long distances, this is never-the-less the first example where a pair-density wave is the strongest SC order on a bipartite lattice. Our results may be relevant to ${rm Na_2IrO_3}$ and $alpha$-${rm RuCl_3}$ upon doping.
In addition to low-energy spin fluctuations, which distinguish them from band insulators, Mott insulators often possess orbital degrees of freedom when crystal-field levels are partially filled. While in most situations spins and orbitals develop long-range order, the possibility for the ground state to be a quantum liquid opens new perspectives. In this paper, we provide clear evidence that the SU(4) symmetric Kugel-Khomskii model on the honeycomb lattice is a quantum spin-orbital liquid. The absence of any form of symmetry breaking - lattice or SU(N) - is supported by a combination of semiclassical and numerical approaches: flavor-wave theory, tensor network algorithm, and exact diagonalizations. In addition, all properties revealed by these methods are very accurately accounted for by a projected variational wave-function based on the pi-flux state of fermions on the honeycomb lattice at 1/4-filling. In that state, correlations are algebraic because of the presence of a Dirac point at the Fermi level, suggesting that the symmetric Kugel-Khomskii model on the honeycomb lattice is an algebraic quantum spin-orbital liquid. This model provides a good starting point to understand the recently discovered spin-orbital liquid behavior of Ba_3CuSb_2O_9. The present results also suggest to choose optical lattices with honeycomb geometry in the search for quantum liquids in ultra-cold four-color fermionic atoms.
Recent proposals for spin-1 Kitaev materials, such as honeycomb Ni oxides with heavy elements of Bi and Sb, have shown that these compounds naturally give rise to antiferromagnetic (AFM) Kitaev couplings. Conceptual interest in such AFM Kitaev systems has been sparked by the observation of a transition to a gapless $U(1)$ spin liquid at intermediate field strengths in the AFM spin-1/2 Kitaev model. However, all hitherto known spin-1/2 Kitaev materials exhibit ferromagnetic bond-directional exchanges. Here we discuss the physics of the spin-1 Kitaev model in a magnetic field and show, by extensive numerical analysis, that for AFM couplings it exhibits an extended gapless quantum spin liquid at intermediate field strengths. The close analogy to its spin-1/2 counterpart suggests that this gapless spin liquid is a $U(1)$ spin liquid with a neutral Fermi surface, that gives rise to enhanced thermal transport signatures.
We study the thermodynamic properties of modified spin-$S$ Kitaev models introduced by Baskaran, Sen and Shankar (Phys. Rev. B 78, 115116 (2008)). These models have the property that for half-odd-integer spins their eigenstates map on to those of spin-1/2 Kitaev models, with well-known highly entangled quantum spin-liquid states and Majorana fermions. For integer spins, the Hamiltonian is made out of commuting local operators. Thus, the eigenstates can be chosen to be completely unentangled between different sites, though with a significant degeneracy for each eigenstate. For half-odd-integer spins, the thermodynamic properties can be related to the spin-1/2 Kitaev models apart from an additional degeneracy. Hence we focus here on the case of integer spins. We use transfer matrix methods, high temperature expansions and Monte Carlo simulations to study the thermodynamic properties of ferromagnetic and antiferromagnetic models with spin $S=1$ and $S=2$. Apart from large residual entropies, which all the models have, we find that they can have a variety of different behaviors. Transfer matrix calculations show that for the different models, the correlation lengths can be finite as $Tto 0$, become critical as $Tto 0$ or diverge exponentially as $Tto 0$. There is a conserved $Z_2$ flux variable associated with each hexagonal plaquette which saturates at the value $+1$ as $Trightarrow0$ in all models except the $S=1$ antiferromagnet where the mean flux remains zero as $Tto 0$. We provide qualitative explanations for these results.
We present a quantu spin liquid state in a spin-1/2 honeycomb lattice with randomness in the exchange interaction. That is, we successfully introduce randomness into the organic radial-based complex and realize a random-singlet (RS) state. All magnetic and thermodynamic experimental results indicate the liquid-like behaviors, which are consistent with those expected in the RS state. These results demonstrate that the randomness or inhomogeneity in the actual systems stabilize the RS state and yield liquid-like behavior.
We investigate the generic features of the low energy dynamical spin structure factor of the Kitaev honeycomb quantum spin liquid perturbed away from its exact soluble limit by generic symmetry-allowed exchange couplings. We find that the spin gap persists in the Kitaev-Heisenberg model, but generally vanishes provided more generic symmetry-allowed interactions exist. We formulate the generic expansion of the spin operator in terms of fractionalized Majorana fermion operators according to the symmetry enriched topological order of the Kitaev spin liquid, described by its projective symmetry group. The dynamical spin structure factor displays power-law scaling bounded by Dirac cones in the vicinity of the $Gamma$, $K$ and $K$ points of the Brillouin zone, rather than the spin gap found for the exactly soluble point.