Design principles for infrared wide-angle perfect absorber based on plasmonic structure

Size: px

Start display at page:

Download "Design principles for infrared wide-angle perfect absorber based on plasmonic structure"

Oswin Chase
6 years ago
Views:

1 Design principles for infrared wide-angle perfect absorber based on plasmonic structure Mingbo Pu, Chenggang Hu, Min Wang, Cheng Huang, Zeyu Zhao, Changtao Wang, Qin Feng, and Xiangang Luo* State Key Laboratory of Optical Technologies for Microfabrication, Institute of Optics and Electronics, Chinese Academy of Science, P.O.box 35, Chengdu 6129, China Abstract: An approach for designing a wide-angle perfect absorber at infrared frequencies is proposed. The technique is based on a perfectly impedance-matched sheet (PIMS) formed by plasmonic nanostructure. It is shown that the effective impedance is more physical meaningful and beneficial than effective medium in describing the electromagnetic properties of metamaterial absorber. As a specific implementation of this technique, a wide-angle polarization-independent dual-band absorber is numerically demonstrated at frequencies of 1THz and 28THz with absorption close to 1% simultaneously. Circuit models are utilized to describe the impedance property of localized plasmon modes and the results show good agreement with that retrieved from reflection coefficient at normal incidence. 211 Optical Society of America OCIS codes: (5.6624) Subwavelength structures; ( ) Metamaterials; (24.668) Surface plasmons; (26.574) Resonance. References and Links 1. N. I. Landy, S. Sajuyigbe, J. J. Mock, D. R. Smith, and W. J. Padilla, Perfect metamaterial absorber, Phys. Rev. Lett. 1(2), 2742 (28). 2. H. Tao, N. I. Landy, C. M. Bingham, X. Zhang, R. D. Averitt, and W. J. Padilla, A metamaterial absorber for the terahertz regime: design, fabrication and characterization, Opt. Express 16(1), (28). 3. H. Tao, C. M. Bingham, A. C. Strikwerda, D. Pilon, D. Shrekenhamer, N. I. Landy, K. Fan, X. Zhang, W. J. Padilla, and R. D. Averitt, Highly flexible wide angle of incidence terahertz metamaterial absorber: design, fabrication, and characterization, Phys. Rev. B 78(24), 24113(R) (28). 4. T. V. Teperik, F. J. García de Abajo, A. G. Borisov, M. Abdelsalam, P. N. Bartlett, Y. Sugawara, and J. J. Baumberg, Omnidirectional absorption in nanostructured metal surfaces, Nat. Photonics 2(5), (28). 5. C. Hu, Z. Zhao, X. Chen, and X. Luo, Realizing near-perfect absorption at visible frequencies, Opt. Express 17(13), (29). 6. Q. Y. Wen, H. W. Zhang, Y. S. Xie, Q. H. Yang, and Y. L. Liu, Dual band terahertz metamaterial absorber: design, fabrication, and characterization, Appl. Phys. Lett. 95(24), (29). 7. M. Diem, T. Koschny, and C. M. Soukoulis, Wide-angle perfect absorber/thermal emitter in the terahertz regime, Phys. Rev. B 79(3), 3311 (29). 8. Y. Avitzour, Y. A. Urzhumov, and G. Shvets, Wide-angle infrared absorber based on a negative-index plasmonic metamaterial, Phys. Rev. B 79(4), (29). 9. N. I. Landy, C. M. Bingham, T. Tyler, N. Jokerst, D. R. Smith, and W. J. Padilla, Design, theory, and measurement of a polarization-insensitive absorber for terahertz imaging, Phys. Rev. B 79(12), (29). 1. Y. Q. Ye, Y. Jin, and S. He, Ominidirectional, polarization-insensitive and broadband thin absorber in the terahertz regime, J. Opt. Soc. Am. B 27(3), (21). 11. J. Hao, J. Wang, X. Liu, W. J. Padilla, L. Zhou, and M. Qiu, High performance optical absorber based on a plasmonic metamaterial, Appl. Phys. Lett. 96(25), (21). 12. X. Liu, T. Starr, A. F. Starr, and W. J. Padilla, Infrared spatial and frequency selective metamaterial with nearunity absorbance, Phys. Rev. Lett. 14(2), 2743 (21). 13. N. Liu, M. Mesch, T. Weiss, M. Hentschel, and H. Giessen, Infrared perfect absorber and its application as plasmonic sensor, Nano Lett. 1(7), (21). 14. J. Hao, L. Zhou, and M. Qiu, Nearly total absorption of light and heat generation by plasmonic metamaterials, Phys. Rev. B 83(16), (211). 15. Y. Li, Y. Xie, H. Zhang, Y. Liu, Q. Wen, and W. Ling, The strong non-reciprocity of metamaterial absorber: characteristic, interpretation and modeling, J. Phys. D Appl. Phys. 42(9), 9548 (29). 16. C. Caloz and T. Itoh, Electromagnetic Metamaterial: Transmission Line Theory and Microwave Applications (Wiley, 25). # $15. USD Received 9 May 211; revised 17 Jul 211; accepted 17 Aug 211; published 19 Aug 211 (C) 211 OSA 29 August 211 / Vol. 19, No. 18 / OPTICS EXPRESS 17413

