Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Tunable mid-infrared selective emitter based on inverse design metasurface for infrared stealth with thermal management

Open Access Open Access

Abstract

Infrared (IR) stealth with thermal management is highly desirable in military applications and astronomy. However, developing selective IR emitters with properties suitable for IR stealth and thermal management is challenging. In this study, we present the theoretical framework for a selective emitter based on an inverse-designed metasurface for IR stealth with thermal management. The emitter comprises an inverse-designed gold grating, a Ge2Sb2Te5 (GST) dielectric layer, and a gold reflective layer. The hat-like function, which describes an ideal thermal selective emitter, is involved in the inverse design algorithm. The emitter exhibits high performance in IR stealth with thermal management, with the low emissivity (ɛ3–5 µm =0.17; ɛ8–14 µm =0.16) for dual-band atmospheric transmission windows and high emissivity (ɛ5–8 µm =0.85) for non-atmospheric windows. Moreover, the proposed selective emitter can realize tunable control of thermal radiation in the wavelength range of 3–14 µm by changing the crystallization fraction of GST. In addition, the polarization-insensitive structure supports strong selective emission at large angles (60°). Thus, the selective emitter has potential for IR stealth, thermal imaging, and mid-infrared multifunctional equipment.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Infrared (IR) stealth technology aims to reduce the IR signature of objects against IR detectors. The application scenarios of IR stealth include mid-infrared (MIR) photodetectors, thermal imaging systems, and heat-seeking missiles [15]. According to the Stefan-Boltzmann law, the radiated thermal energy per unit of a hot surface is related to the temperature of the surface and the infrared emissivity from the surface. Therefore, IR stealth with thermal management can effectively satisfy the requirements of reducing the surface temperature and realizing low IR emissivity. This strategy involves achieving low emissivity in dual-band atmospheric transmission windows (ATWs), which include the mid-wavelength infrared (MWIR; wavelength (λ) ≈ 3–5 µm) and long-wavelength infrared (LWIR; wavelength (λ) ≈ 8–14 µm), IR stealth, and high emissivity in the wavelength ranges of 5–8 µm for radiative cooling. In this case, the ideal thermal selective emitter is a “top-hat-like” function, which shows a unit value of emittance from 5 to 8 µm and zero in dual-band MIR ATWs [6].

Over the past decade, several thermal photonic structures with wavelength-selective abilities have attracted increasing attention. The proposed photonic structures are employed to realize several applications, including solar cells, radiative cooling, and IR stealth [712]. These thermal photonic structures are expected to find a wide range of applications, such as in industry, environment observation, and military. Particularly, in military applications, IR stealth plays an important role in avoiding IR signal tracking and identification. Pan et al. [13] designed and fabricated different photonic grating structures to implement MIR-compatible stealth. Zhu et al. [14] proposed a wavelength-selective emitter based on a multilayer metal structure for applications in high-temperature infrared stealth with efficient thermal management. Song et al. [15] used a metal (insulator) metasurface structure to realize multifunctional thermal stealth. Due to the limitations of traditional structural design methods, it is challenging to implement the top-hat-like function in selective emitters.

The innovative field of inverse design has recently been transforming conventional nanophotonics and discovering digital nanophotonics rather than engineering, primarily by manually selecting parameters or employing semi-analytical models [16,17]. Therefore, the design process of inverse design is different from that of traditional design methods [18]. In some problems, spectral characteristics are directly related to the application of ideal spectral curves, including but not limited to sunlight absorption, radiative cooling, nanophotonics devices, and multi-spectral camouflage problems [1923]. Algorithms that improve the performance of nanophotonics and consider the material distribution for processing technology have emerged, such as direct binary search (DBS) algorithms [24,25], genetic algorithms [26,27], topology optimization [2830], and neural network [3133]. An inverse design that considers the “top-hat-like” function of spectral characteristics as the objective function is expected to achieve IR stealth with thermal management. In addition, tunable thermal emitters with spectral selectivity are required for multifunctional applications such as IR stealth [34] and thermal imaging [35].

Phase-change materials (PCMs) have attracted significant attention, owing to their tunability without changing the structural parameters. Ge2Sb2Te5 (GST) is a PCM widely used in optical storage media [36,37], tunable photonics computing [38], and infrared thermal emission [39]. In the infrared band, GST has a dielectric permittivity of ∼16.0 and ∼34.0 in its amorphous and crystalline states, respectively. Some studies have illustrated that the crystalline state of GST can be converted using methods such as heating, laser pulses, and applied voltage [4042]. Once crystallization or reamorphization is accomplished, it can be maintained for many years at room temperature [43]. In addition, crystalline GST (c-GST) exhibits a high absorption in the MIR region, which can generate resonance at different wavelengths by changing the thickness of GST, whereas amorphous GST (a-GST) is transparent [44]. Based on these merits, GST is an ideal candidate for tuneable MIR multifunctional devices.

In this study, we demonstrate a PCM-based metasurface selective emitter (MSE) that combines the inverse design algorithm with IR stealth for the first time. The DBS algorithm was employed to topologically optimize the top single-layer nanostructured metallic pattern of the MIM metasurface. A GST dielectric layer was introduced to maximize the emittance of the 5–8 µm spectrum and implement the switching function. Similar to the conventional MIM structure, a layer of optically thick gold mirrors was used as the bottom layer. The simulation results reveal a low emissivity (ɛ3–5 µm =0.17; ɛ8–14 µm =0.16) in the dual-band MIR ATWs and a high emissivity (ɛ5–8 µm =0.85) in the non-window band, achieving a performance close to that of the ideal thermal selective emitter. According to the principle of thermal imaging, the proposed MSE can be disguised as an object at a temperature of 70 °C and a real temperature of 200 °C. Compared to other reported thermal stealth structures, this structure has good compatibility between IR stealth and thermal management. In addition, based on the phase-changing characteristics of GST, the proposed metasurface can realize continuous conversion of IR stealth to thermal imaging. The angular and polarization insensitivities of the metasurface were verified. Based on the above advantages, MSE has great potential for applications in IR stealth, thermal imaging, and multifunctional infrared detection.

