## Scaling behavior and giant field-enhancement of the thermal conductivity in the honeycomb antiferromagnet BaCo (AsO ) 2 4 2

Jiayi Hu , Ruidan Zhong 1 2,* , Peter Czajka , Tong Gao 1 1,** , R. J Cava , and N. P. Ong 2 1 1 Department of Physics, Princeton University, Princeton, NJ, 08544, USA 2 Department of Chemistry, Princeton University, Princeton, NJ, 08544, USA

(Dated: October 28, 2024)

The layered honeycomb material BaCo (AsO ) 2 4 2 (BCAO) is of topical interest because its magnetic state is related to that of the Kitaev magnet α -RuCl . 3 Using thermal transport to probe how magnetic excitations interact with phonons in the magnetically disordered regime, we have uncovered an unusually large enhancement of the thermal conductivity κ xx in an in-plane magnetic field H . Just above the N´el temperature e T N , a field of 13 T increases κ xx by a factor ∼ 211, much larger than reported previously in any magnetic insulator. Interestingly, κ xx ( H,T ) exhibits a scaling behavior in the entire magnetically disordered region that surrounds the ordered zigzag state. The ratio ∆ κ xx ( H,T /κ ) xx (13 , T ), measured throughout the disordered region, collapses to a one-parameter scaling function exp( -1 /gx ) (where x = µ B B/k T B and g is a constant).

## I. INTRODUCTION

The family of layered honeycomb spin1 2 magnetic materials has attracted considerable interest. In particular, the ruthenate α -RuCl 3 with local moments derived from 4 d orbitals in Ru has been identified [1, 2] as the closest proximate to the Kitaev Hamiltonian H K [3]. The lowtemperature thermal transport properties have played a central role in illuminating the extent to which predictions based on H K agree with experiment [4-8]. Interest has extended to related layered, honeycomb spin1 2 materials. In the cobaltate BaCo (AsO ) 2 4 2 (BCAO), the local moments based on Co-3 d orbitals undergo a sharp transition to a zig-zag state at the N´el temperature e T N = 5.4 K [9-11]. In its magnetic, thermal and crystallographic properties, BCAO shares many features with α -RuCl . A 3 number of recent reports have compared their magnetic and transport properties [8, 12-20].

regime which covers a large fraction of the H T -plane. We report an unusual scaling behavior of the thermal conductivity κ xx ( H,T ) in BCAO observed throughout the magnetically disordered region below 10 K. In zero H , strong inelastic scattering of phonons from spin excitations severely suppresses phonon conduction. In an in-plane magnetic field H , κ xx ( H,T ) increases dramatically following a one-parameter scaling function. Just above T N , the field-induced enhancement factor ( &gt; 200) exceeds that previously reported in all magnetic insulators to our knowledge.

In the case of BCAO, the theoretical and experimental investigations appear to favor a J 1 -J 3 , XXZ Hamiltonian instead of the Kitaev Hamiltonian [13-18]. Below T N , the spins in BCAO exhibit a helical magnetic texture in the ab -plane. The spins, nominally aligned along the b -axis, exhibit a spiral structure with axis parallel to a ∗ . In an in-plane magnetic field H , the spins are polarized in a ferrimagnetic, colinear pattern above a field H c1 ∼ 0 2 . T. Long-range magnetic order is suppressed at the critical field H c2 = 0 5 T [9][11][13]. . The transition at the lower field H c1 has been investigated by magnetization and terahertz spectroscopy experiments [12, 13, 19]. A recent study combining angular dependence of magnetization and dilatometry suggests an intermediate phase appears before the spins are fully polarized [18]. Another thermal transport experiment measured down to subKelvin temperatures supports the KQSL picture, yet other theoretical proposal attribute it to a Dirac spin liquid phase [19, 21].

Our work here focusses on the magnetically disordered

Underlying the giant field-enhancement of κ xx is the intense scattering of phonons by the spins which leads to 'localization' of the phonons in zero H . Our results suggest that the non-analytic behavior of κ xx ( H,T ) in the limit H → 0, as well as the scaling behavior over a large fraction of the magnetically disordered region, may well be signature hallmarks of strong spin-phonon coupling.

In several topical areas, e.g. spin qubits based on NV centers in diamond [22], spin-phonon coupling is a key topic of interest. The scaling behavior uncovered here suggests an incisive way to probe these couplings in magnetic materials.

## II. EXPERIMENT

To measure the thermal conductivity tensor κ ij of the BCAO crystal, we heat sunk 3 RuO 2 thermometers (matched and calibrated) to its edges using stycast and gold wires (Fig. 1a). One short edge of the crystal is heat sunk to the bath (brass substrate) while the opposite edge is attached to a heater via gold wires. The sample space is kept at a pressure less than 5 × 10 -5 mbar. We label as ˆ x the common axis of the heat current density J and field H .