2 17. L. Fu, H. Schweizer, H. Guo, N. Liu, and H. Giessen, Synthesis of transmission line models for metamaterial slabs at optical frequencies, Phys. Rev. B 78(11), (28). 18. N. Engheta, Circuits with light at nanoscales: optical nanocircuits inspired by metamaterials, Science 317(5845), (27). 19. Q. Y. Wen, Y. S. Xie, H. W. Zhang, Q. H. Yang, Y. X. Li, and Y. L. Liu, Transmission line model and fields analysis of metamaterial absorber in the terahertz band, Opt. Express 17(22), (29). 2. E. D. Palik, Handbook of Optical Constants of Solids (Academic, Orlando, Fla., 1985). 21. T. A. Kelf, Y. Sugawara, R. M. Cole, J. J. Baumberg, M. Abdelsalam, S. Cintra, S. Mahajan, A. Russell, and P. Bartlett, Localized and delocalized plasmons in metallic nanovoids, Phys. Rev. B 74(24), (26). 1. Introduction Metamaterial (MM) absorbers have attracted much attention due to their unique ability to achieve narrowband near-unity absorption in wide frequency region ranging from microwave to terahertz and optical frequencies [1 13]. The sub-wavelength size of unit cell allows many relevant applications including focal plane array (FPA) detectors for imaging [1], thermal emitters [7], and plasmonic sensor [13]. Generally, MM is characterized by its effective permittivity and permeability in the framework of effective medium theory (EMT) [1 3,9,12,14]. However, MM absorber doesn t exactly satisfy the homogeneous-effective limit as it is not uniform and periodic in the direction of wave propagation [15]. In addition, surface plasmons (SPs) can also bring on near-perfect absorption at optical frequency when the period of the surface structure is large enough compared with the wavelength [4,5]. In another theoretical regime, transmission line theory (TLT) is used to study the properties of metamaterials from microwave frequency to optical frequency [16 18]. In this regime, the relation between structure parameters and electromagnetic characters become more obvious. Recently, a transmission line model has been proposed to analyze the threelayer MM absorber in the terahertz band [19]. Although the results are in agreement with simulated ones, the magnetic coupling between two metallic layers is not taken into account and the design methodology is still not clear. In this paper, the proposed absorber consists of a metallic surface structure (metasurface) and a thick metal layer separated by a thin dielectric spacer. The ideal impedance for the metasurface is theoretically deduced to absorb all electromagnetic radiation. Infrared dual band absorption is numerically demonstrated and the conditions of perfect absorption are presented by utilizing equivalent circuit models. The angle-independent absorption is found to be related to the localized plasmon modes, while the SPs induced absorption is eliminated by reducing the period of the surface structure. 2. Perfectly impedance-matched sheet air a as 11 Y s J dielectric spacer exp(-ikd) 1 r m exp(ikd) r m ground plane E, y H, z k, x d Fig. 1. Schematic of the absorber. It composes of a metasurface, dielectric spacer and metallic ground plane. In the dielectric spacer and surrounding space, total fields consist of both the forward and backward going waves. The reflection of thick metal layer and the absorber is denoted as r m and S 11. The amplitude of E field of forward going wave at the ground plane is chosen to be 1 in order to utilize the transfer matrix formulation while the amplitude of incident E field is denoted as a. As shown in Fig. 1, a metasurface is placed in front of a thick metal layer (denoted as ground plane) with a distance d. The metasurface can be thought as an impedance sheet characterized # $15. USD Received 9 May 211; revised 17 Jul 211; accepted 17 Aug 211; published 19 Aug 211 (C) 211 OSA 29 August 211 / Vol. 19, No. 18 / OPTICS EXPRESS 17414

