Abstract
The report of a half-quantized thermal Hall effect and oscillatory structures in the magnetothermal conductivity in the Kitaev material α-RuCl3 have sparked a strong debate on whether it is generated by Majorana fermion edge currents, spinon Fermi surface, or whether other more conventional mechanisms are at its origin. Here, we report low temperature thermal conductivity (κ) of another candidate Kitaev material, Na2Co2TeO6. The application of a magnetic field (B) along different principal axes of the crystal reveals a strong directional-dependent B impact on κ, while no evidence for mobile quasiparticles except phonons can be concluded at any field. Instead, severely scattered phonon transport prevails across the BâT phase diagram, revealing cascades of phase transitions for all B directions. Our results thus cast doubt on recent proposals for significant itinerant magnetic excitations in Na2Co2TeO6, and emphasize the importance of discriminating true spin liquid transport properties from scattered phonons in candidate materials.
Similar content being viewed by others
Introduction
Frustrated magnets refer to a set of materials that do not order down to temperatures well below the corresponding magnetic interaction strength. In some cases, the system is left with a ground state called quantum spin liquid (QSL), which is characterized by entangled spins, topological orders, and fractionalized magnetic excitations1,2,3. Kitaev QSL, proposed for bond-dependent nearest-neighbor interacting spins on the two-dimensional honeycomb lattice4, is among the rare examples of exactly solvable QSL models5,6,7,8. Recent research interest is further fueled by possible achieving Kitaev QSL in realistic materials. In particular, fingerprints of QSL-related excitations were revealed in α-RuCl39,10, both in its zero-field phase above the magnetic order11, and more intriguingly, in its field-induced quantum disordered state12.
However, those sought-after signatures for QSLs observed in α-RuCl3 always come along with alternative explanations. To be more specific on transport experiments, the report of half-integer quantization of thermal Hall conductance was once regarded as a definitive proof for the existence of Majorana fermions in the field-induced quantum disordered state in α-RuCl312. However, following experiments noticed such quantization is less reproducible, and attributed the non-quantized temperature-dependent thermal Hall effect to phonons13, or topological magnons14. On the other hand, quantum oscillations of longitudinal thermal conductivity in the quantum disordered state of α-RuCl3 suggest some fractionalised magnetic excitations forming a Fermi surface15. But such spectacle is also argued to be a result simply due to scattered phonons13,16, in particular taking into account that no itinerant magnetic excitation has been found in a previous low temperature work17. On top of all these inconsistencies are the facts that α-RuCl3 suffers from the notorious stacking fault problem18, from the dechlorination and oxidation after heat treatment19, and more generally from sample dependence20. Thus, it is indispensable to seek for another example that can contribute to clarify to what extent these results are intrinsic to a Kitaev material.
In this paper, we studied the ground state properties of a quantum magnet Na2Co2TeO6 (NCTO), which was recently investigated as a candidate Kitaev QSL21,22,23, through the prism of low temperature thermal conductivity. Similar to its more famous cousin α-RuCl3, NCTO also features an intermediate quantum disordered phase sandwiched between the low-field magnetically ordered state and the high-field polarized state22,23, where a QSL state is prospected. In this work, we report an anisotropic and astonishing large B impact on κ at low temperatures (<1 K). The analysis of our data demonstrates that there is no evidence for mobile gapless QSL excitations (κQSL) contributing to κ. Instead, strongly scattered phonon transport (κph), a quite traditional mechanism, is responsible for all the observed phenomena. Utilizing the high resolution κ(B) data, we reveal a large density of magnetic excitations in the zero-field ground state of NCTO, further promoted excitations by magnetic field perpendicular to the honeycomb plane, and possible quantum criticality driven by in-plane magnetic field.
Results and discussion
κ/T(T) in fixed magnet fields
The residual linear term of thermal conductivity extrapolated to the low temperature limit (κ0/Tââ¡âκ/Tâ£Tâ0) is persistently scrutinized for QSL candidates3,24, since its existence would suggest realization of some sought-after QSL states, including spinons forming a Fermi surface25,26,27,28,29. Figure 1a shows κ/T of five NCTO crystals in zero magnetic field as a function of T. Albeit certain sample dependence, the κ/T(T) profiles of all investigated crystals do not follow the simple power-law behavior expected for common insulators. Instead, regardless of the direction along which the thermal gradient âT is established (the two in-plane principal axes for the Co-honeycomb lattice, denoted as zigzag and armchair directions, indicated by the schematic insets), all κ/T(T) profiles trend to end with a finite κ0/T.
Among them, Sample#1 is especially notable for displaying a pronounced hump which emerges below about 250 mK. κ0/T of Sample#1, extracted from a linear fit to its data above 250 mK, as highlighted by the dashed line in Fig. 1a, is 0.048 mW/Kâ2âcmâ1. At first glance, this finite κ0/T is consistent with the notion that the ground state of NCTO is a QSL, since κQSL contributed by a putative spinon Fermi surface is expected to scale linearly with T at the lowest temperatures24. We are aware of a recent work reporting similar data of NCTO, where these results have been interpreted as evidences for κQSL30. However, as clearly shown in Fig. 1a, the upturn of κ/T(T) for Sample#1 (same for Sample#4 and Sample#5) terminates abruptly below about 100 mK, followed by a sharp drop, which results in a vanishing residual linear term that is incompatible with a gapless QSL ground state. Furthermore, we argue such observation is at odds with a hypothesised gapped QSL state with a gap size of the orderb~100 mK. When a small in-plane magnetic field is applied, κQSL/T(T) at temperatures below the initial gap energy scale is expected to increase due to a proliferation of quasiparticle excitations, while κ/T(T) at higher temperatures should be less sensitive to the field since κph is not directly affected. This is just opposite to our observations. As depicted in Fig. 1b, the field changes κ/T(T) fundamentally above 100 mK, while it has basically no impact on the data at the lowest temperature.
Higher magnetic fields have a drastic impact on the κ(T)/T curves until they finally saturate into a power-law behavior over the whole temperature range at the highest field. A representative case of Sample#1 is shown in Fig. 1c, along with the others displayed in Supplementary Note 1. As can be seen in Fig. 1d, regardless of the (in-plane) field directions, κ(T)/T curves for a given sample saturate at the same value.
Recent high-field magnetization experiments found the field-polarized state is reached for in-plane magnetic field higher than 10 T in NCTO31. As a result, spins are expected to neither contribute directly, nor indirectly through scattering phonons though pertinent fluctuations at the highest field reached in this study. The saturated κ(T)/T thus represents a background phonon heat conductivity κph in the low-temperature limit with κph/TâââTα, of which the exponent α = 1.5. For ballistic phonons condition, usually reached at the low temperature region, α = 2 should be expected, unless the phonon reflection at the sample surface is not perfectly diffuse. As discussed in Supplementary Note 1, there is no sign of phonon men free path exceeding boundaries scattering limit. Thus a more likely explanation of the the deviation from the α = 2 expectation could be attributed ti non-magnetic and/or non-point-like disorders (spatially extended to mm-length that is comparable to phonon wavelength at low temperature) inside the sample. Obviously, the highest field applied in this work is insufficient to affect on this scattering. Thus, the κ(T)/T curves at highest in-plane field are token as the phonon thermal conductivity κph/T(T) free scattering off magnetic fluctuations. We are aware that α < 2 is quite common among frustrated magnet systems17,32,33, and αâââ1.5 is especially frequently observed.
As can be concluded from Fig. 1 and Supplementary Figure 1, there is no evidence for finite κ0/T at any field. Furthermore, κ/T(T) profiles at zero and low fields never exceed the saturated κph/T(T) over the whole temperature range. Hence, any evidence for a direct magnetic thermal transport channel on top of κph, including κQSL, is elusive. If present, it must be negligibly small since the total measured κ at zero field is already at a very small value and at finite fields the saturated κph is never exceeded. Hence, the intricate landscape of κ(T)/T and its field dependence can be reasonably accounted for by strong phonon-spin scattering, which is the main theme of the rest of the paper. We wish to mention that the above-mentioned zero-field low-temperature hump of κ(T)/T can be a result of scattered phonon transport, but some direct magnetic thermal transport contributions can not be excluded from our results.
In-plane magnetic field dependence of κ/T at fixed temperature
The field dependence of κ/T(T) is non-monotonic, see Fig. 1c. A closer inspection can be performed by isothermal κ/T(B) measurements. Figure 2aâc exhibits the representative data of two samples with three different (âT, B) direction combinations. A straightforward conclusion is that κ(B) is dictated by the B direction, irrelevant to the âT direction. That further underlines the dominance of scattered κph in NCTO over the whole field-temperature phase space. More specifically, a magnetic field applied along the zigzag and armchair directions on the same sample results in sharply different κ(B) profiles. For Bââ¥âzigzag (Fig. 2a), there is strictly no hysteresis behavior. The overall feature turns out to be a broad dip centered around 10.2 T, with many minor anomalies marked by different symbols in the figure. For Bââ¥âarmchair (Fig. 2b, c and Supplementary Figure 2eâk), there are two very sharp dips at 7.5 T and 10.2 T. In contrast to the Bââ¥âzigzag case, a large hysteresis between 3 T and 6.5 T is evident. The abrupt increase of κ in the up-ramp isotherm points to a first-order-like transition. Note that the B direction-dependent occurrence of hysteresis was already mentioned by higher temperature investigations23,34. More discussion on the hysteretic behavior can be found in the Supplementary Note 2. With both B directions, κ saturates at high enough fields, as indicated by the dashed lines.
All the salient features in Fig. 2aâc and Supplementary Figure 2 are summarized in Fig. 2d. Beside the saturation field, which increases roughy linearly with temperature, the other anomalies are basically temperature independent. Such rich phase diagram suggests the landscape of successive B induced phase transitions of the Heisenberg-Kitaev model might be inherited by NCTO35. Remarkably, three different finite-temperature anomalies marked by the yellow, red, and blue shaded areas seem to merge at around 10.2 T at 0 K, suggesting a quantum critical point. We expect the presented result to motivate further theoretical studies that account for specified parameters of NCTO in order to demystify these phases.
Some theoretical exploration on the Heisenberg-Kitaev honeycomb model have already revealed the potentially very rich phase diagram tuned by magnetic field35. Of particular interest is a multi-Q phase that can be stabilized in finite fields35. Recent experiments proposed the zero-field ground state of NCTO could be a triple-Q order36,37,38,39,40, potentially stabilized by the proposed ring exchange interactions38. If confirmed, it could be a natural explanation of the extensive scattering of κph at zero field, since a triple-Q ordered state guarantees larger density of state (DOS) of excitations around zero momentum (Πpoint) which interact strongly with phonons38. These magnetic excitations should in principle contribute to specific heat. A decent experiment to exam this property is still pending for NCTO since heat capacity measurement at very low temperature is technically challenging41. Nevertheless, in Supplementary Note 3 we argue our raw data can provide some indirect evidences for huge specific heat in the zero-field ground state. Our data thus corroborate the notion of very large magnetic DOS as is expected for a triple-Q state.
Next we comment on the oscillation-like κ(B) in the intermediate field range when Bââ¥âarmchair, as highlighted in the inset of Fig. 2c. A similar feature exists but is less obvious in another sample and at higher temperatures, see Supplementary Figure 2g, h, and j. Gapless QSL states that foster quantum oscillation of presumably charge neutral magnetic excitations are indeed anticipated for the Kitaev model4,42,43, and was invoked to interpret the κ(B) oscillations in α-RuCl315. Nevertheless, we believe quantum oscillations of κQSL should be irrelevant to NCTO. First of all, the oscillation feature in NCTO exists in the hysteresis region, thus in a magnetically ordered state. Besides, the fact that there is no evidence for itinerant κQSL at any field does not favor a quantum oscillation picture. Additionally, it is hard to imagine a significant direct contribution of κQSL which is exempted from the strong fluctuations that results in a reduction of κph by orders of magnitude. With this in mind, more caution should be taken for uniqueness of possible Kitaev QSL signatures observed in this material class, since thermal transport properties of NCTO share impressively similar features with α-RuCl323,44,45,46. Indeed, a recent careful investigation on the oscillatory structures of κ(B) of α-RuCl3 indeed addressed all the dips to certain magnetic phase transitions16. We stress that our observation in a similar but simpler sample system further supports this conclusion.
Out-of-plane magnetic field impact on κ
Finally, at low temperatures, when B is applied perpendicular to the NCTO honeycomb plane, κ gets further suppressed beyond a critical field of Bâââ3 T, as shown in Fig. 3a. In the language of magnetic scattering, it means the already strong magnetic fluctuations at zero field are even further promoted, contradictory to the intuition that field alignment of the magnetic moments should reduce them.
Indeed, there are some theoretical considerations predicting Bââ¥âc can lead to abundant phase transitions in the original Kitaev model and its variants35,42,47,48,49. Especially, quantum fluctuations in some of these intermediate phases are expected to be rather strong49, in line with our observations that κph gets further suppressed compared to the zero-field phase. As displayed more clearly in Fig. 3b, κ/T(B) experiences a rapid drop between 2 T and 4 T. It is non-monotonic at higher field, featured by a broad hump around 8 T and a broad dip around 15 T, highlighted by the inset. We note the anomaly at about 15 T can be related to phase transitions revealed by a recent magnetization work31, while the rapid drop between 2 T and 4 T do not have a similar correspondence.
To conclude, our high-quality low-temperature thermal transport study on NCTO single crystals established the dominance of phonon thermal transport in its magnetic ground states, and exclude a significant direction contribution from itinerant magnetic excitations (Please be aware magnon excitation contribution to the transverse thermal transport was detected46). The field evolution of κ, which manifests the strength of phonon-spin scattering, was taken as a probe to elucidate the intricate magnetic phases in NCTO. The strongly scattered phonons in the zero-field ground state indicate large DOS of magnetic excitations, in line with the proposal of a triple-Q state in NCTO. The in-plane κ depends significantly of the in-plane field direction, underlines the bond-dependent interactions of a Kitaev material. However, a criticality at 10.2 T is present regardless of the in-plane field direction. Field perpendicular to the honeycomb plane also induces multiple magnetic phases, where phonon-spin scattering is even more prominent than the zero-field phase. We expect similar physics should be quite common among QSL candidate materials, since QSLs are intrinsically at the edge of competing orders. The pervasive phase boundaries should be careful considered when analyzing existing and future experimental results of QSL candidate materials.
Methods
Sample preparation
Single crystals of Na2Co2TeO6 were prepared with a modified flux method. Na2CO3, Co3O4, and TeO2 powders were grounded and loaded into an alumina crucible. Excess amount of TeO2 served as a self-flux. The crucible was heated up to 1050â°C, kept for two days, before cooling down to 600â°C at a rate of 6.5â°C per hour. Ruby-colored hexagonal flakes of typical size ~10âÃâ10âÃâ0.1 mm3 can be mechanically collected out of bluish violet residue. The harvested single crystals were further washed with a NaOH solution. Basic characterizations of the as-grown crystals can be found in ref. 34.
The samples were cut into a size of about ~5âÃâ1âÃâ0.1 mm3, the longest direction along which âT was generated are the armchair or zigzag directions, respectively. In order to make thermal contacts to the fresh surface, the out-most layers of the as-grown crystals were cleaved off by a blade just before the silver paint was glued to it. As a result, the thickness of the samples were reduced to about 60âμm.
Heat transport measurements
The thermal conductivities were measured in a dilution refrigerator, using a standard four-wire steady-state method with two RuO2 chip thermometers, calibrated in situ against a reference RuO2 thermometer. For fixed field κ(T) measurements, each sample was initially cooled down from room temperature without magnetic field. The data were collected with gradually increased magnetic field in order to avoid the complicated situation of involving the hysteresis effect. κ(T) data were collected in a steady-state manner (see Supplementary Figure 3).
To measure the κ(B) isotherms, the system was cooled from the paramagnetic phase in zero field. Then the sample temperature was kept at a set point, and the field was changed from 0 T to the highest field (up-ramp), followed by decreasing to 0 T (down-ramp). The field was changed at a speed no more than 20 mT per minute, in order to minimize its heating effect, and to keep the system in a (quasi-)thermal equilibrium state. In some cases, the field were ramped up and down multiple times, see Supplementary Note 2 for more details.
We wish to point out that according to our experience, NCTO seems to have very strong anisotropic in-plane magnetic susceptibility at very low temperatures, that tend to align the armchair direction to the field. We are aware of a recently posted paper which reports a direct measurement of this magnetic anisotropy that underpins our conjecture50. In order to avoid experimental artifacts caused by such effect, the cooler side of all samples were glued to the heat sink (a gilded silver bulk) directly with a strong epoxy (Wakefield DeltaBond-152 two component adhesives). The samples and the wires attached to them were checked carefully under a microscope after the experiments, to confirm they had not bent during the measurements.
Thermodynamic measurements
The magnetization measurements were carried out using a SQUID magnetometer (iHelium3, MPMS-XL, Quantum Design). The magnetostriction was measured using a commercially available dilatometer (Standard probe, Kuechler) compatible with the Quantum Design PPMS systems. These measurements are based on the capacitance measurement technique (AH2700A, Andeen-Hagerling).
Data availability
The data that support the findings of this study are available from the corresponding authors upon reasonable request.
References
Balents, L. Spin liquids in frustrated magnets. Nature 464, 199â208 (2010).
Savary, L. & Balents, L. Quantum spin liquids: a review. Rep. Prog. Phys. 80, 106502 (2017).
Zhou, Y., Kanoda, K. & Ng, T. K. Quantum spin liquid states. Rev. Mod. Phys. 89, 025003 (2017).
Kitaev, A. Anyons in an exactly solved model and beyond. Ann. Phys. 321, 2â111 (2006).
Moessner, R. & Sondhi, S. L. Resonating valence bond phase in the triangular lattice quantum dimer model. Phys. Rev. Lett. 86, 1881 (2001).
Wen, X.-G. Quantum orders in an exact soluble model. Phys. Rev. Lett. 90, 016803 (2003).
Yao, H., Zhang, S.-C. & Kivelson, S. A. Algebraic spin liquid in an exactly solvable spin model. Phys. Rev. Lett. 102, 217202 (2009).
Ben-Zion, D., Das, D. & McGreevy, J. Exactly solvable models of spin liquids with spinons, and of three-dimensional topological paramagnets. Phys. Rev. B 93, 155147 (2016).
Takagi, H., Takayama, T., Jackeli, G., Khaliullin, G. & Nagler, S. E. Concept and realization of Kitaev quantum spin liquids. Nat. Rev. Phys. 1, 264â280 (2019).
Broholm, C. et al. Quantum spin liquids. Science 367, 263 (2020).
Banerjee, A. et al. Proximate Kitaev quantum spin liquid behavior in a honeycomb magnet. Nat. Mat. 15, 733â740 (2016).
Kasahara, Y. et al. Majorana quantization and half-integer thermal quantum Hall effect in a Kitaev spin liquid. Nature 559, 227â231 (2018).
Lefrançois, Ã. et al. Oscillations in the magnetothermal conductivity of α-RuCl3: evidence of transition anomalies. Phys. Rev. B 107, 064408 (2023).
Czajka, P. et al. Planar thermal Hall effect of topological bosons in the Kitaev magnet α-RuCl3. Nat. Mater. 22, 36â41 (2023).
Czajka, P. et al. Oscillations of the thermal conductivity in the spin-liquid state of α-RuCl3. Nat. Phys. 17, 915â919 (2021).
Bruin, J. A. N. et al. Origin of oscillatory structures in the magnetothermal conductivity of the putative Kitaev magnet α-RuCl3. APL Mater. 10, 090703 (2022).
Yu, Y. J. et al. Ultralow-temperature thermal conductivity of the Kitaev honeycomb magnet α-RuCl3 across the field-induced phase transition. Phys. Rev. Lett. 120, 067202 (2018).
Zhang, H. et al. Stacking disorder and thermal transport properties of α-RuCl3. Phys. Rev. Materials 8, 014402 (2024).
Breitner, F. A., Jesche, A., Tsurkan, A. & Gegenwart, P. Thermal decomposition of the Kitaev material α-RuCl3 and its influence on low-temperature behavior. Phys. Rev. B 108, 045103 (2023).
Kasahara, Y. et al. Quantized and unquantized thermal Hall conductance of the Kitaev spin liquid candidate α-RuCl3. Phys. Rev. B 106, L060410 (2022).
Liu, H. & Khaliullin, G. Pseudospin exchange interactions in d7 cobalt compounds: possible realization of the Kitaev model. Phys. Rev. B 97, 014407 (2018).
Lin, G. et al. Field-induced quantum spin disordered state in spin-1/2 honeycomb magnet Na2Co2TeO6. Nat. Commun. 12, 5559 (2021).
Hong, X. et al. Strongly scattered phonon heat transport of the candidate Kitaev material Na2Co2TeO6. Phys. Rev. B 104, 144426 (2021).
Yamashita, M., Shibauchi, T. & Matsuda, Y. Thermal-transport studies on two-dimensional quantum spin liquids. Chemphyschem 13, 74â78 (2012).
Yamashita, M. et al. Highly mobile gapless excitations in a two-dimensional candidate quantum spin liquid. Science 328, 1246â1248 (2010).
Hartstein, M. et al. Fermi surface in the absence of a Fermi liquid in the Kondo insulator SmB6. Nat. Phys. 14, 166â172 (2018).
Li, N. et al. Possible itinerant excitations and quantum spin state transitions in the effective spin-1/2 triangular-lattice antiferromagnet Na2BaCo(PO4)2. Nat. Commun. 11, 4216 (2020).
Sato, Y. et al. Charge-neutral fermions and magnetic field-driven instability in insulating YbIr3Si7. Nat. Commun. 13, 394 (2022).
Pan, B. et al. Unambiguous experimental verification of linear-in-temperature spinon thermal conductivity in an antiferromagnetic Heisenberg chain. Phys. Rev. Lett. 129, 167201 (2022).
Guang, S. K. et al. Thermal transport of fractionalized antiferromagnetic and field-induced states in the Kitaev material Na2Co2TeO6. Phys. Rev. B 107, 184423 (2023).
Zhang, S. et al. Electronic and magnetic phase diagrams of the Kitaev quantum spin liquid candidate Na2Co2TeO6. Phys. Rev. B 108, 064421 (2023).
Barthélemy, Q. et al. Heat conduction in herbertsmithite: field dependence at the onset of the quantum spin liquid regime. Phys. Rev. B 107, 054434 (2023).
Fauqué, B. et al. Thermal conductivity across the metal-insulator transition in the single-crystalline hyperkagome antiferromagnet Na3+xIr3O8. Phys. Rev. B 91, 075129 (2015).
Yao, W. & Li, Y. Ferrimagnetism and anisotropic phase tunability by magnetic fields in Na2Co2TeO6. Phys. Rev. B 101, 085120 (2020).
Janssen, L., Andrade, E. C. & Vojta, M. Honeycomb-lattice Heisenberg-Kitaev model in a magnetic field: spin canting, metamagnetism, and vortex crystals. Phys. Rev. Lett. 117, 277202 (2016).
Chen, W. et al. Spin-orbit phase behavior of Na2Co2TeO6 at low temperatures. Phys. Rev. B 103, L180404 (2021).
Lee, C. H. et al. Multistage development of anisotropic magnetic correlations in the Co-based honeycomb lattice Na2Co2TeO6. Phys. Rev. B 103, 214447 (2021).
Krüger, W. G. F., Chen, W., Jin, X., Li, Y. & Janssen, L. Triple-Q order in Na2Co2TeO6 from proximity to hidden-SU(2)-symmetric point. Phys. Rev. Lett. 131, 146702 (2023).
Yao, W. et al. Magnetic ground state of the Kitaev Na2Co2TeO6 spin liquid candidate. Phys. Rev. Research 5, L022045 (2023).
Kikuchi, J. et al. Field evolution of magnetic phases and spin dynamics in the honeycomb lattice magnet Na2Co2TeO6: 23Na NMR study. Phys. Rev. B 106, 224416 (2022).
Wilhelm, H., Lühmann, T., Rus, T. & Steglich, F. A compensated heat-pulse calorimeter for low temperatures. Rev. Sci. Instrum. 75, 2700â270 (2004).
Gohlke, M., Moessner, R. & Pollmann, F. Dynamical and topological properties of the Kitaev model in a [111] magnetic field. Phys. Rev. B 98, 014418 (2018).
Hickey, C. & Trebst, S. Emergence of a field-driven U(1) spin liquid in the Kitaev honeycomb model. Nat. Commun. 10, 530 (2019).
Hentrich, R. et al. Unusual phonon heat transport in α-RuCl3: strong spin-phonon scattering and field-induced spin gap. Phys. Rev. Lett. 120, 117204 (2018).
Hentrich, R. et al. Large thermal Hall effect in α-RuCl3: evidence for heat transport by Kitaev-Heisenberg paramagnons. Phys. Rev. B 99, 085136 (2019).
Gillig, M. et al. Phononic-magnetic dichotomy of the thermal Hall effect in the Kitaev material Na2Co2TeO6. Phys. Rev. Res. 5, 043110 (2023).
Sørensen, E. S., Catuneanu, A., Gordon, J. S. & Kee, H.-Y. Heart of entanglement: chiral, nematic, and incommensurate phases in the Kitaev-Gamma ladder in a field. Phys. Rev. X 11, 011013 (2021).
Li, H. et al. Identification of magnetic interactions and high-field quantum spin liquid in α-RuCl3. Nat. Commun. 12, 4007 (2021).
Chern, L. E., Kaneko, R., Lee, H.-Y. & Kim, Y. B. Magnetic field induced competing phases in spin-orbital entangled Kitaev magnets. Phys. Rev. Research 2, 013014 (2020).
Lin, G. et al. Evidence for field induced quantum spin liquid behavior in a spin-1/2 honeycomb magnet. Preprint at https://doi.org/10.21203/rs.3.rs-2034295/v1 (2022).
Acknowledgements
We thank Jeroen van den Brink, Matthias Vojta, Vladislav Kataev, Christoph Wellm, Satoshi Nishimoto, Wenjie Chen, Yang Xu, and Wilhelm Krüger for illuminating discussions and collaboration on related work. We would like to thank Nicolás Pérez for generously sharing the equipment, and thank Danny Baumann and Tino Schreiner for their technical assistance, especially during the pandemic lab-shutdown crisis. This work has been supported by the Deutsche Forschungsgemeinschaft (DFG) through SFB 1143 (Project-ID No.247310070) and the Würzburg-Dresden Cluster of Excellence ct. qmat (EXC 2147, Project-ID No. 390858490). C.H. has received funding from the European Research Council (ERC) under the European Unionâs-Horizon 2020 research and innovation program (Grant Agreement No. 647276-MARS-ERC-2014-CoG). L.J. was supported by the DFG through the Emmy Noether program (JA2306/4-1, Project No. 411750675). V.K. was supported by the Alexander von Humboldt Foundation. Y.L. acknowledges funding support from the National Basic Research Program of China (Grant No. 2021YFA1401901) and the NSF of China (Grant No. 12061131004).
Funding
Open Access funding enabled and organized by Projekt DEAL.
Author information
Authors and Affiliations
Contributions
B.B., C.H., and Y.L. conceived and initiated the study. W.Y. and Y.L. grown the crystals. X.H. and M.G. performed the transport measurements. V.K., S.G., and A.U.B.W. performed the thermodynamic measurements. X.H., C.H., and L.J. analyzed the data. X.H., C.H., L.J., and V.K. wrote the manuscript with input from all authors. C.H. supervised the project.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
41535_2024_628_MOESM1_ESM.pdf
Supplementary Materials for: Phonon thermal transport shaped by strong spin-phonon scattering in a Kitaev material Na2Co2TeO6
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Hong, X., Gillig, M., Yao, W. et al. Phonon thermal transport shaped by strong spin-phonon scattering in a Kitaev material Na2Co2TeO6. npj Quantum Mater. 9, 18 (2024). https://doi.org/10.1038/s41535-024-00628-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41535-024-00628-4