Next Article in Journal
Semantic Mapping for Autonomous Subsea Intervention
Next Article in Special Issue
Metasurface Source Antenna Gain Improvement Using Simple Side Metal Structure
Previous Article in Journal
A Novel Clutter Suppression Method Based on Sparse Bayesian Learning for Airborne Passive Bistatic Radar with Contaminated Reference Signal
Previous Article in Special Issue
A Study of Terahertz-Wave Cylindrical Super-Oscillatory Lens for Industrial Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical and Numerical Analysis of Active Switching for Narrow-Band Thermal Emission with Graphene Ribbon Metasurface

1
Graduate School of Science and Technology, Niigata University, 8050, Ikarashi 2-no-cho, Niigata 950-2181, Japan
2
Department of Mechanical and Production Engineering, Niigata University, 8050, Ikarashi 2-no-cho, Niigata 950-2181, Japan
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(20), 6738; https://doi.org/10.3390/s21206738
Submission received: 17 August 2021 / Revised: 1 October 2021 / Accepted: 6 October 2021 / Published: 11 October 2021
(This article belongs to the Collection Advances in Metamaterials or Plasmonics-Based Sensors)

Abstract

:
Components smaller than the wavelength of electromagnetic waves are called meta-atoms. Thermal emission can be controlled by an artificial structure in which these meta-atoms are arranged on the surface. This artificial structure is called a metasurface, and its optical properties are determined by the materials and shapes of the meta-atoms. However, optical devices may require active control of thermal emission. In the present study, we theoretically and numerically analyze a wavelength-selective emitter using a graphene ribbon metasurface. The graphene ribbon metasurface consists of a graphene ribbon array, potassium bromide thin film, and silver substrate. The geometric parameters of the graphene metasurface are determined based on an equivalent circuit model that agrees well with the results of the electromagnetic field analysis (rigorous coupled-wave analysis). The proposed emitter causes impedance matching depending on the conductivity of the graphene ribbon in a very narrow wavelength range. The conductivity of graphene can be actively controlled by the gate voltage. Therefore, the proposed emitters may realize near-perfect emission with a high quality factor and active controllable switching for various wavelengths. In addition, the quality factor can be changed by adjusting the electron mobility of graphene. The proposed emitter can be used for optical devices such as thermophotovoltaic systems and biosensing.

1. Introduction

Thermal emission is a spontaneous and continuous photon emission from the thermal reservoir [1]. Therefore, it was difficult to realize active switching (i.e., on/off switching) of thermal emission without mechanical shutter. In addition, thermal emission is a broadband light because it originates from the fluctuating current in a material; however, narrow-band thermal emission is a promising technology for the improvement of thermal devices [2].
A meta-atom is a component smaller than the wavelength of an electromagnetic wave. An artificial material that realizes functions that are difficult with natural materials by arranging meta-atoms on the surface is called a metasurface. In recent years, nano/micro-scale artificial materials called metasurfaces were used to control thermal emission. Controlling light with metasurfaces has the possibility to realize active switching of narrow-band thermal emission. Active switching of narrow-band thermal emission is a key technology with applications in thermophotovoltaic systems [3,4], infrared heaters [5], biosensing [6,7,8], microbolometers [9,10], imaging [11], and optical communications [12].
Normally, the optical properties of metasurfaces may be determined by their geometrical shape [13,14,15]. Therefore, thermal emission cannot be actively controlled by conventional metasurfaces. To solve this problem, we focus on graphene, which is one of the carbon allotropes, and has a honeycomb structure where carbon atoms are bonded two-dimensionally [16]. One of the characteristics of graphene is that its electrical conductivity changes when a gate voltage is applied. Thus, active switching of thermal emission may be possible without changing the shape of the structure using graphene metasurfaces [17,18].
In the present study, we propose a graphene ribbon metasurface to design a device for realizing narrow-band emission peak at the target wavelength and active wavelength control. This structure can be designed based on an equivalent circuit model. We computationally demonstrate that the designed graphene ribbon metasurface can exhibit near-perfect narrow-band thermal emission under active switching using electromagnetic wave analysis. In addition, the peak wavelength and intensity of thermal emission can be controlled by adjusting the Fermi energy of graphene.

2. Materials and Methods

