An ultra-broadband terahertz absorber at high terahertz frequency

Broadband absorber at high terahertz frequency is highly required for applications in imaging, detecting, electromagnetic stealth. Although intensive investigations of the broadband absorber have been taken, the challenges still exist both in design and fabrication of an ultra-broadband absorber at high frequency. Herein, a three-layered structure metamaterial has been designed to realize an ultra-broadband terahertz absorber covering 3.94–9.98 THz (6.04 THZ) at the absorption above 80%, and the absorption bandwidth can be increased to 7.64THz (2.34–9.98 THz) by rotating the absorber. Upon on the simulations, the proposed absorber exhibits insensitive to the TM and TE polarization, it means the absorption effect is almost consistent in different polarization modes, the proposed absorber is significant in the practical application. Nevertheless, the absorption bandwidth reduces a little bit to 5.05 THz (absorption > 60%) as β increases to 45°.


Introduction
Terahertz wave is a kind of electromagnetic wave with a frequency range of 0.1-10 THz (Siegel 2002;Yen et al. 2004;Song et al. 2018), which has many unique advantages of strong penetrability, high resolution and low photon energy (Jiang et al. 2018;Pitchappa et al. 2016;Song et al. 2020). It has great potential in imaging, nondestructive detection, spectrum analysis (Chen et al. 2020a(Chen et al. , 2020bAbramov et al. 2021). The discovery of electromagnetic metamaterials has provided a novel approach for preparation of terahertz absorbing material. Since the first experimental demonstration of a near-perfect microwave metamaterial absorber has been conducted (Landy et al. 2008), the metamaterial terahertz absorber has received extensive attentions and researches (Wu et al. 2021;Biabanifard et al. 2020;Bai et al. 2018). However, most attempts are focused on the single absorption frequencies or narrow-band absorption of absorber (Ma et al. 2011;Xu et al. 2020;Huang et al. 2020;Verma et al. 2020;Wang et al. 2020). Actually, the broadband absorber can absorb more electromagnetic waves, it is more important in the fields, such as electromagnetic stealth, terahertz imaging and sensing. For electromagnetic stealth, the broadband absorber can absorb more terahertz wave at different frequency, it means less electromagnetic waves were reflected, finally less reflected electromagnetic waves could be accepted by radar, it's good for evading the search of radar. For terahertz imaging, the reflected THZ waves would be converted to image signal, and broadband absorber can absorb more reflected THZ waves at different frequency, it's good of improving the image accuracy. For terahertz sensing, broadband absorber can absorb more THZ waves which contains more information, it's good of improving the performance of THZ sensor.
In the last decade, many attempts of preparation of broadband absorbers have been made, however, the absorption bands are almost in the low terahertz frequency (0.1-5 THz) (Huang et al. 2018;Zhu et al. 2019;Wang et al. 2014). There is still challenges for designing ultra-broadband absorber at the high terahertz frequency (Liu et al. 2019;Daraei et al. 2020;Mohamadreza et al. 2020).
In this letter, an ultra-broadband terahertz absorber was proposed, which is a metamaterial with a three-layered structure. The absorption of THZ absorber can be above 80% in the high frequency of 3.94-9.98 THz, and almost above 60% at the entire terahertz frequency band. The special designed ten nested aluminum rings with an opening structure will response for such an ultra-broadband absorption. Significantly, the proposed absorber exhibits insensitive to the TM and TE polarization no matter what kinds of incident angles are applied.

Structure and design
The proposed broadband absorber is a kind of periodic metamaterials, and the cell is schematically illustrated in Fig. 1. It contains three layers: a metallic bottom layer, another metallic top layer and a dielectric middle layer. The proposed metamaterial absorber is represented as "NROS" (nested rings with an opening structure) absorber. Specifically, the top layer is consisted of ten nested rings with an opening structure. The diameter of the innermost and outmost circle is represented as x 1 and x 2 , while the width of each circle and inter distance of adjacent circles is recorded as x 3 and x 4 , respectively. The width of opening structure is d, and the period of the unit cell both in x and y directions is p. The thickness of each layer from top to down is defined as l 1 , l 2 and l 3 , respectively. Based on the systematic simulations, the absorption can achieve maximum value when both the relative permittivity (ε 1 ) and permeability (μ 1 ) of metallic layer are about 1. Meanwhile, the ε 2 and μ 2 of middle dielectric layer should be about 3.5 and 1, respectively. Therefore, the aluminum metal is chosen for both of top and bottom layer, while polyethylene is selected to the middle dielectric layer with ε 2 = 3.5 and μ 2 = 1, respectively. All the design parameters are shown in Table 1.
The simulations of the periodic structure are carried out by employing finite difference time domain method based on "ANSYS Electronics 2017". Absorption can be calculated by According to t he matrix transmission principle, the transmission matrix of electromagnetic wave in metamaterials can be calculated as: In the Eq. 1, n is refractive index of uniform tablet, k is wave vector of incident wave, h is thickness of uniform tablet, z is equivalent impedance of medium plate. The structure size of metamaterials is much smaller than the working wavelength, the metamaterials can be consider as homogeneous materials. S 11 and S 21 can be calculated as: Equation 4 shows the relation of refractive index (n), impedance (z), permittivity (ε), and permeability (μ).
According to impedance matching principle, when equivalent impedance (z) is closer to the impedance of metamaterial interface with air, the reflectivity is lower. As Eq. 4 shown, we can change the impedance (z) by controlling permittivity (ε). For the metamaterial absorber, the permittivity (ε) can be selected by controlling the size and geometry of periodic structure. Different from the "CNR" (closed and nested rings) absorber, NROS absorber has an opening structure in the nested rings which makes a better absorption broadband. Figure 2(a-d) are the electric fields of CNR and NROS metamaterial absorber at their resonance frequency, the electric fields are mainly distributed in the horizontal metal line of the closed nested metal ring in the resonant state, it means the horizontal metal line of the closed nested metal ring is equivalent to inductance (L), the vertical metal line is equivalent to capacitance (C) (see Fig. 2(f)). Figure 2 (f) is equivalent circuit of CNR and NROS metamaterial absorber, R in the circuit is the equivalent resistance of energy loss. L and C are the equivalent inductance and capacitance of the metal patch structure. Whether NROS absorber or CNR absorber, each nested ring means a LC resonant circuit, there are 10 parallel LC resonant circuits in the equivalent circuit of CNR and NROS metamaterial absorber. When the metamaterial absorber is in perfect absorption, each LC resonance circuit can make an absorption curve, and the absorption curve is equivalent to the resonance curve, so the frequency of absorption peak is equal to the resonant frequency [see Eq. 9].
T is the oscillation period of LC resonance circuits, it also can be the bandwidth of absorption. There are ten LC resonance circuits, it means there would be ten absorption curves, making ten absorption curves closer by controlling L and C, finally get an broadband absorption[see Fig. 2 Different from the CNR metamaterial absorber, there is an opening structure in the nested rings of NROS metamaterial absorber, it means in each parallel LC resonant circuit, the inductance (Ln, n = 1, 2… 10) has been reduced [see Fig. 2(f)], while the equivalent circuit is a parallel circuit, the equivalent inductance (L) would be increased. The equivalent inductance (L) and equivalent capacitance (C) could be calculated as Eq. 10-11: The equivalent inductance (L) has been reduced [see Fig. 2(f)], while the period of oscillation T is in direct proportion to inductance (L) [see Eq. 9], T of NROS absorber is longer than CNR absorber, it leads the broadband of NROS absorber is longer than it of CNR absorber [see Fig. 2 As Fig. 2(a-d) shown, the electric field of the CNR metamaterial absorber at the resonance frequency (2.5 THZ and 4.3 THZ) is stronger than which of NROS metamaterial absorber (4.3 THZ and 7.2 THZ). In LC resonance circuit, as the electric field intensity increases, the oscillation period will decrease. As Fig. 2(a-d) shown, the electric fields of the CNR metamaterial absorber at its resonance frequency are stronger than which of NROS metamaterial absorber, it means the oscillation period of the NROS metamaterial absorber is longer than which of CNR metamaterial absorber, it also proved that the broadband of NROS absorber is longer than which of CNR absorber. (10) The electric field distribution of the CNR metamaterial absorber at the resonance frequency 2.5 THZ. b The current field distribution of the CNR metamaterial absorber at the resonance frequency 4.3 THZ. c The current field distribution of the NROS metamaterial absorber at the resonance frequency 4.3 THZ. d The current field distribution of the NROS metamaterial absorber at the resonance frequency 7.2 THZ. e Absorption curves with optimized structural parameters of NROS absorber and CNR absorber. f Equivalent circuit diagram of CNR and NROS metamaterial absorber

Results and discussion
Following, the angle between the open structure and the positive direction of the magnetic field was defined as α, which represents the angle between the positive direction of the magnetic field and the positive direction of central axis of opening structure [ Fig. 3(a)]. Both parameters of α and d are proven to have highly effect on the absorption feature of the absorber. As d increases from 1 to 6 μm, the third absorption peak (arrow in Fig. 3(b)) shifts to the high frequency gradually, and resulting in decreasing of bandwidth in the 4-8 THz. Otherwise, the absorption is below 80% in the 4-8 THz as the d is set as 1 μm, while it is below 80% in the 9-10 THz if d is set as 6 μm [ Fig. 3(b)]. As a result, d is optimized as 5 μm in our design. Figure 3(c, d) plot absorption curves with various values of the α in the range of 0-90° at a step of 10°, respectively. The curve shape exhibits two obvious stages [ Fig. 3(c, d)]. As α increases from 0 to 40°, the absorption gradually reduces to below 80% in 5-8 THz, but still above 60% [ Fig. 3(c)]. In contrast, as α increases from 50° to 90°, the absorption broadband shifts to 2-5 THz, but the absorption increases in 2.5-3.5 THz [ Fig. 3(d)]. In this terahertz band, the smallest absorption is still above 70% (3.23 THz, α = 50°), and the absorption broadband is calculated to be 2.8 THz (2.32-5.12 THz) at absorption above 70%. Finally, two isolated absorption broadband can be achieved as α set as 50°-90°, respectively, and the absorption bandwidth can be increased to 7.64THz (2.34-9.98 THz) by controlling the α. Fig. 3 a Schematic diagram of the NROS absorbing structure. b Absorption curves of NROS absorber with various d. Absorption curves of NROS absorber as α is set from 0° to 40° (c), and from 50° to 90° (d). All the absorption results are simulated by the TE polarization mode Furthermore, the structural parameters are systematically investigated by simulations. As the thickness of top layer (l 1 ) increases from 0.02 to 0.10 μm, the absorption curves keep the similar shape but a little decrease of intensity in the 4-10 THz [ Fig. 4(a)]. It indicates a little influence of the thickness of top metals. Therefore, the thickness of top metallic layer is set to 0.02 μm in consideration of the convenience of fabrication process. For the middle dielectric layer, the absorption curves exhibit a big change at the different thicknesses (l 2 ) [ Fig. 4(b)]. In brief, the absorption is below 80% in 5-7 THz as l 2 < 7 μm, while a smaller value is found in 9-10 THz as l 2 > 7 μm [ Fig. 4(b)]. Therefore, l 2 is optimized to be 7 μm in our design.
For the width of each ring (x 3 ), it has a significant influence on the third absorption peak. As x 3 increases from 0.2 to 1.0 μm, the peak shifts to the lower frequency, which can extend absorption broadband efficiently [ Fig. 4(c)]. However, there is part of absorption broadband between the third and fourth peak below 80% [ Fig. 4(c)] as x 3 > 0.4 μm. Meanwhile, there is still a part (7.5-8.5 THz) with absorption below 80% as x 3 = 0.2 μm. Above all, x 3 is optimized to be 0.4 μm herein. In addition, the distance between the adjacent rings (x 4 ) shows a remarkable influence on the absorption peak position and intensity. Figure 4(d) shows a down shift of the third peak as increase of x 4 , leading to extend of absorption regime. Meanwhile, the absorption intensity between the third and fourth peak gradually reduces, even below 80%. For additional reason, there is much difficult and highcost in fabrication process if narrows the adjacent distance further. Above all, the x 4 was set to be 0.2 μm finally.
For the number of rings, the absorption curves [ Fig. 4(e)] show a similar behavior as those of parameter of x 4 . The third peak shifts to the low frequency as increase of number from 7 to 20, thereby extending the absorption broadband efficiently. However, the absorption between the third and fourth peaks will be gradually decreased below 80% if the number is more than 10. To realize the ultra-broadband absorption and reduce the complex of fabrication process, the number of metallic rings is set to 10 in our design.
After optimization analysis, the absorption bandwidth is kept about 6.04 THz when the diameter of the innermost (x 1 ) is about 14.4 μm. Therefore, the x 1 is set to be 14.4 μm in our design. According to the optimized values of x 3 , x 4 and number of rings, the diameter of the outmost circle (x 2 ) is calculated to be 26 μm. In addition, the absorption bandwidth has been kept about 6 THz when the periodicity (p) is around 30.
After optimization of structural parameters of absorber, a broadband of 6.04 THz at high frequency of 3.94-9.98 THz is achieved at incident angle of 0° [ Fig. 4(f) black line]. Notably, all the absorption in this regime is above 80%, with the highest and average absorption of 99.26% and 89.43%, respectively. In contrast, a controlled absorber is developed which has the same structural parameters as the NROS absorber, but with 10 closed nested rings. It is named as "CNR". The red line in Fig. 4(f) exhibits its absorption curve. The absorption bandwidth is obtained cover two frequency regimes of 2.41 THz (from 2.23-4.64 THz) and 1.89 THz (from 7.52-9.41 THz). Nevertheless, the absorption bandwidth is narrower than that of "NROS" absorber, and the absorption is lower than 80% in most of absorption bandwidth too. Figure 4(g) plots absorption curves in the range of 1-10 THz at TM and TE polarized incidence. Two curves are near overlapped, indicating the proposed absorber is insensitive to the TM and TE polarization. The absorption band covers from 3.94 to 9.98 THz at the absorption above 80%. This broadband of 6.04THz at high frequency is rarely achieved before both in the simulation and experiment. Figure 4(h) depicts the absorption curves with various incident angles of 0-60° at a step of 15°. Generally, the absorption curves show the similar shape with gradually decrease of intensity [ Fig. 4(h)]. The largest bandwidth is calculated to be 6.23 THz at absorption above 60% achieved at β = 0°, while it reduces to its smallest value of 5.05 THz at β = 45°. The drop ratio is only 19%, suggesting the proposed NROS absorber is a terahertz polarization insensitive absorber.
To analyse the corresponding background physical mechanism to each simulation result, the electric fields of proposed absorber was shown in Fig. 5-8. Figure 5 depicts the simulated 2D patterns of the electric field distribution of the NROS and CNR absorber at different absorption peaks. Figure 5(a-c) depict the distribution of electric field at absorption peaks of 4.3, 7.2 and 8.8 THz, respectively, when the angle α is 0°. The position of the maximum intensity moves from outer to inner ring, and back to outer one again. Meanwhile, the maximum electric field intensity gradually decreases [ Fig. 5(a-c)], and the stronger field are mainly distributed near the opening structure [ Fig. 5(a-c)].
Because the horizontal metal line of the closed nested metal ring is equivalent to inductance (L n ), the vertical metal line is equivalent to capacitance (C n ). The openning structure of NROS absorber will make the inductance (Ln) decreased, resulting in a decrease in equivalent inductance (L), finally the f will increase, the absorption peak will shift right (Fig. 2e). In addition, the distribution exhibits symmetrical about the central axis of the opening structure at three absorption peaks of 2.5, 4 and 5 THz, and the stronger field are mainly distributed near the opening structure too [ Fig. 5(d-f)]. As a controlled experiment, the same simulations have been carried out on the CNR absorbers to investigate the distribution of electric field. Figure 5(g-i) exhibit the field moves from outer ring to inner one, and then back to outer one again. Meanwhile, the intensity of the electric field is gradually reduced.
As d increases, the electric field intensity of proposed absorber at 3.7THz has been increased ( Fig. 6a-c).
Stronger electric field intensity leads higher absorption, it agrees with the simulation results (Fig. 3b). And according to Eqs. 9 and 10, as d increases, the inductance(L n ) (n = 1,2,3…10) will decrease which leads a increase in equivalent inductance (L), finally the resonant frequency f will be decreased, and the third absorption peak will move to the high frequency (Fig. 3b). As Fig. 6d-f shown, the electric field intensity will be decreased with α increased, which makes a decrease in absorption (Fig. 3c,d). Figure 7a-c illustrate the electric field distribution of different width of metallic ring(x 3 ), absorption is positively correlated with the electric field intensity, when x 3 is 0.4 μm the electric field intensity is strongest (Fig. 7b), it proves the red absorption curve is higher than other curves at 3.7THz (Fig. 4c), the width of metallic ring(x 3 ) should be set as 0.4 μm. In the same way, the distance of adjacent rings (x 4 ) and number of metallic rings should be set as 0.2 μm and 10( Fig. 7d-f). According to Eqs. 9 and 10, the resonant frequency (f) is negatively correlated with the equivalent inductance(L). As x 3 increases, the equivalent inductance(L) will increase with L n increased, which leads a decrease in resonant frequency (f), finally the third absorption peak moves to the left (Fig. 4c). As x 4 increases, the resistance (R n ) will increase which leads the equivalent inductance (L) increased, resulting in a decrease in the resonant frequency (f), finally the third absorption peak moves to the left (Fig. 4d). Fig. 4 a Thickness of top metallic layer (l 1 ). b Thickness of middle dielectric layer (l 2 ). c Width of metallic ring (x 3 ). d Distance of adjacent rings (x 4 ). e Number of metallic rings. f Absorption curves with optimized structural parameters of NROS absorber and its controlled CNR absorber. g Absorption curves with different polarized incidence (TE and TM). h Absorption curves at polarization angles φ = 0-60° at a step of 15°. All the absorption results are simulated by the TE polarization mode except Fig. 4 (g) As Fig. 7g-i shown, when the number of metallic rings is 10, the electric field intensity is strongest (Fig. 7h), it proves the blue absorption curve is higher than other curves at 3.7THz (Fig. 4e). And as the number of metallic rings increase, number of L n and C n will increase, which leads equivalent inductance (L) decreased and equivalent inductance (C) increased(see Eq. 10-11). Because the horizontal metal line of the closed nested metal ring is equivalent to inductance (L n ), the vertical metal line is equivalent to capacitance (C n ), and the horizontal metal line is shorter than vertical metal line, the increments of C is more The absorption peaks at 2.5, 4.0, 4.8 THz at α = 70°, respectively. G-i The absorption peaks at 2.4, 4.3, 7.7 THz of CNR absorber, respectively Fig. 6 2D pattern of distribution of electric field of proposed absorber at 3.7THz. a-c The d is 1 μm, 3 μm and 5 μm, respectively. d-f The α is 30°, 60° and 90°,respectively than the decrements of L. In conclusion, as the number of metallic rings increases, LC has been increased, the resonant frequency (f) will decrease, finally the third absorption peak moves to the left (Fig. 4e). Figure 8a-c illustrate the electric field distribution of different thickness of metallic ring(l 1 ) at 9.0THz, as l 1 increases the electric field intensity has been decreased, the  Distribution of electric field of proposed absorber at 9.0THz. a-c Thickness of metallic ring (l 1 ) is 0.02 μm, 0.04 μm and 0.06 μm, respectively. d-f Thickness of dielectric layer (l 2 ) is 0.1 μm, 0.2 μm, 0.3 μm,respectively absorption at 9.0THz can be decreased (Fig. 4a). When the thickness of l 2 is 7 μm, the electric field intensity of dielectric layer is strongest (Fig. 8e), the absorption at 9.0THz should be highest(green curve Fig. 4b).
Above all, three results can be drawn: 1) the strongest intensity of electric field mainly distributed near the opening structure; 2) the opening structure is response for the absorption peak in NROS absorber, which can also change the distribution of electric field compared to the CNR absorber; 3) the distribution of electric field is symmetrical about the central axis of opening structure. In the proposed structure, the periodic metallic layer is equivalent to an LC (inductance and capacitance) resonance circuit. In this circuit, the intensity of the electric field and the oscillation period (distance between two absorption peaks) are negatively correlated. As a result, the intensity of the electric field of CNR absorber is stronger than that of NROS one [see Fig. 5(g) and Fig. 5(a)], which means the oscillation period of CNR absorber is shorter than that of NROS one. With the reduced oscillation period, the distance between two closed absorption peaks will reduce, thereby the absorption broadband reduces accordingly. Table 2 summarized our findings with the previous results. Significantly, the proposed NROS absorber shows the largest broadband at high terahertz frequency simultaneously. Bandwidth (BW) and Relative absorption bandwidth (RBW) are two important evaluation indexes of the THZ absorber. Compared with the previous results, our designed broadband absorber can cover almost all the high terahertz frequency. This high frequency (equal to short wavelength) absorber can be applied in tiny targets nondestructive detecting and accurate imaging. Meanwhile, the ultra-broadband of 6.04 THz may be used for radar monitoring and terahertz communication.

Conclusion
In conclusion, we present an ultra-broadband terahertz absorber at high terahertz frequency. The absorption bandwidth is up to 6.04 THZ (3.94-9.98 THZ) at the absorption above 80%, while a near entire terahertz regime is achieved at the absorption above 60%. Meanwhile, the highest absorption is up to 99.26%, and the average value is still as high as 89.43%. In addition, the proposed absorber exhibits insensitive to the TM and TE polarization no matter what kinds of incident angles are applied. The opening structure is response for such a wide absorption broadband, and the absorption broadband can be regulated by adjusting the incident angle.