Dual-Mechanism Absorptive Metasurface with Wideband 20 dB RCS Reduction

: This paper presents a dual-mechanism method to design a single-layer absorptive meta-surface with wideband 20 dB RCS reduction by simultaneously combining the absorption and phase cancellation mechanisms. The metasurface comprises two kinds of absorbing unit cells with 10 dB absorption performance but different reﬂection phases. The impedance condition for 20 dB RCS reduction is theoretically analyzed considering both the absorption and the phase cancellation based on the two unit cells, and the relationship between the surface impedance and the reﬂection phase/amplitude is revealed. According to these analyses, two unit cells with absorption performance and different reﬂection phases are designed and utilized to realize the absorptive metasurface. Numerical and experimental results show that the single-layer absorptive metasurface features wideband 20 dB RCS within 11.5–16 GHz with a thickness of only 3 mm.

The phase cancellation mechanism of the metasurface is a spatial domain RCS reduction method that is to properly design the reflection phase of each unit cell, realize phase cancellation or phase gradient, and then control the reflected wave beam distribution.The state-of-the-art approaches are mainly based on the chessboard distribution phase cancellation metasurfaces [11,12,14,15] and the optimized coded metasurfaces [16][17][18][19].The traditional chessboard metasurfaces redirect the reflected wave into four symmetrical beams, making the reflection amplitude highly suppressed in the monostatic direction.The optimized coding metasurfaces shape the reflected wave with approximately uniform echo intensity in all directions, which can also achieve RCS reduction.The absorption mechanism of the metasurface is a method to reduce the RCS in the energy domain, which mainly utilizes the lossy property of the unit cell.The absorbing metasurface can be divided into various forms, including resistive loss, dielectric loss, and magnetic loss [27][28][29][30].
Some researchers use both cancellation and absorption to reduce RCS [31][32][33][34][35][36].In [31], the absorption occurs in low frequency using lumped resistors.The cancellation diffusion occurs in high frequency by constructing random phase distribution, resulting in the dualband 8 dB RCS reduction.In [32], metasurfaces with multiple layers that operate phase cancellation and absorption are proposed to achieve dual-mechanism and obtain 10 dB RCS reduction from 1.55 to 19.2 GHz, respectively.Similarly, Ref. [33] demonstrates a double-layered metasurface to acquire wideband RCS reduction by introducing both phase cancellation and absorption.It should be noted that almost all the works above focus on the expansion of reduction bandwidth using the dual-mechanism but fail to explore the possibility of realizing a lower RCS index.In [36], the absorption is introduced into phase diffusion metasurface to broaden the 10 dB RCS reduction, but fail to achieve lower 20 dB reduction.
Wider band and deeper RCS reduction by lower profile are always the main concerns of stealth material research.Though the above two RCS reduction mechanisms have made significant progress in expanding bandwidth, most researchers design and measure RCS reduction bandwidth with 10 dB standard, that is, the return energy is reduced by 90%.With the rapid development of radar detection technology, the 10 dB RCS reduction can no longer meet the stealth requirements.
A dual-mechanism absorptive metasurface with 20 dB RCS reduction by combining both the phase cancellation and absorption mechanisms is proposed in this paper.Utilizing the reflection phase difference between two absorbing unit cells, we organically combine the absorbing and phase cancellation mechanism to reduce the intensity of the backward echo in both the energy domain and the spatial domain.By optimizing the number of two unit cells, an absorptive metasurface is composed to implement RCS reduction above 20 dB.Numerical and experimental results confirm that the dual-mechanism absorptive metasurface manifests wideband, robustness, and much lower RCS.Compared to current methods, our strategy can broaden the 20 dB reduction bandwidth while keeping that of the original 10 dB almost unchanged, promising great potentials in real applications.

