Surface Impedance Model for Extraordinary Transmission in 1d Metallic and Dielectric Screens

Extraordinary Optical Transmission of TM waves impinging at oblique incidence on metallic or high permittivity dielectric screens with a periodic distribution of 1D slits or any other kind of 1D defects is analyzed. Generalized waveguide theory altogether with the surface impedance concept are used for modeling such phenomena. A numerical analysis based on the mode matching technique proves to be an efficient tool for the characterization of these structures for any angle of incidence and slit or defect apertures. Extraordinary optical transmission through sub-wavelength hole arrays, nanoplasmonic biosensor for direct detection of live viruses from biological media, " Nano Lett. 10, 4962–4969 (2010). Surface plasmons enhance optical transmission through subwavelength holes, " Phys. " Enhanced millimeter-wave transmission through subwavelength hole arrays, " Opt. Optical properties of tungsten thin films perforated with a bidimensional array of subwavelength holes, " Phys. Analytical theory of extraordinary transmission through metallic diffraction screens perforated by small holes, " Opt. Analytical theory of extraordinary optical transmission through realistic metallic screens, " Opt. Extraordinary transmission through dielectric screens with 1D sub-wavelength metallic inclusions, " Opt. Transmission resonances on metallic gratings with very narrow slits, " Phys. Transmission and focusing of light in one-dimensional periodically nanostructured metals, " Phys. Extraordinary transmission through arrays of slits: a circuit theory model, " IEEE Trans. of the dielectric constant and loss tangent of high dielectric-constant materials at terahertz frequencies, " IEEE Trans.


Introduction
The phenomenon of Extraordinary Optical Transmission (EOT) through opaque screens with a periodic array of subwavelength holes, first reported in [1], is a current topic of intensive research due to its potential applications in photonic circuits [2], optical sensing [3] or fabrication of left-handed metamaterials.As nanofabrication techniques improve, efficient numerical models are desirable in order to provide fast characterization of the designs as well as physical insight into the phenomena.The physical background of the phenomena is an issue of controversy.Channeling of surface plasmons (SP) [4] through the holes was first proposed, but could not explain Extraordinary Transmission (ET) in metals at microwave frequencies [5] or in dielectric screens [6] where real part of permittivity is positive and SP do not exist.Subsequently, the surface plasmon concept was rescued to explain ET after considering that plasmon-like surface waves, or spoof plasmons, can be supported by structured metallic or dielectric surfaces [7].Recently, some of the authors presented an analytical theory of EOT through perfect conducting screens based on waveguide analysis [8].In this paper the role of spoof plasmons in extraordinary transmission was also discussed (see e.g.Table 1 in [8]).This theory was later extended to screens made of realistic conductors such as metals at optical frequencies [9] and 1D dielectric configurations [10].These generalizations were based on the surface impedance concept, which is widely used in classical electromagnetism for the analysis of skin effect in imperfect conductors [11].The model in [9,10] has some limitations as it is only valid in the case of normal incidence and for holes smaller than a quarter wavelength.In this letter, we present a unified analysis of extraordinary transmission at oblique incidence through periodic metallic or dielectric screens, loaded with transparent or opaque defects, including thin (sub-skin-depth) screens and wide (∼ λ /2) slits or defects.

Theory
We will consider a metallic or high permittivity dielectric screen with a periodic distribution of slits.These configurations have been analyzed by many authors [12][13][14][15] in order to simplify numerical problems and subsequently have taken on a live of its own.In this section such structure will be analyzed using an approximate mode-matching analysis which is simplified making use of the surface impedance concept.
The unit cell of the structure is depicted in Fig. 1.For an incident TM wave with the magnetic field parallel to the slits and the wave vector in the yz plane, we have a bidimensional problem with no x dependence.The impinging wave imposes phase shift boundary conditions at y = ±a/2.The transverse components of the fields in regions 1, 2 and 3 can be expanded in terms of Bloch and waveguide modes and evaluated in the input and output surfaces of the screen b