2. Structure and model

2.1 Background and MIM structure for IR stealth

As shown in Fig. 1(a), a realistic atmospheric transmittance model plotted for Ref. [45] is exhibited as a light-blue shaded area at 3–14 µm. Owing to the low loss transmission in the wavelength ranges of 3–5 µm and 8–14 µm, dual-band MIR ATWs are used in IR detectors and IR imaging [4648]. For the wavelength range of the non-atmospheric window, it is difficult to implement long distances. Therefore, a non-atmospheric window was adopted to dissipate heat without influencing IR stealth [49]. According to the theory of blackbody radiation, it is approximately 876 W/m2 at a temperature of 473 K in the wavelength range of 5–8 µm, as shown in Fig. 1(a) by the red solid curve. According to Kirchhoff's principle, for objects in a steady state, the absorption spectrum is equivalent to the emission spectrum. Consequently, an ideal “hat-like” function with unit absorption in the non-atmospheric window band and zero absorption in the dual-band MIR ATWs was considered for realizing IR stealth with thermal management. Owing to their superiority in processing and absorption performance, MIM structures are widely used in the design of perfect absorbers. Patterned metal gratings have been developed to excite different surface plasmon (SP) resonances and realize strong spectral absorption, such as strips, crosses, and rings [5052]. Moreover, the MIM structure with PCMs can realize functional conversion and modulation of the working wavelength without changing the structural parameters [53,54].

 figure: Fig. 1.

Fig. 1. (a) Atmospheric transmittance (light-blue shaded area) with detected band (3−5 µm and 8−14 µm) and undetected band (5−8 µm) indicated and radiative heat flux at 473 K. (b) Top view of one unit cell of the design has been divided into 50 × 50 pixels. (The yellow pixel represents gold, and the purple pixel represents air, and the GST is exposed) (c) Schematic of the PCM-based MSE. The MSE is developed with inverse design gold grating coated on the GST and gold mirror. The periodic structural parameters are P = 1 µm. The thicknesses of inverse design gold grating (t1), GST dielectric layer (t2) and gold mirror (t3) are 40, 230 and 300 nm, respectively.

Download Full Size | PDF

2.2 Algorithm optimization and simulation

The top periodic unit is divided into 50 × 50 pixels, as shown in Fig. 1(b), and each is a square of 20 nm × 20 nm, which can be fabricated by electron beam lithography followed by the Au lift-off process. The pixels have a binary state of dielectric properties, i.e., Au or air. To obtain a polarization-insensitive metasurface, we improved the inverse design to obtain an 8-fold symmetrical unit cell. The option of initial structure is discussed in the Fig. S1. And we verify that all-air structure is a appropriate solution to research the locally optimal structure. Therefore, the algorithm optimization process starts from an “all-air” structure comprising a bare thin-film GST and gold mirror. Figure 1(c) shows the schematic of the PCM-based MSE that comprises an inverse-designed gold grating, a GST dielectric layer and a gold reflective layer.

As shown in Fig. 2, we describe all steps for designing the metasurface using the DBS algorithm. In the dashed frame in Fig. 2, we show the top views of the first switching pattern, which chooses the first pixel as the starting point, and switches the pixel state of the first pixel and the final optimized top gold pattern of the proposed metasurface. We switched one pixel state (Au or air) and compared its performance in terms of the objective function with that of the unchanged structure to determine the state of the pixel. In an iteration, each pixel is scanned by a column to determine the state distribution of different pixels. This algorithm-based optimization method provides enhanced flexibility for generating complex and nonintuitive nanostructure features, which enables the designed metasurface to be closer to the objective function and surpass traditional resonant structural elements. Subsequently, a figure-of-merit (FOM) that can approximate the hat function of the ideal emitter is defined as follows:

$$FOM = {\alpha _1} \times \frac{{\int_{{\lambda _1}}^{{\lambda _2}} {\varepsilon (\lambda )d\lambda } }}{{{\lambda _2} - {\lambda _1}}} + {\alpha _2} \times \left( {1 - \frac{{\int_{{\lambda_{\min }}}^{{\lambda_1}} {\varepsilon (\lambda )d\lambda } }}{{{\lambda_1} - {\lambda_{\min }}}}} \right) + {\alpha _3} \times \left( {1 - \frac{{\int_{{\lambda_2}}^{{\lambda_{\max }}} {\varepsilon (\lambda )d\lambda } }}{{{\lambda_{\max }} - {\lambda_2}}}} \right),$$
where λ1, λ2, λmin, and λmax were 5, 8, 3, and 14 µm, respectively. ɛ(λ) represents the emissivity corresponding to different wavelengths. The constants α1, α2, and α3 are the weights of the different wavelength ranges of the selective emitter. To obtain the normalized results, they were 0.5, 0.25, and 0.25 in the simulation. The different weights of normalized α1, α2, and α3 are shown in Fig. S2. We utilized the DBS algorithm to maximize the FOM and determine the optimal topological gold grating pattern. To meet the processing requirements, we also discuss the larger pixels in Fig. S3.

 figure: Fig. 2.

Fig. 2. Designed schematic of the DBS algorithm, initial structure, and optimal metasurface. The steps of the DBS algorithm design include choosing pixel, switching the state of pixel, calculating FOM, judging the spectral characteristic according to the FOM, and outputting the optimal topological structure. The initial structure and the optimal metasurface are given in dashed frames. The initial point position and scanning order are indicated in the dashed frame initial structure diagram.

Download Full Size | PDF

Compared with some other optimization methods, DBS method has the characteristics of orientation, simplicity, and scalability. We also make a comparison of computational cost between the DBS algorithm and genetic algorithm (GA). In the Fig. S4, the comparison illustrate that the DBS method has the advantages of strong orientation and low computational cost. Due to the strict orientation, this makes it challenging to explore the global optimum. There are also some studies to strengthen the DBS method by improving on the starting structure and the scan path.