3 by its impedance. By matching the impedance of ground plane to that of vacuum, perfect absorption can be realized and the corresponding sheet is called as perfectly impedancematched sheet (PIMS). The impedance of PIMS can be obtained by setting reflection coefficient S 11 to be zero. Assuming the dielectric spacer is non-magnetic (µ = 1), the boundary conditions for the sheet at normal incidence can be obtained (L denotes the left side of the impedance sheet and R denotes the right side): L R E = E : a+ as = exp( ikd) + r exp( ikd) 11 L R H = H : Y ( a as ) = Y (exp( ikd ) r exp( ikd)) + J 11 1 ( ) L( R) J = Ys E = Ys a+ as11 m m (1) where Y, Y1 = ε1y and Ym = ε my are the intrinsic admittance of vacuum, dielectric spacer and metal respectively, ε 1 and ε m are permittivities of dielectric and metal. k and k k r = ( Y Y ) / ( Y + Y ) is = ε1 are the wave vector in the vacuum and dielectric spacer. m 1 m 1 m the reflectivity of thick metal layer. J is the current flowing in the sheet while a is the amplitude of incident E field. By eliminating a in Eq. (1), the admittance of PIMS can be obtained: exp( ikd) rm exp( ikd) Ys = 1/ Zs = Y Y1. (2) exp( ikd) + rm exp( ikd) The ideal impedance (Z s = R-iX, where R and X are resistance and reactance) for PIMS at normal incidence is independent of polarization. As depicted in Fig. 2, the impedance should be pure resistive (R = 377Ω, X = ) corresponding to the traditional Salisbury screen when d is near λ/4n ( n= ε1 is the refractive index of the dielectric spacer). In Zone1, the required resistance should be very small and the corresponding reactance should be capacitive (R>, X<). In contrast, the required reactance is inductive (R>, X>) for Zone2. 3. PIMS design Impedance (Ohm) Zone1 R X Salisbury Zone Zone2 Z /2 -Z / nd/lambda Fig. 2. Ideal impedances of PIMS versus the effective thickness of dielectric layer (2nd/λ). n is chosen to be 2 and r m is 1 corresponding to the reflection coefficient of PEC. As the thickness of the absorber is much smaller than the working wavelength in Zone1 as shown in Fig. 2, Zone1 is typically the case for most practical applications. To mimic the impedance of PIMS, lossy metal should be deposited in appropriate geometrical patterns such as square patches, crosses as proposed in the recent designs [11,12,14]. Z # $15. USD Received 9 May 211; revised 17 Jul 211; accepted 17 Aug 211; published 19 Aug 211 (C) 211 OSA 29 August 211 / Vol. 19, No. 18 / OPTICS EXPRESS 17415

