We calculate dynamical spin structure factors (DSFs) of three ab-initio models and one ab-initio-guided model for a honeycomb-lattice magnet
α-RuCl
3 with an exact numerical diagonalization method. We also calculate temperature dependences of the heat capacity, the nearest-neighbor (NN) spin-spin correlation function, and the static spin structure factor by employing thermal pure quantum states. From the obtained results for the four effective models, we find that, even when the magnetic order appears at a low temperature, the intensity at the
Γ point in the DSFs grows as the NN spin correlation grows. In addition, we find that the four models fail in explaining the features of both the inelastic-neutron-scattering experiments and the heat capacity measurements. In the four models, when temperature decreases, the heat capacity shows a prominent peak at a high temperature where the nearest-neighbor spin-spin correlation function grows. However, the peak temperatures in the heat capacity are too low in comparison with the observed values in the experiments. Thus, we propose an effective model with the strong ferromagnetic Kitaev coupling. We show that the proposed model quantitatively reproduces the above two experimental features. To discuss the further adequacy of the proposed model, we calculate the field dependence of the polarized terahertz spectra. We find that the proposed model successfully reproduces the experimental results: the spin-gapped excitation surviving up to an onset field where the magnetic order disappears and the almost linear response in the high field region. From the obtained numerical results, we argue that the feature of the low-energy magnetic excitation in
α-RuCl
3 is mainly characterized by the other interactions rather than the most dominant Kitaev interactions such as the off-diagonal interactions on the nearest-neighbor pairs.
Contrary to the original expectation,
Na2IrO3 is not a Kitaev's quantum spin liquid (QSL) but shows a zigzag-type antiferromagnetic order in experiments. Here, we propose experimental clues and criteria to measure how a material in hand is close to the Kitaev's QSL state. For this purpose, we systematically study thermal and spin excitations of a generalized Kitaev-Heisenberg model studied by Chaloupka et al., Phys. Rev. Lett. 110, 097204 (2013) and an effective ab initio Hamiltonian for
Na2IrO3 proposed by Yamaji et al., Phys. Rev. Lett. 113, 107201 (2014), by employing a numerical diagonalization method. We reveal that closeness to the Kitaev's QSL is characterized by the following properties, besides trivial criteria such as reduction of magnetic ordered moments and N\'eel temperatures. (1) Two peaks in the temperature dependence of specific heat at ${T}_{\ensuremath{\ell}}$ and
Th caused by the fractionalization of spin to two types of Majorana fermions. (2) In between the double peak, a prominent plateau or shoulder pinned at
2Rln2 in the temperature dependence of entropy, where
R is the gas constant. (3) Failure of the linear spin wave approximation at the low-lying excitations of dynamical structure factors. (4) Small ratio ${T}_{\ensuremath{\ell}}/{T}_{h}$ close to or less than 0.03. According to the proposed criteria,
Na2IrO3 is categorized to a compound close to the Kitaev's QSL, and is proven to be a promising candidate for the realization of the QSL if the relevant material parameters can further be tuned by making thin film of
Na2IrO3 on various substrates or applying axial pressure perpendicular to the honeycomb networks of iridium ions. Applications of these characterization to ${({\mathrm{Na}}_{1\ensuremath{-}x}{\mathrm{Li}}_{x})}_{2}{\mathrm{IrO}}_{3}$ and other related materials are also discussed.
Honeycomb-lattice magnet α-RuCl 3 is considered to be a potential candidate of realizing Kitaev spin liquid, although this material undergoes a phase transition to the zigzag magnetically ordered state at T N ∼ 7 K. Quite recently, inelastic neutron-scattering experiments using single crystal α-RuCl 3 have unveiled characteristic dynamical properties. We calculate dynamical spin structure factors of three ab-initio models for α-RuCl 3 with an exact numerical diagonalization method. We also calculate temperature dependences of the specific heat by employing thermal pure quantum states. We compare our numerical results with the experiments and discuss characteristics obtained by using three ab-initio models.
An effective low-energy Hamiltonian of itinerant electrons for iridium oxide Na2IrO3 is derived by an ab initio downfolding scheme. The model is then reduced to an effective spin model on a honeycomb lattice by the strong coupling expansion. Here we show that the ab initio model contains spin-spin anisotropic exchange terms in addition to the extensively studied Kitaev and Heisenberg exchange interactions, and allows us to describe the experimentally observed zigzag magnetic order, interpreted as the state stabilized by the antiferromagnetic coupling of the ferromagnetic chains. We clarify possible routes to realize quantum spin liquids from existing Na2IrO3.
We use an exact numerical diagonalization method to calculate the dynamical spin structure factors of three ab initio models and one ab initio guided model for a honeycomb-lattice magnet $\ensuremath{\alpha}\text{\ensuremath{-}}{\mathrm{RuCl}}_{3}$. We also use thermal pure quantum states to calculate the temperature dependence of the heat capacity, the nearest-neighbor spin-spin correlation function, and the static spin structure factor. From the results obtained from these four effective models, we find that, even when the magnetic order is stabilized at low temperature, the intensity at the $\mathrm{\ensuremath{\Gamma}}$ point in the dynamical spin structure factors increases with increasing nearest-neighbor spin correlation. In addition, we find that the four models fail to explain heat-capacity measurements whereas two of the four models succeed in explaining inelastic-neutron-scattering experiments. In the four models, when temperature decreases, the heat capacity shows a prominent peak at a high temperature where the nearest-neighbor spin-spin correlation function increases. However, the peak temperature in heat capacity is too low in comparison with that observed experimentally. To address these discrepancies, we propose an effective model that includes strong ferromagnetic Kitaev coupling, and we show that this model quantitatively reproduces both inelastic-neutron-scattering experiments and heat-capacity measurements. To further examine the adequacy of the proposed model, we calculate the field dependence of the polarized terahertz spectra, which reproduces the experimental results: the spin-gapped excitation survives up to an onset field where the magnetic order disappears and the response in the high-field region is almost linear. Based on these numerical results, we argue that the low-energy magnetic excitation in $\ensuremath{\alpha}\text{\ensuremath{-}}{\mathrm{RuCl}}_{3}$ is mainly characterized by interactions such as off-diagonal interactions and weak Heisenberg interactions between nearest-neighbor pairs, rather than by the strong Kitaev interactions.
The recently discovered three-dimensional hyperhoneycomb iridate, β -Li 2 IrO 3 , has raised hopes for the realization of the dominant Kitaev interaction between spin-orbit entangled local moments due to its near-ideal lattice structure. If true, this material may lie close to the sought-after quantum spin-liquid phase in three dimensions. Utilizing ab initio electronic structure calculations, we first show that the spin-orbit entangled basis, , correctly captures the low-energy electronic structure. The effective spin model derived in the strong-coupling limit supplemented by the ab initio results is shown to be dominated by the Kitaev interaction. We demonstrated that the possible range of parameters is consistent with a non-coplanar spiral magnetic order found in a recent experiment. All of these analyses suggest that β -Li 2 IrO 3 may be the closest among known materials to the Kitaev spin-liquid regime.