Transparent and Ultra-Thin Flexible Checkerboard Metasurface for Radar–Infrared Bi-Stealth

This paper proposed a single-layer checkerboard metasurface with simultaneous wideband radar cross-section (RCS) reduction characteristics and low infrared (IR) emissivity. The metasurface consists of an indium tin oxide (ITO)-patterned film, a polyethylene terephthalate (PET) substrate and an ITO backplane from the top downwards, with a total ultra-thin thickness of 1.6 mm. This design also allows the metasurface to have good optical transparency and flexibility. Based on phase cancellation and absorption, the metasurface can achieve a wideband RCS reduction of 10 dB from 10.6 to 19.4 GHz under normal incidence. When the metasurface is slightly cylindrically curved, an RCS reduction of approximately 10 dB can still be achieved from 11 to 19 GHz. The polarization and angular stability of the metasurface have also been verified. The filling rate of the top ITO-patterned film is 0.81, which makes the metasurface have a low theoretical IR emissivity of 0.24. Both simulation and experimental results have verified the excellent characteristics of the proposed checkerboard metasurface, demonstrating its great potential application in radar–IR bi-stealth.


Introduction
The development of modern radar detection technology has brought great challenges to the stealth characteristics of military targets.Radar cross-section (RCS) is an important indicator to measure the stealth capability of targets, and reducing RCS can effectively improve the survivability of targets [1,2].Applying radar-absorbing materials (RAMs) and redirecting scattering fields are two commonly used methods for RCS reduction [3,4].RAMs can convert incident energy to heat, and have been widely studied for their advantage of convenience and efficiency.Although some broadband RAMs have been proposed, they have the disadvantage of being too thick or too complex in structure [5,6].Changing the shape of targets to keep the scattered field away from the direction is another feasible way to reduce RCS, but this may seriously affect the mechanical performance of targets.
Metasurfaces have the powerful ability to manipulate electromagnetic (EM) waves, which provides new strategies to achieve RCS reduction [7][8][9][10].Through different design methods, metasurfaces can suppress scattering field through absorption, phase cancellation, and polarization conversion [11][12][13].Among them, the checkerboard metasurface, based on phase interference cancellation technology, makes it easier to effectively reduce RCS in a wide frequency band, while avoiding thermal accumulation [14,15].The main principle of a checkerboard metasurface is to use different types of artificial magnetic conductor (AMC) patterns and to make their reflection phase difference equal approximately 180 • .The simplest structure is to use two types of AMC units with a reflection phase difference within 180 • ± 37 • , and many structures based on this principle have been designed [16,17].Y. Modi et al. designed a blended checkerboard metasurface using square rings and circular patch elements, which could achieve 10 dB RCS reduction in 3.9-9.45GHz [18].F. El-Sewedy et al. proposed a thin metasurface composed of triangular patches with different sizes, which had a thickness of only 0.8 mm but could achieve a 10 dB RCS reduction in 23.71-33.52GHz [19].Random and coded checkerboard metasurfaces have also been extensively studied, and through the algorithmic optimization of the layout, the frequency range of RCS reduction can be greatly expanded [20][21][22].In [23], a binary-coded metasurface was designed, based on machine learning optimisation techniques, to achieve 10 dB RCS reduction in 14-20 GHz.Qu et al. proposed an optimized multi-element phase cancellation (OMEPC) method to design checkerboard metasurfaces [24].By optimizing the shape and size of various unit cells, a 10 dB RCS reduction can be achieved in the ultra-wideband frequency band of 7.4-64.8GHz.
With the increasing application of multiple detection, materials and structures with radar-IR bi-stealth are gradually attracting attention.Some composites with both microwave absorption and thermal insulation ability have been proposed to achieve radar-IR integrated stealth, such as Ni-MXene/Melamine [25], Fe/Fe 2 O 3 porous carbon [26].However, these composite materials are difficult to fabricate and have narrow absorption bands.Since metals usually have low IR emissivity, covering metamaterial absorbers with high filling rate metal patches is an effective method to achieve radar-IR bi-stealth, and many structures based on this method have been proposed [27][28][29][30][31].It should be noted that covering metal patches will change the impedance matching of metamaterial absorbers, thus changing the original absorbing effect.For special applications such as aircraft windows and satellite solar panels, radar-IR bi-stealth metamaterial absorber with simultaneous optical transparency is designed [32][33][34].However, metamaterial absorbers convert absorbed energy into thermal energy, which is not conducive to IR stealth.In addition, these structures are usually multi-layered, which suffers from the disadvantages of being too thick and less flexible.Some studies have proposed using checkerboard metasurface structures with high metal filling rates to achieve radar-IR stealth, but these structures have the disadvantages of large IR emissivity and thickness [35][36][37].
In this work, a single-layer checkerboard metasurface with simultaneous optical transparency, ultra-thin, flexible, wideband RCS reduction, and low IR emissivity is proposed.The proposed structure consists of a top indium tin oxide (ITO)-patterned film, a middle polyethylene terephthalate (PET) substrate, and an ITO backplane.The total structural thickness is only 1.6mm, which allows for good flexibility.The top ITO-patterned film contains two types of units with high ITO filling rates and, therefore, has lower IR emissivity.Finally, a sample of the proposed checkerboard metasurface was fabricated, and its RCS reduction performance and IR emissivity were verified through microwave darkroom RCS testing and thermal IR imaging.