Figure 1a shows the schematic of the proposed graphene ribbon metasurface emitter. It consists of a graphene ribbon array, a dielectric thin film, and a metallic substrate. Such ribbon-shaped graphene structures can be feasible via top-down structuring such as plasma CVD [19]. The graphene ribbon is periodic in the x-direction and extended infinitely in the y-direction. The metallic substrate was sufficiently thick. This structure has period Λ [µm], graphene ribbon width W [µm], and dielectric thickness d [µm]. The optical conductivity of graphene σg [S] is given as [20]
σ g = ln 2 | E F | ( ω i τ 1 ) ħ 2 | E F | + ( ω i τ 1 ) ħ + i e 2 k B T π ħ 2 ( ω i τ 1 ) { E F k B T + 2 ln ( e E F / k B T + 1 ) } ,
where EF [eV] is the Fermi energy of graphene, ω [rad/s] is the angular frequency in vacuum, τ [ps] is the relaxation time of graphene, ħ [m2kg/s] is the reduced Planck constant, e [C] is the electron charge, kB [m2kg/s2K] is the Boltzmann constant, and T [K] is the temperature of graphene. The Fermi energy of graphene can be tuned by applying a gate voltage.
The permittivity of graphene εg is given as [21]
ε g = i σ g ε 0 ω Δ ,
where ε0 [F/m] is the permittivity of vacuum and [nm] is the thickness of the graphene.
The equivalent circuit model is useful for designing graphene ribbon metasurfaces and elucidating the underlying mechanism of emissivity enhancement [22]. Figure 1b shows the equivalent circuit model corresponding to Figure 1a. Y0 [S] is the admittance of free space, Yin [S] is the admittance of graphene ribbon metasurface, Yd [S] is the admittance of the dielectric layer, Ym [S] is the admittance of the metallic substrate, and Y m tr [S] is the admittance of structure below graphene ribbon. R1 [Ω], L1 [H], and C1 [F] are the resistance, inductance, and capacitance, respectively, corresponding to the first-order mode of graphene plasmon. In this equivalent circuit model, the influence of higher-order modes of graphene plasmon is negligible compared to its first-order mode [23]. The mechanism of emissivity enhancement in this equivalent circuit model is based on the impedance matching theory. A reflected wave is canceled under the impedance matching condition between the graphene ribbon metasurface and vacuum, and thus, nearly perfect emission can be achieved. Since admittance is the reciprocal of impedance, a nearly perfect emission can be obtained when the admittance of the graphene ribbon metasurface is equal to that of the vacuum ( Y in = Y 0 ).
Here, Yin and Y0 are given as follows:
Y in = Y m tr + Y G ,
Y G = 1 R 1 + i ω L 1 + 1 i ω C 1 ,
Y m tr = Y d Y m + i Y d tan ( ω n d / c 0 ) Y d + i Y m tan ( ω n d / c 0 ) ,
Y 0 = 1 / η 0 ,
where YG [S] is the admittance of graphene, n is the refractive index of the dielectric layer, c0 [m/s] is the speed of light in vacuum, and η0 [Ω] is the impedance of the vacuum (=120π). Each coefficient representing the first-order mode of graphene plasmon is given by the following equations in the wavelength range, where the first term on the right side of Equation (1) may be negligible.
R 1 = Λ S 1 2 π ħ 2 e 2 E F τ ,
L 1 = Λ S 1 2 π ħ 2 e 2 E F = τ R 1 ,
C 1 = S 1 2 Λ 2 W ε e f f π r 1 ,
S 1 2 8 9 W ,
where εeff is the average relative permittivity of the upper and lower media of graphene, and in this structure ε e f f = ε 0 ( 1 + n 2 ) / 2 . r1 is a function of the fill factor and takes a value corresponding to W/Λ [22]. It is derived from the surface current density in the graphene plasmon first-order mode obtained using perturbation theory. In this case, the non-perturbative term is the surface current density when light is incident on a single graphene ribbon. The perturbation term is the effect of adjacent ribbons on the surface current density when the ribbons are arranged in a periodic array.
Thus, the conditional expression for admittance matching for a target wavelength λt [µm] is given as
R 1 = η 0 ,
L 1 C 1 = 1 / ( 2 π f t ) 2 ,
ω t n d / c 0 = π / 2 ,
Y m ,
Here, f t = c 0 / λ t [Hz] is the target frequency, and ω t = 2 π f t [rad/s] is the target angle frequency. The absorptivity predicted by the equivalent circuit model used in this design is given as:
A ( λ ) = 1 | Y in ( λ ) Y 0 Y in ( λ ) + Y 0 | 2 .
Since the silver (Ag) substrate is opaque, spectral emissivity (ελ) can be calculated from the Kirchhoff’s law, i.e., ε λ = A λ .
Structural admittance α = Re ( Y in ) η 0 and β = n R 1 2 / 8 η 0 f t L 1 are defined for the real and imaginary parts, respectively. These parameters are useful for determining the bandwidth of the emitter. α indicates the value of the real part of the structural admittance at the target wavelength. When β is close to unity, the bandwidth becomes narrower and the emissivity increases. As α increases from unity, the bandwidth becomes broader and the emissivity decreases. β indicates the slope of the imaginary part of the structural admittance at the target wavelength. When β is closer to zero, the bandwidth becomes narrower. As β increases from zero, the band becomes broader. Therefore, α = 1 and β 0 are appropriate for the optimized design of a narrow-band thermal emitter.
The undetermined structural parameters are calculated using the above equations. According to Y in = Y 0 , Equations (6) and (11) and the definition of α,
Re ( Y in ) = 1 / R 1 = α / η 0 .
Rearranging Equation (13) gives the form:
d = c 0 4 n f t .
If it is desirable that the slope of the structural admittance imaginary part is 0,
d d λ Im ( Y in ) = 0 .
The admittance imaginary part below graphene and the admittance imaginary part of graphene are expressed from Equations (5), (12)–(14) and (17) in the following.
Im ( Y m tr ) | f f t π n 2 η 0 ( f f t 1 ) ,
Im ( Y G ) | f f t 2 π L 1 R 1 2 f ( f 2 f t 2 ) .
Differentiating Equations (19) and (20) and substituting them into Equations (3) and (18), lead to
π n 2 η 0 f t 4 π L 1 R 1 2 = 0 .
Define β as follows:
β = n R 1 2 8 η 0 f t L 1 .
Substituting Equations (8) and (16) into the Equation (22), lead to
π = n 8 α β f t .
According to the relation between relaxation time and Fermi energy in graphene,
E F = π e v F 2 / μ .
where v F = 10 6 [m/s] is the Fermi velocity of graphene, µ [m2/Vs] is the electron mobility of graphene. According to Equations (7) and (10),
W Λ = 9 π ħ 2 8 e 2 E F τ R 1 .
According to Equations (8), (9) and (12),
W = r 1 e 2 E F 2 ħ 2 ω t 2 ε e f f .
Λ is determined by Equations (25) and (26). Moreover, there are two restrictions in determining the structural parameters: (1) the graphene ribbon width must not exceed the period; (2) because of the effectiveness of the equivalent circuit model, the period must not exceed the wavelength. Due to the above restrictions,
W < Λ s 1 ,
Λ < λ 0 n s 2 ,
where s1 and s2 are safety factor parameters, and s 1 = 0.9 and s 2 = 0.8 are employed for the calculation.
Based on the equivalent circuit model, narrow-band emitters at target wavelengths λt = 6, 8, and 10 µm were designed. A potassium bromide (KBr) layer was used as the dielectric layer, and an Ag substrate was used as the metallic substrate. Here, the phenomenon can be simplified if the refractive index of the dielectric layer is constant. In addition, unintended emission appears when the dielectric layer has the extinction coefficient. Therefore, KBr was used for the dielectric layer, which has a substantially constant refractive index and almost no extinction coefficient in the near-infrared region. The refractive index and dielectric function of KBr was obtained from the tabulated data from Palik’s databook [24]. The dielectric function of Ag was obtained using a Drude model [1]: The permittivity of graphene εAg is given as ε Ag = ε ω p 2 / ( ω ( ω + i γ ) ) with a high-frequency dielectric constant ε = 3.40, plasma frequency ω p = 1.39 × 10 16 rad/s, and scattering rate γ = 2.70 × 10 13 rad/s. The temperature of graphene was set to be T = 300 K, the electron mobility of graphene is µ = 2.0 m2/Vs, and the structural parameters are determined using the equivalent circuit model. The rigorous couple-wave analysis (RCWA) method [25], which is a semi-analytical method, is employed as an electromagnetic wave analysis. In the RCWA method, the structure is treated as a grid with a uniform dielectric constant distribution in the depth direction. First, the electromagnetic waves in each layer are expressed by Fourier expansion, and the general solution of the electromagnetic field that can exist in the layers is obtained. Next, the solution of the electromagnetic field in the entire region can be obtained by imposing continuous conditions of the electromagnetic field at the boundary of each layer. As a result, the distribution of electromagnetic fields and diffraction efficiency can be calculated. The diffraction order was set to 200, and only transverse magnetic waves were perpendicularly incident on the surface. Note that transverse electric waves were not considered. In the case of transverse magnetic waves, the incident wavevector does not have y-component. Therefore, the incident wavevector can be expressed as k inc = k x x ^ + k z z ^ = k 0 sin θ x ^ + k 0 cos θ z ^ , where k0 is the wavevector in vacuum and θ is the incident angle. The surface roughness was not considered in the RCWA calculations in this study. However, previous studies have shown that the effect of surface roughness is significant in plasmonic absorbers [26,27]. Although the present calculations dealt with ideal conditions, our future works can deal with this problem by considering periodic structures that approximate the surface roughness in the RCWA calculations. To evaluate the performance of the designed structures, the following quality factor was employed:
Q = f t f 2 f 1 ( f 1 < f t < f 2 ) ,
where f1 [Hz] and f2 [Hz] are the frequencies at which the emissivity is half the peak emission value on the low-frequency side and high-frequency side of the peak frequency, respectively.

