Five-Band Tunable Terahertz Metamaterial Absorber Using Two Sets of Different-Sized Graphene-Based Copper-Coin-like Resonators

: In this paper, a five-band metamaterial absorber with a tunable function in a terahertz band is proposed, which consists of a gold grounding layer, a polyimide dielectric layer, and a periodic patterned graphene layer. The patterned graphene layer is constructed from two sets of copper-coin-shaped structures of different sizes. The designed absorber achieves absorptions of 96.4%, 99.4%, 99.8%, 98.4%, and 99.9% at 4.62 THz, 7.29 THz, 7.70 THz, 8.19 THz, and 8.93 THz, respectively, with an average absorption intensity of 98.78%. The physical mechanism of this five-band absorber was explained by the impedance matching principle and electric field distribution. The absorption performance of the five-band absorber can be effectively tuned by changing the geometry of the patterned graphene array and the thickness of the dielectric layer. Given that the resonant frequency of the absorber varies in proportion to the Fermi level, by varying the Fermi level of the graphene hypersurface, we can achieve the continuous tuning of the absorption performance over a wide frequency range. The five-band absorber has a stable absorption performance over a wide incidence angle of 0–65 ◦ , and by combining the merits of high absorption, dynamic adjustability, and a large number of absorption peaks, the given absorber could have great potential for applications in nondestructive testing


Introduction
Terahertz (THz) is an electromagnetic wave with a frequency range from 0.1 THz to 10 THz [1].Due to the special position of the terahertz band between the microwave and infrared bands of the electromagnetic spectrum, terahertz bands share the properties of electrons and photons.Consequently, terahertz exist vast applications in imaging technology [2], spectral resolution [3], and communication technology [4].Metamaterials are artificially structured materials consisting of subwavelength arrays of periodically or non-periodically arranged resonant crystal cells [5], which are regarded as equivalent media, and the equivalent permittivity ε(ω) and equivalent permeability µ(ω) can be manipulated according to the actual needs to realize the absorption of electromagnetic waves of different frequencies [6].It has bizarre properties beyond natural materials, including a negative refractive index [7,8], ultra-thin perfect absorption [9,10], an inverse Doppler effect [11,12], and so on.In view of the physical properties of metamaterials, in line with the current design needs for the miniaturization and integration of electronic components, a variety of metamaterial functional devices applicable to different scenarios has been widely proposed, which are of great significance for the research of stealth technology, perfect absorption [13], plasma-induced transparency [14,15], and sensors [16,17].
Among the above functional devices, metamaterial absorbers have received much attention for their advantages of small sizes, light weight, and high absorption.Landy et al. realized perfect absorption in the microwave band for the first time based on the resonant structure of metal surfaces in 2008 [18].With the continuous development of micro and nanotechnology, perfect absorption has been gradually introduced into visible, infrared, and terahertz bands [19][20][21].Among them, perfect absorbers (PAs), which produce a strong response in the terahertz band, have become quite attractive in recent years [22,23].To expand the potential applications of PAs, there is an urgent need for PAs with tunable performance or multi-band perfect absorbers.According to the number of resonant absorption peaks, multi-band metamaterial absorber devices can be categorized as dual-band absorbers, triple-band absorbers, quadruple-band absorbers, or even five-band absorbers.For example, Wen et al. developed a dual-band polarization-insensitive electromagnetic wave absorber by exploiting the LC resonance of metamaterials to produce a strong absorption point [24].A three-band terahertz absorber utilizing three concentric square ring resonators was reported to achieve triple-frequency absorption with both high sensitivity and polarization insensitivity [25].Fu et al. realized a multiwavelength absorber based on a metallic subwavelength structure using hexagonal filled nanoring-nanodisc composite arrays to introduce an asymmetric structure and realize new absorption functions [26].However, the above proposed multi-band absorbers are usually based on metallic subwavelength structures, and once the original geometrical structure is determined, their optical properties are difficult to change, with the typical defect of a limited tuning range [27][28][29].
Graphene, as a two-dimensional nanomaterial, has excellent optical and electrical properties such as high carrier mobility, low loss, and an ultrafast optical response time [30][31][32].When a bias voltage is added to the surface of graphene, the Fermi level of graphene varies with the external voltage, which results in the continuous regulation of conductivity over a wide frequency range [33][34][35][36][37][38][39][40][41].Graphene has been widely used to construct tunable PAs, mainly including multi-band and broadband absorption.The realization of multi-band absorption performance is mainly classified into vertical stacked and coplanar types, which can bring about significant improvements in absorption performance by stacking multilayers of graphene, designing graphene gratings, and changing the periodic array pattern to increase the number of graphene resonators.For example, in 2019, Liu et al. proposed a dynamically and independently tunable absorber based on multilayered metal-graphene metamaterials, which can be arbitrarily tailored to dual-band or even multi-band absorption efficacy by stacking multiple metal-graphene layers [41][42][43][44][45][46][47][48][49].However, the vertical stacking type often brings more complicated structures, and various multilayer designs lead to thicker structure dimensions, which makes practical production difficult [50][51][52].Most of the absorbers also have limitations, such as a low number of absorption peaks, poor absorption quality, and insensitivity of polarization, so it is still feasible to explore a simple structure of a high-quality multi-band absorber based on a single-layer coplanar design.
In this paper, a tunable five-band terahertz absorber based on a graphene metasurface is proposed.It is a typical three-layer structure consisting of a bottom gold grounding layer, an intermediate polyimide dielectric layer, and a top patterned graphene layer.Five near-unity absorption could be realized, which could be characterized by polarization and angle sensitivity and wide-angle and flexible tunability.The results show that the peak absorbance reaches 96.4%, 99.4%, 99.8%, 98.4%, and 99.9% at 4.62 THz, 7.29 THz, 7.70 THz, 8.19 THz, and 8.93 THz, respectively.The resonance frequency of the absorber varies in direct proportion to the graphene's Fermi level.By adjusting the Fermi level of graphene, the absorber performance can be continuously adjusted within the working range, and the average absorption can be increased from 96.38% to 98.78%.By changing the angle of incidence, the absorber can be considered to have stable absorption performance within 0-65 • , characterized by wide-angle incidence.Due to the excellent properties of the structure, such as high absorbance, strong tunability, and wide-angle absorption, it has a broad application prospect in the fields of terahertz sensing, nondestructive testing, imaging, stealth, and detectors.

