High-Performance Terahertz Coherent Perfect Absorption with Asymmetric Graphene Metasurface

: In this work, we introduce a novel coherent perfect absorber, accentuating its novelty by emphasizing the broad bandwidth, reduced thickness, tunable property


Introduction
High technology advancements have made electromagnetic absorption a major research focus [1].Coherent perfect absorption (CPA) is a promising phenomenon where two coherent light beams, illuminating from opposite sides of a medium or structure, interfere in such a way that the medium or structure absorbs all the incoming light, resulting in no reflection or transmission [2][3][4][5].This process hinges on the precise control of the phase and amplitude of the light beams to achieve perfect absorption through coherent interference.CPA can be realized through various structures, including thin films and waveguides, offering more absorption tunability through interference interactions compared to perfect absorbers that rely on strong resonance [6,7].Numerous CPA structures have been designed and experimentally validated for a wide range of frequencies, from microwave to optical [8][9][10][11][12][13].This makes them highly promising for applications in electromagnetic modulation, switching, and signal processing.
The simplest structure for achieving CPA is a thin film or flat plate.Two-dimensional (2D) materials [14], such as graphene, are ideal for this purpose due to their ultrathin profile and excellent electrical properties [15][16][17].Graphene, a 2D mono-layer of carbon atoms, has exceptional optical and electrical properties [18,19].It has been widely used in various optical devices, including polarizers, phase shifters, and absorbers, due to its ability to manipulate the interaction between light and matter [20][21][22].However, the absorption rate of single-layer graphene is only 2.3% [18,23], which is determined by its thin-layer conductivity.To improve the bandwidth and reduce the thickness of the CPA structure, researchers have explored the theory of CPA in patterned graphene structures, such as square graphene ring-and T-shaped graphene structures [24][25][26][27].The graphene metasurface-based CPA system enhances incident light absorption and widens the absorption bandwidth [25,28,29].An asymmetric graphene metasurface structure has been developed to achieve CPA, and it has a 90% absorption bandwidth of 1.5 THz and a thickness of only 2 um [30].Additionally, the performance of the graphene metasurface can be dynamically controlled by varying the graphene Fermi level through using electrostatic gating or chemical doping, resulting in a wide frequency tunable range [31,32].Moreover, the graphene metasurface structure is also suitable for small-scale integration.
In this study, an ultrabroadband THz CPA structure is proposed based on an asymmetric graphene metasurface, which consists of both square and circular graphene patches arranged on opposite sides of a silicon substrate.The results show that with a graphene Fermi level of 0.8 eV, this designed asymmetric structure can achieve an over 90% absorption bandwidth of 2.84 THz, ranging from 1.65 to 4.49 THz.In comparison to previously reported CPA structures, the absorption performance of this new metamaterial absorber is significantly enhanced, while its structure is relatively simpler.The effects of phase difference between two coherent lights, the structural geometry, and the incident angle on the CP performance are demonstrated in detail.Importantly, we achieve CPA tunability by varying the graphene Fermi level.This ultrathin, wideband coherent perfect absorption device is expected to have a significant impact on the development of on-chip THz communication technologies such as light modulators, sensors, and photodetectors.