3. Results and Discussion

Figure 2 represents the spectral normal emissivity calculated by the RCWA method (solid line) and the equivalent circuit model (dashed line). Table 1 shows the structural parameters designed for each target wavelength. As shown in Figure 2, the multiple wavelength-selective thermal emissions could be successfully obtained by the RCWA analysis and the equivalent circuit model, and their results agreed well at each target wavelength. Furthermore, the emissivity with EF = 0.01 eV is nearly zero in each test case, which implies that active thermal emission switching was also successfully achieved. Therefore, the equivalent circuit theory is significantly effective for designing the graphene ribbon metasurface.
Next, we investigated the relationship between the quality factor and the parameter β. Table 1 shows the quality factor of each emission spectrum in Figure 2. The quality factor increases as the target wavelength becomes shorter. Figure 3 shows the contour plot of the quality factor with varying target wavelengths and β. The white region indicates absence of data because the constraint conditions in Equations (27) and (28) are not satisfied. The quality factor depends only on β and is independent of the target wavelengths. Therefore, the design is performed using the smallest β among the β satisfying Equations (27) and (28). As the target wavelength increases, the minimum value of β increases due to the constraint of Equation (28). Therefore, the minimum value of β increases as the target wavelength increases, and the quality factor decreases.
In addition, the electric field distribution is shown to investigate the mechanism of the emissivity peak. Figure 4a is a contour diagram of electric field distribution and Figure 4b of x component of the electric field from the x-z plane around the graphene ribbon at the peak wavelength of 10 µm. Figure 5 shows that the electric field is enhanced at the center and edges of the graphene ribbon. Therefore, graphene plasmons are excited in the graphene ribbon [28]. Furthermore, this structure forms an asymmetric FP cavity in which a dielectric layer with a low refractive index is sandwiched between a graphene ribbon array and a metal substrate. In this structure, the excitation of graphene plasmons causes FP resonance [29]. Here, the admittance of the graphene ribbon in the equivalent circuit model is based on the first-order mode of graphene plasmon. Therefore, it supports the excitation of graphene plasmons in Figure 5. As a result, in the proposed structure, the excitation of graphene plasmons in the asymmetric FP cavity causes FP resonance and increases emissivity. The phenomenon can be explained as follows: The emissivity enhancement may be caused by graphene plasmon [22]. When the frequency of the incident photons matches the first localized mode of the graphene ribbons, some photons are absorbed, and others are reflected and transmitted. The transmitted photons are reflected by the metallic substrate, some of which change the phase and exit the structure, resulting in destructive interference. In addition, since the metallic substrate is sufficiently thick, transmission does not occur. Therefore, the incident waves are completely absorbed by the graphene ribbons. On the other hand, when EF = 0.01 eV, the normal emissivity is almost zero at any wavelength. This is because graphene plasmon cannot be excited in this situation. Hence, the active switching of emissivity becomes possible by adjusting the Fermi energy in the proposed structure.
Moreover, the dispersion relation is shown in Figure 5 to confirm that the asymmetric FP resonance is excited in the proposed structure. Figure 5 shows that the emission curve is divided into two peaks when the wavevector is increased. This phenomenon is consistent with the characteristics of asymmetric FP resonance [29]. In addition, the division of this emission peak is due to Rabi splitting analogues [29,30,31]. Rabi splitting analogues are a phenomenon in which when two energy levels are combined, a new split energy state appears due to the coupling between the two energy levels and the new eigenstate. The Fabry–Perot mode and surface plasmon polariton (SPP) mode of graphene ribbon are excited in the proposed structure. Therefore, Rabi splitting analogues occur when these modes are coupled. The horizontal mode indicates the FP mode, and the tilted mode indicates the SPP mode in Figure 5.
To investigate the effect of Fermi energy, normal emissivity spectra at several Fermi energies were calculated, as shown in Figure 6. In this calculation, the structural parameters of the graphene ribbon metasurface were the same as the test case in which the target wavelength was λt = 10 µm, and the Fermi energy of the graphene was 0.5–1.5 times EF = 1.37 eV. As the Fermi energy was reduced, the peak wavelength of the emissivity shifted to a longer wavelength. This is because the resonance wavelength of the graphene plasmon is red-shifted as the Fermi energy decreases. Moreover, as the Fermi energy decreases, the intensity of the emissivity peak decreases. This can be explained using the equivalent circuit model. As the Fermi energy decreases, the conductivity of graphene decreases. Thereby, the resistance and reactance of graphene, and the admittance of graphene, increase. As a result, emissivity decreases because impedance matching is less likely to occur. Since the proposed structure can change the peak wavelength of emissivity by adjusting the Fermi energy of graphene, it can also be used for wavelength control of the thermal emitter without changing the structure. It is important to note that if the Fermi energy is too high for the experiment, the substrate might be destroyed. In this study, we have determined the Fermi energy of graphene based on previous studies with numerical simulations [32,33,34], thus confirming its theoretical performance.
To investigate the effect of changes in the electron mobility of graphene on the structure, the emission spectra at several electron mobilities were calculated, as shown in Figure 7a. Table 2 shows the structural parameters and quality factors for each spectrum. These results show that the quality factor decreases as the electron mobility of graphene decreases. The reason for this can be understood from Figure 7b. Figure 7b shows the quality factor of the emissivity peak of the structure designed by determining the electron mobility of graphene and β. The region without data is the region that does not satisfy Equations (27) and (28). The quality factor depends only on β and hardly on the electron mobility of graphene (Figure 7b). Similar to the relationship between the target wavelength and the quality factor, the minimum value of β increases owing to the constraint of Equation (28) as the electron mobility of graphene increases. Therefore, β increases as the electron mobility of graphene increases, and the quality factor decreases. Since the electron mobility of graphene varies depending on the fabrication method, it is important to select the fabrication method according to the desired bandwidth.
In the proposed structure, there are many parameters that can affect the emissivity. Therefore, Figure 8 shows the emissivity spectrum when the period, graphene ribbon width, and dielectric thickness are changed. In this calculation, the structural parameters of the graphene ribbon metasurface were the same as the test case in which the target wavelength was λt = 10 µm. In Figure 8a, the period was 0.9–1.1 times Λ = 5.23 µm, in Figure 8b, the graphene ribbon width was 0.9–1.1 times W = 0.355 µm, and in Figure 8c, the dielectric thickness was 0.9–1.1 times d = 1.64 µm. Figure 8a,c show that the normal emissivity is nearly independent of the period and dielectric thickness, respectively. However, Figure 8b shows that the peak wavelength of the emissivity shifts to a longer wavelength as the graphene ribbon width is increased. This is because graphene plasmons are excited at the ends of the graphene ribbon (Figure 4), and the resonance wavelength shifts to the longer wavelength as the graphene ribbon width increases.

