Numerical Investigation of Graphene and STO Based Tunable Terahertz Absorber with Switchable Bifunctionality of Broadband and Narrowband Absorption

A graphene metamaterial and strontium titanate (STO)-based terahertz absorber with tunable and switchable bifunctionality has been numerically investigated in this work. Through electrically tuning the Fermi energy level of the cross-shaped graphene, the bandwidth of the proposed absorber varies continuously from 0.12 THz to 0.38 THz with the absorptance exceeding 90%, which indicates the functionality of broadband absorption. When the Fermi energy level of the cross-shaped graphene is 0 eV, the proposed absorber exhibits the other functionality of narrowband absorption owing to the thermal control of the relative permittivity of STO, and the rate of change of the center frequency is 50% ranging from 0.56 THz to 0.84 THz. The peak intensity of the narrowband absorption approximates to nearly 100% through adjusting the Fermi energy level of the graphene strips. The calculated results indicate that it is not sensitive to the polarization for wide incidence angles. The proposed absorber can realize tunable bifunctionality of broadband absorption with a tunable bandwidth and narrowband absorption with a tunable center frequency, which provides an alternative design opinion of the tunable terahertz devices with high performance for high-density integrated systems.


Introduction
Metamaterials, artificially engineered by subwavelength electromagnetic materials, show some optical properties that differ from natural materials. The resonance frequency can be arbitrarily customized in a large frequency range in the microwave [1,2], terahertz (THz) [3][4][5], and near-infrared regions [6]. Metamaterial-based absorbers (MMAs) have attracted great interest due to their scalable properties and a wide variety of potential applications, for example solar and thermophotovoltaic energy conversion [7,8], cloaking [9,10], sensors [11], and thermal emitters [12,13]. Conventional MMAs are composed of a sandwich structure with a dielectric spacer, which acts as a divider between the subwavelength metallic patterns and a continuous metallic ground plane. Since the first perfect MMA was proposed and investigated [14], MMAs with single- [15,16], dual- [17], multi- [18,19], and broad-band [20][21][22] absorption have been proposed and investigated. However, the absorption performance of MMAs influenced by the geometric parameters of unit cells cannot be dynamically adjusted anymore once the fabrication is completed.
To enrich the functionalities, active materials are introduced into metamaterial devices, such as vanadium dioxide [23], phase-change materials [24], liquid crystals [25], liquid metals [26], and especially graphene [27,28]. Due to its extraordinary physical properties [29,30], such as high electron mobility, flexible tunability [31], relatively low loss, In this study, a graphene metamaterial and STO-based absorber with switchable and tunable bifunctionality is proposed in the terahertz regime. Through controlling the Fermi energy level of the cross-shaped graphene, a tunable broadband absorption of the bifunctional absorber can be obtained, and its bandwidth can realize a variation in the range of frequencies from 0.12 to 0.38 THz with excellent absorptance of over 90%. The relative impedance as well as the field analyses are investigated to reveal its absorbing mechanism. When the Fermi energy level of the cross-shaped graphene is 0 eV, a narrowband absorption with tunable center frequencies can be achieved utilizing thermal control of the relative permittivity of STO. Furthermore, the effects of the geometrical parameters and various incidence angles for different polarization on absorption performances are also discussed. The bifunctional absorber inspires the design of dynamically tunable devices with multiple functionalities in the terahertz regime.

