arxiv: v2 [cond-mat.str-el] 20 Nov 2017

Transcription

1 Unusual thermal Hall effect in a Kitaev spin liquid candidate α-rucl 3 Y. Kasahara 1, K. Sugii,, T. Ohnishi 1, M. Shimozawa, M. Yamashita, N. Kurita 3, H. Tanaka 3, J. Nasu 3, Y. Motome, T. Shibauchi 5, and Y. Matsuda 1 1 Department of Physics, Kyoto University, Kyoto 66-5, Japan Institute for Solid State Physics, University of Tokyo, Kashiwa 77-51, Japan 3 Department of Physics, Tokyo Institute of Technology, Meguro, Tokyo , Japan Department of Applied Physics, University of Tokyo, Bunkyo, Tokyo , Japan and 5 Department of Advanced Materials Science, University of Tokyo, Chiba , Japan arxiv:179.16v [cond-mat.str-el] Nov 17 The Kitaev quantum spin liquid (QSL) displays the fractionalization of quantum spins into Majorana fermions. The emergent Majorana edge current is predicted to manifest itself in the form of a finite thermal Hall effect, a feature commonly discussed in topological superconductors. Here we report on thermal Hall conductivity κ xy measurements in α-rucl 3, a candidate Kitaev magnet with the two-dimensional (D) honeycomb lattice. In a spin-liquid (Kitaev paramagnetic) state below the temperature characterized by the Kitaev interaction J K/k B K, positive κ xy develops gradually upon cooling, demonstrating the presence of highly unusual itinerant excitations. Although the zero-temperature property is masked by the magnetic ordering at = 7K, the sign, magnitude, and T-dependence of κ xy/t at intermediate temperatures follows the predicted trend of the itinerant Majorana excitations. Quantum spin liquids (QSLs) are novel states of matter lacking magnetic order while possessing some special patterns of quantum mechanical entanglement [1]. Among the long-standing experimental challenges associated with the identification of these exotic states, is the detection of fractionalized excitations, which are signatures of topological order inherent to the QSL. Recently, the Kitaev spin model of insulating magnets on a D honeycomb lattice has attracted interest, as it hosts a quantum spin liquid (QSL) where spins are fractionalized into two types of Majorana fermions; one is localized and composes localized Z fluxes with finite gap of G.6J K, and the other is itinerant (mobile) and gapless at zero field [Fig.1(a)] [, 3]. A spin-orbital-assisted Mott insulator α-rucl 3, in which local j eff = 1/ spins are aligned in twodimensional (D) honeycomb layers, has been recently considered to be a Kitaev candidate material [Fig. 1]. Although α-rucl 3 exhibits an antiferromagnetic ordering with zigzag spin structure at 7K [] due to non-kitaev interactions, such as Heisenberg exchange and off-diagonal interactions [Fig.1(c)], J K is predominant among the magnetic interactions [5 1]. The signature of the spin liquid, in which spins are not ordered but strongly entangled due to Kitaev interaction, is expected to appear in the temperature range between and J K /k B. The unusual broad magnetic excitations reported by Raman and neutron inelastic scattering measurements have been suggested to reflect the proximity to the Kitaev model [ 11]. However, the interpretation of these excitations in terms of itinerant Majorana fermions is controversial [3, 1]. Thus, the further investigation of excitations is strongly desired. Thermal transport experiments constitute sensitive probes to study low-energy itinerant quasiparticle excitations, not influenced by localized ones. In conventional insulating magnets, the longitudinal thermal conductivity is observed due to bosonic magnon excitations in the ordered state, but it is strongly suppressed in the conventional paramagnetic state. On the other hand, in the disordered QSL states realized in 1D chain, and D triangular and 3D pyrochlore systems, large thermal conductivity is observed, which has been attributed to fractionalized exotic quasiparticles, such as spinons and monopoles [13 15]. Recently, the thermal Hall effect, thermal analogue of the electronic Hall effect, has aroused great attention as a unique probe of non-trivial topological excitations [16 ]. In particular, it has been predicted that in the spin-liquid state of Kitaev magnets under an external field thermal Hall effect can be generated by Majorana fermion excitations characterized by the non-zero topological Chern number. Therefore, thermal Hall effect is a powerful probe in the search for Majorana fermions as commonly discussed in quantum Hall systems and chiral p-wave superconductors [ ]. Here we report the observation of a large κ xy with positive sign in the spin-liquid (Kitaev paramagnetic) state of α-rucl 3 at < T < J K /k B. We find that the observed thermal Hall response is essentially different from the other insulating magnets where finite κ xy has been reported, demonstrating the presence of highly unusual itinerant excitations inherent to the spin-liquid state of α-rucl 3. High-quality single crystals of α-rucl 3 were grown by a vertical Bridgman method [5]. The crystal quality has been checked by magnetic susceptibility and specific heat measurements [Fig.1(d)]. For thermal transport measurements, we carefully picked up thin crystals with plate-like shape. In-plane thermal conductivity and thermal Hall conductivity were measured simultaneously on the same crystal by the standard steady state method, using the experimental setup illustrated in the inset of