4. Conclusions

In summary, we proposed a graphene ribbon metasurface for active thermal emission switching. Near perfect thermal emission at the targeted wavelengths could be achieved by adjusting the structural parameters, and emissivity switching was achieved by adjusting the Fermi energy of graphene in those structures. In addition, the proposed structure can be easily designed using an equivalent circuit model based on the impedance matching theory. FP cavities are formed in the proposed structure, and graphene plasmons cause an asymmetric FP resonance. The wavelength and intensity of resonance can be changed by adjusting the Fermi energy of graphene. Furthermore, the quality factor may also be affected by the electron mobility of graphene. This study may be applied to the initial design of a structure for active control, and facilitate the development of graphene-based thermal devices.

Author Contributions

Investigation, K.Y., T.S., H.O. and A.S.; writing—original draft preparation, K.Y.; writing—review and editing, K.Y. and A.S.; supervision, A.S.; project administration, A.S.; funding acquisition, A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by JSPS KAKENHI Grant Numbers JP15K17985, JP18K03974.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Modest, M.F. Radiative Heat Transfer, 3rd ed.; Academic Press: Cambridge, MA, USA, 2013. [Google Scholar]
  2. Sakurai, A.; Yada, K.; Simomura, T.; Ju, S.; Kashiwagi, M.; Okada, H.; Nagao, T.; Tsuda, K.; Shiomi, J. Ultranarrow-Band Wavelength-Selective Thermal Emission with Aperiodic Multilayered Metamaterials Designed by Bayesian Optimization. ACS Cent. Sci. 2019, 5, 319–326. [Google Scholar] [CrossRef]
  3. Lenert, A.; Bierman, D.; Nam, Y.; Chan, W.; Celanovic, I.; Soljacic, M.; Wang, E. A nanophotonic solar thermophotovoltaic device. Nat. Nanotechnol. 2014, 9, 126–130. [Google Scholar] [CrossRef] [PubMed]
  4. Rephaeli, E.; Fan, S. Absorber and emitter for solar thermophotovoltaic systems to achieve efficiency exceeding the Shockley-Queisser limit. Opt. Express 2009, 17, 15145–15159. [Google Scholar] [CrossRef] [PubMed]
  5. Totani, T.; Sakurai, A.; Kondo, Y. A wavelength control emitter for drying furnace. In Proceedings of the Asian Conference on Thermal Sciences ACTS-P00423, Jeju, Korea, 27–29 March 2017. [Google Scholar]
  6. Liu, N.; Mesch, M.; Weiss, T.; Hentschel, M.; Giessen, H. Infrared Perfect Absorber and Its Application As Plasmonic Sensor. Nano Lett. 2010, 10, 2342–2348. [Google Scholar] [CrossRef] [PubMed]
  7. Wu, C.; Khanikaev, A.; Adato, R.; Arju, N.; Yanik, A.; Altug, H.; Shvets, G. Fano-resonant asymmetric metamaterials for ultrasensitive spectroscopy and identification of molecular monolayers. Nat. Mater. 2012, 11, 69–75. [Google Scholar] [CrossRef]
  8. Luo, S.; Zhao, J.; Zuo, D.; Wang, X. Perfect narrow band absorber for sensing applications. Opt. Express 2016, 24, 9288–9294. [Google Scholar] [CrossRef]
  9. Liu, X.; Wang, L.; Zhang, Z. Wideband Tunable Omnidirectional Infrared Absorbers Based on Doped-Silicon Nanowire Arrays. J. Heat Transf.-Trans. ASME 2013, 135, 061602. [Google Scholar] [CrossRef]
  10. Du, K.; Li, Q.; Zhang, W.; Yang, Y.; Qiu, M. Wavelength and Thermal Distribution Selectable Microbolometers Based on Metamaterial Absorbers. IEEE Photonics J. 2015, 7, 1–8. [Google Scholar] [CrossRef] [Green Version]
  11. Landy, N.; Bingham, C.; Tyler, T.; Jokerst, N.; Smith, D.; Padilla, W. Design, theory, and measurement of a polarization-insensitive absorber for terahertz imaging. Phys. Rev. B 2009, 79, 125104. [Google Scholar] [CrossRef] [Green Version]
  12. Miyoshi, Y.; Fukazawa, Y.; Amasaka, Y.; Reckmann, R.; Yokoi, T.; Ishida, K.; Kawahara, K.; Ago, H.; Maki, H. High-speed and on-chip graphene blackbody emitters for optical communications by remote heat transfer. Nat. Commun. 2018, 9, 1279. [Google Scholar] [CrossRef] [Green Version]
  13. Wang, L.; Zhang, Z. Phonon-mediated magnetic polaritons in the infrared region. Opt. Express 2011, 19, A126–A135. [Google Scholar] [CrossRef] [PubMed]
  14. Smaali, R.; Omeis, F.; Moreau, A.; Taliercio, T.; Centeno, E. A universal design to realize a tunable perfect absorber from infrared to microwaves. Sci. Rep. 2016, 6, 32589. [Google Scholar] [CrossRef] [PubMed]
  15. Matsuno, Y.; Sakurai, A. Electromagnetic resonances of wavelength-selective solar absorbers with film-coupled fishnet gratings. Opt. Commun. 2017, 385, 118–123. [Google Scholar] [CrossRef]
  16. Geim, A.; Novoselov, K. The rise of graphene. Nat. Mater. 2007, 6, 183–191. [Google Scholar] [CrossRef] [PubMed]
  17. Cai, Y.; Zhu, J.; Liu, Q.; Lin, T.; Zhou, J.; Ye, L.; Cai, Z. Enhanced spatial near-infrared modulation of graphene-loaded perfect absorbers using plasmonic nanoslits. Opt. Express 2015, 23, 32318–32328. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Zhao, B.; Zhang, Z. Strong Plasmonic Coupling between Graphene Ribbon Array and Metal Gratings. Acs Photonics 2015, 2, 1611–1618. [Google Scholar] [CrossRef]
  19. Suzuki, H.; Kaneko, T.; Shibuta, Y.; Ohno, M.; Maekawa, Y.; Kato, T. Wafer-scale fabrication and growth dynamics of suspended graphene nanoribbon arrays. Nat. Commun. 2016, 7, 11797. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Hanson, G. Dyadic Green’s functions and guided surface waves for a surface conductivity model of graphene. J. Appl. Phys. 2008, 103, 064302. [Google Scholar] [CrossRef] [Green Version]
  21. Vakil, A.; Engheta, N. Transformation Optics Using Graphene. Science 2011, 332, 1291–1294. [Google Scholar] [CrossRef] [Green Version]
  22. Khavasi, A. Design of ultra-broadband graphene absorber using circuit theory. J. Opt. Soc. Am. B-Opt. Phys. 2015, 32, 1941–1946. [Google Scholar] [CrossRef]
  23. Khavasi, A.; Rejaei, B. Analytical Modeling of Graphene Ribbons as Optical Circuit Elements. IEEE J. Quantum Electron. 2014, 50, 397–403. [Google Scholar] [CrossRef]
  24. Palik, E.D. Handbook of Optical Constants of Solids II; Academic Press: Cambridge, MA, USA, 1991. [Google Scholar]
  25. Georgia Tech Nanoscale Thermal Radiation Laboratory. Rigorous Coupled-Wave Analysis (RCWA). Available online: http://zhang-nano.gatech.edu. (accessed on 10 May 2021).
  26. Nishijima, Y.; Balcytis, A.; Naganuma, S.; Seniutinas, G.; Juodkazis, S. Tailoring metal and insulator contributions in plasmonic perfect absorber metasurfaces. ACS Appl. Nano Mater. 2018, 1, 3557–3564. [Google Scholar] [CrossRef]
  27. To, N.; Juodkazis, S.; Nishijima, Y. Detailed Experiment-Theory Comparison of Mid-Infrared Metasurface Perfect Absorbers. Micromachines 2020, 11, 409. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Hou, H.; Teng, J.; Palacios, T.; Chua, S. Edge plasmons and cut-off behavior of graphene nano-ribbon waveguides. Opt. Commun. 2016, 370, 226–230. [Google Scholar] [CrossRef]
  29. Chen, D.; Yang, J.; Zhang, J.; Huang, J.; Zhang, Z. Section 1 Tunable broadband terahertz absorbers based on multiple layers of graphene ribbons. Sci. Rep. 2017, 7, 15836. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Zhu, Y.; Gauthier, D.J.; Morin, S.; Wu, Q.; Carmichael, H.; Mossberg, T. Vacuum Rabi splitting as a feature of linear-dispersion theory: Analysis and experimental observations. Phys. Rev. Lett. 1990, 64, 2499. [Google Scholar] [CrossRef] [PubMed]
  31. Liu, F.; Cubukcu, E. Tunable omnidirectional strong light-matter interactions mediated by graphene surface plasmons. Phys. Rev. B 2013, 88, 115439. [Google Scholar] [CrossRef] [Green Version]
  32. Wan, Y.; Deng, L.; Wang, L.; Yang, M.; Wang, Y. Modulation of Visible and Near-Infrared Surface Plasmon Resonance of Au Nanoparticles Based on Highly Doped Graphene. Plasmonics 2017, 12, 1317–1324. [Google Scholar] [CrossRef]
  33. Baqir, M.; Choudhury, P.; Fatima, T.; Ibrahim, A.-B. Graphene-over-graphite-based metamaterial structure as optical filter in the visible regime. Optik 2019, 180, 832–839. [Google Scholar] [CrossRef]
  34. Zeng, L.; Huang, T.; Liu, G.-B.; Zhang, H.-F. A tunable ultra-broadband linear-to-circular polarization converter containing the graphene. Opt. Commun. 2019, 436, 7–13. [Google Scholar] [CrossRef]
