Finite-Temperature Phase Transition to a Quantum Spin Liquid in a Three-Dimensional Kitaev Model on a Hyperhoneycomb Lattice

We present numerical evidence for the presence of a finite-temperature ($T$) phase transition separating paramagnet and quantum spin liquid in a three-dimensional variant of the Kitaev model defined on a hyperhoneycomb lattice in the limit of strong anisotropy; the model is mapped onto an effective Ising-type model, where elementary excitations consist of closed loops of flipped Ising-type variables. Analyzing this effective model by Monte Carlo simulation, we find a phase transition from quantum spin liquid to paramagnet at a finite critical temperature $T_c$. We also compute the magnetic properties in terms of the original quantum spins. We find that the magnetic susceptibility exhibits a broad hump above $T_c$, while it obeys the Curie law at high $T$ and approaches a nonzero Van Vleck-type constant at low $T$. Although the susceptibility changes continuously at $T_c$, its $T$ derivative shows critical divergence at $T_c$. We also clarify that the dynamical spin correlation function is momentum independent but shows quantized peaks corresponding to the discretized excitations. Although the phase transition accompanies no apparent symmetry breaking in terms of the Ising-type variables as well as the original quantum spins, we characterize it from a topological viewpoint. We find that, by defining the flux density for loops of the Ising-type variables, the transition is interpreted as the one occurring from the zero-flux quantum spin liquid to the nonzero-flux paramagnet; the latter has a Coulombic nature due to the local constraints. The role of global constraints on the Ising-type variables is examined in comparison with the results in the two-dimensional loop model. A correspondence of our model to the Ising model on a diamond lattice is also discussed. A possible relevance of our results to the recently-discovered hyperhoneycomb compound, $beta$-Li$_2$IrO$_3$, is mentioned.