Tuning the extraordinary transmission in a metallic / dielectric CDC hole array by changing the temperature

Tunable extraordinary transmission via changing temperature of a porous metallic layer on top of a thin layer of dielectric strontium titanate (STO) was studied. The metallic layer has a through-hole array and each hole has a circular converging-diverging channel (CDC) shape, which induces the excitation of surface plasmon polaritons (SPPs) and then results in a controllable extraordinary optical transmission in the terahertz (THz) frequency range. We used a three-dimensional (3D) finite element method to analyze the transmission characteristics of the structure. Location and magnitude of the transmission peaks can be adjusted by hole size, converging angle, and thicknesses of metal and STO layers. Remarkably, the suggested structure presents a strong transmission dependency on temperature, which offers a new approach to actively and externally tune the transmission. This new design could lead to a family of temperaturesensitive devices working in the THz frequency range, promising in many applications including photonics, nanolithography, imaging, and sensing. ©2010 Optical Society of America OCIS codes: (000.2700) General science. References and links 1. M. Tonouchi, “Cutting-edge terahertz technology,” Nat. Photonics 1(2), 97–105 (2007). 2. R. Piesiewicz, T. Kleine-Ostmann, N. Krumbholz, D. Mittleman, M. Koch, J. Schoebel, and T. Kurner, “Shortrange ultra-broadband terahertz communications: Concepts and perspectives,” IEEE Ant. Prop. Mag. 49(6), 24– 39 (2007). 3. R. M. Woodward, B. E. Cole, V. P. Wallace, R. J. Pye, D. D. Arnone, E. H. Linfield, and M. Pepper, “Terahertz pulse imaging in reflection geometry of human skin cancer and skin tissue,” Phys. Med. Biol. 47(21), 3853–3863 (2002). 4. C. J. Strachan, P. F. Taday, D. A. Newnham, K. C. Gordon, J. A. Zeitler, M. Pepper, and T. Rades, “Using terahertz pulsed spectroscopy to quantify pharmaceutical polymorphism and crystallinity,” J. Pharm. Sci. 94(4), 837–846 (2005). 5. M. Nagel, P. H. Bolivar, M. Brucherseifer, H. Kurz, A. Bosserhoff, and R. Büttner, “Integrated THz technology for label-free genetic diagnostics,” Appl. Phys. Lett. 80(1), 154–156 (2002). 6. 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). 7. Q. Cao, and P. Lalanne, “Negative role of surface plasmons in the transmission of metallic gratings with very narrow slits,” Phys. Rev. Lett. 88(5), 057403 (2002). 8. P. Lalanne, J. P. Hugonin, and J. C. Rodier, “Theory of surface plasmon generation at nanoslit apertures,” Phys. Rev. Lett. 95(26), 263902 (2005). 9. H. J. Lezec, and T. Thio, “Diffracted evanescent wave model for enhanced and suppressed optical transmission through subwavelength hole arrays,” Opt. Express 12(16), 3629–3651 (2004). 10. D. X. Qu, D. Grischkowsky, and W. L. Zhang, “Terahertz transmission properties of thin, subwavelength metallic hole arrays,” Opt. Lett. 29(8), 896–898 (2004). 11. 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). 12. F. Miyamaru, and M. Hangyo, “Finite size effect of transmission property for metal hole arrays in subterahertz region,” Appl. Phys. Lett. 84(15), 2742–2744 (2004). 13. A. K. Azad, Y. Zhao, and W. Zhang, “Transmission properties of terahertz pulses through an ultrathin subwavelength silicon hole array,” Appl. Phys. Lett. 86(14), 141102 (2005). #129836 $15.00 USD Received 9 Jun 2010; revised 20 Jun 2010; accepted 20 Jun 2010; published 7 Jul 2010 (C) 2010 OSA 19 July 2010 / Vol. 18, No. 15 / OPTICS EXPRESS 15553 14. J. G. Rivas, C. Schotsch, P. H. Bolivar, and H. Kurz, “Enhanced transmission of THz radiation through subwavelength holes,” Phys. Rev. B 68, (2003). 15. C. Janke, J. G. Rivas, C. Schotsch, L. Beckmann, P. H. Bolivar, and H. Kurz, “Optimization of enhanced terahertz transmission through arrays of subwavelength apertures,” Phys. Rev. B 69(20), 205314 (2004). 16. 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). 17. J. B. Pendry, L. Martín-Moreno, and F. J. Garcia-Vidal, “Mimicking surface plasmons with structured surfaces,” Science 305(5685), 847–848 (2004). 18. A. Degiron, H. J. Lezec, W. L. Barnes, and T. W. Ebbesen, “Effects of hole depth on enhanced light transmission through subwavelength hole arrays,” Appl. Phys. Lett. 81(23), 4327–4329 (2002). 19. A. K. Azad, and W. L. Zhang, “Resonant terahertz transmission in subwavelength metallic hole arrays of subskin-depth thickness,” Opt. Lett. 30(21), 2945–2947 (2005). 20. 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). 21. J. M. Steele, Z. W. Liu, Y. Wang, and X. Zhang, “Resonant and non-resonant generation and focusing of surface plasmons with circular gratings,” Opt. Express 14(12), 5664–5670 (2006). 22. K. J. K. Koerkamp, S. Enoch, F. B. Segerink, N. F. van Hulst, and L. Kuipers, “Strong influence of hole shape on extraordinary transmission through periodic arrays of subwavelength holes,” Phys. Rev. Lett. 92(18), 183901 (2004). 23. A. Battula, Y. L. Lu, R. J. Knize, K. Reinhardt, and S. C. Chen, “Tunable transmission at 100 THz through a metallic hole array with a varying hole channel shape,” Opt. Express 15(22), 14629–14635 (2007). 24. W. Wang, Y. L. Lu, R. J. Knize, K. Reinhardt, and S. C. Chen, “Tunable and polarization-selective THz range transmission properties of metallic rectangular array with a varying hole channel shape,” Opt. Express 17(9), 7361–7367 (2009). 25. A. Battula, S. Chen, Y. Lu, R. J. Knize, and K. Reinhardt, “Tuning the extraordinary optical transmission through subwavelength hole array by applying a magnetic field,” Opt. Lett. 32(18), 2692–2694 (2007). 26. P. Kuzel, and F. Kadlec, “Tunable structures and modulators for THz light,” C. R. Phys. 9(2), 197–214 (2008). 27. A. K. Tagantsev, V. O. Sherman, K. F. Astafiev, J. Venkatesh, and N. Setter, “Ferroelectric materials for microwave tunable applications,” J. Electroceram. 11(1/2), 5–66 (2003). 28. E. D. Palik, Handbook of Optical Constants of Solids (Academic, 1985). 29. A. Lavrinenko, P. I. Borel, L. H. Frandsen, M. Thorhauge, A. Harpøth, M. Kristensen, T. Niemi, and H. M. H. Chong, “Comprehensive FDTD modelling of photonic crystal waveguide components,” Opt. Express 12(2), 234–248 (2004). 30. H. A. Bethe, “Theory of diffraction by small holes,” Phys. Rev. 66(7-8), 163–182 (1944).


