Light transmission by subwavelength square coaxial aperture arrays in metallic films

Using the Fourier modal method, we study the enhanced transmission exhibited by arrays of square coaxial apertures in a metallic film. The calculated transmission spectrum is in good agreement w ith FDTD calculations. We show that the enhanced transmission can be explained when we consider a few guided modes of a coaxial waveguide. © 2003 Optical Society of America OCIS codes:(050.1120) Apertures; (240.6680) Surface plasmons; (230.7 37 ) Waveguides References and links 1. T. W. Ebbesen, H. J. Lezec, H. F. Ghaemi, T. Thio, and P. A. Wol f, “Extraordinary optical transmission through sub-wavelength hole arrays,” Nature 391, 667 (1998). 2. J. A. Porto, F. J. Garcia-Vidal, and J. B. Pendry, “Transmis ion resonances on metallic gratings with very narrow slits,” Phys. Rev. Lett. 83, 2845 (1999). 3. L. Martin-Moreno, F. J. Garcia-Vidal, H. J. Lezec, K. M. Pe llerin, T. Thio, J. B. Pendry, and T. W. Ebbesen, “Theory of extraordinary optical transmission through subw avelength hole arrays,” Phys. Rev. Lett. 86, 114 (2001). 4. Q. Cao and P. Lalanne, “Negative role of surface plasmons in the transmission of metallic gratings with very narrow slits,” Phys. Rev. Lett. 88, 057403 (2002). 5. S. Astilean, P. Lalanne, and M. Palamaru, “Light transmissi on through metallic channels much smaller than the wavelength,” Opt. Commun. 175, 265 (2000). 6. M. M. J. Treacy, “Dynamical diffraction in metallic optical gratings,” Appl. Phys. Lett. 75, 606 (1999). 7. M. M. J. Treacy, “Dynamical diffraction explanation of the anomalous transmission of light through mettallic gratings,” Phys. Rev. B66, 195105 (2002). 8. S. Collin, F. Pardo, R. Teissier, and J.-L. Pelouard, “Hor izontal and vertical surface resonances in transmission metallic gratings,” J. Opt. A:Pure Appl. Opt. 4, S154 (2002). 9. L. Li, “New formulation of the Fourier modal method for corsse d surface-relief gratings,” J. Opt. Soc. Am. A 11, 2758-2767 (1997). 10. E. Popov, M. Nevì ere, S. Enoch, and R. Reinisch, “Theory of light transmission through subwavelength periodic hole arrays,” Phys. Rev. B 62, 16100 (2000). 11. F. I. Baida and D. Van Labeke, “Light transmission by subwa velength annular aperture arrays in metallic films,” Opt. Commun.209, 17-22 (2002). 12. G. Granet and J. P. Plumey, “Parametric formulation of the Fo urier modal method for crossed surface-relief gratings,” J. Opt. A: Pure Appl. Opt. 4, S145 (2002).


Introduction
Nowadays, opticians are greatly interested in structures that exhibit anomalous effects, because they have potential applications in novel photonic devices.The extraordinary enhanced transmission by subwavelength metallic hole arrays is one such phenomenon.Since the publication of by Ebbesen et al., 1 many experimental and theoretical studies were carried out in order to determine the physical origin of the observed enhanced transmission.Three kinds of explanations have been proposed since, relating the enhanced transmission to the excitation of surface plasmons, 2,3 to a Fabry-Pérot cavity behavior of the holes, 4,5 or explaining the transmission in terms of dynamical 6,7 diffraction.It is now established that both horizontal and vertical resonances play a role 8 in the extraordinary transmission.It is then of importance to characterize and to understand the electromagnetic behavior of the channel through which the light propagates inside the metallic film. 10Recently, numerical simulations have shown that a transmission as high as 80% can be obtained with anular 11 apertures.The aim of the present communication is to study the spectral response of metallic films with a square coaxial aperture.Those structures are similar to the above-mentioned ones from the electromagnetic point of view.Since the aperture dimensions are of the order of magnitude of the wavelength, a rigorous electromagnetic theory is necessary to analyze the behavior of such structures.Although the FDTD method allows us to calculate rigorously the reflection and transmission of a plane wave by a periodical structure in the resonance domain, the Fourier modal method gives a more physical insight in the present resonant phenomenon. 12The diffraction problem is reduced to the searching of eigenvalues and eigenvectors of a particular matrix.It permits us to calculate the effective index of the modes of the coaxial aperture and the coupling of these modes with the reflected and transmitted order.