CPA Theory and Design of Ultrabroadband CPA with Graphene Metasurfaces
Figure 1a depicts the configuration of a unit cell of the graphene metasurface, which consists of arrays of square and circular graphene patches positioned on opposite sides of a silicon substrate.It is illuminated by two counter-propagating and coherently modulated optical beams, designated as E in1 and E in2 , with E out1 and E out2 representing the corresponding output waves.By optimizing the proposed structure to have specific reflection and transmission coefficients from singe-side illumination, the CPA conditions can be satisfied.We need to deduce the CPA conditions via the scattering matrix S by linking the outgoing waves E out1 and E out2 with the incident waves E in1 and E in2 .The formula can be described as follows [33,34]: where r 1 , r 2 , t 1 , and t 2 are the reflection coefficients and transmission coefficients of incident waves E in1 and E in2 , excited on the front and back sides of the metasurface structure, respectively.Since this configuration does not break the reciprocity, the transmission coefficient is equal from both incident directions, i.e., t = t 1 = t 2 .For symmetric coherent illumination, we have E in1 = E in2 e iθ , where θ is the phase difference between the two incident waves.To quantitatively investigate CPA, we define the absorption coefficients as: By deriving from Equation (2), the absorption coefficient can be simplified to the following: The CPA operation implies that the absorption coefficients are equal to 1.As a result, the conditions to achieve CPA should be r 1 = −te −iθ and r 2 = −te iθ .In order to achieve CPA, we need to optimize the proposed structure to meet the CPA conditions.This configuration in Figure 1a induces considerable structural asymmetry, which, in turn, broadens the bandwidth of the CPA response.Design optimization was conducted through numerical simulations using COMSOL Multiphysics.In the numerical simulations, a 3D model was constructed, and periodic boundary conditions were applied along the x and y directions, reflecting the finite periodicity of the proposed structure in these dimensions.Port boundary conditions are placed in the z direction to create the incident plane wave, which is assumed to be polarized along the x-axis, i.e., transverse magnetic (TM) polarization.The parameters of the optimized design are detailed in Figure 1a; the square graphene patch length l is 5 um, the circular graphene patch radius r is 1.5 um, the silicon substrate thickness t is 1.5 um, and the structure's period p is 5.2 um.The proposed CPA structure can be easily fabricated in potential experiments.Single-layer graphene can be cultivated on copper substrates via chemical vapor deposition.Graphene in square and circular forms can be achieved through the use of photoresist techniques, electron beam exposure, and plasma etching processes [35].Furthermore, silicon films of micron thickness can be acquired by mechanically grinding, polishing, and etching thick monocrystalline silicon [36].To modulate the graphene metasurface, an ion gel with a low relative permittivity of 1.82 can be employed as the gate dielectric material atop the graphene metasurfaces.Two gold gate contacts can be fabricated on the ion gel layers [37][38][39][40], as depicted in Figure 1b, to apply voltage, dope the graphene, and thereby tune the absorption response.Owing to the thin thickness and low permittivity of the ion gel layers, their impact on the system's performance was determined to be negligible (not illustrated here) and was therefore disregarded in our simulations.
The CPA operation implies that the absorption coefficients are equal to 1.As a result, the conditions to achieve CPA should be 1  This configuration in Figure 1a induces considerable structural asymmetry, which, in turn, broadens the bandwidth of the CPA response.Design optimization was conducted through numerical simulations using COMSOL Multiphysics.In the numerical simulations, a 3D model was constructed, and periodic boundary conditions were applied along the x and y directions, reflecting the finite periodicity of the proposed structure in these dimensions.Port boundary conditions are placed in the z direction to create the incident plane wave, which is assumed to be polarized along the x-axis, i.e., transverse magnetic (TM) polarization.The parameters of the optimized design are detailed in Figure 1a; the square graphene patch length l is 5 um, the circular graphene patch radius r is 1.5 um, the silicon substrate thickness t is 1.5 um, and the structure's period p is 5.2 um.The proposed CPA structure can be easily fabricated in potential experiments.Single-layer graphene can be cultivated on copper substrates via chemical vapor deposition.Graphene in square and circular forms can be achieved through the use of photoresist techniques, electron beam exposure, and plasma etching processes [35].Furthermore, silicon films of micron thickness can be acquired by mechanically grinding, polishing, and etching thick monocrystalline silicon [36].To modulate the graphene metasurface, an ion gel with a low relative permittivity of 1.82 can be employed as the gate dielectric material atop the graphene metasurfaces.Two gold gate contacts can be fabricated on the ion gel layers [37][38][39][40], as depicted in Figure 1b, to apply voltage, dope the graphene, and thereby tune the absorption response.Owing to the thin thickness and low permittivity of the ion gel layers, their impact on the system's performance was determined to be negligible (not illustrated here) and was therefore disregarded in our simulations.
Due to graphene's planar nature, it is modeled as a surface current in our simulations, Due to graphene's planar nature, it is modeled as a surface current in our simulations, described as J = σ g (w)E, where σ g is the conductivity of the graphene, and E is the electric field along the graphene surface.Graphene conductivity can be modulated by chemical doping or external electric field bias, and it is defined as follows [41]: where σ intra (ω) and σ inter (ω) represent the intra-and inter-band transition conductivities, e is the electron charge, k B is the Boltzmann constant, is the Plank's constant, T is the Kelvin temperature, τ is the relaxation time, and E F is the Fermi level.In the low-THz range, intra-band transitions dominate the graphene response, and the graphene conductivity can be expressed by use of the Drude model [42]: Here, the relaxation time is assumed to be 10 −13 s.