Introduction
The terahertz (THz) frequency band, which lies between domains of microwave electronics and mid infrared optics, has received increasing attentions recently because of its possible applications in fundamental science, new imaging and sensing modalities, and high bandwidth signal processing [1].Examples of such applications include communication [2], medical imaging [3], pharmaceutical quality control [4], and chemical and biochemical sensing [5].Utilizing the light in the THz band, however, requires a precise control of the wave propagation as well as the electromagnetic field distribution.For example, in the applications as THz filters, photonic circuits, and sensors, control of the transmission spectrum, especially a capability of active and precise frequency tuning, is highly desirable.
The extraordinary optical transmission (EOT) phenomenon was first reported by T. W. Ebbesen et. al. in 1998 [6].This work has generated considerable interests in recent years and led to the development of a new family of optical devices based on subwavelength hole arrays with extraordinary transmission spectra.This unique phenomenon is believed to be due to the SPPs excited at the metal/dielectric interface [7][8][9].Although metal surfaces in the THz frequency region act as perfect conductors making SPPs excitation unexpected, EOT has been observed recently in both metallic [10][11][12] and semiconductor [13][14][15] hole arrays, and it has been proven both experimentally and theoretically that the resonantly-excited SPP-like mode plays an important role in the much enhanced transmissions [16,17].Excited by Ebbesen's pioneer work, intensive studies were done to analyze the tunability of the transmission.By changing the geometric variables (film thickness, interface media, periodicity, hole shape, aperture size, channel shape, etc.), the transmission properties (magnitude, location of the peaks, the full width at half maximum (FWHM), and selectivity of different polarization) can be tuned correspondingly [10,11,[18][19][20][21][22][23][24].
However, tuning by changing the geometry variables indicates either deforming or replacing the device will be required, and is not active.This actually limits their future applications in real-time systems.Previous work has shown that the transmission of light through a thin Ag film with a specifically designed periodic subwavelength hole array can be influenced by the presence of the externally applied magnetic field [25], which indeed indicates an initial move toward the active frequency tuning direction.Considering the fact that SPPs strongly rely on the metal/dielectric interface, modification of the interface would become an effective route for realizing the active tuning, and recent development of dielectric materials with a strong non-linear response [26] will actually support the idea.Here, the term "non-linear response" means their optical properties can be controlled by external parameters (temperature, electric or magnetic field, light pulse, etc.).Among these non-linear materials, so called incipient ferroelectric materials [27] such as strontium titanate (SrTiO3, STO) or potassium tantalate (KTaO3, KTO) exhibit a potential to realize the large frequency range and active tuning.For example, between 0.3 THz to 0.9 THz, i.e. wavelengths between 300 µm to 1 mm, both real and imaginary parts of permittivity of STO undergo an increase as temperature decreases [26].This unique property offers a new way to design strong temperature-dependent devices working in the THz band.
In this work, we propose a new structure with a metallic layer on top of a thin layer of the STO.The metallic layer has a through-hole array with a circular, converging-diverging channel (CDC) shape.Use of the CDC shape has been proven with similar EOT effect, in addition, with an extra degree of freedom in a geometric variable to further fine-tune the transmission spectrum [23].Simulations were performed to study the EOT properties under different temperature conditions.By careful design of the converging angle and STO layer thickness, the proposed structure can function well at certain THz bands with a realistic transmission spectrum profile.More importantly, the strong transmission dependency on temperature provides us the capability to modulate the transmission without using deformation or replacement-a true active frequency tuning.