Design Principle of Metasurface
Suppose a metasurface, which consists of M × N unit cells, is under the normal incidence of a plane wave, as shown in Figure 1.The conception of RCS reduction can be explained by the principle of planar array theory [32].For a metasurface composed of two different unit cells, which is the simplest case for cancellation layout, the monostatic RCS reduction of metasurface can be represented by: If "1" and "2" have no absorption and all energy is reflected, then A1 = A2 = 1.The relationship between RCS reduction and parameter ,   can be expressed by Figure 2. The maximum range of   appears at  = 0.5.In order to obtain 10 dB RCS reduction, Here, A 1 and A 2 represents the amplitude of the reflected EM waves of unit cells "1" and "2".Parameter α is the ratio of two unit cells and α = n 1 n 1 +n 2 , 1 − α = n 2 n 1 +n 2 , where n 1 is the number of unit cell "1" in whole metasurface while n 2 is that of unit cell "2".The phase difference of two unit cells is ϕ d , which can be expressed as: ϕ 1 and ϕ 2 are the reflection phases of "1" and "2", respectively.If "1" and "2" have no absorption and all energy is reflected, then A 1 = A 2 = 1.The relationship between RCS reduction and parameter α, ϕ d can be expressed by Figure 2. The maximum range of ϕ d appears at α = 0.5.In order to obtain 10 dB RCS reduction, the phase difference needs to be at (180    ), which is consistent with the conclusion of AMC chessboard [14].For RCS reduction above 20 dB, the phase difference needs to be at (180    ).However, for the resonance phase, it is an extremely harsh constraint for two unit cells to meet the phase difference of 180 • ± 11 • .Therefore, the reduction effect using phase cancellation, which could get more than 20 dB RCS reduction, is only for individual frequency points.If "1" has 10 dB absorption property while "2" is lossless (A 1 = 0.316, A 2 = 1), the relationship can be expressed by Figure 3a.The 20 dB maximum range of ϕ d appears at α = 0.78 from 180 • − 24 • to 180 • + 24 • , which is wider than 180 • ± 11 • .If "1" and "2" both have 10 dB absorption property (A 1 = A 2 = 0.316) (see Figure 3b), the 20 dB maximum range of ϕ d appears at α = 0.5 and the range from 180 • −       To sum up, absorption can help to expand the range of   of specific RCS reduction value, and the  is 0.5 when absorption ability of "1" and "2" are the same.When the absorption coefficient is unequal, the  will move to an off-center position.The degree of   To sum up, absorption can help to expand the range of   of specific RCS reduction value, and the  is 0.5 when absorption ability of "1" and "2" are the same.When the absorption coefficient is unequal, the  will move to an off-center position.The degree of deviation increases with the increase in absorption coefficient difference.To sum up, absorption can help to expand the range of ϕ d of specific RCS reduction value, and the α is 0.5 when absorption ability of "1" and "2" are the same.When the absorption coefficient is unequal, the α will move to an off-center position.The degree of deviation increases with the increase in absorption coefficient difference.

