## Transmission of light through periodic arrays of sub-wavelength slits in metallic hosts

Optics Express, Vol. 14, Issue 14, pp. 6400-6413 (2006)

http://dx.doi.org/10.1364/OE.14.006400

Acrobat PDF (487 KB)

### Abstract

Using Bloch modes to study the extraordinary transmission of light through a periodic array of slits in a metallic host, we discuss the differing roles of surface plasmon polaritons and Wood’s anomalies in the observed behavior of such structures. Under certain circumstances, the first few excited modes appear to play a decisive role in determining the transmission efficiency of the array. Surface plasmon excitations tend to reduce the transmissivity of a semi-infinitely thick slit array, yet, paradoxically, the same reduction can account for enhanced transmission in an array of finite thickness *τ*, provided that *τ* is tuned to a Fabry-Perot-like resonance between the entrance and exit facets of the slit array. At the Wood anomaly, power redistribution produces sharp peaks in the diffraction efficiencies of various reflected and transmitted orders of the semi-infinite structure. With skew incidence, the degenerate states split, resulting in two peaks and two valleys, as observed by Wood in his 1902 experiments.

© 2006 Optical Society of America

## 1. Introduction

1. T. W. Ebbesen, H. J. Lezec, H. F. Ghaemi, T. Thio, and P. A. Wolff, “Extraordinary optical transmission through subwavelength hole arrays,” Nature **39**, 667–669 (1998). [CrossRef]

11. Y. Xie, A.R. Zakharian, J. V. Moloney, and M. Mansuripur, “Transmission of light through slit apertures in metallic films,” Opt. Express **12**, 6106 (2004). [CrossRef] [PubMed]

2. M. M. J. Treacy, “Dynamical diffraction explanation of the anomalous transmission of light through metallic gratings,” Phys. Rev. B **66**, 195105–11 (2002). [CrossRef]

12. R. W. Wood, “On a remarkable case of uneven distribution of light in a diffraction grating spectrum,” Proc. Phys. Soc. London **18**, 269–275 (1902). [CrossRef]

13. R. W. Wood, “Anomalous diffraction gratings,” Phys. Rev. **48**, 928–937 (1935). [CrossRef]

*E*-and

*H*-fields uniquely specifies the common (tangential) component of the

*k*-vector for these two beams as

*k*

_{o}=2

*π*/λ

_{o}is the propagation constant in vacuum, and

*ε*

_{m}and

*ε*

_{d}are the (relative) permittivities of the metal and the dielectric, respectively. (The surface plasmon resonance at

*ω*

_{sp}. While the resonant SP, a collective oscillation of all surface electrons, has a group velocity of zero, the SPP is a localized excitation that propagates with a finite group velocity as a guided surface mode.)

*modes*(instead of Bloch

*waves*) is that, inside the periodic structure, such modes exist and propagate independently of all the other modes. In a semi-infinite slit array located in the half-space

*z*>0 and illuminated by a plane-wave arriving from the free space region

*z*<0, the strength of the various excited modes can be determined by matching the boundary conditions at the entrance facet of the slit array, namely, in the

*xy*-plane located at

*z*=0. This is done in Section 3, where we minimize (simultaneously) the difference between the tangential

*E*-fields, as well as that between the tangential

*H*-fields, across the interface. The (complex) amplitudes of all the excited modes are thus determined and, with enough modes included in the calculation, the boundary conditions at

*z*=0 are satisfied.

*p*of the array (all other parameters being kept constant), the excitation of SPP and the occurrence of Wood’s anomaly will be seen to be vastly different phenomena insofar as the number and strength of excited Bloch modes are concerned. In Section 4, we analyze the modal content of reflected and transmitted beams, identify dominant modes, and explain their roles in enhancing or weakening the transmission under conditions that either favor the excitation of SPP or usher in the Wood anomalies. Section 5 provides a brief analysis of the SPP and Wood anomalies at oblique incidence. Final remarks and conclusions appear in Section 6.

## 2. Bloch modes of a periodic array of slits in a semi-infinite metallic host

*z*>0. The array is illuminated by a plane-wave from the free-space region above the surface, which surface is the

*xy*-plane at

*z*=0. The host material is silver, the slits are empty (i.e., vacuum or air-filled), the incident wavelength is

*λ*

_{o}, the array periodicity along the

*y*-axis is

*p*, and the slit-width is

*w*. We shall consider only the case of incident plane-waves whose propagation vector

*k*lies in the

*yz*-plane; these will be assumed to be linearly polarized within the plane of incidence; in other words, the incident plane-waves are transverse-magnetic (TM) or

*p*-polarized. Although, for the most part, this paper will be concerned with a normally incident plane-wave polarized along the

*y*-axis, (

*E*

_{y},

*H*

_{x}), generalizing the method of analysis to oblique incidence is fairly straightforward; some numerical results pertaining to the case of oblique incidence are presented in Section 5.