2 FIG. 1. (a) Schematic of Kitaev model on honeycomb lattice. Blue, green, and red bonds represent Ising-like bonddirectional interactions between x-, y-, and z-components of the 1/ spins at the apexes of the hexagons, respectively. The quantum spins are fractionalized into itinerant Majorana fermions (purple spheres) and Z fluxes with W p = ±1 (orange hexagons are excited fluxes for W p = 1). Crystal structure of α-rucl 3 in the ab plane. Edge-sharing RuCl 6 octahedra forms D honeycomb lattice of Ru ions [Ru 3+ (d 5 ), j eff = 1/]. (c) Schematic phase diagram for the D Kitaev model as a function of g, where g is the ratio of the Kitaev interaction J K and non-kitaev interactions, such as Heisenberg and off-diagonal terms (in the pure Kitaev model, g = ). At low g, the ground state is a QSL. α-rucl 3 corresponds to finite g with magnetic ordering. There are two characteristic temperatures, J K/k B and G/k B.6J K/k B. At T > J K/k B (regime-iii), the system is in a conventional paramagnetic state. At T < J K/k B, the spin-liquid (Kitaev paramagnetic) state appears. Above G/k B (regime-ii), the itinerant Majorana fermions are scattered by the thermally excited Z fluxes. Below G/k B (regime-i), the numbers of the excited Z fluxes are strongly reduced. (d) Temperature dependence of magnetic susceptibility χ (left axis) measured in an in-plane magnetic field of H = 1Oe and specific heat C in zero field (right axis) for α-rucl 3 single crystals. Fig.(a). Heat current q and magnetic field H were applied along the ab plane (q x) and c axis (H c z), respectively. The temperature gradients x T x and y T y were measured by carefully calibrated thermometers. For the measurements of the thermal Hall effect, we removed the longitudinal response due to misalignment of the contacts by anti-symmetrizing the measured y T as y T asym (H) = [ y T(H) y T( H)]/at each temperature. κ xx and κ xy were obtained from longitudinal thermal resistivity w xx = x T/q and thermal Hall resistivity w xy = y T asym /q as κ xx = w xx /(wxx + wyy ) and κ xy = w xy /(wxx +w yy ). To avoid background Hall signal, a LiF heat bath and non-metallic grease were (a) κ xx (W/Km) κ xx (H)/κ xx () Sample H q Ty µ H (T) heat bath 3 heater Tx thermometer Sample 1 T 1 T Sample T 1 1 T 1 T 5 1 N K 35 K 7 K K 1 K FIG.. (a) Temperature dependence of κ xx at T and 1T for two different samples. Lower inset shows expanded view of κ xx(t) at low temperatures. Solid arrows indicate the AFM transition temperature 7K. No anomaly is observed at 1 K at which the transition associated with the stacking fault occurs (dashed arrow). Upper inset illustrates a schematic of the measurement setup for the in-plane thermal conductivity and thermal Hall conductivity. Field dependence of magneto-thermal conductivity κ xx(h)/κ xx()[ κ(h) κ() κ xx() ] in the spin-liquid state. used [19]. We note that both κ xx (T) and κ xy (T) measured on two different single crystals and by using setups with different calibrated thermometers well coincide with each other. Figure (a) depicts the in-plane thermal conductivity κ xx (T) of α-rucl 3 single crystals in zero and finite magnetic fields (H) perpendicular to the D planes. The overall T-dependence is similar to the previous reports [6 ]. At, κ xx exhibits a sharp kink anomaly. Although the magnetic susceptibility shows a tiny anomaly at 1K due to the transition associated with staking faults in the crystals [Fig.1(d)], κ xx and specific heat exhibit no corresponding anomaly [Figs. 1(d) and (a)], indicating that this effect is negligible in the present analysis. These results, along with the fact that the magnitude ofκ xx in the presentcrystalsis-3times largerthan the previously reported values, indicate the high quality of the present crystals. 7 16

