Transmission of electromagnetic waves through sub-wavelength channels

We propose a method of tunneling electromagnetic (EM) waves through a channel with sub-wavelength cross section. By filling the channel with high-ε isotropic material and implementing two matching layers with uniaxial metamterial substrates, the guided waves can go through the narrow channel without being cut off, as if it has just passed through the original empty waveguide. Both the magnitude and phase information of the EM fields can be effectively restored after passing this channel, regardless of the polarization of the incoming wave. The performance of this subwavelength channel, which is designed with coordinate transformation methodology, is studied theoretically and numerically. ©2010 Optical Society of America OCIS codes: (260.2110) Electromagnetic optics; (160.3918) Metamaterials; (350.4010) Microwaves. References and links 1. H. T. Miyazaki, and Y. Kurokawa, “Squeezing visible light waves into a 3-nm-thick and 55-nm-long plasmon cavity,” Phys. Rev. Lett. 96(9), 097401 (2006). 2. M. Silveirinha, and N. Engheta, “Tunneling of electromagnetic energy through subwavelength channels and bends using ε-near-zero materials,” Phys. Rev. Lett. 97(15), 157403 (2006). 3. R. Liu, Q. Cheng, T. Hand, J. J. Mock, T. J. Cui, S. A. Cummer, and D. R. Smith, “Experimental demonstration of electromagnetic tunneling through an epsilon-near-zero metamaterial at microwave frequencies,” Phys. Rev. Lett. 100(2), 023903 (2008). 4. B. Edwards, A. Alù, M. E. Young, M. Silveirinha, and N. Engheta, “Experimental verification of epsilon-nearzero metamaterial coupling and energy squeezing using a microwave waveguide,” Phys. Rev. Lett. 100(3), 033903 (2008). 5. M. G. Silveirinha, and N. Engheta, “Transporting an Image through a Subwavelength Hole,” Phys. Rev. Lett. 102(10), 103902 (2009). 6. J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science 312(5781), 1780–1782 (2006). 7. U. Leonhardt, “Optical conformal mapping,” Science 312(5781), 1777–1780 (2006). 8. J. Li, and J. B. Pendry, “Hiding under the carpet: a new strategy for cloaking,” Phys. Rev. Lett. 101(20), 203901 (2008). 9. R. Liu, C. Ji, J. J. Mock, J. Y. Chin, T. J. Cui, and D. R. Smith, “Broadband ground-plane cloak,” Science 323(5912), 366–369 (2009). 10. J. Valentine, J. Li, T. Zentgraf, G. Bartal, and X. Zhang, “An optical cloak made of dielectrics,” Nat. Mater. 8(7), 568–571 (2009). 11. L. H. Gabrielli, J. Cardenas, C. B. Poitras, and M. Lipson, “Silicon Nanostructure Cloak Operating at Optical Frequencies,” Nat. Photonics 3(8), 461–463 (2009). 12. H. Chen, and C. T. Chan, “Transformation Media that Rotate Electromagnetic Fields,” Appl. Phys. Lett. 90(24), 241105 (2007). 13. J. Zhang, Y. Luo, H. Chen, and B. I. Wu, “Manipulating the directivity of antennas with metamaterial,” Opt. Express 16(15), 10962–10967 (2008), http://www.opticsinfobase.org/oe/abstract.cfm?uri=oe-16-15-10962. 14. Y. Luo, J. Zhang, L. Ran, H. Chen, and J. A. Kong, “New Concept Conformal Antennas Utilizing Metamaterial and Transformation Optics,” IEEE Antennas Wirel. Propag. Lett. 7, 509–511 (2008). 15. M. Rahm, D. A. Roberts, J. B. Pendry, and D. R. Smith, “Transformation-optical design of adaptive beam bends and beam expanders,” Opt. Express 16(15), 11555–11567 (2008), http://www.opticsinfobase.org/abstract.cfm?URI=oe-16-15-11555. 16. D. H. Kwon, and D. H. Werner, “Polarization splitter and polarization rotator designs based on transformation optics,” Opt. Express 16(23), 18731–18738 (2008), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-1623-18731. #120736 $15.00 USD Received 1 Dec 2009; revised 22 Jan 2010; accepted 4 Feb 2010; published 11 Feb 2010 (C) 2010 OSA 15 February 2010 / Vol. 18, No. 4 / OPTICS EXPRESS 3864 17. H. Ma, S. Qu, Z. Xu, and J. Wang, “General method for designing wave shape transformers,” Opt. Express 16(26), 22072–22082 (2008), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-26-22072. 18. Y. Luo, J. Zhang, H. Chen, J. Huangfu, and L. Ran, “High-directivity Antenna with Small Antenna Aperture,” Appl. Phys. Lett. 95(19), 193506 (2009). 19. Y. Luo, J. Zhang, H. Chen, L. Ran, B.-I. Wu, and J. A. Kong, “A Rigorous Analysis of Plane-transformed Invisibility Cloaks,” IEEE Trans. Antenn. Propag. 57(12), 3926–3933 (2009). 20. J. Pendry, “Taking the Wraps off Cloaking,” Physics 2, 95 (2009).