The finite-difference time domain (FDTD) is a method of approximately solving Maxwell's equations, and has been used for numerical simulations. The commercial software Lumerical FDTD Solutions was employed in a simulation where IR plane waves propagated to the metasurface along the z-direction. Both the boundary conditions in the x- and y-directions are periodic layers. The upper and lower boundary conditions of the z-direction are perfectly matched layers and metal layers, respectively.

The spectral transmittance T(λ) and reflectance R(λ) of the proposed metasurface were monitored during the simulation. The absorptivity A(λ) was calculated as follows: A(λ) 1 − T(λ) − R(λ). According to Kirchhoff’s law, absorptivity A(λ) equals emissivity ɛ(λ). Consequently, the emissivity spectrum ɛ(λ) of the proposed metasurface was obtained from the simulation results. Combining the FDTD methods and DBS algorithm, it takes about 70 hours (6 iteration) to get the optimized metasurface on a computer with a 6-core central processing unit (Intel Core i7-9750H). The total memory of the computer is 32 GB (SAMSUNG DDR4 2666 MHz).

2.3 Material optical characterization

As shown in Fig. 3, the dielectric constants of the GST in the crystalline and amorphous states are plotted for Ref. [55,56]. In the IR band, amorphous GST is a low-loss material, whereas crystalline GST is a lossy material, which allows it to achieve intrinsic absorption in the MIR spectrum with changes in thickness [39]. In addition, GST materials have several intermediate phases that can be controlled by various methods such as heating time, laser irradiation, and gating voltage. The effective dielectric constants of GST with different degrees of crystallinity were calculated using the Lorentz–Lorenz relation [5759]:

$$\frac{{{\varepsilon _{GST}}({\lambda ,C} )- 1}}{{{\varepsilon _{GST}}({\lambda ,C} )+ 2}} = C \times \frac{{{\varepsilon _{cGST}}({\lambda ,C} )- 1}}{{{\varepsilon _{cGST}}({\lambda ,C} )+ 2}} + ({1 - C} )\times \frac{{{\varepsilon _{aGST}}({\lambda ,C} )- 1}}{{{\varepsilon _{aGST}}({\lambda ,C} )+ 2}},$$
where ɛaGST and ɛcGST are the dielectric constants of a-GST and c-GST, respectively. The constant C represents the doping coefficient of c-GST, which ranges from 0 to 1. The relative permittivity of Au was obtained from Palik’s handbook [60].

 figure: Fig. 3.

Fig. 3. Dielectric permittivity (ɛGST) of crystalline and amorphous GST.

Download Full Size | PDF

3. Results and discussion

3.1 Simulation result and physical mechanism

First, we consider the case in which IR plane waves are normally illuminated on an MSE, as illustrated in Fig. 1(c). Figure 4(a) shows the emissivity of the proposed PCM-based MSE in the IR band. Owing to the intrinsic absorption of GST and the coupling mode of the patterned gold grating, the absorption of the proposed structure was 85.3% in the wavelength range of 5–8 µm. Low emissivity is achieved in the dual-band atmospheric windows (ɛ3–5 µm =0.17; ɛ8–14 µm =0.16). The IR stealth with thermal management of the proposed inverse-designed metasurface was fully verified by comparison with a dual-layer stacked structure (GST thin layer and gold optical layer). The blue line in Fig. 4(a) shows the emission spectrum of the stacked structure without the patterned gold grating. The ideal hat-like function, which is compatible with IR stealth with thermal management, is shown in Fig. 4(a) by the red dashed line. It is evident that the emissivity spectrum of the metasurface is closer to the ideal hat-like function than that of the bare GST thin film and gold mirror.

 figure: Fig. 4.

Fig. 4. (a) Simulated spectral emissivity at normal incidence and spectral emissivity of ideal thermal selective emitter. (b) Electric field distribution in the indicated xy- and xz- cross sections at the two peak wavelengths.

Download Full Size | PDF

To explain the physical mechanism of the high average emissivity in the wavelength range of 5–8 µm, we provide the electric field distribution of the two resonance peaks of the proposed metasurface, as shown in Fig. 4(b). The emission peak at 5.5 µm wavelength is mainly due to the SP resonance between the patterned gold gratings. The other emission peak at 7.0 µm wavelength is realized by the coupling between the SP resonance and Fabry-Pérot cavity of the c-GST and Au thin film [39]. It is worth noting that the high emissivity of wavelengths ranging from 5 to 8 µm has no effect on thermal imaging because of high atmospheric absorption. Therefore, this part of the energy was quickly dissipated through the high emissivity of the metasurface. For a 473 K object, the radiation power of this part of the spectrum was as high as 876 W/m2.

3.2 Performance of IR stealth with thermal management

Most natural objects have a blackbody-like radiation spectrum. For the IR detector, the temperature of the object is determined by the detected radiation power [13,14]:

$${T_r} = {P^{ - 1}}({\varepsilon _{IR}},T).$$
where Tr represents the temperature obtained using the inverse function of the blackbody radiation spectrum, P represents the radiation power detected by the IR detector, which corresponds to the radiation power radiated by the black body at temperature T. ɛIR is the emissivity of the IR detector in the working wavelength (7.5–14 µm), and ɛIR = 1 for the infrared camera. The detection power of an infrared camera is generally determined by two factors: the black-body radiation power and the environmental reflection power of the object [1315] as
$$P(\varepsilon ,T) = {P_{rad}}(\varepsilon ,T) + {P_{ref}}(\varepsilon ,{\varepsilon _a},{T_a}) = \varepsilon (\lambda ){I_{BB}}(T) + [{1 - \varepsilon (\lambda )} ]{\varepsilon _a}(\lambda ){I_{BB}}({T_a}),$$
where ɛ and ɛa represent the emissivity of the designed structure and ambient environment, respectively. T and Ta represent the temperature of the object and ambient environment, respectively. IBB represents the blackbody irradiance at the corresponding temperature, and the blackbody radiation threshold is determined by the working wavelength of the infrared camera. This indicates that when the emissivity of the object at the working wavelength of the infrared camera is lower, it will make objects with the same temperature as the environment look “colder” than the environment. It is easy to detect a hot object by an infrared camera in the low ambient temperature. To realize infrared stealth, the low emissivity in the working wavelength of the infrared camera makes the hot object conceal in the low ambient temperature.