3 3 (a) 1 κ xy (1-3 W/Km) Sample -1 1 µ H (T) 37 K K 15 K 1 K 7 K 1 κ xy (1-3 W/Km) κ xy (1-3 W/Km) µ H (T) 1 K 7 K µ H = 1 T Sample 1 Sample 1 15 (a) 6 κ xy /T (1 - W/K m) LRO - - Sample 1 (6 T) Sample 1 (1 T) Sample (15 T) 6 π/ (κ xy D /T)/(k B /h) I II III h*/j K = T/(J K /k B ) π/ (κ xy D /T)/(k B /h) FIG. 3. (a) Field dependence of κ xy for sample. Temperature dependence of κ xy near. Inset shows the field dependence of κ xy below and above. As shown in Fig., κ xx (H) decreases with H. In the present system, the heat is carried by spin excitations and phonons; κ xx = κ sp xx+κ ph xx. As the magnetic field enhances the phonon mean free path due to the alignment ofspins, suppressingthe spin-phononscatteringrate, κ ph xx is expected to increase with H. Thus the observed reduction of κ xx (H) with H is caused by κ sp xx. However, lacking detailed information on spin-phonon scattering, which is suggested to be present in this material [6], the quantitative contribution of spin excitations is difficult to be extracted from κ xx (H). Thermal Hall effect is, therefore, a more suitable probe to unveil the nature of spin excitations. As shown in Fig.3(a), κ xy, obtained by antisymmetrizing thermal response with respect to the field direction, is clearly resolved, establishing a finite thermal Hall effect. Since the magnitude of κ xy is three orders of magnitude smaller than κ xx even at 1T, special care was taken to detect the intrinsic thermal Hall signal [19]. The observed κ xy is slightly non-linear in H for all measured temperatures. Figure 3 and its inset depict the T- and H-dependences of κ xy, respectively, at low temperatures below and above. In the spin-liquid state above, κ xy is positive and displays slight upward curvature in the field dependence, while it becomes negative and exhibits slight downward curvature or H-linear dependence upon entering the ordered state. Figure (a) displays κ xy /T in a wide temperature range. Above, κ xy /T increases steeply with T and decreases gradually after reaching a maximum at around K. No discernible Hall signal is observed within our experimental resolution in the conventional paramagnetic state above K. The changes in Hall sign and the sign of curvature in the H-dependence of κ xy when crossing around FIG.. (a) Temperature dependence of κ xy/t of α-rucl 3 in several magnetic fields in the spin-liquid and magnetically ordered states. The right scale represents (κ D xy/t)/(kb/ ), where κ D xy = κ xyd (d = 5.7Å is the interlayer distance []). Violet dashed line represents the quantized thermal Hall effect; (κ D xy/t)/(kb/ ) = π/1. Numerical results of (κ D xy/t)/(kb/ ) plotted as a function of T/(J K/k B) for the pure Kitaev model in the presence of the effective magnetic field h proportional to H 3, which is introduced by treating H as a perturbation [1]. Inset shows the same plot on a semi-logarithmic scale of T/(J K/k B). Regimes I and II represent the T-regime below and above G/k B in the spin-liquid state, respectively, and regime III represents the conventional paramagnetic state. suggest that the thermal Hall effect in the ordered and spin-liquid states is different in origin. This is supported by the following argument. The finite κ xy in the ordered state appears as a result of the magnon Hall effect in the presence of the spin chirality arising from the Dzyaloshinsky-Moriya (DM) interaction, which is allowed by symmetry for the second neighbor bonds [16, 17,, 9, 3]. The magnitude of DM interaction of α-rucl 3 is, however, reported to be D/k B.6J K /k B 5K [31], while on the other hand the prominent positive Hall signal is found up to a much higher temperature scale. We stress that in usual Heisenberg systems thermal Hall effect is absent in spin-liquid states in the absence of DM interaction []; the finite Hall effect is predicted in exotic spin-liquid state with emergent spinon Fermi surface, but it is unlikely that such a situation is realized in the present system [1, 16, 19]. Moreover, the phonon thermal Hall conductivity is two orders of magnitude smaller than the observed κ xy /T in the spin-liquid state and shows essentially different T-dependence[3, 33]. Finite thermal Hall signal in the spin systems has been reported in Lu V O 7 [17], Cu(1,3-bdc) [1], Tb Ti O 7 [3], and volborthite Cu 3 V O 7 (OH) H O [19]. However, we point out that