Materials and Methods
Based on the design needs of a multi-band absorber, we designed a five-band absorber by combining two different shapes of split-ring structures based on the patterns of circular and square graphene rings.The designed five-band PA structure is shown in Figure 1, where Figure 1a,b show the side view of the structure and the top view of the patterned layer, respectively.The absorber system has a classical sandwich structure, where the bottom metal layer is separated from the top patterned layer by a layer of polyimide with a relative dielectric constant of 2.89 at a thickness of t 2 = 4 µm, and the patterned layer consists of two sets of different-sized graphene-based copper-coin-like resonators.Gold was chosen as the all-reflective metallic ground plane, whose conductivity is 4.56 × 10 7 S/m with a thickness t 1 = 0.

Results and Discussion
The absorption curve of the five-band metamaterial absorber is given in Figure 2. The absorption rate can be expressed as A(ω) denotes the absorptivity, , are the scattering coefficients, the reflectivity is = | | , and the transmittance is = | | .In the terahertz range, the thickness of the gold is much larger than the skin depth.It can be approximated that = 0, so the absorptivity can be simplified to The reflection is highly dependent on the impedance, and when the impedance is adjusted to match perfectly with free space, the incident wave can be absorbed significantly.From the impedance matching principle, the absorption rate can be expressed as follows [60][61][62]: where Z and Z0 are the effective impedance of the absorber and free space, respectively, and the relative impedance of the absorber can be calculated by parametric inversion [63]: The reflection coefficient is minimized when the absorber equivalent impedance matches the free space impedance.From the given design parameters, perfect absorption up to 96.4%, 99.4%, 99.8%, 98.4%, and 99.9% is achieved at five discrete points (labeled as The surface dielectric constant of graphene can be expressed as follows [53,54]: where the thickness of graphene is 0.334 nm, the free space impedance is Z 0 = 377 ohm, k 0 denotes the wave vector number of free space electromagnetic wave, and σ gra denotes the surface conductivity of graphene.In the model, the surface conductivity of monolayer graphene can be described by Kubo's formula, which includes intraband and interband terms [55,56]: where e, T, k B , and h are the charge of the electron, the absolute temperature (T = 300 K), Boltzmann's constant, and the approximate Planck's constant, respectively, and ω, τ, and E f are the angular frequency, the electron-phonon relaxation time, and the Fermi level, respectively.At a room temperature of 300 K, hω ≪ 2E f within the terahertz frequency according to the Pauli exclusion principle.Intraband leaps in the terahertz region dominate compared to interband leaps.Therefore, the above equations can be simplified to a Drudelike form [57][58][59]: In the following study, the Fermi level is set as 1.2 eV, and the electron-phonon relaxation time is 1 ps.The absorption properties and physical mechanisms of this absorbing system are investigated in detail using the COMSOL Multiphysics 6.1 software based on the finite element method.The periodic boundary conditions are set in both x-axis and y-axis directions, and Floquet ports are assigned in the direction perpendicular to the graphene plane.In addition, the incident light source is a plane wave, incident in the x direction.

Results and Discussion
The absorption curve of the five-band metamaterial absorber is given in Figure 2. The absorption rate can be expressed as where A(ω) denotes the absorptivity, S 11 , S 21 are the scattering coefficients, the reflectivity is R(ω) = |S 11 | 2 , and the transmittance is T(ω) = |S 21 | 2 .In the terahertz range, the thickness of the gold is much larger than the skin depth.It can be approximated that T(ω) = 0, so the absorptivity can be simplified to The reflection is highly dependent on the impedance, and when the impedance is adjusted to match perfectly with free space, the incident wave can be absorbed significantly.From the impedance matching principle, the absorption rate can be expressed as follows [60][61][62]: where Z and Z 0 are the effective impedance of the absorber and free space, respectively, and the relative impedance of the absorber can be calculated by parametric inversion [63]:   The reflection coefficient is minimized when the absorber equivalent impedance matches the free space impedance.From the given design parameters, perfect absorption up to 96.4%, 99.4%, 99.8%, 98.4%, and 99.9% is achieved at five discrete points (labeled as f 1 , f 2 , f 3 , f 4 , and f 5 ) at 4.62 THz, 7.29 THz, 7.70 THz, 8.19 THz, and 8.93 THz, respectively.Figure 3 shows the relative impedance of the absorber at the operating frequencies.The relative impedance of the absorber is highly matched to free space at five absorption frequencies, with the real part close to 1 and the imaginary part close to 0.  At f3 = 7.70 THz, the electric field is mainly concentrated at the cleavage of the r square ring, where the graphene surface couples with the magnetic field to enhance absorption effect.The spacing between the upper two resonators is pulled in compare the lower two resonators, and there is a strong coupling effect, and the electromagn In order to further elaborate the absorption mechanism of five-band high absorption, the electric field distributions, E Z , and the current density distributions corresponding to the five absorption peaks of the metamaterial absorber are given in Figures 4 and 5.For f 1 = 4.62 THz, consider the square resonator as an inductor, L, and the opening as a capacitor, C. According to the principle of equivalent circuit, the induced charges are generated by an electromagnetic wave at an incidence converge at the slit to form an electric field.Due to the narrow spatial position at the square ring slit, the movement of the charge forms a displacement current.As shown in Figure 4c, at this location, the current forms a loop on the graphene surface, producing an LC resonance mode.We can find that the electric field is partially concentrated at the top and bottom of the lower side of the copper-coin-like resonator due to the excitation of the surface plasmon excitations.It can be assumed that the high absorption at f 1 = 4.62 THz is produced by the combined effect of two different sizes of copper-coin-like resonators.For the absorption peak at f 2 = 7.29 THz, the electric field is mainly concentrated above the left-side square ring resonance cavity splits and below the right-side square ring resonance cavity splits, while some electric field distributions of relatively low intensity also exist around the outer ring.The coupling between the surface electronic oscillations of graphene and the incident electromagnetic wave enhances the localized electric field at the cleavage.This localization of the electromagnetic field leads to the localization of the optical field, which significantly enhances the absorption of light.Figure 4b,e represent the distribution of the z-component of the electric field in the corresponding modes, both of which are in good agreement with the distribution of the electric field, corroborating the analysis.
peak at f5 = 8.93 THz, the absorption pattern is similar to the peak at f3 = 7.70 THz, and the electric field is concentrated at the slit of the square annular resonant cavity.The electromagnetic wave induces localized surface plasmon resonance at the above location, which significantly enhances the absorption effect.The surface currents in Figure 5e,f,i also correspond well with the electric field, and the currents are mainly concentrated at the slit or in the narrow region of the space between the two resonators.electric field is concentrated at the slit of the square annular resonant cavity.The electromagnetic wave induces localized surface plasmon resonance at the above location, which significantly enhances the absorption effect.The surface currents in Figure 5e,f,i also correspond well with the electric field, and the currents are mainly concentrated at the slit or in the narrow region of the space between the two resonators.At f 3 = 7.70 THz, the electric field is mainly concentrated at the cleavage of the right square ring, where the graphene surface couples with the magnetic field to enhance the absorption effect.The spacing between the upper two resonators is pulled in compared to the lower two resonators, and there is a strong coupling effect, and the electromagnetic wave induces a localized surface plasma resonance at this slit.At the same time, the uppermost edge of the outer ring is also accompanied by an electric field distribution, which is conducive to increased absorption.At f 4 = 8.19 THz, the electric field is mainly concentrated in the lower part of the square ring resonator cavity at the slit.Figure 5e shows the E z distribution at f 4 = 8.19 THz.For the lower side opening of the left square ring resonator, the positive charges are concentrated at the left end, and the negative charges are concentrated at the right end, forming a dipole resonance.Similarly, a dipole resonance is formed on the right-side square ring resonator, but with opposite polarity to the left side.For the peak at f 5 = 8.93 THz, the absorption pattern is similar to the peak at f 3 = 7.70 THz, and the electric field is concentrated at the slit of the square annular resonant cavity.The electromagnetic wave induces localized surface plasmon resonance at the above location, which significantly enhances the absorption effect.The surface currents in Figure 5e,f,i also correspond well with the electric field, and the currents are mainly concentrated at the slit or in the narrow region of the space between the two resonators.
Figure 6 illustrates the effect of the variation in the geometrical parameters of the device on the absorption effect.As shown in Figure 6a, when the edge length a increases, the frequency of the f 1 undergoes a slight redshift, and the f 2 -f 5 absorption intensities undergo significant changes.With the increase in a, the effective area of graphene increases, the distance between the two sub-resonant cavities above gradually draws closer, the mutual coupling interaction is enhanced, and the absorption capacity increases.However, after increasing to a certain value, an interaction that is too strong causes the absorption rate to decrease.Figure 6b shows the absorption spectra under different parameters above the hollow square sheet b.With the change in b, the overall frequency range is slightly redshifted and then significantly blue-shifted, and there is a tendency to cleave into two peaks at f 1 , and the overall absorption quality increases and then weakens.With an increasing b value, the internal arrangement of graphene becomes denser, and the increase in surface plasmon exciton polarizers leads to an increase in absorbance.However, after reaching the peak, the resonance peaks excited between different regions interfere with each other, leading to a decrease in the absorption rate.Similarly, the change in the absorption pattern of parameter c is similar to that of b.The effect of the parameter r of the copper-coin-shaped resonator on the absorption rate is similar to that of a, which can be attributed to the increase in the effective area of graphene and the generation of interactions.The Fermi level of graphene can be continuously tuned by applying a bias voltage or chemical doping on the graphene surface.The absorption spectra at different Fermi levels are shown in Figure 7a.The theoretical calculations show that the five absorption frequency points exhibit an obvious blue shift with the increase in the Fermi level, the absorption intensity at f1 decreases, and the absorption intensities at f2-f5 increase and then decrease.Meanwhile, to clearly illustrate the frequency tunability of the designed graphene metamaterial absorber, we plotted the absorptivity as a function of the Fermi levels and resonance frequencies, as shown in Figure 7b.The relationship between the Fermi level and the graphene-based absorber resonance frequency can be expressed as follows [64]: where 0 = 2 ℏ is the structural constant of graphene, L is the length of the graphene strip, and is the speed of light in a vacuum.Based on the dynamic tunable properties of gra- Based on the need for a compact design, in the process of translating each of the upper two square rings to the center axis, we found interesting phenomena.According to the absorption curves and trends, the number of absorption peaks mutates from three peaks to five peaks.When d = 0 µm, the overall pattern is horizontally symmetric about the central axis, at which time there are three absorption peaks, but the center front absorption quality is low.With the increase in d, the upper two square rings are gradually close to each other, and the coupling between them increases significantly, improving the absorption quality.After reaching the peak, the operating band undergoes a significant blue shift.In order to improve the overall absorption performance of the absorber, based on the original design, the slits are skeletonized at different locations to form splits. From the absorption mechanism analyzed earlier, the absorption strength depends largely on the plasma resonance at the slits.When δ = 0 µm, both splits are located at the centrosymmetric position of the associated geometrical structure, but at this time, the absorption mode is three-peaked and there are oscillating peaks with poor absorption quality.With the change in the position of δ, the absorption intensity reaches the optimal value at δ = 0.15 µm.The thickness of the dielectric layer has a significant effect on the absorption profile.
The Fermi level of graphene can be continuously tuned by applying a bias voltage or chemical doping on the graphene surface.The absorption spectra at different Fermi levels are shown in Figure 7a.The theoretical calculations show that the five absorption frequency points exhibit an obvious blue shift with the increase in the Fermi level, the absorption intensity at f 1 decreases, and the absorption intensities at f 2 -f 5 increase and then decrease.Meanwhile, to clearly illustrate the frequency tunability of the designed graphene metamaterial absorber, we plotted the absorptivity as a function of the Fermi levels and resonance frequencies, as shown in Figure 7b.The relationship between the Fermi level and the graphene-based absorber resonance frequency can be expressed as follows [64]: where α 0 = e 2 ℏ is the structural constant of graphene, L is the length of the graphene strip, and υ is the speed of light in a vacuum.Based on the dynamic tunable properties of graphene, we can flexibly adjust the absorption intensity and operating frequency of the absorber by changing the Fermi level of graphene without changing the geometrical parameters of the absorber.In addition, as shown in Figure 7c, changes in the thickness of the dielectric layer bring about significant changes in the number of absorption peaks.A dielectric layer that is too thin can cause electromagnetic waves that penetrate the entire structure, while a dielectric layer that is too thick can cause multiple oscillations of electromagnetic waves in the cavity.Thus, t 2 reaches the optimum value at 4 µm.results show that the absorber can maintain the absorption stability under a wide incidence angle of 0-65° and exhibits wide-angle incidence characteristics.

Conclusions
In summary, a multi-band tunable terahertz perfect absorber based on a graphene metasurface, which has a typical three-layer sandwich structure, is proposed in this paper.The numerical simulation results confirm the existence of five perfect absorption peaks in Combined with practical engineering applications, we need to consider the incident wave under oblique incidence.In the COMSOL Multiphysics 6.1 software, we parameterized the incidence elevation angle and finally obtained the absorption spectra at different incidence angles, as shown in Figure 8.As the incident angle increases in 10 • steps, different frequency points have different angular dependencies.The absorption performance of f 1 is stable with the increase in the incidence angle, and the intensity remains unchanged between 0 and 70 • .The f 2 and f 3 values decrease in absorption intensity at the absorption angle of 65 • , and the position of the frequency point is almost unchanged.The f 4 and f 5 values undergo a significant decrease in intensity at the incidence angle of more than 60 • , and the frequency point is shifted slightly in comparison with the center position.The results show that the absorber can maintain the absorption stability under a wide incidence angle of 0-65 • and exhibits wide-angle incidence characteristics.

Conclusions
In summary, a multi-band tunable terahertz perfect absorber based on a graphene metasurface, which has a typical three-layer sandwich structure, is proposed in this paper.The numerical simulation results confirm the existence of five perfect absorption peaks in the terahertz operating range of this structure, which can realize high-quality absorption with an average absorption rate of 98.78%.By precisely adjusting the structural parameters of the patterned graphene layer and the thickness of the dielectric layer, the number and intensity of absorption peaks can be adjusted to the desired values.The numerical results show that the absorber has stable absorption performance in the range of 0~65° in the incidence angle with wide-angle incidence.When the Fermi level is adjusted in the range of 1.0-1.4eV, the average peak absorption increases from 96.38% to 98.78%.Additionally, we can realize the continuous tuning of the absorption performance in a wide frequency range by changing the Fermi level of the graphene.Given the excellent properties of the designed absorber, such as high absorption, dynamic tuning, and wide incidence angle, our design has potential applications in the fields of nondestructive testing, imaging, communication, sensing, and detectors.
3 µm.The top view of the patterned graphene cell structure is shown in Figure 1b, and the specific optimized geometrical parameters are as follows: a = 2.0 µm, b = 0.5 µm, c = 1.2 µm, r = 1.0 µm, d = 0.15 µm, and δ = 0.15 µm.The period of the unit cell in the x and y directions is P = 4 µm.

Figure 4 .
Figure 4.The electric field distributions of the absorber at (a) 4.62 THz and (d) 7.29 THz, while (b) and (e) show the respective EZ maps.The current density distributions at (c) 4.62 THz and (f) 7.29 THz.

Figure 4 .
Figure 4.The electric field distributions of the absorber at (a) 4.62 THz and (d) 7.29 THz, while (b) and (e) show the respective E Z maps.The current density distributions at (c) 4.62 THz and (f) 7.29 THz.

Figure 4 .
Figure 4.The electric field distributions of the absorber at (a) 4.62 THz and (d) 7.29 THz, while (b) and (e) show the respective EZ maps.The current density distributions at (c) 4.62 THz and (f) 7.29 THz.

Figure 6 .
Figure 6.Absorption spectra with different values of (a) parameter a, (b) parameter b, (c) parameter c, (d) parameter r, (e) parameter d, and (f) parameter r.The other default parameters are the same as those illustrated in Figure 1.

Figure 7 .
Figure 7. (a) Absorption spectra and (b) frequency shifts at different Fermi levels.The absorption spectra with (c) different dielectric layer thickness t2.

Figure 8 .
Figure 8. Absorption as a function of frequency and angle of incidence.

Figure 7 .
Figure 7. (a) Absorption spectra and (b) frequency shifts at different Fermi levels.The absorption spectra with (c) different dielectric layer thickness t 2 .

Figure 7 .
Figure 7. (a) Absorption spectra and (b) frequency shifts at different Fermi levels.The absorption spectra with (c) different dielectric layer thickness t2.

Figure 8 .
Figure 8. Absorption as a function of frequency and angle of incidence.