Design and Analysis
In order to have good optical transparency and flexibility, the designed metasurface is only composed of a top ITO-patterned film, a PET substrate, and an ITO backplane.The IR radiation intensity of the metasurface can be calculated by the Stefan-Boltzmann equation as j = εσT 4 , where ε is the IR emissivity, σ is the Stefan constant, and T is the temperature.As the temperature T of the metasurface is difficult to control, the IR radiation can be reduced by reducing ε.The IR emissivity ε ITO of ITO increases with the increase of sheet resistance R. When R is less than 10 Ω/sq , ε ITO is about 0.09 [31].Since the operating wavelength range of the infrared band is only 3-12 µm, which is much lower than the radar wavelength, the IR characteristics of the metasurface mainly depend on the top ITO film pattern.The IR emissivity of PET is ε PET = 0.9.Then, the total IR emissivity can be calculated as , where f ITO is the filling rate of ITO.As the IR emissivity of PET is much larger than that of ITO, ε TOTAL decreases with the increase of f ITO .But when f ITO = one, the incident radar wave is almost fully reflected, which cannot achieve the effect of reducing RCS.Therefore, the top ITO pattern can be made into a bandpass periodic array with a high filling rate.Figure 1 shows the basic unit cell of the metasurface.The ITO pattern consists of four triangular patches with a sheet resistance of R 1 = 5 Ω/sq, and the period size is p.This design allows the pattern to have multiple gaps, which can generate more resonant modes.The relative permittivity of PET is approximately three, and its loss on the incident wave can be ignored [37].It can be calculated that f ITO and ε TOTAL are 0.81 and 0.24, respectively.The function of ITO backplane with a sheet resistance of R 2 is to reduce the transmittance of radar waves.It can be inferred that the radar wave transmittance decreases with the decrease of R 2 .When R 2 is zero, the ITO backplane is equivalent to an perfect electric conductor (PEC) surface, and there will be no transmission at this time.The sheet resistance of ITO backplane is select as 5 Ω/sq, which can make the transmittance lower than −20 dB.
Figure 2 shows the reflection amplitude and phase of the basic unit at different period sizes p.The frequency range of the incident wave is 0-23 GHz, and the thickness of PET substrate is h PET = 1.6 mm.These results were obtained using CST Studio Suite 2022 software.The frequency domain solver is adopted, where the "unit cell" boundaries were set in the xy plane and the "add open space" boundaries were set in the z direction.The reflection amplitude value is S 11 .Due to the absorption effect of ITO, it can be seen that the reflectivity is not one.As the increase of p, the reflection amplitude curve generates more resonances, with the lowest reflection amplitude of 0.54 at p = 10 mm.Therefore, a 10 dB RCS reduction will not be achieved by relying on the absorption of one particular size of unit.From Figure 2b, it can be seen that there are significant differences in the reflection phase of units with different period sizes.Therefore, units of different sizes can be arranged in a checkerboard shape, and the principle of reflect phase cancellation can be used to achieve wideband RCS reduction.To make the metasurface have better polarization stability, two basic units of different sizes are used.Figure 3 shows a period of the proposed checkerboard metasurface, and also illustrates its working mechanism, which can scatter radar waves in multiple directions to reduce RCS in a specific direction, while having low IR emissivity and being able to penetrate visible light.The metasurface contains two different basic units A1 and A2, with periodic sizes of p 1 and p 2 , respectively.For a checkerboard metasurface containing two types of units, an approximated expression for RCS reduction is presented as in [26], as follows: where E r1 and E r2 are the reflected amplitudes of A1 and A2, while φ r1 and φ r2 are their phases, respectively.By substituting the relevant parameters into (1), the final optimized p 1 and p 2 are 2 mm and 8 mm, respectively.The total metasurface has a period length of P = 80 mm.The period size of A2 is four times that of A1, so the 4 × 4 size A1 unit array can also be regarded as a basic unit.A1 and A2 form 20 × 20 and 5 × 5 arrays respectively, and are distributed crosswise.Figure 4a shows the reflectivity R(ω) and absorptivity A(ω) of the two units A1 and A2, while Figure 4b shows their corresponding reflection phases and their phase difference.The reflectivity is calculated by R = |S 11 | 2 and the absorptivity is A = 1 − R. It can be seen that both A1 and A2 have a certain absorption effect, where A1 generates a resonance point with a maximum absorption rate of 43% at 18.2 GHz, and A2 has two resonance points with a maximum absorption rate of 88% at 10.5 GHz.The phase difference between A1 and A2 at 11.6 GHz and 16.9 GHz is close to 180 • , so it is expected to generate two reflection reduction peaks.Under normal incidence, Figure 5a shows the RCS reduction obtained by a full-wave simulation and the approximate RCS reduction obtained by substituting the reflection coefficient and phase into (1).The full-wave simulation is calculated using the time domain solver of CST.The excitation source is a Gaussian pulse, and the RCS results in a wide frequency band are obtained through a time-frequency conversion.The time domain solver uses hexahedral meshes for spatial discretization, which cannot perfectly fit the geometric characteristics of the metasurface.Therefore, there are also certain errors in the simulation results.The metasurface is placed in the xy plane, and the angles between the incident wave direction with the z-axis and x-axis are θ and φ, respectively.The backscattering of the metasurface and the PEC surface with the same size is simulated, respectively, and the RCS reduction of the metasurface is obtained by subtracting the two results.It can be seen that the full-wave simulation results are basically consistent with the theoretical results obtained by (1).Considering the numerical error in the full-wave simulation, the difference between the two results is acceptable.The simulated RCS has a reduction of over 10 dB in the bandwidth of 10.6-19.4GHz, with a relative bandwidth of 58.6%.As predicted, there are two RCS reduction peaks, reaching 18.9 dB and 19.1 dB at 11.8 GHz and 16.8 GHz, respectively.Figure 5b shows the energy distribution ratio of the proposed checkerboard metasurface under normal incidence, where the energy ratio of absorption A(ω), phase cancellation PC(ω), and reflection R(ω) are calculated as follows: When the total energy ratio of the absorption and phase interference exceeds 90%, the reflected energy ratio will be less than 10%, and the RCS reduction will exceed 10 dB.It can be concluded that the achievement of broadband RCS reduction is the result of the combined effect of absorption and phase cancellation, and the contribution of phase cancellation to RCS reduction is greater than that of absorption.Figure 6 shows the three-dimensional (3D) scattering patterns of the checkerboard metasurface and equal-sized PEC surface at frequencies of 11.8 and 16.8 GHz under normal incidence, as well as the corresponding bistatic RCS in the xz plane.Due to the symmetry of the metasurface, the scattering patterns at this time are not affected by the polarization mode.It can be seen that the main scattering lobe of the PEC surface is concentrated in the normal incident direction.Four scattering lobes appear on the metasurface, and the scattering intensity is much lower than that of PEC surface.The bistatic RCS of the metasurface in the xz plane is lower than that of the PEC surface.In the direction of θ = 0 • , there is a significant decrease in the scattering intensity of the metasurface.The backscattering RCS of the PEC surface and the metasurface at 11.8 GHz are −0.9 dB and −19 dB, and the backscattering RCS of the PEC surface and metasurface at 16.8 GHz are 2.2 dB and −25 dB.These results indicate that the metasurface has good RCS reduction performance at normal incidence.Since the metasurface has an excellent flexibility, it can be cylindrically curved with a curvature angle of γ. Figure 7 shows the RCS reduction of the metasurface under different polarizations with γ values of 5 • , 10 • , 15 • , and 20 • at 15 GHz.It can be seen that for TM and TE polarization, as the curvature angles γ increases, the RCS reduction effect gradually decreases.When γ is less than 15 • , an RCS reduction of approximately 10 dB can still be achieved in the frequency range of 11-19 GHz.This is because the curvature changes the impedance match of the base unit, thereby changing the reflection amplitude and phase.For TM polarization, 10 dB RCS reduction can be achieved in a wide frequency band when γ is less than 15 • , which has a better performance than TE polarization.This is because, under TM polarization, the surface current induced by the magnetic field is still parallel to the metasurface, and a more stable scattering pattern is generated [34].Figure 8 shows the 3D scattering patterns of the PEC surface and the checkerboard metasurface with γ values of 10 • and 20 • , as well as the bistatic RCS in the xz plane.The incident frequency is 15 GHz and the polarization mode is TM polarization.It can be seen that the main scattering lobe of the PEC surface is still concentrated in the normal incident direction, and four scattering lobes appear on the metasurface.In addition, it can be seen that the main scattering lobes have a wider range of azimuths when γ = 20 • than when γ = 10 • .At this time, the bistatic RCS of the metasurface in the xz plane is still lower than the PEC surface.The backscattering RCS of the PEC surface and metasurface at γ = 10 • are 0.9 dB and −12.8 dB, and the backscattering RCS at γ = 20 • are −0.4 dB and −25 dB.Therefore, the designed checkerboard metasurface still has excellent RCS reduction performance after proper bending.The RCS reduction effect of the checkerboard metasurface at different oblique incidence angles are also studied.When electromagnetic waves are obliquely incident on the PEC surface, the reflected lobes will mainly concentrate in the specular direction, so the RCS reduction in the specular direction is calculated here.Under oblique incidence, the surface impedance and phase difference of A1 and A2 are functions of the incidence angle.Figure 9 shows the RCS reduction results under different polarizations at incidence angles θ of 15 • , 30 • , 45 • , and 60 • .It can be seen that, when θ is less than 45 • , the checkerboard metasurface can still achieve 10 dB RCS reduction in a wide frequency band under TE and TM polarization.For TM polarization, broadband 10 dB RCS reduction can still be achieved when θ is greater than 45 • .However, for TE polarization, the RCS reduction performance deteriorates sharply when θ is greater than 45 • .This is also because, under TM polarization, the magnetic field is parallel to the metasurface and, thus, generates a sufficiently large induced current, which is more conducive to the absorption of the incident wave.Also, for TE polarization, the magnetic field component parallel to the metasurface decreases with an increasing incident angle, ultimately leading to a decrease in absorption and RCS reduction effects.Figure 10 shows the 3D scattering patterns of PEC surface and the checkerboard metasurface with θ values of 30 • and 60 • , as well as the bistatic RCS in the xz plane.The incident frequency is 15 GHz and the polarization mode is TM polarization.It can also be seen that the PEC surface generates one main scattering lobe in the specular direction, while the metasurface exhibits four scattering lobes.Compared to the PEC surface, the bistatic RCS of the metasurface in the xz plane is also significantly reduced.

