On the extraordinary optical transmission in parallel plate waveguides for non-TEM modes

On the extraordinary optical transmission in parallel plate waveguides for non-TEM modes MIGUEL CAMACHO,1,* RAFAEL R. BOIX,2FRANCISCO MEDINA,2 ALASTAIR P. HIBBINS,1AND J. ROY SAMBLES1 1Electromagnetic and Acoustic Materials Group, Department of Physics and Astronomy, University of Exeter, Stocker Road, Exeter EX4 4QL, UK 2Department of Electronics and Electromagnetism, College of Physics, University of Seville, Avda Reina Mercedes s/n, Seville, Spain *[email protected] Abstract: Extraordinary transmission has been recently measured in a parallel plate waveguide (PPWG) through a metal strip with a patterned 1-D periodic array of circular holes, the metal strip being embedded inside the PPWG. Wood’s anomaly and the extraordinary transmission peak (EOT) were detected for transverse electric (TE) mode excitation at frequencies higher than those found for TEM mode excitation. In this paper we provide an explanation for this frequency shift by decomposing the problem of a TE mode impinging on the 1-D array of holes into two problems of plane waves impinging obliquely on 2-D periodic arrays of holes. By then solving the integral equation for the electric field on the surface of the holes, the origin of the frequency shift is proved both mathematically and physically in terms of the symmetries present in the system. Published by The Optical Society under the terms of the Creative Commons Attribution 4.0 License. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI. OCIS codes: (230.7370) Waveguides; (050.1755) Computational electromagnetic methods, (260.2110) Electromagnetic optics. References and links 1. 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, 667–669 (1998). 2. H. A. Bethe, “Theory of diffraction by small holes,” Phys. Rev. 66, 163–182 (1944). 3. 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, 6779–6782 (1998). 4. M. Beruete, M. Sorolla, I. Campillo, J. S. Dolado, L. Martín-Moreno, J. Bravo-Abad, and F. J. García-Vidal, “Enhanced millimeter-wave transmission through subwavelength hole arrays,” Opt. Lett. 29, 2500 (2004). 5. J. B. Pendry, L. Martín-Moreno, and F. J. Garcia-Vidal, “Mimicking surface plasmons with structured surfaces.” Science 305, 847–848 (2004). 6. F. J. García De Abajo, and J. J. Sáenz, “Electromagnetic surface modes in structured perfect-conductor surfaces.” Phys. Rev. Lett. 95 233901 (2005), 7. F. J. Garcia-Vidal, L. Martin-Moreno, T. W. Ebbesen, and L. Kuipers, “Light passing through subwavelength apertures,” Rev. Mod. Phys. 82, 729–787 (2010). 8. C. Genet and T. W. Ebbesen, “Light in tiny holes,” Nature 445, 39–46 (2007). 9. F. J. Garcia De Abajo, “Colloquium: Light scattering by particle and hole arrays,” Rev. Mod. Phys. 79, 1267–1290 (2007). 10. J. Bravo-Abad, A. I. Fernández-Domínguez, F. J. García-Vidal, and L. Martín-Moreno, “Theory of extraordinary transmission of light through quasiperiodic arrays of subwavelength holes,” Phys. Rev. Lett. 99, 203905 (2007). 11. M. Camacho, R. R. Boix, and F. Medina, “Comparative study between resonant transmission and extraordinary transmission in truncated periodic arrays of slots,” in “2017 IEEE MTT-S International Conference on Numerical Electromagnetic and Multiphysics Modeling and Optimization for RF, Microwave, and Terahertz Applications (NEMO),” (IEEE, 2017), pp. 257–259. 12. Y. Pang, A. Hone, P. So, and R. Gordon, “Total optical transmission through a small hole in a metal waveguide screen,” Opt. Express 17, 4433–4441 (2009). 13. F. Medina, J. A. Ruiz-Cruz, F. Mesa, J. M. Rebollar, J. R. Montejo-Garai, and R. Marqués, “Experimental verification of extraordinary transmission without surface plasmons,” Appl. Phys. Lett. 95, 071102 (2009). Vol. 25, No. 20 | 2 Oct 2017 | OPTICS EXPRESS 24670 #304027 https://doi.org/10.1364/OE.25.024670 Journal © 2017 Received 2 Aug 2017; revised 22 Sep 2017; accepted 22 Sep 2017; published 27 Sep 2017 14. F. Medina, F. Mesa, and R. Marqués, “Extraordinary transmission through arrays of electrically small holes from a circuit theory perspective,” IEEE Trans. Microw. Theory Techn. 56, 3108–3120 (2008). 15. K. S. Reichel, P. Y. Lu, S. Backus, R. Mendis, and D. M. Mittleman, “Extraordinary optical transmission inside a waveguide: spatial mode dependence,” Opt. Express 24, 28221 (2016). 16. D. M. Pozar, Microwave Engineering (John Wiley & Sons Inc, 2005). 17. M. Camacho, R. R. Boix, and F. Medina, “Computationally efficient analysis of extraordinary optical transmission through infinite and truncated subwavelength hole arrays,” Phys. Rev. E 93, 063312 (2016). 18. M. Camacho, R. R. Boix, F. Medina, A. P. Hibbins, and J. R. Sambles, “Theoretical and experimental exploration of finite sample size effects on the propagation of surface waves supported by slot arrays,” Phys. Rev. B 95, 245425 (2017). 19. R. F. Harrington, Field Computation by Moment Methods (Wiley-IEEE, 1993). 1. Introduction With the experimental discovery of extraordinary optical transmission (EOT) in the late nineteen nineties [1], it was unveiled that the transmission through small holes could be accomplished by taking advantage of periodicity, largely overcoming the predictions of Bethe’s single aperture theory [2]. It was found that an enhanced transmission peak was always present at a frequency lower than that of the first onset of diffraction, even when the intrinsic resonance of the hole was well above the latter. In the optical regime, this was explained through the coupling of the free space radiation to bound modes commonly known as surface plasmon polaritons (SPPs) [3]. The finding of enhanced transmission at frequencies at which metals do not support SPPs, such as microwaves [4], made it clear that genuine SPPs were not essential to explain EOT. The EOT theory was later unified across a large range of frequencies thanks to the existence of surface modes supported by periodic arrays of holes at microwave frequencies, which were shown in some cases to behave in a similar dispersive manner to the optical SPPs [5,6]. A large body of work has been undertaken by both the optics and microwave engineering communities on this topic, collected into extensive review articles [7–9] and it still remains important due to the rich physics behind it [10,11]. Although most of the studies of the transmission through small apertures have been based on free-space wave propagation, the EOT phenomenon is also found in the transmission through metallic waveguide diaphragms near the cut-offof higher order waveguide modes, which play the role of the onset of diffracted modes in periodic structures [12,13]. Moreover, the use of transmission line theory to explain the EOT as an impedance matching problem is based on the transverse electromagnetic (TEM) mode supported by a parallel plate waveguide (PPWG) [14]. Only recently has the EOT effect been studied for transmission through a 1-D periodic array of circular holes embedded in a PPWG by Reichel et al. [15]. In particular, these researchers measured and computed the transmission coefficient when the array was excited both by the TEM mode and by the PPWG TE1mode. They found that the EOT peak appears at higher frequencies in the latter case than in the former case. In this paper we present an approach to understand the origin of the EOT peak frequency shift in PPWG under TE1mode excitation. This theory is based on the decomposition of the problem into two problems of EOT through a 2-D periodic array of symmetric holes under oblique plane wave excitation. These EOT problems are in turn accurately analyzed by solving an integral equation involving a periodic Green’s function. The steps followed in the solution of the integral equation provide a way to predict the appearance of Wood’s anomalies and the associated extraordinary transmission through both the 2-D array of holes and the equivalent 1-D array of holes in PPWG under TE1mode excitation. Our approach provides new physical insight into the origin of Wood’s anomalies and proposes a highly efficient analysis tool for the numerical study of structures supporting EOT phenomena. It only requires the solution of a small system of equations to obtain accurate results as opposite to other numerical methods such as coupled-mode method (also known as mode-matching). And this is because it adequately Vol. 25, No. 20 | 2 Oct 2017 | OPTICS EXPRESS 24671 accounts for the singular behavior of the field quantities at the edges of the metal screen holes. This singular behavior cannot be handled by the expansion modes used in the coupled-mode method presented in [15]. 2. Theoretical model a b ls w s x y z x y z Fig. 1. Perspective view of a parallel plate waveguide. The two plates are connected through a negligible thickness PEC strip perforated with an infinite 1-D periodic array of slots. A zoomed view of the hole array with the definition of the geometry parameters is also shown on the bottom-left of the figure. Let us consider the problem depicted in Fig. 1, where one pure TEz 1mode propagating in an infinite PPWG impinges on a negligible thickness conducting strip filling the space between the waveguide plates. The strip, located at z=0, is perforated with a 1-D periodic array of holes. The holes will be assumed to be rectangular slots for the sake of mathematical simplicity, although the physics governing the problem would remain identical if circular holes were to be chosen as in [15]. Both the parallel plates and the strip containing the holes will be assumed to be perfect electric conductors (PECs). In the following, a time-harmonic dependence of the physical quantities of the type ejωtwill be assumed and suppressed throughout. Under these conditions, the complex electric field of the TEz 1mode incident on the array of slots can be written as [16] EPPWG 1=EPPWG 0e−jk1,−−e−jk1,+ˆ x(1) where k1,−=−(π/b)ˆ y+k2 0−(π/b)2ˆ zand k1,+=+(π/b)ˆ y+k2 0−(π/b)2ˆ z(k0=ω√μ00= 2π/λ0,λ0being the free space wavelength). Equation (1) indicates that the TEz 1mode of the PPWG can be seen as the sum of two plane waves, the first propagating with wavenumber vector k1,−along the direction given by the angular spherical coordinates (φ− inc,θ− inc)=(3π/2,θ 0), and the second propagating with wavenumber vector k1,+along the direction given by the angular spherical coordinates (φ+ inc,θ + inc)=(π/2,θ 0), where θ0=arctan ⎛⎜⎜⎜⎜⎜⎜⎜⎜⎜⎝ π k2 0b2−π2⎞⎟⎟⎟⎟⎟⎟⎟⎟⎟⎠(2) Imagine each of the two aforementioned plane waves separately propagates in free space and is obliquely incident on a 2-D periodic array of rectangular slots perforated in a PEC screen at z=0 shown in Fig. 2, where the unit cell is exactly equal to that of the 1-D periodic structure of Fig. 1. If we superpose these two waves in the way shown in (1), the resulting wave will fulfill the electric wall (EW) boundary conditions imposed by the PPWG of Fig. 1 at y=0 and y=b. This means the problem of the transmission of the TEz 1mode through the 1-D array of holes can Vol. 25, No. 20 | 2 Oct 2017 | OPTICS EXPRESS 24672 be viewed as the superposition of two problems of transmission of plane waves propagating in free space which are obliquely incident on a 2-D periodic array of holes in a conducting screen, and the superposition principle can be invoked. b inc ϕinc ki Fig. 2. Perspective view of part of the two dimensional periodic array of holes perforated into a negligible thickness PEC screen. The array is illuminated by an obliquely incident plane wave along the direction given by the spherical angular coordinates (φinc,θinc). A zoomed view of the hole array with the definition of the geometry parameters is also shown. Now let us concentrate on the problem of a plane wave obliquely incident on the 2-D array of rectangular slots in a PEC screen that is shown in Fig. 2. Let us assume that the propagation direction of the plane wave is given by the angular spherical coordinates (φinc,θ inc) shown in Fig. 2. The complex electric field of this wave can be written as Ei=E0ej(kx0x+ky0y+kz0z)ˆ u=E0e−jki·rˆ u(3) where kx0=−k0sinθinc cosφinc,ky0=−k0sinθinc sinφinc,kz0=−k0cosθinc, and ˆ uis a unit vector such that ˆ u·ki=0. In accordance with the derivations shown in [17] (only valid for the case of normal incidence), in order to determine the scattered tangential electric field, Esc t(x,y,z=0), in the rectangular slot δ00 ≡{(a−w)/2<x<(a+w)/2; (b−l)/2<y<(b+l)/2}of the reference unit cell C00 ≡{0≤x≤a;0≤y≤b}of the 2-D periodic array of slots shown in Fig. 2, we need to solve the integral equation Jas +a 0b 0Gper M(x−x,y−y)·Esc t(x,y,z=0) dxdy=0(x,y)∈δ00 (4) where Jas is the surface current density that would be induced by the obliquely incident wave in the PEC screen in the absence of rectangular slots, and Gper M(x−x,y−y) is the periodic dyadic Green’s function given by [18] Gper M(x,y)= +∞  m=−∞ +∞  n=−∞ GM(x−ma,y−nb)e j(kx0ma+ky0nb).(5) In Equation (5), GM(x,y) stands for the dyadic Green’s function defined by GM(x,y)=⎛⎜⎜⎜⎜⎜⎜⎜⎝k2 0+∂2 ∂y2g(x,y)−∂2g(x,y) ∂x∂y −∂2g(x,y) ∂x∂yk2 0+∂2 ∂x2g(x,y)⎞⎟⎟⎟⎟⎟⎟⎟⎠(6) Vol. 25, No. 20 | 2 Oct 2017 | OPTICS EXPRESS 24673 where g(x,y)=−je−jk0√x2+y2 πk0Z0x2+y2(7) and where k0and Z0represent the free space wavenumber and impedance (Z0=μ0/0) respectively. In order to solve the integral equation of (4), the unknown quantity Esc t(x,y,z=0) is expanded in terms of basis functions as shown below Esc t(x,y,z=0) ≈ Nb  j=1 fjbj(x,y)(x,y)∈δ00 (8) where bj(x,y)(j=1,...,Nb) are chosen to be Chebyshev polynomials multiplied by the edge behavior expected for each of the two possible polarizations (see [17] for details). When the approximation of (8) is substituted in (4) and the Galerkin’s version of the Method of Moments (MoM) [19] is applied, the following system of linear equations is obtained for the unknown coefficients fj Nb  j=1 Γ ij fj=Ci(i=1,...,Nb) (9) The coefficients of the MoM matrix of (9) can be obtained in the Fourier transform domain as [17] Γ ij =ab ∞  m,n=−∞  b∗ i(kxm,kyn)· GM(kx=kxm,ky=kyn)· bj(kxm,kyn)(10) where bj(kxm,kyn)(j=1,...,Nb) are the 2-D discrete Fourier transforms of bj(x,y), kxm = 2πm/a+kx0and kyn =2πn/b+ky0. Also,  GM(kx=kxm,ky=kyn) is the 2-D continuous Fourier transform of GM(x,y) sampled at kx=kxm and ky=kyn, which is given by  GM(kx=kxm,ky=kyn)=−2 k0Z0k2 0−k2 xm −k2 yn ×k2 0−k2 yn kxmkyn kxmkyn k2 0−k2 xm(11) The constant terms Ci(i=1,...,Nb) of the system of equations (9) can be obtained as Ci=−a 0b 0bi(x,y)e−j(kx0x+ky0y)dxdy∗ ·Jas =−ab b∗ i(kx0,ky0)·Jas.(12) Note that the spectral dyadic Green’s function  GM(kx=kxm,ky=kyn) has poles at frequencies for which k2 0−k2 xm −k2 yn =0⇒k2 xm +k2 yn =k2 0(13) When condition (13) is fulfilled, the elements of  GM(kx=kxm,ky=kyn) tend to infinity, the MoM matrix entries Γ ij usually tend to infinity and the unknown coefficients fjtend to zero Vol. 25, No. 20 | 2 Oct 2017 | OPTICS EXPRESS 24674 to keep Cibounded in (9). As a consequence of this, Esc t(x,y,z=0) becomes zero in (8), which means the tangential electric field in the holes of the 2-D array vanishes, i.e., the PEC screen with holes behaves as a solid PEC and reflects all the incoming power, which is the origin of a Wood’s anomaly. Therefore, Equation (13) provides the condition for Wood’s anomalies and marks the onset of the m−n-th diffracted mode (also called m−n-th grating lobe). At frequencies just below this divergence, the fields on the slot are largely enhanced, providing a mechanism for the appearance of EOT peaks. This has also been explained in terms of circuit theory reported in [14], showing that an EOT peak always appears at a frequency slightly lower than a Wood’s anomaly owing to impedance matching provided by a resonance created by a frequency dependent capacitance that becomes singular at the Wood’s anomaly. So, a Wood anomaly always presents an associated EOT peak at slightly lower frequencies. 3. Results and discussion Transmission Coefficient T00 Transmission Coefficient T00 b) TEz a) TMz Einc Hinc ki bb lsls z yx Einc Hinc ki aa wsws z y x Fig. 3. Transmission spectra of a plane wave impinging on a 2-D array of rectangular slots in a PEC screen. Two polarizations TMz(a) and TEz(b) are considered, and different angles of incidence (θinc). In both cases, the direction of the incident electric field is contained in the x−zplane so as to excite the fields in the slots. The dimensions of the unit cell were chosen to be ls/a=0.4, ws/a=0.05 and a=b. For a PEC screen with a 2-D array of holes under normal incidence conditions (θinc =0), the first Wood’s anomaly is expected to appear at a frequency for which λ0=max(a,b). This corresponds to the Wood’s anomaly experimentally detected in [15] for a PPWG under TEM excitation since this problem is equivalent to the EOT problem of a 2-D periodic array of holes under normal incidence [17]. However, for the 2-D array of holes of Fig. 2 under oblique incidence, Wood’s anomalies and the associated EOT peaks in principle tend to appear at frequencies above and below the frequency for which λ0=max(a,b). For example, let us consider the transmission curves Vol. 25, No. 20 | 2 Oct 2017 | OPTICS EXPRESS 24675 shown in Fig. 3(a) for a 2-D periodic array of rectangular slots illuminated by obliquely incident plane waves propagating along different directions given by (φinc =0,θ inc) (all the curves in Fig. 3 have been obtained by means of MoM in accordance with Eqns. (4) to (12)). For these particular incidence directions, we would expect Wood’s anomalies for the m=±1,n= 0diffracted modes (see (13)) when a/λ0=1/1±sinθinc. Also, we would expect Wood’s anomalies for the m=0,n=±1 modes when b/λ0=1/cosθinc. When a=b, which is the particular case treated in Fig. 3(a), the Wood’s anomalies appear for a/λ0<1 and a/λ0>1, i.e., above and below the frequency of the Wood anomaly for normal incidence (a/λ0=1). Note that total transmission does not occur in the EOT peaks associated with the Wood’s anomalies for which a/λ0>1, which is due to the fact that some diffracted modes have already been launched at those frequencies and are capturing part of the energy of the original incident wave. However, Wood’s anomalies and EOT peaks may not always appear for oblique incidence at frequencies below the frequency for which λ0=max(a,b). In particular, Fig. 3(b) shows the transmission curves of a2-D periodic array of rectangular slots illuminated by obliquely incident plane waves propagating along different directions given by (φinc =π/2,θ inc). For these particular incidence directions, we should expect Wood’s anomalies for the m=±1,n=0diffracted modes when a/λ0=1/cosθinc, and for the m=0,n=±1diffracted modes when b/λ0=1/1±sinθinc. In the case shown in Fig. 3(b) where a=b, the Wood anomalies and EOT peaks for the m=±1,n=0 modes appear when a/λ0>1, but those for the m=0,n=±1 modes do not appear when a/λ0<1. The explanation for this latter behavior is hidden in Eqns. (10) and (11). Under the conditions studied in Fig. 3(b), it turns out that kx0=0. This means that when Eqn. (13) is fulfilled for the m=0,n=±1 modes, there is a zero-pole cancellation in all the elements of the matrix  GM(kx=kx0,ky=ky,±1) except for the y−yelement. Also, under the illumination conditions studied in Fig. 3(b), the functions bj(x,y)·ˆ yare odd functions of xsince the planes x=ma +a/2(m=...,−1,0,1,...) of the 2-D periodic array are all EWs [17]. As a consequence of this, it turns out that  bj(kx0=0,ky,±1)·ˆ y=0. So, when  GM(kx=kx0,ky=ky,±1) and  bj(kx0=0,ky,±1) are both substituted in (10), there is an additional zero-pole cancellation that prevents Γ ij from being infinite when m=0,n=±1 and condition (13) is fulfilled. Note that total transmission does not occur in the EOT peaks associated with the m=±1,n=0 modes in Fig. 3(b) when θinc 0, which indicates that the m=0,n=±1diffracted modes have been launched, even though the Wood’s anomaly and EOT peaks for these modes have not appeared. MoM MoM MoM TE1 Einc Fig. 4. Transmission spectra obtained when the TE1mode of the PPWG of Fig. 1 impinges on the 1-D array of rectangular slots. Our results (MoM) are compared with HFSS results. The dimensions of the unit cell were chosen to be ws/a=0.05 and a=b. The MoM formulation of (4) to (12) and the superposition principle have both been used to compute the transmission through the 1-D array of rectangular slots of Fig. 1 under TEz 1mode Vol. 25, No. 20 | 2 Oct 2017 | OPTICS EXPRESS 24676 excitation. The results are plotted in Fig. 4. For validation purposes, our MoM results have been compared with results obtained by the commercial software HFSS, and excellent agreement has been found. It should be pointed out that our in-house software is about 225 times faster than HFSS for the generation of Fig. 3, which justifies the efforts carried out to develop inhouse software in addition to the physical insight provided for the understanding of the Wood’s anomaly phenomena. Taking into account the Wood’s anomalies detected in Fig. 3(b) for the diffracted modes m=±1,n=0, we have inserted these two values of mand nin (13) plus the incidence angle values (φinc =π/2|3π/2,θ inc =θ0) (where θ0is defined in (2)) for the TEz 1mode in the PPWG, and we have found that a Wood’s anomaly and the associated EOT peak is to be expected when a/λ0=b/λ0=√5/2>1(a=bhas been assumed in Fig. 4). This exactly coincides with the results obtained in Fig. 4, and it can be verified that this frequency estimation for the Wood’s anomaly roughly matches those experimentally obtained by Reichel et al. [15] under TEz 1mode excitation. There is another way to explain the appearance of the Wood’s anomaly for a/λ0=b/λ0= √5/2>1, which is based on the fact that the planes x=ma and x=ma +a/2(m= ...,−1,0,1,...) act as EWs and the plane y=b/2 acts as a magnetic wall (MW) in Fig. 1 under TEz 1mode excitation. This means the periodic unit cell acts as an equivalent rectangular waveguide excited by a TE01 mode. The symmetry conditions imposed by the EWs and the MWs mean that the first higher order TM mode that can be excited in this equivalent waveguide is a TM21 mode. Since TM modes in rectangular waveguides are precisely those which give rise to Wood’s anomalies and EOT peaks [12,14], it is clear that the first Wood’s anomaly of the PPWG under TEz 1mode excitation will show up at the cutofffrequency of the TM21 mode, which occurs for a/λ0 = b/λ0 = √5/2 when a = b. This is in agreement with previous studies of EOT in hollow metal waveguides [12]. When a/λ0 = b/λ0 > √5/2, a combination of TE1 and TM1 modes will be excited in the PPWG propagating at oblique angles from the z direction since TM21 modes cannot propagate in a PPWG. 4. Conclusion In summary, we have presented a theory based on the application of the Method of Moments that explains why the EOT peaks of PPWG with 1-D periodic arrays of symmetric holes under TE1z mode excitation appear at frequencies higher than those found for TEM mode excitation. The theory indicates that EOT peaks for TE1z excitation should also be expected at frequencies lower than those of TEM mode excitation, but they do not appear because of the symmetry of the holes. The symmetry plane perpendicular to the x axis eliminates the divergent response of the fields on the first onset of diffraction, therefore finding the first peak associated to extraordinary transmission at higher frequencies than at normal incidence. Funding Engineering and Physical Sciences Research Council (EPSRC) (EP/L015331/1) and “Junta de Andalucía” (P12-TIC-1435). Acknowledgment All data created during this research are openly available from the University of Exeter’s institutional repository at https://ore.exeter.ac.uk/. Vol. 25, No. 20 | 2 Oct 2017 | OPTICS EXPRESS 24677