Results of Broadband Absorber
Here, we initially investigate single illumination on the proposed asymmetric metasurface.The spectral responses, including reflection, transmission, and absorption, under normal incidence for both front and back illuminations are depicted in Figures 2a and Figure 2b, respectively.The Fermi level of graphene E F is 0.8 eV.It is clear that the transmission from both sides is identical because reciprocity is not broken with this designed configuration.However, reflection and absorption are different due to the structural asymmetry.As illustrated in Figure 2a,b, at a frequency of 2.5 THz, the absorption coefficient for the front illumination is A 1 = 45.2%.For the back illumination, the absorption coefficient is A 2 = 52.4%.At a frequency of 3.7 THz, the absorption coefficients for the front and back illuminations are A 1 = 44.1% and A 2 = 55.5%, respectively.We find that the 50% absorption bound is broken for the back illumination due to the proposed extremely asymmetric structural design.Furthermore, across a broad frequency range, the absorption for a single incident wave from either direction approaches 50%, which lays an important groundwork for realizing the broadband CPA response.

 
where int ( ) ra σ ω and int ( ) er σ ω represent the intra-and inter-band transition conductivities, e is the electron charge, kB is the Boltzmann constant,  is the Plank's constant, T is the Kelvin temperature, τ is the relaxation time, and EF is the Fermi level.In the low-THz range, intra-band transitions dominate the graphene response, and the graphene conductivity can be expressed by use of the Drude model [42]: Here, the relaxation time is assumed to be 10 −13 s.

