- Advanced Photonics
- Vol. 6, Issue 4, 046006 (2024)
Abstract
1 Introduction
The demand for green and sustainable energy-saving technologies has extensively promoted the development of refrigeration technology. Radiative cooling is emerging as a cutting-edge technology, attracting widespread attention for its unique working principle.1
Passive radiative cooling relies on the emission of heat from terrestrial objects to the 3K outer space through the thermal infrared transparency window ( to ), commonly referred to as the “atmospheric window.” To effectively dissipate heat, materials for passive radiative cooling are typically designed with minimal solar absorptance () in the wavelength range of 0.3 to and maximal thermal emittance () within the atmospheric window. Such a property has been realized in various materials such as herbs,1,12 paints,13
In this work, we designed and prepared a temperature-adaptive metasurface radiative cooling device (ATMRD). The ATMRD is a simple sandwich structure consisting of silver (Ag), hafnium oxide (), and patterned . Through optical simulation and experimental verification, the size range of the optimal micro-structure that achieves synergy between solar absorptance and emittance was clarified. Based on the results of simulation optimization, the ATMRD with the best microstructure was fabricated using sample mask-filling engineering. The prepared ATMRD exhibited excellent radiative cooling and solar reflectivity performance, automatically switching the atmospheric window from 0.13 to 0.85 when the surface temperature exceeds a pre-set transition temperature; meanwhile, achieving a satisfactory of 27.71% at 20°C (or 35.44% at 90°C). The potential applications were further evaluated using simulated solar heating tests and infrared (IR) thermography. This scheme represents a significant advancement in balancing conflicting requirements to enhance radiative cooling performance. Simultaneously, the adoption of a simple and cost-effective micro–nano-preparation technology offers a practical solution for the fabrication of micro-nanostructures.
Sign up for Advanced Photonics TOC. Get the latest issue of Advanced Photonics delivered right to you!Sign up now
2 Results
2.1 Design of Radiative Cooling Device and Numerical Modeling
Figure 1(a) illustrates a typical temperature-adaptive radiative device (ATRD) based on the F–P resonator. A highly reflective metal layer is placed at the bottom, followed by a middle layer of highly IR transparent material and a continuous film covers the top. In the insulating (I) state of at (metal–insulator transition temperature), allows high transmission of IR light, and incident IR light is reflected by the bottom metal layer, achieving a low IR absorption. The heat inside the object cannot radiate away, preventing the overcooling phenomenon. By contrast, the ATRD becomes highly absorptant to IR light when switches to the metal (M) state at . The absorption is further amplified in the atmospheric window through a design of a 1/4 wavelength cavity.3 According to Kirchhoff’s law, the IR absorptance () equals emittance () for IR opaque materials.33 Hence, the emittance achieves a temperature-adaptive switch in the atmospheric window from low to high (or high to low) depending on the I–M phase transition of , as illustrated in Fig. 1(b).
Figure 1.(a), (c) Schematics of the ATRD and the ATMRD. (b) Working diagram of the ideal temperature-adaptive radiative cooling. (d) Surface pattern dimensions and cross-sectional diagram of designed ATMRD. (e), (f) Comparison of simulated solar absorptance and thermal emissivity between ATRD and ATMRD. The blue line represents low temperature, and the red line represents high temperature.
Here, the reflective layer uses an Ag metal with the highest reflection compared with other common metals [Figs. S1(a) and S1(b) in the Supplementary Material], which will contribute to achieving a high solar reflection [Fig. S1(c) in the Supplementary Material]. is selected as the dielectric layer due to its high IR transparency at a broad bandwidth.34 The simulated optimal emittance tunability () can be achieved in ATRD by optimizing the thickness of the and layer (Fig. S2 in the Supplementary Material), whereas the is 33.48% at low temperature [the dotted horizontal line in Fig. 2(h)] and 43.51% at high temperature [the dotted horizontal line in Fig. 2(i)]. High solar absorptance, especially at high temperatures, will cause the object’s temperature to increase and weaken the radiation heat dissipation ability to a certain extent during daytime.
Figure 2.(a), (b) Simulated solar absorptance spectra for ATMRD
Aiming at the goal of achieving low and high , we employ a periodic array pattern to replace the continuous film, as depicted in Figs. 1(c) and 1(d), namely, constructing an ATMRD. The incorporation of metasurface can obviously reduce and enhance through rational structural design to tune the resonant frequency to the atmospheric window, as shown in Figs. 1(e) and 1(f). The resonant frequency of the metasurface strongly depends on the structure parameters of the top metallic patterns,35 including geometric shape, arrangement period (), feature size (), and gap () with the relation , as shown in Fig. 1(c). We create simple squares and circular patterns to consider the effects of and on and . For convenience, the samples with different and for squares and circular patterns are labeled as and , respectively, where represents the variation of or from 1 to , and represents the ranging from 0.5 to . The simulated solar absorptance results for ATMRD with a square array structure are presented in Figs. 2(a) and 2(b) and Fig. S3 in the Supplementary Material. The amplitude of the solar absorptance curve for all samples gradually weakens as increases at a fixed or decreases at fixed . The same trend is found in ATMRD with a circular array structure shown in Fig. S4 in the Supplementary Material, indicating that significantly depends on the coverage area of in ATMRD. The integration results in Figs. 2(c) and 2(d) using Eq. (1) (see Sec. 4.2) suggest a proportional relationship between and the coverage area of : the larger the coverage area, the stronger .
Figures 3(a)–3(g) show the simulated emissivity spectra of ATMRD with the square patterns of different and in the frequency range from 2.5 to . Compared with the emissivity spectra of ATRD, additional peaks are induced by the periodic patterns besides the resonant absorption peak. Variations in feature size and gap significantly impact the intensity and position of the emittance peak as is in the metallic phase. When is fixed, both the emissivity intensity and bandwidth increase as becomes narrower. As increases at fixed , the intensity of the enhanced emissivity peak gradually decreases, whereas the long wave absorption edge shifts toward a longer wavelength, as shown in Fig. 3(g). A comparison is carried out for different or in the same period (), as shown in Fig. S5 in the Supplementary Material. It is found that the emissivity intensity and bandwidth obviously increase as increases, or decreases at fixed when is in the metallic phase. The thermal emittances () at different and are calculated with Eq. (2) (presented in Sec. 4.2) by referencing the blackbody radiation spectra and exhibited in Fig. 3(h). At low temperatures, the thermal emittance () of ATMRD is lower than that of ATRD, which is more beneficial in preventing the overcooling of objects at low temperatures. As the feature size increases in ATMRD, the slightly increases due to the increased area of . In comparison, the emittance at high temperature () is significantly enhanced. The of ATMRD increases rapidly with increasing as and reaches the maximum value at for the G0.5 and G1 ATMRDs, followed by a slight decline with a further increase in . As , the of G1 and G0.5 surpasses that of ATRD. However, the of ATMRD with a larger gap () is lower than that of ATRD. The emittance tunability () of ATMRD outperforms ATRD over a broader pattern range when and , with the peak value located at , as shown in Fig. 3(i). The results of ATMRD with a circular array structure (Fig. S6 in the Supplementary Material) also revealed consistent regularity, with slightly lower than square patterns [Fig. 3(i)] because the rounded corners of circular patterns increase the gap . The above results indicate that the appropriate values of feature size and gap (or period) are crucial to achieving emittance enhancement.
Figure 3.(a)–(g) Simulated emissivity spectra at different
2.2 Preparation and Characterization of ATMRD
Based on the above simulated results, we selected three representative feature sizes (, 6, and ) and three gap values (, 2, and ) to fabricate ATMRD using a simple mask-filling engineering. The experimental section describes the fabrication process in detail. Figure 4(a) exhibits a photograph of the ATMRD. It presents a rainbow spectrum under visible light due to the structural color induced by the micro–nanostructure.36 The cross-sectional scanning electron microscopy (SEM) image and energy dispersive spectroscopy (EDS) mapping of the ATMRD in Fig. 4(b) exhibit good adhesion among different layers, providing a fundamental guarantee for the device’s performance. The thickness of the Ag, , and layers is 100, 800, and 50 nm, respectively. Figure 4(c) shows the SEM morphology of the ATMRD’s top view. All samples present well-designed geometries over the range of feature sizes from 4 to and from 1 to , attributed to advanced photolithography and magnetron sputtering deposition technology. This method can significantly reduce manufacturing costs and achieve scale preparation compared with electron beam lithography.
Figure 4.(a) Photograph of ATMRD L4-G1, and the middle colored area is the micro–nanostructure. (b) Cross-sectional SEM image and EDS mappings of ATMRD. (c) SEM morphology of the ATMRD’s top view with different
Figures 5(a) and 5(b) show the X-ray diffractometer (XRD) and Raman results of ATMRD films for phase identification. The XRD pattern in Fig. 5(a) primarily exhibits the crystalline structure of , along with diffraction peaks of Ag. No other impurities are observed. The Raman results in Fig. 5(b) confirm the existence of the M1-phase film.37 In addition, the film reveals an excellent IR switch from high transmittance and high reflectance at , as shown in Figs. 5(c) and 5(d).
Figure 5.(a) XRD pattern of the ATMRD. (b) Raman spectrum at room temperature for
2.3 Radiative Cooling Performance of ATMRD
The emissivity spectra of ATMRD are detected under different temperatures, as depicted in Figs. 6(a)–6(d). The emissivity of the devices changes significantly from low to high values with the increase in temperature. The ATRD exhibits strong absorption for wavelengths larger than , with an emissivity peak at when heated to 90°C [Fig. 6(a)], corresponding to the blackbody radiation peak. For ATMRD at 90°C, the peak intensity of L4-G1 and L6-G1 devices is significantly enhanced and almost reaches 100% [Figs. 6(b) and 6(c)], whereas L8-G1 shows a slight decrease in peak intensity [Fig. 6(d)]. The increasing causes the peak to shift toward longer wavelengths, as indicated in Figs. 6(b)–6(d) and Fig. 7(a). The temperature-dependent thermal emittances are obtained from Figs. 6(a)–6(d) using Eq. (2), intuitively evaluating the , as shown in Figs. 6(e)–6(h). A rapid switch from to is observed when the temperature exceeds . The values for L4-G1, L6-G1, and L8-G1 are 0.85, 0.89, and 0.84, respectively, representing a more than 13% improvement over ATRD’s value of 0.75. As the temperature gradually declines, the can maintain over a broad temperature range and then switch to at due to a typical hysteresis effect of phase transition, contributing to more effective heat dissipation at high temperatures. However, the thermal emissivity gradually decreases as broadens from 1 to at a fixed , as shown in Fig. S7 in the Supplementary Material. Only achieves significantly improved emittance compared with ATRD, as indicated in Fig. 7(b), which is consistent with the simulated results. The of ATMRD-L4-G1 () is enhanced by 20% compared with the ATRD ().
Figure 6.Thermal emissivity spectra at different temperatures, emittance-dependent temperature during heating and cooling, and solar absorptance at low and high temperatures for (a), (e), (i) ATRD, (b), (f), (j) ATMRD L4-G1, (c), (g), (k) ATMRD L6-G1, and (d), (h), (l) ATMRD L8-G1, respectively.
Figure 7.(a) Comparison of the thermal emittance. (b) Comparison of emittance tunability between ATMRD
Figures 6(i)–6(m) and Fig. S8 in the Supplementary Material present the measured solar absorptance spectra at low and high temperatures. The relationship of solar absorptance, feature size, and gap size at 20°C and 90°C is illustrated in Fig. 7(c), respectively. Solar absorptance at low temperatures shows a noticeable decline in ATMRD compared with ATRD (35.25%). When ATMRD switches to the state, solar absorptance rises slightly but remains significantly lower than ATRD’s (46.55%). Solar absorptance of the L4-G1 is 27.71% at low temperatures and rises to 35.44% at high temperatures. ATMRDs with more significant gaps achieve smaller solar absorptance but lower emissivity. Hence, the ATMRD L4-G1 achieves optimal performance with at 20°C ( at 90°C) and (), as shown in Fig. 7(d). Figure 7(e) compares the emittance performance of dynamic radiation cooling devices reported in recent years. The performance of the radiative cooling device designed in this work is better than most literature results, affirming the viability of the design scheme and preparation method presented in this study.
2.4 Application Potential Assessment
We carried out verification experiments to evaluate the practical application potential of ATMRD devices. A solar simulator was employed to conduct a preliminary test of the sunlight reflection capability; the schematic diagram is shown in Fig. 8(a). The 3M tape (), silicon (Si) wafer (), and aluminum (Al) sheet () are selected as reference samples. The results in Fig. 8(b) show that the temperature of ATRD quickly rises and then gradually levels off to 37°C within 120 s, but the temperature is much lower than that of the 3M tape (53°C) and Si wafer (44°C). The ATMRD presents a lower temperature than the ATRD, and the temperature further declines toward that of the Al sheet with decreasing . The temperature of ATMRD L4-G1 is only 33°C, slightly higher than that of the Al sheet (28°C). Such a low solar absorptance significantly reduces solar radiation interference to the cooling device.
Figure 8.(a) Schematic of a simulated solar heating test and (c) IR thermal emittance test, and the
The emittance switching of ATMRD, along with two reference samples, was examined by IR imaging with a thermal IR camera, as shown in Fig. 8(c). One reference sample is an Al sheet, which features a consistently low thermal emittance of 0.10, and the other is the commercial 3M tape, possessing a steady high thermal emittance of 0.96, pasted over the substrate used in ATMRD. Compared with the two references, the ATMRD shows a significant thermal radiation change over the phase transition [Fig. 8(d)]. Before the phase transition, the thermal radiation of ATMRD is consistent with that of the Al sheet and then switches to a high emittance state after the phase transition. Moreover, the emittance exhibits weak dependence on the detection angle, as displayed in Fig. 8(e). The result demonstrates that the ATMRD has excellent solar reflectance and considerable heat dissipation.
2.5 Mechanism Analysis
To elucidate the underlying mechanism responsible for the metasurface-enhanced absorption based on the F–P resonant cavity, electromagnetic field distributions in the cross section of one array cell [Fig. 9(a)] were investigated for ATMRDs with different feature size with in the metal phase. Coupled with the F–P resonances in ATRD [Fig. 9(c)], the multimode polariton resonances are excited in ATMRD, as shown in Figs. 9(d)–9(k) and Figs. S9–S11 in the Supplementary Material, producing the corresponding sharp peaks in the emissivity spectra in Fig. 3. The results show that the polariton resonance strongly depends on the structural parameters of feature size and gap size or period of the metasurface.
Figure 9.(a) Schematic diagram of simulated ATMRD using FDTD, the
The surface plasmon polariton (SPP) resonance emerges on the surface of the units when the feature size (or period ) matches the wavelength of the incident wave. For ATMRD L4-G1, the SPP resonance occurs at due to the strong electric polarization above the units in Fig. 9(d). For L6-G1, L8-G1, and L10-G1, the SPP resonance appears at 7, 9, and , respectively, as shown in Figs. S10C and S11E in the Supplementary Material and Fig. 9(k).
As the incident wavelength , the electric polarization in the medium between the metallic and the Ag substrate gradually increases, as characterized by the polarized electric field vectors, which indicate electric current loops in Figs. 9(e)–9(g). The alternatively polarized electric field induces strong localized magnetic polaritons (MPs).47 The electromagnetic polariton resonances can be understood by employing an equivalent inductor–capacitor (LC) circuit model.48 The sandwiched spacer serves as a capacitor (), whereas the top metallic microstructure and the bottom Ag substrate function as inductors (), as illustrated in Fig. 9(b), forming a resonant LC circuit. The LC circuit induces the magnetic field and develops magnetic dipoles according to Lenz’s law.49 Eventually, the localized strong magnetic field is induced at the resonance frequency of the magnetic and electric dipoles, where the total impedance equals zero.35 In other words, the MP resonance is excited at the resonance frequency, resulting in an enhanced absorption effect. At , an intense MP resonance centers below the gaps of units of L4-G1 in Fig. 9(e). As the wavelength of incident light becomes longer, the induced electric polarization period also changes for a longer period of time. The strong polariton resonances occur at and , with the enhanced magnetic response centering below the units, as shown in Figs. 9(f) and 9(g). Further increasing the incident wavelength, the polarization becomes weak. Hence, the emissivity spectrum experiences gradual suppression at a longer wavelength.
For , the polarized electric field period is smaller than the feature size . Figures 9(i) and 9(j) indicate two localized SPP resonances on one unit for L8-G1 and L10-G1. This is because two periodic polarization electric fields are formed on the surface of each metallic unit when is half of , thus exciting two localized SPP resonances. Meanwhile, the multiple MP resonances can also be excited between layers and Ag layers at resonant frequencies. However, the polarization strength gradually weakens as the becomes much smaller than . Therefore, the metastructure loses its effects, as is much smaller than .
3 Discussion
We have successfully designed and fabricated the temperature-adaptive metasurface radiative device—ATMRD by sample mask-filling engineering. The fabrication of a metasurface resolves the issue of balancing high emittance and low solar absorptance. The metasurface stimulates the localized polariton resonances as in the metallic state, enhancing the absorption of electromagnetic radiation in the thermal IR and declining through reducing the area of . The fabricated ATMRD L4-G1 demonstrates the ability to switch the atmospheric window thermal emittance from 0.13 to 0.85 over the phase transition. Meanwhile, it achieves a commendable of 27.71%. The simulation-based analysis reveals the relationship between the structural parameters of the metasurface and performance, providing a guide for the design of metasurface radiative cooling devices. The solar absorptance is directly proportional to the coverage area of , independent of the pattern shape. However, in terms of emissivity, the structural parameters (feature size and gap size or period) of the metasurface play a pivotal role in triggering the polariton resonances of different manners at specific frequencies. Reduced pattern gaps facilitate the enhancement of IR absorption, whereas the feature sizes of 4 and perform the optimal thermal emittance because they can ultimately excite the polariton resonances in the thermal IR range.
4 Appendix: Materials and Methods
4.1 Fabrication of Radiative Cooling Device
Radiative cooling devices featuring a triple-layer film structure of were fabricated on a substrate using magnetron sputtering. Initially, a reflective layer of Ag was deposited on the substrate using a high-purity Ag target. Subsequently, an 800-nm layer was deposited utilizing a ceramic target. For ATRD, a 50-nm-thick layer of vanadium is directly deposited on the top of the layer using a high-purity metal vanadium target at room temperature. The vanadium film underwent thermal treatment to transform into the M1 phase using a rapid annealing furnace (RTP-500, Beijing East-star Research Office of Applied Physics, Beijing, China). For ATMRD, the microstructure was fabricated through mask-filling engineering, as illustrated in Fig. S12 in the Supplementary Material.
The pre-prepared device underwent plasma cleaning initially to remove surface contaminants. Photoresist (L300, Seoul, Republic of Korea) was spin-coated on the layer at a speed of 4000 r/min, and the patterned potholes on the layer were achieved through exposure in a photolithography machine. Subsequently, pure vanadium was deposited into the patterned potholes, and the photoresist mask, along with excess vanadium, was eliminated through immersion and ultrasonic delamination in an acetone solution. Following this step, the patterned vanadium microstructure underwent additional thermal treatment to transform the (M) phase, thereby producing the final ATMRD.
4.2 Characterization
The morphology was observed by field emission scanning electron microscopy (FE-SEM, Hitachi S-4800, Tokyo, Japan). The phase structure of the films was measured utilizing an XRD (Rigaku SmartLab SE, Tokyo, Japan) and a Raman spectrometer (Horiba LabRAM HR Evolution, Paris, France) equipped with a 532-nm laser. An ultraviolet–visible–near-IR spectrophotometer (UV-VIS-NIR, Agilent UV3600, Santa Clara, California, United States) fitted with an integrating sphere measured the transmittance and reflectivity from 300 to 2500 nm at different temperatures using a homemade heating table. Thermal spectral reflectance was characterized by a Fourier transform IR spectrometer (Nicolet iS50, Madison, Wisconsin, United States) equipped with mid-IR IntegratIR spheres (PIKE Technologies, Inc., Madison, Wisconsin, United States) and mercury cadmium telluride (MCT) detectors over the to wavelength range at 12-deg angle of incidence. According to Kirchhoff’s law of radiation, in a state of thermodynamic equilibrium, the spectral emittance equals the spectral absorptance .33 As a radiative cooling device comprises a thick Ag layer, the transmittance () in the solar and thermal IR wavelength ranges was considered zero. Thus, its thermal spectral emittance in this range was computed as .
The solar absorptance () and thermal emittance () can be calculated from the corresponding spectral data using the following equations:
The solar heating simulation test was carried out using a solar simulator (PLS-SXE300, Beijing Perfectlight Technology Co., Ltd., Beijing, China), and a temperature recorder fitted with a Pt temperature sensor recorded the temperatures. Before measuring each type of material surface, the shutter of the solar simulator remained closed, and the sample was first stabilized at room temperature. Then, the shutter was opened, and the surface temperature was recorded as a function of time. Thermal IR images were captured by an FLIR A310 infrared camera (IR camera) working at a 7.5- to wavelength. To minimize reflections in the surrounding environment, the experiments were conducted in an open, outdoor environment. Before testing, the IR camera was calibrated using a thermocouple and 3M tape (a commercial product with a constant emittance of 0.96) as a reference to ensure that the displayed surface temperature matches the reference sample’s actual temperature. The sample was heated using a heating stage, and the IR camera recorded video images at a fixed detection angle.
4.3 Simulation Design
We used finite-difference time-domain (FDTD, Ansys Lumerical, Canonsburg, Pennsylvania, United States) modeling for numerical simulations of the radiative cooling device. A three-dimensional unit cell model was established, including the submicrometer square or circle, layer, and Ag substrate. The materials’ optical properties were obtained from Ref. 24 for and Ref. 34 for . The optical properties of Ag employed in the simulation were calculated from Palik’s handbook.50 The simulation region had a size of referring to the period of the micrometer disks. The perfect matching layer boundary conditions were applied along the propagation of electromagnetic waves ( direction), and periodic boundary conditions were applied along the direction and direction. A broadband linearly polarized plane wave was directed toward the unit cell from above the structure at normal incidence. Positioned above the plane wave source was a frequency–domain power monitor designed to capture the reflected waves. The simulation exclusively focuses on the transverse magnetic (TM) wave, where metamaterial particles are excited by the component of the magnetic field. Due to the symmetry of the considered patterns, the reflectance remains consistent for both transverse electric and TM waves. All simulations detailed in this paper are conducted within a three-dimensional computational domain, employing a nonuniform structured mesh with a minimum mesh size of 1 nm.35
Junlin Yang is a PhD candidate at the Beijing Institute of Technology, under the supervision of professor Jingbo Li. His research primarily focuses on the structure and performance regulation of thermochromic phase change materials and the design of radiative thermal control devices.
Haibo Jin is a professor at the Beijing Institute of Technology. His main research interests focus on functional materials and energy storage and conversion, with a primary emphasis on key materials and technologies for secondary batteries and hydrogen energy development. He has published over 100 papers, which have been cited more than 10,000 times, and holds 15 authorized invention patents. He was once selected as one of the “top 1% highly cited Chinese authors” by a journal of the Royal Society of Chemistry.
Jingbo Li is a professor at the Beijing Institute of Technology. His research focuses on novel electromagnetic functional materials, multi-field regulation, and device design and fabrication, with a primary emphasis on key materials and technologies for broadband intelligent electromagnetic modulation and the development of energy-saving windows. He has led and completed more than 10 major projects, published over 100 papers, held 12 national invention patents, and edited 2 textbooks.
Biographies of the other authors are not available.
References
[1] T. Li et al. A radiative cooling structural material. Science, 364, 760-763(2019).
[7] X. Zhao et al. A solution-processed radiative cooling glass. Science, 382, 684-691(2023).
[19] D. Zhao, H. Tang. Staying stably cool in the sunlight. Science, 382, 644-645(2023).
[26] S. Abdollahramezani et al. Dynamic hybrid metasurfaces. Nano Lett., 21, 1238-1245(2021).
[33] J. Agassi. The Kirchhoff-Planck radiation law. Science, 156, 30-37(1967).
[50] E. D. Palik. Handbook of Optical Constants of Solids, 3(1998).
Set citation alerts for the article
Please enter your email address