Processing math: 52%

Nonreciprocal thermal metamaterials: Methods and applications

Zhengjiao Xu, Chuanbao Liu, Xueqian Wang, Yongliang Li, Yang Bai

图(6)

分享

计量
  • 文章访问数:  1270
  • HTML全文浏览量:  212
  • PDF下载量:  52
  • 被引次数: 0

目录

Cite this article as:

Zhengjiao Xu, Chuanbao Liu, Xueqian Wang, Yongliang Li, and Yang Bai, Nonreciprocal thermal metamaterials: Methods and applications, Int. J. Miner. Metall. Mater., 31(2024), No. 7, pp.1678-1693. https://dx.doi.org/10.1007/s12613-023-2811-6
Zhengjiao Xu, Chuanbao Liu, Xueqian Wang, Yongliang Li, and Yang Bai, Nonreciprocal thermal metamaterials: Methods and applications, Int. J. Miner. Metall. Mater., 31(2024), No. 7, pp.1678-1693. https://dx.doi.org/10.1007/s12613-023-2811-6
引用本文 PDF XML SpringerLink
特约综述

非互易热超材料:方法和应用

文章亮点

(1) 归纳总结了实现非互易热超材料的不同方法 (2) 详细阐述了非互易热超材料(热传导和热辐射)的发展和应用前景 (3) 介绍了非互易超材料在其他拉普拉斯物理场中的应用
在天然材料中,由于受限于热流固有的各向同性扩散属性,难以实现非互易传热行为。而利用人工构建的热超材料,则可以打破系统互易性,实现对热流大小和方向的任意调控,在隔离保护、单向传输和能量收集等方面具有重要应用前景。本文首先从昂萨格互易定理出发,介绍了非互易热超材料设计的基本理论。其次,探讨了实现非互易传热的不同方法,包括使用非线性材料或结构、时空调制和引入角动量偏差等。进一步,讨论了实现非互易热辐射方面的实现方法。最后,介绍了非互易性在其他拉普拉斯物理场的潜在应用,并展望了非互易热超材料的未来研究方向。

 

Invited Review

Nonreciprocal thermal metamaterials: Methods and applications

Author Affilications
    Corresponding author:

    Chuanbao Liu      E-mail: cbliu@ustb.edu.cn

    Yang Bai      E-mail: baiy@mater.ustb.edu.cn

  • Received: 13 August 2023; Revised: 25 November 2023; Accepted: 13 December 2023; Available online: 14 December 2023