4 Here, we choose circular patches arrayed in hexagonal lattice due to its high azimuthal symmetry. As illustrated in Fig. 3, the radius of cylindrical patch is r and the center distance between adjacent patches is a. The thicknesses of patch and dielectric spacer are t and d, respectively. The thickness of thick metal layer h is far larger (1nm in our case) than the penetration depth to suppress transmission. All metallic parts are chosen as gold and experimental value of permittivity is used [2]. Al 2 O 3 is used as dielectric spacer and the dielectric constant is 2.92 at wavelengths around 3µm (Al 2 O 3 is lossless in this frequency regime [2]). (a) Gold Gold a Al 2 O 3 y r 6 y d h x z z x t Fig. 3. (a) Front view and side view of the absorber. The rectangular region in (a) is the unit cell used in our simulations. We simulated the reflection coefficient S 11 using frequency domain Finite Element method (FEM). The unit cell is shown in Fig. 3 (a) with dashed line and periodic boundary condition is used to mimic infinite cell array. The absorption is calculated as A = 1-S 11 2 since there is no transmittance. For TE (TM) polarization, the E (H) field direction is along y-axis. As an example, the geometrical parameters are adjusted to optimize the absorption at frequency of 1THz (λ = 3µm). Following the same strategy as Eq. (1) and Eq. (2), the effective admittance of the surface structure can be retrieved from the simulated reflection coefficient S 11 : 1 S exp( ikd ) rm exp( ikd) Ysretri = 1/ Zsretri = Y Y. 1+ S exp( ikd) + r exp( ikd) In fact, the ideal admittance given in Eq. (2) is obtained on a more specific condition (S 11 = ). When r = 277nm, a = 693nm, t = 32nm, d = 3nm, the retrieved impedance for TM polarization overlap with the ideal one at f = 1THz and near-perfect absorption is realized as shown in Fig. 4 (a). Interestingly, another absorption peak occurs at the right side of the resonant frequency f = 245THz. The simulated absorption coefficients are 99.96% and 98.89% for 1THz and 28THz while the bandwidths are very narrow for both peaks. It is worth noting that the two absorption peaks both occur at the intersection points of the retrieved and ideal impedances. In fact, the retrieved impedance is similar with that of a typical parallel RLC circuit with resonance frequency at f = 245THz. Although the resonant frequency of the retrieved impedance does not correspond to any absorption peak, the resonant characteristic is the key of impedance matching for the two peaks. To obtain an accurate circuit model, the field distributions in the structure should be explored at different frequencies. As shown in Figs. 5 (a, c), two different resonances in the voids between the patches and ground plane are found at 1THz and 28THz. Using the classical equivalent circuit theory [16,17] to analyze the electromagnetic properties of the two resonant modes, the effective impedance in the whole simulated frequency region can be written as m (3) # $15. USD Received 9 May 211; revised 17 Jul 211; accepted 17 Aug 211; published 19 Aug 211 (C) 211 OSA 29 August 211 / Vol. 19, No. 18 / OPTICS EXPRESS 17416

5 1 1 Zseff ( ω) = +, f < f 1/ ( R+ jωl) + jωc jωc Zseff ( ω) = + +, f > fc, 1/ ( R + jωl ) + jωc 1/ ( R + jωl ) + jωc jωc c (4) (a) Absorption (c) Absorption Sim. Absorption Retr. R Retr. X Ideal R Ideal X Retr. R Retr. X Ideal R Ideal X Cal. Absorption Cal. R Cal. X Ideal R Ideal X Frequency (THz) Fig. 4. Absorption of the dual-band absorber for TM polarization wave and the corresponding impedances versus the frequency. The simulated absorption and retrieved impedance are shown in (a), while that calculated from equivalent circuits are shown in (c). Scaled impedances retrieved from S parameters are shown in. Near perfect absorption peaks are achieved at the intersection points of retrieved impedances and the ideal ones. The 9% absorption bandwidths are 5.5% and 1.59% for 1THz and 28THz, respectively. The results for TE polarization are the same due to high azimuthal symmetry of the structure. Impedance (Ohm) Impedance (Ohm) (a) E H k R L C C 1THz R L C C E normal H tangential (c) E H k (d) 28THz R 1 R L1 2 L 2 C 3 C 1 C 2 E normal H tangential E H J -Max Max Fig. 5. (a)(c) Side view of the field distribution and the corresponding circuit models at 1THz and 28THz. (d) Front view of the maximal normal electric field and tangential magnetic field at the center of dielectric spacer. There is 9 phase shift between E and H fields. # $15. USD Received 9 May 211; revised 17 Jul 211; accepted 17 Aug 211; published 19 Aug 211 (C) 211 OSA 29 August 211 / Vol. 19, No. 18 / OPTICS EXPRESS 17417