Results of Broadband Absorber
Here, we initially investigate single illumination on the proposed asymmetric metasurface.The spectral responses, including reflection, transmission, and absorption, under normal incidence for both front and back illuminations are depicted in Figure 2a and Figure 2b, respectively.The Fermi level of graphene EF is 0.8 eV.It is clear that the transmission from both sides is identical because reciprocity is not broken with this designed configuration.However, reflection and absorption are different due to the structural asymmetry.As illustrated in Figure 2a,b, at a frequency of 2.5 THz, the absorption coefficient for the front illumination is A1 = 45.2%.For the back illumination, the absorption coefficient is A2 = 52.4%.At a frequency of 3.7 THz, the absorption coefficients for the front and back illuminations are A1 = 44.1% and A2 = 55.5%, respectively.We find that the 50% absorption bound is broken for the back illumination due to the proposed extremely asymmetric structural design.Furthermore, across a broad frequency range, the absorption for a single incident wave from either direction approaches 50%, which lays an important groundwork for realizing the broadband CPA response.Next, two counter-propagating TM-polarized incident waves illuminate the proposed design normally from the front and back sides simultaneously.Controlling phase difference is very important in achieving the CPA effect; since identical-amplitude coherent illumination was utilized in this work.By varying the phase difference θ between the two incident waves, we aimed to satisfy the CPA condition [34] and achieve perfect absorption.Figure 3a illustrates the calculated absorption coefficient [26] as a function of the operating frequency and phase difference θ between the two input waves.It can be found that over a wide range of frequencies, nearly perfect absorption is obtained when the phase difference θ is smaller than 30 • .Importantly, by tuning the phase difference between the incident waves, the absorption coefficients can be modulated from 0 to 99.7% at 2.5 THz and to 99.8% at 3.7 THz.This implies that a coherent beam traversing the CPA structure can be modulated from near-complete absorption to near-complete transparency.The modulation depth (MD) [4,30,43] is a convenient figure of merit for evaluating the performance of the proposed structure for use in devices such as modulators, switches, or detectors.It can be calculated as MD = (Amax − Amin)/Amax, where Amax and Amin are the maximum and minimum of absorption from phase modulation or Fermi level modulation.Here, the MD can reach as high as 0.997 at both 2.5 THz and 3.7 THz through phase modulation.These findings could enable a host of potential applications, including planar terahertz modulators and graphene-based coherent photodetectors [44].To effectively illustrate the excellent absorption capability of our design, we have set the phase difference constant to compute the relevant output coefficient across various frequencies.Here, we select a specific angle, θ = 0 • , as it aligns with the best absorption performance, marked by the destructive interference of the two incoming waves.The outcomes are graphically represented in Figure 3b, showcasing a distinct broadband CPA response with an impressive ≥90% absorption coefficient spanning a broad frequency spectrum from 1.65 THz to 4.49 THz.Note the same broadband CPA response can be achieved with this structure under TE-polarized incidence, since the structure is symmetric along the x-axis and y-axis.
sorption.Figure 3a illustrates the calculated absorption coefficient [26] as a function of the operating frequency and phase difference θ between the two input waves.It can be found that over a wide range of frequencies, nearly perfect absorption is obtained when the phase difference θ is smaller than 30°.Importantly, by tuning the phase difference between the incident waves, the absorption coefficients can be modulated from 0 to 99.7% at 2.5 THz and to 99.8% at 3.7 THz.This implies that a coherent beam traversing the CPA structure can be modulated from near-complete absorption to near-complete transparency.The modulation depth (MD) [4,30,43] is a convenient figure of merit for evaluating the performance of the proposed structure for use in devices such as modulators, switches, or detectors.It can be calculated as MD = (Amax − Amin)/Amax, where Amax and Amin are the maximum and minimum of absorption from phase modulation or Fermi level modulation.Here, the MD can reach as high as 0.997 at both 2.5 THz and 3.7 THz through phase modulation.These findings could enable a host of potential applications, including planar terahertz modulators and graphene-based coherent photodetectors [44].To effectively illustrate the excellent absorption capability of our design, we have set the phase difference constant to compute the relevant output coefficient across various frequencies.Here, we select a specific angle, θ = 0°, as it aligns with the best absorption performance, marked by the destructive interference of the two incoming waves.The outcomes are graphically represented in Figure 3b, showcasing a distinct broadband CPA response with an impressive ≥90% absorption coefficient spanning a broad frequency spectrum from 1.65 THz to 4.49 THz.Note that the same broadband CPA response can be achieved with this structure under TE-polarized incidence, since the structure is symmetric along the x-axis and y-axis.Our analysis suggests that the wideband CPA response arises due to the hybridization and coupling of narrowband resonant CPA responses from differently shaped graphene patch arrays on the front and back sides.To verify this analysis, we calculated the absorption coefficients for a structure featuring a single-sided graphene patch array.Figure 4a shows the calculated absorption of the structure with a square graphene patch array solely on the front side of the silicon substrate.It can be observed that a single resonance occurs at 2.34 THz.By optimizing the length l of the square graphene patches, the resonance shifts to 2.47 THz at l = 4.9 um, and the absorption coefficient approaches 1.With a solely circular graphene patch array, a single resonance persists at 3.77 THz, and the absorption coefficient reaches 70.5%, as shown in Figure 4b.We further optimized the structure by adjusting the radius r of the circular graphene patches.The absorption coefficient can reach 99.7% with r = 2.5 um, yet it remains a single resonance.It is implied that broadband response cannot be achieved without efficient coupling between two patch arrays on both sides of the substrate.To further illustrate the impact of structural asymmetry on the broadband absorption effect, here, we compare symmetric and asymmetric bifacial graphene metasurfaces.The result is shown in Figure 4c,d.Figure 4c illustrates the absorption coefficient for a design with square graphene patch arrays on both sides of the substrate.The red line represents the absorption coefficient with a patch length of 5 µm (a value used for the square patch in the asymmetric design), while the blue line shows the result with an optimized radius of 4 µm for the symmetric design.This figure indicates that identical square patch arrays on both sides fail to achieve a broadband CPA effect, even with further optimization.Figure 4d displays the absorption coefficient for a design with identical circular graphene patch arrays on both sides of the substrate.The blue line, symbolizing the optimized symmetric structure, demonstrates that the symmetric design fails to meet the criteria for broadband coherent absorption.These results indicate that significant structural asymmetry is essential to achieving a wideband CPA response, given that a graphene metasurface design exhibiting minimal asymmetry will exhibit CPA only at a single resonant frequency.