Introduction
Due to its potential applications in a wide range of fields, the transmission of electromagnetic energy through sub-wavelength channels has intrigued an extensive investigation in the recent years. The technology of exciting surface plasmon polaritons on both sides of the channel explores the possibility of overcoming the barrier of miniaturization and designing subwavelength waveguides [1]. Another inventive approach to the problem of squeezing electromagnetic energies through sub-wavelength channels was also introduced [2], making use of ε-near zero (ENZ) materials. This idea was later experimentally verified in microwave frequencies by two groups. One uses a planar waveguide setup, where enhanced transmission was measured when the narrow channel connecting waveguides was filled with patterned complementary split ring resonators [3]. The other achieved energy squeezing and tunneling through an ultra-narrow rectangular waveguide channel that mimicked zero-permittivity properties [4].The scheme of utilizing ENZ materials to improve the energy transmission relies on the principle that there are no relevant reflection losses at abrupt junctions or bends where the permittivity of the materials is extremely small, and has the main advantage that the reflection reduction effect is independent of the shapes of the junctions or bends. However, this method cannot be used to guide EM fields through a channel where EM wave is inherently cut off. In addition, it is not effective for tunneling waves in high modes. In Ref [5], this method was improved when arrays of metallic wires were incorporated and was applied to transport an electromagnetic image through a tiny hole in the sub-wavelength spatial scale. In this paper, we suggest another possible solution to the problem of sub-wavelength tunneling, which is based on the coordinate transformation methodology. By exploiting a transformation which compresses the original space into a much narrower region, EM waves which are supposed to transmit in the original space will be confined in a channel whose cross-section is much smaller than the free-space wavelengths. This transformation gives rise to a set of materials parameters which could be achieved using natural positive materials as well as uniaxial metamaterials. Compared with the ENZ material approach in Ref [2]. which only works for transverse electromagnetic (TEM) mode, the channel designed with transformation methodology can potentially squeeze and support all the incoming waves regardless of the polarizations and the transmission modes. Not only the magnitude but also the phase information of EM waves will be restored at the output interface of the channel. Both theoretical analysis and numerical simulations demonstrate the effectiveness of this subwavelength tunneling scheme.