Unit Cell Design
Based on the equivalent transmission line method, a single-layer reflected unit cell, which is an infinite periodic structure in free space (see Figure 4a), can be equivalent to a single-port transmission-line network.As shown in Figure 4b, Zr and Zd denote the equivalent impedance of the periodic structure and the substrate layer, respectively.Zr is a complex value that can be expressed by Zr = R + jX.Zin is the input impedance of this network.The reflection coefficient of the unit cell can be expressed as And =  0 /√  (5) where Z0 is the intrinsic impedance of a vacuum, d and   are the thickness and relative permittivity of the substrate layer, respectively.Thus, the reflection amplitude and phase can be expressed as Here, we fix d = 3 mm,   = 3, and the center frequency is set to be 13 GHz.The relationship among ,  , R, and X can be shown in Figure 5.The black line indicates the phase of the unit cell, and the color map is the reflection amplitude.White dotted line represents reflection amplitude  = −10 dB.Therefore, we can select the points inside the white dotted circle to design the unit cells with proper phase and absorption accordingly.As shown in Figure 4b, Z r and Z d denote the equivalent impedance of the periodic structure and the substrate layer, respectively.Z r is a complex value that can be expressed by Z r = R + jX.Z in is the input impedance of this network.The reflection coefficient of the unit cell can be expressed as And where Z 0 is the intrinsic impedance of a vacuum, d and ε r are the thickness and relative permittivity of the substrate layer, respectively.Thus, the reflection amplitude and phase can be expressed as Here, we fix d = 3 mm, ε r = 3, and the center frequency is set to be 13 GHz.The relationship among A, ϕ, R, and X can be shown in Figure 5.The black line indicates the phase of the unit cell, and the color map is the reflection amplitude.White dotted line represents reflection amplitude A = −10 dB.Therefore, we can select the points inside the white dotted circle to design the unit cells with proper phase and absorption accordingly.
However, for two different unit cells "1" and "2", there are eight interrelated parameters, namely It is complex and time-consuming to optimize and select the eligible two unit cells.Thus, to simplify the process, we firstly design one absorption unit cell "1" and then choose the other counterpart.
According to the pre-defined parameters, the optimized topology of the absorbed "1" is shown in Figure 6a.It is a single-layer structure composed of top compound resistive film, a middle dielectric spacer, and a backed metallic ground.The top compound resistive film patterns are two square rings, which can offer two resonance points at the given frequency band.The resistive film will add resistive loss at these two frequencies which ensure the wideband absorption.In addition, the fourfold rotation symmetry makes the unit cells polarization insensitive.Additionally, the stable wideband reflection phase difference by varied dimensions has been verified in ref. [16].The resistive film can offer R and X simultaneously for 10 dB absorption, as shown in Figure 5.In the current design, the F4BM-1 substrate with permittivity of 3, a 1 = 5 mm, b 1 = a 1 − 3 mm, the height of h = 3 mm and dielectric loss of tan δ = 0.00081 is utilized as the dielectric spacer.For "1", the sheet resistance of resistive film patterns is 10 Ω/sq, and the reflection phase and amplitude curves are shown in Figure 6b.The 10 dB absorption bandwidth is 9.9-17.5GHz, but that of 20 dB is only 10.9-12.2GHz.However, for two different unit cells "1" and "2", there are eight interrelated parameters, namely A1, A2, P1, P1, X1, X2, R1, R2.It is complex and time-consuming to optimize and select the eligible two unit cells.Thus, to simplify the process, we firstly design one absorption unit cell "1" and then choose the other counterpart.
According to the pre-defined parameters, the optimized topology of the absorbed "1" is shown in Figure 6a.It is a single-layer structure composed of top compound resistive film, a middle dielectric spacer, and a backed metallic ground.The top compound resistive film patterns are two square rings, which can offer two resonance points at the given frequency band.The resistive film will add resistive loss at these two frequencies which ensure the wideband absorption.In addition, the fourfold rotation symmetry makes the unit cells polarization insensitive.Additionally, the stable wideband reflection phase difference by varied dimensions has been verified in ref. [16].The resistive film can offer R and X simultaneously for 10 dB absorption, as shown in Figure 5.In the current design, the F4BM-1 substrate with permittivity of 3, a1 = 5 mm, b1 = a1 − 3 mm, the height of h = 3 mm and dielectric loss of tan δ = 0.00081 is utilized as the dielectric spacer.For "1", the sheet resistance of resistive film patterns is 10 Ω/sq, and the reflection phase and amplitude curves are shown in Figure 6b.The 10 dB absorption bandwidth is 9.9-17.5GHz, but that of 20 dB is only 10.9-12.2GHz.Actually, the double square ring pattern is not the only choice.Any pattern has above features can play the role as unit cell.In this paper, our research attempts to provide a design recipe of whole metasurface instead of the unit cells.The unit cells we utilized are as an example to verify the effectiveness of our design recipe.Thus, we have not done too much of a complex design on the structure pattern and optimization on the performance of the unit cell.
Here, we need to confirm the design priority between absorption ability and the desired phase because these two indicators are difficult to be satisfied meanwhile for "2".Actually, the double square ring pattern is not the only choice.Any pattern has above features can play the role as unit cell.In this paper, our research attempts to provide a design recipe of whole metasurface instead of the unit cells.The unit cells we utilized are as an example to verify the effectiveness of our design recipe.Thus, we have not done too much of a complex design on the structure pattern and optimization on the performance of the unit cell.
Here, we need to confirm the design priority between absorption ability and the desired phase because these two indicators are difficult to be satisfied meanwhile for "2".Suppose "2" has 10 dB absorption but no phase difference with "1", a metasurface composed by "1" and "2" will only have 10 dB RCS reduction.While if "2" has no absorption but 180 • phase difference with "1", just like the situation shown in Figure 3a, the metasurface composed by these two cells will reach 20 dB RCS reduction and more.Thus, the phase difference is the key indicator than absorption for "2".It does not mean that the absorption ability is of no use.The absorption ability will expand the phase range of 20 dB RCS reduction, which could get the wideband bandwidth.
In order to obtain the structural parameters of "2" that meet the phase difference and absorption conditions, the impedance analysis of two square resistive film rings is carried out utilizing the CST de-embedding simulation.The curves of real part R and imaginary part X of Z r are shown in Figure 7.When the sheet resistance value remains unchanged and only the length of the outer square ring varies, the curves of the real part R are almost overlapped, but the imaginary part X varies obviously.In another case, when the parameter a remains unchanged, and only the sheet resistance value varies, an opposite result will happen.Thus, according to the rule shown in Figure 7 and the relationship shown in Figure 5, unit cell "2" that meets the phase and absorption requirements can be designed.Through the above analysis and further parameter optimization and scanning, unit cell "2" is obtained, and the topology is shown in Figure 6c.It has a similar structure to "1" but its sheet resistance and the dimension of resistive film patterns are different.The sheet resistance of resistive film patterns in "2" is 30 Ω/sq.The reflection phase and amplitude curves of "2" are shown in Figure 6d.As shown in Figure 6d, the unit cell 2 shows a relatively narrow absorption bandwidth, and the 10 dB bandwidth is only 11.8-16.3GHz.On the other hand, the phase response curve is approximately linear, which can get a stable phase difference with "1".The phase difference between "1" and "2" maintains about 180° from 11.5-15.5GHz, see Figure 8. Through the above analysis and further parameter optimization and scanning, unit cell "2" is obtained, and the topology is shown in Figure 6c.It has a similar structure to "1" but its sheet resistance and the dimension of resistive film patterns are different.The sheet resistance of resistive film patterns in "2" is 30 Ω/sq.The reflection phase and amplitude curves of "2" are shown in Figure 6d.As shown in Figure 6d, the unit cell 2 shows a relatively narrow absorption bandwidth, and the 10 dB bandwidth is only 11.8-16.3GHz.On the other hand, the phase response curve is approximately linear, which can get a stable phase difference with "1".The phase difference between "1" and "2" maintains about 180 • from 11.5-15.5GHz, see Figure 8.

