Abstract
A plasmonic metamaterial is proposed for which an array of subwavelength apertures is pierced into a metallic foil whose one flat surface has been made with periodic rectangle holes of finite depth. Designed surface plasmons sustained by the holes are explored when the size and spacing of the holes are much smaller than those of the apertures. The transmission property of electromagnetic waves through the metamaterials is analyzed. Results show that the designed surface plasmons characterized by the holes could support the transmission resonances of the incident wave passing through the subwavelength apertures, and that the peak transmission wavelengths could be designed by controlling the geometrical and optical parameters of the holes. Example is taken at THz regime. Our work proposes a method to design the peak wavelengths, and may affect further engineering of surface plasmon optics, especially in THz to microwave regimes.
©2009 Optical Society of America
1. Introduction
Coupling between light and the collective oscillations of the electron plasma at a metal–dielectric interface gives rise to surface plasmons (SPs) [1]. These are the SP states, responsible for a host of phenomena unique to metals, for which the phenomenon of enhanced transmission of electromagnetic waves through subwavelength apertures in metal films or foils has stimulated significant interest in recent years [2] since the pioneering work by Ebbesen and his coauthors [3]. It was believed that the interaction of the incident radiation with the SPs leads to an enhancement of the transmission [4]. The resonant excitation of SPs creates huge electric fields at the surface that force light through the apertures, giving very high transmission coefficients. Such a viewpoint meets challenge at low frequency regime, in which metals can be regarded as perfect conductors, and excitation of the SPs cannot be expected, and consequently the SP excitation does not contribute to the enhanced transmission property. To overcome this problem, Pendry and his coauthors presented a new idea that there should be surface plasmon polariton-like bound surface states, namely mimicking surface plasmons (MSPs), on structured surfaces of a perfect conductor [5]. The existence of MSPs was indirectly verified using microwave reflectivity measurements [6]. MSPs promise the ability to engineer a SP at almost any frequency, thus MSP has a name designer surface plasmon (DSP) [6]. For a semi-infinite perfect conductor perforated with a two-dimensional (2D) square array of square holes whose size and spacing are much smaller than the incident wavelength, by taking a method in which a perfect conductor was optically described in the long wavelength limit as an effective medium, Pendry and his coauthors demonstrated this symmetric structure to be of a dielectric function of plasmon form characterized by its structural parameters [5]. By using the effective-medium method and the well-known method of modal expansion, they repeated their results and analyzed the other two structures made by 1D arrays of grooves and 2D arrays of square holes of infinite or finite depth [7]. In fact, the effective-medium method had been used to consider 1D periodic slits for designing metallic metamaterials with a high index of refraction [8]. Besides, the field-averaging method could provide a useful means of characterizing the effective-medium parameters of an arbitrary metamaterial [9]. It was believed that a measured enhanced transmission phenomenon at THz regime should originate from the crucial contribution of DSPs [10]. DSP propagation on structured metal surfaces was directly observed and the results were quantitatively consistent with theoretical predictions [11]. Based on a modal expansion of electromagnetic fields, a rigorous method for analyzing DSPs on a periodically corrugated metal surface has been formulated [12]. The existence of DSPs has been confirmed by rigorous numerical modeling (such as using finite element [13] and finite time domain [14] numerical methods) and by means of both analytical techniques and rigorous numerical solution of Maxwell’s equations [15]. The explicit analytical expressions for determining the dispersion of DSPs sustained by a two-dimensional array of square holes were derived by employing a modified modal matching approach [16]. The DSP theory was used to study the slow-light characteristics of structured metal surfaces with a one-dimensional periodic rectangle [17] or with a graded grating of varying depth in THz domain [18].
DSPs promise the ability to engineer a SP on structured metal surfaces. It is the character distinguishes DSPs from SPs. This character allows us to actively design some relevant experiments, such as microwave reflectivity measurements [6] and highly confined guiding THz surface wave [11]. However, based on the DSP theory to actively design the transmission property of electromagnetic waves through subwavelength apertures have not been considered so far. This mainly originates from the fact that the analytical form of the dispersion relationship for DSPs has not been found when the DSPs are sustained by subwavelength apertures whose size and spacing are comparing to the incident wavelength. As we know that the analytical form is available for DSPs sustained by arrays of holes whose size and spacing are much smaller than the incident wavelength, we naturally address an idea as the following: Can we construct a kind of new plasmonic metamaterials by which we can use the DSP theory to actively design an experiment of the transmission resonances of electromagnetic waves through subwavelength apertures in THz domain? Such metamaterials are made of metallic foils. At first, some rectangle holes of finite depth are made on one flat surface of the foil with the size and spacing of the holes being much smaller than the incident wavelength. And then arrays of subwavelength apertures are pierced on the foil. We imagine that the DSP states sustained by the holes of finite depth could support the transmission resonances of incident radiation through the subwavelength apertures. In this paper, based on the DSP’s theory [5,7], we obtain an analytical form of the dispersion relationship for the anisotropic DSP states sustained by the holes, and then analyze theoretically the transmission property of radiation through the subwavelength apertures. Results indicate that the DSP states characterized by the holes could support the transmission resonances of THz radiation through the subwavelength apertures. The formula of the peak transmission wavelengths is obtained, which involves the geometrical and optical parameters of the metamaterials, including those of the holes, indicating that the peak wavelengths could be designed by controlling the parameters of the holes. The shape of holes is chosen as rectangle because it involves three geometrical parameters, i.e., long side, short side and depth, more than that of square or circle shape, thus allowing us to adjust more parameters to design the peak wavelengths. Numerical examples are taken at THz regime. Our research propose a method to use the DSP theory to design experiments of transmission resonances, and may affect further engineering of surface plasmon optics, especially in THz to microwave regimes.
2. Theoretical model
At first, we pay our attention on the structure of a metallic foil pierced with periodic rectangle holes of infinite depth and are devoted to find the effective medium of this structure. A set of a × b rectangle holes arranged on a d × d lattice is made into a square metallic foil, as shown in Fig. 1 (a) . A cartesian coordinate is built for which the x- and y-axis are parallel to the sides a and b, respectively, and the z-axis is normal to the surface of the foil. For such a structure, if , where λ is the wavelength of radiation, externally incident radiation is insensitive to the details of the holes, which it can see only an average response that can be described by an effective anisotropic homogeneous medium with effective permittivity , , and , and effective permeability ,, and , for which and were determined as and by observing that the dispersion of the waveguide mode in the holes is unaffected by parallel momentum [5,7]. In what follows, we will endeavor to find the expressions of ,, , and for the asymmetric structure.
For the waveguide modes excited by an incident wave, both electric and magnetic fields are zero inside the conductor but in the rectangle holes the electric fields of transverse electric (TE) waveguide modes TEmn have the form
where ,,, m, n = 0,1,2,3⋯, is vacuum permeability, is constant, x, y and z are cartesian coordinates, ω is the angular frequency, t is time, andwhere is the free space wave vector, , and are the permittivity and permeability, respectively, of any material that may be filling the rectangle holes. The magnetic fields follow from Eq. (1), where the superscript T denotes the operation of transposition.For these waveguide modes, only those modes TEm0 and TE0n (m, n ≥1) can sustain DSP states because the average fields at the surface are equal to zero for the modes TE mn with m≠0 and n ≠0. As we are interested in the limit, we assume that the fundamental modes will dominate because they are the least strongly decaying. We thereby only consider two modes: TE10 and TE01. For the mode TE10, we have and with
where . We then suppose that there exist effective homogeneous fields in the effective medium as and within which is the same as in the waveguide, i.e., , and is defined by the incident direction. If this effective field is to match to the incident and reflected fields external to the surface, and must give the same average fields at the surface. Hence, by the use of the two equations above, after matching to incident and reflected waves, we haveBy taking the same approach, we have = 0 because and ≠0, indicating that the effective homogeneous fields only have two components and , thus we have and . We then consider the instantaneous flows of energy across the surface, and . They must be the same inside and outside the surface, both for the real and effective media, thus we haveFrom the two equations above, we haveSubstitute and Eq. (7) into Eq. (2) with m = 1 and n = 0, we have where c is the velocity of wave in a vacuo.For the mode TE01, we have and with
where . We then suppose that the effective homogeneous fields in the medium are and within which , too, and is defined by the incident direction. For this case, similarly, we have . By the use of the same approach as above, we obtain The results above show that ≠. This is easy understood for the asymmetric structure. Taking a = b, we can get the results for the case of square holes [5,7].We then consider the DSP states sustained by a structured surface with rectangle holes of finite depth. As shown Fig. 1 (b), a set of a × b rectangle holes with depth h arranged on a d × d lattice is made into one surface of a square metallic foil. In the effective medium approximation the structure displayed in Fig. 1 (b) behaves as an homogeneous but anisotropic layer of thickness h on top of a perfect conductor, as shown in Fig. 1 (c). For an interface between the effective medium layer and an isotropic homogeneous dielectric medium with , there should be an anisotropic bound surface state, i.e., DSP state. Let β be the wave vector of the DSP state, which has two components and in the x and y directions, respectively, thus
where and are the unit vectors in the x and y directions, respectively.After some straightforward algebra the specular reflection coefficient, , for a plane wave whose electric and magnetic fields have the form and impinging at the surface of the homogeneous but anisotropic layer of thickness h with and given by Eqs. (7) and (8) can be written as
here with being the refractive index of the dielectric medium. By extending this formula to the case and looking at the zeros of the denominator of and note that for this case, we can obtainwhere for, and for . For the plane wave whose electric and magnetic fields have the form and impinging at the surface of the homogeneous but anisotropic layer of thickness h with and given by Eqs. (12) and (13), the specular reflection coefficient, , can be obtained ashere . Note that for this case, similarly, we havewhere for , and for .From Eqs. (17) and (19), we can obtain the dispersion curves of the anisotropic DSP waves, as shown in Fig. 2 . Notably, the asymptotic frequencies corresponding to and are different, in contrast to an array of square holes where such two asymptotic frequencies are the same [5].
For a foil with arrays of holes of finite depth made on its one surface, when an incident wave with impinges at the surface, this wave will see only an average response that can be described by the effective anisotropic homogeneous medium layer as described above. If the foil is pierced with arrays of subwavelength apertures with size A and spacing D, as shown in Fig. 3 , the incident wave will passes through the apertures. The DSP states sustained by the holes could support the transmission resonances of incident wave through the subwavelength apertures. The dispersion relation for DSPs associated with this case is not well established. Nevertheless, it has been shown experimentally that the SP dispersion relation for a plane metal-dielectric interface without subwavelength apertures allows one to determine the locations of the transmission peaks when the interface is perforated at optical frequencies to good approximation. Therefore, we assume that such an approximation is equally valid here, thus the conservation of momentum is then given by [4]
for radiation normally incident of the interface, where j and l are zero and integers, and are the reciprocal lattice vectors for a square lattice with . From the two equations above, we find that the locations of transmission peak wavelengths are given bywhere , and . Note that the peak wavelengths depend on the geometrical and optical parameters of the pierced foils, especially on those of the holes, thus allowing us to design the peak wavelengths by adjusting the parameters of the holes.It is necessary to point out that Eq. (21) is obtained following the approach presented in Ref. 4, in which only the momentum-matching condition was considered. This equation hence does not take into account the boundary condition between the transmission region and inside the aperture. As a consequence, it neglects some physical processes [19], such as the interference that gives rise to a resonance shift [20].
3. Numerical results
In order to establish a concept of magnitude, we carried out the calculations of , and for , a, b, h, and d with and D = 500 μm. Figures 4 and 5 show the parameter dependence of these peak wavelengths. In Fig. 4, we use , and to normalize the values of , and , respectively. In fact, these parameters are those ones of , and at , corresponding to the peak wavelengths of the subwavelength square apertures without the rectangle holes. As can be seen from Fig. 4 (a), will increase from 708.30 μm to 708.48 μm for increasing from 1 to 3, having a wavelength shift , while the corresponding values for and being and , respectively, meaning that the higher the peak wavelength’s order is, the more sensitively the peak wavelength varies with . If we took a nematic liquid crystal as the material to fill the holes, we could adjust the peak wavelength via controlling the refractive index of the crystal by applying a magnetic field [21]. The biggest advantage lies in here is that there is no any materials filling the subwavelength square apertures for which the THz wave will pass through. For the other parameters (that is, a, b, h, and d), the higher the order of the peak wavelength is, the more sensitively the peak wavelength varies with these parameters, which is the same as that of the parameter. These numerical examples indicate that the peak wavelengths could be designed in a certain extent by adjusting the geometrical and optical parameters of the array of the rectangle holes, indicating that the DPS states sustained by the rectangle holes would affect the transmission property of electromagnetic wave passing through the subwavelength apertures.
4. Discussion and conclusion
The present experiments at THz regime cannot be used to check the theoretical results above because only those subwavelength apertures with the size and spacing being hundred-micron dimension were cut in freestanding 75 µm thick stainless steel foils [22]. Some suitable experiments are expected for which the samples should be pierced with subwavelength apertures and drilled with holes whose size and spacing should be taken as ten-micron dimension or less. Such experiments not only can be used to check the results above, but also would provide evidences for the existence of DSP states at THz regime. Our research proposes a way to connect the DSP theory to the engineering of surface plasmon optics, for which the creation of designer SPs with expectant resonant peaks could readily be realized by controlling the parameters of the holes, especially in THz to microwave regimes.
Acknowledgments
The Research Foundation of Wuhan National Laboratory under Grant No. P080008 and National “973” Project under Grant No. 2007CB310403 have supported this research.
References and links
1. S. A. Maier, Plasmonics—Fundamentals and Applications (Springer, New York, 2007).
2. C. Genet and T. W. Ebbesen, “Light in tiny holes,” Nature 445(7123), 39–46 (2007). [CrossRef] [PubMed]
3. T. W. Ebbesen, H. J. Lezec, H. F. Ghaemi, T. Thio, and P. A. Wolff, “Extraordinary optical transmission through sub-wavelength hole arrays,” Nature 391(6668), 667–669 (1998). [CrossRef]
4. H. F. Ghaemi, T. Thio, D. E. Grupp, T. W. Ebbesen, and H. J. Lezec, “Surface plasmons enhance optical transmission through subwavelength holes,” Phys. Rev. B 58(11), 6779–6782 (1998). [CrossRef]
5. J. B. Pendry, L. Martín-Moreno, and F. J. Garcia-Vidal, “Mimicking surface plasmons with structured surfaces,” Science 305(5685), 847–848 (2004). [CrossRef] [PubMed]
6. A. P. Hibbins, B. R. Evans, and J. R. Sambles, “Experimental verification of designer surface plasmons,” Science 308(5722), 670–672 (2005). [CrossRef] [PubMed]
7. F. J. Garcia-Vidal, L. Martin-Moreno, and J. B. Pendry, “Surfaces with holes in them: New plasmonic metamaterials,” J. Opt. A 7, S97–S101 (2005). [CrossRef]
8. J. T. Shen, P. B. Catrysse, and S. Fan, “Mechanism for designing metallic metamaterials with a high index of refraction,” Phys. Rev. Lett. 94(19), 197401–197404 (2005). [CrossRef] [PubMed]
9. D. R. Smith and J. B. Pendry, “Homogenization of metamaterials by field averaging,” J. Opt. Soc. Am. B 23(3), 391–403 (2006). [CrossRef]
10. M. Tanaka, F. Miyamaru, M. Hangyo, T. Tanaka, M. Akazawa, and E. Sano, “Effect of a thin dielectric layer on terahertz transmission characteristics for metal hole arrays,” Opt. Lett. 30(10), 1210–1212 (2005). [CrossRef] [PubMed]
11. C. R. Williams, S. R. Andrews, S. A. Maier, A. I. Fernández-Domínguez, L. Martín-Moreno, and F. J. García-Vidal, “Highly confined guiding of terahertz surface plasmon polaritons on structured metal surfaces,” Nat. Photonics 2(3), 175–179 (2008). [CrossRef]
12. L. Shen, X. Chen, and T. J. Yang, “Terahertz surface plasmon polaritons on periodically corrugated metal surfaces,” Opt. Express 16(5), 3326–3333 (2008). [CrossRef] [PubMed]
13. Y. C. Lan and R. L. Chern, “Surface plasmon-like modes on structured perfectly conducting surfaces,” Opt. Express 14(23), 11339–11347 (2006). [CrossRef] [PubMed]
14. S. A. Maier and S. R. Andrews, “Terahertz pulse propagation using plasmon-polariton-like surface modes on structured conductive surfaces,” Appl. Phys. Lett. 88(25), 251120 (2006). [CrossRef]
15. F. J. García de Abajo and J. J. Sáenz, “Electromagnetic surface modes in structured perfect-conductor surfaces,” Phys. Rev. Lett. 95(23), 233901 (2005). [CrossRef]
16. E. Hendry, A. P. Hibbins, and J. R. Sambles, “Importance of diffraction in determining the dispersion of designer surface plasmons,” Phys. Rev. B 78(23), 235426 (2008). [CrossRef]
17. Z. Ruan and M. Qiu, “Slow electromagnetic wave guided in subwavelength region along one-dimensional periodically structured metal surface,” Appl. Phys. Lett. 90(20), 201906 (2007). [CrossRef]
18. Q. Gan, Z. Fu, Y. J. Ding, and F. J. Bartoli, “Ultrawide-bandwidth slow-light system based on THz plasmonic graded metallic grating structures,” Phys. Rev. Lett. 100(25), 256803 (2008). [CrossRef] [PubMed]
19. C. Genet and T. W. Ebbesen, “Light in tiny holes,” Nature 445(7123), 39–46 (2007). [CrossRef] [PubMed]
20. C. Genet, M. P. van Exter, and J. P. Woerdman, “Fano-type interpretation of red shifts and red tails in hole array transmission spectra,” Opt. Commun. 225(4-6), 331–336 (2003). [CrossRef]
21. C.-L. Pan, C.-F. Hsieh, R.-P. Pan, M. Tanaka, F. Miyamaru, M. Tani, and M. Hangyo, “Control of enhanced THz transmission through metallic hole arrays using nematic liquid crystal,” Opt. Express 13(11), 3921–3930 (2005). [CrossRef] [PubMed]
22. H. Cao and A. Nahata, “Influence of aperture shape on the transmission properties of a periodic array of subwavelength apertures,” Opt. Express 12(16), 3664–3672 (2004). [CrossRef] [PubMed]