Materials and Methods
The structural schematic and polarization configuration of the incident wave of the bifunctional absorber is shown in Figure 1a, which is composed of a cross-shaped graphene metamaterial layer, an insulating spacer (polyethylene cyclic olefin copolymer, Topas) [48], a 1.9-µm-thick STO layer covered by strip-shaped graphene metamaterial, and a bottom gold film with conductivity of 4.56 × 10 7 S/m used as the continuous metallic reflector. The thickness of the gold film is 0.5 µm, which is much greater than the maximum skin depth conductivity of the gold film. Topas, a transparent and stiff amorphous thermoplastic copolymer, shows superior optical properties such as high stability, excellent heat resistance, negligible absorption coefficient, and constant refractive index in the THz range [49]. The relative permittivity of Topas is 2.35 with negligible loss and dispersion in this work [50]. Thicknesses of the insulating spacer h g and h d are 40 µm and 6.4 µm, respectively. Figure 1b,c shows the upper and lower metamaterial of the bifunctional absorber, respectively. The branch length and width of the cross-shaped graphene is l s = 36 µm and w s = 6 µm, respectively, and the width of the graphene strips is w r = 42 µm. The periods of the unit cell for the bifunctional absorber are P x = P y = 60 µm. The side view is depicted in Figure 1d, where the cross-shaped graphene metamaterial is covered by a 20-µm-thick electrolyte (poly-(ethylene oxide)/LiClO 4 ) layer with refractive index 1.7 [51].
impedance as well as the field analyses are investigated to reveal its absorbing mechanism. When the Fermi energy level of the cross-shaped graphene is 0 eV, a narrowband absorption with tunable center frequencies can be achieved utilizing thermal control of the relative permittivity of STO. Furthermore, the effects of the geometrical parameters and various incidence angles for different polarization on absorption performances are also discussed. The bifunctional absorber inspires the design of dynamically tunable devices with multiple functionalities in the terahertz regime.

Materials and Methods
The structural schematic and polarization configuration of the incident wave of the bifunctional absorber is shown in Figure 1a, which is composed of a cross-shaped graphene metamaterial layer, an insulating spacer (polyethylene cyclic olefin copolymer, Topas) [48], a 1.9-μm-thick STO layer covered by strip-shaped graphene metamaterial, and a bottom gold film with conductivity of 4.56 × 10 7 S/m used as the continuous metallic reflector. The thickness of the gold film is 0.5 μm, which is much greater than the maximum skin depth conductivity of the gold film. Topas, a transparent and stiff amorphous thermoplastic copolymer, shows superior optical properties such as high stability, excellent heat resistance, negligible absorption coefficient, and constant refractive index in the THz range [49]. The relative permittivity of Topas is 2.35 with negligible loss and dispersion in this work [50]. Thicknesses of the insulating spacer hg and hd are 40 μm and 6.4 μm, respectively. Figure 1b,c shows the upper and lower metamaterial of the bifunctional absorber, respectively. The branch length and width of the cross-shaped graphene is ls = 36 μm and ws = 6 μm, respectively, and the width of the graphene strips is wr = 42 μm. The periods of the unit cell for the bifunctional absorber are Px = Py = 60 μm. The side view is depicted in Figure 1d, where the cross-shaped graphene metamaterial is covered by a 20μm-thick electrolyte (poly-(ethylene oxide)/LiClO4) layer with refractive index 1.7 [51]. The CVD-grown graphene layer is transferred onto the multilayer substrate by a transfer technique using polymethylmethacrylate (PMMA) supporting layers and is subsequently patterned by photolithography and oxygen plasma etching [52]. The graphene The CVD-grown graphene layer is transferred onto the multilayer substrate by a transfer technique using polymethylmethacrylate (PMMA) supporting layers and is subsequently patterned by photolithography and oxygen plasma etching [52]. The graphene layer can be modeled as a surface current defined as J = σ g E t according to Ohm's law, where E t is the tangential electric field, and σ g is the complex conductivity of graphene. The surface conductivity of graphene from Kubo formula is determined by the combination of intraband and interband contributions [53]: where ω is the incident angular frequency, k B ≈ 1.381 × 10 −23 J/K is the Boltzmann constant, e ≈ 1.602 × 10 −19 C is electron charge,h ≈ 1.055 × 10 −34 , J·s is the reduced Planck constant, and T is temperature in kelvin. The values of chemical potential µ c and Fermi energy level E f are equal when k B < µ c . The relaxation time τ can be given by τ = µE f e −1 υ F −2 , where the Fermi velocity υ F is 10 6 m/s, µ is the carrier mobility assumed as 1310 cm 2 V −1 s −1 for the graphene strips and 26,250 cm 2 V −1 s −1 for the crossshaped graphene, and the corresponding relaxation time is 0.105 ps and 2.1 ps with E f = 0.8 eV, respectively. The complex relative permittivity of STO material is sensitive to the temperature, which can be expressed as follows [40,54] where k is the wave number of the incident wave, ε ∞ = 9.6 is the high-frequency bulk permittivity, k 0 = 31.2(T − 42.5) cm −1 is the wave number of the soft mode, f = 2.3 × 10 6 cm −2 represents the oscillator strength, and γ = −3.3 + 0.094T cm −1 is the damping parameter of the soft mode. The external temperature T is set as 400 K in this work except when stated otherwise. In order to gain deep insight into the temperature-dependent property of STO, the real and imaginary parts of the permittivity as a function of frequency are calculated with various temperatures, as shown in Figure 2a,b, respectively. Both the real part and the imaginary part decreases with the increasing temperature. However, the value of the imaginary part is much smaller than that of the real part. The resonance frequency is mainly affected by the real part, and the imaginary part characterizes the losses. Therefore, the resonance frequency caused by the STO material can be significantly shifted by adjusting the environment temperature, while the intensity of the absorption changes little.
where Et is the tangential electric field, and σg is the complex conductivity of graphene. The surface conductivity of graphene from Kubo formula is determined by the combination of intraband and interband contributions [53]: where ω is the incident angular frequency, kB ≈ 1.381 × 10 −23 J/K is the Boltzmann constant, e ≈ 1.602 × 10 −19 C is electron charge, ħ ≈ 1.055 × 10 −34 , J·s is the reduced Planck constant, and T is temperature in kelvin. The values of chemical potential μc and Fermi energy level Ef are equal when kB < μc. The relaxation time τ can be given by Fermi velocity υF is 10 6 m/s, μ is the carrier mobility assumed as 1310 cm 2 V −1 s −1 for the graphene strips and 26,250 cm 2 V −1 s −1 for the cross-shaped graphene, and the corresponding relaxation time is 0.105 ps and 2.1 ps with Ef = 0.8 eV, respectively. The complex relative permittivity of STO material is sensitive to the temperature, which can be expressed as follows [40,54] where k is the wave number of the incident wave, ε∞ = 9.6 is the high-frequency bulk per- represents the oscillator strength, and γ = −3.3 + 0.094T cm −1 is the damping parameter of the soft mode. The external temperature T is set as 400 K in this work except when stated otherwise.
In order to gain deep insight into the temperature-dependent property of STO, the real and imaginary parts of the permittivity as a function of frequency are calculated with various temperatures, as shown in Figure 2a,b, respectively. Both the real part and the imaginary part decreases with the increasing temperature. However, the value of the imaginary part is much smaller than that of the real part. The resonance frequency is mainly affected by the real part, and the imaginary part characterizes the losses. Therefore, the resonance frequency caused by the STO material can be significantly shifted by adjusting the environment temperature, while the intensity of the absorption changes little.

