We experimentally demonstrate the ability to create additional transmission resonances in a double-layer aperture array by varying the interlayer gap spacing. In the case of periodic aperture arrays, these additional resonances are sharply peaked, while for random aperture arrays the resonances are broad. Surprisingly, these additional resonances only occur when the interlayer gap spacing is greater than half the aperture spacing on a single array. Since there is no corresponding periodicity in the random arrays, these resonances occur regardless of how small the gap spacing is made. This phenomenon can be accurately modeled only if the correct frequency-dependent complex dielectric function of a metal film perforated with subwavelength apertures is used. Using THz time-domain spectroscopy, we are able to directly obtain the complex dielectric response function from the THz experimental transmission measurements. We conclude by demonstrating several passive free-space THz filters using multilayer aperture arrays. Importantly, we show that the magnitude of the lowest order resonance can be approximately maintained, while the background transmission can be significantly suppressed leading to a significant improvement in the optical filter fidelity.
© 2011 OSA
Technology development within the terahertz (THz) spectral range (1 THz = 1012 Hz) has lagged severely behind that of other regions of the electromagnetic spectrum [1,2]. While reasonable progress has been made in developing sources and detectors, other device technologies are still at a rather rudimentary stage [1,2]. A significant reason for this lies in the fact that dielectric materials typically used for passive and active device applications in microwave and optical applications tend to exhibit high loss at THz frequencies. Metals, on the other hand, are highly conductive at these frequencies , allowing for low propagation losses. The field of plasmonics , and the use of surface plasmon polaritons (SPPs), is therefore particularly attractive for THz device development. As an example, single layer periodic  and aperiodic  arrays of subwavelength apertures fabricated in metal films have been shown to exhibit resonantly enhanced transmission approaching unity [7,8], even for relatively small fractional aperture areas. Here we show that when two or more aperture arrays are placed in close proximity to one another, new transmission resonances arise that are associated with the gap spacing between the arrays, offering an additional mechanism for tailoring their optical properties. Importantly, proper design of these additional resonances requires that the complex dielectric response of the single aperture array be determined. Furthermore, with appropriate layer spacing in multilayer structures, we obtained extremely high transmission at the primary resonance frequency, along with substantial suppression of the background transmission increasing the filtering fidelity.
The resonant optical properties of subwavelength aperture arrays (also known as ‘plasmonic lattices’) have been explained in terms of SPPs that mediate the metal-light interaction via coupling to corrugated metal surfaces , where the electromagnetic field associated with the SPPs decays exponentially from the metal-dielectric interface. While the nearly exclusive focus of both theoretical and experimental work has been on the transmission response of single layer plasmonic lattices , there have been a number of experimental and theoretical studies of double-layer aperture arrays [10–15]. It was found that when the gap spacing, d between the two plasmonic lattices was large, the resulting transmission properties could be obtained by considering each array independently. However, when d becomes sufficiently small, it is expected that the coupling between the SPP waves on the inner, adjacent surfaces of the two layers would yield new and interesting spectral transmission characteristics. Here we show that this coupling forms additional resonances for small values of d that may be used in designing THz filters.
2. Experimental details
Two-dimensional (2D) aperture arrays were fabricated in 75 μm thick, ~5x5 cm2 area of free-standing stainless steel metal foils by laser cutting using a frequency tripled Nd:YAG laser. Two different aperture array pairs, based on periodic and random hole patterns respectively, were designed and fabricated to have the same fractional aperture area of ~19%. One aperture array pair consisted of 2D periodically spaced circular holes on a square lattice, while the other pair consisted of 2D randomly distributed circular holes. In Fig. 1 , we show a schematic diagram of the double-layer aperture array with the relevant dimensions.
We used THz time-domain spectroscopy (THz-TDS) to measure the optical transmission spectra, t(ω) of the perforated metal films, where the THz frequency ν = ω/2π . A pair of off-axis paraboloidal mirrors was used to collect and collimate the THz radiation from the emitter, which was then normally incident on the aperture array structures. Each array was attached to a solid metal plate with a 5 cm x 5 cm opening that was placed in the path of the collimated THz beam and placed on a translation stage to vary the distance between the arrays. The 1/e THz beam diameter was smaller than the aperture opening in the metal holders and the spatial extent of the aperture array structures, thereby minimizing edge effects due to the finite size of the array sample. Reference transmission spectra were taken with just the metal holders in the THz beam path using the same setup. The metallic foils were completely opaque prior to fabrication of the arrays. Therefore transmission through the perforated metallic films was uniquely due to the apertures. The detected transient photocurrent was then Fourier transformed and normalized to the reference transmission, yielding the electric field transmission spectrum, t(ω). THz TDS allows for the direct measurement of the THz electric field, yielding both amplitude and phase information. By transforming the time-domain data to the frequency domain, we are able to determine independently both the magnitude and phase of the amplitude transmission coefficient, t(ω), using the relation
In this expression, Eincident and Etransmitted are the incident and transmitted THz fields, respectively, |t(ω)| and φ(ω) are the magnitude and phase of the amplitude transmission coefficient, respectively, and ω /2π is the THz frequency.
3. Experimental results, simulation and discussion
Figure 2 summarizes our experimental studies on the periodically perforated metallic films. We begin by characterizing the properties of a single periodic aperture array in which the center-to-center aperture spacing, a = 1.5 mm. This plasmonic lattice forms a lowest order ( ± 1,0) anti-resonance (AR) frequency, AR1 = 0.2 THz and second lowest ( ± 1, ± 1) AR frequency, AR2 = 0.28 THz, as shown in the measured transmission spectrum in Fig. 2(a). We have shown previously that the AR feature (i.e. the sharp dip immediately to the right of the resonance peak) is directly related with the underlying structure factor . We note that although the fractional aperture area is only ~19%, the absolute resonant transmission amplitude at the lowest resonance, R1 is >80%. In Fig. 2(b), we show the obtained THz transmission spectra for a double-layer plasmonic lattice structure versus the interlayer spacing, d. The basic features of the transmission spectra can be summarized as follows. For d < d1 = 0.75 mm (where d1 = a/2), t(ω) for the double layer aperture array is similar to that of the corresponding single aperture array. Specifically, the resonance and anti-resonance frequencies are identical for the single and double arrays (Figs. 2(a) and 2(b)), although the double layer structure exhibits lower overall transmissivity. For values of d > d1, however, we obtain additional transmission resonances in the double layer array that are not present in the case of a single layer. We label these new resonances, meta-Fabry-Perot (MFP) resonances, as explained below. The new resonance, labeled MFP1 appears at frequencies ω1/2π < ω(R1)/2π, and moves progressively to lower frequencies with increasing d. However when d > d2 = 1.5 mm (where d2 = a), we obtain a second new resonance, labeled MFP2 (Fig. 2(b)) whose frequency, ω2/2π also decreases with d. The various new resonance frequencies are summarized in Fig. 2(c).
In order to determine the origin of the additional MFP resonances we performed numerical simulations based on the experimental parameters (Fig. 1) and the interlayer spacing, d. Since the MFP resonances emerge at array distances associated with integer values of a/2, namely d1 = (a/2) and d2 = 2(a/2), it is reasonable to consider Fabry-Perot reflections as the underlying mechanism. These reflections, however, need be modified according to the dielectric response of the underlying plasmonic lattice, which is actually a ‘metamaterial’. In the study of plasmonic lattices, a common approach for determining the transmission resonance frequencies is to use the dielectric properties of the (unperforated) metal film, although this typically yields only approximate resonance frequencies. In Fig. 3 we calculate the resulting transmission spectra taking all Fabry-Perot reflections between the two aperture array films into account, assuming that the dielectric properties of the aperture arrays can be approximated by that of an unstructured metal film at THz frequencies having a dielectric response ε ~-3x104 + i106 . The poor agreement with the experimental results demonstrates that a more sophisticated model for ε(ω) response of the plasmonic lattice needs be considered.
Using the THz-TDS complex transmission spectra associated with a single aperture array, we can simultaneously extract the real and imaginary components of the effective ε(ω) using the amplitude transmittance, |t(ω)|, and the phase spectra, φ(ω) , assuming that the effective permeability is frequency independent and of the order of 1 [17,18]. As shown in Fig. 4(a) , with fit parameters given in Table 1 , the complex ε(ω) response of the plasmonic lattices can be modeled using the general form:Eq. (2) corresponds to a smooth ‘plasma-type’ response, and is the only term used in modeling the random arrays (discussed below), while the last (summation) term in Eq. (2) corresponds to the discrete resonances associated with the periodic aperture arrays. The calculated ε(ω) along with the modeled ε(ω) using Eq. (2) are shown in Fig. 4(a) of a periodic aperture array having a = 1.5 mm and D = 750 µm. In fitting the data for the periodic arrays, we used the same effective plasma frequency () that we will use later in describing ε(ω) of the corresponding random aperture array. The TO frequencies () in the resonant contribution correspond to AR1 and AR2 frequencies (Fig. 2(a)), which can also be derived from the standard surface plasmon equation [4,5].
We calculated the transmission spectra of the double layer array using the obtained dielectric response, taking multiple reflections into account and assuming a single cycle broadband THz pulse at normal incidence. While contributions from Fabry-Perot reflections typically need to take all four surfaces of the double layer structure into account, for d > d1 we found that reflections between the inner array surfaces dominate. In Fig. 4(b), we show the resulting simulated spectra as a function of d; they are in good agreement with the experimental data. In Fig. 4(c), we plot the calculated frequencies associated with the R1, MFP1 and MFP2 resonances as a function of d, and compare these to the experimental frequencies; once again, the excellent agreement between the simulation and experimental data validates our approach.
In order to further demonstrate that the emergence of the MFP transmission resonances can be explained by using the correct complex dielectric function of the perforated metal film rather than the unperforated metallic response, we measured the transmission properties of a commensurate pair of random aperture arrays. First, we show the transmission spectrum of a single random array (Fig. 5(a) ). As expected, the transmission spectrum does not show evidence of any resonance, which is consistent with the fact that the Fourier transform of this geometry does not contain any discrete Fourier components. When two such arrays are separated by d, however we observe broad resonances in the otherwise smoothly varying transmission spectrum, as shown in Fig. 5(b). Using the magnitude and phase of the transmission spectrum associated with a single random array, we again extract the real and imaginary components of the effective ε(ω), as shown in Fig. 5(c). ε(ω) is essentially identical to that of the periodic array, but includes only the lossy plasma contribution (the first term of Eq. (2)) . The fit parameters are given in Table 2 . Using this extracted ε(ω) response we repeated the simulation mentioned above to obtain t(ω) as a function of d, as shown in Fig. 5(d), and summarized in Fig. 5(e). We again see excellent agreement between our experimental data and numerical computations. In contrast to what was observed with the periodic arrays, there does not appear to be a minimum d associated with the appearance of MFP resonances. This is consistent with the lack of periodicity in this double layer structure.
Based on the existence of MFP resonances between closely spaced aperture arrays, we are now able to extend this general idea for demonstrating highly efficient free-space THz filters, as summarized in Fig. 6 . In Fig. 6(a), we show the transmission properties of two 0.75 mm thick periodic arrays having d = 0 and d = 0.75 mm, respectively. For d = 0, the double-layer aperture array is equivalent to a single, 1.5 mm thick periodic aperture array. In fact the transmission spectrum is identical to what is observed for a single layer array (Fig. 2(a)), albeit with somewhat reduced amplitude. For d = 0.75 mm, no MFPs are expected. However, since d = 0.8 mm corresponds to a local maxima in t(ω) vs. d for the ( ± 1,0) resonance, we expect larger transmission compared with that for d = 0. In fact, the R1 transmission amplitude is nearly identical to that of the single layer array, although t(ω) at other frequencies above and below the R1 resonance is dramatically reduced. Similar behavior is expected for other array spacings that correspond to local maxima (i.e. the reflections yield constructive interference). The notion of maintaining high transmission response for the R1 resonance while reducing the background transmission can be implemented also using multiple arrays. In this case, the spacing between the arrays is significantly larger (~1 cm), as shown in Fig. 6(b). We introduce a small 5° rotation in the center array to further minimize the transmission away from the R1 resonance . While the R1 transmission for the triple array is slightly smaller than that of the single array, the resonance quality factor is higher, and the transmission away from the R1 resonance is dramatically minimized compared with a single array.
Our results demonstrate that carefully designed multilayer aperture arrays can allow for two important device capabilities: (1) the introduction of additional resonances associated with the layer spacing, which may be used to generate more complex THz spectral filtering properties, and (2) the ability to create high quality narrow bandpass THz filters. These capabilities were illustrated primarily using double-layer aperture arrays, although multilayer structures may allow for additional refinement of the transmission properties. In addition, the use of aperiodic geometries  as well as conductive non-metallic materials  may allow for greater control over the MFP resonance frequencies, as well as the potential for active THz optoelectronic devices applications. Finally, it should be noted that interesting “three dimensional” structures using stacks of aperture arrays layers in a variety of geometries may give rise to new and interesting optical interference phenomena.
This work was supported in part by the National Science Foundation under Grant # ECCS- 0801965 and the CEMRI program at the University of Utah. We gratefully acknowledge funding from the University of Utah Open Access Publishing Fund.
References and Links
1. M. Tonouchi, “Cutting-edge terahertz technology,” Nat. Photonics 1(2), 97–105 (2007). [CrossRef]
3. M. A. Ordal, L. L. Long, R. J. Bell, S. E. Bell, R. R. Bell, R. W. Alexander Jr, and C. A. Ward, “Optical properties of the metals Al, Co, Cu, Au, Fe, Pb, Ni, Pd, Pt, Ag, Ti, and W in the infrared and far infrared,” Appl. Opt. 22(7), 1099–20 (1983). [CrossRef] [PubMed]
5. T. W. Ebbesen, H. J. Lezec, H. F. Ghaemi, T. Thio, and P. Wolff, “Extraordinary optical transmission through sub-wavelength hole arrays,” Nature 391(6668), 667–669 (1998). [CrossRef]
8. J. W. Lee, M. A. Seo, J. Y. Sohn, Y. H. Ahn, D. S. Kim, S. C. Jeoung, Ch. Lienau, and Q.-H. Park, “Invisible plasmonic meta-materials through impedance matching to vacuum,” Opt. Express 13(26), 10681–10687 (2005). [CrossRef] [PubMed]
9. F. J. Garcia-Vidal, T. W. Ebbesen, and L. Kuipers, “Light passing through subwavelength apertures,” Rev. Mod. Phys. 82(1), 729–787 (2010). [CrossRef]
11. H. B. Chan, Z. Marcet, K. Woo, D. B. Tanner, D. W. Carr, J. E. Bower, R. A. Cirelli, E. Ferry, F. Klemens, J. Miner, C. S. Pai, and J. A. Taylor, “Optical transmission through double-layer metallic subwavelength slit arrays,” Opt. Lett. 31(4), 516–518 (2006). [CrossRef] [PubMed]
12. F. Miyamaru and M. Hangyo, “Anomalous terahertz transmission through double-layer metal hole arrays by coupling of surface plasmon polaritons,” Phys. Rev. B 71(16), 165408 (2005). [CrossRef]
14. H. Li, S. Xie, R. Zhou, Q. Liu, X. Zhou, and M. Yuan, “Two different transmission tunnels of light through double-layer gold nanohole arrays,” J. Phys. Condens. Matter 20(41), 415223 (2008). [CrossRef]
15. R. Ortuño, C. García-Meca, F. J. Rodríguez-Fortuño, J. Martí, and A. Martínez, “Role of surface plasmon polaritons on optical transmission through double layer metallic hole arrays,” Phys. Rev. B 79(7), 075425 (2009). [CrossRef]
18. F. J. Garcia-Vidal, L. Martin-Moreno, and J. B. Pendry, “Surfaces with holes in them: new plasmonic metamaterials,” J. Opt. A, Pure Appl. Opt. 7(2), S97–S101 (2005). [CrossRef]
19. T. D. Nguyen, A. Nahata, and Z. V. Vardeny, “THz anomalous transmission in plasmonic lattices: incidence angle dependence,” Proc. SPIE 7394, 73940H, 73940H-7 (2009). [CrossRef]
20. T. Matsui, Z. V. Vardeny, A. Agrawal, A. Nahata, and R. Menon, “Resonantly-enhanced transmission through a periodic array of subwavelength apertures in heavily-doped conducting polymer films,” Appl. Phys. Lett. 88(7), 071101 (2006). [CrossRef]