Typically, we measure the field profile of κ xx ( H,T )

FIG. 1: Panel (a): Schematic of the crystal attached to the heater, 3 thermometers and the bath by gold wires. The applied H direction is indicated. Panel (b): The two-up, two-down zigzag phase inferred from neutron scattering [11]. Panel (c): Comparison of the zeroH κ T ( ) (square blue symbols, multiplied by a factor of 30) with κ xx vs. T measured in an in-plane 13-T magnetic field (red symbols). The faint blue curve passing through the blue symbols is an interpolation supplemented by measurements in zero H in which T is continuously varied. The red curve is a fit to the 13-T data using the standard Callaway expression. More details can be found in the Supplement [23]. Panel (d): The inferred phase diagram of BCAO in the H T -plane. The magnetically ordered zigzag state is shaded red and coral. Throughout the disordered region (shaded green), κ xx ( H,T ) obeys a one-parameter scaling behavior. Deviation from this scaling occurs in the blue region. The fields H c1 and H c2 are obtained from data presented in Figure 2a.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00873v2-figure-79f4398f-11cd-4af5-88e6-eb0394a91589.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00873v2-figure-b0600a88-d8f4-432d-8025-a70fa60289f8.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00873v2-figure-36512e5a-5517-4592-95f1-d4f8b10141df.jpg)

keeping the bath temperature T bath regulated at a set point as H is swept from 0 to 13 T and back. At low T , the applied heat current invariably leads to a finite shift δT = T -T bath between the sample's temperature T and T bath (the shift is usually nearly H independent). Here, however, the anomalously large (211-fold) increase in κ xx ( H,T ) causes δT H,T ( ) to be strongly field dependent. In all field profiles reported here, we have recorded δT H,T ( ) to correct for these shifts. We also adopted steps to avoid errors arising from large magnetocaloric effects which lead to long relaxation times for both the sample and bath to reach equilibrium whenever H undergoes a step-change (Secs. 1, 3 and 4 of Supplement [23]). Finally, the tendency of the crystal to fracture under strain exerted by magnetic torque restricts measurements so far to in-plane H .

( ≡ κ xx (0 , T ), plotted as blue symbols in Fig 1c) slowly decreases to an anomalously low value ( ∼ 0 6 W/mK) as . T is lowered to 5.4 K. At the N´el temperature e T N (5.4 K), the spins undergo a phase transition to an antiferromagnetic state with zig-zag order. In the ordered state, κ rises to a broad peak that reflects contributions from phonons and magnons in the ordered state.

Motivated by the observation of an unusual planar thermal Hall conductivity κ xy in the Kitaev magnet α -RuCl 3 [6, 7], we have searched for a planar thermal Hall signal in BCAO focusing on the disordered region close to the ordered phase. Despite a careful search we did not detect a finite κ xy signal in BCAO within our resolution (See Sec. 7 of Supplement [23]).

## III. FIELD AND TEMPERATURE PROFILES OF THERMAL CONDUCTIVITY

The thermal conductivity κ xx ( H,T ) of BCAO displays a rich set of features as a function of temperature T and an in-plane magnetic field H . In the absence of H , κ T ( )

In our experiment, we observe that long-range order is abruptly suppressed by an in-plane field H ∼ 0.5 T, consistent with previous studies [12, 13, 15, 20]. As described below, phonons remain strongly coupled to spin excitations throughout the disordered region surrounding the ordered state (green region in Fig. 1d). The strong coupling is steadily weakened as H increases. Eventually, at large H (13 T), scattering of phonons by the spins becomes negligible, and κ xx attains large values (red symbols in Panel c). The field enhancement, expressed by the ratio κ xx ( H,T /κ T ) ( ) ( ∼ 211 just above T N ) exceeds that in all other magnetic insulator reported to date. Our key result is that, over the entire magnetically disordered region (green region in Panel d), κ xx ( H,T ) obeys a oneparameter scaling function.

Figure 2a displays field profiles of κ xx within the ordered state with T fixed at the values indicated. At the lowest T (0.49 K), the profile of κ xx is virtually flat. When T is raised (1.64 to 2.75 K) the field profiles feature a dome that is abruptly truncated on the largeH side by a vertical cliff as H approaches H c2 (0.5 T), and long-range order is suppressed. On the low H side, the dome displays a notch that corresponds to the lower critical field H c1 ( ∼ 0 12 T). In previous reports, . H c1 was in-