E
(2) where T and R are the transmission and reflection coefficients; T n , R n , S ± 0 and S ± m are the amplitudes of the different modes excited in the screen; k y,0 is the transversal component of the wave vector of the incident wave and are the propagation constants of the modes inside the slits which may be filled with a dielectric of relative permittivity ε r .Expressions for magnetic fields x are obtained by just multiplying each mode in Eqs. ( 1)-( 4) by its corresponding wave admittance.Expansions in Eqs. ( 2) and (3) correspond to a parallel plate waveguide with PEC plates at y = ±b/2, so losses or effects of field penetration inside the screen in the inner walls of the slits are neglected.This approximation is plausible only in the case of ε r << ε s , which is also the meaningful case for EOT.A mode-matching analysis is employed to calculate all the coefficients in Eqs. ( 1)-( 4).In this mode-matching analysis, specific relations between the electromagnetic fields at both sides of the screen (z = ±t/2) will be considered.In the area of the slits, continuity of transverse electromagnetic fields is imposed whereas in the area of the screen, each mode of the transverse electromagnetic fields in regions 1 and 3 evaluated in z = ±t/2 are related through a surface impedance matrix [9] x,n + H where and Z (2) are the surface impedances corresponding to each mode, are the propagation constants and wave admittances of the modes inside the metallic or high permittivity screen and ε s is the relative permittivity of the screen.The effect of the surface impedance in Eq. ( 8) is to provide an additional coupling between both sides of the screen, which may be dominant for thin sub skin-depth metallic screens, or in dielectric screens (see discussion in Sec. 3, at the end).It also accounts for losses in the screen.Matrix equation Eq. ( 7) is exact in the case of an infinite and homogenous plain slab; however, despite the presence of the slits, it proves to be a powerful approximation and allows us to considerably simplify the mode matching algorithm that strictly should also account for modal expansion inside the screen.When slits are very thin (t ≈ a/100) the model is less accurate as edge effects become more relevant.In any case, 1D configurations do not yield relevant transmission peaks when the screens are so thin.Now, multiplying Eq. ( 6) by and Eq. ( 7) by exp −i k y,0 + 2nπ a y , n = −N, ..., 0..., N; (11) and integrating in the corresponding ranges of y, we obtain a set of overdetermined equations which can be numerically solved using a least square algorithm.

Results
In Figs.2-6 the transmission coefficients obtained with the reported model are compared with electromagnetic simulations using CST Microwave Studio.In order to obtain convergent results, the resolution of the higher modes inside the slits (M) and in the input and output regions (N) must be similar.Convergent results were obtained employing N = 60 modes in Eq. ( 6) and M ≈ N(b/a) modes in Eq. ( 7) with a CPU time of ∼ 0.5s vs CPU time of ∼ 10s with the electromagnetic solver.In Figs. 2 and 3 the metallic screen is modeled by a finite conductivity σ = 59.6 × 10 6 S/m (corresponding to copper) at frequencies around 1 THz: and in Fig. 4 by a Drude-Lorentz permittivity for silver [16] at frequencies around 300 THz Tr ansmittance Tr ansmittance with plasma frequency ω p = 2π ×2175 THz and frequency of collision 35 THz, where a correction due to finite thickness of the screen has been taken into account [9].In the abscissas, frequencies are normalized to the Wood's anomalies frequencies for each angle of incidence for n < 0 (14) corresponding to the divergence of the TM 2n modes excited in the generalized waveguide.In the case of metallic screens, the numerical model works quite well for any angle of incidence and different Wood's anomalies resonances (n = −1, n = 1 and n = −2).At a given angle of incidence (θ = 19.47 o ) the Wood's anomaly frequency corresponding to the divergence of the scattered TM 2 and TM −4 modes is the same.For higher angles of incidence the peaks corresponding to the divergence of TM 2 mode vanish (Fig. 2(b)) and those corresponding to the divergence of TM −4 grow.In Fig. 2 the numerical model presented here is compared with the model reported in [9].It can be seen how this previous numerical model, in which variations of orthogonal functions in the region of the small slits are neglected, fails to characterize screens with wide slits.The proposed numerical model also applies to screens made of high permittivity dielectric (Figs. 5  and 6).However, in this case there are some limitations.When the Wood's anomaly frequency is close to a Fabry-Perot resonance and background transmission is close to unity, the numerical model does not provide physically meaningful results.In any case, when background transmission is high, EOT peaks lack practical interest.In Fig. 5 the slits are empty and in Fig. 6 the slits are filled with a PEC (Z n = 0), in which case we find that transmission peaks are higher and background transmission lower than in the case of empty slits.This effect was previously reported in [10] for normal incidence.Here is confirmed for all angles of incidence.In computations not shown here, we have checked the validity of the model for lower values of permittivity (ε s ≈ 10).However, for such low values of the permittivity, background transmission is high and EOT peaks lack interest.
In order to gain physical insight into the phenomenon, the longitudinal real part of the Poynting vector (ℜ(S z )) at the frequencies of some of the peaks, computed with the electromagnetic simulator are plotted in Fig. 7 and 8.In Fig. 7 the screens are modeled as lossy copper.The  graphics are complex, with vortexes of power bouncing forward and backward along the structure; a fact which express the important role played by evanescent modes near the screen.A high concentration of power appears at the edges of the slits, where strong currents are excited.Figure 7(a) corresponds to the divergence of the TM −2 mode and Fig. 7(b) to the divergence of the TM −4 mode excited in the screen for an angle of incidence θ = 20 o .In both cases, the power flow pattern corresponding to each resonance can be clearly appreciated.In Fig. 8 the screens are modeled as zirconium-tin-titanate (ε = 92.7(1+ 0.005i)).In Fig. 8(a) the slits are empty and in Fig. 8(b) the slits are filled with PEC.In both cases the patterns correspond to the divergence of the TM −2 mode for an angle of incidence θ = 20 o ; and show that most power is transmitted through the dielectric screens, thus suggesting that EOT in these structures is obtained as a consequence of the periodic distribution of defects, and no necessarily associated to the presence of transparent slits.Even more, in Fig. 8(a) the empty slits seem to act as a drain of power, which is first transmitted through the dielectric from left to right and then partially put back in the left space through the slits.