According to the emissivity of the MSE, we plotted the detection power under an infrared camera working in the wavelength range of 7.5–14 µm. Figure 5(a) shows that when our structure is at 200 °C, the observation power under the infrared camera is 406 W/m2, which corresponds to the black body radiation power at 70 °C. Therefore, our metasurface is expected to be disguised as an object at a temperature of 70–200 °C. The measured radiation temperatures (Tr) of the objects are significantly lower than their real temperatures (T), and the radiation temperatures increase with increasing real temperatures. In addition, we provide the radiative heat flux curve of the MSE at 473 K. Compared with an ideal blackbody at 473 K, the proposed metasurface has obvious wavelength-selective emission, which has a high heat dissipation in the wavelength range of 5–8 µm, as shown in Fig. 5(b).

 figure: Fig. 5.

Fig. 5. (a) The relation between radiation temperatures and real temperatures of the MSE: the lines are calculated with the integrated power. (b) The spectral radiation intensity of the blackbody and MSE at 473 K, the black and red lines represent the blackbody and the MSE, respectively.

Download Full Size | PDF

To evaluate the performance of our designed metasurface structure, we defined the efficiency of the thermal stealth compatible with thermal management. The performance of the MSE should be determined in two parts: IR stealth, including dual-band ATWs, and thermal management. Therefore, we define the efficiency as the product of the dual-band MIR efficiency and the out-of-band efficiency.

$$eff = eff_{MWIR}^{in} \cdot ef{f^{out}} \cdot eff_{LWIR}^{in}.$$
$$eff_{MWIR}^{in} = 1 - \frac{{\int_{{\lambda _{\min }}}^{{\lambda _1}} {\varepsilon (\lambda )d\lambda } }}{{{\lambda _1} - {\lambda _{\min }}}},ef{f^{out}} = \frac{{\int_{{\lambda _1}}^{{\lambda _2}} {\varepsilon (\lambda )d\lambda } }}{{{\lambda _2} - {\lambda _1}}},$$
$$eff_{LWIR}^{in} = 1 - \frac{{\int_{{\lambda _2}}^{{\lambda _{\max }}} {\varepsilon (\lambda )d\lambda } }}{{{\lambda _{\max }} - {\lambda _2}}}.$$

Here, $eff_{MWIR}^{in}$, $ef{f^{out}}$, and $eff_{LWIR}^{in}$ represent dual-band MIR and out-of-band efficiencies, respectively, and ɛ(λ) is the spectral emissivity of the optimized emitter. The emissivity in the dual-band MIR ATWs reflects the stealth performance of the proposed structure under infrared imaging. Out-of-band emissivity determines the thermal management performance of the proposed selective emitter, which is related to radiative cooling.

We also provide a comparison of these results with other works [13,14,49,61], as shown in Table 1. It can be seen that the designed structural performance has certain advantages, and it has good compatibility with the IR stealth and thermal management performance. For the FOM evaluation, we used equal weights for normalized α1, α2, and α3(α1= 0.5, α2= α3= 0.25). The definition of compatible efficiency (eff), which is another evaluation standard for the compatibility between IR stealth and thermal management inspired by a previous study [19], is given by Eqs. (5) and (6), respectively. We also compare the designed metasurface with the initial structure (dual-layer film). The comparison shows that the optimized metasurface significantly improves the performance of infrared stealth with thermal management.

Tables Icon

Table 1. Performance of the selective emitters for IR stealth with thermal management under two types of evaluations

3.3 Switchable thermal management and thermal imaging

In the previous study, a dual-layer PCM-based structure has been fabricated and realized the tunable emissivity in the LWIR [62]. Notably, the PCM-based structure is expected to achieve functional transformation without changing the structural parameters [63]. Inspired by these, we discuss the tunable properties of the designed metasurface in the MWIR. This will enable the designed selective emitter to achieve a light-dark-like transition that can transform IR stealth into a thermal imaging configuration that has the potential to be used in the aerospace and military fields. Figure 6 shows the emissivity spectral characteristics under different crystalline fraction parameters. According to the Eq. (2), we obtained the dielectric permittivity of GST under different crystallization fractions. With an increase in the crystalline fraction, the selective emission spectrum exhibited a red shift, which was related to the resonance of the MIM structure and the properties of the GST material. We also provide a color map of the metasurface with MWIR ATWs (the integral transmittance is 67% in the wavelength range of 3–5 µm) as the background. Owing to the IR stealth performance of the MSE, it is difficult to distinguish the environment from the proposed MSE with c-GST. As the crystalline fraction decreases, the proposed PCM-based MSE applied to the IR stealth transforms into a thermal imaging structure. Furthermore, the PCM-based MSE can achieve tunable absorption for the range of approximately 2 µm in the wavelength range of 3–8 µm, including selective emission of 3–5 µm, which has potential applications in the field of thermal imaging [38], encryption with IR color [64], and adaptive thermal management [5,65]. In addition, because the phase change material GST can achieve rapid switching within the nanosecond/picosecond range [66,67], the proposed MSE can achieve rapid functional switching.

 figure: Fig. 6.

Fig. 6. Simulated spectral emissivity of the plasmonic thermal emitter gradually changes as crystallization fraction of GST increases from 0% to 100%. The colormaps of PCM-based MSE with different crystallization fractions comparing the MWIR ATWs (the integral transmittance is 67% in the wavelength range of 3–5 µm) have been marked in each emissivity spectrum.

Download Full Size | PDF

3.4 Dependence of the IR stealth performance with polarization and incident angle

It is necessary for the emitters to support a wide range of angles so that light waves entering from different directions can be well emitted by the nanostructure [6870]. In order to obtain a polarization-insensitive MSE configuration, an 8-fold symmetrical unit cell was designed using the improved DBS algorithm. The emissivity values for different polarization angles are shown in Fig. 7(a). Because all polarized light waves can be decomposed into the sum of the x-direction and y-direction light waves that respond to a symmetrical structure, the proposed metasurface exhibits perfect polarization insensitivity. In addition, we discuss the incident-angle insensitivity of the designed metasurface. As shown in Fig. 7(b), when large-angle light (60°) is incident, the designed metasurface can still maintain good performance.

 figure: Fig. 7.