FIG. 2: Profiles of κ xx vs. field and temperature. Panel (a): Field profiles of κ xx ( H,T ) at low H µ H &lt; ( 0 1 2 T). In each . curve, T is fixed at the values indicated. Curves in the temperature interval 1 6 . &lt; T &lt; 2 8 K display a vertical cliff that . signals the abrupt collapse of long-range magnetic order as H approaches H c2 = 0.5 T. At low H , a second sharp feature H c1 ( ∼ 0.1 T) is observed. Panel (b): Field profiles of κ xx extended to 11 T, with T fixed at the values indicated. At each T , κ xx initially increases exponentially with H . At large H , it gradually saturates to its high-field limit. Profiles shown in Panel (a) are visible as the small mesa below 0.5 T. Panel (c): Curves of κ xx ( H,T ) plotted versus T with H fixed at the values indicated. The curve κ xx (13 , T ) at the maximum applied H (13 T) serves as an envelope that caps the maximum value attained at each temperature T . Each curve (e.g. the one at H = 4 T) initially tracks closely κ xx (13 , T ) as T is raised above 1 K, but deviates away from κ xx (13 , T ) as T is raised (at ∼ 2 5 K). Curves taken at larger . H (e.g. 10 T) deviate from κ xx (13 , T ) at a higher temperature (4 K). This pattern reflects the one-parameter scaling behavior discussed below. Note the anomalously low values of κ xx in the curve at 0.6 T. Panel (d): First hints of scaling behavior above T N . When we plot the field-induced change ∆ κ xx ( H,T ) vs. the scaling variable µ H/T 0 , we find that all data points above T N converge to the same curve when the scaling variable µ H/T 0 falls below ∼ 0 4 T/K. The temperature of each data point is given by the color scale. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00873v2-figure-bdaaea08-be2a-45e5-a16c-4afbafde005f.jpg)

terpreted as a transition between either an incommensurate spiral state or a double zigzag pattern to the cAFM state [11, 20]. Unlike in Refs. [12, 19], we observe only a single sharp transition within the field interval 0.1-0.2 T. In the disordered phase above H c2 , κ xx increases steeply with H . This is a small fraction of the scaling pattern uncovered below.

axis extended to 11 T (the profiles in the ordered phase are now reduced to small mesas visible below 0.5 T). At each T (e.g. 3 K), κ xx increases from a small value ( ∼ 1 W/mK) before slowly saturating to a maximum value ( ∼ 42 W/mK) as H increases to 13 T. All the profiles appear to follow the same functional form, with a pattern amenable to simple analysis.

Figure 2b shows the field profiles of κ xx with the H

To this end, we turn to the orthogonal cuts in which

κ xx ( H,T ) is plotted vs. T with H fixed (Fig. 2c). At each fixed value of H , the curve of κ xx vs. T is now non-monotonic, exhibiting a broad maximum at the temperature T p ( H ). As shown in the inset, T p ( H ) increases roughly linearly with H . This suggests that, in the disordered region, κ xx depends on H solely through the scaling variable

$$x = \mu _ { \text{B} } B / k _ { \text{B} } T, \text{ \quad \ \ } ( 1 ) \text{ \quad \ } \text{pl}.$$

where µ B and k B are, respectively, the Bohr magneton and the Boltzmann constant, and B = µ H 0 is the magentic flux density with µ 0 the vacuum permeability.

To test this preliminary hint of scaling behavior, we plot in Fig. 2d the field-induced change ∆ κ xx ( H,T ) = κ xx ( H,T ) -κ xx (0 , T ) versus x in log-log scale for T between 5.4 and 9.5 K. We find that the data merge to define a scaling curve that extends over nearly 3 decades in ∆ κ xx ( H,T ).

Another important feature in Fig. 2c is the role of κ xx (13 , T ) (black squares) relative to κ xx measured at lower fields throughout the disordered region. Combining the perspectives afforded by the fixedT and fixedH cuts, we infer that when T is fixed (e.g. at 4 K), κ xx increases towards κ xx (13 , T ) but does not exceed it.

Physically, this implies that, in the high-field limit, phonon conduction asymptotes to the regime in which ℓ , restricted by scattering from crystalline disorder only, becomes independent of T . In this limit, the T dependence observed in κ xx ( H,T ) = 1 3 c gs ( T v ℓ ) s arises from the heat capacity c gs ( T ) of the spin-phonon system ( v s is the group velocity of the acoustic branch). As inferred from the field profiles in Fig. 2b, κ xx (13 , T ) is a good approximation to this high-field limit below T = 6 K. The approximation is less satisfactory above 6 K as saturation is approached at increasingly large H . The dashed curve in Fig. 2c is a fit to κ xx (13 , T ) using the Callaway model in this limit (see Sec. 2 of Supplement [23]).