Simulation considerations
Figure 1 shows a side view of the proposed structure used for the simulation where the incident light transmits through the thin layers of silver and STO.Their thicknesses are t 1 and t 2 , respectively.There is a circular CDC hole array with a semi converging angle θ in the silver layer [23].The radii of the holes at surface are fixed at 20µm.The transmission was calculated for different wavelengths, temperatures and geometric variables.Figure 2 shows the top view of the silver film.The circular hole array is periodic in both the x and y directions, and the periodicity for the two directions are the same (50µm).The inset shows the details of a single circular hole.
Fig. 1.The side view of the device structure.The thickness for the metallic layer is t1, the thickness for the dielectric layer with a temperature dependent dielectric constant is t2, and a circular CDC hole array with a semi-angle θ will be made inside the metallic layer.
The converging angle, θ, and the metal and STO thickness, t 1 and t 2 , were changed accordingly in the simulation, while keeping all other parameters unchanged.The goal is to study the effect of such varying parameters on the relationship among transmission, wavelength, and temperature.For this study we considered the frequency to be around 0.6THz or 450µm wavelength and the temperature range to be from 270K to 305K.The dielectric constant of STO was described by a damped harmonic oscillator model: [26].The dielectric constant of silver used in the simulation was described by the Drude model ε = ε ∞ -ω p 2 / (ω 2 + ݅γ s ω), where ε ∞ = 1, ω p = 1.44 × 10 16 s −1 , and γ s = 4 × 10 13 s −1 [28].Within the suggested temperature range in our calculation, the dielectric constant of silver is set to be temperature independent.The electromagnetic fields were assumed to be time harmonic and the resulting governing equations for the steady-state distribution was solved using commercially available three-dimensional (3D) finite element software (COMSOL 3.3).The computational domain considered is a single unit cell surrounded either by periodic boundary conditions or by perfectly matching layers (PML) [29].The light is incident normal to the film surface and the transmittance is calculated from the obtained electromagnetic field distributions.