Experimental Validation
A sample of the metasurface consisting of a 2 × 2 period was fabricated and measured.First, a 180 nm thick ITO coating was deposited on a 0.05 mm thick PET film by magnetron sputtering technology, to obtain an ITO-PET film, which can be directly used as the ITO backplane of the metasurface.The ITO coating thickness had little effect on the total thickness of the film.Then, the ITO-PET film obtained in the previous step was processed using laser etching technology to obtain the ITO pattern on the top layer of the metasurface.The actual processed ITO film resistance is approximately 4-6 Ω/sq.The surface current of ITO pattern and the charge accumulation on the edge will generate inductive reactance and capacitive reactance, respectively.According to the method in [38], it can be calculated that the inductive reactance and capacitive reactance of the ITO pattern are much larger than the surface resistance, so the ITO film resistance error has little impact on the results.In addition, the full wave simulation results also verified that the difference in surface resistance had almost no impact on the results.Finally, we used transparent adhesive to stick the upper and lower layers of ITO-PET film onto a 1.5 mm thick PET substrate, to obtain the sample.From Figure 11a, the scene behind can be clearly seen through the sample, which verifies the good optical transparency of the checkerboard metasurface.The sample has good flexibility and can be easily cylindrically curved, as shown in Figure 11b.Figure 11c shows the magnified photograph image of the fabricated sample.The transparency of the metasurface in the visible wavelength range of 400-800 nm is measured using a ultraviolet (UV)/visible (Vis)/near-infrared (NIR) spectrometer, and Figure 11d shows the measurement results.Due to its ultra-thin thickness of 1.6 mm, the sample exhibits a transmittance of over 80% in the visible light band.Figure 12a shows the experimental set up for RCS measurement in a microwave darkroom, in which two standard horn antennas are connected to a vector network analyzer (N9918A), and the operating frequency range of antennas is 8-18 GHz.The transmitting antenna emits signals and converts them into plane waves through a metal reflective screen.The scattered waves are reflected by the metal screen and received by the receiving antenna.Due to limitations in experimental conditions, only two cases under TE polarization at normal incidence were tested: when the metasurface was not curved and the metasurface has a cylindrically curved angle of γ = 20 • .Figure 12b shows the comparison of RCS reduction results obtained from experimental testing and simulation for the two cases.It can be seen that, although there are differences between the experimental results and the simulation results, the overall trend of change is basically consistent.The reason for the difference is the inability to precisely control the orientation and curvature angle of the metasurface.In addition, there are numerical errors in the full-wave simulation results.The IR stealth ability of the proposed metasurface was tested through IR imaging.We imaged the fabricated sample using a thermal IR camera (Thermoview T72) (Huashengchang, Shenzhen, China) operating at the wavelength of 7.5-14 µm and compared the results with the IR imaging of a metal plate and a PET plate.Figure 13a shows the image of the three samples.Figure 13b is the IR image at a room temperature of 25 • C, and Figure 13c is the IR image of the three samples after being heated on a 50 • C heating plate for ten minutes.The higher the temperature displayed in the IR image, the stronger the IR radiation of the object.It can be found that the IR radiation of the checkerboard metasurface is higher than the metal plate, but is significantly lower than PET, which confirms the good IR stealth performance of the metasurface.Based on the imaging results, it can be estimated that the IR emissivity of the checkerboard metasurface is as follows [39]: where T r is the temperature displayed in the IR image, T o is the ambient temperature of the metasurface, and T a is the atmospheric temperature.The ambient temperature is 50 • C, the atmospheric temperature is 25 • C, and the checkerboard metasurface in Figure 13c shows a temperature of about 32 • C.These three temperatures need to be converted into Kelvin temperatures when substituting into (3), and the result can be eventually calculated as ϵ exp = 0.256.This is very close to the theoretical result, proving that the metasurface has an extremely low IR emissivity.Table 1 shows a comprehensive comparison of the proposed checkerboard metasurface with some previously reported structures.The comparison shows that although the proposed metasurface has a relatively low 10 dB RCS reduction bandwidth, it has low thickness, optical transparency, low infrared emissivity, and flexibility at the same time.The thickness of the metasurface is only 0.08λ L , where λ L is the wavelength at the center frequency of the 10 dB RCS reduction band.Therefore, compared with previously pub-lished radar-IR bi-stealth structures, the proposed metasurface has significant advantages in thickness and flexibility.