## IV. SCALING BEHAVIOR

The functional form of the scaling curve in Fig. 2d becomes evident when we plot ∆ κ xx vs. 1 /x in semilog plot. As shown in Fig. 3a, all the data above 5.4 K collapse onto the straight line corresponding to the function ∆ κ xx ∼ e -1 /gx with g a constant.

The simple scaling based on x accounts only for the H dependence of ∆ κ xx ( T, H ). The further dependence on T requires more care. This becomes apparent if we attempt to apply the same analysis to the data below 5.4 K in the disordered region ( H &gt; H c2 ).

The semilog plot of κ xx vs. 1 /x shown in Fig. 3b reveals that the data segregate into two distinct groups. Data in the ordered phase (not relevant to scaling) follow flat trajectories that are nearly independent of x . By contrast, the data in the disordered region fall along diagonal curves that mimic the scaling function except that they are now displaced vertically by a factor of the form

T β with exponent 2 &lt; β &lt; 3, which we infer arises from the heat capacity c gs ( T ). From the discussion on Fig. 2c, we anticipate that the offsets can be cancelled if we divide κ xx ( H,T ) by κ xx (13 , T ).

This inference is confirmed when we replot κ xx ( H,T /κ ) xx (13 , T ) vs. 1 /x in Fig. 3c (upper panel). Now all data throughout the disordered region below 5.4 K collapse onto one scaling curve. A similar plot shows that data above 5.4 K also collapse onto the same curve (lower panel). As noted above, deviation from scaling appears when x exceeds ∼ 0 3. .

Hence, throughout the disordered region below and above 5.4 K, the ratio ∆ κ xx ( H,T /κ ) xx (13 , T ) is described by the one-parameter scaling form

$$\frac { \Delta \kappa _ { x x } ( H, T ) } { \kappa _ { x x } ( 1 3, T ) } = A \exp [ - \frac { 1 } { g x } ] \, \quad \quad ( 2 )$$

where x is given by Eq. 1 and A and g are constants.

The function e -1 /gx implies that ∆ κ xx initially increases from zero at a very slow rate (the non-analytic limit H → 0 requires all derivatives to vanish). As discussed, the ultra-slow initial growth may be seen in Fig. 2a. However, once B rises above 1 Tesla, ∆ κ xx undergoes an exponential increase before ultimately saturating to a very large value consistent with phonons in a crystal with a low density of lattice and spin disorder.

We sketch the physical picture suggested by our results. The single-parameter scaling form constrains the candidate scenarios for describing the field-enhanced κ xx .

First, the anomalously low value of κ T ( ) in the fluctuating interval 5 4 . &lt; T &lt; 10 K at H = 0 implies strong spin-phonon coupling which leads to repeated absorption or emission of phonons accompanied by spin flips. Secondly, the non-analytic nature in Eq. 2 in the limit x → 0 implies that, in weak H , deviation of ∆ κ xx ( H,T ) from zero is initially nearly unresolvable until H exceeds a field scale of ∼ 1 T. Then κ xx increases exponentially. This strongly non-linear increase is not amenable to a perturbative treatment. (We also point out that Eq. 2 is very different from the familiar expression e -∆ b /k B T describing thermal activation of charge carriers across a band gap ∆ .) b

Throughout the disordered region, we infer that the phonons and spin excitations constitute a tightly coupled system in which ℓ 0 (the phonon mean-free path at H = 0) is limited to less than 30 nm by intense scattering. Although the scattering may be elastic or inelastic, the scaling form implies that inelastic processes are dominant in the disordered region. By analogy with experiments on the honeycomb system α -RuCl 3 [24-27] we may expect that spin excitations in the disordered state exhibit a broad, featureless spectrum that is gapless in zero H . The low-lying spin excitations interact strongly with phonons at a select population of sites, which absorb and emit phonons. At temperature T , the phonon energy distribution peaks sharply at ℏ ω ≃ k B T .

An applied Zeeman field opens a gap ∆ = gµ B H in the spin spectrum where g is the effective gyromagnetic

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00873v2-figure-bcabee15-2129-4cb2-9a31-8018cd6236f5.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00873v2-figure-b323bae9-540f-4739-b7d1-04367e1994b3.jpg)