6 where f c 2THz is the frequency where mode begins to change. The corresponding equivalent circuit models are shown in the insets of Figs. 5 (a) and (c). At frequency below f c, only the lowest resonance mode is responsible for the circuit model. However, the applicability of the single resonance model breaks down and higher order modes must be taken into account at higher frequencies. Although the circuit elements in Eq. (4) are frequency dependent in general [18], they can be assumed to be constant in a narrow frequency band. In general, the resistors (R, R 1, R 2 ) come from the intrinsic loss in metal (the dielectric spacer is lossless here) and become larger for smaller patch thickness t. The capacitors (C, C, C 1, C 2, C 3 ) result from the electric field distribution in dielectric between two metallic elements and tend to be larger for smaller distance between them. The inductors (L, L 1, L 2 ) are related with the current distribution in the metallic patches. These relations between the circuit elements and geometrical parameters are used in the optimization of absorptance. To obtain the circuit parameters, one must fit the impedance calculated using Eq. (4) to the retrieved one. Because there are so many parameters to be determined, the relations between these parameters and impedances must be taken advantage to solve the problem. As an example, the resonance frequency and bandwidth are determined by parallel connected R (R 1, R 2 ), L (L 1, L 2 ) and C (C 1, C 2 ) while the overall capacitance can be tuned by the series connected capacitor (C, C 3 ). The fitted parameters for the dual-band band absorber are: R = 1.7Ω, L = 2fH, C = 2aF, C = 27aF, R 1 = 4.1Ω, L 1 = 75fH, C 1 = 5.7aF, R 2 = 4Ω, L 2 = 8fH, C 2 = 2.3aF, C 3 = 12aF. The corresponding impedances and absorption are depicted in Fig. 4 (c), in good agreement with that in Fig. 4 (a). Thus, the impact of cavity modes on effective impedance is successfully modeled through parallel connection of RLC circuits. 4. Localized and delocalized plasmons 3 (a) Frequency(THz) Frequency(THz) SPs Angle of incidence (deg) Angle of incidence (deg) Fig. 6. Absorption of the dual-band absorber as a function of frequency and the angle of incidence for (a) TE and TM polarizations. The angle-invariant absorption results from the localized plasmon, while the SPs induced absorption is shown at the top right corner of. The angle dependences of the absorption of the dual-band absorber are shown in Fig. 6. The absorption for TM (TE) polarization at 1THz is larger than 9% at incidence angle of 7 (55 ), while the absorption is still larger than 8% at incidence angle of 3 (45 ) at 28THz. This kind of wide-angle absorption has been found at other nano-structured metal surfaces and is attributed to the localized void plasmons [4,21]. The field distribution in Fig. 5 demonstrates that most energy is localized in the cavity formed by the patches and ground plane. As shown at the top right corner of Fig. 6, there is an angle dependent absorption curve, which is related to SPs. The SPs cannot be excited when the period of the structure is in sub-wavelength scale. However, this effect will become much more obvious when the period is comparable with the working wavelength, as discussed later. At large incidence angles, other absorption peaks also occur at frequencies around 18THz for TM polarization due to excitation of horizontal resonance modes. # $15. USD Received 9 May 211; revised 17 Jul 211; accepted 17 Aug 211; published 19 Aug 211 (C) 211 OSA 29 August 211 / Vol. 19, No. 18 / OPTICS EXPRESS 17418