Nonreciprocity of thermal metamaterials has significant application prospects in isolation protection, unidirectional transmission, and energy harvesting. However, due to the inherent isotropic diffusion law of heat flow, it is extremely difficult to achieve nonreciprocity of heat transfer. This review presents the recent developments in thermal nonreciprocity and explores the fundamental theories, which underpin the design of nonreciprocal thermal metamaterials, i.e., the Onsager reciprocity theorem. Next, three methods for achieving nonreciprocal metamaterials in the thermal field are elucidated, namely, nonlinearity, spatiotemporal modulation, and angular momentum bias, and the applications of nonreciprocal thermal metamaterials are outlined. We also discuss nonreciprocal thermal radiation. Moreover, the potential applications of nonreciprocity to other Laplacian physical fields are discussed. Finally, the prospects for advancing nonreciprocal thermal metamaterials are highlighted, including developments in device design and manufacturing techniques and machine learning-assisted material design.

 

  • With the miniaturization of devices and chip integration in communication technology and electronic packaging, there is an increasing demand for materials with superior thermal properties. Based on the second law of thermodynamics, heat spontaneously flows from high temperature to low temperature in homogenous materials when a temperature gradient is present. This process of energy transfer is ubiquitous in nature. However, to guarantee the stability of the target system, manipulating the heat flux at will becomes a significant task. Because of the isotropic material parameters and fixed geometry shape, traditional thermal materials have a limited ability to tune thermal properties and adapt to complex heat flux distributions. Thus, it is urgent to examine novel thermal management methods. In this regard, thermal metamaterials inspired by electromagnetic metamaterials have attracted considerable attention to overcome the limitations of natural materials [12].

    Metamaterials are artificially engineered media comprising a periodic or nonperiodic array of unit structures that demonstrate remarkable physical properties and applications such as electromagnetic invisibility cloaks [37], three-dimensional holograms [8], and achromatic metalens [910]. Previous studies on metamaterials were mainly focused on manipulating electromagnetic waves and investigated through the use of identical structural units arranged in a periodic layout to realize exotic physical phenomena such as negative refractive [11] and reversed Doppler effect [12]. Later, Pendry et al [4]. and Alù et al [13]. proposed transformation optics and scattering elimination theory, respectively, allowing metamaterials to control electromagnetic waves more effectively through the spatial sequence design of differentiated structural units. These metamaterial design approaches are not only applicable to wave fields (e.g., electromagnetic fields) but also to Laplace physical fields (e.g., thermal fields), which satisfy the second-order partial differential equation at the steady state 2φ=0, where 2 is the Laplacian operator and φ is the potential function, which can be a physical field, a combination of physical fields, or a nonlinear term of the original field T [1415]. In 2008, Fan et al. [1] pioneered the concept of an electromagnetic invisibility cloak into the field of heat and theoretically predicted the thermal invisibility cloak. In 2012, Narayana et al. [16] proved by experiments the thermal invisibility cloak for the first time. Later, research on thermal invisibility cloaks has increased, and transient thermal invisibility cloaks and multifunction thermal invisibility cloaks have been explored successively [1622]. With the proposal and implementation of the thermal invisibility cloak, several thermal metamaterials have emerged, including heat flow concentrators for collecting heat energy [17,2325], heat flow rotators for rotating heat fields [16,2627], thermal illusion devices for camouflage [2831], and thermal encoding [3233].

    The abovementioned thermal metamaterials are reciprocal and present symmetry in energy transfer in opposite directions. However, several scenarios require more precise and adaptive heat field regulation to overcome the inherent isotropic diffusion of heat flow and achieve directional control over both the magnitude and orientation of heat flux. Nonreciprocal thermal metamaterials provide a promising solution for achieving greater degrees of freedom in thermal field regulation, resulting in numerous novel applications. For example, researchers have proposed thermal diodes to allow unidirectional heat transfer [34], developed a geometric heat pump that pumps extra heat ably diffusing from cold to hot [35], and demonstrated a nonreciprocal infrared thermal emitter caused by a spatiotemporal modulation grating [3638].

    The design of nonreciprocal thermal metamaterials is mainly based on the Onsager reciprocity theorem. In this review, we offer a detailed introduction to this theory and examine many typical implementation methods, including nonlinearity, spatiotemporal modulation, and angular momentum bias, for realizing nonreciprocity in the thermal fields. Nonreciprocal thermal radiation using magnetic response, time-variant systems, and optical nonlinearity are also discussed. Finally, it outlooks several future directions for the development of nonreciprocal thermal metamaterials.

    Reciprocity requires the system to exhibit symmetry in response to energy transfers that occur in opposite directions. Specifically, the transmission channel should respond symmetrically to the input source when the transmitter and receiver positions are interchanged. The notion of reciprocity was first discussed theoretically by Stokes and Helmholtz for light waves, followed by successive proposals of Lorentz reciprocity and Onsager reciprocity. Among these, the Onsager reciprocity can be applied to different irreversible physical processes such as acoustic, electromagnetic, mechanical, thermoelectric, and diffusion phenomena. In this section, we outline how the Onsager principle of microscopic reversibility is derived in a heat field.

    At the microscopic level, most physical processes show time-reversal symmetry. On the basis of microscopic reversibility, Onsager explained the basic relationship between the time-reversal invariance of the microscopic dynamic equation and reciprocity. The Onsager reciprocity principle is derived from four fundamental assumptions [39]: time-reversal symmetry of microscopic equations, linearity, causality (an irreversible process that increases entropy), and thermodynamic quasiequilibrium (a system reaches a stable state after interacting with its surroundings for a sufficient amount of time).

    Onsager reciprocity emphasizes that irreversible processes originate from a generalized flow J (the amount of heat flux and concentration passing through a unit area per unit time) driven by a generalized force X (temperature gradient, chemical potential gradient). In the thermal quasiequilibrium state, these two quantities follow a linear phenomenological relationship, denoted as Eq. (1).

    Ji=jLijXj (1)

    where Lij refers to the phenomenological coefficient, indicating that the i-th generalized flow Ji is affected by the j-th generalized force Xj. Based on the limitation of the microscopic reversibility hypothesis, the linear phenomenological coefficients have symmetry of Lij=Lji, which is the Onsager reciprocity relation.

    Specifically, for a three-dimensional anisotropic heat conduction process, the law of heat conduction (i.e., Fourier’s law) can also be expressed in accordance with the phenomenological relationship of Eq. (2).

    qi=LijXj,(i,j=x,y,z) (2)

    where qi refers to the heat flux, and Xj refers to the temperature gradient. In Eq. (2), the phenomenological coefficient Lij is the thermal conductivity as a tensor, with its value range in the Cartesian coordinate system (x, y, z) forming an ellipsoid characterized by three principal axes (a, b, c) of heat conduction. Based on the limitation of the linear assumption, the heat flux through any point on the surface of an object varies linearly with the temperature gradient. The phenomenological coefficientsLij (i,j=x,y,z) in the Cartesian coordinate system can be expressed as phenomenological coefficients Mnn (n=a,b,c) along the principal axes, denoted as Eq. (3).

    Lij=n=a,b,ceinejnMnn = eiaejaMaa+eibejbMbb+eicejcMcc (3)

    where eia, eja, eib, ejb, eic, and ejc(i, j=x,y,z) is the cosine of the angle between the Cartesian coordinate system and the direction corresponding to the principal axes. According to matrix arithmetic, the phenomenological coefficients Lij in Eq. (3) is symmetric and can be a diagonalized matrix. Thus, the Onsager reciprocity is satisfied in the heat conduction process.

    To realize nonreciprocity in a system, it is required to invalidate at least one of the four basic assumptions of Onsager reciprocity [39]. With regard to macroscopic heat conduction, causality is described by the second law of thermodynamics, which states that heat flows from high temperature to low temperature. This is a well-established law of nature and is hard to violate in everyday applications and most engineering systems. Albeit negative heat transfer phenomena have been achieved in certain extreme and specialized situations [1,16,27,4041], the causality or the second law of thermodynamics is not violated essentially. For instance, by using a thermal rotator in a system with a fixed temperature difference, the local heat flux can be completely reversed. However, the heat flux still takes place from the high-temperature boundary to the low-temperature boundary from the angle of the whole system [27,40]. Thus, breaking causality in the overall heat transfer process is demanding.

    Nonreciprocal heat transfer can be realized by using dissipative structures far from the distance equilibrium, which violates the thermodynamic quasiequilibrium assumption. Dissipative structures are macroscopic structures that appear from the exchange of energy and matter between the system and its surroundings when the system is in a nonequilibrium state. The formation and maintenance of such structures are accompanied by dissipation. It was demonstrated that optical nonreciprocity and unidirectional energy transmission can be attained through system energy dissipation [4243]. In heat conduction, the introduction of surface convection or radiant heat transfer allows the exchange of matter and energy with the surroundings [44]. This introduces nonlinear coupling and energy conversion into the heat transfer process, leading to a nonlinear relationship between the temperature gradient and heat flux, thus allowing nonreciprocal heat transfer. However, the principle of minimum entropy production ensures that the nonequilibrium state always converges to a fixed state within the system region over time, irrespective of external interference or internal parameter fluctuations. The disturbed dynamics will show higher entropy production than the original stationary state and eventually go back to its initial equilibrium [45]. Thus, thermodynamic quasiequilibrium is inherently stable, i.e., unbalanced processes naturally tend toward equilibrium.

    In comparison to causality and thermodynamic quasiequilibrium, breaking the linear response and time-reversal symmetry makes it easier to realize nonreciprocal heat transfer, which has triggered immense attention in this area. In the following, a detailed discussion is presented regarding achieving thermal nonreciprocity according to these two assumptions. One approach is to use nonlinear thermal materials or structures in which physical transport properties are not constrained by the reciprocity theorem, which can break the linear response (Section 3.1). The second is to employ the speed generated by the linear spatiotemporal modulation (Section 3.2) or the angular momentum bias resulting from the rotation (Section 3.3), which can break the time-reversal symmetry.

    The best way to realize thermal nonreciprocity using nonlinear thermal materials or structures lies in the realization of the nonlinearity of thermal conductivity. For passive conditions and without convection, the one-dimensional heat conduction equation is:

    ρcTt=κ2Tx2 (4)

    where T refers to the temperature, x refers to the coordinates of the heat transfer direction, t refers to the time of heat transfer, κ refers to the thermal conductivity, and ρ and c refer to the density and specific heat capacity of the material, respectively. Generally, the thermal conductivity remains constant irrespective of temperature, thus making the governing equation a linear partial differential equation. As a consequence, Green’s function can be utilized to solve the temperature distribution in heat conduction processes [4648]. Based on the linear superposition principle, the field distribution of the heat source under the same boundary conditions can be acquired using Green’s function (field generated by the point heat source). In this case, Green’s function is symmetric, i.e., the field is symmetric with respect to the source [49]. Thus, the heat conduction process satisfies the reciprocity theorem on both local and global scales [50]. However, if the thermal conductivity of the material is temperature-dependent, the governing equation transforms into a nonlinear partial differential equation:

    ρcTt=κ(T)2Tx2 (5)

    Consequently, obtaining a solution via Green’s function becomes unattainable and makes the reciprocity relation invalid.

    To realize a nonreciprocal heat conduction process, nonlinear thermal materials with temperature-dependent thermal conductivity can be employed. Considering nonlinear isotropic thermal materials 1 and 2, material 1 demonstrates a positive correlation between thermal conductivity and temperature, while material 2 has a negative correlation (the top left of Fig. 1(a)). First, the nonreciprocity of heat conduction can be attained using a single nonlinear thermal material 1 combined with the asymmetry of space geometry. Because more parts of the system can be heated from the right side, when the heat source and cold source are exchanged, the amount of heat in the forward direction should be small compared with that in the backward direction. Therefore, the heat conduction process in the two directions is nonreciprocal and is confirmed (the bottom left of Fig. 1(a)). From this phenomenon, this structure can be utilized as a thermal diode [51] or a thermal rectifier [52]. However, the thermal rectification ratio (η=|Q(ΔT)|/|Q(ΔT)||Q(ΔT)||Q(ΔT)| (where |Q(ΔT)| and |Q(ΔT)| denote the heat fluxes in the forward and backward directions, respectively) of the structure is as low as 100.48% when ΔT = 25 K, leading to suboptimal performance as a thermal diode or rectifier. To further enhance the thermal rectification ratio, asymmetric nonlinearity is proposed as a possible method by simultaneously utilizing nonlinear thermal materials 1 and 2 (the bottom right of Fig. 1(a)). The heat flux through the asymmetric geometry and asymmetric nonlinearity system is displayed in the upper right corner of Fig. 1(a). Due to the contrasting temperature dependence of materials 1 and 2, both regions show high (low) thermal conductivity concurrently during forward (backward) heat conduction, giving rise to an outstandingly high thermal rectification ratio of up to 184.32% at ΔT = 25 K. These results imply that realizing thermal rectification in one-dimensional thermal diffusion requires the inseparability of spatial and temperature-dependent thermal conductivities [53].

    Figure  1.  Nonreciprocity resulting from nonlinear thermal materials: (a) thermal diode made of nonlinear material. Two nonlinear materials (top left), steady-state temperature distribution in geometric asymmetry composed of material 1 (bottom left) and geometric symmetry structure spliced by materials 1 and 2 (bottom right), and heat flow across the two systems (top right). Reprinted figure with permission from [50] as follows: Y. Li, J.X. Li, M.H. Qi, C.W. Qiu, and H.S. Chen, Phys. Rev. B, 103, 014307 (2021). Copyright 2021 by the American Physical Society; (b) nonlinear structure-based macroscopic thermal diode. Diode in the off-state (first column) or on-state (second column). Reprinted figure with permission from [54] as follows: Y. Li, X. Shen, Z. Wu, et al., Phys. Rev. Lett., 115, 195503 (2015). Copyright 2015 by the American Physical Society; (c) isolation thermostat constructed with a nonlinear structure. Left: structure of the isolation thermostat and material properties of two shape memory alloys; right: effect of the isolation thermostat structure at different heat source temperatures. Reprinted figure with permission from [55] as follows: X. Shen, Y. Li, C. Jiang, and J. Huang, Phys. Rev. Lett., 117, 055501 (2016). Copyright 2016 by the American Physical Society.

    Besides using nonlinear thermal materials directly, a structure with a specific response to temperature can also be used to attain the nonlinear equivalent thermal conductivity of the system. For instance, a macroscopic thermal diode can be fabricated with shape memory alloys that deform with temperature (Fig. 1(b)) [54]. The left-hand alloy demonstrates flattening behavior at high temperatures and bending behavior at low temperatures, whereas the right-hand alloy has flattening behavior at low temperatures and bending behavior at high temperatures. Temperature-dependent variations in alloy shape indirectly affect the equivalent thermal conductivity of the system. When the left container (blue) is filled with cold water and the right container (red) is filled with hot water, both shape memory alloys simultaneously experience warping (the top left of Fig. 1(b)). During this process, the system shows an extremely low equivalent thermal conductivity, inhibiting heat transfer from the heat source to the cold one, leading to an off-state for the diode (the bottom left of Fig. 1(b)). However, switching the locations of the two containers, both shape memory alloys are straight (the top right of Fig. 1(b)). At this stage, there is a significantly high effective thermal conductivity within the system, facilitating heat flow throughout and giving rise to an on-state for the diode (the bottom right of Fig. 1(b)).

    Similarly, sandwich structures of isolation thermostats can be realized using shape memory alloys (SMA) to achieve nonreciprocal isolation protection [5556]. In the top left of Fig. 1(c), the middle section (light yellow) comprises copper, whereas the terminals are composed of distinct shape memory alloys (gray) and phosphor copper (orange) with identical phase transition temperatures but opposite phase transition characteristics: Type-A flattens at high temperature and bends at low temperature, whereas Type-B has the opposite behavior. The effective thermal conductivities of Type-A and Type-B as a function of temperature are shown at the bottom of Fig. 1(c). When connected to the heat source and the cold source, both shape memory alloys simultaneously experience warping, leading to a simultaneous decrease in effective thermal conductivity at both ends and forming an isolated protection effect in the middle region (as shown in the right panel of Fig. 1(c)). The central temperature of the experimental group remains relatively stable compared with that of the control group. In addition, the nonlinear heat transfer process can be constructed by using the variation of thermal conductivity with ambient temperature during phase transformation processes, such as solid–liquid phase transition [57] and liquid–solid composites [58].

    Because the thermal conductivity of nonlinear thermal materials depends on temperature, the constructed thermal nonreciprocal devices are limited to specific operating temperatures. To overcome this limitation, a spatiotemporal medium is introduced, which breaks the symmetry of temporal and spatial inversion in linear systems, allowing the realization of thermal nonreciprocity.

    Without convection and under passive conditions, nonreciprocal heat transport can be achieved by introducing traveling-wave spatiotemporal modulation (resembling a traveling wave-like change over time and space) in both volumetric heat capacity (the product of mass density ρ and specific heat capacity c) and thermal conductivity. In this scenario, the one-dimensional heat conduction equation can be reformulated as

    C(xv0t)Tt=x[κ(xv0t)Tx] (6)

    where C(xv0t) refers to the volumetric heat capacity and κ(xv0t) denotes the thermal conductivity. Considering a simple modulation case for these two parameters,

    C(xv0t)=C0[1+ΔCcos(β(xv0t))] (7)
    κ(xv0t)=κ0[1+Δκcos(β(xv0t)+δ)] (8)

    Both demonstrate periodic variations between the minimum value (C0(1Δρ), κ0(1Δκ)) and the maximum value (C0(1+Δρ), κ0(1+Δκ)), where β is the wavenumber, v0 is the speed of movement, C0 and κ0 are the volumetric heat capacity and thermal conductivity of the initial homogeneous material, respectively, and ΔC and Δκ are their modulation amplitudes.

    The modulation of volumetric heat capacity encompasses two distinct scenarios: one is the modulation of mass density, and the other is the modulation of heat capacity. In continuous media, mass conservation is maintained by the continuity equation ρ/t+(ρv0)=0, incorporating an additional convection term of (CTv0) that is not included in Eq. (6). Therefore, without an external energy or mass input, a physical system keeps a constant mass density over time as governed by the principle of mass conservation [59]. Thus, rather than relying on density modulation, a tunable specific heat capacity can serve as a viable alternative since the product of density and heat capacity is what ultimately matters.

    According to Bloch’s theorem, the temperature distribution of heat wave diffusion is

    T=ϕ(xv0t)ei(βxωt) (9)

    where β and ω are the wavenumber and frequency, respectively, and ϕ(xv0t) is an amplitude-modulated function that shares the same period as C(xv0t) and κ(xv0t) in Eq. (6).

    By substituting Eqs. (7)–(9) into Eq. (6) and adopting the homogenization theory of differential equations [60], we can obtain the homogeneous equation with a constant coefficient that has the same solution as Eq. (6) [6162]:

    C˜Tt+V˜Txκ2˜Tx2S2˜Txt=0 (10)

    Compared with the one-dimensional diffusion equation, the spatiotemporal modulated thermal diffusion equation has two terms: the convection-like term V˜Tx and the thermal equivalent term S2˜Txt (to be discussed below). The convection-like term (V˜Tx) can break the spatial inversion symmetry, facilitating the realization of nonreciprocal thermal diffusion. The appearance of this convection-like term represents the effect of “medium bias” due to the spatial modulation of the parameters. Thus, the spatial direction preference of thermal diffusion will emerge [63]. Both the convection-like and convection terms can achieve nonreciprocal heat transfer. Compared with convection terms involving actual mechanical or physical flow [6465], the convection-like term resulting in nonreciprocal heat transfer includes only the “motion” of material parameters. The homogenization parameters in Eq. (10) are described as follows:

    κκ0[1(Δκ)2211+Γ2] (11)
    CC0[1(ΔC)22Γ21+Γ2] (12)
    Vv0κ0βΔCΔκ211+Γ2(cosδ+Γsinδ) (13)
    S1v0C0ΔCΔκ2βΓ21+Γ2(cosδ+Γ1sinδ) (14)

    where Γ=v0C0/βκ0 denotes the dimensionless modulation speed and ˜T refers to the envelope of the actual temperature wave.

    When the modulation speed Γ=0, the homogenization parameters V and S are both zero, and the homogenization material parameters show that the average volumetric heat capacity C=C0 and the equivalent thermal conductivity κ=κ0[1(Δκ)2], leading to reciprocal thermal diffusion. Similarly, reciprocal thermal diffusion can also be realized by setting the modulation speed Γ+ where homogenization parameters V and S approach zero and the homogenization material parameters show equivalent volumetric heat capacity C=C0[1(ΔC)2] and average thermal conductivity κ=κ0. However, when the modulation speed Γ0, a nonzero value of V will introduce a convective term, leading to a biased temperature distribution that exhibits nonreciprocity of thermal diffusion, as displayed in Fig. 2(a) [66].

    Figure  2.  Nonreciprocal thermal diffusion in spatiotemporal modulated materials: (a) nonreciprocal temperature distributions in spatiotemporal modulated metamaterials under forward (left) and backward (right) propagation. Γ is the space–time modulation moving speed. Γ = 0 denotes no modulation. The nonreciprocal effect shows an initial increase and then a decrease as the modulation speed varies. Reprinted figure with permission from [66] as follows: D. Torrent, O. Poncelet, and J.C. Batsale, Phys. Rev. Lett., 120, 125501 (2018). Copyright 2018 by the American Physical Society; (b) forward propagation of temperature waves in metamaterials (top) and two backward propagation scenarios: changing the source position but maintaining the modulation direction (middle) and changing the modulation direction but maintaining the source position (bottom). u is modulation speed and equivalent to the speed of movement v0; (c) phase difference δ measures the extent of nonreciprocity. When δ = –π/4, nonreciprocity disappears, and when δ = π/4, heat wave propagation shows nonreciprocal behavior. Reprinted figure with permission from [62] as follows: L.J. Xu, J.P. Huang, and X.P. Ouyang, Phys. Rev. E, 103, 032128 (2021). Copyright 2021 by the American Physical Society.

    It is clear that both C and κ are independent of the phase difference δ in Eqs. (11)–(14), while parameters V and S show dependence on the phase difference δ. The phase difference δ provides a flexible and tunable parameter to control the degree of nonreciprocity by affecting V and S. By manipulating the position of the heat source and the modulation direction of the parameters, the direction of heat flow transmission can possibly be controlled, as illustrated in Fig. 2(b) [62]. When keeping a fixed position for the heat source, modifying the direction of parameter modulation leads to a corresponding reversal in the direction of heat flow transmission. Conversely, if we exchange the position of the heat source while keeping the direction of parameter modulation unchanged, the opposite effect is observed in heat flow transmission. To depict the effect of phase difference δ on the degree of nonreciprocity, the forward and backward-1 cases (Fig. 2(b)) with the same modulation direction but opposite source position are taken as examples. When δ = –π/4, the evolution of temperature and heat flux (the first column of Fig. 2(c)) for these two cases is the same, demonstrating reciprocal propagation. Meanwhile, when δ = π/4, the temperature amplitudes and heat flux are different, demonstrating nonreciprocal propagation (the second column of Fig. 2(c)).

    In Eq. (10), we find that the homogenization equation has an additional term of S2˜Txt, which is the thermal equivalent term of the Willis coefficient in the elastic dynamics of an inhomogeneous medium [67]. Thermal Willis coupling in thermal diffusion is the interplay between heat flux and temperature change rate. The time-independent convection-like term (V(˜T/x)) can break the spatial inversion symmetry and realize asymmetric diffusion in both transient and quasisteady states. While the thermal equivalent term (S2˜Txt) is time-dependent and relies on the rate of temperature change, nonreciprocal thermal diffusion by Willis coupling takes place only in transient processes [63,68].

    A Willis thermal metamaterial, developed with spatiotemporal modulation, has precise control over the directional propagation of the temperature field. Its propagation direction can be reversed precisely around the critical point of the modulation speed (Fig. 3(a)) [69]. In porous media, the Fizeau drag in thermal diffusion can be observed in the Willis thermal metamaterial (Fig. 3(b)), indicating distinct propagation velocities of the temperature field in opposite directions and giving rise to a nonreciprocal process of thermal diffusion [70]. The spatiotemporal modulation-based Willis coupling not only realizes nonreciprocal thermal diffusion but also offers ideas for controlling nonequilibrium mass and energy transport.

    Figure  3.  Spatiotemporally modulated nonreciprocal thermal metamaterials: (a) asymmetric thermal diffusion in Willis thermal metamaterials induced by spatiotemporal modulation parameters. Reprinted figure with permission from [69] as follows: L.J. Xu, G.Q. Xu, J.X. Li, Y. Li, J.P. Huang, and C.W. Qiu, Phys. Rev. Lett., 129, 155901 (2022). Copyright 2022 by the American Physical Society; (b) thermal Fizeau drag resulting from the thermal Willis coupling mechanism. Reprinted figure with permission from [70] as follows: L.J. Xu, G.Q. Xu, J.P. Huang, and C.W. Qiu, Phys. Rev. Lett., 128, 145901 (2022). Copyright 2022 by the American Physical Society; (c) thermal wave diodes show rectification of the heat flux in forward and backward propagation (left). The amplitude profiles of the temperature field are shown in the center, and the right panel displays the rectification coefficient. The arrows and their amplitudes are the direction and magnitude of the heat wave, respectively. Reprinted figure with permission from [68] as follows: J. Ordonez-Miranda, Y.Y. Guo, J.J. Alvarado-Gil, S. Volz, and M. Nomura, Phys. Rev. Appl., 16, L041002 (2021). Copyright 2021 by the American Physical Society.

    The material’s thermal parameters with wave-like modulation that achieve a nonreciprocal heat transfer process can be determined by switching between low and high thermal conductivity and heat capacity states. When material parameters present sensitivity to external stimuli (such as magnetic fields, electric fields, and light), periodically arranged layered media can be produced by applying periodic external stimuli. For example, upon exposure to ultraviolet (375 nm) and green light (530 nm), the azobenzene polymer shows a rapid reversal of thermal conductivity at room temperature within seconds [71]. The layered structure, comprising p-azobenzene polymer layers, uses ultraviolet and green period-lighting to efficiently modulate thermal conductivity and volumetric heat capacity. On the basis of this method, a thermal diode with a rectification factor (R = ||QF(x=d)||QB(x=0)||max, where {Q_F}(x = d) and {Q_B}(x = 0) are the forward and backward heat fluxes at x = d and x = 0, respectively) above 86% has been attained (Fig. 3(c)) [68]. Compared with traditional nonlinear steady-state diodes, this diode is implemented by propagating heat waves driven by Willis coupling produced through spatiotemporal modulation of thermal parameters.

    Another way to realize nonreciprocity in linear systems is based on the Zeeman effect caused by angular momentum bias. For example, in a three-port annular cavity with rotating airflow, the circular motion of the airflow induces an angular momentum bias that results in the acoustic Zeeman effect, thereby allowing nonreciprocal acoustic behavior of the device [72]. In convection–diffusion systems, the Zeeman effect originates from the angular momentum bias induced by the volume force that acts on the fluid, which can be gravity, centrifugal force, or other forms of force [73].

    In a one-dimensional scenario with convection velocity v, the diffusion equation is

    \rho c\frac{{\partial T}}{{\partial t}} = \kappa \frac{{{\partial ^2}T}}{{\partial {x^2}}} - \rho cv\frac{{\partial T}}{{\partial x}} (15)

    where ρ, c, κ, and v are the temperature-independent mass density, specific heat capacity, thermal conductivity, and convective velocity of the fluid, respectively. Fluid refers to a steady, incompressible creeping flow between parallel plates, with its velocity demonstrating a parabolic distribution along the vertical direction of the flow. For a narrow plate of heighth, the convective velocity is represented by v = - {h^2}{{(\nabla P - f)} \mathord{\left/ {\vphantom {{(\nabla P - f)} {(12\mu )}}} \right. } {(12\mu )}} , where μ is the dynamic viscosity, P is the pressure, and f is the volume force applied to the fluid. As shown on the left side of Fig. 4(a), the fluid velocity + {v_0}( - {v_0}) is produced by the fluid pressure in the + x( - x) direction. The application of a volume force f to the fluid induces an additional convective velocity in the direction of the force. The fluid velocity in the + x( - x) direction is {v_ + }({v_ - }), while \left| {{v_ + }} \right| \ne \left| {{v_ - }} \right|.

    Figure  4.  Thermal nonreciprocity induced by angular momentum bias: (a) thermal Zeeman effect. Schematic of the frequency splitting principle (left) and degree of splitting of the real part of the frequency as a function of the volume force (right); (b) nonreciprocal heat transport due to the thermal Zeeman effect. An angular momentum bias ring exhibiting thermal nonreciprocity and isolation in which the fluid is water. Temperature distribution at 600 s under different volume forces (first row), temperature changes at ports 2 and 3 under different volume forces for 500 to 600 s (second row), and convective velocity distribution (third row). Reprinted with permission from L.J. Xu, J.P. Huang, and X.P. Ouyang, Appl. Phys. Lett., 118, 221902 (2021) [73]. Copyright 2021 AIP Publishing LLC; (c) topological protection in thermal convection–diffusion systems. Two thermal spin structures (first row) consisting of counterclockwise (red unit) and clockwise convection (blue unit), topological protection in the unidirectional hot edge state (second line), and unidirectional heat transfer at the topological thermal interface state (third and fourth lines). Reprinted from L.J. Xu and J.P. Huang, EPL Europhys. Lett., 134, 60001 (2021) [75]. © IOP Publishing. Reproduced with permission. All rights reserved.

    By applying a periodic temperature input aligned with the fluid flow direction T = {T_0} + A\cos (\beta x - \omega t), where {T_0}, A, β, and ω are the reference temperature, temperature amplitude, wavenumber, and frequency, respectively. By taking the counterclockwise direction as the + x direction (indicated by the blue arrow in Fig. 4(a)) and substituting the periodic temperature into the convection–diffusion equation, a frequency of \omega_0=\beta v_0-\mathrm{i}\beta^2D with the thermal diffusion coefficient D = {\kappa \mathord{\left/ {\vphantom {\kappa {(\rho c)}}} \right. } {(\rho c)}} is obtained for both clockwise and counterclockwise directions when no volume force exists. \mathrm{Re}(\omega_0) and -\mathrm{Im}(\omega_0) are the frequency and attenuation factor of temperature waves during propagation, respectively [74].

    If volume force f = 0, the counterclockwise and clockwise frequencies are no longer equal and can be mathematically expressed as follows:

    {\omega _ \pm } = \beta {v_ \pm } - {\text{i}}{\beta ^2}D (16)

    where convective velocities {v_ \pm } = {v_0} \pm {{{h^2}f} }/ {(12\mu )} . This phenomenon of frequency splitting caused by the angular momentum bias created by the volume force is called the thermal Zeeman effect, and the degree of frequency splitting increases with increasing volume force (the right picture in Fig. 4(a)).

    In Fig. 4(b), Xu et al. [73] demonstrated nonreciprocal heat transport due to the thermal Zeeman effect in a fluid-filled three-port ring cavity. The fluid flow within the cavity can be considered as the flow between the narrow plates. Port 1 acts as the input port (high voltage P_{\mathrm{H}} ), while ports 2 and 3 are output ports (low voltage {P_{\text{L}}}). A periodic temperature source is placed on port 1, with ports 2 and 3 set as open borders. For zero-volume forces, the ring exhibits two symmetric velocities ({v_{1 \to 2}} = {v_{1 \to 3}}) with identical temperature amplitudes at ports 2 and 3. When a counterclockwise volume force is applied, {v_{1 \to 2}} increases while port {v_{1 \to 3}} decreases, leading to an increase in temperature amplitude at port 2 and complete temperature isolation at port 3, which achieves the nonreciprocal propagation of heat waves and the temperature isolation of port 3. The utilization of thermal nonreciprocity caused by angular momentum bias can be utilized to realize topological protection [75], as depicted in Fig. 4(c). A thermal spin structure with counterclockwise or concurrent flow is arranged in a honeycomb lattice. When the heat source is positioned in the lower left corner, a unidirectional thermal edge state appears within the lattice [76]. This phenomenon arises because of the nonreciprocal thermal isolation induced by the angular momentum bias resulting from the spin of the fluid in the structure, allowing the lattice surface to only support temperature propagation in the same direction as the spin of the structural unit (the second row of Fig. 4(c)). Combining two types of thermal spin structures results in the formation of topological interface states at their interface [77]. Because of the distinct convection directions of the two thermal spin unit structures, the attenuation rate is minimized at the unique boundary surface with a matching convection direction. This allows the heat wave to propagate along the flow field while preventing its propagation along the countercurrent field, thus realizing nonreciprocal propagation of the heat wave (the third and fourth lines of Fig. 4(c)) [75]. The results presented here provide potential for examining the topological properties related to diffusion dynamics, in particular those related to topological heat.

    Besides the above nonreciprocal thermal transfer, recently, nonreciprocal thermal radiation is also a rapid development direction for nonreciprocal thermal metamaterials.

    Thermal radiation processes are governed by Kirchhoff’s law [78]:

    e(\omega ,\theta ,\phi ) = \alpha (\omega ,\theta ,\phi ) (17)

    where e and \alpha refer to the directional spectral emissivity and absorptivity, respectively, \omega refers to the frequency, and \phi and \theta refer to the azimuth angle and pitch angle, respectively. From Eq. (17), the structures or devices with high/low absorption also have high/low emission characteristics, which limits the energy conversion efficiency in solar energy harvesting [79] and radiative cooling [80]. Thus, researchers are striving to overcome this balance limitation to realize nonreciprocal emission/absorption.

    The essence of the thermal radiation process is that an object with a temperature above absolute zero radiates electromagnetic waves outward. This phenomenon follows the Lorentz reciprocity theorem, which acts as a case of the application of the Onsager reciprocal relations to electromagnetic processes [39]. Kirchhoff’s law, in accordance with the Lorentz reciprocity theorem, is applicable only to nonmagnetic, time-invariant, and linear materials. As a consequence, it becomes feasible to realize nonreciprocal emission and absorption by removing any of the three conditions, such as magnetic response, time-variant systems, or optical nonlinearity [81].

    One approach to breaking the Lorentz reciprocity is to utilize materials with an asymmetric permittivity tensor. For example, under an external magnetic field, the permittivity of a magneto-optical material is characterized by an asymmetric tensor. Among these materials, InAs, which is a prominent magneto-optical material, has been extensively applied to achieve heightened nonreciprocity. For example, Zhu et al. [78] first proposed a grating structure madding of InAs to show the violation of Kirchhoff’s law. The one-dimensional grating of n-InAs atop aluminum is shown at the top of Fig. 5(a). Under external magnetic field B = 3T in the z-direction, the absorptivity and emissivity no longer overlap, and Kirchhoff’s law of thermal radiation is almost completely broken (the bottom of Fig. 5(a)). Later, many structures based on InAs with different physical mechanisms have been introduced into thermal emitters to achieve strong nonreciprocal radiation [8285].

    Figure  5.  Nonreciprocal thermal radiation in different mechanisms: (a) grating of n-InAs atop aluminum (top). The structure is periodic in the x direction, and has the following geometry parameters: p = 7.24 μm, w = 3.2 μm, t1 = 1.981 μm, and t2 = 0.485 μm. Emissivity spectra (e) and absorptivity spectra (α) at θ = 162.8° (bottom). Reprinted figure with permission from [78] as follows: L.X. Zhu and S.H. Fan, Phys. Rev. B, 90, 220301 (2014). Copyright 2014 by the American Physical Society; (b) Weyl semimetallic photonic crystals for nonreciprocal thermal emission(top). Emissivity and absorptivity spectra in the θ = 80° direction [86]. Reprinted with permission from B. Zhao, C. Guo, C.A.C. Garcia, P. Narang, and S. Fan, Nano Lett., 20, 1923 (2020) [86]. Copyright 2020 American Chemical Society; (c) grating design and the corresponding unit cell (left) with period Λ = 1.7 μm. The dashed region undergoes a time-dependent modulation of the refractive index. Band structure of grating (middle) and unequal absorptivity and emissivity (right). Reprinted with permission from A. Ghanekar, J.H. Wang, S.H. Fan, and M.L. Povinelli, ACS Photonics, 9, 1157 (2022) [37]. Copyright 2022 American Chemical Society; (d) nonreciprocal thermal radiation based on Kerr nonlinearity. Schematic of the three resonators at equal temperature supporting modes (middle). Heat extraction power Hex from the emitter and unconverted power Hup compared with the temperature-driven radiative heat transfer between two vacuum-separated SiC plates held at 300 K (light blue) and 1 K (dark blue) temperature differences [108].

    However, the utilization of nonreciprocal thermal emitters fabricated from InAs is typically confined to a large applied magnetic field, which restricts its potential applications. Without an external magnetic field, magnetic Weyl semimetals (WSMs) can violate Lorentz reciprocity and time-reversal symmetry, which is due to the significant anomalous Hall effect, which is linked to an improved Berry curvature at the Weyl nodes [86]. The vector 2b that separates the two Weyl cones in momentum space functions similarly to an internal magnetic field, as shown at the top of Fig. 5(b). WSM-based emitters can break Kirchhoff’s law of thermal radiation even without a magnetic field (the bottom of Fig. 5(b)). Over the past few years, many photonic structures involving magnetic Weyl semimetals, including one-dimensional nanowire arrays [87] and planar structures [8889], have been developed to breach the reciprocity between absorbance and emittance.

    The above structures rely on the external/internal magnetic response and can be summarized as grating structures [8287,9092], thin film structures [89,93], and multilayer structures [94100], which refer to a number of physical effects, including guided mode resonance, surface plasmon resonance, Tamm plasmon resonance, Fabry–Perot resonance, topological interface state, and bound state in the continuum. Very recently, it is heartening to note that the experimental demonstration of strong nonreciprocal radiation has been realized [101103].

    Time modulation provides a magnet-free alternative to breaking the Lorentz reciprocity [3738]. Nonreciprocal emission and absorption based on time-variant systems have been experimentally demonstrated at radio frequencies [104], demonstrating the potential of spatiotemporal modulation methods for thermal management and energy harvesting at infrared frequencies. The left view of Fig. 5(c) presents a nonreciprocal infrared thermal emitter that depends on dynamic modulation using the spatiotemporal grating refractive index to drive photon transitions between guided resonance modes [37]. It includes a dielectric grating coated with a perfect electrical conductor on both the top and bottom surfaces, and a detail of the unit cell is shown below the grating. The region marked by the dashed line undergoes a spatiotemporal permittivity modulation represented by ε = εb + \Delta\varepsilon\cos(\mathit{\Omega}t+Kx) , where \varepsilon_{\mathrm{b}} refers to the base permittivity of the dielectric, \Delta \varepsilon refers to the amplitude of modulation, \mathit{\Omega} refers to the modulation frequency, and K refers to the spatial frequency of modulation. \tau=\Delta\varepsilon\mathord{\left/\vphantom{\Delta\varepsilon\varepsilon_b}\right.}\varepsilon_{\mathrm{b}} is the modulation depth. For zero modulation depth, the band structure of the grating plots is in the middle of Fig. 5(c). Modes excited by {S_{1 + }} and {S_{2 + }} are separated by modulation frequency Ω = 47.93 GHz and normalized wave vector K\mathit{\Lambda}/2\text{π }=\pm1/3 . The forward (modes excited by obliquely incident wave {S_{1 + }} can couple to outgoing waves {S_{2 - }} and {S_{3 - }}) and inverse (an outgoing wave {S_{1 - }} can result only from input wave {S_{3 + }}) problems are defined. This demonstrates that the grating exhibits perfect nonreciprocal reflection (the right of Fig. 5(c)). It turns out that the time modulation of the Fermi energy of graphene can be used to generate nonreciprocity [38]. These works offer an idea for the development of nonreciprocal thermal emitters using electro-optical devices. Time modulation can also be employed to achieve photon refrigeration [105] and heat radiation pumps [106107].

    The violation of Lorentz reciprocity can take place in nonlinear materials, as exemplified by the Kerr effect used in thermal radiation control [108109]. The Kerr effect, a typical nonlinear optical effect, originates from the three-boundary nonlinear coefficients of materials. When an intense laser beam impinges on a material, changes in its refractive index occur, which in turn influences the transmission characteristics of the laser.

    As shown in the middle of Fig. 5(d), a three-resonator system is considered [108], involving resonators that support modes at {\omega _1}, {\omega _3}, and the resonator with resonant frequency {\omega _2} = {{({\omega _3} - {\omega _1})} \mathord{\left/ {\vphantom {{({\omega _3} - {\omega _1})} 2}} \right. } 2} is excited by laser irradiance and placed in a nonlinear material with {\chi ^{(3)}} nonlinearity. Through a four-wave mixing process, photons in resonators 1 and 3 are coupled through interaction with excited photons in resonator 2. In such a four-wave mixing process, the total photon number is conserved, and the coupling between resonance modes 1 and 3 is nonreciprocal. The physical construction of such a scheme involves an indium tin oxide absorber (resonator 3) and a silicon carbide (SiC) emitter (resonator 1), with a nonlinear {\chi ^{(3)}} spacer in between (the left of Fig. 5(d)). Graphene nanosheets with fermi energy EF = 0.7 eV are employed to produce strong local plasmon resonances to improve the nonlinear response. Calculations based on coupled-mode theory show that heat is extracted from the SiC emitter (the right of Fig. 5(d)) even when the two substrates are at the same temperature. The heat extraction study confirmed the significant potential of refrigeration using Kerr nonlinearity.

    In the earlier sections, we extensively examined nonreciprocity phenomena in terms of the thermal field, as well as the design and applications of nonreciprocal thermal metamaterials. However, it is crucial to acknowledge that nonreciprocity phenomena extend beyond the thermal field. It also attracts a broader exploration of other physical fields governed by the Laplace equation, such as the electric charge diffusion process. This section is dedicated to elucidating the possibility of nonreciprocity in the thermal field extending to other Laplace physical fields.

    For heat, the relevant modulation physical quantities are heat capacity and thermal conductivity. For the diffusion of electric charge, the corresponding physical quantities are capacitance and electrical conductivity. Electric charge diffusion can be macroscopically expressed by Fick’s diffusion equation [63]:

    \frac{{\partial q(x,t)}}{{\partial t}} = D\frac{{{\partial ^2}}q(x,t)}{{\partial {x^2}}} (18)

    where D=\sigma g is the diffusion coefficient, \sigma and g are the conductivity and the inverse of the capacity of the medium, respectively, and q(x,t) is charge density. It can be seen that the electric charge diffusion equation (Eq. (18)) has a form similar to that of the heat conduction equation (Eqs. (4) and (6)). Thus, by introducing traveling-wave spatiotemporal modulation to dynamically change the material parameters such as conductivity and capacitance, a nonreciprocal electric charge diffusion process can be realized [63,110]. In Fig. 6(a), a constant voltage source is applied at the left end of the sample, while the right end is kept electrically insulated. The sample consists of 50 disk capacitors that rotate at a constant angular velocity to generate spatiotemporal variation in conductivity and capacitance, which induces an asymmetric transfer of charge.

    Figure  6.  Metamaterials in the electric charge and mass diffusion process: (a) experimental setup consisting of 50 discs (left), theoretical prediction of the time-average voltage for each disk at each rotational period (middle), and experimental measurement (right) [63]; (b) mass diffusion metamaterials that manipulate the behavior of ion diffusion can realize ion cloaking, concentration, and ion selection in liquid solvents [114].

    Besides charge diffusion, Fick’s diffusion equation also holds for the mass diffusion process, which is driven by the potential gradient. For instance, metamaterial-based design methods have been designed to realize mass diffusion invisibility cloaks [111114], concentrators [114], rotators, and material separation devices [115117]. Fig. 6(b) illustrates a device that integrates a mass diffusion metamaterial with different functions for manipulating ion diffusion in liquids, including a bilayer cloak, a concentrator, and an ion selector [114]. By adjusting the material parameters, it is feasible to realize nonreciprocal material diffusion, allowing the selective transmission of specific substances. Such developments find valuable applications in tasks such as separating mixtures and filtering contaminants, for example, in areas such as wastewater treatment, air filtration, and fractionation.

    This study presents a comprehensive overview of nonreciprocal thermal transfer. First, we present an elaborate explanation of the Onsager reciprocity relation, which is a basic theory in the design of nonreciprocal thermal metamaterials. By considering four key assumptions, we derive the Onsager reciprocal relations and show their application to heat conduction processes. Next, we explore different approaches to breaking reciprocity and extensively discuss three commonly applied methods: nonlinear systems, spatiotemporal modulation, and angular momentum bias.

    Nonlinear thermal materials or structures can break the linear response and realize nonreciprocal heat transfer. However, its application is limited by temperature-dependent material parameters. By adjusting the material parameter properties, spatiotemporal modulation induces the breaking of both temporal and spatial inversion symmetry in linear systems. Angular momentum bias leads to a difference in the distribution of angular momentum by introducing rotation, which realizes nonreciprocal heat transfer. These methods give us a higher degree of freedom of thermal field regulation so that thermal metamaterials exhibit several novel applications, such as thermal isolation [73], geometric heat pumps [35], and topological protection [75].

    Future research directions are outlined. First, it is imperative to perform extensive studies on the design and preparation methods of nonreciprocal thermal metamaterial devices to address real-life thermal regulation challenges. For example, the development of nonreciprocal thermal radiation devices should progress from the initial single-channel to dual-channel and multichannel systems. To address the demands of more intricate operating environments, dynamically tunable nonreciprocal radiation should also be considered. Integration with other material design approaches should be considered, such as leveraging machine learning for parametric optimization [118]. In addition, synergistic combinations of nonreciprocal thermal metamaterials with advanced materials and technologies such as nanotechnology [119] and photonics [81] hold promise for realizing improved levels of thermal field regulation and application. Finally, it is crucial to examine the potential application of the design methodology applied in nonreciprocal thermal metamaterials to other Laplacian fields, including current and material diffusion fields. To sum up, investigating nonreciprocity in thermal metamaterials provides innovative insights and methodologies for controlling and regulating heat transfer.

    This work was supported by grants from the National Natural Science Foundation of China (No. 52325208), the Fundamental Research Funds for the Central Universities (No. 06500174), and National Key Research and Development Program of China (No. 2022YFB3807401).

  • [1]

    C.Z. Fan, Y. Gao, and J.P. Huang, Shaped graded materials with an apparent negative thermal conductivity, Appl. Phys. Lett., 92(2008), No. 25, art. No. 251907. DOI: 10.1063/1.2951600

    [2]

    T. Chen, C.N. Weng, and J.S. Chen, Cloak for curvilinearly anisotropic media in conduction, Appl. Phys. Lett., 93(2008), No. 11, art. No. 114103. DOI: 10.1063/1.2988181

    [3]

    U. Leonhardt, Optical conformal mapping, Science, 312(2006), No. 5781, p. 1777. DOI: 10.1126/science.1126493

    [4]

    J.B. Pendry, D. Schurig, and D.R. Smith, Controlling electromagnetic fields, Science, 312(2006), No. 5781, p. 1780. DOI: 10.1126/science.1125907

    [5]

    S.A. Cummer, B.I. Popa, D. Schurig, D.R. Smith, and J. Pendry, Full-wave simulations of electromagnetic cloaking structures, Phys. Rev. E, 74(2006), No. 3, art. No. 036621. DOI: 10.1103/PhysRevE.74.036621

    [6]

    D. Schurig, J.J. Mock, B.J. Justice, et al., Metamaterial electromagnetic cloak at microwave frequencies, Science, 314(2006), No. 5801, p. 977. DOI: 10.1126/science.1133628

    [7]

    W.S. Cai, U.K. Chettiar, A.V. Kildishev, and V.M. Shalaev, Optical cloaking with metamaterials, Nat. Photonics, 1(2007), p. 224. DOI: 10.1038/nphoton.2007.28

    [8]

    G.X. Zheng, H. Mühlenbernd, M. Kenney, G.X. Li, T. Zentgraf, and S. Zhang, Metasurface holograms reaching 80% efficiency, Nat. Nanotechnol., 10(2015), p. 308. DOI: 10.1038/nnano.2015.2

    [9]

    W.T. Chen, A.Y. Zhu, V. Sanjeev, et al., A broadband achromatic metalens for focusing and imaging in the visible, Nat. Nanotechnol., 13(2018), p. 220. DOI: 10.1038/s41565-017-0034-6

    [10]

    S.M. Wang, P.C. Wu, V.C. Su, et al., A broadband achromatic metalens in the visible, Nat. Nanotechnol., 13(2018), p. 227. DOI: 10.1038/s41565-017-0052-4

    [11]

    R.A. Shelby, D.R. Smith, and S. Schultz, Experimental verification of a negative index of refraction, Science, 292(2001), No. 5514, p. 77. DOI: 10.1126/science.1058847

    [12]

    S.H. Lee, C.M. Park, Y.M. Seo, and C.K. Kim, Reversed Doppler effect in double negative metamaterials, Phys. Rev. B, 81(2010), No. 24, art. No. 241102. DOI: 10.1103/PhysRevB.81.241102

    [13]

    A. Alù and N. Engheta, Achieving transparency with plasmonic and metamaterial coatings, Phy. Rev. E, 72(2005), art. No. 016623. DOI: 10.1103/PhysRevE.72.016623

    [14]

    S. Yang, L.J. Xu, G.L. Dai, and J.P. Huang, Omnithermal metamaterials switchable between transparency and cloaking, J. Appl. Phys., 128(2020), No. 9, art. No. 095102. DOI: 10.1063/5.0013270

    [15]

    T. Qu, J. Wang, and J.P. Huang, Manipulating thermoelectric fields with bilayer schemes beyond Laplacian metamaterials, EPL Europhys. Lett., 135(2021), No. 5, art. No. 54004. DOI: 10.1209/0295-5075/ac1648

    [16]

    S. Narayana and Y. Sato, Heat flux manipulation with engineered thermal materials, Phys. Rev. Lett., 108(2012), No. 21, art. No. 214303. DOI: 10.1103/PhysRevLett.108.214303

    [17]

    S. Guenneau, C. Amra, and D. Veynante, Transformation thermodynamics: Cloaking and concentrating heat flux, Opt. Express, 20(2012), No. 7, p. 8207. DOI: 10.1364/OE.20.008207

    [18]

    R. Schittny, M. Kadic, S. Guenneau, and M. Wegener, Experiments on transformation thermodynamics: Molding the flow of heat, Phys. Rev. Lett., 110(2013), No. 19, art. No. 195901. DOI: 10.1103/PhysRevLett.110.195901

    [19]

    S. Narayana, S. Savo, and Y. Sato, Transient heat flux shielding using thermal metamaterials, Appl. Phys. Lett., 102(2013), No. 20, art. No. 201904. DOI: 10.1063/1.4807744

    [20]

    H.Y. Xu, X.H. Shi, F. Gao, H.D. Sun, and B.L. Zhang, Ultrathin three-dimensional thermal cloak, Phys. Rev. Lett., 112(2014), No. 5, art. No. 054301. DOI: 10.1103/PhysRevLett.112.054301

    [21]

    T.C. Han, X. Bai, D.L. Gao, J.T.L. Thong, B.W. Li, and C.W. Qiu, Experimental demonstration of a bilayer thermal cloak, Phys. Rev. Lett., 112(2014), No. 5, art. No. 054302. DOI: 10.1103/PhysRevLett.112.054302

    [22]

    D.M. Nguyen, H.Y. Xu, Y.M. Zhang, and B.L. Zhang, Active thermal cloak, Appl. Phys. Lett., 107(2015), No. 12, art. No. 121901. DOI: 10.1063/1.4930989

    [23]

    T.C. Han, J.J. Zhao, T. Yuan, D.Y. Lei, B.W. Li, and C.W. Qiu, Theoretical realization of an ultra-efficient thermal-energy harvesting cell made of natural materials, Energy Environ. Sci., 6(2013), No. 12, p. 3537. DOI: 10.1039/c3ee41512k

    [24]

    C. Fei and Y.L. Dang, Experimental realization of extreme heat flux concentration with easy-to-make thermal metamaterials, Sci. Rep., 5(2015), art. No. 11552. DOI: 10.1038/srep11552

    [25]

    W. Liu, C. Lan, M. Ji, et al., A flower-shaped thermal energy harvester made by metamaterials, Global Challenges, 1(2017), No. 6, art. No. 1700017. DOI: 10.1002/gch2.201700017

    [26]

    S. Guenneau and C. Amra, Anisotropic conductivity rotates heat fluxes in transient regimes, Opt. Express, 21(2013), No. 5, p. 6578. DOI: 10.1364/OE.21.006578

    [27]

    F.B. Yang, B.Y. Tian, L.J. Xu, and J.P. Huang, Experimental demonstration of thermal chameleonlike rotators with transformation-invariant metamaterials, Phys. Rev. Appl., 14(2020), No. 5, art. No. 054024. DOI: 10.1103/PhysRevApplied.14.054024

    [28]

    T.C. Han, X. Bai, J.T.L. Thong, B.W. Li, and C.W. Qiu, Full control and manipulation of heat signatures: Cloaking, camouflage and thermal metamaterials, Adv. Mater., 26(2014), No. 11, p. 1731. DOI: 10.1002/adma.201304448

    [29]

    T.Z. Yang, Y.S. Su, W.K. Xu, and X.D. Yang, Transient thermal camouflage and heat signature control, Appl. Phys. Lett., 109(2016), No. 12, art. No. 121905. DOI: 10.1063/1.4963095

    [30]

    S. Hong, S. Shin, and R.K. Chen, An adaptive and wearable thermal camouflage device, Adv. Funct. Mater., 30(2020), No. 11, art. No. 1909788. DOI: 10.1002/adfm.201909788

    [31]

    R. Hu, W. Xi, Y.D. Liu, et al., Thermal camouflaging metamaterials, Mater. Today, 45(2021), p. 120. DOI: 10.1016/j.mattod.2020.11.013

    [32]

    R. Hu, S.Y. Huang, M. Wang, L.L. Zhou, X.Y. Peng, and X.B. Luo, Binary thermal encoding by energy shielding and harvesting units, Phys. Rev. Appl., 10(2018), No. 5, art. No. 054032. DOI: 10.1103/PhysRevApplied.10.054032

    [33]

    M. Lei, C.R. Jiang, F.B. Yang, J. Wang, and J.P. Huang, Programmable all-thermal encoding with metamaterials, Int. J. Heat Mass Transf., 207(2023), art. No. 124033. DOI: 10.1016/j.ijheatmasstransfer.2023.124033

    [34]

    M. Kasprzak, M. Sledzinska, K. Zaleski, et al., High-temperature silicon thermal diode and switch, Nano Energy, 78(2020), art. No. 105261. DOI: 10.1016/j.nanoen.2020.105261

    [35]

    Z. Wang, J. Chen, and J. Ren, Geometric heat pump and no-go restrictions of nonreciprocity in modulated thermal diffusion, Phys. Rev. E, 106(2022), No. 3, art. No. L032102. DOI: 10.1103/PhysRevE.106.L032102

    [36]

    Z.J. Coppens and J.G. Valentine, Spatial and temporal modulation of thermal emission, Adv. Mater., 29(2017), No. 39, art. No. 1701275. DOI: 10.1002/adma.201701275

    [37]

    A. Ghanekar, J.H. Wang, S.H. Fan, and M.L. Povinelli, Violation of Kirchhoff’s law of thermal radiation with space–time modulated grating, ACS Photonics, 9(2022), No. 4, p. 1157. DOI: 10.1021/acsphotonics.1c01350

    [38]

    A. Ghanekar, J.H. Wang, C. Guo, S.H. Fan, and M.L. Povinelli, Nonreciprocal thermal emission using spatiotemporal modulation of graphene, ACS Photonics, 10(2022), No. 1, p. 170.

    [39]

    V.S. Asadchy, M.S. Mirmoosa, A. Díaz-Rubio, S.H. Fan, and S.A. Tretyakov, Tutorial on electromagnetic nonreciprocity and its origins, Proc. IEEE, 108(2020), No. 10, p. 1684. DOI: 10.1109/JPROC.2020.3012381

    [40]

    L.L. Zhou, S.Y. Huang, M. Wang, R. Hu, and X.B. Luo, While rotating while cloaking, Phys. Lett. A, 383(2019), No. 8, p. 759. DOI: 10.1016/j.physleta.2018.11.041

    [41]

    L.J. Xu and J.P. Huang, Negative thermal transport in conduction and advection, Chin. Phys. Lett., 37(2020), No. 8, art. No. 080502. DOI: 10.1088/0256-307X/37/8/080502

    [42]

    X.Y. Huang, C.C. Lu, C. Liang, H.G. Tao, and Y.C. Liu, Loss-induced nonreciprocity, Light. Sci. Appl., 10(2021), No. 1, art. No. 30. DOI: 10.1038/s41377-021-00464-2

    [43]

    X.Y. Huang and Y.C. Liu, Perfect nonreciprocity by loss engineering, Phys. Rev. A, 107(2023), No. 2, art. No. 023703. DOI: 10.1103/PhysRevA.107.023703

    [44]

    Y.S. Su, Y. Li, M.H. Qi, S. Guenneau, H.G. Li, and J. Xiong, Asymmetric heat transfer with linear conductive metamaterials, Phys. Rev. Appl., 20(2023), No. 3, art. No. 034013. DOI: 10.1103/PhysRevApplied.20.034013

    [45]

    W. Muschik, Fundamentals of Nonequilibrium Thermodynamics, , Springer, Vienna, 1993.

    [46]

    W.J. Mansur, C.A.B. Vasconcellos, N.J.M. Zambrozuski, and O.C. Rotunno Filho, Numerical solution for the linear transient heat conduction equation using an explicit Green’s approach, Int. J. Heat Mass Transf., 52(2009), No. 3-4, p. 694. DOI: 10.1016/j.ijheatmasstransfer.2008.07.036

    [47]

    T.M. Chen, A hybrid Green’s function method for the hyperbolic heat conduction problems, Int. J. Heat Mass Transf., 52(2009), No. 19-20, p. 4273. DOI: 10.1016/j.ijheatmasstransfer.2009.04.026

    [48]

    M. Leindl, E.R. Oberaigner, and T. Antretter, Solution of a time-dependent heat conduction problem by an integral-equation approach, Comput. Mater. Sci., 52(2012), No. 1, p. 178. DOI: 10.1016/j.commatsci.2011.07.033

    [49]

    A. Mandelis, Diffusion-wave Fields : Mathematical Methods and Green Functions, Springer Science & Business Media, Berlin, 2013.

    [50]

    Y. Li, J.X. Li, M.H. Qi, C.W. Qiu, and H.S. Chen, Diffusive nonreciprocity and thermal diode, Phys. Rev. B, 103(2021), No. 1, art. No. 014307. DOI: 10.1103/PhysRevB.103.014307

    [51]

    G. Wehmeyer, T. Yabuki, C. Monachon, J.Q. Wu, and C. Dames, Thermal diodes, regulators, and switches: Physical mechanisms and potential applications, Appl. Phys. Rev., 4(2017), No. 4, art. No. 041304. DOI: 10.1063/1.5001072

    [52]

    C.W. Chang, D. Okawa, A. Majumdar, and A. Zettl, Solid-state thermal rectifier, Science, 314(2006), No. 5802, p. 1121. DOI: 10.1126/science.1132898

    [53]

    D.B. Go and M. Sen, On the condition for thermal rectification using bulk materials, J. Heat Transf., 132(2010), No. 12, art. No. 1.

    [54]

    Y. Li, X. Shen, Z. Wu, et al., Temperature-dependent transformation thermotics: From switchable thermal cloaks to macroscopic thermal diodes, Phys. Rev. Lett., 115(2015), No. 19, art. No. 195503. DOI: 10.1103/PhysRevLett.115.195503

    [55]

    X. Shen, Y. Li, C. Jiang, and J. Huang, Temperature trapping: Energy-free maintenance of constant temperatures as ambient temperature gradients change, Phys. Rev. Lett., 117(2016), No. 5, art. No. 055501. DOI: 10.1103/PhysRevLett.117.055501

    [56]

    J. Wang, J. Shang, and J.P. Huang, Negative energy consumption of thermostats at ambient temperature: Electricity generation with zero energy maintenance, Phys. Rev. Appl., 11(2019), No. 2, art. No. 024053. DOI: 10.1103/PhysRevApplied.11.024053

    [57]

    S.D. Lubner, J. Choi, G. Wehmeyer, et al., Reusable bi-directional 3ω sensor to measure thermal conductivity of 100-μm thick biological tissues, Rev. Sci. Instrum., 86(2015), No. 1, art. No. 014905. DOI: 10.1063/1.4905680

    [58]

    R.T. Zheng, J.W. Gao, J.J. Wang, and G. Chen, Reversible temperature regulation of electrical and thermal conductivity using liquid–solid phase transitions, Nat. Commun., 2(2011), art. No. 289. DOI: 10.1038/ncomms1288

    [59]

    J.X. Li, Y. Li, P.C. Cao, et al., Reciprocity of thermal diffusion in time-modulated systems, Nat. Commun., 13(2022), art. No. 167. DOI: 10.1038/s41467-021-27903-3

    [60]

    M.A. Biot, Thermoelasticity and irreversible thermodynamics, J. Appl. Phys., 27(1956), No. 3, p. 240. DOI: 10.1063/1.1722351

    [61]

    P.A. Huidobro, M.G. Silveirinha, E. Galiffi, and J.B. Pendry, Homogenization theory of space-time metamaterials, Phys. Rev. Appl., 16(2021), No. 1, art. No. 014044. DOI: 10.1103/PhysRevApplied.16.014044

    [62]

    L.J. Xu, J.P. Huang, and X.P. Ouyang, Tunable thermal wave nonreciprocity by spatiotemporal modulation, Phys. Rev. E, 103(2021), No. 3, art. No. 032128. DOI: 10.1103/PhysRevE.103.032128

    [63]

    M. Camacho, B. Edwards, and N. Engheta, Achieving asymmetry and trapping in diffusion with spatiotemporal metamaterials, Nat. Commun., 11(2020), art. No. 3733. DOI: 10.1038/s41467-020-17550-5

    [64]

    J. Li, Y. Li, P.C. Cao, et al., A continuously tunable solid-like convective thermal metadevice on the reciprocal line, Adv. Mater., 32(2020), No. 42, art. No. e2003823. DOI: 10.1002/adma.202003823

    [65]

    J. Li, Y. Li, W. Wang, L. Li, and C.W. Qiu, Effective medium theory for thermal scattering off rotating structures, Opt. Express, 28(2020), No. 18, p. 25894. DOI: 10.1364/OE.399799

    [66]

    D. Torrent, O. Poncelet, and J.C. Batsale, Nonreciprocal thermal material by spatiotemporal modulation, Phys. Rev. Lett., 120(2018), No. 12, art. No. 125501. DOI: 10.1103/PhysRevLett.120.125501

    [67]

    A.N. No. ris, A.L. Shuvalov, and A.A. Kutsenko, Analytical formulation of three-dimensional dynamic homogenization for periodic elastic systems, Proc. R. Soc. A, 468(2012), No. 2142, p. 1629. DOI: 10.1098/rspa.2011.0698

    [68]

    J. Ordonez-Miranda, Y.Y. Guo, J.J. Alvarado-Gil, S. Volz, and M. Nomura, Thermal-wave diode, Phys. Rev. Appl., 16(2021), No. 4, art. No. L041002. DOI: 10.1103/PhysRevApplied.16.L041002

    [69]

    L.J. Xu, G.Q. Xu, J.X. Li, Y. Li, J.P. Huang, and C.W. Qiu, Thermal Willis coupling in spatiotemporal diffusive metamaterials, Phys. Rev. Lett., 129(2022), No. 15, art. No. 155901. DOI: 10.1103/PhysRevLett.129.155901

    [70]

    L.J. Xu, G.Q. Xu, J.P. Huang, and C.W. Qiu, Diffusive fizeau drag in spatiotemporal thermal metamaterials, Phys. Rev. Lett., 128(2022), No. 14, art. No. 145901. DOI: 10.1103/PhysRevLett.128.145901

    [71]

    S.D. Sun, S.F. Liang, W.C. Xu, G.F. Xu, and S. Wu, Photoresponsive polymers with multi-azobenzene groups, Polym. Chem., 10(2019), No. 32, p. 4389. DOI: 10.1039/C9PY00793H

    [72]

    R. Fleury, D.L. Sounas, C.F. Sieck, M.R. Haberman, and A. Alù, Sound isolation and giant linear nonreciprocity in a compact acoustic circulator, Science, 343(2014), No. 6170, p. 516. DOI: 10.1126/science.1246957

    [73]

    L.J. Xu, J.P. Huang, and X.P. Ouyang, Nonreciprocity and isolation induced by an angular momentum bias in convection-diffusion systems, Appl. Phys. Lett., 118(2021), No. 22, art. No. 221902. DOI: 10.1063/5.0049774

    [74]

    Y. Li, Y.G. Peng, L. Han, et al., Anti-parity-time symmetry in diffusive systems, Science, 364(2019), No. 6436, p. 170. DOI: 10.1126/science.aaw6259

    [75]

    L.J. Xu and J.P. Huang, Robust one-way edge state in convection-diffusion systems, EPL Europhys. Lett., 134(2021), No. 6, art. No. 60001. DOI: 10.1209/0295-5075/134/60001

    [76]

    H. Hu, S. Han, Y. Yang, et al., Observation of topological edge states in thermal diffusion, Adv. Mater., 34(2022), No. 31, art. No. e2202257. DOI: 10.1002/adma.202202257

    [77]

    Y. Hatsugai, Edge states in the integer quantum Hall effect and the Riemann surface of the Bloch function, Phys. Rev. B: Condens. Matter., 48(1993), No. 16, p. 11851. DOI: 10.1103/PhysRevB.48.11851

    [78]

    L.X. Zhu and S.H. Fan, Near-complete violation of detailed balance in thermal radiation, Phys. Rev. B, 90(2014), No. 22, art. No. 220301. DOI: 10.1103/PhysRevB.90.220301

    [79]

    Z. Chen, L.X. Zhu, W. Li, and S.H. Fan, Simultaneously and synergistically harvest energy from the Sun and outer space, Joule, 3(2019), No. 1, p. 101. DOI: 10.1016/j.joule.2018.10.009

    [80]

    L.X. Zhu, A.P. Raman, and S.H. Fan, Radiative cooling of solar absorbers using a visibly transparent photonic crystal thermal blackbody, Proc. Natl. Acad. Sci. USA, 112(2015), No. 40, p. 12282. DOI: 10.1073/pnas.1509453112

    [81]

    Z.N. Zhang and L.X. Zhu, Nonreciprocal thermal photonics for energy conversion and radiative heat transfer, Phys. Rev. Appl., 18(2022), No. 2, art. No. 027001. DOI: 10.1103/PhysRevApplied.18.027001

    [82]

    J. Wu, F. Wu, T.C. Zhao, and X.H. Wu, Tunable nonreciprocal thermal emitter based on metal grating and graphene, Int. J. Therm. Sci., 172(2022), art. No. 107316. DOI: 10.1016/j.ijthermalsci.2021.107316

    [83]

    J. Wu and Y.M. Qing, The enhancement of nonreciprocal radiation for light near to normal incidence with double-layer grating, Adv. Compos. Hybrid Mater., 6(2023), No. 3, art. No. 87. DOI: 10.1007/s42114-023-00671-y

    [84]

    J. Wu and Y.M. Qing, Strong nonreciprocal radiation for extreme small incident angle, Int. Commun. Heat Mass Transf., 144(2023), art. No. 106794. DOI: 10.1016/j.icheatmasstransfer.2023.106794

    [85]

    J. Wu and Y.M. Qing, Near-perfect nonreciprocal radiation for extremely small incident angle based on cascaded grating structure, Int. J. Therm. Sci., 190(2023), art. No. 108340. DOI: 10.1016/j.ijthermalsci.2023.108340

    [86]

    B. Zhao, C. Guo, C.A.C. Garcia, P. Narang, and S. Fan, Axion-field-enabled nonreciprocal thermal radiation in Weyl semimetals, Nano Lett., 20(2020), No. 3, p. 1923. DOI: 10.1021/acs.nanolett.9b05179

    [87]

    Y. Tsurimaki, X. Qian, S. Pajovic, F. Han, M.D. Li, and G. Chen, Large nonreciprocal absorption and emission of radiation in type-I Weyl semimetals with time reversal symmetry breaking, Phys. Rev. B, 101(2020), No. 16, art. No. 165426. DOI: 10.1103/PhysRevB.101.165426

    [88]

    X.H. Wu, H.Y. Yu, F. Wu, and B.Y. Wu, Enhanced nonreciprocal radiation in Weyl semimetals by attenuated total reflection, AIP Adv., 11(2021), No. 7, art. No. 075106. DOI: 10.1063/5.0055418

    [89]

    J. Wu, Z.M. Wang, H. Zhai, Z.X. Shi, X.H. Wu, and F. Wu, Near-complete violation of Kirchhoff’s law of thermal radiation in ultrathin magnetic Weyl semimetal films, Opt. Mater. Express, 11(2021), No. 12, art. No. 4058. DOI: 10.1364/OME.444308

    [90]

    J. Wu and Y.M. Qing, Nonreciprocal thermal emitter for near perpendicular incident light with cascade grating involving weyl semimetal, Mater. Today Phys., 32(2023), art. No. 101025. DOI: 10.1016/j.mtphys.2023.101025

    [91]

    J. Wu, Y.S. Sun, B.Y. Wu, Z.M. Wang, and X.H. Wu, Extremely wide-angle nonreciprocal thermal emitters based on Weyl semimetals with dielectric grating structure, Case Stud. Therm. Eng., 40(2022), art. No. 102566. DOI: 10.1016/j.csite.2022.102566

    [92]

    J. Wu and Y.M. Qing, Tunable near-perfect nonreciprocal radiation with a Weyl semimetal and graphene, Phys. Chem. Chem. Phys., 25(2023), No. 13, p. 9586. DOI: 10.1039/D2CP05945B

    [93]

    M.Q. Liu, C. Zhao, Y.X. Zeng, Y. Chen, C.Y. Zhao, and C.W. Qiu, Evolution and nonreciprocity of loss-induced topological phase singularity pairs, Phys. Rev. Lett., 127(2021), No. 26, art. No. 266101. DOI: 10.1103/PhysRevLett.127.266101

    [94]

    J. Wu and Y.M. Qing, Strong nonreciprocal radiation with topological photonic crystal heterostructure, Appl. Phys. Lett., 121(2022), No. 11, art. No. 112101. DOI: 10.1063/5.0107022

    [95]

    J. Wu, Z.M. Wang, B.Y. Wu, Z.X. Shi, and X.H. Wu, The giant enhancement of nonreciprocal radiation in Thue-morse aperiodic structures, Opt. Laser Technol., 152(2022), art. No. 108138. DOI: 10.1016/j.optlastec.2022.108138

    [96]

    J. Wu, F. Wu, T.C. Zhao, M. Antezza, and X.H. Wu, Dual-band nonreciprocal thermal radiation by coupling optical Tamm states in magnetophotonic multilayers, Int. J. Therm. Sci., 175(2022), art. No. 107457. DOI: 10.1016/j.ijthermalsci.2022.107457

    [97]

    J. Wu and Y.M. Qing, Strong multi-band nonreciprocal radiation with Fibonacci multilayer involving Weyl semimetal, Results Phys., 51(2023), art. No. 106642. DOI: 10.1016/j.rinp.2023.106642

    [98]

    J. Wu and Y.M. Qing, Multichannel nonreciprocal thermal radiation with Weyl semimetal and photonic crystal heterostructure, Case Stud. Therm. Eng., 48(2023), art. No. 103161. DOI: 10.1016/j.csite.2023.103161

    [99]

    J. Wu and Y.M. Qing, A multi-band nonreciprocal thermal emitter involving a Weyl semimetal with a Thue–Morse multilayer, Phys. Chem. Chem. Phys., 25(2023), No. 16, p. 11477. DOI: 10.1039/D3CP00492A

    [100]

    J. Wu, B.Y. Wu, Z.M. Wang, and X.H. Wu, The enhanced nonreciprocal radiation with topological interface states, Opt. Laser Technol., 158(2023), art. No. 108907. DOI: 10.1016/j.optlastec.2022.108907

    [101]

    K.J. Shayegan, B. Zhao, Y. Kim, S. Fan, and H.A. Atwater, Nonreciprocal infrared absorption via resonant magneto-optical coupling to InAs, Sci. Adv., 8(2022), No. 18, art. No. eabm4308. DOI: 10.1126/sciadv.abm4308

    [102]

    M.Q. Liu, S. Xia, W.J. Wan, et al., Broadband mid-infrared non-reciprocal absorption using magnetized gradient epsilon-near-zero thin films, Nat. Mater., 22(2023), p. 1196. DOI: 10.1038/s41563-023-01635-9

    [103]

    K.J. Shayegan, S. Biswas, B. Zhao, S.H. Fan, and H.A. Atwater, Direct observation of the violation of Kirchhoff’s law of thermal radiation, Nat. Photonics, 17(2023), p. 891. DOI: 10.1038/s41566-023-01261-6

    [104]

    Y. Hadad, J.C. Soric, and A. Alu, Breaking temporal symmetries for emission and absorption, Proc. Natl. Acad. Sci. USA, 113(2016), No. 13, p. 3471. DOI: 10.1073/pnas.1517363113

    [105]

    S. Buddhiraju, W. Li, and S.H. Fan, Photonic refrigeration from time-modulated thermal emission, Phys. Rev. Lett., 124(2020), No. 7, art. No. 077402. DOI: 10.1103/PhysRevLett.124.077402

    [106]

    L.J. Fernández-Alcázar, H.N. Li, M. Nafari, and T. Kottos, Implementation of optimal thermal radiation pumps using adiabatically modulated photonic cavities, ACS Photonics, 8(2021), No. 10, p. 2973. DOI: 10.1021/acsphotonics.1c00896

    [107]

    H.N. Li, L.J. Fernández-Alcázar, F. Ellis, B. Shapiro, and T. Kottos, Adiabatic thermal radiation pumps for thermal photonics, Phys. Rev. Lett., 123(2019), No. 16, art. No. 165901. DOI: 10.1103/PhysRevLett.123.165901

    [108]

    C. Khandekar and A.W. Rodriguez, Near-field thermal upconversion and energy transfer through a Kerr medium, Opt. Express, 25(2017), No. 19, p. 23164. DOI: 10.1364/OE.25.023164

    [109]

    C. Khandekar, R. Messina, and A.W. Rodriguez, Near-field refrigeration and tunable heat exchange through four-wave mixing, AIP Adv., 8(2018), No. 5, art. No. 055029. DOI: 10.1063/1.5018734

    [110]

    J. Li, Z. Zhang, G. Xu, et al., Tunable rectification of diffusion-wave fields by spatiotemporal metamaterials, Phys. Rev. Lett., 129(2022), No. 25, art. No. 256601. DOI: 10.1103/PhysRevLett.129.256601

    [111]

    L.W. Zeng and R.X. Song, Controlling chloride ions diffusion in concrete, Sci. Rep., 3(2013), art. No. 3359. DOI: 10.1038/srep03359

    [112]

    Y. Li, C.B. Liu, Y. Bai, L.J. Qiao, and J. Zhou, Ultrathin hydrogen diffusion cloak, Adv. Theory Simul., 1(2018), No. 1, art. No. 1700004. DOI: 10.1002/adts.201700004

    [113]

    Y. Li, C.B. Liu, P. Li, et al., Scattering cancellation by a monolayer cloak in oxide dispersion-strengthened alloys, Adv. Funct. Mater., 30(2020), No. 36, art. No. 2003270. DOI: 10.1002/adfm.202003270

    [114]

    Y. Li, C.Y. Yu, C.B. Liu, et al., Mass diffusion metamaterials with “plug and switch” modules for ion cloaking, concentrating, and selection: Design and experiments, Adv. Sci., 9(2022), No. 30, art. No. 2201032. DOI: 10.1002/advs.202201032

    [115]

    J.M. Restrepo-Flórez and M. Maldovan, Mass separation by metamaterials, Sci. Rep., 6(2016), art. No. 21971. DOI: 10.1038/srep21971

    [116]

    X. Zhou, G.Q. Xu, and H.Y. Zhang, Binary masses manipulation with composite bilayer metamaterial, Compos. Struct., 267(2021), art. No. 113866. DOI: 10.1016/j.compstruct.2021.113866

    [117]

    Z. Zhang, L. Xu, and J. Huang, Controlling chemical waves by transforming transient mass transfer, Adv. Theor. Simul., 5(2022), No. 3, art. No. 2100375. DOI: 10.1002/adts.202100375

    [118]

    B. Liu, L.J. Xu, and J.P. Huang, Thermal transparency with periodic particle distribution: A machine learning approach, J. Appl. Phys., 129(2021), No. 6, art. No. 065101. DOI: 10.1063/5.0039002

    [119]

    S.A.M. Loos, S. Arabha, A. Rajabpour, A. Hassanali, and É. Roldán, Nonreciprocal forces enable cold-to-hot heat transfer between nanoparticles, Sci. Rep., 13(2023), No. 1, art. No. 4517. DOI: 10.1038/s41598-023-31583-y

Relative Articles

Lei Li, Zhi Wang, Ming-yong Wang, Yi Zhang. Modulation of active Cr(III) complexes by bath preparation to adjust Cr(III) electrodeposition [J]. 矿物冶金与材料学报(英文版). DOI: 10.1007/s12613-013-0813-5

View details

Jingyin Fei, Guozheng Liang, Wenli Xin, Weikang Wang, Jianghong Liu. Corrosion behaviors of zinc and Zn-Ni alloy compositionally modulated multilayer coatings [J]. 矿物冶金与材料学报(英文版).

View details

Li Peng, Qing Li, Wen Li, Shufen Li. NN adaptive pole placement method for sintering finish point control [J]. 矿物冶金与材料学报(英文版).

View details

Li Peng, Wen Li, Guohuan Lou, Xuyan Tu. An integrated scheme of neural network and optimal predictive control [J]. 矿物冶金与材料学报(英文版).

View details

Int. J. Miner. Metall. Mater., 1999, 6(4): 302-305.

PDF View details

Jianping Wei, Huade Li, Datai Yu. Inverse Modeling of Nonlinear Systems and Application [J]. 矿物冶金与材料学报(英文版).

View details

Jing Wang, Hui Chen. Identification of Nonlinear Dynamic Systems Using Diagonal Recurrent Neural Networks [J]. 矿物冶金与材料学报(英文版).

View details

Int. J. Miner. Metall. Mater., 1994, 1(1-2): 45-52.

PDF View details

邢宏, 邹家祥, 钱桦. 非线性随机振动隔离系统的参数优化 [J]. 矿物冶金与材料学报(英文版).

View details

/

返回文章
返回