Results of Tunable Absorber
Here, we investigate the impact of different geometric parameters on the absorptivity of the asymmetric graphene metasurface.It is imperative to examine the effects of the silicon thickness t, the radius r of the circular patches, and the length l of the square

Results of Tunable Absorber
Here, we investigate the impact of different geometric parameters on the absorptivity of the asymmetric graphene metasurface.It is imperative to examine the effects of the silicon thickness t, the radius r of the circular patches, and the length l of the square patches on absorption characteristics.The corresponding simulated results are depicted in Figure 5a-c.Figure 5a illustrates the absorption spectra generated by varying the silicon thickness with the parameters held constant at r = 1.5 um and l = 5 um.As the silicon thickness increases, the first resonance at the lower frequency remains virtually unchanged, yet there is a slight decrease in absorption amplitude.Conversely, the second resonance at the higher frequency exhibits a blueshift, accompanied by a significant decrease in the absorption coefficient.Optimal absorption bandwidth is attained at a silicon thickness of 1.5 um.Subsequently, the absorption at both resonances initially increases and subsequently decreases as we adjust the radius of the circular patches, starting from 1.1 um, with the silicon thickness and square patches length fixed at t = 1.5 um and l = 5 um, respectively, as shown in Figure 5b. Figure 5c presents the absorption spectra for varying square patch lengths with t = 1.5 um and r = 1.5 um.Increasing the length l causes the first resonance to redshift, enhancing the maximum absorption coefficient.The structure demonstrates optimal performance when the length reaches 5 um.When the square graphene patch entirely covers the silicon substrate, namely, when l equals 5.2 um, a single resonance is observed, as indicated by the light cyan line in Figure 5c, which resembles the performance of a single-side circular graphene patch.This result implies that a complete graphene sheet does not resonate with the excited electromagnetic waves.It is worth noting that in practical applications, small relative shifts between the centers of the rectangular and circular graphene components are inevitable during fabrication.We set out to clarify this issue by performing more simulations, and it turns out that the relative shift shows a limited effect on the broadband absorption response.the performance of a single-side circular graphene patch.This result implies that a complete graphene sheet does not resonate with the excited electromagnetic waves.It is worth noting that in practical applications, small relative shifts between the centers of the rectangular and circular graphene components are inevitable during fabrication.We set out to clarify this issue by performing more simulations, and it turns out that the relative shift shows a limited effect on the broadband absorption response.Tunability is a unique property of graphene-based structures which enables dynamic control of the coherent absorption in the proposed device by adjusting the graphene's Fermi level E F .To demonstrate the tunability of the proposed CPA design, we calculated the output coefficients by sweeping E F from 0.1 eV to 1 eV.The results are shown in Figure 6.It can be observed that the CPA response can be achieved across a broad frequency range when the E F is greater than 0.3 eV, and as the E F increases, a blueshift in the resonant frequency is observed.This phenomenon can be attributed to graphene exhibiting more metallic behavior with the elevation of its Fermi level, which aligns with the established correlation between the resonant frequency and the Fermi level in graphene patches, which can be concisely represented as f r ∝ √ E F /L [21,30,45], with L denoting the side length of the graphene patch or the structural geometry.Moreover, we computed the MD by manipulating the graphene Fermi level, with the MD reaching a maximum of 0.88 at 2.5 THz and 0.94 at 3.7 THz.These results suggest that the graphene Fermi level is critical to the tunability of the absorption response, indicating that the structure functions as an efficient light modulator.In practical applications, wide-angle absorption represents a crucial fea sorbers.We investigated the impact of various incident angles on the absorpt mance.Figure 7 illustrates the coherent absorption of the proposed asymmetri metasurface as a function of the operating frequency and incident angle φ.C absorption magnitude and bandwidth are well maintained until the inciden creases to 45°, indicating that the absorption resonance is insensitive to small i gles.When the incident angle reaches 60°, the absorption sharply decreases weak coupling between the structure and the oblique incident waves.In practical applications, wide-angle absorption represents a crucial feature of absorbers.We investigated the impact of various incident angles on the absorption performance.Figure 7 illustrates the coherent absorption of the proposed asymmetric graphene metasurface as a function of the operating frequency and incident angle ϕ.Clearly, the absorption magnitude and bandwidth are well maintained until the incident angle increases to 45 • , indicating that the absorption resonance is insensitive to small incident angles.When the incident angle reaches 60 • , the absorption sharply decreases, implying weak coupling between the structure and the oblique incident waves.
To fully appreciate the advantages of this new device, it is essential to analyze its performance in comparison to previous reports on broadband absorption.Our comparison is shown in Table 1.Previous reports on broadband absorption have shown promising results, but they achieve enhanced absorption performance at the expense of structural simplicity.Our new device's ultrathin and wideband characteristics make it a more appealing option for on-chip THz communication technologies, as illustrated in Table 1.This will help in understanding the potential impact and value of this technology in the context of on-chip technologies.It can be found that this ultrathin, wideband CPA device has demonstrated its potential to outperform existing solutions in terms of absorption bandwidth and efficiency.
mance. Figure 7 illustrates the coherent absorption of the proposed asymmetric gr metasurface as a function of the operating frequency and incident angle φ.Clea absorption magnitude and bandwidth are well maintained until the incident a creases to 45°, indicating that the absorption resonance is insensitive to small incid gles.When the incident angle reaches 60°, the absorption sharply decreases, im weak coupling between the structure and the oblique incident waves.To fully appreciate the advantages of this new device, it is essential to ana performance in comparison to previous reports on broadband absorption.Our c son is shown in Table 1.Previous reports on broadband absorption have shown pr results, but they achieve enhanced absorption performance at the expense of st simplicity.Our new device's ultrathin and wideband characteristics make it a m pealing option for on-chip THz communication technologies, as illustrated in Tabl will help in understanding the potential impact and value of this technology in the of on-chip technologies.It can be found that this ultrathin, wideband CPA dev demonstrated its potential to outperform existing solutions in terms of absorptio width and efficiency.
Table 1.Comparison of the proposed device's broadband absorption with that of devices f previous reports.