Fig. 7. (a) Emissivity of the MSE in the range of 3–14 µm, and (b) MSE is insensitive to the angle and polarization.

Download Full Size | PDF

4. Conclusion

In conclusion, we demonstrated a PCM-based MSE for IR stealth using thermal management. The MSE is based on an MIM structure, which includes an inverse-designed gold grating, GST dielectric layer, and gold mirror. Simulation results show that low emissivity in the dual-band MIR ATWs and high emissivity in the non-window wavelength range of 5–8 µm were achieved. By combining inverse design and ideal thermal selective emitter theory, the proposed selective emitter exhibits better performance than traditional IR stealth with thermal management. Owing to the introduction of PCM materials, the proposed metasurface can realize tunable control of thermal radiation in the wavelength range of 3–14 µm. In addition, the proposed selective emitter has good insensitivity to polarization and supports large-angle incidence (60°), which is crucial to IR stealth equipment. Our study provides an approach and a reference for inverse-design IR stealth with thermal management.

Funding

National Natural Science Foundation of China (60907003, 61805278); Foundation of NUDT (JC13-02-13, ZK17-03-01); Natural Science Foundation of Hunan Province (13JJ3001); Program for New Century Excellent Talents in University (NCET-12-0142).

Acknowledgements

The authors gratefully thank Yifei Xiao for her help with the schematics of the configurations, thank Jie Huang and Zheng Peng for their useful advice.

Disclosures

The authors declare that there are no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Supplemental document

See Supplement 1 for supporting content.

References

1. N. Baranwal and S. P. Mahulikar, “Infrared Signature of Aircraft Engine with Choked Converging Nozzle,” J. Thermophys. Heat Transf. 30(4), 854–862 (2016). [CrossRef]  

2. L. Phan, W. G. Walkup, D. D. Ordinario, E. Karshalev, J. Jocson, A. M. Burke, and A. A. Gorodetsky, “Reconfigurable Infrared Camouflage Coatings from a Cephalopod Protein,” Adv. Mater. 25(39), 5621–5625 (2013). [CrossRef]  

3. X. Xie, X. Li, M. Pu, X. Ma, K. Liu, Y. Guo, and X. Luo, “Plasmonic Metasurfaces for Simultaneous Thermal Infrared Invisibility and Holographic Illusion,” Adv. Funct. Mater. 28(14), 1706673 (2018). [CrossRef]  

4. C. Y. Xu, G. T. Stiubianu, and A. A. Gorodetsky, “Adaptive infrared-reflecting systems inspired by cephalopods,” Science 359(6383), 1495–1500 (2018). [CrossRef]  

5. O. Salihoglu, H. Uzlu, O. Yakar, S. Aas, O. Balci, N. Balci, S. Olcum, S. Süzer, and C. Kocabas, “Graphene Based Adaptive Thermal Camouflage,” Nano Lett. 18(7), 4541–4548 (2018). [CrossRef]  

6. N. Lee, T. Kim, J. Lim, I. Chang, and H. H. Cho, “Metamaterial-Selective Emitter for Maximizing Infrared Camouflage Performance with Energy Dissipation,” ACS Appl. Mater. Interfaces 11(23), 21250–21257 (2019). [CrossRef]  

7. S. Wu, Y. Ye, Z. Jiang, and T. Yang, “Large-Area, Ultrathin Metasurface Exhibiting Strong Unpolarized Ultrabroadband Absorption,” Adv. Opt. Mater. 7(24), 1901162 (2019). [CrossRef]  

8. X. Jiang, Z. Zhang, D. Chen, K. Wen, and J. Yang, “Tunable multilayer-graphene-based broadband metamaterial selective absorber,” Appl. Opt. 59(35), 11137 (2020). [CrossRef]  

9. Z. Chen, Y. Weng, J. Liu, N. Guo, Y. Yu, and L. Xiao, “Dual-band perfect absorber for a mid-infrared photodetector based on a dielectric metal metasurface,” Photonics Res. 9(1), 27 (2021). [CrossRef]  

10. A. P. Raman, M. A. Anoma, L. Zhu, E. Rephaeli, and S. Fan, “Passive radiative cooling below ambient air temperature under direct sunlight,” Nature 515(7528), 540–544 (2014). [CrossRef]  

11. L. Zhao, H. Liu, Z. He, and S. Dong, “All-metal frequency-selective absorber/emitter for laser stealth and infrared stealth,” Appl. Opt. 57(8), 1757–1764 (2018). [CrossRef]  

12. J. Zhang, J. Shi, D. Zhao, Q. Wang, and C. Wang, “Realization of compatible stealth material for infrared, laser and radar based on one-dimensional doping-structure photonic crystals,” Infrared Phys. Technol. 85, 62–65 (2017). [CrossRef]  

13. M. Pan, Y. Huang, Q. Li, H. Luo, H. Zhu, S. Kaur, and M. Qiu, “Multi-band middle-infrared-compatible camouflage with thermal management via simple photonic structures,” Nano Energy 69, 104449 (2020). [CrossRef]  

14. H. Zhu, Q. Li, C. Zheng, Y. Hong, Z. Xu, H. Wang, W. Shen, S. Kaur, P. Ghosh, and M. Qiu, “High-temperature infrared camouflage with efficient thermal management,” Light: Sci. Appl. 9(1), 60 (2020). [CrossRef]  

15. J. Song, S. Huang, Y. Ma, Q. Cheng, R. Hu, and X. Luo, “Hierarchical metamaterials for laser-infrared-microwave compatible camouflage,” Opt. Express 28, 675 (2020). [CrossRef]  

16. J. Huang, H. Ma, D. Chen, H. Yuan, J. Zhang, Z. Li, J. Han, J. Wu, and J. Yang, “Digital nanophotonics: the highway to the integration of subwavelength-scale photonics,” Nano Photonics 10(3), 1011–1030 (2021). [CrossRef]  