FIG. 3: Uncovering the scaling behavior. (Panel a) The analytical form of the scaling behavior implied in Fig. 2d becomes apparent in a semilog plot of ∆ κ xx vs. 1 /x with x defined in Eq. 1. All data points taken above T N fall on the dashed line given by Eqs. 2 and 1. The equation of the fit (dashed line) is ln ∆ κ xx = -0 628 . x -1 +5 476, which gives . g = 1 6. . The scatter as x → 0 reflects uncertainties in ∆ κ xx . (Panel b) Finding the missing pre-factor in the scaling form. At temperatures below 5.4 K, the data points κ xx ( H,T ) segregate into two distinct groups when plotted vs. 1 /x . In the ordered phase ( H &lt; H c2 ), the data follow nearly flat trajectories. By contrast, data in the disordered region ( H &gt; H c2 ) fall on diagonal curves that match the dashed line in Panel (a), except that they are systematically displaced downwards as T → 0. This suggests that κ xx should be divided by a T -dependent factor T β , which we approximate by the measured curve κ xx (13 , T ). (Panel c): Semilog plots of κ xx /κ xx (13 , T ) vs. 1 /x for T &lt; 5.4 K (upper panel) and T &gt; 5.4 K (lower). Once κ xx is divided by κ xx (13 , T ), the diagonal curves in Panel (b) collapse to a single curve, which represents the universal curve given in Eq. 2. The data above 5.4 K also collapse onto the same curve (lower panel). Hence, over the disordered region (green region in Fig. 1d), scaling behavior prevails. Departure from scaling becomes observable when 1 /x decreases below ∼ 3 (dark blue region in Fig. 1d).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00873v2-figure-72f0ee2b-5b19-453e-ba2b-ba44f882a814.jpg)

parameter for the low-lying spin excitations. Phonon annihilation is blocked when the gap exceeds k B T , i.e.

$$g \mu _ { \text{B} } B > k _ { \text{B} } T. \quad \quad \ ( 3 ) \quad 2 \text{a}$$

In effect, the identified absorption sites act as escape hatches that are propped open by H . Rendering these sites transparent to phonons allows them to escape the trapping well. Despite their small population, these sites play an outsize role in determining the initial growth of κ xx ( H,T ).

hatches far exceeds ℓ 0 . Each liberated phonon does not travel far before becoming retrapped. Hence the initial increase in κ xx is nearly unresolvable, as shown in Figures 2a and 2b, and implied by the non-analyticity of e -T/B in the limit B → 0. Crucially, by Eq. 3, the field at which phonons escape at the identified sites is proportional to T . Hence, the field-induced increase in ℓ automatically satisfies the scaling in Eq. 2.

When H is weak, the average spacing between escape

Increasing H weakens the amplitude of the fluctuating magnetic order by increasing the polarization of spins, which causes the population of escape hatches to pro-

liferate. In analogy with percolation conductance, the largest enhancement of ℓ occurs when connected paths form skeletal spines that span the crystal. When H exceeds ∼ 1 T, rapid proliferation of the hatches enables the connectivity of phonon propagation paths to exceed the percolation limit, and κ xx undergoes the exponential growth observed.

## V. COMPARISON WITH OTHER MAGNETIC INSULATORS

Experiments showing large field-induced enhancement of κ xx have been reported in several layered magnetic insulators. We briefly survey the results beginning with the largest previously reported (Fig. 4 and Table S1 in Supplement [23]). Zhao et al. found that, in Co V O , 3 2 8 κ xx increases 100 fold in a field of 14 T [28] at 7 5 . K . They interpret the increase as caused by a series of magnetic phase transitions occuring in a narrow window (5 transitions between 6 and 12 K). While there are fewer known phase transitions in the case of BCAO, recent dilatometry experiment indeed suggests more intermediate phases [18]. One difference we want to point out though is that BCAO has only one critical temperature, and it has been relatively robust even when applied pressure [15]. More importantly, the field-induced enhancement persists up to as high as 13T well beyond the critical field H c 2 = 0 5 . T , and the scaling behavior also covers a wide regime in the phase diagram surrounding the ordering states. As a result, a succession of phase transitions doesn't seem to be the whole story, though intense spin fluctuations do seem to be one of the key components.

A recent experiment on another honeycomb Kitaev magnet candidate Na Co TeO 2 2 6 found that κ xx increases by a factor of 70 in a field of 16 T [29]. The authors attribute this intense scattering to the abundance of magnetic excitations above a triple-Q ground state. The authors do not report any scaling behavior in the observed κ xx vs. H .