Conclusion
A new model for oblique incidence EOT in metallic or dielectric screens with 1D periodic arrays of slits or other kind of defects has been provided.The model is based on a combination of the generalized waveguide modal analysis with the surface impedance concept.A numerical technique has been developed based on such model, which allows for the rapid characterization of EOT in metallic or high permittivity dielectric screens loaded with 1D transparent or opaque defects.This technique provides a unified treatment of all these examples, revealing that the same physics is behind them and that periodicity, instead of the specific mechanism of power transmission, is the relevant physical fact.The analyzed structures could be useful as spatial filters for certain angles of incidence and frequencies.In the case of dielectric screens, metallic inclusions in the slits result in enhanced transmission peaks, suggesting the possibility of interesting phenomena, such as optically induced EOT in photoconductive semiconductor screens.

Fig. 1 .
Fig. 1.Front (a) and side (b) view of the screen with a periodic distribution of slits.Unit cell of the structure (c) with the incident TM wave.The impinging wave imposes periodic boundary conditions with a constant phase shift ∆φ = k y,0 a along y.

Fig. 2 .
Fig.2.Transmission through an array of slits in a lossy copper screen (σ = 59.6×10 6 S/m) at normal incidence and different sizes of the slits.Periodicity is a = 300µm and thickness of the screen is t = a/20.Continuous lines correspond to mode matching model, dashed lines to CST simulations and dot-dashed lines to our previous numerical model[9].

Fig. 3 .Fig. 4 .Fig. 5 .
Fig. 3. Transmission through an array of slits in a lossy copper screen (σ = 59.6×10 6 S/m) for different angles of incidence and Wood's anomalies resonances.Periodicity is a = 300µm, size of slits is b = a/4 and thickness of the screen is t = a/20.Continuous lines correspond to mode matching model and dashed lines to CST simulations.

Fig. 6 .
Fig.6.Transmission through an array of PEC inclusions in a zirconium-tin-titanate (ε = 92.7(1+ 0.005i)[17]) screen for different angles of incidence and Wood's anomalies resonances.Periodicity is a = 3mm, size of the inclusions is b = a/6 and thickness of the screen is t = a/25.Continuous lines correspond to mode matching model and dashed lines to CST simulations.

Fig. 7 .
Fig. 7. Propagating component of the real part of the Poynting vector at the frequencies of the EOT peaks for the lossy copper screen analyzed in Fig.2.In (a) EOT is associated to the divergence of the T M −2 mode and in (b) to the divergence of the T M −4 mode.Angle of incidence is θ = 20 o in both cases.

Fig. 8 .
Fig.8.Propagating component of the real part of the Poynting vector at the frequencies of the EOT peaks for the zirconium-tin-titanate screen analyzed in Figs.5 and 6.In (a) the slits are empty and in (b) filled with PEC.Angle of incidence is θ = 20 o and EOT is associated to the divergence of the T M −2 mode in both cases.