Metasurface Design and Simulation
According to Equation (1) and Figure 3a, the varied absorption ability of unit cells at different frequencies will make the maximum range of ϕ d deviate.At each specific frequency, there is a maximum range of ϕ d related to the reflection amplitude of "1" and "2", see Figure 9a-c.To find the optimal ratio α of the whole metasurface in the wideband frequency, we need to overlap the −20 dB ranges of different frequency points.Here, an Crystals 2022, 12, 493 7 of 14 example is given in Figure 9d that the overlap area is calculated by three frequency points, 12, 13 and 14 GHz.The ordinate range corresponding to the blue overlapping area is the value range of α to meet the 20 dB reduction, see shades of gray in Figure 9d.The ordinate corresponding to the maximum phase difference range is the optimal α value of these three frequency points, see the red point in Figure 9d.Through the above analysis and further parameter optimization and scanning, unit cell "2" is obtained, and the topology is shown in Figure 6c.It has a similar structure to "1" but its sheet resistance and the dimension of resistive film patterns are different.The sheet resistance of resistive film patterns in "2" is 30 Ω/sq.The reflection phase and am plitude curves of "2" are shown in Figure 6d.As shown in Figure 6d, the unit cell 2 shows a relatively narrow absorption bandwidth, and the 10 dB bandwidth is only 11.8-16.3GHz.On the other hand, the phase response curve is approximately linear, which can ge a stable phase difference with "1".The phase difference between "1" and "2" maintains about 180° from 11.5-15.5GHz, see Figure 8.