Principle
Transformation optics offers the flexibility of controlling the paths of EM waves at will and has recently been widely used to design many EM devices with extraordinary capabilities, such as invisibility cloaks, field rotators, conformal antennas, highly directive antennas, beam expanders, etc [6][7][8][9][10][11][12][13][14][15][16][17][18][19]. Here we consider using the knowledge of the transformation to obtain the values of material parameters required to compress and guide the EM waves in a subwavelength channel. For simplification, we confine our discussion to the two-dimensional (2D) case, where the fields are invariant in z direction. As depicted in Fig. 1, a trapeziform region with the height h 2 , the short side 2l, and the long side 2(l + d), has been compressed into the blue-gridded polygonal region. The upper and lower boundaries of the original and the transformed space are all perfect electric conducting (PEC) boundaries. The following equations describe such a transformation: In fact, by applying this transformation, the curved perfect reflecting surface marked in red in Fig. 1 will behave as a conducting plane (denoted by green solid line) for EM wave illuminations. This is similar to the idea of the "carpet cloak" [8][9][10][11]19], where an obstacle hidden under a conducting surface would be perceived as a flat conducting sheet. Therefore, after the waves travel through the blue-gridded region with conducting boundaries, they will go on propagating as if they just pass a rectangular region with reflecting surfaces. This explains why the waves transmitted along the x axis can still be supported even when the value of h 2 −h 1 is far less than the free-space wavelength, and all the transmission modes will be restored at the far side of the transformed region. The corresponding permittivity and permeability tensors of the transformed region can be expressed as [ ε , µ appear on the same footing because of the symmetry between electric and magnetic fields. The parameter tensors with elements in Eqs. (2-a) and (2-b) can be easily diagonized by rotating the coordinate system, and can be realized with uniaxial materials. After the diagonization process, all the components of the material parameters become positive constants. In addition, by compressing the space in the central tunneling region along the x axis by a factor of (h 2 −h 1 )/h 2 [20], we can further obtain a set of even simpler parameters: For a certain TE (or TM) incidence, only ε zz , µ xx , µ yy (or µ zz , ε xx , ε yy ) enter into the Maxwell's equations, so the material in the channel can be replaced by the isotropic material with parameters ε = [h 2 /(h 2 −h 1 )] 2 , µ = 1 (or ε = 1, µ = [h 2 /(h 2 −h 1 )] 2 ), while still allowing for the tunneling effect.

Numerical simulations
The COMSOL Multiphysics finite element-based electromagnetics solver is used to verify the above idea and illustrate the wave propagation. In all the simulations, the working frequency of the EM waves is for simplicity 2GHz. We first consider a case where only the first transverse-electric (TE1) mode is supported by a 0.1m (h 2 = 0.1m) wide planar waveguide, as shown in Fig. 2(a). When a channel with a width of only 0.02m (h 2 −h 1 = 0.01m) is embedded in the waveguide (the boundaries of the channel is PEC, and l = 0.1m, d = 0.05m for all the following simulations), the EM wave is cut off and the fields vanish in the wedgy junction region between the waveguide and the channel, as depicted by Fig. 2(b). Next we fill the narrow channel as well as the junction regions on the two sides with the materials characterized by Eq. (2). The result in Fig. 2(c) shows that all the fields are guided through the narrow channel to the planar waveguide on the other side. Compared with the electric field distributions shown in Fig. 2(a), we see that the EM waves continue propagating after passing the channel as if they have never met any blockage in the original waveguide. In Fig. 2(d), all the simulation setups are the same as those in Fig. 2(c) except that we use the material characterized by Eq. (3) for the channel, and similar results could be observed. To understand the tunneling effect from a physical perspective, we may need to check the dispersion relation 2 2 2 x y r r k k ω ε µ + = , µ = 1, k y = mπ/d, m = 1, 2, 3…. where d is the width of the planar waveguide. When ε r increases to a value large enough, k x will change from imaginary to real. In other words, the cut-off frequency for the waveguide will decrease if it is stuffed with materials with high permittivity or permeability. In the simulation shown in Fig. 2(d), the material in the channel has an ε equal to 100, so the propagation of waves would still be supported even though the width of the channel is only 1/15 of the freespace wavelength. Next, we address whether we can achieve the tunneling effect by simply inserting highpermittivity medium in the channel. Figure 2(e) gives the answer, showing that the EM wave still cannot be transmitted in this case. This is because the impedance is strongly mismatched at the junction of waveguides with different cross-sections. In fact, the layers on both sides of the channel (-(l + d) < x < -l and l < x < (l + d)) are used for the purpose of impedance matching. The performance of the matching layers can be studied theoretically as follows.
To simplify our calculation, we first suppose there is no reflection from the interface between the air and the matching layer and only consider the interface between the matching layer and the channel.
Then the problem will be simplified to the case shown in Fig. 3(a) where a TE wave is incident upon the interface between an isotropic medium and an anisotropic medium. The incident electric field, reflected electric field and the transmitted electric field can be expressed as , .
Since the permeability tensor in our case is symmetric (µ xy = µ yx ), we have It is easy to find that as long as µ xx ε zz = µ 1 ε 1 , Eq. (7) is independent of k y and can be further simplified to 2 2 1 xx yy xy µ µ µ µ − = .
Since in the case we discussed the condition µ 1 = 1 holds, the total transmission condition can now be summarized as 2 1 1, .
xx yy xy xx zz The fact that Eq. (8) is satisfied for material parameters of the matching layer at the interface implies there is no wave reflection at this interface. The interface between the matching layer and the air can be similarly treated except that we need to rotate the coordinate system ([x, y, z) by an angle -φ to a new system (x', y', z'), where  From the above analysis, we can see that our proposed approach is able to realize TE mode wave tunneling. In fact, for the TEM wave case, which can be dealt with using ENZ material, our approach is also applicable. Figure 4(a) is the electric field distribution for a waveguide supporting TEM wave, and Fig. 4(b) shows the electric field distribution when a narrow channel is built inside the waveguides. Here h 2 = 0.06m, h 1 = 0.05 m. The reflection introduced by the discontinuity of the boundary can be observed by comparing the magnitude of electric fields on the two sides of the channel. In Fig. 4(c) and 4(d), we simulate the tunneling performance of the ENZ material method and our transformation method, respectively. However, different with the ENZ material method which can only be applied to TEM polarization cases, our approach can also treat EM waves in higher modes.
In the following simulation, we increase the width of the planar waveguides to 0.3m so that TM3 mode can be supported. In order to observe the high mode clearly, we only let TM3 mode be excited in our simulation, as displayed in Fig. 5(a). When a 0.02m wide channel is constructed inside such a planar waveguide, the TM3 mode will be cut off at the narrow channel and some lower mode will be excited at the discontinuous boundary. Therefore, some energy is transmitted through the channel without maintaining the original mode, which can be observed in Fig. 5(b). After we fill the channel with high-ε material as well as adding the two matching layers, the waves can be guided smoothly through the channel and continue propagating in the planar waveguide with TM3 mode [see Fig. 5(c)]. As a comparison, we also simulate how the ENZ material scheme would operate in this case. The result in Fig. 5(d) demonstrates that the ENZ material cannot achieve the tunneling effect for the waves of higher modes.

Conclusion
A new method which allows for the EM wave squeezing and transmission trough narrow channels and waveguide junctions, was introduced in this paper. The validity of this approach was theoretically analyzed and numerically tested in different polarization conditions, using waveguides with different cross sections. Instead of being cut off, all the energies can be totally tunneled through the narrow channel independent of the polarization of the incoming waves. The information of transmission modes will also be persevered after the tunneling operation. Transformation optics provides the theoretical background of this scheme and uniaxial metamaterials provide the means of achieving the prescribed parameters.