Results and Discussion
The electromagnetic absorptance(A), i.e., the intensity of absorption, can be defined as A = 1 − T − R, in which transmittance T = |S 21 | 2 and reflectance R = |S 11 | 2 can be obtained from S-parameters in the simulation results calculated by COMSOL Multiphysics. When the Fermi energy level of graphene strips E f2 is set as 0.8 eV and that of the crossshaped graphene E f1 is adjusted from 0.4 to 0.8 eV, the proposed absorber shows a tunable broadband absorption, as depicted in Figure 2. The bandwidth with the absorptance above 90% is 0.38 THz in the frequency range from 0.78 to 1.16 THz. As E f1 decreases to 0.4 eV, the bandwidth grows narrower and reaches the minimum 0.12 THz. The absorption broadband disappears and divides to become two absorption bands when E f1 continues to decrease, as depicted in Figure 3 by the dashed curves. Therefore, the bandwidth can Nanomaterials 2021, 11, 2044 5 of 12 be dynamically and continually tuned in the range of frequencies from 0.12 to 0.38 THz through controlling E f1 .
as A = 1−T−R, in which transmittance T = |S21| 2 and reflectance R = |S11| 2 can be obtained from S-parameters in the simulation results calculated by COMSOL Multiphysics. When the Fermi energy level of graphene strips Ef2 is set as 0.8 eV and that of the cross-shaped graphene Ef1 is adjusted from 0.4 to 0.8 eV, the proposed absorber shows a tunable broadband absorption, as depicted in Figure 2. The bandwidth with the absorptance above 90% is 0.38 THz in the frequency range from 0.78 to 1.16 THz. As Ef1 decreases to 0.4 eV, the bandwidth grows narrower and reaches the minimum 0.12 THz. The absorption broadband disappears and divides to become two absorption bands when Ef1 continues to decrease, as depicted in Figure 3 by the dashed curves. Therefore, the bandwidth can be dynamically and continually tuned in the range of frequencies from 0.12 to 0.38 THz through controlling Ef1. To understand how the absorption broadbands are formed, the influences on absorption spectra of the individual upper and lower metamaterials are calculated, respectively. The absorption performance of the individual cross-shaped graphene with Ef1 varying from 0.1 eV to 0.8 eV is shown in Figure 4a. When Ef1 increases, the center frequency shows a monotonous blue-shift. The absorption enhancement starts to increase, then decreases, and achieves the maximum absorptance with Ef1 = 0.3 eV. Thus, Ef1 not only affects the location of the center frequency, but also the intensity of absorption. As shown in Figure  4b, the absorption intensity of the individual lower layer of STO combined with graphene strips is changed through adjusting Ef2, while the variation of the center frequency is very small, which indicates that Ef2 only influences the intensity of the absorption. The impedance-matching theory can be used as a physical explanation the effect of Fermi energy level on the absorption intensity of the individual upper and lower metamaterials. Consequently, due to the impedance matching between the free space and the proposed absorber, the upper and lower metamaterials combined effectively, which contributes to the stable high absorptance and the wide operating bandwidth, as depicted in Figure 2 by the solid curves. To understand how the absorption broadbands are formed, the influences on absorption spectra of the individual upper and lower metamaterials are calculated, respectively. The absorption performance of the individual cross-shaped graphene with E f1 varying from 0.1 eV to 0.8 eV is shown in Figure 4a. When E f1 increases, the center frequency shows a monotonous blue-shift. The absorption enhancement starts to increase, then decreases, and achieves the maximum absorptance with E f1 = 0.3 eV. Thus, E f1 not only affects the location of the center frequency, but also the intensity of absorption. As shown in Figure 4b, the absorption intensity of the individual lower layer of STO combined with graphene strips is changed through adjusting E f2 , while the variation of the center frequency is very small, which indicates that E f2 only influences the intensity of the absorption. The impedance-matching theory can be used as a physical explanation the effect of Fermi energy level on the absorption intensity of the individual upper and lower metamaterials. Consequently, due to the impedance matching between the free space and the proposed absorber, the upper and lower metamaterials combined effectively, which contributes to the stable high absorptance and the wide operating bandwidth, as depicted in Figure 2 by the solid curves.  The distributions of the amplitude of electric field |E| and the z-component of electric field Ez for the individual upper metamaterials are depicted in Figure 5a,b, respectively, at the resonance frequency of 0.78 THz with Ef1 = 0.3 eV. The amplitude of |E| concentrates mainly around the ends and edges of the horizontal branches for the cross-shaped graphene. The distribution of Ez shows that opposite charges accumulate at both ends of the horizontal branches, which indicates a typical electric dipole resonance. Figure 5c describes the distributions of magnetic field |H| as well as the surface current (black arrows)  The distributions of the amplitude of electric field |E| and the z-component of electric field E z for the individual upper metamaterials are depicted in Figure 5a,b, respectively, at the resonance frequency of 0.78 THz with E f1 = 0.3 eV. The amplitude of |E| concentrates mainly around the ends and edges of the horizontal branches for the cross-shaped graphene.
The distribution of E z shows that opposite charges accumulate at both ends of the horizontal branches, which indicates a typical electric dipole resonance. Figure 5c describes the distributions of magnetic field |H| as well as the surface current (black arrows) for the individual lower metamaterial at the resonance frequency 0.85 THz with E f2 set as 0.7 eV. The directions of the currents of the unit cell without the bottom gold block are marked by black arrows, while the magnetic field is concentrated below the STO layer, which demonstrates that there are induced currents on the gold layer. Then, a strong magnetic resonance is caused by these parallel surface currents with opposite direction. Therefore, the resonance absorptions of the individual upper and lower metamaterials can be attributed to the electric dipole resonance and the magnetic resonance, respectively. The distributions of the amplitude of electric field |E| and the z-component of electric field Ez for the individual upper metamaterials are depicted in Figure 5a,b, respectively, at the resonance frequency of 0.78 THz with Ef1 = 0.3 eV. The amplitude of |E| concentrates mainly around the ends and edges of the horizontal branches for the cross-shaped graphene. The distribution of Ez shows that opposite charges accumulate at both ends of the horizontal branches, which indicates a typical electric dipole resonance. Figure 5c describes the distributions of magnetic field |H| as well as the surface current (black arrows) for the individual lower metamaterial at the resonance frequency 0.85 THz with Ef2 set as 0.7 eV. The directions of the currents of the unit cell without the bottom gold block are marked by black arrows, while the magnetic field is concentrated below the STO layer, which demonstrates that there are induced currents on the gold layer. Then, a strong magnetic resonance is caused by these parallel surface currents with opposite direction. Therefore, the resonance absorptions of the individual upper and lower metamaterials can be attributed to the electric dipole resonance and the magnetic resonance, respectively. The amplitude distributions of electric field and magnetic field at the resonance frequencies under normal TM polarized incidence are shown in Figure 6 to further investigate the merging of its upper and lower layer of metamaterials. When E f1 is 0.7 eV and E f2 is 0.8 eV, the first resonance frequency is 0.88 THz, and the second resonance frequency is 1.06 THz. The electric dipole resonance and the magnetic resonance coexist to originate the broadband absorption. For the cross-shaped graphene, the electric dipole resonance at the second resonance is much stronger than that at the first resonance, as shown in Figure 6a,b. Figure 6c,d shows that the magnetic resonance below the STO layer at 0.88 THz is more concentrated than that at 1.06 THz. Consequently, the upper metamaterial contributes more to the broadband absorbing performance at the second resonance, while the lower metamaterial contributes more to the broadband absorption at the first resonance.
1.06 THz. The electric dipole resonance and the magnetic resonance coexist to originate the broadband absorption. For the cross-shaped graphene, the electric dipole resonance at the second resonance is much stronger than that at the first resonance, as shown in Figure 6a,b. Figure 6c,d shows that the magnetic resonance below the STO layer at 0.88 THz is more concentrated than that at 1.06 THz. Consequently, the upper metamaterial contributes more to the broadband absorbing performance at the second resonance, while the lower metamaterial contributes more to the broadband absorption at the first resonance. To comprehend the influence of the geometrical parameters on the broadband absorbing properties, the absorption spectra that vary with various parameters are simulated. The absorption spectra that vary with the thickness hd of the Topas layer between STO and the metallic ground is shown in Figure 7a. The location of the first resonance clearly shows a red-shift and the bandwidth grows wider gradually, which can be attributed to the fact that the increasing hd will affect the magnetic resonance below the STO layer. It is evident that the red shift of the second resonance frequency changes slightly and the bandwidth remains almost unchanged with the increase of the insulating spacer hg because the upper metamaterial contributes more to the broadband absorption at the second resonance, as shown in Figure 7b. When the thickness of the STO material layer hSTO varies from 1.7 μm to 2.1 μm, the first absorption resonance has red shift, which results in the increment of bandwidth, as shown in Figure 7c. The absorption spectra of the branch length ls increasing from 34 μm to 38 μm is shown in Figure 7d. It is obvious that the second resonance has red shift with the increment of ls. That is because the resonance frequency varies inversely to the effective length of the patterned structure on the basis of the LC circuit model. Figure 7e shows that the bandwidth with an absorptance over 90% decreases slightly when the width of the graphene strip wr changes from 31 μm to 53 μm. The absorption spectra show quite a change when the temperature is lower than 375 K, which can be attributed to the influence of temperature on the conductivity of the STO material, as the impedance matching is destroyed with the drop in temperature, as shown in Figure 7f. To comprehend the influence of the geometrical parameters on the broadband absorbing properties, the absorption spectra that vary with various parameters are simulated. The absorption spectra that vary with the thickness h d of the Topas layer between STO and the metallic ground is shown in Figure 7a. The location of the first resonance clearly shows a red-shift and the bandwidth grows wider gradually, which can be attributed to the fact that the increasing h d will affect the magnetic resonance below the STO layer. It is evident that the red shift of the second resonance frequency changes slightly and the bandwidth remains almost unchanged with the increase of the insulating spacer h g because the upper metamaterial contributes more to the broadband absorption at the second resonance, as shown in Figure 7b. When the thickness of the STO material layer h STO varies from 1.7 µm to 2.1 µm, the first absorption resonance has red shift, which results in the increment of bandwidth, as shown in Figure 7c. The absorption spectra of the branch length l s increasing from 34 µm to 38 µm is shown in Figure 7d. It is obvious that the second resonance has red shift with the increment of l s . That is because the resonance frequency varies inversely to the effective length of the patterned structure on the basis of the LC circuit model. Figure 7e shows that the bandwidth with an absorptance over 90% decreases slightly when the width of the graphene strip w r changes from 31 µm to 53 µm. The absorption spectra show quite a change when the temperature is lower than 375 K, which can be attributed to the influence of temperature on the conductivity of the STO material, as the impedance matching is destroyed with the drop in temperature, as shown in Figure 7f.
Furthermore, Figure 8a shows that the first absorption peak decreases with the decline of E f2 , while the location of the absorption peak remains unchanged. The Fermi energy level of the graphene strips only influences the flatness of the broadband absorption, and the impedance matching theory can be used to explained this phenomenon. The intensity of the absorption can reach the maximum when the relative impedance is 1. Figure 8b,c denotes the real and imaginary parts of the relative impedance Z r with E f1 changing in the range of 0.4-0.8 eV. When E f2 is set as 0.8 eV, the real part of Z r approaches 1, and the imaginary part of Z r approaches 0, respectively, in the frequency range from 0.78 to 1.16 THz, which indicates that the impedance matching of the proposed absorber and the free space, i.e., the maximum absorption intensity, can be achieved. The range of matching frequencies falls as the Fermi energy level decreases. The absorption bandwidth reaches the minimum as the Fermi energy level decreases to 0.4 eV.
In addition, the broadband absorber can be transformed into a narrowband absorber when the E f1 is fixed as zero, as depicted in Figure 9. Attributed to the thermal control relative permittivity of STO, the center frequency is clearly blue-shifted from 0.58 THz to 0.88 THz when the external temperature rises from 200 to 400 K with E f2 = 0.8 eV, as shown by the dashed curves in Figure 9. The absorptance can be further improved through tuning E f2 . The peak absorption intensity is over 96% at various temperatures with E f2 = 0.3 eV, as shown by the solid curves. Thus, a perfect narrowband absorber with a dynamically tunable center frequency can be achieved. Furthermore, Figure 8a shows that the first absorption peak decreases with the decline of Ef2, while the location of the absorption peak remains unchanged. The Fermi energy level of the graphene strips only influences the flatness of the broadband absorption, and the impedance matching theory can be used to explained this phenomenon. The intensity of the absorption can reach the maximum when the relative impedance is 1. Figure  8b,c denotes the real and imaginary parts of the relative impedance Zr with Ef1 changing in the range of 0.4-0.8 eV. When Ef2 is set as 0.8 eV, the real part of Zr approaches 1, and the imaginary part of Zr approaches 0, respectively, in the frequency range from 0.78 to 1.16 THz, which indicates that the impedance matching of the proposed absorber and the free space, i.e., the maximum absorption intensity, can be achieved. The range of matching frequencies falls as the Fermi energy level decreases. The absorption bandwidth reaches the minimum as the Fermi energy level decreases to 0.4 eV. In addition, the broadband absorber can be transformed into a narrowband absorber when the Ef1 is fixed as zero, as depicted in Figure 9. Attributed to the thermal control relative permittivity of STO, the center frequency is clearly blue-shifted from 0.58 THz to 0.88 THz when the external temperature rises from 200 to 400 K with Ef2 = 0.8 eV, as shown by the dashed curves in Figure 9. The absorptance can be further improved through tuning Ef2. The peak absorption intensity is over 96% at various temperatures with Ef2 = 0.3 eV, as shown by the solid curves. Thus, a perfect narrowband absorber with a dynamically tunable center frequency can be achieved. The effects of oblique incidence under both transverse electric (TE) and transverse magnetic (TM) waves on the absorption performance of the bifunctional absorber are finally investigated. The refraction angle θ m in layer m with the incident angle θ i can be expressed as θ m = arcsin(n 1 sin(θ i )/n m ) based on Snell's law, where n 1 and n m are the refractive index of the first and m-th layers. It is noteworthy that the refraction angle is a complex value when the m-th medium is a lossy material, which reflects that the refracted wave in the lossy material is a non-uniform plane wave [55]. Therefore, the overall reflection of the proposed absorber is then the superposition of the multiple reflections at different interfaces.
In addition, the broadband absorber can be transformed into a narrowband absorber when the Ef1 is fixed as zero, as depicted in Figure 9. Attributed to the thermal control relative permittivity of STO, the center frequency is clearly blue-shifted from 0.58 THz to 0.88 THz when the external temperature rises from 200 to 400 K with Ef2 = 0.8 eV, as shown by the dashed curves in Figure 9. The absorptance can be further improved through tuning Ef2. The peak absorption intensity is over 96% at various temperatures with Ef2 = 0.3 eV, as shown by the solid curves. Thus, a perfect narrowband absorber with a dynamically tunable center frequency can be achieved. As shown in Figure 10a, when the bifunctional absorber acts as a broadband absorber, for the TE wave, more than 80% absorptance can be achieved for the incidence angle is up to 55 • , and the center operating frequency remains unchanged. The closer the incidence angle to 90 • , the smaller the absorptance. For TM polarization, the absorptance remains over 80% for the incidence angle up to 80 • , as depicted in Figure 10b. When the incidence angle is over 40 • , the center frequency shows blue shift clearly. The influence on absorption spectra is attributed to the change in zero-reflection condition under various incidence angles. For the narrowband absorption, the temperature is fixed at 400 K, and E f2 is set as 0.3 eV. As shown in Figure 10c,d, the absorptance remains above 80% with the incidence angle smaller than 80 • for both TE and TM waves. With the increase in the incidence angle, the center frequency remains unchanged for the TE wave, while a clear blue shift is observed for TM polarization, which is consistent with the results in Figure 10a,b. Consequently, the bifunctional absorber is not sensitive to the polarization for wide oblique incidence angles. wave in the lossy material is a non-uniform plane wave [55]. Therefore, the overall reflection of the proposed absorber is then the superposition of the multiple reflections at different interfaces. As shown in Figure 10a, when the bifunctional absorber acts as a broadband absorber, for the TE wave, more than 80% absorptance can be achieved for the incidence angle is up to 55°, and the center operating frequency remains unchanged. The closer the incidence angle to 90°, the smaller the absorptance. For TM polarization, the absorptance remains over 80% for the incidence angle up to 80°, as depicted in Figure 10b. When the incidence angle is over 40°, the center frequency shows blue shift clearly. The influence on absorption spectra is attributed to the change in zero-reflection condition under various incidence angles. For the narrowband absorption, the temperature is fixed at 400 K, and Ef2 is set as 0.3 eV. As shown in Figure 10c,d, the absorptance remains above 80% with the incidence angle smaller than 80° for both TE and TM waves. With the increase in the incidence angle, the center frequency remains unchanged for the TE wave, while a clear blue shift is observed for TM polarization, which is consistent with the results in Figure 10a,b. Consequently, the bifunctional absorber is not sensitive to the polarization for wide oblique incidence angles.