Metasurface Design and Simulation
According to Equation (1) and Figure 3a, the varied absorption ability of unit cells a different frequencies will make the maximum range of   deviate.At each specific frequency, there is a maximum range of   related to the reflection amplitude of "1" and "2", see Figure 9a-c.To find the optimal ratio  of the whole metasurface in the wideband frequency, we need to overlap the −20 dB ranges of different frequency points.Here an example is given in Figure 9d that the overlap area is calculated by three frequency points, 12, 13 and 14 GHz.The ordinate range corresponding to the blue overlapping area is the value range of  to meet the 20 dB reduction, see shades of gray in Figure 9d.The ordinate corresponding to the maximum phase difference range is the optimal  value of these three frequency points, see the red point in Figure 9d.In the case that the overlap area of -20 dB exists, substitute the amplitude values of the two unit cells at each frequency point from 11.5-15.5GHz (step with 0.1 GHz), which is the frequency range that two unit cells have about 180° phase difference, and finally get the optimal  = 0.74.For simplicity, we set  = 0.75, that is to say, the ratio of the number of "1" to the number of "2" is 3:1.Thus, the simplest subarray is designed, as shown in Figure 10a.The whole metasurface is composed of 6 × 6 subarrays (12 × 12 unit cells), see Figure 10b.In the case that the overlap area of -20 dB exists, substitute the amplitude values of the two unit cells at each frequency point from 11.5-15.5GHz (step with 0.1 GHz), which is the frequency range that two unit cells have about 180 • phase difference, and finally get the optimal α = 0.74.For simplicity, we set α = 0.75, that is to say, the ratio of the number of "1" to the number of "2" is 3:1.Thus, the simplest subarray is designed, as shown in Figure 10a.The whole metasurface is composed of 6 × 6 subarrays (12 × 12 unit cells), see Figure 10b.
In the case that the overlap area of -20 dB exists, substitute the amplitude values of the two unit cells at each frequency point from 11.5-15.5GHz (step with 0.1 GHz), which is the frequency range that two unit cells have about 180° phase difference, and finally get the optimal  = 0.74.For simplicity, we set  = 0.75, that is to say, the ratio of the number of "1" to the number of "2" is 3:1.Thus, the simplest subarray is designed, as shown in Figure 10a.The whole metasurface is composed of 6 × 6 subarrays × 12 unit cells), see Figure 10b.To characterize the further RCS reduction behavior of the metasurface, full-wave simulations are performed in CST Microwave Studio.For sharp comparison, the equalsize PEC plane is also simulated.The RCS curves of the designed metasurface and PEC plane are shown in Figure 11a.It is easy to see the RCS of the metasurface is much lower than the equal-size PEC plane from 11-17 GHz.The RCS reduction curve is shown in To characterize the further RCS reduction behavior of the metasurface, full-wave simulations are performed in CST Microwave Studio.For sharp comparison, the equal-size PEC plane is also simulated.The RCS curves of the designed metasurface and PEC plane are shown in Figure 11a.It is easy to see the RCS of the metasurface is much lower than the equal-size PEC plane from 11-17 GHz.The RCS reduction curve is shown in Figure 11b to further illustrate the reduction ability of the designed metasurface.From Figure 11b, we can see that the above 20 dB reduction bandwidth is from 11.5-16.4GHz, which is much wider than the surface composed only by initial unit cell "1" (10.9-12.2GHz).Moreover, the 10 dB bandwidth is from 10.2-18 GHz, which is almost the same compared with the initial "1" (9.9-17.5 GHz).Not only the 20 dB bandwidth but also the 25 dB and even 30 dB are widened, see from Figure 11b.Without destroying the 10 dB bandwidth, the designed metasurface can get a broadband 20 dB RCS reduction.Due to the fourfold rotation symmetry feature of two unit cells, the metasurface is polarization insensitive.For verification, different polarizations (TE, TM) of incident wave are simulated in CST and the results are shown in Figure 11a.Obviously, it can be observed that the curves of different polarizations are overlapping.Thus, for the sake of simplicity, only the simulation results of TE polarization mode are given later.To illustrate that the further lower RCS is caused by two mechanisms: energy domain absorption and spatial domain diffusion, the 3D far-field patterns are shown in Figure 12.The patterns of the surface composed only by "1" (only has absorption but no cancellation diffusion, Figure 12b), patterns of traditional chessboard cancellation metasurface (only has cancellation diffusion but no absorption, composed by two 180°-phase-difference unit cells with total reflection, Figure 12c) and patterns of PEC plane (Figure 12a) are shown for comparison.From Figure 12a,b,d, we can easily see the designed metasurface would diffuse the reflected wave in the spatial domain, which inherits the feature of phase cancellation metasurface.Additionally, the designed metasurface would reduce the amplitude of reflected wave in the energy domain, which inherits the absorption feature of "1" and "2", see Figure 12a,c,d.These are strong evidence for the dual-mechanism work char- To illustrate that the further lower RCS is caused by two mechanisms: energy domain absorption and spatial domain diffusion, the 3D far-field patterns are shown in Figure 12.The patterns of the surface composed only by "1" (only has absorption but no cancellation diffusion, Figure 12b), patterns of traditional chessboard cancellation metasurface (only has cancellation diffusion but no absorption, composed by two 180 • -phase-difference unit cells with total reflection, Figure 12c) and patterns of PEC plane (Figure 12a) are shown for comparison.From Figure 12a,b,d, we can easily see the designed metasurface would diffuse the reflected wave in the spatial domain, which inherits the feature of phase cancellation metasurface.Additionally, the designed metasurface would reduce the amplitude of reflected wave in the energy domain, which inherits the absorption feature of "1" and "2", see Figure 12a,c,d.These are strong evidence for the dual-mechanism work characteristics.When the ratio  = 0.75, there is more than one arrangement of the whole metasurface.As shown in Figure 13, the subarray can also be 2 × 2, 3 × 3, 6 × 6, or random, etc.According to Equation (2), monostatic RCS at normal incidence is only related to the ratio between "1" and "2", regardless of the specific layout.For verification, the simulated results of different layouts but the same ratio are shown in Figure 13.As we can see from Figure 13f, the value of RCS reduction is almost the same.Only the 3-D far-field patterns are different, see Figure 13a-e.Once the ratio is fixed, the monostatic RCS at normal incidence is determined, and the different layout will only affect the scattering patterns.
The scattering patterns have no influence on monostatic RCS reduction in these cases but can affect the bi-static RCS.Thus, we choose the metasurface with 1 × 1 subarrays as the final design after comprehensively comparing the performance of the above layouts.This layout has the simplest topology and satisfactory mono/bi-static RCS reduction.It should be pointed out that the current layout is caused by the different absorption abilities of "1" and "2".For different cases, the dual-mechanism layout will vary with the different ratios .When the ratio α = 0.75, there is more than one arrangement of the whole metasurface.As shown in Figure 13, the subarray can also be 2 × 2, 3 × 3, 6 × 6, or random, etc.According to Equation (2), monostatic RCS at normal incidence is only related to the ratio between "1" and "2", regardless of the specific layout.For verification, the simulated results of different layouts but the same ratio are shown in Figure 13.As we can see from Figure 13f, the value of RCS reduction is almost the same.Only the 3-D far-field patterns are different, see Figure 13a-e.Once the ratio is fixed, the monostatic RCS at normal incidence is determined, and the different layout will only affect the scattering patterns.For illustrating the operating characteristics of oblique incidence, the simulated results of metasurface with 1 × 1 subarrays are shown in Figure 14 when illuminated by TE and TM polarized waves at 15°, 30°, 45° angles of incidence.When the incident angle is narrow, the 20 dB reduction remains obvious.As the incidence angle increases, the reduction effect becomes worse.However, the TE mode is always better than the TM mode under the same incident angle, which is consistent with the reference [37].The scattering patterns have no influence on monostatic RCS reduction in these cases but can affect the bi-static RCS.Thus, we choose the metasurface with 1 × 1 subarrays as the final design after comprehensively comparing the performance of the above layouts.This layout has the simplest topology and satisfactory mono/bi-static RCS reduction.It should be pointed out that the current layout is caused by the different absorption abilities of "1" and "2".For different cases, the dual-mechanism layout will vary with the different ratios α.
For illustrating the operating characteristics of oblique incidence, the simulated results of metasurface with 1 × 1 subarrays are shown in Figure 14 when illuminated by TE and TM polarized waves at 15 • , 30 • , 45 • angles of incidence.When the incident angle is narrow, the 20 dB reduction remains obvious.As the incidence angle increases, the reduction effect becomes worse.However, the TE mode is always better than the TM mode under the same incident angle, which is consistent with the reference [37].For illustrating the operating characteristics of oblique incidence, the simula sults of metasurface with 1 × 1 subarrays are shown in Figure 14 when illuminated and TM polarized waves at 15°, 30°, 45° angles of incidence.When the incident a narrow, the 20 dB reduction remains obvious.As the incidence angle increases, the tion effect becomes worse.However, the TE mode is always better than the TM under the same incident angle, which is consistent with the reference [37].

Metasurface Fabrication and Measurement
For experimental verification, a sample of the dual-mechanism absorptive m face is fabricated and experimentally characterized, see the photograph shown in 15a.According to the standard of the test system, a sample with the size of 300 mm mm is fabricated.In this paper, we use CAPITON-900 as the resistive film metasurface.CAPITON-900′s is a carbon polymer thick film ink.The two carbo

Metasurface Fabrication and Measurement
For experimental verification, a sample of the dual-mechanism absorptive metasurface is fabricated and experimentally characterized, see the photograph shown in Figure 15a.According to the standard of the test system, a sample with the size of 300 mm × 300 mm is fabricated.In this paper, we use CAPITON-900 as the resistive film in the metasurface.CAPITON-900 s is a carbon polymer thick film ink.The two carbon inks with different sheet resistance values are printed on one side of the F4BM-1 substrate by overprint technology, see Figure 15c.The RCS measurement is conducted in an anechoic chamber to avoid the disturbance of the background environment, see Figure 15b.Meanwhile, The TRL calibration technique and the time-domain gating feature of the network analyzer are used to eliminate the effects of undesirable multiple reflections.
In addition, there is a roughly 0.3 mm position error between two prints, which is the systematic error of the overprint system we used.Thus, we consider the overprint deviation during sample manufacturing and simulate the related results, see Figure 16.Suppose the "2" will be printed with 0.5 mm deviation, shown as the yellow shaded area, which is larger than 0.3 mm overprint error, and the simulated results are shown in Figure 16.We can see the RCS curves almost coincide with each other.Therefore, the errors within a certain range of fabrication cannot influence the performance of the metasurface sample.Meanwhile, it is powerful evidence that our design method has a degree of robustness.
with different sheet resistance values are printed on one side of the F4BM-1 substrate by overprint technology, see Figure 15c.The RCS measurement is conducted in an anechoic chamber to avoid the disturbance of the background environment, see Figure 15b.Meanwhile, The TRL calibration technique and the time-domain gating feature of the network analyzer are used to eliminate the effects of undesirable multiple reflections.In addition, there is a roughly 0.3 mm position error between two prints, which is the systematic error of the overprint system we used.Thus, we consider the overprint deviation during sample manufacturing and simulate the related results, see Figure 16.Suppose the "2" will be printed with 0.5 mm deviation, shown as the yellow shaded area, which is larger than 0.3 mm overprint error, and the simulated results are shown in Figure 16.We can see the RCS curves almost coincide with each other.Therefore, the errors within a certain range of fabrication cannot influence the performance of the metasurface sample.Meanwhile, it is powerful evidence that our design method has a degree of robustness.During the monostatic RCS testing, the transmitter and the receiver must have a 14° angle in our test system, see Figure 17.For the sake of rigor, the 14° angle case is also simulated and compared with the measured 14° "monostatic RCS" results.As shown in Figure 16, the 14° simulated result curve has a similar trend with the 0° curve, only is lower from 12-16 GHz.Additionally, their 20 dB bandwidths are almost the same, from 11.7-16 GHz.The measured curve, similarly, has the approximately coincident 20 dB bandwidth compared with the 14° simulated curve.The 10 dB bandwidth of the 14° measured result is wider than the simulated one.The slight deviations at low frequency are probably attributable to the finite size effect and tolerances that are inherent in fabrications and measurements.At high frequency, the main radiation lobe of double-ridge horns antennas we used start to split.In that case, the incident wavefront is not the plane, which can generate the phase gradient on the metasurface and make the reflected energy deflected, causing the wider 10 dB bandwidth at high frequency.In addition, the differences in the numerical and experimental results are mainly caused by phase shift of metasurface During the monostatic RCS testing, the transmitter and the receiver must have a 14 • angle in our test system, see Figure 17.For the sake of rigor, the 14 • angle case is also simulated and compared with the measured 14 • "monostatic RCS" results.As shown in Figure 16, the 14 • simulated result curve has a similar trend with the 0 • curve, only is lower from 12-16 GHz.Additionally, their 20 dB bandwidths are almost the same, from 11.7-16 GHz.The measured curve, similarly, has the approximately coincident 20 dB bandwidth compared with the 14 • simulated curve.The 10 dB bandwidth of the 14 • measured result is wider than the simulated one.The slight deviations at low frequency are probably attributable to the finite size effect and tolerances that are inherent in fabrications and measurements.At high frequency, the main radiation lobe of double-ridge horns antennas we used start to split.In that case, the incident wavefront is not the plane, which can generate the phase gradient on the metasurface and make the reflected energy deflected, causing the wider 10 dB bandwidth at high frequency.In addition, the differences in the numerical and experimental results are mainly caused by phase shift of metasurface sample.Since the phase is more sensitive than amplitude, a small processing error can cause a large phase change.The error in resistance film value will bring about not only impedance mismatch, but phase shift.All of the above factors along with mentioned manuscript cause the differences in the numerical and experimental results.Otherwise, we have also tested the oblique incidence cases utilizing the arch, see Figure 18.Due to limited experimental conditions, we only measure oblique incidence in TE mode.Inevitably, the 20 dB bandwidth is deteriorating with the oblique angles increasing, but the RCS reduction still keeps the satisfactory level.Additionally, the measured results are in good agreement with the simulation results in Figure 14.