Conclusions
In summary, a checkerboard metasurface with simultaneous optical transparency, ultra-thin, flexible, wideband RCS reduction, and low IR emissivity is designed.Based on the hybrid mechanism of phase cancellation and absorption, the designed metasurface can effectively achieve 10 dB RCS reduction in the frequency band of 10.6-19.4GHz, and has a low IR emissivity of 0.24.The total thickness of the metasurface is only 1.6 mm, which allows for good flexibility.The metasurface can still achieve a good RCS reduction effect under normal incidence when it is cylindrically curved.A sample of the checkerboard metasurface was fabricated and measured, and the experimental results also verified its excellent performance.The proposed checkerboard metasurface shows a good radar-IR bi-stealth capability, and has great potential for applications such as aircraft canopies and stealth glass.

Figure 1 .
Figure 1.Unit cell design of the checkerboard metasurface.

Figure 2 .
Figure 2. (a) Reflection amplitude and (b) reflection phase of basic units under different period sizes.

Figure 3 .
Figure 3.One period of the checkerboard metasurface and the schematic of RCS reduction, low IR emissivity, and optical transparency.

Figure 4 .
Figure 4. (a) Simulated reflectivity and absorptivity of A1 and A2.(b) Simulated reflection phases of A1, A2, and the phase difference between them.The vertical axis corresponding to the black dashed line is 180 • .

