• Photonics Research
  • Vol. 6, Issue 5, 443 (2018)
Zi-Lan Deng1, Yaoyu Cao1, Xiangping Li1、3、*, and Guo Ping Wang2、4、*
Author Affiliations
  • 1Guangdong Provincial Key Laboratory of Optical Fiber Sensing and Communications, Institute of Photonics Technology, Jinan University, Guangzhou 510632, China
  • 2College of Electronic Science and Technology, Shenzhen University, Shenzhen 518060, China
  • 3e-mail: xiangpingli@jnu.edu.cn
  • 4e-mail: gpwang@szu.edu.cn
  • show less
    DOI: 10.1364/PRJ.6.000443 Cite this Article Set citation alerts
    Zi-Lan Deng, Yaoyu Cao, Xiangping Li, Guo Ping Wang, "Multifunctional metasurface: from extraordinary optical transmission to extraordinary optical diffraction in a single structure," Photonics Res. 6, 443 (2018) Copy Citation Text show less

    Abstract

    We show that a metasurface composed of a subwavelength metallic slit array embedded in an asymmetric dielectric environment can exhibit either extraordinary optical transmission (EOT) or extraordinary optical diffraction (EOD). The cascaded refractive indices of the dielectrics can leverage multiple decaying passages into variant subsections with different diffraction order combinations according to the diffraction order chart in the k-vector space, providing a flexible mean to tailor resonance decaying pathways of the metallic slit cavity mode by changing the wavevector of the incident light. As a result, either the zeroth transmission or 1st reflection efficiencies can be enhanced to near unity by the excitation of the localized slit cavity mode, leading to either EOT or EOD in a single structure, depending on the illumination angle. Based on this appealing feature, a multifunctional metasurface that can switch its functionality between transmission filter, mirror, and off-axis lens is demonstrated. Our findings provide a convenient way to construct multifunctional miniaturized optical components on a single planar device.

    1. INTRODUCTION

    Integrating various functionalities on a single miniaturized photonic device can largely simplify the overall optical system design, and therefore is a continuing trend in the nanophotonic area. There are many ways to realize multifunctional devices, such as the so-called “Janus” transformation optics device [1], acting as either a lens or a beam-shifter; parity-time-symmetric multilayers, which can simultaneously support extraordinary transmission and reflection [2]; and a single parity-time-symmetric cavity that can support both lasing and anti-lasing [3]. Those demonstrations require large coherent interaction lengths during light’s propagation in such devices.

    Recently, exploring various wavefront-shaping functionalities in a single ultrathin layer, i.e., metasurface, has been widely studied in the nanophotonic field due to the easy fabrication of a flat monolayer instead of a three-dimensional bulky structure [46]. Resorting to the abrupt phase change of each resonant element instead of accumulated interaction with large coherent length between light and structure, the ultrathin metasurface can provide various applications, such as beam deflection [711], focusing [1216], complex beam shaping [1723], holographic display [2429], and multiplexed hologram recording [3039]. In particular, single metasurfaces with multiple functionalities represent a newly emerging research frontier and receive onrushing developments [40,41]. Generally, the previous metasurface with high efficiencies can only work on either transmission mode [7] or reflection mode [9,27] depending on the limited number of deterministic decaying pathways. For transmission metasurfaces [42,43], only the transmission decaying channel is used for wavefront manipulation, while the reflection decaying channel is abandoned. For reflection metasurfaces [44,45], a metallic mirror is usually placed below the metasurface nanostructure to block all transmission channels, and only the reflection channel is tailored on demand. A metasurface with versatile functionalities demands multiple decaying pathways of its resonance, which inevitably leads to the tradeoff between functionalities and efficiencies owing to the competition among different decaying passages. In this paper, we present a metallic slit metasurface that harnesses both the transmission and reflection decaying channels of different diffraction orders, and demonstrate angularly switchable functionalities with nearly unitary efficiencies.

    It is well known that, for a metallic slit array with subwavelength period (p<λ), the localized resonance mode in the slit area leads to the enhanced transmission of incident light, which is known as extraordinary optical transmission (EOT) [46,47]. To obtain unitary diffraction efficiency of a grating, the usual method is to mimic a blazed grating by making a gradient phase profile with a discrete set of different nanoresonators in a large supercell [7,48,49]. Very recently, it has been shown that metasurfaces with identical unit cells, rather than the gradient phase approach, are also capable to obtain near-unitary diffraction efficiencies by harnessing the decaying pathways of the resonance cavity mode into high diffraction orders [20,50,51], which is referred to as extraordinary optical diffraction (EOD) [20]. Now, we consider a metallic subwavelength slit array embedded in an asymmetric dielectric environment exhibiting cascaded refractive indices without a metallic ground plate, and demonstrate that both EOT and EOD can happen in a single structure. According to the diffraction order chart in k-vector space, the resonance decaying pathways of the slit cavity mode can be tuned to either zeroth transmission or 1st reflection channel by changing the parallel wavevector of incident light. As a result, the metasurface exhibits resonantly enhanced transmission for small incident angles, acting as a transmission filter; exhibits total internal reflection (TIR) for moderate incident angles, acting as a refection mirror; and exhibits EOD for large incident angles, which can be applied for arbitrary wavefront shaping, such as deflecting, focusing, beam shaping, and hologram displaying. Such metasurfaces with versatile functionalities harnessing of multiple decaying channels are fundamentally different from the previous work based on single diffraction order enhancement approaches [20,40,50,52]. The demonstrated capability to control both the transmission and reflection channels of impinging light of our proposed metasurface provides a convenient and efficient way to realize multifunctional miniaturized optical components on a single surface.

    2. MODE EXPANSION THEORY MODEL

    Figure 1(a) illustrates the schematic of the multifunctional metasurface, which is a metallic subwavelength slit array (with slit width w, slit height h, and period p) embedded in an asymmetric environment. The refractive index of the medium in the incident side is larger than that of the medium in the transmission side (n1>n3) to make sure that there are intersection areas that are not only higher than zeroth reflection Woods anomaly (RWA) but also lower than zeroth transmission Woods anomaly (TWA) [green and red patches in Fig. 1(b)]. The refractive index in the slit area n2 is chosen properly so that the resonance wavelength of the slit cavity mode can cross all three kinds of areas with truncated diffraction order combinations, which are marked as blue, green, and red patches in Fig. 1(b). For small, moderate, and large incident angles, the beam will undergo EOT, TIR, and EOD, respectively, as the decaying pathways of the slit cavity mode are composed by different diffraction order combinations for different angles. Assume a plane wave with wavenumber k0 and parallel wavevector kx0 illuminating the structure from the superstrate; the possible propagating diffraction orders for both reflection (r0,r±1,) and transmission (t0,t±1,) can be found in the diffraction order chart as shown in Fig. 1(b). The diffraction order chart is constructed by RWA (solid) and TWA (dashed) lines, which are governed by the equations kx0+m2π/p=±n1k0 and kx0+m2π/p=±n3k0 (m is an integer), respectively. Above the mth (m=0,±1,±2,) RWA/TWA, the propagating reflection/transmission of the mth diffraction order (rm/tm) begins to appear, while below the zeroth RWA, which is also known as the light line, there are only evanescent waves without any radiation. For simplicity, we only consider the lowest areas of the diffraction order chart, where only a few diffraction orders (0 and ±1) exist. In the lowest sidelobe areas (green patches), there is only a r0 channel, where TIR happens; in the lowest central area (blue patch), there are only two propagating channels of r0 and t0, which may give rise to EOT; and in the higher sidelobe areas (red patches), there are only two propagating channels of r0 and r1/r1, allowing for the possibility of EOD. As the diffraction order chart is symmetric with respect to the k0 axis, we only consider the kx0>0 case in the following discussions.

    (a) Schematic of the metasurface structure that can exhibit EOT, TIR, or EOD, depending on the incident angle. The index of the upper layer, slit area, and lower layer is n1=2, n2=1.4, and n3=1, respectively. (b) Diffraction order chart in k-vector space of the structure in (a). Solid lines represent the RWAs with orders ±0, ±1, and ±2, while the dashed lines represent TWAs with orders ±0, ±1.

    Figure 1.(a) Schematic of the metasurface structure that can exhibit EOT, TIR, or EOD, depending on the incident angle. The index of the upper layer, slit area, and lower layer is n1=2, n2=1.4, and n3=1, respectively. (b) Diffraction order chart in k-vector space of the structure in (a). Solid lines represent the RWAs with orders ±0, ±1, and ±2, while the dashed lines represent TWAs with orders ±0, ±1.

    To obtain the transmission and reflection spectra of different diffraction orders, we first employ the analytical model expansion theory by considering the metal as a perfect electric conductor (PEC) [20,53]. When a plane wave with parallel wavevector kx0=n1k0sinθ and transverse magnetic (TM) polarization illuminates on the metallic slit array, by matching the boundary conditions, which require that the tangential components of both the electric field (Ex) and the magnetic field (Hy) are continuous at interfaces, we can analytically obtain the zeroth reflectance R0, zeroth transmittance T0, and 1st reflectance R1 with closed forms as follows (see Appendix A): where, ϵ1, ϵ2, and ϵ3 are the permittivities of the upper layer, slit area, and lower layer, respectively; a0+ and a0 are forward and backward complex amplitude coefficients of the cavity mode in the subwavelength slit area, which can be obtained by solving linear equations with where kx0=k0sinθ, kz0(1)=n12k02kx02 is the wavevector of the incident plane wave, and kxm=kx0+2mπ/p, kzm(j)=nj2k02kxm2 (j=1,3) is the wavevector of the mth diffraction order in the superstrate and substrate, respectively.

    3. RESULTS AND DISCUSSION

    Figures 2(a)2(c) show the analytically obtained R0, T0, and R1 for varying wavelengths and incident angles, by Eqs. 1(a)1(c), respectively. The parameters of the metallic slit array are chosen as p=1  μm, w=0.1  μm, h=1  μm; the refractive indices of the upper layer, the slit area, and the lower layer are n1=2, n2=1.4, and n3=1, which can be realized by silicon nitride, silica, and air, respectively. In Fig. 2(a), there are three abrupt color changing contours, which correspond to the zeroth TWA (critical angle for TIR), 1st RWA, and 1st TWA in Fig. 1(b), respectively. The subareas ① (upper-left), ② (upper-middle), and ③ (upper-right) in Fig. 2(a) separated by those WA lines correspond to the blue patch, green patch, and red patch in Fig. 1(b), respectively. Thus, there are only r0 and t0 channels in ①, only r0 in ②, and only r0, r1 channels in ③, respectively. We see that R0 exhibits dips in subareas ① and ③, at the resonant wavelength of about 3.3 μm, which corresponds to the fundamental cavity mode in the metallic silts. For different incident angles, the resonant wavelength stays the same, because the slit cavity mode that leads to the resonant dip of R0 is a localized resonant mode, whose excitation is independent of the incident angles. Correspondingly, T0 exhibits a resonant peak approaching unity at the same position of the R0 dip in subarea ① [θ0<asin(n3/n1)=30°] as shown in Fig. 2(b). Note that the subarea ① is the widely studied zeroth-order-diffraction zone in the metamaterial community; many phenomena such as EOT [46], the generalized Brewster effect [5456], and perfect absorption [57] are found in this zone, while the influence of the resonance behavior on areas other than subarea ① receives scarce attention. When θ0>30°, T0 is always zero because the incident angle is larger than the critical angle for TIR. Figure 2(c) shows that R1 also exhibits a peak at the same position of the R0 dips in subarea ③. Therefore, the slit cavity mode indeed can enhance the diffraction efficiency of the 1st order to near unity, as long as the propagating channels are restricted to only r0 and r1. We note that, between the EOT and EOD, there is a transitional area (subarea ②), where R0 is constantly unity [Fig. 2(a)]. This is due to the fact that r0 is the only propagating channel in this subarea; the resonant behavior of the slit cavity can only lead to a phase change of the reflection coefficient, while the reflection amplitude always stays unity. However, when the absorption loss channel is considered in the practical case, there may be a resonant absorption phenomenon that can also lead to the R0 dip in subarea ②.

    Analytical results of (a) zeroth reflectance R0, (b) zeroth transmittance T0, and (c) −1st reflectance R−1, respectively, of the metasurface for varying incident angle θ and wavelength λ. The three abrupt color changing contours in (a) automatically show the critical angles for TIR, the −1st RWA, and the −1st TWA, respectively, as marked in the figure, and they split the whole phase map into different subareas that contain combinations of different propagating channels (also denoted in the figure). In subarea ①, there are r0 and t0 channels; R0 exhibits a dip in (a), while T0 exhibits a peak in (b) under the resonance condition of the slit cavity mode. In subarea ②, there is only a r0 channel; R0 is unity in the whole subarea. In subarea ③, there are r0 and r−1 channels; R0 exhibits a dip in (a), while R−1 exhibits a peak in (c) under the resonance condition of the slit cavity mode.

    Figure 2.Analytical results of (a) zeroth reflectance R0, (b) zeroth transmittance T0, and (c) 1st reflectance R1, respectively, of the metasurface for varying incident angle θ and wavelength λ. The three abrupt color changing contours in (a) automatically show the critical angles for TIR, the 1st RWA, and the 1st TWA, respectively, as marked in the figure, and they split the whole phase map into different subareas that contain combinations of different propagating channels (also denoted in the figure). In subarea ①, there are r0 and t0 channels; R0 exhibits a dip in (a), while T0 exhibits a peak in (b) under the resonance condition of the slit cavity mode. In subarea ②, there is only a r0 channel; R0 is unity in the whole subarea. In subarea ③, there are r0 and r1 channels; R0 exhibits a dip in (a), while R1 exhibits a peak in (c) under the resonance condition of the slit cavity mode.

    According to the theory we discussed above, now let us consider the practical case of real metallic structures. Silver is chosen as the metal material due to its small ohmic losses. The permittivity of silver is obtained by fitting the experimental data [58] to the Drude model. Figure 3 shows the R0, T0, and R1 spectra by the finite element method (FEM), which is implemented by the commercial package COMSOL. At incident angles of 0° [Fig. 3(a)] and 15° [Fig. 3(b)] that belong to the subarea ①, R0 (green) exhibits resonant dips below 0.1, while T0 (blue) exhibits peaks exceeding 0.9 at the wavelength of about 3.3 μm, which is the typical EOT phenomenon. At the incident angle of 35° [Fig. 3(c)] that belongs to the subarea ②, both R1 and T0 are suppressed as zero, while R0 is the only nonzero value after the Wood’s anomaly wavelength (λ=3.0  μm). At incident angles of 45° [Fig. 3(d)] and 60° [Fig. 3(e)] that belong to the subarea ③, R0 (green) also exhibits resonant dips at the wavelength around 3.2 μm, while R1 (red) exhibits resonant peaks exceeding 0.9 at the corresponding wavelength, which is the EOD phenomenon. We note that there are sharp peaks for R1 and sharp dips for R0 with Fano lineshape at the short wavelength side in Figs. 3(b)3(d), which is due to the Wood’s anomalies near those wavelengths. Compared with the model expansion theory results [dashed lines in Fig. 3], the peak values of T0 and R1 of the FEM simulation results are slightly lower, and the resonant wavelength shifts a little to shorter wavelengths, because the real metal inevitably has absorption loss and finite negative permittivity, which give rise to overall power loss and the decrease of the effective wavelength of the slit mode. Nevertheless, the FEM simulation results for real metal (solid lines in Fig. 3) are generally consistent with the theoretical results (dashed lines in Fig. 3). The appearance of the EOT and EOD relies on the different resonant decaying pathways (zeroth or 1st diffraction order channels) of the local slit cavity mode; therefore we can tailor the working wavelength and bandwidth of EOT and EOD through the slit parameters in the wavelength range (n1+n3)p<λ<2n1p, where all three kinds of subareas ①, ②, and ③ can be reached.

    Reflection (R0,R−1) and transmission (T0) spectra of the metasurface composed of real metal (silver), for different incident angles (a) θ=0°, (b) θ=15°, (c) θ=35°, (d) θ=45°, and (e) θ=60°, respectively. Dashed curves show the corresponding analytical results of the PEC metasurface for comparison.

    Figure 3.Reflection (R0,R1) and transmission (T0) spectra of the metasurface composed of real metal (silver), for different incident angles (a) θ=0°, (b) θ=15°, (c) θ=35°, (d) θ=45°, and (e) θ=60°, respectively. Dashed curves show the corresponding analytical results of the PEC metasurface for comparison.

    In general, the resonant wavelength is linearly dependent on the slit height h, and the bandwidth is positively correlated to the slit width w. Thus, the proposed structure can tailor the working condition of the multiple functionalities on demand.

    To clearly illustrate the evolution from EOT to EOD, we fix the working wavelength at 3.3 μm, and plot the R0, T0, and R1 as a function of the incident angle as shown in Fig. 4(a). In the practical experimental setup, the apparatus to steer the beam angle can be facilely realized by many methods, such as scanning mirror, rotational stage under the metasurface, or circular rail under the illumination source [38,59]. It indicates that, as we gradually increase the incident angle from 0° to 90°, the beam first undergoes EOT (T0>0.8; R0,R1<0.1). And after the critical angle of 30°, the beam experiences TIR, where R0 is the only nonzero component. To further increase the incident angle to cross the 1st RWA line (θ>asin[λ/(n1p)1]=40°), the EOD begins to happen, where R1 becomes the largest component, while R0 and T0 are effectively suppressed. Figures 4(b)4(f) show the field patterns (|Hz|) for a Gaussian beam illuminating the metasurface with different incident angles, which are calculated by the COMSOL package. We can clearly see that, for small incident angles (0°, 15°), the beam passes through the metallic slit metasurface [Figs. 4(b) and 4(c)] even if the slit width is in deep subwavelength range (w=0.1  μmλ=3.3  μm), and the relation between the transmission angle θt and incidence angle θ obeys the Snell law (n1sinθ=n3sinθt). For incident angle 35° in subarea ②, the incident beam undergoes ordinary TIR, with the reflection angle equal to the incident angle [Fig. 4(d)]. For large incident angles (45°, 55.6°) in subarea ③, the incident beam undergoes negative reflection, with the reflected beam on the same side of the incident beam with respect to the normal [Figs. 4(e) and 4(f)], which is the typical behavior of EOD. Retroreflection occurs when the incident angle is 55.6° (Littrow mount angle); that is, the beam is reflected back along its original path, even for oblique incidence [60]. Note that the beam width of the incident beam is set as w0=7.5λ, which gives the divergence angle of 5°; the above simulation results reveal that such a small divergence angle has little influence on the angular tailoring of multiple functionalities. In a realistic experiment, the footprint of the proposed metasurface can be designed on the order of 200 μm by 200 μm, which can be compatible with a weakly focused Gaussian beam with a small divergence angle limited to less than 5° (for example, beam width 10 mm together with a lens focal length f=300  mm gives a divergence angle of about 2°).

    (a) Reflection (R0, R−1) and transmission (T0) as a function of the incident angle θ at the working wavelength λ=3.3 μm for the periodic slit array of which the zoom-in structure is shown in the upper panel. Dashed curves show the corresponding analytical results of PEC metasurface for comparison. (b)–(f) Field patterns (|Hz|) of a Gaussian beam with λ=3.3 μm illuminating the metasurface (the area marked by the red dashed frame) with different incident angles: (b) 0°, (c) 15°, (d) 35°, (e) 45°, and (f) 55.6°, respectively.

    Figure 4.(a) Reflection (R0, R1) and transmission (T0) as a function of the incident angle θ at the working wavelength λ=3.3  μm for the periodic slit array of which the zoom-in structure is shown in the upper panel. Dashed curves show the corresponding analytical results of PEC metasurface for comparison. (b)–(f) Field patterns (|Hz|) of a Gaussian beam with λ=3.3  μm illuminating the metasurface (the area marked by the red dashed frame) with different incident angles: (b) 0°, (c) 15°, (d) 35°, (e) 45°, and (f) 55.6°, respectively.

    Based on EOD, we can design arbitrary wavefront-shaping metasurfaces with the binary computer-generated hologram approach [61]. As an example, we design a multifunctional metasurface that can behave as an off-axis lens in the EOD domain by modulating the metallic slit array with a quadratic phase profile φ(x,y)=k0x2/(2f). The modulated structure is based on the periodic grating structure with geometric parameters w=0.1  μm, p=1  μm, and h=1  μm; after the modulation, both the slit width and the period of the grating will exhibit linear gradient profiles as shown in Fig. 5(a). For small incident angles [Figs. 5(b) and 5(c)], the beam still can pass through the nonperiodic subwavelength slit array, exhibiting the EOT phenomenon. For the incident angle of 35° [Fig. 5(d)], the metasurface acts as an ordinary mirror, making the incident beam undergo TIR. For large incident angles [Figs. 5(e) and 5(f)], the metasurface becomes an off-axis lens, which can focus the incident beam into a single line in the 1st diffraction direction, while the specular reflection in the right side is near-totally suppressed. Therefore, for different incident angles, the metasurface indeed can act as an optical device with totally different functionalities. This angle-dependent multifunctional device may be potentially useful in panoramic viewing systems. Although we only designed 1D functional devices, arbitrary 2D wavefront-shaping devices can also be obtained by modulating the structure in both x and y directions according to a predefined phase profile, which typically results in a curved grating profile [20].

    (a) Schematic of the modulated metasurface with a gradient grating profile with a zoomed-in structure on the upper panel. (b)–(f) Field patterns (|Hz|) of a Gaussian beam with λ=3.3 μm illuminating the modulated metasurface (the area marked by the red dashed frame) with different incident angles: (b) 0°, (c) 15°, (d) 35°, (e) 45°, and (f) 55.6°, respectively. In (b) and (c), the incident beam also undergoes EOT. In (d), the incident beam still undergoes TIR. In (e) and (f), the incident beam undergoes the focusing capability, with near-total suppression of the zeroth diffraction order.

    Figure 5.(a) Schematic of the modulated metasurface with a gradient grating profile with a zoomed-in structure on the upper panel. (b)–(f) Field patterns (|Hz|) of a Gaussian beam with λ=3.3  μm illuminating the modulated metasurface (the area marked by the red dashed frame) with different incident angles: (b) 0°, (c) 15°, (d) 35°, (e) 45°, and (f) 55.6°, respectively. In (b) and (c), the incident beam also undergoes EOT. In (d), the incident beam still undergoes TIR. In (e) and (f), the incident beam undergoes the focusing capability, with near-total suppression of the zeroth diffraction order.

    4. CONCLUSION

    In conclusion, we show that a metallic subwavelength slit array embedded in an asymmetric environment can support EOT, TIR, and EOD simultaneously. The cascaded dielectric environment can leverage multiple decaying passages into variant subsections with truncated diffraction order combinations according to the diffraction order chart in the k-vector space, which allows the wavevector-dependent tunneling of the slit cavity mode to different diffraction orders. For small incident angles, the resonance decaying pathway goes to the zeroth transmission channel, giving rise to enhanced transmission to near unity, while for large incident angles, the resonance decaying pathway switches to the 1st reflection channel, leading to enhanced diffraction to near unity. Between the EOT and EOD areas, there is also a transitional area with moderate incident angles, where ordinary specular reflection happens. Based on the above theory, a multifunctional metasurface that can behave as a transmission filter, a mirror, and an off-axis lens simultaneously is demonstrated. Our proposed metasurface can be used to realize multifunctional optical components on a single planar device.

    APPENDIX A

    To deduce Eqs. (1)–(3) in the main text, let us consider that a plane wave with an arbitrary incident angle (or parallel wavevector kx0), where TM polarization (magnetic field is always parallel to the slit direction) illuminates the structure from the superstrate. Owing to the diffraction effect of the grating, the electromagnetic (EM) field in the superstrate and the substrate can be expanded as plane waves of all diffraction orders (including both propagating and evanescent components). As a result, the transverse EM field components in the superstrate and the substrate have the following forms: H1y=H0(eikx0x+ikz0(1)z m= rneikxmx ikzm(1)z),E1x=kz0(1)H0ω 1 0(eikx0x+ikz0(1)z+m= rnkzm(1)kz0(1)eikxmx ikzm(1)z),H3y=H0m= tmeikxmx+ikzm(3)z,E3x=kz0(3)H0ω 3 0m= tnkzm(3)kz0(3)eikxmx+ikzm(3)z,where kx0, kz0(1)=n12k02 kx02 is the wavevector of the incident plane wave, and kxm=kx0+2mπ/p, kzm(j)=nj2k02 kxm2 (j=1,3) is the wavevector of the mth diffraction order in the superstrate and the substrate, respectively. tm and rm are the complex transmission and reflection coefficients for the mth diffraction order, respectively. 0 is the permittivity of vacuum; ω is the angular frequency, and H0 is the amplitude of the incident wave.

    In the slit area, there are generally multiple waveguide modes in each metallic slit. However, for the ultra-narrow slit in the deep subwavelength scale, the higher order modes are cut off, and hence only the fundamental mode needs to be considered. We assume a0+ and a0 are the forward and backward amplitude coefficients in each slit, respectively; then the transverse EM field components in the metallic slit area are Hy=H0 rect(x/w)(a0+ein2k0z a0 e in2k0z),for 0<z<h,Ex=n2k0H0ω 2 0 rect(x/w)(a0+ein2k0z+a0 e in2k0z),for 0zh,where rect(·) is the rectangular function, and the amplitude coefficients a0+ and a0 are unknown variables that can be determined by the boundary condition requiring that the parallel EM field components should be continuous at all interfaces.

    At the interface z=0, the boundary condition yields eikx0x m= rmeikxmx=rect(x/w)(a0+ a0 ),eikx0x+m= rmkzm(1)kz0(1)eikxmx= 1 2n2k0kz0(1) rect(x/w)(a0++a0 ),and at the interface z=h, the boundary condition yields m= tmeikxmx+ikzm(3)h=rect(x/w)(a0+ein2k0h a0 e in2k0h),m= tmkzm(3)kz0(3)eikxmx+ikzm(3)h= 3 2n2k0kz0(3) rect(x/w) (a0+ein2k0h+a0 e in2k0h).

    Considering the orthogonality of different diffraction orders, we multiply Eqs. (A4b) and (A5b) by e ikxnx and then integrate over p/2<x<p/2, respectively. This yields pδm0+kzm(1)kz0(1)prm= 1 2n2k0kz0(1)w sin c(kxmw2)(a0++a0 ),kzm(3)kz0(3)ptmeikzm(3)h= 3 2n2k0kz0(3)w sin c(kxmw2) (a0+ein2k0h+a0 e in2k0h).

    Thus, tm, rm can be expressed in terms of ak+, ak as rm= kz0(1)kzm(1)δm0+ 1 2n2k0kzm(1)wp sin c(kxmw2)(a0++a0 ),tm=e ikzm(3)h 3 2n2k0kzm(3)wp sin c(kxmw2)(a0+e in2k0h+a0 ein2k0h).

    Multiplying Eqs. (A4a) and (A5a) by rect(x/w) and then integrating in terms of x in region p/2<x< p/2, respectively, yields sin c(kx0w2)=m= rm sin c(kxmw2)+a0+ a0 ,0=m= tmeikzm(3)h sin c(kxmw2) (a0+ein2k0h a0 e in2k0h).

    Substituting Eqs. (A7a) and (A7b) into Eqs. (A8a) and (A8b), respectively, yields 2sin c(kx0w2)= 1 2wp(a0++a0 )m= n2k0kzm(1) sin c2(kxmw2)+(a0+ a0 ),0= 3 2wp(a0+ein2k0h+a0 e in2k0h)m= n2k0kzm(3) sin c2(kxmw2) (a0+ein2k0h a0 e in2k0h).

    From Eq. (A9), we can build a 2×2 matrix M with elements as expressed in Eq. (3) in the main text.

    Then, the amplitude coefficients of cavity modes in metallic slits a=(a0+,a0 )T can be solved by a linear matrix equation, a=M 1b, where b=(2 sin c(kx0w2),0)T. Finally, substituting the solved amplitude coefficients a0+ and a0 back into Eq. (A7), we can readily obtain the reflection and transmission coefficients for arbitrary diffraction orders.

    In the main text, we focus on the zeroth reflection r0, zeroth transmission t0, and 1st reflection r 1. They can be analytically expressed as r0= 1 2n2k0kz0(1)wp sin c(kx0w2)(a0++a0 ) 1,t0=e ikz0(3)h 3 2n2k0kz0(3)wp sin c(kx0w2)(a0+e in2k0h+a0 ein2k0h),r 1= 1 2n2k0kz, 1(1)wp sin c(kx, 1w2)(a0++a0 ).

    The corresponding power reflectance and transmittance are R0=|r0|2, T0=|t0|2, and R 1=|r 1|2, respectively, which are exactly the expressions in Eq. (1) in the main text.

    References

    [1] T. Zentgraf, J. Valentine, N. Tapia, J. Li, X. Zhang. An optical ‘Janus’ device for integrated photonics. Adv. Mater., 22, 2561-2564(2010).

    [2] S. Ding, G. P. Wang. Extraordinary reflection and transmission with direction dependent wavelength selectivity based on parity-time-symmetric multilayers. J. Appl. Phys., 117, 023104(2015).

    [3] Z. J. Wong, Y.-L. Xu, J. Kim, K. O’Brien, Y. Wang, L. Feng, X. Zhang. Lasing and anti-lasing in a single cavity. Nat. Photonics, 10, 796-801(2016).

    [4] A. V. Kildishev, A. Boltasseva, V. M. Shalaev. Planar photonics with metasurfaces. Science, 339, 1232009(2013).

    [5] N. Yu, F. Capasso. Flat optics with designer metasurfaces. Nat. Mater., 13, 139-150(2014).

    [6] N. Meinzer, W. L. Barnes, I. R. Hooper. Plasmonic meta-atoms and metasurfaces. Nat. Photonics, 8, 889-898(2014).

    [7] N. Yu, P. Genevet, M. A. Kats, F. Aieta, J.-P. Tetienne, F. Capasso, Z. Gaburro. Light propagation with phase discontinuities: generalized laws of reflection and refraction. Science, 334, 333-337(2011).

    [8] F. Aieta, P. Genevet, N. Yu, M. A. Kats, Z. Gaburro, F. Capasso. Out-of-plane reflection and refraction of light by anisotropic optical antenna metasurfaces with phase discontinuities. Nano Lett., 12, 1702-1706(2012).

    [9] S. Sun, Q. He, S. Xiao, Q. Xu, X. Li, L. Zhou. Gradient-index meta-surfaces as a bridge linking propagating waves and surface waves. Nat. Mater., 11, 426-431(2012).

    [10] L. Huang, X. Chen, H. Mühlenbernd, G. Li, B. Bai, Q. Tan, G. Jin, T. Zentgraf, S. Zhang. Dispersionless phase discontinuities for controlling light propagation. Nano Lett., 12, 5750-5755(2012).

    [11] D. Lin, P. Fan, E. Hasman, M. L. Brongersma. Dielectric gradient metasurface optical elements. Science, 345, 298-302(2014).

    [12] F. Aieta, P. Genevet, M. A. Kats, N. Yu, R. Blanchard, Z. Gaburro, F. Capasso. Aberration-free ultrathin flat lenses and axicons at telecom wavelengths based on plasmonic metasurfaces. Nano Lett., 12, 4932-4936(2012).

    [13] X. Chen, L. Huang, H. Mühlenbernd, G. Li, B. Bai, Q. Tan, G. Jin, C.-W. Qiu, S. Zhang, T. Zentgraf. Dual-polarity plasmonic metalens for visible light. Nat. Commun., 3, 1198(2012).

    [14] X. Li, S. Xiao, B. Cai, Q. He, T. J. Cui, L. Zhou. Flat metasurfaces to focus electromagnetic waves in reflection geometry. Opt. Lett., 37, 4940-4942(2012).

    [15] M. Khorasaninejad, W. T. Chen, R. C. Devlin, J. Oh, A. Y. Zhu, F. Capasso. Metalenses at visible wavelengths: diffraction-limited focusing and subwavelength resolution imaging. Science, 352, 1190-1194(2016).

    [16] F. Zhou, W. Cao, B. Dong, T. Reissman, W. Zhang, C. Sun. Additive manufacturing of a 3D terahertz gradient-refractive index lens. Adv. Opt. Mater., 4, 1034-1040(2016).

    [17] I. Dolev, I. Epstein, A. Arie. Surface-plasmon holographic beam shaping. Phys. Rev. Lett., 109, 203903(2012).

    [18] P. Genevet, J. Lin, M. A. Kats, F. Capasso. Holographic detection of the orbital angular momentum of light with plasmonic photodiodes. Nat. Commun., 3, 1278(2012).

    [19] J. Lin, P. Genevet, M. A. Kats, N. Antoniou, F. Capasso. Nanostructured holograms for broadband manipulation of vector beams. Nano Lett., 13, 4269-4274(2013).

    [20] Z.-L. Deng, S. Zhang, G. P. Wang. A facile grating approach towards broadband, wide-angle and high-efficiency holographic metasurfaces. Nanoscale, 8, 1588-1594(2016).

    [21] B. G. Cai, Y. B. Li, W. X. Jiang, Q. Cheng, T. J. Cui. Generation of spatial Bessel beams using holographic metasurface. Opt. Express, 23, 7593-7601(2015).

    [22] W. Zhu, M. Jiang, H. Guan, J. Yu, H. Lu, J. Zhang, Z. Chen. Tunable spin splitting of Laguerre-Gaussian beams in graphene metamaterials. Photon. Res., 5, 684-688(2017).

    [23] T. Cao, Y. Li, X. Zhang, Y. Zou. Theoretical study of tunable chirality from graphene integrated achiral metasurfaces. Photon. Res., 5, 441-449(2017).

    [24] L. Huang, X. Chen, H. Mühlenbernd, H. Zhang, S. Chen, B. Bai, Q. Tan, G. Jin, K.-W. Cheah, C.-W. Qiu, J. Li, T. Zentgraf, S. Zhang. Three-dimensional optical holography using a plasmonic metasurface. Nat. Commun., 4, 2808(2013).

    [25] X. Ni, A. V. Kildishev, V. M. Shalaev. Metasurface holograms for visible light. Nat. Commun., 4, 2807(2013).

    [26] Y.-H. Chen, L. Huang, L. Gan, Z.-Y. Li. Wavefront shaping of infrared light through a subwavelength hole. Light Sci. Appl., 1, e26(2012).

    [27] G. Zheng, H. Mühlenbernd, M. Kenney, G. Li, T. Zentgraf, S. Zhang. Metasurface holograms reaching 80% efficiency. Nat. Nanotechnol., 10, 308-312(2015).

    [28] K. E. Chong, L. Wang, I. Staude, A. R. James, J. Dominguez, S. Liu, G. S. Subramania, M. Decker, D. N. Neshev, I. Brener, Y. S. Kivshar. Efficient polarization-insensitive complex wavefront control using Huygens’ metasurfaces based on dielectric resonant meta-atoms. ACS Photon., 3, 514-519(2016).

    [29] Z.-L. Deng, G. Li. Metasurface optical holography. Mater. Today Phys., 3, 16-32(2017).

    [30] W. T. Chen, K.-Y. Yang, C.-M. Wang, Y.-W. Huang, G. Sun, I. D. Chiang, C. Y. Liao, W.-L. Hsu, H. T. Lin, S. Sun, L. Zhou, A. Q. Liu, D. P. Tsai. High-efficiency broadband meta-hologram with polarization-controlled dual images. Nano Lett., 14, 225-230(2013).

    [31] C. Min, J. Liu, T. Lei, G. Si, Z. Xie, J. Lin, L. Du, X. Yuan. Plasmonic nano-slits assisted polarization selective detour phase meta-hologram. Laser Photon. Rev., 10, 978-985(2016).

    [32] D. Wen, F. Yue, G. Li, G. Zheng, K. Chan, S. Chen, M. Chen, K. F. Li, P. W. H. Wong, K. W. Cheah, E. Yue, B. Pun, S. Zhang, X. Chen. Helicity multiplexed broadband metasurface holograms. Nat. Commun., 6, 8241(2015).

    [33] M. Khorasaninejad, A. Ambrosio, P. Kanhaiya, F. Capasso. Broadband and chiral binary dielectric meta-holograms. Sci. Adv., 2, e1501258(2016).

    [34] M. Khorasaninejad, W. T. Chen, J. Oh, F. Capasso. Super-dispersive off-axis meta-lenses for compact high resolution spectroscopy. Nano Lett., 16, 3732-3737(2016).

    [35] W. Wan, J. Gao, X. Yang. Full-color plasmonic metasurface holograms. ACS Nano, 10, 10671-10680(2016).

    [36] X. Li, L. Chen, Y. Li, X. Zhang, M. Pu, Z. Zhao, X. Ma, Y. Wang, M. Hong, X. Luo. Multicolor 3D meta-holography by broadband plasmonic modulation. Sci. Adv., 2, e1601102(2016).

    [37] M. Q. Mehmood, S. Mei, S. Hussain, K. Huang, S. Y. Siew, L. Zhang, T. Zhang, X. Ling, H. Liu, J. Teng, A. Danner, S. Zhang, C.-W. Qiu. Visible-frequency metasurface for structuring and spatially multiplexing optical vortices. Adv. Mater., 28, 2533-2539(2016).

    [38] L. Huang, H. Mühlenbernd, X. Li, X. Song, B. Bai, Y. Wang, T. Zentgraf. Broadband hybrid holographic multiplexing with geometric metasurfaces. Adv. Mater., 27, 6444-6449(2015).

    [39] H. Ren, X. Li, Q. Zhang, M. Gu. On-chip noninterference angular momentum multiplexing of broadband light. Science, 352, 805-809(2016).

    [40] M. Khorasaninejad, F. Capasso. Broadband multifunctional efficient meta-gratings based on dielectric waveguide phase shifters. Nano Lett., 15, 6709-6715(2015).

    [41] D. Wen, S. Chen, F. Yue, K. Chan, M. Chen, M. Ardron, K. F. Li, P. W. H. Wong, K. W. Cheah, E. Y. B. Pun, G. Li, S. Zhang, X. Chen. Metasurface device with helicity-dependent functionality. Adv. Opt. Mater., 4, 321-327(2016).

    [42] M. Pu, X. Li, X. Ma, Y. Wang, Z. Zhao, C. Wang, C. Hu, P. Gao, C. Huang, H. Ren, X. Li, F. Qin, J. Yang, M. Gu, M. Hong, X. Luo. Catenary optics for achromatic generation of perfect optical angular momentum. Sci. Adv., 1, e1500396(2015).

    [43] D. Tang, C. Wang, Z. Zhao, Y. Wang, M. Pu, X. Li, P. Gao, X. Luo. Ultrabroadband superoscillatory lens composed by plasmonic metasurfaces for subdiffraction light focusing. Laser Photon. Rev., 9, 713-719(2015).

    [44] X. Luo. Principles of electromagnetic waves in metasurfaces. Sci. China Phys. Mech. Astron., 58, 594201(2015).

    [45] Z. Zhao, M. Pu, H. Gao, J. Jin, X. Li, X. Ma, Y. Wang, P. Gao, X. Luo. Multispectral optical metasurfaces enabled by achromatic phase transition. Sci. Rep., 5, 15781(2015).

    [46] T. W. Ebbesen, H. J. Lezec, H. F. Ghaemi, T. Thio, P. A. Wolff. Extraordinary optical transmission through sub-wavelength hole arrays. Nature, 391, 667-669(1998).

    [47] J. A. Porto, F. J. García-Vidal, J. B. Pendry. Transmission resonances on metallic gratings with very narrow slits. Phys. Rev. Lett., 83, 2845-2848(1999).

    [48] F. Aieta, M. A. Kats, P. Genevet, F. Capasso. Multiwavelength achromatic metasurfaces by dispersive phase compensation. Science, 347, 1342-1345(2015).

    [49] Y. F. Yu, A. Y. Zhu, R. Paniagua-Domínguez, Y. H. Fu, B. Luk’yanchuk, A. I. Kuznetsov. High-transmission dielectric metasurface with 2π phase control at visible wavelengths. Laser Photon. Rev., 9, 412-418(2015).

    [50] Y. Ra’di, D. L. Sounas, A. Alù. Metagratings: beyond the limits of graded metasurfaces for wave front control. Phys. Rev. Lett., 119, 067404(2017).

    [51] Z.-L. Deng, S. Zhang, G. P. Wang. Wide-angled off-axis achromatic metasurfaces for visible light. Opt. Express, 24, 23118-23128(2016).

    [52] N. M. Estakhri, V. Neder, M. W. Knight, A. Polman, A. Alù. Visible light, wide-angle graded metasurface for back reflection. ACS Photon., 4, 228-235(2017).

    [53] Z.-L. Deng, N. Yogesh, X.-D. Chen, W.-J. Chen, J.-W. Dong, Z. Ouyang, G. P. Wang. Full controlling of Fano resonances in metal-slit superlattice. Sci. Rep., 5, 18461(2015).

    [54] A. Alù, G. D’Aguanno, N. Mattiucci, M. Bloemer. Plasmonic Brewster angle: broadband extraordinary transmission through optical gratings. Phys. Rev. Lett., 106, 123902(2011).

    [55] R.-H. Fan, R.-W. Peng, X.-R. Huang, J. Li, Y. Liu, Q. Hu, M. Wang, X. Zhang. Transparent metals for ultrabroadband electromagnetic waves. Adv. Mater., 24, 1980-1986(2012).

    [56] R. Paniagua-Domínguez, Y. F. Yu, A. E. Miroshnichenko, L. A. Krivitsky, Y. H. Fu, V. Valuckas, L. Gonzaga, Y. T. Toh, A. Y. S. Kay, B. Luk’yanchuk, A. I. Kuznetsov. Generalized Brewster effect in dielectric metasurfaces. Nat. Commun., 7, 10362(2016).

    [57] Y. Cui, Y. He, Y. Jin, F. Ding, L. Yang, Y. Ye, S. Zhong, Y. Lin, S. He. Plasmonic and metamaterial structures as electromagnetic absorbers. Laser Photon. Rev., 8, 495-520(2014).

    [58] P. B. Johnson, R. W. Christy. Optical constants of the noble metals. Phys. Rev. B, 6, 4370-4379(1972).

    [59] X. Zhang, J. Jin, M. Pu, X. Li, X. Ma, P. Gao, Z. Zhao, Y. Wang, C. Wang, X. Luo. Ultrahigh-capacity dynamic holographic displays via anisotropic nanoholes. Nanoscale, 9, 1409-1415(2017).

    [60] A. Arbabi, E. Arbabi, Y. Horie, S. M. Kamali, A. Faraon. Planar metasurface retroreflector. Nat. Photonics, 11, 415-420(2017).

    [61] W.-H. Lee. Binary computer-generated holograms. Appl. Opt., 18, 3661-3669(1979).

    Zi-Lan Deng, Yaoyu Cao, Xiangping Li, Guo Ping Wang, "Multifunctional metasurface: from extraordinary optical transmission to extraordinary optical diffraction in a single structure," Photonics Res. 6, 443 (2018)
    Download Citation