Conclusions
In summary, a bifunctional absorber based on a graphene metamaterial and STO is designed with dynamically tunable and switchable properties. Compared with some other's previous work, the bifunctional absorber presents excellent tunable ability. When Ef1 varies from 0.4 eV to 0.8 eV, a broadband absorption is achieved, and its bandwidth

Conclusions
In summary, a bifunctional absorber based on a graphene metamaterial and STO is designed with dynamically tunable and switchable properties. Compared with some other's previous work, the bifunctional absorber presents excellent tunable ability. When E f1 varies from 0.4 eV to 0.8 eV, a broadband absorption is achieved, and its bandwidth varies from 0.12 THz to 0.38 THz with the intensity exceeding 90%. The results indicate that the broadband absorption benefits from the combination of graphene metamaterial and STO material. For the broadband absorption, the rate of change of the bandwidth is 216.6%, and that of the center frequency is 21%. The broadband absorber can be transformed into a narrowband absorber, when E f1 is set as 0. The center frequency of the narrowband absorption can be adjusted from 0.56 THz to 0.84 THz by controlling the temperature of STO. The peak intensity of the narrowband absorption is approximately 100% when E f2 is set as 0.3 eV, and the rate of change of the center frequency is 50%. The bifunctional absorber is not sensitive to polarization for large incidence angles. Hence, the proposed absorber may be suitable for many potential applications, for example sensing, optical cloaking, and some other switchable devices.

Conflicts of Interest:
The authors declare no conflict of interest.