Results and discussion
Figure 3 shows the transmission vs. temperature and wavelength for the structure having various geometric variables.To better evaluate the transmission property, all figures in Fig. 3 were plotted at the same scale.All the figures shown are characterized by their strong transmission peaks which are much larger than (A/λ) 4 , the Bethe-Bouwkamp power law model [30].In this model, 'A' is the aperture size and 'λ' is the wavelength of the light.This indicates the existence of EOT inside the proposed structure.Further examination reveals that the transmission spectrum changes as both temperature and geometry change.Let's first look at the transmission's dependency on temperature.Each figure in Fig. 3 shows the full transmission spectrum under certain geometric conditions.Since a higher temperature results in both smaller imaginary and real parts of the dielectric constant of the STO material, the magnitude of the transmission peak increases and the peak blue-shifts as the temperature increases.
To better study the role of changing both temperature and converging angle, we compare Fig. 3a-c altogether.Responses from both CDC and the straight-hole array due to a temperature change of 36 K are re-plotted in Fig. 4. When changing the temperature 36 K around the room temperature as shown in Fig. 4a, the transmission peak's location and magnitude among all cases change about 33 µm, which is about 3% central to the transmission wavelength.This indicates a fairly large frequency tunable range, and is actually practical for many long wavelength photonics application.Another remarkable benefit is that these transmission properties are nearly linearly dependent on temperature.Moreover, the full width at half maximum (FWHM) is relatively independent of temperature change, which remains to be around 60 µm.
For the hole-arrayed metal/dielectric structure, an increase in the angle θ actually gives rise to a decrease in transmission, this is slightly different from the case without the dielectric layer [23].However, it shows a blue-shifted transmission when increasing the angle θ, indicating an additional manipulation over the transmission by changing the angle is possible.Figure 4b shows the different transmission spectra at different converging angles at 269 K.As the apertures in the metal layer change from straight-channel to the CDC shape with a converging angle of 60°, the peak transmission reduces from 20.1% to 15.8% with a blue shift of ~32 µm, and the FWHM of those peaks is actually becoming narrower too.The results here indicate, for the suggested metal/dielectric structure, introducing CDC to the design is not beneficial to the transmission amplitude, but does add a freedom in fine-adjusting the peak position and in modifying the FWHM of the transmission.
The comparison of Fig. 3b and 3d reveals the transmission dependency on the metal film thickness.With the same converging angle and temperature, the magnitude of the transmission peak is much lower for a thicker metal layer.Also, as the metal layer thickness increases, the location of the transmission peak red-shifts at the same time.It is important to point out that Fig. 3c (using a larger converging angle) and 3d (using a thicker metal layer thickness) show very closed transmission spectrum, which indicates that adding the CDC design to the structure can provide not only an extra degree of freedom to fine-tune the transmission properties, but also as an important substitution to change other geometric variables, for example, here to use less amount of metallic material.
The thickness of STO the layer also plays an important role in the transmission properties.Figure 3b, 3e and 3f show the transmission spectrum with a same converging angle of 45° for different STO thicknesses.Although all the transmission peaks in the three plots are still much larger than those predicted by the Bethe-Bouwkamp power law model, the magnitude of each peak is much lower for a thicker layer of STO.Along with magnitude reduction, the location of the transmission peak red-shifts at the same time.Since STO is not a lossless material, the reduction of transmission results from the increased energy loss in longer travel distance in the absorbing material.In the case of our simulation, as shown in Fig. 3, when the STO layer thickness increases from 1.5 µm to 2.5 µm, a steep magnitude drop of the transmission peaks of around 10% with a red shift of nearly 100 µm is observed.The dramatic transmission change suggests that the STO layer thickness is a crucial parameter on which special attention is needed to meet the desired performance requirement.A further investigation over the case having holes all the way through the metal/dielectric structure will be also interesting, and this is an ongoing effort.

Conclusion
We have investigated the transmission spectra of metallic hole arrays with different converging-diverging channels combined with a layer of STO.Besides extraordinarily-high subwavelength transmission, the simulation results also show strong tunable transmission characteristics.Both the location and magnitude of the transmission peak can be tuned by temperature, converging angle, and both of the metal and STO film thickness.Moreover, adding a STO layer offers a new approach to actively control the transmission without mechanically changing the optical device.As a result, this proposed structure could lead to extraordinary transmission at different wavelengths and can be used to develop temperaturetuning THz filters.

Fig. 2 .
Fig. 2. The top view of the structure.The circular hole array is periodic in both the x and y directions, and the periodicity for the two directions are the same (50µm).The inset shows the top view of a single circular hole.

Fig. 3 .
Fig. 3. Transmittance spectrum with different geometric variables and temperatures.(a)-(c) are plots of transmission versus wavelength and temperature with converging angle (a) 0°, (b) 45°, and (c) 60°.Thickness of the metal and STO layer are both fixed at 2µm.(d) shows the case at which the thickness of the metal and STO are 3µm and 2µm, respectively.The converging angle is 45°.(e) and (f) show the transmittance spectrums with the same metal film thicknesses (2µm) and converging angle (45°), but with different STO layer thicknesses: (e) t2 = 1µm and (f) t2 = 2.5µm.For all cases above, the hole area is unchanged.The transmission peak vs. temperature curve is shown at the bottom in each figure (red curve).

Fig. 4 .
Fig. 4. Transmittance spectrum curves.(a) shows transmission spectrum with different temperatures.The converging angles are fixed at 0° (straight hole) for all cases.(b) shows transmission spectrum with different converging angles.The temperatures are fixed at 269 K for all cases.