Conclusions
In summary, a dual-mechanism absorptive metasurface is proposed and investigated.The design method can further reduce the RCS in the energy domain by absorbing Otherwise, we have also tested the oblique incidence cases utilizing the arch, see Figure 18.Due to limited experimental conditions, we only measure oblique incidence in TE mode.Inevitably, the 20 dB bandwidth is deteriorating with the oblique angles increasing, but the RCS reduction still keeps the satisfactory level.Additionally, the measured results are in good agreement with the simulation results in Figure 14.Otherwise, we have also tested the oblique incidence cases utilizing the arch, see Figure 18.Due to limited experimental conditions, we only measure oblique incidence in TE mode.Inevitably, the 20 dB bandwidth is deteriorating with the oblique angles increasing, but the RCS reduction still keeps the satisfactory level.Additionally, the measured results are in good agreement with the simulation results in Figure 14.

Conclusions
In summary, a dual-mechanism absorptive metasurface is proposed and investigated.The design method can further reduce the RCS in the energy domain by absorbing

Figure 1 .
Figure 1.The dual-mechanism absorptive meta surface of M × N.

Figure 4 .
Figure 4.The (a) structure and (b) equivalent transmission line model of the unit cell.

Figure 4 .
Figure 4.The (a) structure and (b) equivalent transmission line model of the unit cell.