4 the origin of κ xy in α-rucl 3 is essentially different from these systems. In fact, κ xy is observed only in the ferromagnetically ordered state in Lu V O 7, and a peak-like featureofκ xy (T)/T andchangeinhallsignappearbelow ferromagnetic ordering temperature in Cu(1,3-bdc). In Tb Ti O 7, finite κ xy appears in the paramagnetic state far above the characteristic temperature scale of the exchange interaction. In volborthite, although similar peak structureofκ xy (T)/T isobservedinthespin-liquid state, the Hall sign is negative and the amplitude of κ xy /T is one order of magnitude smaller than in α-rucl 3. These results lead us to conclude that the observed thermal Hall effect above reflects unusual quasiparticle excitations inherent to the spin-liquid state of α-rucl 3. In the Kitaev magnets itinerant Majorana fermions give rise to the finite thermal Hall effect. The sign of κ xy is determined by the topological Chern number and is predicted to be positive[], which is consistent with the present results. Moreover, the non-linear H-dependence with upward curvature of κ xy (H) is predicted [1], which is also consistent with the observedκ xy (H) above. In the zero-temperature limit, the T-linear coefficient of the thermal Hall conductivity per D layer κ D xy/t is quantized as κ D xy/t = (π/1)(kb / ) [, 1]. The right axis of Fig.(a) plots (κ D xy /T)/(k B / ), where κd xy = κ xyd (d is the interlayer distance). The result is highlighted by the fact that (κ D xy/t)/(kb / ) at K reaches about a half of π/1 at T. Considering the occurrence of the magnetic ordering and the presence of non-kitaev terms, this observation is rather remarkable and is quite telling of the nature of the spin-liquid phase of this material. Figure displays the numerical results of κ xy (T)/T for the D pure Kitaev model calculated by the quantum Monte Carlo method [1]. There are three distinct temperature regimes. In the high-temperature conventional paramagnetic state (regime-iii), κ xy /T vanishes. In the spin-liquid state below J K /k B (regime-ii), κ xy (T)/T is generated by itinerant Majorana fermions, which are scattered by the thermally excited Z fluxes. A broad peak of κ xy /T at T.1J K /k B appears as a result of the interaction between the two excitations; itinerant Majorana fermions and Z fluxes [Fig.1(a)]. At low temperatures at T < G /k B (regime-i), the rapid reduction of the number of the excited Z fluxes leads to a steep increase of κ xy /T, approaching the quantized value at T. The numerical calculations do not take into account the non-kitaev terms, and thus the occurrence of magnetic ordering. Apart from the low-temperature range where the effects of magnetic order take place, the experimental curve of κ xy (T)/T follows the trend of the calculated results. Antiferromagnetic fluctuations above would suppress the Kitaev excitations, which may be partly responsible for the reduction of κ xy (T)/T near. However, whether such an effect is significant up to much higher temperature than is not clear. Although it should be scrutinised by other probes, the fact that the thermal conductivity, magnetic susceptibility and specific heat show no critical features associated with the transition except for very near, suggests that the experimentally observedpeak ofκ xy /T at much higher temperature K (.5J K /k B ) likely has the same origin as the peak feature found in the calculation. Even if the antiferromagnetic fluctuations originating from the ordered state below influence the thermal Hall signal, we can safely conclude that κ xy /T in the high temperature regime above K, where κ xy /T increases with decreasing T, is most likely to be attributed to unusual spin excitations inherent to the Kitaev spin liquid. Although the zero-temperature property is masked by the magnetic ordering, the magnitude of thermal Hall signal as well as its sign and its T-dependence above, including a peak structure are well reproduced by the Monte Carlo calculation. Therefore, the observed thermal Hall signal appear to provide a signature of the Majorana fermion excitations in the spin-liquid (Kitaev paramagnetic) state, where itinerant topological Majorana fermions propagate in the background of thermallyexcited Z fluxes. In summary, we have investigated the low-energy itinerant quasiparticle excitations in the Kitaev spin liquid candidate α-rucl 3 by thermal Hall effect. Our measurements reveal a large thermal Hall signal as well as the anomalous T-dependence in the spin-liquid state of α- RuCl 3, which is markedly different from the other insulating magnets. Although the anomalous behavior is masked by the magnetic ordering below, the sign, magnitude, and T-dependence of the observed κ xy /T above are well reproduced by the calculated κ xy /T for the D pure Kitaev model, suggesting that the observed thermal Hall response arises from quasiparticle excitations inherent to the Kitaev spin liquid. We thank S. Fujimoto, H. Katsura, E.-G. Moon, H. Lee, K. Nomura, Y. Tokiwa, K. Totsuka, and M. Udagawa, for useful discussions. This work was supported by Grants-in-Aid for Scientific Research (KAKENHI) (No. 571, 15H1, 15H16, and 15H557) and Grants-in-Aid for Scientific Research on innovative areas Topological Materials Science (No. JP15H55) from Japan Society for the Promotion of Science (JSPS), and by Yamada Science Foundation and Toray Science Foundation. These authors contributed equally to this work. [1] L. Balents, Nature 6, 199 (1). [] A. Kitaev, Ann. Phys. 31, (6). [3] S. Trebst, arxiv: [] R. D. Johnson, S. C. Williams, A. A. Haghighirad, J. Singleton, V. Zapf, P. Manuel, I. I. Mazin, Y. Li, H. O. Jeschke, R. Valentí, R. Coldea, Phys. Rev. B 9, 3119