In the Kitaev magnet α -RuCl , 3 κ xx has been reported to increase 13-fold compared to the value at 7 T , the critical field of its zig-zag state [30]. The ratio is smaller if compared to its zero field value (See Table 1 in the Supplement [23]). The authors are able to fit κ ab with a modified Callaway model, adding a relaxation term that takes into account the resonance between phonon and magnetic gap, whose size is linear to the applied field. Motivated by this study, We also tried fitting our data at finite field with the modified Callaway model, with spin-phonon resonance terms. However, we have not succeeded in fitting our data with physically meaningful parameters (the exception is the curve at 13 T where ℓ is limited by disorder scattering). In fact, the non-analyticity inherent in the H → 0 limit in Eq.3 rules out a Callaway-type description in the scaling regime in BCAO. Because it is not possible to reproduce the non-analyticity in a Callowaytype formulation, our conclusion is that the empirical fits

FIG. 4: Comparison of the field enhancement of κ xx (expressed as a ratio) with the ratio reported for other magnetic insulators[28-34]. BCAO is marked with a red star. In the ratio for HoMnO ∗ 3 the denominator is κ xx (2.2T). For RuCl ∗∗ 3 , the denominator is κ xx (7.5T). See Sec. 6 in Supplement [23] for detailed information.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00873v2-figure-0f8a6349-0c2b-477b-9294-b1e5f33fa686.jpg)

used in α -RuCl 3 cannot be applied to BCAO.

***************

- * Present address of R.Z. : Tsung-Dao Lee Institute and School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 201210, China
- ** Present address of T.G. : Building 3, West Area, Xibeiwang East Street, Haidian District, Beijing 100193, China

Acknowledgement The thermal transport experiments were supported by the US Department of Energy (Basic Energy Sciences) under contract DE-SC0017863 to N.P.O. An award from the U.S. National Science Foundation (DMR 2011750) supported the materials growth (R.J.C.) and acquisition of cryogenic equipment (N.P.O.). N.P.O. acknowledges support from the Gordon and Betty Moore Foundation EPiQS initiative (grant GBMF9466).

## Author contributions

J.H., P.C. and N.P.O. conceptualized the experiment. J.H. performed all measurements with initial guidance from P.C. Data in an earlier experiment by T.G. provided substantial insights. The present data were analyzed by

- [1] Banerjee, A., Bridges, C.A., Yan, J.Q., Aczel, A.A., Li, L., Stone, M.B., Granroth, G.E., Lumsden, M.D., Yiu, Y., Knolle, J. and Bhattacharjee, S., Proximate Kitaev quantum spin liquid behaviour in a honeycomb magnet. Nature Materials . 15 , 733-740 (2016).
- [2] Jackeli, G. and Khaliullin, G. Mott Insulators in the Strong Spin-Orbit Coupling Limit: From Heisenberg to a Quantum Compass and Kitaev Models. Phys. Rev. Lett. . 102 , 017205 (2009,1).
- [3] Kitaev, A. Anyons in an exactly solved model and beyond. Annals Of Physics . 321 , 2-111 (2006), January Special Issue.
- [4] Leahy, I., Pocs, C., Siegfried, P., Graf, D., Do, S., Choi, K., Normand, B. and Lee, M. Anomalous thermal conductivity and magnetic torque response in the honeycomb magnet α -RuCl . 3 Physical Review Letters . 118 , 187203 (2017).
- [5] Czajka, P., Gao, T., Hirschberger, M., Lampen-Kelley, P., Banerjee, A., Yan, J., Mandrus, D., Nagler, S. and Ong, N. Oscillations of the thermal conductivity in the spin-liquid state of α -RuCl . 3 Nature Physics . 17 , 915-919 (2021).
- [6] Kasahara, Y., Ohnishi, T., Mizukami, Y., Tanaka, O., Ma, S., Sugii, K., Kurita, N., Tanaka, H., Nasu, J., Motome, Y. and Others, Majorana quantization and halfinteger thermal quantum Hall effect in a Kitaev spin liquid. Nature . 559 , 227-231 (2018).
- [7] Czajka, P., Gao, T., Hirschberger, M., Lampen-Kelley, P., Banerjee, A., Quirk, N., Mandrus, D., Nagler, S. and Ong, N. Planar thermal Hall effect of topological bosons in the Kitaev magnet α -RuCl . 3 Nature Materials . 22 , 3641 (2023).
- [8] Liu, H. and Khaliullin, G. Pseudospin exchange interactions in d 7 cobalt compounds: Possible realization of the Kitaev model. Phys. Rev. B . 97 , 014407 (2018,1).

J.H. and N.P.O. The crystals were grown and characterized by R.Z. and R.J.C. The manuscript was written by J.H. and N.P.O. with input from all authors.

## Competing financial interests

The authors declare no competing financial interests.

## Additional Information

Supplementary Materials is available in the online version of the paper.

Correspondence and requests for materials should be addressed to J.H. (jiayih@princeton.edu) or N.P.O. (npo@princeton.edu).