17. L. Ma, J. Li, Z. Liu, Y. Zhang, N. Zhang, S. Zheng, and C. Lu, “Intelligent algorithms: new avenues for designing nanophotonic devices,” Chin. Opt. Lett. 19(1), 011301 (2021). [CrossRef]  

18. S. Molesky, Z. Lin, A. Y. Piggott, W. Jin, J. Vucković, and A. W. Rodriguez, “Inverse design in nanophotonics,” Nat. Photonics 12(11), 659–670 (2018). [CrossRef]  

19. Z. A. Kudyshev, A. V. Kildishev, V. M. Shalaev, and A. Boltasseva, “Machine-learning-assisted metasurface design for high-efficiency thermal emitter optimization,” Appl. Phys. Rev. 7(2), 021407 (2020). [CrossRef]  

20. Y. Shi, W. Li, A. Raman, and S. Fan, “Optimization of Multilayer Optical Films with a Memetic Algorithm and Mixed Integer Programming,” ACS Photonics 5(3), 684–691 (2018). [CrossRef]  

21. D. Chae, H. Lim, S. So, S. Son, S. Ju, W. Kim, J. Rho, and H. Lee, “Spectrally Selective Nanoparticle Mixture Coating for Passive Daytime Radiative Cooling,” ACS Appl. Mater. Interfaces 13(18), 21119–21126 (2021). [CrossRef]  

22. H. Ma, J. Huang, K. Zhang, and J. Yang, “"Inverse-designed arbitrary-input and ultra-compact 1×N power splitters based on high symmetric structure,” Sci. Rep. 10(1), 11757 (2020). [CrossRef]  

23. H. Zhu, Q. Li, C. Tao, Y. Hong, Z. Xu, W. Shen, S. Kaur, P. Ghosh, and M. Qiu, “Multispectral camouflage for infrared, visible, lasers and microwave with radiative cooling,” Nat. Commun. 12(1), 1805 (2021). [CrossRef]  

24. Y. Liu, K. Xu, S. Wang, W. Shen, H. Xie, Y. Wang, S. Xiao, Y. Yao, J. Du, Z. He, and Q. Song, “Arbitrarily routed mode-division multiplexed photonic circuits for dense integration,” Nat. Commun. 10(1), 3263 (2019). [CrossRef]  

25. H. Ma, J. Huang, K. Zhang, and J. Yang, “Ultra-compact and efficient 1×2 mode converters based on rotatable direct-binary-search algorithm,” Opt. Express 28(11), 170109 (2020). [CrossRef]  

26. X. Jiang, H. Yuan, D. Chen, Z. Zhang, T. Du, H. Ma, and J. Yang, “Metasurface Based on Inverse Design for Maximizing Solar Spectral Absorption,” Adv. Opt. Mater. 9(19), 2100575 (2021). [CrossRef]  

27. Z. Yu, H. Cui, and X. Sun, “Genetically optimized on-chip wideband ultracompact reflectors and Fabry–Perot cavities,” Photonics Res. 5(6), B15–B19 (2017). [CrossRef]  

28. H. Zhang, L. Qian, C. Wang, C. Ji, Y. Liu, J. Chen, and C. Lu, “Topological rainbow based on graded topological photonic crystals,” Opt. Lett. 46(6), 1237 (2021). [CrossRef]  

29. J. Huang, J. Yang, D. Chen, W. Bai, J. Han, Z. Zhang, J. Zhang, X. He, Y. Han, and L. Liang, “Implementation of on-chip multi-channel focusing wavelength demultiplexer with regularized digital metamaterials,” Nano Photonics 9(1), 159–166 (2019). [CrossRef]  

30. J. Huang, J. Yang, D. Chen, X. He, Y. Han, J. Zhang, and Z. Zhang, “Ultra-compact broadband polarization beam splitter with strong expansibility,” Photon. Res. 6(6), 574 (2018). [CrossRef]  

31. S. So, T. Badloe, J. Noh, J. B. Abad, and J. Rho, “Deep learning enabled inverse design in nanophotonics,” Nanophotonics 9(5), 1041–1057 (2020). [CrossRef]  

32. S. So and J. Rho, “Designing nanophotonic structures using conditional deep convolutional generative adversarial networks,” Nanophotonics 8(7), 1255–1261 (2019). [CrossRef]  

33. S. So, D. Lee, T. Badloe, and J. Rho, “Inverse design of ultra-narrowband selective thermal emitters designed by artificial neural networks,” Opt. Mater. Express 11(7), 1863 (2021). [CrossRef]  

34. H. Ji, D. Liu, C. Zhang, and H. Cheng, “VO2ZnS core-shell nanoparticle for the adaptive infrared camouflage application with modified color and enhanced oxidation resistance,” Sol. Energ. Mat. Sol. C. 176, 1–8 (2018). [CrossRef]  

35. S. Park, Y. Kim, Y. Moon, J. Cho, and S. Kim, “Tuning of polarized room-temperature thermal radiation based on nanogap plasmon resonance,” Opt. Express 28(10), 15472 (2020). [CrossRef]  

36. M. Wuttig and N. Yamada, “Phase-change materials for rewriteable data storage,” Nat. Mater. 6(11), 824–832 (2007). [CrossRef]  

37. Y. Hu, H. Zou, J. Zhang, J. Xue, and Y. Sui, “Ge2Sb2Te5/Sb superlattice-like thin film for high speed phase change memory application,” Appl. Phys. Lett. 107(26), 263105 (2015). [CrossRef]  

38. J. Feldmann, N. Youngblood, C. D. Wright, H. Bhaskaran, and W. H. P. Pernice, “All-optical spiking neurosynaptic networks with self-learning capabilities,” Nature 569(7755), 208–214 (2019). [CrossRef]  

39. K. Du, Q. Li, Y. Lyu, J. Ding, Y. Lu, Z. Cheng, and M. Qiu, “Control over emissivity of zero-static-power thermal emitters based on phase-changing material GST,” Light-Sci. Appl. 6(1), e16194 (2017). [CrossRef]  