5 5 (15). [5] K. W. Plumb, J. P. Clancy, L. J. Sandilands, V. V. Shankar, Y. F. Hu, K. S. Burch, H.-Y. Kee, Y.-J. Kim, Phys. Rev. B 9, 111(R) (1). [6] H.-S. Kim, V. V. Shankar, A. Catuneanu, H.-Y. Kee, Phys. Rev. B 91, 111(R) (15). [7] G. Jackeli, G. Khaliullin, Phys. Rev. Lett. 1, 175 (9). [] A. Banerjee, C. A. Bridges, J.-Q. Yan, A. A. Aczel, L. Li, M. B.Stone, G. E. Granoroth, M. D.Lumsden, Y.Yiu, J. Knolle, S. Bhattacharjee, D. L. Kovrizhin, R. Moessner, D. A. Tennant, D. G. Mandrus, S. E. Nagler, Nat. Mat. 15, 733 (16). [9] Phys. Rev. Lett. 11, 171 (15). [1] J. Nasu, J. Knolle, D. L. Kovrizhin, Y. Motome, R. Moessner, Nat. Phys. 1, 91 (16). [11] J. Yoshitake, J. Nasu, Y. Motome, Phys. Rev. Lett. 117, 1573 (16). [1] S. Winter, K. Riedl, A. Honecker, R. Valentí, Nat. Commun., 115 (17). [13] A. V. Sologubenko, K. Giannò, H. R. Ott, A. Vietkine, A. Revcolevschi, Phys. Rev. B 6, 51 (1). [1] M. Yamashita, N. Nakata, Y. Senshu, M. Nagata, H. M. Yamamoto, R. Kato, T. Shibauchi, Y. Matsuda, Science 3, 16 (1). [15] Y. Tokiwa, T. Yamashita, M. Udagawa, S. Kittaka, T. Sakakibara, D. Terazawa, Y. Shimoyama, T. Terashima, Y. Yasui, T. Shibauchi, Y. Matsuda, Nat. Commun. 7, 17 (16). [16] H. Katsura, N. Nagaosa, P. A. Lee, Phys. Rev. Lett. 1, 663 (1). [17] Y. Onose, T. Ideue, H. Katsura, Y. Shiomi, N. Nagaosa, Y. Tokura, Science 39, 97 (1). [1] M. Hirschberger, R. Chisnell, Y. S. Lee, N. P. Ong, Phys. Rev. Lett. 115, 1663 (15). [19] D. Watanabe, K. Sugii, M. Shimozawa, Y. Suzuki, T. Yajima, H. Ishikawa, Z. Hiroi, T. Shibauchi, Y. Matsuda, M. Yamashita, Proc. Natl. Acad. Sci. USA 113, 653 (16). [] J. H. Han, H. Lee, J. Phys. Soc. Jpn. 6, 117 (16). [1] J. Nasu, J. Yoshitake, Y. Motome, Phys. Rev. Lett. 119, 17 (17). [] N. Read, D. Green, Phys. Rev. B 61, 167 (). [3] K. Nomura, S. Ryu, A. Furusaki, N. Nagaosa, Phys. Rev. Lett. 1, 6 (1). [] H. Sumiyoshi, S. Fujimoto, J. Phys. Soc. Jpn., 36 (13). [5] Y. Kubota, H. Tanaka, T. Ono, Y. Narumi, K. Kindo, Phys. Rev. B 91, 9 (15). [6] R. Hentrich, A. U. B. Wolter, X. Zotos, W. Brenig, D. Nowak, A. Isaeva, T. Doert, A. Banerjee, P. Lampen- Kelly, D. G. Mandrus, S. E. Nagler, J. Sears, Y.-J. Kim, B. Büchner, C. Hess, arxiv: [7] D. Hirobe, M. Sato, Y. Shiomi, H. Tanaka, E. Saitoh, Phys. Rev. B 95, 111(R) (17). [] I. A. Leahy, C. A. Pocs, P. E. Siegfried, D. Graf, S.-H. Do, K.-Y. Choi, B. Normand, M. Lee, Phys. Rev. Lett. 11, 173 (17). [9] S. A. Owerre, J. Appl. Phys. 1, 393 (16). [3] V. A. Zyuzin and A. A. Kovalev, Phys. Rev. Lett. 117, 173 (16). [31] S. M. Winter, Y. Li, H. O. Jeschke, R. Valentí, Phys. Rev. B 93, 131 (16). [3] C. Strohm, G. L. J. A. Rikken, P. Wyder, Phys. Rev. Lett. 95, (5). [33] K. Sugii, M. Shimozawa, D. Watanabe, Y. Suzuki, M. Halim, M. Kimata, Y. Matsumoto, S. Nakatsuji, M. Yamashita, Phys. Rev. Lett. 11, 159 (17). [3] M. Hirschberger, J. W. Krizan, R. J. Cava, N. P. Ong, Science 3, 16 (15).