- [9] Regnault, L., Burlet, P. and Rossat-Mignod, J. Magnetic ordering in a planar X - Y model: BaCo (AsO ) . 2 4 2 Physica B+C . 86-88 pp. 660-662 (1977).
- [10] Regnault, L., Rossat-Mignod, J., Villain, J. and Combarieu, A. Specific heat of the quasi 2D-XY helimagnets BaCo (AsO ) 2 4 2 and BaCo (PO ) . 2 4 2 Journal De Physique Colloques . 39 , C6-759-C6-761 (1978)
- [11] Regnault, L., Boullier, C. and Lorenzo, J. Polarizedneutron investigation of magnetic ordering and spin dynamics in BaCo (AsO ) 2 4 2 frustrated honeycomb-lattice magnet. Heliyon . 4 , e00507 (2018).
- [12] Zhong, R., Gao, T., Ong, N. and Robert J. Cava Weakfield induced nonmagnetic state in a Co-based honeycomb. Science Advances . 6 , eaay6953 (2020).
- [13] Halloran, T., Desrochers, F., Zhang, E., Chen, T., Chern, L., Xu, Z., Winn, B., Graves-Brook, M., Stone, M., Kolesnikov, A. and Others, Geometrical frustration versus Kitaev interactions in BaCo 2 (AsO ) . 4 2 Proceedings Of The National Academy Of Sciences . 120 , e2215509119 (2023).
- [14] Maksimov, P., Ushakov, A., Pchelkina, Z., Li, Y., Winter, S. and Streltsov, S. Ab initio guided minimal model for the 'Kitaev' material BaCo (AsO ) : Importance of 2 4 2 direct hopping, third-neighbor exchange, and quantum fluctuations. Physical Review B . 106 , 165131 (2022)
- [15] Huyan, S., Schmidt, J., Gati, E., Zhong, R., Cava, R., Canfield, P. and Bud'Ko, S. Hydrostatic pressure effect on the Co-based honeycomb magnet BaCo (AsO ) . 2 4 2 Physical Review B . 105 , 184431 (2022)
- [16] Das, S., Voleti, S., Saha-Dasgupta, T. and Paramekanti, A. XY magnetism, Kitaev exchange, and long-range frustration in the J eff = 1 2 honeycomb cobaltates. / Physical Review B . 104 , 134425 (2021)
- [17] Winter, S. Magnetic couplings in edge-sharing high-spin d 7 compounds. Journal Of Physics: Materials . 5 , 045003

(2022)