Conclusions
In conclusion, the research on CPA using graphene metasurfaces has made significant progress in attaining ultrabroadband absorption capabilities within the THz spectrum.By harnessing the distinctive attributes of graphene and investigating its patterned configurations, including the asymmetric arrangement that merges square and circular graphene patches, we have succeeded in boosting absorption efficiency while preserving structural simplicity.The adjustability of the CPA system via the manipulation of the graphene Fermi level opens up a plethora of opportunities for on-chip THz communication technologies.The remarkable 90% absorption bandwidth of 2.84 THz and the MD value reaching as high as 0.997 demonstrate the potential of this advanced methodology for the development of sophisticated devices such as light modulators, photodetectors, and more.In summary, this study underscores the bright prospects of graphene-based metasurfaces in realizing efficient and tunable CPA over a broad frequency spectrum.
− .In order to achieve CPA, we need to optimize the proposed structure to meet the CPA conditions.

Figure 1 .
Figure 1.(a) Schematic for a unite cell of an asymmetric graphene metasurface illuminated by two counter-propagating and coherently modulated input waves, Ein1 and Ein2.(b) Schematic of an ion gel-gating approach used to dope the graphene patch array.

Figure 1 .
Figure 1.(a) Schematic for a unite cell of an asymmetric graphene metasurface illuminated by two counter-propagating and coherently modulated input waves, E in1 and E in2 .(b) Schematic of an ion gel-gating approach used to dope the graphene patch array.