Figure 5 .
Figure 5. (a) RCS reduction obtained by full wave simulation and approximate RCS reduction obtained by (1) under normal incidence.(b) The energy distribution ratio of the metasurface, and the vertical axis corresponding to the black solid line is 0.9.

Figure 6 .
Figure 6.The 3D scattering patterns of the PEC surface and the checkerboard metasurface, as well as the corresponding bistatic RCS in the xz plane at different frequencies.(a-c) 11.8 GHz.(d-f) 16.8 GHz.

Figure 8 .
Figure 8.The 3D scattering patterns of the PEC surface and the checkerboard metasurface, as well as the corresponding bistatic RCS in the xz plane at different curvature angles.(a-c) γ = 10 • .(d-f) γ = 20 • .

Figure 9 .
Figure 9. Simulated RCS reduction in specular direction of the checkerboard metasurface under different incident angles θ and polarizations.(a) TM polarization.(b) TE polarization.

Figure 10 .
Figure 10.The 3D scattering patterns of the PEC surface and the checkerboard metasurface, as well as the corresponding bistatic RCS in the xz plane at different oblique incidence angles.(a-c) θ = 30 • .(d-f) θ = 60 • .

Figure 11 .
Figure 11.(a) Photograph of the fabricated sample showing optical transparency and (b) photograph of the fabricated sample showing flexibility.(c) Magnified photograph image of the fabricated sample.(d) The transmittance of the sample measured in the visible light wavelength range.

Table 1 .
A comprehensive comparison of the proposed checkerboard metasurface with some previously reported structures.Relative bandwidth; b the wavelength corresponding to the center frequency of the 10 dB RCS reduction band. a