- [18] Mukharjee, P., Shen, B., Erdmann, S., Jesche, A., Kaiser, J., Baral, P., Zaharko, O., Gegenwart, P. and Tsirlin, A. Intermediate field-induced phase of the honeycomb magnet BaCo (AsO ) . 2 4 2 ArXiv Preprint arXiv:2403.04466 . (2024)
- [19] Tu, C., Dai, D., Zhang, X., Zhao, C., Jin, X., Gao, B., Dai, P. and Li, S. Evidence for gapless quantum spin liquid in a honeycomb lattice. ArXiv Preprint arXiv:2212.07322 . (2022).
- [20] Zhang, X., Xu, Y., Halloran, T., Zhong, R., Broholm, C., Cava, R., Drichko, N. and Armitage, N. A magnetic continuum in the cobalt-based honeycomb magnet BaCo (AsO ) . 2 4 2 Nature Materials . 22 , 58-63 (2023).
- [21] Bose, A., Routh, M., Voleti, S., Saha, S., Kumar, M., Saha-Dasgupta, T. and Paramekanti, A. Proximate Dirac spin liquid in the honeycomb lattice J 1 -J 3 XXZ model: Numerical study and application to cobaltates. Physical Review B . 108 , 174422 (2023)
- [22] Wolfowicz, G., Heremans, F., Anderson, C., Kanai, S., Seo, H., Gali, A., Galli, G. and Awschalom, D. Quantum guidelines for solid-state spin defects. Nature Reviews Materials . 6 , 906-925 (2021).
- [23] See Supplemental Information.
- [24] Ponomaryov, A., Schulze, E., Wosnitza, J., LampenKelley, P., Banerjee, A., Yan, J., Bridges, C., Mandrus, D., Nagler, S., Kolezhuk, A. and Zvyagin, S. Unconventional spin dynamics in the honeycomb-lattice material α -RuCl : High-field electron spin resonance studies. 3 Physical Review B . 96 , 241107 (2017).
- [25] Ponomaryov, A., Zviagina, L., Wosnitza, J., LampenKelley, P., Banerjee, A., Yan, J., Bridges, C., Mandrus, D., Nagler, S. and Zvyagin, S. Nature of magnetic excitations in the high-field phase of α -RuCl . 3 Physical Review Letters . 125 , 037202 (2020).
- [26] Wellm, C., Zeisner, J., Alfonsov, A., Wolter, A., Roslova, M., Isaeva, A., Doert, T., Vojta, M., B¨ uchner, B. and Kataev, V. Signatures of low-energy fractionalized excitations in α -RuCl 3 from field-dependent microwave absorption. Physical Review B . 98 , 184408 (2018).
- [27] Balz, C., Lampen-Kelley, P., Banerjee, A., Yan, J., Lu,
- Z., Hu, X., Yadav, S., Takano, Y., Liu, Y., Tennant, D., Lumsden, M., Mandrus, D. and Nagler, S. Finite field regime for a quantum spin liquid in α -RuCl . 3 Physical Review B . 100 , 060405 (2019).
- [28] Zhao, X., Wu, J., Zhao, Z., He, Z., Song, J., Zhao, J., Liu, X., Sun, X. and Li, X. Heat switch effect in an antiferromagnetic insulator Co V O . 3 2 8 Applied Physics Letters . 108 , 242405 (2016).
- [29] Hong, X., Gillig, M., Yao, W., Janssen, L., Kocsis, V., Gass, S., Li, Y., Wolter, A., B¨chner, B. and Hess, C. u Phonon thermal transport shaped by strong spin-phonon scattering in a Kitaev material Na 2 Co TeO . 2 6 Npj Quantum Materials . 9 , 18 (2024).
- [30] Hentrich, R., Wolter, A., Zotos, X., Brenig, W., Nowak, D., Isaeva, A., Doert, T., Banerjee, A., LampenKelley, P., Mandrus, D., Nagler, S., Sears, J., Kim, Y., B¨ uchner, B. and Hess, C. Unusual Phonon Heat Transport in α -RuCl : 3 Strong Spin-Phonon Scattering and Field-Induced Spin Gap. Phys. Rev. Lett. . 120 , 117204 (2018,3).
- [31] Li, Q., Zhao, Z., Fan, C., Zhang, F., Zhou, H., Zhao, X. and Sun, X. Phonon-glass-like behavior of magnetic origin in single-crystal Tb 2 Ti O . 2 7 Phys. Rev. B . 87 , 214408 (2013,6).
- [32] Wang, X., Fan, C., Zhao, Z., Tao, W., Liu, X., Ke, W., Zhao, X. and Sun, X. Large magnetothermal conductivity of HoMnO 3 single crystals and its relation to the magnetic-field-induced transitions of magnetic structure. Phys. Rev. B . 82 , 094405 (2010,9).
- [33] Berggold, K., Baier, J., Meier, D., Mydosh, J., Lorenz, T., Hemberger, J., Balbashov, A., Aliouane, N. and Argyriou, D. Anomalous thermal expansion and strong damping of the thermal conductivity of NdMnO 3 and TbMnO 3 due to 4f crystal-field excitations. Physical Review B-Condensed Matter And Materials Physics . 76 , 094418 (2007).
- [34] Pocs, C., Leahy, I., Zheng, H., Cao, G., Choi, E., Do, S., Choi, K., Normand, B. and Lee, M. Giant thermal magnetoconductivity in CrCl 3 and a general model for spinphonon scattering. Phys. Rev. Res. . 2 , 013059 (2020,1).
Jiayi Hu Ruidan Zhong Peter Czajka Tong Gao R. J Cava N. P. Ong
Scaling behavior and giant field-enhancement of the thermal conductivity in the honeycomb antiferromagnet BaCo2(AsO4)2
The layered honeycomb material BaCo$_2$(AsO$_4$)$_2$ (BCAO) is of topical
interest because its magnetic state is related to that of the Kitaev magnet
$\alpha$-RuCl$_3$. Using thermal transport to probe how magnetic excitations
interact with phonons in the magnetically disordered regime, we have uncovered
an unusually large enhancement of the thermal conductivity $\kappa_{xx}$ in an
in-plane magnetic field ${\bf H}$. Just above the N\'{e}el temperature $T_{\rm
N}$, a field of 13 T increases $\kappa_{xx}$ by a factor $\sim 211$, much
larger than reported previously in any magnetic insulator. Interestingly,
$\kappa_{xx}(H,T)$ exhibits a scaling behavior in the entire magnetically
disordered region that surrounds the ordered zigzag state. The ratio $\Delta
\kappa_{xx}(H,T)/\kappa_{xx}(13,T)$, measured throughout the disordered region,
collapses to a one-parameter scaling function ${\rm exp}(-1/gx)$ (where $x =
\mu_{\rm B}B/k_{\rm B}T$ and $g$ is a constant).