Figure 2 .
Figure 2. The calculated reflection, transmission, and absorption of the proposed asymmetric metasurface when excited by a normal incident wave with (a) front-side illumination and (b) backside illumination.

Figure 2 .
Figure 2. The calculated reflection, transmission, and absorption of the proposed asymmetric metasurface when excited by a normal incident wave with (a) front-side illumination and (b) backside illumination.

Figure 3 .
Figure 3. (a) The computed absorption coefficient of the proposed asymmetric metasurface design as a function of the operating frequency and phase difference between the two incident waves.(b) The computed absorption coefficient as a function of frequency a with fixed phase difference, θ = 0 • .Maximum absorption coefficients are obtained at frequencies of 2.5 THz and 3.7 THz.The ≥90% absorption coefficient spans a broad frequency spectrum ranging from 1.65 THz to 4.49 THz.

Figure 4 .
Figure 4. (a,b) The computed coherent absorption coefficients of the graphene metasurface design with (a) a square graphene patch array only on the front side of the silicon substrate and (b) a circular graphene patch array only on the back side of the silicon substrate as a function of operating frequency.The blue line presents the optimized results for the corresponding structural design.(c,d) The computed absorption coefficients of the symmetric graphene metasurface design with (c) an identical square graphene patch array on both sides of the silicon substrate and (d) an identical circular graphene patch array on both sides of the silicon substrate as a function of operating frequency.The black line in (c,d) depicts the absorption coefficient of the asymmetric graphene metasurface design.

Figure 4 .
Figure 4. (a,b) The computed coherent absorption coefficients of the graphene metasurface design with (a) a square graphene patch array only on the front side of the silicon substrate and (b) a circular graphene patch array only on the back side of the silicon substrate as a function of operating frequency.The blue line presents the optimized results for the corresponding structural design.(c,d) The computed absorption coefficients of the symmetric graphene metasurface design with (c) an identical square graphene patch array on both sides of the silicon substrate and (d) an identical circular graphene patch array on both sides of the silicon substrate as a function of operating frequency.The black line in (c,d) depicts the absorption coefficient of the asymmetric graphene metasurface design.

Figure 5 .
Figure 5.The effect of the geometric parameters of the proposed asymmetric design on the absorption performance: (a) the silicon thickness t, (b) the radius r of the circular graphene patches, and (c) the length l of the square graphene patches.

Figure 5 .
Figure 5.The effect of the geometric parameters of the proposed asymmetric design on the absorption performance: (a) the silicon thickness t, (b) the radius r of the circular graphene patches, and (c) the length l of the square graphene patches.

Figure 6 .
Figure 6.The calculated coherent absorption of the proposed design as a function of the operating frequency and Fermi level of graphene.

Figure 7 .
Figure 7.The calculated coherent absorption of the proposed design as a function of the o frequency and incident angles.

Figure 7 .
Figure 7.The calculated coherent absorption of the proposed design as a function of the operating frequency and incident angles.

Table 1 .
Comparison of the proposed device's broadband absorption with that of devices from previous reports.