*z*>0, while ordinary homogeneous (i.e., propagating) and inhomogeneous (i.e., evanescent) plane-waves reside in the free-space region above the slit array, where

*z*<0. Our construction of the Bloch modes of the slit-array ensures that Maxwell’s equations hold both within the slits and in the metallic host, with boundary conditions satisfied on the slits’ vertical walls. Matching the boundary conditions at the top surface,

*z*=0, leads to a determination of the (complex) amplitudes of the various excited modes that reside above and below the interface.

*et al*. [15

15. Ph. Lalanne, C. Sauvan, J. P. Hugonin, J. C. Rodier, and P. Chavel, “Perturbative approach for surface plasmon effects on flat interfaces periodically corrugated by subwavelength apertures,” Phys. Rev. B **68**, 125404 (2003). [CrossRef]

16. Ph. Lalanne, J. P. Hugonin, and J. C. Rodier, “Theory of Surface Plasmon Generation at Nanoslit Apertures,” Phys. Rev. Lett. **95**, 263902 (2005). [CrossRef]

*ε*

_{m}, and will include as many Bloch modes as needed to achieve convergence. In this section, we express the

*E*-and

*H*-fields inside the slits and within the metallic host as superpositions of inhomogeneous plane-waves, and solve the transcendental equation that results from the continuity of

*H*

_{x}and

*E*

_{z}on the slit walls. At normal incidence, the equations are somewhat simplified due to symmetry.

### 2.1 Transcendental equation for Bloch modes of the slit array

*z*<0 above the slit array are propagating as well as evanescent plane-waves (i.e., Rayleigh expansion, similar to RCWA method [17

17. M. G. Moharam, E. B. Grann, D. A. Pommet, and T. K. Gaylord, “Formulation for stable and efficient implementation of the rigorous coupled-wave analysis of binary gratings,” J. Opt. Soc. Am. A **12**, 1068–76 (1995). [CrossRef]

18. Ph. Lalanne and G. M. Morris, “Highly improved convergence of the coupled-wave method for TM polarization,” J. Opt. Soc. Am. A **13**, 779–84 (1996). [CrossRef]

*E*

_{x},

*H*

_{y},

*H*

_{z}) and TM (

*H*

_{x},

*E*

_{y},

*E*

_{z}) modes. The sole concern of the present paper is TM modes, as TE modes do not propagate through slits narrower than a half-wavelength. We express the

*E*-and

*H*-fields inside the slits–which are either empty or filled with a transparent material of dielectric constant

*ε*

_{s}-and also the fields in the metallic regions between adjacent slits as superpositions of two (generally inhomogeneous) plane-waves bouncing back and forth between the vertical walls. The

*H*-field transmitted below the surface at

*z*=0 is thus written

*n*is the mode index (1, 2, 3, …),

*k*

_{o}=2

*π*/

*λ*

_{o}is the vacuum wave-number,

*n*

^{th}mode’s propagation constant along the

*z*-axis, the subscript

*s*denotes the slit region, and the subscript

*m*denotes the metallic host medium (i.e., cladding for the slit waveguides). Inside the slits (

^{2}+(

^{2}=

*ε*

_{s}, where

*ε*

_{s}is the relative permittivity of the filling material, while in the metal (

^{2}+(

^{2}=

*ε*

_{m}, where

*ε*

_{m}is the relative permittivity of the host material. The horizontal component

*E*

_{y}and the vertical component

*E*

_{z}of the

*E*-field are proportional to the partial derivatives of

*H*

_{x}with respect to

*z*and

*y*, respectively:

*ε*

_{o}is the vacuum permittivity, and

*ω*=2

*πc*/

*λ*

_{o}is the angular frequency of the incident beam. Defining Λ=exp(i

*k*

_{o}

*p*sin

*θ*), where

*θ*is the incident plane-wave’s angle of incidence, the Bloch condition requires that

*E*-and

*H*-fields at the vertical walls of the slits yields:

*H*-field amplitudes (

*h*

_{1s},

*h*

_{2s},

*h*

_{1m},

*h*

_{2m}) of the

*n*

^{th}mode, namely,

*a*=exp(i

*k*

_{o}

*b*=exp[i

*k*

_{o}

*σ*

_{ym}(

*p*-

*w*)]. Existence of a non-trivial solution for (

*h*

_{1s},

*h*

_{2s},

*h*

_{1m},

*h*

_{2m}) requires that the determinant of the coefficient matrix be zero. For each mode, therefore, the propagation constant

*h*

_{1s}=

*h*

_{2s}and

*h*

_{1m}=

*h*

_{2m}, and the characteristic equation simplifies as follows:

### 2.2 Bloch mode profiles

*σ*

_{z}. We enumerate the solutions in such a way that lower-order modes will lose less power as they propagate along the

*z*-axis; that is,

*σ*

_{z}will be ordered in accordance with the strength of its imaginary part. If the solutions include guided modes, these will constitute the least lossy modes of all.

*p*=0.9µm,

*w*=0.1µm at

*λ*

_{o}=1.0µm (normal incidence), we find

*λ*

_{o}) along the

*z*-axis. For this system, Fig. 2(a) shows the locations of

*σ*

_{z}for modes 2 through 30 in the complex plane; the corresponding plots for

*σ*

_{ys}and

*σ*

_{ym}are shown in Figs. 2(b, c); the arrows identify the direction of increasing mode number, from 2 to 30.

*σ*

_{z}is determined, the mode profiles

*H*

_{x}(

*y*,

*z*),

*E*

_{y}(

*y*,

*z*),

*E*

_{z}(

*y*,

*z*) are readily calculated. Figure 3 shows profiles |

*H*

_{x}(

*y*,

*z*=0

^{+})| and |

*E*

_{y}(

*y*,

*z*=0

^{+})| of the first 10 Bloch modes for the

*p*=0.9 µm,

*w*=0.1 µm slit array under normally-incident plane-wave (

*λ*

_{o}=1.0 µm). Each mode, of course, will decay exponentially along the propagation direction

*z*in accordance with the factor exp(i

*k*

_{o}

*z*). Similar profiles can be computed for other sets of parameter values.

*n*+1)

^{st}mode typically has one more wiggle than the

*n*

^{th}mode in the metal on either side of the slit. Inside the slit, however, an additional wiggle appears only once every several modes. With an increasing period

*p*, determination of the modes becomes more difficult as more roots

*σ*

_{z}of the transcendental Eq. (6) become available per unit area of the complex plane. When

*w*<

*λ*

_{o}/2, the first mode can propagate deep into the slits, while all other modes exhibit large absorption along

*z*. The mode amplitudes thus drop rapidly after propagating even a short distance; for higher-order modes this distance is much smaller than the skin-depth of the metallic host. The modes typically extend along

*y*in the host material with fairly small attenuation, as can be inferred from the complex-plane distribution of

*σ*

_{ym}depicted in Fig. 2(c). These modes build up the charges and currents on the top surface of the metal, and contribute to the accumulated charges at the corners of the slits.

### 2.3 Orthogonality of Bloch modes

*not*mutually orthogonal [20]. This can be seen in Fig. 4, which shows the magnitudes of the inner products of the first 60 modes for both

*H*

_{x}and

*E*

_{y}fields. (Inner products are defined as ∫

*y*, 0)

*H*

^{m}

^{*}

_{x}(

*y*, 0) dy and ∫

*y*, 0)

*E*

^{m}

^{*}

_{y}(

*y*, 0) d

*y*, with individual modes doubly normalized, i.e., ∫ |

*y*, 0)|2 dy=1 and ∫ |

*y*, 0)|

^{2}d

*y*=1.) Lack of orthogonality is particularly striking in the case of

*E*

_{y}profiles, where nearby modes exhibit a substantial degree of overlap. The apparent regularity of the inner products might signify a deep connection with the abnormal roots depicted in Fig. 2, although their relationship is not understood at present.

## 3. Excitation of Bloch modes by an incident plane wave

*z*<0. These include both homogeneous and inhomogeneous plane-waves, which are expressed by the general form:

*n*is the mode number (0, ±1, ±2, ±3, ……),

*θ*)+

*nλ*

_{o}/

*p*, and (

^{2}+(

*σ*

_{y},

*σ*

_{z}from similar symbols used for the slit array in Section 2.) The

*E*-field can be obtained in the same way as before, from Eq. (2). At normal incidence, where

*θ*=0°, symmetry allows only even modes, namely,

*z*=0,

*H*

_{x}and

*E*

_{y}in the free-space region

*z*<0–where the incident plane-wave (

*H*

_{xI},

*E*

_{yI},

*E*

_{zI}) must be added to the sum of homogeneous and evanescent reflected modes–are set equal to the corresponding fields resulting from a superposition of the excited Bloch modes in the

*z*>0 region of the slit array, that is,

*n*

^{th}reflected and transmitted modes, respectively. Equation (10) may be further simplified by moving all the terms to the left-hand-side, that is,

*C*

_{n}}={-

*p*] interval along the

*y*–axis. The total number of modes, 2

*N*, included in these calculations is equally divided between reflected and transmitted modes.

### 3.1 The method of least squares

*C*

_{n}} by minimizing the integral of the squared modulus of the mismatch at the interface. The mismatch arises from the difference between the

*E*

_{y}-field profiles just above (

*z*=0

^{-}) and just below (

*z*=0

^{+}) the interface, as well as the difference between the

*H*

_{x}-field profiles in these two (adjacent) regions. We fix the number of modes in each region at some desired value

*N*. (The mismatch should decline as more modes are incorporated.) The least-squares error may be written

*C*

_{n}=

*C*

_{n}

*′*+i

*C*

_{n}

*″*is the (initially unknown) coefficient of the

*n*

^{th}excited mode. Since in the MKSA system of units the

*E*-field magnitude is generally greater than that of the

*H*-field, we have used the weighting factor

*N*modes, there are 4

*N*real-valued unknowns {

*C*

_{n}

*′*,

*C*

_{n}

*″*}. Setting to zero the derivative of the

*Error*function with respect to each unknown, namely,

*∂Error*/

*∂C*

_{n}

*′*=0 and

*∂Error*/

*∂C*

_{n}

*″*=0, we find the following 2

*N*equations for the 2

*N*complex unknowns {

*C*

_{n}}:

*j*, ranging from 1 to 2

*N*, is the index of the unknown coefficient,

*C*

_{j}, with respect to which the partial derivatives were evaluated. Equation (13), a set of 2

*N*linear equations in 2

*N*unknowns, may be written in matrix form, then solved by inverting the coefficients matrix.

*R*

_{n}and

*T*

_{n}, as well as the total reflected and transmitted optical power,

*R*and

*T*, may be computed using the corresponding component

*S*

_{z}of the Poynting vector

*along the*

**S***z*-axis, that is,

*R*

_{n}=0 for the evanescent modes, while, for propagating modes,

*R*

_{n}is the

*n*

^{th}order diffraction efficiency. Orthogonality of the modes in this region allows one to obtain the total reflectance

*R*by summing the various

*R*

_{n}. In transmission, however, the various Bloch modes are not mutually orthogonal, and the total transmissivity must be computed directly by evaluating the

*z*-component of the Poynting vector, as given by Eq. 14(c). While the transmitted guided mode loses its energy slowly to the vertical slit walls, the non-propagating modes of the slit array deposit their entire energy in a thin layer near the top surface of the metallic structure.

### 3.2. Convergence

*p*=0.9µm,

*w*=0.1µm,

*λ*

_{o}=1.0µm, normal incidence), Fig. 5 shows the

*E*

_{y}and

*H*

_{x}field profiles (on both sides of the

*z*=0 interface) with

*N*=80 modes included in the computation. The match is excellent, both in the slit region and over the metal surface. The difference between the incident power and the computed

*R*+

*T*was found to be less than 0.3%. Although the high-frequency oscillations of

*E*

_{y}within the slit may, in part, be a manifestation of the Gibbs phenomenon [21], the sharp peaks at the slit edges are real and represent the significant accumulation of electrical charge on these sharp corners.

*C*

_{n}| as functions of

*N*(the number of modes on each side of the

*z*=0 interface used to minimize the mismatch). Figure 6(a) shows the amplitude variations of the first five Bloch modes of the slit-array versus

*N*, indicating that the mode coefficients stabilize after incorporating ~10–15 modes on each side of the interface. Figure 6(b) shows the magnitudes of the first 25 Bloch modes of the slit-array when a total of

*N*=80 modes are used to match the boundary conditions. It is readily observed that the mode amplitudes are strong for the low-order modes (i.e., dominant modes), while higherorder modes are fairly weak. The solution converges rather slowly after the first few tens of modes, which is reminiscent of the RCWA behavior [17

17. M. G. Moharam, E. B. Grann, D. A. Pommet, and T. K. Gaylord, “Formulation for stable and efficient implementation of the rigorous coupled-wave analysis of binary gratings,” J. Opt. Soc. Am. A **12**, 1068–76 (1995). [CrossRef]

18. Ph. Lalanne and G. M. Morris, “Highly improved convergence of the coupled-wave method for TM polarization,” J. Opt. Soc. Am. A **13**, 779–84 (1996). [CrossRef]

## 4. Transmittance of the slit array

*w*=0.1µm at

*λ*

_{o}=1.0 µm (normal incidence) versus the period

*p*(ranging from 0.4µm to 3.2µm) are shown in Fig. 7. The blue curve represents the transmittance

*T*

_{1}of the guided mode only, whereas the red curve corresponds to the total transmission efficiency

*T*, which includes contributions by all the Bloch modes of the slit array. The non-guided modes of the array are, of course, absorbed within the first few nanometers upon entering the top facet of the metallic medium.

*p*=0.4µm, both

*T*and

*T*

_{1}gradually decline to almost zero at

*p*=0.99µm, which is just below the Wood anomaly at

*p*=1.0µm. The transmission then climbs rapidly to a local maximum at

*p*=1.0µm, and falls immediately afterward. Subsequently,

*T*(

*p*) and

*T*

_{1}(

*p*) repeat this behavior (decline followed by a sharp rise and a local peak) for

*p*=1.0-2.0µm and

*p*=2.0-3.0µm. The minima of

*T*occur at those values of

*p*that are generally associated with the excitation of SPP, i.e., when

*p*becomes an integer-multiple of the SPP wavelength

*λ*

_{spp}=

*λ*

_{o}

*Real*

*λ*

_{spp}≈0.99µm.

*T*

_{1}) is extremely small, although not exactly zero, at

*p*=

*mλ*

_{spp}(

*m*an integer). For

*m*=1, 2, 3, we examined the neighborhood of each such point carefully, and found that the minimum of

*T*

_{1}in each case occurs at

*p*=

*mλ*

_{spp}, with

*T*

_{1}rising rapidly when the period

*p*deviates, ever so slightly, from

*mλ*

_{spp}.

### 4.1. Modal analysis of the anomalies

*p*=

*λ*

_{o}and

*p*=2

*λ*

_{o}. For this discussion, the mode-numbering scheme is as follows: In transmission, modes 1 through 5 are the top five modes depicted in Fig. 3 (albeit for the relevant value of the period

*p*); these are the five least lossy modes, starting with the guided mode as #1. In reflection, the first mode is the plane-wave directly reflected from the entrance facet of the slit array; this corresponds to the mode described by Eq. (8) with

*n*=0. The second mode is the sum of the two plane-waves in Eq. (8) having

*n*=±1; the left-right symmetry at normal incidence allows us to combine these two plane-waves into a single mode, as in Eq. (9). The third reflected mode is the sum of the two plane-waves having

*n*=±2, and so on.

*p*=

*λ*

_{o}. The case of

*p*=0.8µm represents the normal behavior of the slit array, where a good mix of modes is excited on both sides of the interface. As

*p*approaches

*λ*

_{spp}=0.99µm, however, most of the modes diminish, leaving the normally reflected planewave (#1) plus a single evanescent mode (#2) in the region above the slit array, and essentially nothing but two non-propagating Bloch modes (#2 and #3) within the slit array. The single evanescent mode above the array (#2 in reflection) is, of course, characteristic of an excited SPP at a flat metal surface (with no slits), and the two Bloch modes of the slit array (transmitted modes #2 and #3) in combination resemble the requisite SPP excitation beneath a flat metal surface (again with no slits).

*H*

_{x}(

*y*,

*z*=0

^{+})| and |

*E*

_{y}(

*y*,

*z*=0

^{+})| over a full period at the

*z*=0 interface (

*p*=0.99µm). The green curve contains only the first 3 transmitted modes, while the blue curve is a superposition of all 80 modes. The net effect is that, at

*p*=0.99µm, the few excited modes of the structure conspire to “hide” the presence of the slits, making the slit array behave as if the slits did not exist. With the slits all but wiped out from the picture, the handful of excited modes (plus the incident beam, of course) manage to satisfy the boundary conditions at the

*z*=0 interface. The guided mode (#1 in transmission) is excited only weakly (

*T*≈3%), the incident beam is nearly fully reflected (

*R*≈97%).

12. R. W. Wood, “On a remarkable case of uneven distribution of light in a diffraction grating spectrum,” Proc. Phys. Soc. London **18**, 269–275 (1902). [CrossRef]

22. Lord Rayleigh, “On the dynamic theory of gratings”, Proc. R. Soc. A **79**, 399–416 (1907). [CrossRef]

13. R. W. Wood, “Anomalous diffraction gratings,” Phys. Rev. **48**, 928–937 (1935). [CrossRef]

*p*=

*λ*

_{o}is caused by a redistribution of energy due to the passing-off of the evanescent mode at grazing angle along the interface [23]. Table I shows that at

*p*=

*λ*

_{o}many modes are once again necessary to yield a good match at the

*z*=0 interface, and the graphs in the second row of Fig. 8 confirm the inadequacy of the first few modes to accomplish this task all by themselves. In going from

*p*=

*λ*

_{spp}to

*p*=

*λ*

_{o}, the guided mode (#1 in transmission) suddenly rises to significance, the specularly reflected mode (#1 in reflection) weakens substantially, and the ±1

^{st}order diffracted beams (#2 in reflection), while remaining strong, switch from evanescence to propagating status. This upsets the balance of the

*E*-and

*H*-fields at the interface, thus requiring the excitation of many higher-order modes to restore the balance.

*p*=2

*λ*

_{o}. In general, as in the case of

*p*<

*λ*

_{spp}, many modes are needed to achieve a good match at the

*z*=0 interface. However, As

*p*approaches 2

*λ*

_{spp}=1.98µm from below, the guided mode (#1 in transmission) is extinguished, the specularly reflected mode (#1 in reflection) gains prominence, the ±1

^{st}order diffracted beams (#2 in reflection) nearly vanish, and the evanescent mode with the resonant SPP frequency (#3 in reflection) rises high above all other evanescent modes in the free-space region

*z*<0. The first few Bloch modes of the slit array (#2, 3, 4 in transmission) are then needed to match the boundary conditions at the

*z*=0 interface, as can be readily observed in the |

*H*

_{x}(

*y*,

*z*=0

^{+})| and |

*E*

_{y}(

*y*,

*z*=0

^{+})| plots of Fig. 8, third row. Once again, as in the case of

*p*=

*λ*

_{spp}, it is seen that at

*p*=2

*λ*

_{spp}the SPP excitation (most readily recognized in the dominance of the reflected mode #3), masks the presence of the slits, thus preventing the excitation of the guided mode (#1 in transmission), and returning nearly all the incident optical power via specular reflection (

*R*≈97%).

*p*increases beyond 2

*λ*

_{spp}, the arrival of the second Wood anomaly at

*p*=2

*λ*

_{o}=2.0µm once again wreaks havoc on the balance of the fields at the

*z*=0 boundary, throwing a large fraction of the incident optical power into the guided mode (#1 in transmission) and the ±1

^{st}order diffracted modes (#2 in reflection). The first few modes are no longer sufficient to satisfy the boundary conditions at the

*z*=0 interface, thus necessitating the excitation of higher order modes, as shown in Fig. 8, bottom row.

## 5. Oblique incidence

*p*=1.2µm,

*w*=0.1µm) in a semi-infinite silver host in the wavelength range

*λ*

_{o}=0.4-0.7µm at oblique incidence. Silver’s permittivity in the considered range is interpolated linearly from known discrete values [24], and the incidence angles are fixed at

*θ*=0°, 2°, and 4°. Computed plots of transmittance

*T*, specular reflectance

*R*

_{o}, and the ±1

^{st}order reflected diffraction efficiencies

*R*

_{±1}are plotted versus

*λ*

_{o}in Fig. 9. A remarkable feature of these curves is that the transmittance

*T*and the ±1

^{st}order diffraction efficiencies

*R*

_{±1}show very similar behavioral patterns. The anomalies are sharp and strong in the vicinity of

*λ*

_{o}=

*p*/2=0.6µm, and less so around

*λ*

_{o}=

*p*/3=0.4µm, indicating that a larger distance (in units of

*λ*

_{o}) between adjacent slits tends to reduce the effects of interference among the slits. We concentrate our attention where the anomalies are strong, i.e., at

*λ*

_{o}~0.6µm.

*h*

_{1m}|>|

*h*

_{2m}|, the other with |

*h*

_{1m}|<|

*h*

_{2m}|; see Eq. (1). Denoting the incidence wavelengths corresponding to the Wood and SPP anomalies by

*λ*

_{Wood}and

*λ*

_{spp}, their values may be calculated from:

*m*=±2 and

*n*

_{spp}=

*Real*

*θ*=2°,

*λ*

_{Wood}=0.6±0.02µm,

*λ*

_{spp}=0.62±0.02µm, and in the case of

*θ*=4°,

*λ*

_{Wood}=0.6±0.04µm,

*λ*

_{spp}=0.62±0.04µm, in agreement with the numerical computation results depicted in Fig. 9. We will present a more detailed discussion of oblique incidence in a forthcoming paper.

## 6. Concluding remarks

7. Q. Cao and Ph. Lalanne, “Negative role of surface plasmons in the transmission of metallic gratings with very narrow slits,” Phys. Rev. Lett. **88**, 057403(4) (2002). [CrossRef]

11. Y. Xie, A.R. Zakharian, J. V. Moloney, and M. Mansuripur, “Transmission of light through slit apertures in metallic films,” Opt. Express **12**, 6106 (2004). [CrossRef] [PubMed]

15. Ph. Lalanne, C. Sauvan, J. P. Hugonin, J. C. Rodier, and P. Chavel, “Perturbative approach for surface plasmon effects on flat interfaces periodically corrugated by subwavelength apertures,” Phys. Rev. B **68**, 125404 (2003). [CrossRef]

*τ*, an altogether different kind of resonance can take place in the cavity formed within the slit waveguide. The guided mode bounces back and forth between the entrance and exit facets of the slit (i.e., at the top and bottom facets of the metal film, separated by a distance

*τ*), building strength in exactly the same way as the light trapped between the mirrors of a Fabry-Perot resonator gains strength by multiple reflections. The low transmissivity of the semi-infinite slit array (i.e.,

*τ*=∞) is, therefore, a necessary condition for the strong build-up of trapped light in a slit array of finite thickness. As the period

*p*of a given array moves closer to

*λ*

_{spp}, the semi-infinite array becomes less transmissive, but, by the same token, the transmissivity of a finite thickness array (having the same values of

*p*and

*w*) can become extremely large if the Fabry-Perot resonance condition is satisfied, namely,

*φ*

_{r}+2

*πτ*/

*λ*=

*mπ*, where

*φ*

_{r}is the phase acquired by the guided mode upon internal reflection from the top (or bottom) of the slit,

*λ*is the guided mode’s propagation wavelength, and

*m*is an integer. This behavior, which has been extensively studied by Lalanne

*et al*. [15

15. Ph. Lalanne, C. Sauvan, J. P. Hugonin, J. C. Rodier, and P. Chavel, “Perturbative approach for surface plasmon effects on flat interfaces periodically corrugated by subwavelength apertures,” Phys. Rev. B **68**, 125404 (2003). [CrossRef]

16. Ph. Lalanne, J. P. Hugonin, and J. C. Rodier, “Theory of Surface Plasmon Generation at Nanoslit Apertures,” Phys. Rev. Lett. **95**, 263902 (2005). [CrossRef]

8. Y. Xie, A. R. Zakharian, J. V. Moloney, and M. Mansuripur, “Transmission of light through a periodic array of slits in a thick metallic film,” Opt. Express **13**, 4485 (2005) [CrossRef] [PubMed]

## Acknowledgments

## References and links

1. | T. W. Ebbesen, H. J. Lezec, H. F. Ghaemi, T. Thio, and P. A. Wolff, “Extraordinary optical transmission through subwavelength hole arrays,” Nature |

2. | M. M. J. Treacy, “Dynamical diffraction explanation of the anomalous transmission of light through metallic gratings,” Phys. Rev. B |

3. | H. J. Lezec and T. Thio, “Diffracted evanescent wave model for enhanced and suppressed optical transmission through subwavelength hole arrays,” Opt. Express |

4. | G. Gay, O. Alloschery, B. Viaris de Lesegno, C. O’Dwyer, J. Weiner, and H. J. Lezec, “The optical response of nanostructured surfaces and the composite diffracted evanescent wave model,” Nature Phys. |

5. | S. H. Chang, S. Gray, and G. Schatz, “Surface plasmon generation and light transmission by isolated nanoholes and arrays of nanoholes in thin metal films,” Opt. Express |

6. | 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. Communications |

7. | Q. Cao and Ph. Lalanne, “Negative role of surface plasmons in the transmission of metallic gratings with very narrow slits,” Phys. Rev. Lett. |

8. | Y. Xie, A. R. Zakharian, J. V. Moloney, and M. Mansuripur, “Transmission of light through a periodic array of slits in a thick metallic film,” Opt. Express |

9. | J. A. Porto, F. J. García-Vidal, and J. B. Pendry, “Transmission resonance on metallic gratings with very narrow slits,” Phys. Rev. Lett. |

10. | J. Bravo-Abad, L. Martín-Moreno, and F. J. García-Vidal, “Transmission properties of a single metallic slit: from the subwavelength regime to the geometrical-optics limit,” Phys. Rev. E |

11. | Y. Xie, A.R. Zakharian, J. V. Moloney, and M. Mansuripur, “Transmission of light through slit apertures in metallic films,” Opt. Express |

12. | R. W. Wood, “On a remarkable case of uneven distribution of light in a diffraction grating spectrum,” Proc. Phys. Soc. London |

13. | R. W. Wood, “Anomalous diffraction gratings,” Phys. Rev. |

14. | H. Raether, Surface Plasmons on smooth and rough surfaces and on gratings, (Springer-Verlag, Berlin, 1986). |

15. | Ph. Lalanne, C. Sauvan, J. P. Hugonin, J. C. Rodier, and P. Chavel, “Perturbative approach for surface plasmon effects on flat interfaces periodically corrugated by subwavelength apertures,” Phys. Rev. B |

16. | Ph. Lalanne, J. P. Hugonin, and J. C. Rodier, “Theory of Surface Plasmon Generation at Nanoslit Apertures,” Phys. Rev. Lett. |

17. | M. G. Moharam, E. B. Grann, D. A. Pommet, and T. K. Gaylord, “Formulation for stable and efficient implementation of the rigorous coupled-wave analysis of binary gratings,” J. Opt. Soc. Am. A |

18. | Ph. Lalanne and G. M. Morris, “Highly improved convergence of the coupled-wave method for TM polarization,” J. Opt. Soc. Am. A |

19. | J. D. Jackson, |

20. | A. W. Snyder and J. D. Love, |

21. | R. N. Bracewell, |

22. | Lord Rayleigh, “On the dynamic theory of gratings”, Proc. R. Soc. A |

23. | E. Noponen, “Electromagnetic Theory of Diffractive Optics,” dissertation, Dept. of Technical Physics, Helsinki University of Technology, Finland (1994). |

24. | P. Edward, |

**OCIS Codes**

(050.1220) Diffraction and gratings : Apertures

(050.1960) Diffraction and gratings : Diffraction theory

(240.6680) Optics at surfaces : Surface plasmons

(260.3910) Physical optics : Metal optics

(310.2790) Thin films : Guided waves

**ToC Category:**

Diffraction and Gratings

**History**

Original Manuscript: March 27, 2006

Revised Manuscript: June 22, 2006

Manuscript Accepted: June 22, 2006

Published: July 10, 2006

**Citation**

Y. Xie, A. R. Zakharian, J. V. Moloney, and M. Mansuripur, "Transmission of light through periodic arrays of sub-wavelength slits in metallic hosts," Opt. Express **14**, 6400-6413 (2006)

http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-14-14-6400

Sort: Year | Journal | Reset

### References

- T. W. Ebbesen, H. J. Lezec, H. F. Ghaemi, T. Thio, P. A. Wolff, "Extraordinary optical transmission through subwavelength hole arrays," Nature 39, 667-669 (1998). [CrossRef]
- M. M. J. Treacy, "Dynamical diffraction explanation of the anomalous transmission of light through metallic gratings," Phys. Rev. B 66,195105-11 (2002). [CrossRef]
- H. J. Lezec and T. Thio, "Diffracted evanescent wave model for enhanced and suppressed optical transmission through subwavelength hole arrays," Opt. Express 12, 3629-3651 (2004). [CrossRef] [PubMed]
- G. Gay, O. Alloschery, B. Viaris de Lesegno, C. O'Dwyer, J. Weiner, H. J. Lezec, "The optical response of nanostructured surfaces and the composite diffracted evanescent wave model," Nature Phys. 264, 262 - 67 (2006). [CrossRef]
- S. H. Chang, S. Gray, and G. Schatz, "Surface plasmon generation and light transmission by isolated nanoholes and arrays of nanoholes in thin metal films," Opt. Express 13, 3150-3165 (2005). [CrossRef] [PubMed]
- 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. Communications 225, 331 (2003). [CrossRef]
- Q. Cao and Ph. Lalanne, "Negative role of surface plasmons in the transmission of metallic gratings with very narrow slits," Phys. Rev. Lett. 88, 057403(4) (2002). [CrossRef]
- Y. Xie, A. R. Zakharian, J. V. Moloney, and M. Mansuripur, "Transmission of light through a periodic array of slits in a thick metallic film," Opt. Express 13, 4485 (2005) [CrossRef] [PubMed]
- J. A. Porto, F. J. García-Vidal, J. B. Pendry, "Transmission resonance on metallic gratings with very narrow slits," Phys. Rev. Lett. 83, 02845(4) (1999). [CrossRef]
- J. Bravo-Abad, L. Martín-Moreno, F. J. García-Vidal, "Transmission properties of a single metallic slit: from the subwavelength regime to the geometrical-optics limit," Phys. Rev. E 69, 26601(6) (2004). [CrossRef]
- Y. Xie, A.R. Zakharian, J. V. Moloney, M. Mansuripur, "Transmission of light through slit apertures in metallic films," Opt. Express 12, 6106 (2004). [CrossRef] [PubMed]
- R. W. Wood, "On a remarkable case of uneven distribution of light in a diffraction grating spectrum," Proc. Phys. Soc. London 18, 269-275 (1902). [CrossRef]
- R. W. Wood, "Anomalous diffraction gratings," Phys. Rev. 48, 928-937 (1935). [CrossRef]
- H. Raether, Surface Plasmons on smooth and rough surfaces and on gratings, (Springer-Verlag, Berlin, 1986).
- Ph. Lalanne, C. Sauvan, J. P. Hugonin, J. C. Rodier, and P. Chavel, "Perturbative approach for surface plasmon effects on flat interfaces periodically corrugated by subwavelength apertures," Phys. Rev. B 68, 125404 (2003). [CrossRef]
- Ph. Lalanne, J. P. Hugonin, and J. C. Rodier, "Theory of Surface Plasmon Generation at Nanoslit Apertures," Phys. Rev. Lett. 95, 263902 (2005). [CrossRef]
- M. G. Moharam, E. B. Grann, D. A. Pommet, and T. K. Gaylord, "Formulation for stable and efficient implementation of the rigorous coupled-wave analysis of binary gratings," J. Opt. Soc. Am. A 12, 1068-76 (1995). [CrossRef]
- Ph. Lalanne and G. M. Morris, "Highly improved convergence of the coupled-wave method for TM polarization," J. Opt. Soc. Am. A 13, 779-84 (1996). [CrossRef]
- J. D. Jackson, Classical Electrodynamics, Chapter 8, 3rd edition, Wiley, New York, 1999.
- A. W. Snyder and J. D. Love, Optical Waveguide Theory, Chapman and Hall, London, 1983.
- R. N. Bracewell, The Fourier Transform and its Applications, McGraw-Hill, New York, 1978.
- Lord Rayleigh, "On the dynamic theory of gratings", Proc. R. Soc. A 79, 399-416 (1907). [CrossRef]
- E. Noponen, "Electromagnetic Theory of Diffractive Optics," dissertation, Dept. of Technical Physics, Helsinki University of Technology, Finland (1994).
- P. Edward, Handbook of optical constants of solids, 1st edition, Academic press, 1997.

## Cited By |
Alert me when this paper is cited |

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

« Previous Article | Next Article »

OSA is a member of CrossRef.