7 To further exploit the influence of localized modes and delocalized SPs on the absorption properties, another absorber at 1THz with increased geometrical parameters (Sample2) is designed. Compared with the absorption of the first design (Fig. 4), two additional absorption peaks around 27THz are observed (Fig. 7). The circuit parameters are not calculated here because the retrieved impedances are sufficient for the optimization. As can be seen from the retrieved impedances, there are three resonances at frequencies around 155THz, 255THz and 26THz, corresponding to the drastic variations of the curves. Absorption Absorption Retr. R Retr. X Ideal R Ideal X Impedance (Ohm) Frequency (THz) -4 Fig. 7. Absorption of Sample2 at normal incidence and the corresponding effective impedances. The geometric parameters are r = 34nm, a = 1385nm, t = 4nm, d = 6nm, h = 1nm. The absorption coefficient at 1THz is 99.84% and 9% absorption bandwidth is 3.54%, respectively. The absorption around 27THz, however, is rather complicated due to higher order resonances and SPs coupling. 3 (a) Frequency(THz) Frequency(THz) Angle of incidence (deg) Angle of incidence (deg) Fig. 8. Absorption of Sample2 as a function of frequency and the angle of incidence for (a) TE and TM polarizations. The red dashed lines illustrate the lowest plasmon mode excited by the reciprocal vector along z-axis. The angle dependences of the absorption for Sample2 are shown in Fig. 8. The absorption for TM (TE) polarization at 8 (55 ) is still 9% at 1THz, which is less angle dependent compared with the first absorber. Absorptions at frequencies above 15THz are characterized by very bad angular stability. Obviously, the absorption peak at 1THz is mainly related with the localized plasmon resonance, while absorption peaks at higher frequencies are related with the interplay of localized and delocalized plasmon resonances. As wide-angle perfect absorber is needed in most applications, the influence of SPs should be eliminated in the design. In general, the coupling of propagating wave to delocalized SPs in the nanostructured periodic surfaces is realized by a reciprocal vector added to the wave vector of the incident wave and this mechanism can be explained using Bragg scattering theory [21]: # $15. USD Received 9 May 211; revised 17 Jul 211; accepted 17 Aug 211; published 19 Aug 211 (C) 211 OSA 29 August 211 / Vol. 19, No. 18 / OPTICS EXPRESS 17419

8 sin β = neff k = k θ+ q mn (5) 1 1 where neff 1/ ε m + ε d is an approximation of the effective refractive index of SPs mode and q = 4 π / a is the lowest reciprocal vector along z-axis. The theoretical result is 11(11) shown by red dashed lines in Fig. 8. According to this theory, reciprocal vector will decreases and the coupling frequency will shifts to lower frequency as the patch separation a increases. For the first absorber, the patch separation a is 693nm and there is only a SPs tail at the top right corner in Fig. 6. For sample2 (a = 1385nm), however, the SPs coupling is obvious and the angle-stability of the high frequency absorption is dramatically decreased. In the design process, the patch separation a should be chosen to be less than 2λ/(n eff + 1) in order to remove the impact of SPs. This sub-wavelength scale is just one of the conditions for the homogenization of MMs. 5. Conclusion In conclusion, we have proposed a methodology based on plasmonic circuit theory to design metamaterials absorber at infrared frequency. Based on this methodology, a dual-band polarization-insensitive absorber at 1THz and 28THz with absorption of 99.96% and 98.89% is numerically demonstrated. Circuit models are constructed to describe the frequency-dependent impedance property. The influence of localized and delocalized plasmonic modes on the angle dependence of absorption is discussed. It is shown that wideangle absorption can be realized by utilizing the properties of localized plasmonic modes and eliminating the impact of surface plasmons. Although the design procedure is devoted to perfect absorber at infrared frequencies, our conclusions can be extended to other frequency bands. Acknowledgments This work was supported by 973 Program of China (No. 211CB318). # $15. USD Received 9 May 211; revised 17 Jul 211; accepted 17 Aug 211; published 19 Aug 211 (C) 211 OSA 29 August 211 / Vol. 19, No. 18 / OPTICS EXPRESS 1742

Engineering heavily doped silicon for broadband absorber in the terahertz regime

Engineering heavily doped silicon for broadband absorber in the terahertz regime Mingbo Pu, Min Wang, Chenggang Hu, Cheng Huang, Zeyu Zhao, Yanqin Wang, and Xiangang Luo * State Key Laboratory of Optical