40. A. Tittl, A. U. Michel, M. Schäferling, X. Yin, B. Gholipour, L. Cui, M. Wuttig, T. Taubner, F. Neubrech, and H. Giessen, “A Switchable Mid-Infrared Plasmonic Perfect Absorber with Multispectral Thermal Imaging Capability,” Adv. Mater. 27(31), 4597–4603 (2015). [CrossRef]  

41. Y. Chen, T. Kao, B. Ng, X. Li, X. Luo, B. Luk’yanchuk, A. Maier, and M. Hong, “Hybrid phase-change plasmonic crystals for active tuning of lattice resonances,” Opt. Express 21(11), 13691 (2013). [CrossRef]  

42. A. K. U. Michel, P. Zalden, D. N. Chigrin, M. Wuttig, A. M. Lindenberg, and T. Taubner, “Reversible Optical Switching of Infrared Antenna Resonances with Ultrathin Phase-Change Layers Using Femtosecond Laser Pulses,” ACS Photonics 1(9), 833–839 (2014). [CrossRef]  

43. C. Ríos, M. Stegmaier, P. Hosseini, D. Wang, and T. Scherer, “Integrated all-photonic non-volatile multi-level memory,” Nat. Photonics 9(11), 725–732 (2015). [CrossRef]  

44. D. Zhu, W. Shen, and H. Zhen, “Anisotropic optical constants of in-plane oriented polyfluorene thin films on rubbed substrate,” J. Appl. Phys. 106(8), 084504 (2009). [CrossRef]  

45. IR Transmission Spectra and Gemini Observatory. http://www.gemini.edu/?q=node/10789.

46. M. N. Julian, C. Williams, S. Borg, S. Bartram, and H. Kim, “Reversible optical tuning of GeSbTe phase-change metasurface spectral filters for mid-wave infrared imaging,” Optica 7(7), 746 (2020). [CrossRef]  

47. Z. Wang, T. S. Luk, Y. X. Tan, D. X. Ji, M. Zhou, Q. Q. Gan, and Z. F. Yu, “Tunneling enabled spectrally selective thermal emitter based on flat metallic films,” Appl. Phys. Lett. 106(10), 101104 (2015). [CrossRef]  

48. J. S. Tyo, D. L. Goldstein, D. B. Chenault, and J. A. Shaw, “Review of passive imaging polarimetry for remote sensing applications,” Appl. Opt. 45(22), 5453–5469 (2006). [CrossRef]  

49. L. Peng, D. Liu, H. Cheng, S. Zhou, and M. Zu, “A Multilayer Film Based Selective Thermal Emitter for Infrared Stealth Technology,” Adv. Optical Mater. 6(23), 1801006 (2018). [CrossRef]  

50. H. Lu, G. Wang, and X. Liu, “Manipulation of light in MIM plasmonic waveguide systems,” Sci. Bull. 58(30), 3607–3616 (2013). [CrossRef]  

51. Y. Qu, Q. Li, K. Du, L. Cai, J. Lu, and M. Qiu, “Dynamic Thermal Emission Control Based on Ultrathin Plasmonic Metamaterials Including Phase-Changing Material GST,” Laser Photonics Rev. 11(5), 1700091 (2017). [CrossRef]  

52. X. Liu, T. Tyler, T. Starr, A. F. Starr, N. M. Jokerst, and W. J. Padilla, “Taming the blackbody with infrared metamaterials as selective thermal emitters,” Phys. Rev. Lett. 107(4), 045901 (2011). [CrossRef]  

53. X. Jiang, D. Chen, Z. Zhang, J. Huang, K. Wen, J. He, and J. Yang, “Dual-channel optical switch, refractive index sensor and slow light device based on a graphene metasurface,” Opt. Express 28(23), 34079–34092 (2020). [CrossRef]  

54. X. Jiang, Z. Zhang, K. Wen, G. Li, J. He, and J. Yang, “A triple-band hybridization coherent perfect absorber based on graphene metamaterial,” Appl. Sci. 10(5), 1750 (2020). [CrossRef]  

55. K. Shportko, S. Kremers, M. Woda, D. Lencer, J. Robertson, and M. Wuttig, “Resonant bonding in crystalline phase-change materials,” Nat. Mater. 7(8), 653–658 (2008). [CrossRef]  

56. J. Park, S. Eom, and H. Lee, “Optical properties of pseudobinary GeTe, Ge2Sb2Te5, GeSb2Te4, GeSb4Te7, and Sb2Te3 from ellipsometry and density functional theory,” Phys. Rev. B 80(11), 115209 (2009). [CrossRef]  

57. J. Tian, H. Luo, Y. Yang, F. Ding, Y. Qu, D. Zhao, and S. I. Bozhevolnyi, “Active control of anapole states by structuring the phase-change alloy Ge2Sb2Te5,” Nat. Commun. 10(1), 396 (2019). [CrossRef]  

58. H. Li, “Refractive index of alkaline earth halides and its wavelength and temperature derivatives,” J. Phys. Chem. Ref. Data 9(1), 161–290 (1980). [CrossRef]  

59. Z. Zhang, J. Yang, W. Bai, Y. Han, X. He, J. Hyang, D. Chen, S. Xu, and W. Xie, “All-optical switch and logic gates based on hybrid silicon-Ge2Sb2Te5 metasurfaces,” Appl. Opt. 58(27), 7392 (2019). [CrossRef]  

60. E. D. Palik, Handbook of Optical Constants of Solids (Academic Press, 1985).

61. M. Li, D. Liu, H. Cheng, L. Peng, and M. Zu, “Manipulating metals for adaptive thermal camouflage,” Sci. Adv. 6(22), eaba3494 (2020). [CrossRef]  

62. Y. Qu, Q. Li, L. Cai, M. Pan, P. Ghosh, K. Du, and M. Qiu, “Thermal camouflage based on the phase-changing material GST,” Light Sci Appl 7(1), 1–10 (2018). [CrossRef]  

63. Y. Qu, Q. Li, L. Cai, and M. Qiu, “Polarization switching of thermal emissions based on plasmonic structures incorporating phase-changing material Ge2Sb2Te5,” Opt. Mater. Express 8(8), 2312–2320 (2018). [CrossRef]  