Figure 6 .
Figure 6.Schematic of the two unit cells (a) unit cell "1" and (c) unit cell "2".The amplitude and phase curves of reflectivity of the two unit cells, (b) unit cell "1" and (d) unit cell "2".

17 Figure 7 .
Figure 7.The impedance analysis of two square resistive film rings.(a) Sheet resistance = 10 Ω/sq, parameter a varies from 4.5 to 6 with 0.5 step.(b) Parameter a = 5 mm, sheet resistance varies from 5 to 25 mm with 5 step.

Figure 7 .
Figure 7.The impedance analysis of two square resistive film rings.(a) Sheet resistance = 10 Ω/sq, parameter a varies from 4.5 to 6 with 0.5 step.(b) Parameter a = 5 mm, sheet resistance varies from 5 to 25 mm with 5 step.

5
to 25 mm with 5 step.

Figure 9 .
Figure 9.The 20 dB reduction area of (a) 12 GHz (b) 13 GHz and (c) 14 GHz.(d) The overlap area of these three frequency points.

Figure 9 .
Figure 9.The 20 dB reduction area of (a) 12 GHz (b) 13 GHz and (c) 14 GHz.(d) The overlap area of these three frequency points.

Figure 10 .
Figure 10.(a) The simplest subarray of ratio  = 0.75.(b) The schematic diagram of the whole metasurface.
Figure11bto further illustrate the reduction ability of the designed metasurface.From Figure11b, we can see that the above 20 dB reduction bandwidth is from 11.5-16.4GHz, which is much wider than the surface composed only by initial unit cell "1" (10.9-12.2GHz).Moreover, the 10 dB bandwidth is from 10.2-18 GHz, which is almost the same compared with the initial "1" (9.9-17.5 GHz).Not only the 20 dB bandwidth but also the 25 dB and even 30 dB are widened, see from Figure11b.Without destroying the 10 dB bandwidth, the designed metasurface can get a broadband 20 dB RCS reduction.Due to the fourfold rotation symmetry feature of two unit cells, the metasurface is polarization insensitive.For verification, different polarizations (TE, TM) of incident wave are simulated in CST and the results are shown in Figure11a.Obviously, it can be observed that the curves of different polarizations are overlapping.Thus, for the sake of simplicity, only the simulation results of TE polarization mode are given later.

Figure 11 .
Figure 11.(a) Simulated RCS curves of designed metasurface with TE and TM polarization, and equal-size PEC plane.(b) Simulated RCS reduction curves of designed metasurface and surface composed only by unit cell "1".

Figure 11 .
Figure 11.(a) Simulated RCS curves of designed metasurface with TE and TM polarization, and equal-size PEC plane.(b) Simulated RCS reduction curves of designed metasurface and surface composed only by unit cell "1".

Crystals 2022 , 17 Figure 13 .
Figure 13.(a-e) The different arrangements of the whole metasurface with the same ratio  and their 3D far-field patterns.(f) The RCS reduction curves of these metasurface arrangements.

Figure 13 .
Figure 13.(a-e) The different arrangements of the whole metasurface with the same ratio α and their 3D far-field patterns.(f) The RCS reduction curves of these metasurface arrangements.

Figure 13 .
Figure 13.(a-e) The different arrangements of the whole metasurface with the same ratio their 3D far-field patterns.(f) The RCS reduction curves of these metasurface arrangements

Figure 14 .
Figure 14.The operating characteristics when illuminated by TE and TM polarized wave 30°, 45° angles of incidence.

Figure 14 .
Figure 14.The operating characteristics when illuminated by TE and TM polarized waves at 15 • , 30 • , 45 • angles of incidence.

Figure 15 .
Figure 15.(a) The sample of dual-mechanism absorptive metasurface.(b) The testing platform of the anechoic chamber.(c) The production process of the overprinted metasurface.

Figure 15 .
Figure 15.(a) The sample of dual-mechanism absorptive metasurface.(b) The testing platform of the anechoic chamber.(c) The production process of the overprinted metasurface.Crystals 2022, 12, x FOR PEER REVIEW 14 of 17

Crystals 2022 , 17 Figure 17 .
Figure 17.The comparison of simulated and measured results.

Figure 18 .
Figure 18.The measured results of oblique incidence (TE mode).

Figure 17 .
Figure 17.The comparison of simulated and measured results.

Crystals 2022 , 17 Figure 17 .
Figure 17.The comparison of simulated and measured results.

Figure 18 .
Figure 18.The measured results of oblique incidence (TE mode).