Abstract
Breaking of spatial inversion symmetry induces unique phenomena in condensed matter. In particular, by combining this symmetry with magnetic fields or another type of time-reversal symmetry breaking, noncentrosymmetric materials can be made to exhibit nonreciprocal responses, which are responses that differ for rightward and leftward stimuli. However, the effect of spatial inversion symmetry breaking on thermal transport in uniform media remains to be elucidated. Here, we show nonreciprocal thermal transport in the multiferroic helimagnet TbMnO3. The longitudinal thermal conductivity depends on whether the thermal current is parallel or antiparallel to the vector product of the electric polarization and magnetization. This phenomenon is thermal rectification that is controllable with external fields in a uniform crystal. This discovery may pave the way to thermal diodes with controllability and scalability.
INTRODUCTION
Unconventional physical properties in noncentrosymmetric materials have been investigated for long (1, 2). Besides the classic examples such as natural optical activity and piezoelectricity, the spin-orbit interaction further enriches the effect of spatial inversion symmetry (SIS) breaking, as exemplified by the Rashba effect (3). When SIS and time-reversal symmetry (TRS) are broken simultaneously, there emerge distinct phenomena denoted as nonreciprocal responses, which are responses that differ for rightward and leftward stimuli (4). Nonreciprocal responses were first studied in terms of optical properties (5, 6). When TRS and SIS are simultaneously broken, the optical properties such as absorption and luminescence become dependent on the sign of the optical wave vector irrespective of the polarization. Similar directional dependences were observed for magnon and phonon propagation (7–10). Regarding electronic transport, resistivity becomes directionally dependent; that is to say, rectification appears in materials without TRS and SIS (11–14). The symmetrical rule for the unidirectional response is common among these nonreciprocities. In chiral materials, in which any mirror symmetry is broken, directional dependence is observed along the magnetization M for any nonreciprocities. On the other hand, for polar materials, which have finite electric polarization P, the directional dependence is observed along P × M. In this paper, we study nonreciprocal thermal transport in a polar system.
We chose the multiferroic spiral magnet TbMnO3 as a sample (15). TbMnO3 shows transverse spiral magnetic order with the propagation vector along the b axis, accompanied by P along the c axis. We measured nonreciprocal thermal transport in magnetic fields parallel to the a axis, in which a conical magnetic structure is realized (Fig. 1A). Let us discuss the phononic state in this magnetic structure to study the thermal transport microscopically. As Hamada and co-workers theoretically suggested, the breaking of SIS induces phonon angular momentum L (16), similarly to the Rashba and Dresselhaus effects in electronic states (17). In polar systems, the induced L at momentum k is along k × P. Because L is caused by circular or ellipsoidal phonon polarization in the plane perpendicular to L, it originates from the mixing of two linearly polarized phonon modes in the plane. In the low-energy acoustic branches, these two modes respectively correspond to the longitudinal mode and one of the transverse modes with phonon polarization along P. The solid lines in Fig. 1B show the expected acoustic phonon modes for ferroelectric TbMnO3 at zero field. The longitudinal mode and one of the transverse modes have the phonon angular momenta L0 and −L0, respectively, along k × P as a result of mixing. On the other hand, the other transverse mode has linear polarization along k × P, and therefore, the phonon angular momentum is zero. When k changes its sign, L is reversed, but the energy is unchanged. By the application of the magnetic field along the a axis, the energies of phonon modes shift depending on L (dotted lines in Fig. 1B), owing to the so-called Raman-type spin-phonon interaction HRaman = K m∙L (18, 19), where m is a magnetic moment. Thus, the phonon energy dispersion becomes asymmetric in the magnetic field. This asymmetry of phononic dispersion is the origin of thermal nonreciprocity.
(A) Illustration of magnetic structure in a magnetic field H parallel to the a axis. (B) Expected acoustic phonon modes along the b axis at H = 0 (solid lines) and finite H (dotted lines). (C) Illustration of the experimental setup for nonreciprocal thermal transport measurement. The distance between two thermometers is 1.05 mm. The sample thickness and width are 0.4 and 1.15 mm, respectively.
RESULTS
Experimental setup
Thermal conductivity is usually measured in an experimental setup where one side of the sample is attached to a thermal bath and the other side to a heater (Fig. 1C). Because one needs to modify the experimental setup to reverse the heat current, it is difficult to estimate small rectification precisely. Here, we demonstrate thermal nonreciprocity in the multiferroic spiral magnet TbMnO3 by a comparison of the thermal conductivities for ±P and ±M states, without changing the direction of the thermal current. The thermal conductivity is expected to be rectified along the b axis (∣∣P × M). Reversal of thermal current is equivalent to P (M) reversals and 180° rotation of experimental setup around M (P). The 180° rotation does not change any physical responses, and therefore, the reversal of P or M corresponds to reversal of the thermal current. Thus, we can estimate a small thermal rectification by comparative measurements with positive and negative P and M. For this purpose, we attached a pair of electrodes to control the electric polarization in addition to the conventional thermal transport setup, as shown in Fig. 1C.
Magnetic, electric, and thermal transport properties of TbMnO3
TbMnO3 shows several magnetic phase transitions in the low-temperature region. In Fig. 2A, we show the temperature dependence of magnetization at 0.5 T for TbMnO3. The anomalies at 7, 27, and 43 K indicate the magnetic transitions (15, 20). The manganese magnetic moments order incommensurately below 43 K. While the magnetic moments are almost collinear above 27 K, a cycloidal magnetic structure emerges, and ferroelectric polarization is induced along the c axis below 27 K (Fig. 2B). The increase at 7 K corresponds to the incommensurate order of terbium moments. Note that terbium oxides frequently show gigantic magnetoelastic phenomena such as phonon Hall effect (21, 22) and acoustic Faraday effect (23). The thermal nonreciprocity in the present study also seems to be related to the strong magnetoelastic coupling of terbium moments, as discussed later. The temperature dependence of thermal conductivity is shown in Fig. 2C. The thermal conductivity is relatively small and monotonically decreases with lowering temperature. A previous study suggests that these are caused by strong scattering due to terbium moments (24).
(A) Temperature dependence of magnetization at H = 0.5 T parallel to the c axis. Inset shows the low-temperature and large magnetization region. (B) Temperature dependence of electric polarization parallel to the c axis at H = 0. (C) Temperature dependence of thermal conductivity at H = 0. The vertical dashed lines indicate the magnetic transition temperatures. The kink temperature in the magnetization is slightly higher than that in electric polarization around 27 K, owing to the application of the magnetic field.
Nonreciprocal thermal transport
Next, let us discuss the magnetic field dependence of the magnetic and thermal properties at 4.2 K to demonstrate the thermal nonreciprocity. Figure 3A shows the magnetization as a function of the magnetic field parallel to the a axis at 4.2 K. The magnetization shows a step-like increase around 2 T. Because the c-axis polarization due to the spiral spin structure of manganese moments is robust in this magnetic field range (25), the transition is relevant to the ferromagnetic ordering of terbium moments. In Fig. 3B, we show the magnetic field dependence of thermal conductivity at 4.2 K. The thermal conductivity shows spike-like anomalies at the magnetic transition fields. A small peak is also observed at 0 T. Small hysteresis is observed in the low-field region below the transition field. In Fig. 3C, we plot the magnetic field dependence of thermal conductivity above the metamagnetic transition field on a magnified scale for the ±P states. It shows an asymmetric magnetic field dependence, and the asymmetry is reversed when the electric polarization is reversed. To scrutinize this phenomenon, we show in Fig. 3D the asymmetric part of the magnetic field–dependent thermal conductivity
(A) Magnetization curve along the a axis at 4.2 K. (B and C) Magnetic field variation of thermal conductivity at 4.2 K for H || a. In (C), the thermal conductivity scale is magnified around κ(H)/κ(0) = 0.83. (D) Magnetic field dependence of thermal nonreciprocity
In Fig. 3E, we show the thermal current dependence of
The nonreciprocity of thermal conductivity
(A to H) Magnetic field variation of thermal nonreciprocity Δκ/κ for ±P and H || a at various temperatures. (I) Temperature dependence of averaged thermal nonreciprocity Δκav/κ for ±P. The error bars are estimated as the error of thermal conductivity from the temperature fluctuation during the measurement divided by the square of number of averaged data points. Solid lines are merely guides for the eyes.
DISCUSSION
Last, let us discuss the origin of the nonreciprocal thermal transport. As mentioned above, the nonreciprocal signal correlates with the terbium magnetic state. Nevertheless, the terbium magnetic excitation hardly contributes to the thermal transport because it has a dispersionless local nature (26) and cannot carry heat. Therefore, inelastic scattering of phonons induced by the terbium moments seems to be the most plausible origin. When a phonon with angular momentum L is scattered by a terbium moment with total angular momentum J, the sum of L and J must be conserved through the entire scattering process. Therefore, the transition matrix is expressed as
To understand how the terbium magnetic scattering affects the thermal transport properties, let us calculate the nonreciprocal thermal transport for a simplified model. This is just for a qualitative understanding of nonreciprocal thermal transport in this system, while a more elaborate theory is needed for the quantitative understanding. We assume one-dimensional asymmetric phonon bands as indicated by the dotted lines in Fig. 1B and an induced phonon momentum magnitude ∣L0∣ = 1. Therefore, the phonon angular momenta for three acoustic branches are L = 0, ± 1. The phonon density for the angular momentum L at momentum k under a temperature gradient is expressed as g(k, L) = g0(k, L) + g′(k, L), where g0 and g′ are the phonon density in thermal equilibrium and the deviation, respectively. Assuming the relaxation time approximation, we get
Here, ϵ is phonon energy vk. Similarly, for
Then, the nonreciprocity is
This relation is finite when v+, + ≠ v+, − or v−, + ≠ v−, −. Thus, the scattering by ferromagnetic terbium moments and simultaneous breaking of SIS and TRS induce the nonreciprocal thermal transport. This is independent of the magnitude of the thermal current density. If the magnitude of the induced phonon angular momentum L0 is less than 1, the nonreciprocity is reduced but still finite unless L0 becomes zero. This simple model is consistent with the experimental observations; the nonreciprocity is fairly enhanced by the ferromagnetic alignment of terbium moments and Δκ/κ is independent of thermal current density.
In summary, we have demonstrated nonreciprocal thermal transport in multiferroic TbMnO3; the thermal transport is rectified along the vector product of the polarization and magnetization. Similar thermal rectification phenomena were already studied theoretically and experimentally. Zhang et al. (27) theoretically suggest that thermal rectification can be induced by a magnetic field in a nanoscale three terminal junctions. Some theories also suggest that nonlinear lattice dynamics induces the thermal rectification (28, 29). Experimentally, thermal rectifications were previously realized by using combinations of two different materials or asymmetrically shaped media (30–35). On the other hand, we achieved the thermal rectification in a uniform bulk crystal that is induced by the phononic band asymmetry that originates from the symmetry breaking. The present result suggests that the thermal rectification should be emergent in any materials without both the time-reversal and spatial inversion symmetries. In some of these materials, the symmetries are broken even at room temperature. Therefore, thermal rectification at room temperature is also expected. Much larger nonreciprocity may also show up in some of these materials. A promising way to increase the magnitude of rectification is to use the material with the large electric polarization. The electric polarization in TbMnO3 is orders of magnitude smaller than that in typical ferroelectrics because the ferroelectricity in TbMnO3 is induced by the magnetic structure. Materials with larger polarization are expected to show larger thermal nonreciprocity. The thermal rectification seems useful for the efficient management of heat in various circumstances (30). In the cases of previously realized thermal rectification (30–35), the forward and reverse directions are determined by structural asymmetries, and the directionality cannot be reversed by external fields. In addition, an asymmetric structure limits the scale of the device. In the present case, the thermal rectification is realized in a uniform crystal and can be reversed by the external electric or magnetic fields. If the magnitude of thermal nonreciprocity becomes large enough to make practical applications feasible, the scalability of the uniform crystal and the ability to control the phenomenon with the external fields are substantial advantages over the previous approaches. While thermal current can be controlled by the application of electric current with the use of thermoelectric effect as exemplified by a recent paper (36), constant application of electric fields is not needed after the alignment of P, and the energy consumption is much less in the present case. Thus, the present results suggest that the exploration of thermal properties in symmetry-broken materials is a new research strategy that should be investigated toward the realization of thermal diode devices.
MATERIALS AND METHODS
A TbMnO3 single crystal was grown by means of the floating zone technique. The conditions were almost the same as those in previous reports (15, 25). Magnetization was measured with the use of a superconducting quantum interference device. Thermal conductivity measurements were performed by means of the steady-state method in a superconducting magnet. The measurement atmosphere was evacuated down to 3 × 10−3 Pa. Heat current was generated by chip resistance. The thermal gradient was measured by two thermometers (CERNOX, Lake Shore Cryotronics Inc.). The sample temperature is estimated as the average of the measured temperatures. To control the electric polarization, a pair of Au/Ti electrodes were formed on the surface of the thermal transport sample. The thicknesses of the Au and Ti were 20 and 2 nm, respectively. Before the thermal conductivity measurement, we performed the following poling procedure to obtain the full alignment of ferroelectric polarization (15, 25). We applied an electric field between the electrodes as large as 750 V/mm in high-temperature paraelectric state and slowly decreased the temperature to the measurement temperature. Before the thermal conductivity measurement, we turned off the external electric field and wait a long time so that the relative temperature fluctuation δT/T is comparable to 10−4.
SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/6/40/eabd3703/DC1
This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.
REFERENCES AND NOTES
- Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).