Statement of the problem
Let us consider a metallic film deposed on a glass substrate with an engraved periodic structure of square coaxial apertures (see Fig. 1).The refractive index of the metal is described by a simple free-electron Drude model with a plasma frequency ω p = 1.374 × 10 16 s −1 and a relaxation time τ = 0.3 × 10 −14 s.The periods are d x in the x direction and d y in the y direction.The width and the position of an aperture are controlled by two parameters w 1 and w 2 (see Fig. 1).Finally, the thickness is denoted by h.The structure is illuminated in vacuum, under normal incidence, by a monochromatic linearly polarized plane wave, with a wavelength λ , a wavenumber k = 2π/λ , and a time dependence exp(iωt).Our goal is to calculate and to understand the reflection and transmission spectra of this structure with the help of the Fourier modal method.In the layer, any component F of the electric or magnetic field can indeed be expressed as a superposition of eigenmodes: where m and n are integers such that −M ≤ m ≤ M and −N ≤ n ≤ N. The integers M and N describe the truncation scheme.The matrix from which eigenvalues and eigenvectors are calculated is then of rank 2(2M + 1)(2N + 1).A + q and A − q are the unknown complex amplitudes of the upward and downward propagating or decaying waves.Our numerical code includes the correct factorization rules derived by Li, 9 our personal parametric formulation, and the S matrix approach for writing the boundary conditions.It should be emphasized that the abovementioned numerical tools are of great importance for obtaining reliable and converged results, although this point is beyond the scope of our paper.To compare the Fourier modal method the FDTD that was used by Baida and Van Labeke, 11 we have calculated the transmission spectrum of a structure with the following parameters: It can be seen that both methods give resonances at the same place even though a small difference is observed in their intensity.

Analysis of the mode
Our goal is to analyze the enhanced transmission by using the guided modes of the coaxial apertures.Since we consider a metallic medium, an aperture is not coupled with its neighbors.A mode for the entire structure thus corresponds exactly to a mode of a sole aperture, and thus no distinction is made in this discussion between them.Indeed, the eigenvalues and the fields inside the apertures corresponding to an eigenmode do not change when the distance between holes varies.As a consequence, the eigenvalues γ q give an immediate access to the effective index for each guided mode.
Because of the metal, all the propagating constants are complex but some of them can be considered as guided modes with low losses.For the considered structure we have found that there were three such modes, two of them being degenerated as a result of the square symmetry.The numerically obtained dispersion relations are plotted in Figs. 3 and 4  Figures 5 and 6 show a map of the modulus of the transverse electric field of the first and the second modes.The mode whose effective index has the largest real part and the lowest imaginary part corresponds to the TEM mode of the same coaxial structure with perfect conducting walls.This mode is characterized by an electric field normal to the walls and has no cut-off.In the present case, it is not strictly speaking a TEM mode since its effective index is greater than one.However, when the width of the aperture becomes larger, the coupling between the opposite sides of the coaxial waveguide diminishes resulting in a lower effective index.The two other guided modes have a cut-off ∼λ = 845 nm.

Analysis of the coupling of the modes to free radiation
We have shown the existence of attenuated guided modes.In the present section we are interested in the way they can be excited by an incident plane wave.The S matrix approach is a very appropriate tool for such an analysis.Let us consider the particular wavelength λ = 558 nm where a resonance occurs.Figures 7 and 8 show the 21st calculated modal coefficients corresponding to the upward and downward waves inside the aperture when the film is illuminated by an x-polarized plane wave under normal incidence.For convenience, sorting is displayed in decreasing order.To obtain some physical information from this spectrum analysis, we have carefully normalized all the eigenvectors.It should be noted that the above coefficients are calculated on the interface where the corresponding wave has been excited.By considering the eigenvalues, i.e., the normalized propagating constants, one can easily deduce which kind of mode is excited.Figure 9 represents the location in the complex plane of the propagating constant associated to the modal amplitude of Fig. 7.In Fig. 7, the first and the third modes are a degenerated mode whose the imaginary part of the effective index is as high as 22.8.In the present case, we can conclude that the mode responsible for the resonant transmission is the attenuated guided mode that matches the polarization of the incident wave.This mode has an effective index of 1.39 − 0.006i. .The ten first complex propagating constants associated to the ten first modes that are exited on the upper face inside the coaxial waveguide.The red one corresponds to an attenuated guide wave; its value is γ = 1.39 − 0.006i.

Conclusion
We have numerically studied the spectral response of subwavelength coaxial apertures.We have calculated the propagating constants of the modes supported by a square coaxial waveguide.Some of them correspond to attenuated guided modes.However, the excitation of such modes is possible only when the incident wave matches the mode profile.Owing to the electric properties of metals at optical wavelengths, the dispersion relations of the modes of the transmission channel are very specific and very different from those of the same channel with perfectly conducting walls.This preliminary study paves the way for future investigations in order to engineer the modes and their excitation for applications.

Fig. 2 .
Fig.2.Transmission of a square coaxial aperture calculated with the FDTD (blue curve) and the Fourier modal method (red curve).

Fig. 4 .
Fig. 4. Dispersion curves of the second mode.Blue curve, real part; red curve, imaginary part.

Fig. 5 .Fig. 6 .
Fig. 5. Modulus of the transverse electric field of the first guided mode.

Fig. 7 .
Fig. 7. Twenty-first modal amplitude coefficients inside the coaxial on the upper face.The red bar corresponds to an attenuated guided wave.

Fig. 8 .
Fig. 8. Twenty-first modal amplitude coefficients inside the coaxial on the lower face.The red bar corresponds to an attenuated guided wave.

Fig. 9
Fig.9.The ten first complex propagating constants associated to the ten first modes that are exited on the upper face inside the coaxial waveguide.The red one corresponds to an attenuated guide wave; its value is γ = 1.39 − 0.006i.
. Fig. 3. Dispersion curves of the first mode.Blue curve, real part; red curve, imaginary part.The presence of dips is probably due to the right angle corners.