64. M. Song, D. Wang, Z. A. Kudyshev, Y. Xuan, Z. Wang, A. Boltasseva, V. M. Shalaev, and A. V. Kildishev, “Enabling Optical Steganography, Data Storage, and Encryption with Plasmonic Colors,” Laser. Photonics Rev. 15(3), 2000343 (2021). [CrossRef]  

65. S. Chandra, D. Franklin, J. Cozart, A. Safaei, and D. Chanda, “Adaptive Multispectral Infrared Camouflage,” ACS Photonics 5(11), 4513–4519 (2018). [CrossRef]  

66. Y. Zhang, J. B. Chou, J. Li, H. Li, Q. Du, A. Yadav, S. Zhou, M. Y. Shalaginov, Z. Fang, H. Zhong, C. Roberts, P. Robinson, B. Bohlin, C. Ríos, H. Lin, M. Kang, T. Gu, J. Warner, V. Liberman, K. Richardson, and J. Hu, “Broadband transparent optical phase change materials for high-performance nonvolatile photonics,” Nat. Commun. 10(1), 4279 (2019). [CrossRef]  

67. K. Du, L. Cai, H. Luo, Y. Lu, J. Tian, Y. Qu, P. Ghosh, Y. Lyu, Z. Cheng, M. Qiu, and Q. Li, “Wavelength-tunable mid-infrared thermal emitters with nonvolatile phase changing material,” Nanoscale 10(9), 4415–4420 (2018). [CrossRef]  

68. T. Cao and M. J. Cryan, “Study of incident angle dependence for dual-band double negative-index material using elliptical nanohole arrays,” J. Opt. Soc. Am. A 29(3), 209 (2012). [CrossRef]  

69. T. Cao, L. Zhang, R. E. Simpson, and M. J. Cryan, “Mid-infrared tunable polarization-independent perfect absorber using a phase-change metamaterial,” J. Opt. Soc. Am. B 30(6), 1580 (2013). [CrossRef]  

70. J. A. Bossard, L. Lin, S. Yun, L. Liu, D. H. Werner, and T. S. Mayer, “Near-ideal optical metamaterial absorbers with super-octave bandwidth,” Acs Nano 8(2), 1517–1524 (2014). [CrossRef]  

Supplementary Material (1)

NameDescription
Supplement 1       Discussion of inverse design and initial structure

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1.
Fig. 1. (a) Atmospheric transmittance (light-blue shaded area) with detected band (3−5 µm and 8−14 µm) and undetected band (5−8 µm) indicated and radiative heat flux at 473 K. (b) Top view of one unit cell of the design has been divided into 50 × 50 pixels. (The yellow pixel represents gold, and the purple pixel represents air, and the GST is exposed) (c) Schematic of the PCM-based MSE. The MSE is developed with inverse design gold grating coated on the GST and gold mirror. The periodic structural parameters are P = 1 µm. The thicknesses of inverse design gold grating (t1), GST dielectric layer (t2) and gold mirror (t3) are 40, 230 and 300 nm, respectively.
Fig. 2.
Fig. 2. Designed schematic of the DBS algorithm, initial structure, and optimal metasurface. The steps of the DBS algorithm design include choosing pixel, switching the state of pixel, calculating FOM, judging the spectral characteristic according to the FOM, and outputting the optimal topological structure. The initial structure and the optimal metasurface are given in dashed frames. The initial point position and scanning order are indicated in the dashed frame initial structure diagram.
Fig. 3.
Fig. 3. Dielectric permittivity (ɛGST) of crystalline and amorphous GST.
Fig. 4.
Fig. 4. (a) Simulated spectral emissivity at normal incidence and spectral emissivity of ideal thermal selective emitter. (b) Electric field distribution in the indicated xy- and xz- cross sections at the two peak wavelengths.
Fig. 5.
Fig. 5. (a) The relation between radiation temperatures and real temperatures of the MSE: the lines are calculated with the integrated power. (b) The spectral radiation intensity of the blackbody and MSE at 473 K, the black and red lines represent the blackbody and the MSE, respectively.
Fig. 6.
Fig. 6. Simulated spectral emissivity of the plasmonic thermal emitter gradually changes as crystallization fraction of GST increases from 0% to 100%. The colormaps of PCM-based MSE with different crystallization fractions comparing the MWIR ATWs (the integral transmittance is 67% in the wavelength range of 3–5 µm) have been marked in each emissivity spectrum.
Fig. 7.
Fig. 7. (a) Emissivity of the MSE in the range of 3–14 µm, and (b) MSE is insensitive to the angle and polarization.

Tables (1)

Tables Icon

Table 1. Performance of the selective emitters for IR stealth with thermal management under two types of evaluations

Equations (7)

Equations on this page are rendered with MathJax. Learn more.

F O M = α 1 × λ 1 λ 2 ε ( λ ) d λ λ 2 λ 1 + α 2 × ( 1 λ min λ 1 ε ( λ ) d λ λ 1 λ min ) + α 3 × ( 1 λ 2 λ max ε ( λ ) d λ λ max λ 2 ) ,
ε G S T ( λ , C ) 1 ε G S T ( λ , C ) + 2 = C × ε c G S T ( λ , C ) 1 ε c G S T ( λ , C ) + 2 + ( 1 C ) × ε a G S T ( λ , C ) 1 ε a G S T ( λ , C ) + 2 ,
T r = P 1 ( ε I R , T ) .
P ( ε , T ) = P r a d ( ε , T ) + P r e f ( ε , ε a , T a ) = ε ( λ ) I B B ( T ) + [ 1 ε ( λ ) ] ε a ( λ ) I B B ( T a ) ,
e f f = e f f M W I R i n e f f o u t e f f L W I R i n .
e f f M W I R i n = 1 λ min λ 1 ε ( λ ) d λ λ 1 λ min , e f f o u t = λ 1 λ 2 ε ( λ ) d λ λ 2 λ 1 ,
e f f L W I R i n = 1 λ 2 λ max ε ( λ ) d λ λ max λ 2 .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.