Figure 1. (a) Schematic of graphene ribbon metasurface. This structure consists of Ag substrate, KBr layer, and graphene ribbon layer. The period is Λ [µm], graphene ribbon width is W [µm], and dielectric thickness is d [µm]. Transverse magnetic wave is used as the incident wave. (b) Equivalent circuit model based on graphene plasmon.
Figure 1. (a) Schematic of graphene ribbon metasurface. This structure consists of Ag substrate, KBr layer, and graphene ribbon layer. The period is Λ [µm], graphene ribbon width is W [µm], and dielectric thickness is d [µm]. Transverse magnetic wave is used as the incident wave. (b) Equivalent circuit model based on graphene plasmon.
Sensors 21 06738 g001
Figure 2. Normal emissivity spectrum for RCWA (colored solid line) and the equivalent circuit model circuit model (colored dashed line) of (a) switch-on state and (b) switch-off state. Target wavelengths are 6, 8, and 10 µm, and their Fermi energies of switch-on state are 1.64, 1.48, and 1.37 eV, respectively. Their Fermi energies of switch-off state are 0.01 eV.
Figure 2. Normal emissivity spectrum for RCWA (colored solid line) and the equivalent circuit model circuit model (colored dashed line) of (a) switch-on state and (b) switch-off state. Target wavelengths are 6, 8, and 10 µm, and their Fermi energies of switch-on state are 1.64, 1.48, and 1.37 eV, respectively. Their Fermi energies of switch-off state are 0.01 eV.
Sensors 21 06738 g002
Figure 3. Quality factor for different target wavelengths, and bandwidth parameter β. Normal emissivity is calculated from the equivalent circuit model to obtain the quality factor.
Figure 3. Quality factor for different target wavelengths, and bandwidth parameter β. Normal emissivity is calculated from the equivalent circuit model to obtain the quality factor.
Sensors 21 06738 g003
Figure 4. Contour plots of (a) Electric field intensity and (b) x-component of electric field contours around graphene ribbon at the target wavelength λt = 10 µm.
Figure 4. Contour plots of (a) Electric field intensity and (b) x-component of electric field contours around graphene ribbon at the target wavelength λt = 10 µm.
Sensors 21 06738 g004
Figure 5. Emissivity contours in terms of the wavevector and x-component of wavevector for the graphene ribbon metasurface. The horizontal mode is the FP mode, and the tilted mode is the SPP mode.
Figure 5. Emissivity contours in terms of the wavevector and x-component of wavevector for the graphene ribbon metasurface. The horizontal mode is the FP mode, and the tilted mode is the SPP mode.
Sensors 21 06738 g005
Figure 6. Normal emissivity spectrum for graphene ribbon metasurface with different Fermi energies of graphene. The target wavelength is 10 µm, and the Fermi energy of graphene is 0.5–1.5 times of EF = 1.37 eV.
Figure 6. Normal emissivity spectrum for graphene ribbon metasurface with different Fermi energies of graphene. The target wavelength is 10 µm, and the Fermi energy of graphene is 0.5–1.5 times of EF = 1.37 eV.
Sensors 21 06738 g006
Figure 7. (a) Normal emissivity spectrum for graphene ribbon metasurface with different electron mobilities of graphene. The target wavelength is 10 µm, and the electron mobilities of graphene are µ = 0.10, 0.20, 0.60, 2.0, and 6.0 m2/Vs. (b) Quality factor for different mobilities of graphene, and bandwidth parameter β. Normal emissivity is calculated from RLC circuit model to obtain the quality factor.
Figure 7. (a) Normal emissivity spectrum for graphene ribbon metasurface with different electron mobilities of graphene. The target wavelength is 10 µm, and the electron mobilities of graphene are µ = 0.10, 0.20, 0.60, 2.0, and 6.0 m2/Vs. (b) Quality factor for different mobilities of graphene, and bandwidth parameter β. Normal emissivity is calculated from RLC circuit model to obtain the quality factor.
Sensors 21 06738 g007
Figure 8. Normal emissivity spectrum for graphene ribbon metasurface with different period, graphene ribbon width, and dielectric thickness. (a) The target wavelength is 10 µm, and the period is 0.9–1.1 times of Λ = 5.23 µm. (b) The target wavelength is 10 µm, and the graphene ribbon width is 0.9–1.1 times of W = 0.355 µm. (c) The target wavelength is 10 µm, and the dielectric thickness is 0.9–1.1 times of d = 1.64 µm.
Figure 8. Normal emissivity spectrum for graphene ribbon metasurface with different period, graphene ribbon width, and dielectric thickness. (a) The target wavelength is 10 µm, and the period is 0.9–1.1 times of Λ = 5.23 µm. (b) The target wavelength is 10 µm, and the graphene ribbon width is 0.9–1.1 times of W = 0.355 µm. (c) The target wavelength is 10 µm, and the dielectric thickness is 0.9–1.1 times of d = 1.64 µm.
Sensors 21 06738 g008aSensors 21 06738 g008b
Table 1. Design parameters and obtained quality factors of the graphene ribbon metasurface.
Table 1. Design parameters and obtained quality factors of the graphene ribbon metasurface.
Target Wavelengthλt [µm]6810
PeriodΛ [µm]2.854.015.23
Width of graphene ribbonW [µm]0.1340.2310.355
Thickness of dielectricd [µm]0.8911.251.64
Relaxation time of grapheneτ [ps]0.3290.2960.273
Fermi energy of grapheneEF [eV]1.641.481.37
Quality factorQ56.040.025.9
Table 2. Design parameters and obtained quality factors of graphene ribbon metasurface used for the calculation in Figure 7a.
Table 2. Design parameters and obtained quality factors of graphene ribbon metasurface used for the calculation in Figure 7a.
Mobility of Grapheneµ [m2Vs]0.100.200.602.06.0
Target wavelengthλt [µm]10
PeriodΛ [µm]5.245.235.245.245.23
Width of graphene ribbonW [µm]0.4470.3550.2460.1650.114
Thickness of dielectricd [µm]1.64
Relaxation time of grapheneτ [ps]0.1720.2370.5691.272.64
Fermi energy of grapheneEF [eV]1.721.370.9480.6350.440
Quality factorQ15.525.952.2114220
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yada, K.; Shimojo, T.; Okada, H.; Sakurai, A. Theoretical and Numerical Analysis of Active Switching for Narrow-Band Thermal Emission with Graphene Ribbon Metasurface. Sensors 2021, 21, 6738. https://doi.org/10.3390/s21206738

AMA Style

Yada K, Shimojo T, Okada H, Sakurai A. Theoretical and Numerical Analysis of Active Switching for Narrow-Band Thermal Emission with Graphene Ribbon Metasurface. Sensors. 2021; 21(20):6738. https://doi.org/10.3390/s21206738

Chicago/Turabian Style

Yada, Kyohei, Takashi Shimojo, Hideyuki Okada, and Atsushi Sakurai. 2021. "Theoretical and Numerical Analysis of Active Switching for Narrow-Band Thermal Emission with Graphene Ribbon Metasurface" Sensors 21, no. 20: 6738. https://doi.org/10.3390/s21206738

APA Style

Yada, K., Shimojo, T., Okada, H., & Sakurai, A. (2021). Theoretical and Numerical Analysis of Active Switching for Narrow-Band Thermal Emission with Graphene Ribbon Metasurface. Sensors, 21(20), 6738. https://doi.org/10.3390/s21206738

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop