## Abstract

We suggest a principle for isolation of circularly polarized waves in magnetically active extreme-parameter metamaterials. Using theoretical analysis and numerical simulations, we show that metamaterials with extreme parameters, such as epsilon-near-zero materials (ENZ), when merged with magneto-optical materials, become transparent for forward circularly polarized waves of a given handedness and opaque for backward propagating waves of the same handedness. We theoretically study two possible implementations of such hybrid materials: (1) the case of metal-dielectric stacks; and (2) rectangular waveguide near its cut-off frequency. We prove that these structures can be utilized as compact isolators for circularly polarized waves.

© 2013 OSA

## 1. Introduction

Many objects in nature ranging from molecules and biological species to astrophysical bodies exhibit different sensitivity to left- and right-handed circular polarizations of light [1–3]. Therefore the ability to filter and isolate circular polarizations is of practical importance for science and engineering. Progress with chiral heterostructures [4–8], paved the way for the design and engineering of ultra-compact composite structures distinguishing left- and right-handed polarization states. Nevertheless, any chiral metastructure satisfying the Lorentz reciprocity theorem [9] is reciprocal, implying that interaction of the backward and forward propagation of the same circularly polarized states with such structures are identical and can not be distinguished. Imagine, for example, a right-handed circularly polarized (RCP) light incident normally in +*z* (forward) direction on a slab of chiral material. Its transmission and reflection may differ from left-handed circularly polarized (LCP) light incident from the same direction. At the same time due to reciprocity of the system RCP light propagating backward (in −*z* direction) experiences exactly the same transmission and reflection as forward propagating RCP, see schematic in Fig. 1(a). Hence, clearly the isolation and separation of forward and backward waves is not possible in such structures. It is worth mentioning here the works on asymmetric transmission with some types of chiral heterostructures [10–12], however these structures still satisfy the Lorentz reciprocity conditions and can not be considered as truly nonreciprocal systems with asymmetric transmission.

Difference in forward and backward light flow directions is closely related with the time reversal symmetry breaking of Maxwell’s equations [13,14] and nonreciprocity of electromagnetic laws [9, 15]. The latter is possible with either time noninvariant, gyrotropic, or nonlinear materials. In particular, recent advances with magneto-photonic crystals [16, 17] and nonlinear dielectric stacks [18, 19] demonstrate the possibility for one-way transmission of linearly polarized waves in such systems. Furthermore, an optical isolation of circularly polarized waves was investigated in both nonlinear chiral metastructures [20] and in a thin magneto-photonic crystal slab [21]. However, both systems have limitations in their applicability; for example, the nonlinear chiral metastructure operates at microwave frequencies [20], and the extension of this interesting idea to higher frequencies would be challenging due to the weak optical nonlinearities. On the other hand, the elegant idea of optical isolation in perforated magneto-optical slabs [21] might experience difficulties with fabrication due to the mechanical stability of magneto-optical materials.

Merging the concepts of metamaterial design with magneto-optical (MO) response might open possibilities for engineering novel nonreciprocal systems for optical isolation and signal processing. In particular, recent advances with nonreciprocal magneto-active metastructures [22–25], including plasmonic crystals [26–28], have shown a significant enhancement of Faraday and Kerr effects. Metamaterials – artificially engineered heterostructures with elements much smaller than operating wavelengths [29, 30] – can be treated macroscopically as a homogeneous media exhibiting exotic values for permittivity and permeability, such as negative refractive index [31–33] naturally not observed in nature. Recently, materials with epsilon-near-zero (ENZ) permittivity were proposed and fabricated [34–36]. One of the most intriguing properties of these materials is the nearly zero phase progression of propagating waves, which can be employed for the efficient energy tunneling [34,36] and for the design of the next generation of waveguides [37]. What is more important in these materials is that a small variation of their parameters may allow tuning their properties form being opaque “metals” (with Re(*ε*) < 0) to transparent dielectrics (when Re(*ε*) > 0).

In this paper, we suggest a design approach for optical isolation for circularly polarized waves by utilizing both the epsilon-near-zero materials or structures and the magneto-optical materials. We show theoretically that magnetically active ENZ media is transparent for one of the circularly polarized waves, whereas it is opaque for another polarization handedness. We propose two realistic examples where such a scenario is feasible. Specifically, we study the case of a stack of alternating metal-dielectric slabs, and the case of a waveguide filled with magneto-optical media operating near its cut-off frequency. We demonstrate numerically that the isolation of circularly polarized states is possible in both of these example structures.

The paper is organized as follows. In Section II we give a generalized concept for designing the ENZ-MO isolator for circularly polarized waves. We discuss two particular examples of such isolator in Section III. In particular, we consider a metal-dielectric stack and a waveguide near its cut-off, and show that both structures exhibit effective ENZ properties combined with gyrotropy allowing isolation of one of the circularly polarized states when magneto-optical activity is present. Section IV provides the conclusion.

## 2. Heuristic approach for the ENZ isolator deign

We begin our analysis with a brief insight into the properties of a conventional magneto-optical materials. Consider, for example, a slab of magneto-optical medium with a magnetization direction normal to the slab interface, see Fig. 1(b). In this case the MO relative permittivity can be described by an antisymmetric tensor [9]:

*ε*and

_{mo}*ε*

_{⊥}are diagonal components of relative permittivity tensor,

*α*> 0 is the off-diagonal component of the permittivity tensor responsible for the “strength” of magneto-optical activity of the medium. At optical frequencies, the parameter

*α*is often very small and usually of the order of 10

^{−2}or even smaller.

For plane waves propagating along the magnetization direction (i.e. in Faraday configuration) the eigenmodes in magneto-optical media are two circularly polarized waves: *E _{x}* ±

*iE*. Solving Maxwell’s equations for the circularly polarized waves propagating in the magneto-optical media in forward direction we get a well known relation for the corresponding wavevectors [9]: ${k}_{LCP}^{+}=\left(\omega /c\right)\sqrt{{\epsilon}_{mo}-\alpha}$ for left-handed circular polarized waves and ${k}_{RCP}^{+}=\left(\omega /c\right)\sqrt{{\epsilon}_{mo}+\alpha}$ for right-handed circular polarized waves. To account for the backward propagating waves the time reversal operation (

_{y}*t*→ −

*t*) needs to be performed. In this case the change of sign in front of

*α*has to be accounted for in the wave dispersions: ${k}_{LCP}^{-}=-\left(\omega /c\right)\sqrt{{\epsilon}_{mo}+\alpha}$ and ${k}_{RCP}^{-}=-\left(\omega /c\right)\sqrt{{\epsilon}_{mo}-\alpha}$ for backward LCP and RCP waves, respectively. Clearly, backward and forward propagating circular polarized waves of the same handedness have different wave numbers with differing dispersion (|

*k*

^{+}| ≠ |

*k*

^{−}|) yielding difference in their transmittance and reflectance [9].

Exploiting the metamaterial concept, one can engineer metamaterials and structures in order to deliberately tailor the value of permittivity function [30]. In Refs. [36–38] it was shown that in properly designed structures, one can achieve *ε* → 0, and consequently very small phase progression (*k* → 0). Examples of such structures include waveguides near their cut-off frequency for the TE modes [36], layered metal-dielectric stacks [38] and natural materials at frequencies for which the real part of their permittivity is near zero (while its imaginary par is small) [39]. By integrating such epsilon-near zero structures with magneto-optical materials it may be possible to design effective magneto-active media with nulled diagonal components in the permittivity tensor, i.e. *ε _{mo}* → 0 in Eq.(1). Assuming now that such a case is possible, we obtain following dispersion relations for circularly polarized wave eigenmodes:

Assuming a lossless condition, as evident from Eq.(2), for *α* > 0 the forward propagating LCP waves and backward propagating RCP waves have imaginary wavevectors in such a medium, implying that the medium is opaque for these waves, while for the two other eigen-modes, the wave numbers are real, and thus the medium is transparent to these waves. Consequently, it would be possible to achieve two effects simultaneously: 1) isolate forward and backward propagating waves with the same handedness (see schematic in Fig. 1(c)), and 2) isolate LCP from RCP waves propagating in the same direction.

In Fig. 1(d) we plot the transmittance of the RCP wave in forward and backward directions through a slab made of such material with the variation of slab thickness. To be more specific we consider here the operating wavelength to be *λ* = 1.5*μ*m, ambient medium is air and *α* = 0.1. Note that, *α* is usually smaller than this value, for example, for bismuth iron garnet, a commonly used magneto-optical material, *α* = 0.06 at *λ* = 1.55*μ*m. However, here we intentionally consider slightly higher values of magneto-optical activity in order to stress and highlight the observed effects. We see that the transmission of backward RCP wave exponentially decays as the slab thickness increases, and for the 2*μ*m thick slab it already falls below 1%. At the same time the slab acts as a Fabry-Perot resonator for the forward propagating RCP with a first transmission resonance at ∼ 2.3*μ*m thick slab. Analogously the forward propagating LCP would be suppressed similarly to backward RCP. Hence, engineering a magnetically active ENZ slab would allow to achieve isolation of circularly polarized waves with a slab thickness comparable with the incident wavelength. In the following section we propose two possible realizations of such magneto-active ENZ isolators.

## 3. Design of ENZ isolator

In order to engineer the magneto-active ENZ media we follow the two methods developed in the area of ENZ metamaterials. In the first approach we consider a metal-dielectric stack [38], shown schematically in Fig. 2(a). For the second method, a rectangular waveguide operating at the cutoff frequency of its TE mode is utilized [36], see Fig. 2(b).

#### Metal-dielectric stack

We consider a stack of periodically alternating metal and magneto-optical slabs with thicknesses *a* and *b*, respectively (Fig. 2(a)). When the thickness of the slabs is much smaller than the wavelength of operation the effective-medium approach may be used [38]. Using the Cartesian coordinate systems shown in Fig. 2, the *x* and *y* components of the effective medium permittivity for an unbiased system, i.e. when *α* = 0, are given by: *ε _{xx}* =

*ε*= (

_{yy}*aε*+

_{met}*bε*)/(

_{MO}*a*+

*b*). By properly choosing these parameters, we can achieve near-zero permittivity for

*ε*and

_{xx}*ε*. Then according to the predictions made in Section II we expect that such a stack, when magnetically biased, behaves as an isolator for circularly polarized waves. To prove our predictions we perform full-wave numerical simulations with the CST Microwave Studio

_{yy}*[40].*

^{TM}For the purpose of numerical demonstration here we consider the magneto-optical layers to be made of bismuth iron garnet with *ε _{MO}* = 6.25 and

*α*= 0.1, and the metal to be lossless with relative permittivity

*ε*= −12.5. In this case choosing

*a*= (1/2)

*b*we obtain zero effective epsilon for

*ε*and

_{xx}*ε*. A plane wave with free space wavelength of 1.5

_{yy}*μ*m is incident normally on a stack containing 64 unit cells. We choose the thickness of the metal layer to be

*a*= 20nm and that of magneto-optical slab to be

*b*= 40nm, so that the overall stack thickness of the entire structure is 3.84

*μ*m. Such choice of slab thicknesses allows the applicability of the effective medium approximation, i.e.

*a*+

*b*≪

*λ*.

Figures 3(a) and 3(b) show instantaneous distributions of the electric field vector for the RCP and LCP forward propagating plane waves at the specific time instant (*ωt* = 0). We observe that RCP wave propagates through the stack and preserves its amplitude, Fig. 3(a), whereas LCP is almost completely reflected and only small fraction of it is tunneled trough the stack, Fig. 3(b). Note that, the incident LCP and its reflection from the stack form a standing circularly polarized wave. Instantaneously it looks like linearly polarized wave, however the polarization plane rotates with time (not shown here). To analyze the degree of polarization we plot the front view of the electric field vector, i.e. its projection on the plane transverse to the direction of propagation (*x* − *y* plane) at the same specific time instant (*ωt* = 0), Figs. 3(c) and 3(d). Clearly, RCP preserves its polarization handedness with propagation, i.e. stays circular after the transmission through the stack, Fig. 3(c). At the same time reflected LCP, as was noted, forms an instantaneous linearly polarized state (which is rotating with time), Fig. 3(d). We note that the backward propagating RCP behaves similarly to the forward LCP, i.e., in particular it is reflected from the stack, indicating that the stack acts as an isolator for circularly polarized waves.

#### Waveguide at cut-off frequency

It is known that a metallic rectangular waveguide with the perfectly conducting walls, near its cut-off frequency for TE10 mode effectively behaves as an ENZ strcture [36]. We show that such a waveguide when filled with a magneto-optical material may also act as an ENZ isolator.

We consider a rectangular waveguide with a square cross-section (10*μ*m×10*μ*m) filled with bismuth iron garnet (*ε _{mo}* = 6.25 and

*α*= 0.1), see Fig. 2(b). Note that the direction of magnetization in this case is along the waveguide, i.e. along the direction of mode propagation. The square cross section of the waveguide leads to the degeneracy of the lower order TE modes (

*TE*

_{01}and

*T E*

_{10}). For the given waveguide dimensions both modes exhibit cut-off at 6THz (corresponding wavelength of 50

*μ*m).

When magneto-optical activity is introduced these degenerate states couple and form two new conjugate eigen-states *ψ*_{1} ± *iψ*_{2}. Note that such coupling is possible only between orthogonal states existing simultaneously in the unbiased waveguide, for example, a single mode rectangular waveguide does not allow any rotating states. These new conjugate modes can be considered as analogues of LCP and RCP plane waves as the eigen-modes of the bulk media. Using this analogy we define these states as left-rotating and right-rotating eigen-modes respectively. The dispersion relations defining the eigenvalues of right and left rotating states are given as:
${k}_{\pm}=\sqrt{\left({\omega}^{2}/{c}^{2}\right)\left(\epsilon \pm \alpha \right)-{\pi}^{2}/{a}^{2}}$, where *a* is the waveguide’s edge dimension. Hence near the cut-off frequency, when
$\left(\omega /c\right)\sqrt{\epsilon}=\pi /a$, we arrive to the conditions similar to that discussed in Section II. In particular, we see that near the cut-off frequency, one of the rotating states is not propagating, whereas the same rotating state propagates in the reverse direction.

To demonstrate this property of electromagnetic wave isolation inside the rectangular waveguide filled with the MO material, we perform a series of full-wave numerical calculations. In our simulations we assume that the waveguide walls are made of perfect electric conductor (PEC). We excite the waveguide with *TE*_{10} mode and study the corresponding field distribution and transmission characteristics. For magneto-optically active case the *TE*_{10} is decomposed into the counter-rotating eigen-modes supported by magnetized waveguide. Above the cut-off frequency both counter-rotating states exist, and due to the difference between their wavenumbers the phase between the *TE*_{10} and *TE*_{01} is accumulated with propagation along the waveguide, similarly to the Faraday effect in bulk media [9]. Below cut-off only one of the counter-rotating states exists, hence *TE*_{10} launched into the waveguide would generate this rotating state.

We trace the time variation of electric field at several different positions along the waveguide, i.e. along the central line, see Figs. 4(a) and 4(b). Fig. 4(a) shows the electric field variation above cut-off at 6.1 THz. We observe that in this case the polarization is highly elliptical, almost linear. Furthermore, we notice rotation of the polarization axis as we move along the waveguide as predicted earlier. Below the cut-off at 5.98THz, see Fig. 4(b), we observe propagation of purely circularly polarized state, which clearly manifests the isolation of one of the rotating states. Finally, we note that the small deviation from the linear polarization on Fig. 4(a) and modulation of the amplitude on Fig. 4(b) is due to small reflection appearing due to the numerical mismatch at the input and output ports of the waveguide.

## 4. Conclusion

We showed that magnetically active epsilon-near-zero structures acts as compact isolator for circularly polarized waves. We have demonstrated two cases for such media: a metal-dielectric stack and a rectangular waveguide near the cut-off frequency. Using full-wave numerical simulation, we have shown that in both structures the complete suppression of one of the circularly polarized states is possible.

## Acknowledgments

This work is supported in part by the US Air Force Office of Scientific Research (AFOSR) grant number FA9550-10-1-0408.

## References and links

**1. **W. B. Sparks, J. Hough, T. A. Germer, F. Chen, S. DasSarma, P. DasSarma, F. T. Robb, N. Manset, L. Kolokolova, N. Reid, F. D. Macchetto, and W. Martin, “Detection of circular polarization in light scattered from photosynthetic microbes,” PNAS **106**, 7816–7821 (2009). [CrossRef] [PubMed]

**2. **V. S. Degtjarev and L. O. Kolokolova, “Possible application of circular polarisation for remote sensing of cosmic bodies,” Earth, Moon and Planets **57**, 213–223 (1992). [CrossRef]

**3. **L. Kolokolova, W. Sparksm, and D. Mackowski, “Astrobiological remote sensing with circular polarization” in *Polarimetric Detection, Characterisation and Remote Sensing*, M. I. Mishchenko, Y. S. Yatskiv, V. K. Rosenbush, and G. Videen, eds. (Springer, Netherlands, 2011) pp. 277–294. [CrossRef]

**4. **M. Decker, M. W. Klein, M. Wegener, and S. Linden, “Circular dichroism of planar chiral magnetic metamaterials,” Opt. Lett. **32**, 856–858 (2007). [CrossRef] [PubMed]

**5. **E. Plum, V. A. Fedotov, and N. I Zheludev, “Optical activity in extrinsically chiral metamaterial,” Appl. Phys. Lett. **93**, 191911–3 (2008). [CrossRef]

**6. **C. Helgert, E. Pshenay-Severin, M. Falkner, C. Menzel, C. Rockstuhl, E.-B. Kley, A. Tunnermann, F. Lederer, and T. Pertsch, “Chiral metamaterial composed of three-dimensional plasmonic nanostructures,” NanoLett. **11**, 4400–4404 (2011). [CrossRef]

**7. **Y. Zhao, M. A. Belkin, and A. Alu, “Twisted optical metamaterials for planarized ultrathin broadband circular polarizers,” Nature Commun. **3**, 1–7 (2012). [CrossRef]

**8. **M. Schferling, D. Dregely, M. Hentschel, and H. Giessen, “Tailoring enhanced optical chirality: Design principles for chiral plasmonic nanostructures,” Phys. Rev. X **2**, 031010 (2012). [CrossRef]

**9. **L. D. Landau, L. P. Pitaevskii, and E. M. Lifshitz, *Electrodynamics of Continuous Media, Second Edition: Volume 8 (Course of Theoretical Physics)* (Butterworth-Heinemann, 1984).

**10. **V. A. Fedotov, P. L. Mladyonov, S. L. Prosvirnin, A. V. Rogacheva, Y. Chen, and N. I. Zheludev, “Asymmetric propagation of electromagnetic waves through a planar chiral structure,” Phys. Rev. Lett. **97**, 167401 (2006). [CrossRef] [PubMed]

**11. **A. S. Schwanecke, V. A. Fedotov, V. V. Khardikov, S. L. Prosvirnin, Y. Chen, and N. I. Zheludev, “Nanostructured metal film with asymmetric optical transmission,” NanoLett **8**, 2940–2943 (2008). [CrossRef]

**12. **C. Menzel, C. Helgert, C. Rockstuhl, E.-B. Kley, A. Tunnermann, T. Pertsch, and F. Lederer, “Asymmetric transmission of linearly polarized light at optical metamaterials,” Phys. Rev. Lett. **104**, 253902 (2010). [CrossRef] [PubMed]

**13. **F. D. M. Haldane and S. Raghu, “Possible realization of directional optical waveguides in photonic crystals with broken time-reversal symmetry,” Phys. Rev. Lett. **100**, 013904 (2008). [CrossRef] [PubMed]

**14. **Z. Wang, Y. Chong, J. D. Joannopoulos, and M. Soljacic, “Observation of unidirectional backscattering-immune topological electromagnetic states,” Nature **461**, 772 (2009). [CrossRef] [PubMed]

**15. **R. J. Potton, “Reciprocity in optics,” Rep. Prog. Phys. **67**, 717754 (2004). [CrossRef]

**16. **Z. Yu, Z. Wang, and S. Fan, “One-way total reflection with one-dimensional magneto-optical photonic crystals,” Appl. Phys. Lett. **90**, 121133–3 (2007). [CrossRef]

**17. **A. B. Khanikaev, S. H. Mousavi, G. Shvets, and Y. S. Kivshar, “One-way extraordinary optical transmission and nonreciprocal spoof plasmons,” Phys. Rev. Lett. **105**, 126804 (2010). [CrossRef] [PubMed]

**18. **M. D. Tocci, M. J. Bloemer, M. Scalora, J. P. Dowling, and C. M. Bowden, “Thin-film nonlinear optical diode,” Appl. Phys. Lett. **66**, 23246 (1995). [CrossRef]

**19. **A. E. Miroshnichenko, E. Brasselet, and Y. S. Kivshar, “Reversible optical nonreciprocity in periodic structures with liquid crystals,” Appl. Phys. Lett. **96**, 063302–3 (2010). [CrossRef]

**20. **I. V. Shadrivov, V. A. Fedotov, D. A. Powell, Y. S. Kivshar, and N. I. Zheludev, “Electromagnetic wave analogue of an electronic diode,” New J. Phys. **13**, 033025–8 (2011). [CrossRef]

**21. **K. Fang, Z. Yu, V. Liu, and S. Fan, “Ultracompact non-reciprocal optical isolator based on guided resonance in a magneto-optical photonic crystal slab,” Opt. Lett. **36**, 4254–4256 (2011). [CrossRef] [PubMed]

**22. **D. J. Bergman and Y. M. Strelniker, “Anisotropic ac electrical permittivity of a periodic metal-dielectric composite film in a strong magnetic field,” Phys. Rev. Lett. **80**, 857–860 (1998). [CrossRef]

**23. **F. J Rachford, D. N. Armstead, V. G. Harris, and C. Vittoria, “Simulations of ferrite-dielectric-wire composite negative index materials,” Phys. Rev. Lett. **99**, 057202–4 (2007). [CrossRef] [PubMed]

**24. **W. Li, Z. Liu, X. Zhang, and X. Jiang, “Switchable hyperbolic metamaterials with magnetic control,” Appl. Phys. Lett. **100**, 161108–4 (2012). [CrossRef]

**25. **P. K. Jain, Y. Xiao, R. Walsworth, and A. E. Cohen, “Surface plasmon resonance enhanced magneto-optics (SuPREMO): Faraday rotation enhancement in gold-coated iron oxide nanocrystals,” NanoLett. **9**, 1644–1650 (2009). [CrossRef]

**26. **A. B. Khanikaev, A. V. Baryshev, A. A. Fedyanin, A. B. Granovsky, and M. Inoue, “Anomalous Faraday effect of a system with extraordinary optical transmittance,” Opt. Express **15**, 6612–6622 (2007). [CrossRef] [PubMed]

**27. **V. I. Belotelov, I. A. Akimov, M. Pohl, V. A. Kotov, S. Kasture, A. S. Vengurlekar, A. V. Gopal, D. R. Yakovlev, A. K. Zvezdin, and M. Bayer, “Enhanced magneto-optical effects in magnetoplasmonic crystals,” Nature Nanotech. **6**, 370–376 (2011). [CrossRef]

**28. **J. Y. Chin, T. Steinle, T. Wehlus, D. Dregely, T. Weiss, V. I. Belotelov, B. Stritzker, and H. Giessen, “Nonreciprocal nanoplasmonics: Giant enhancement of thin film Faraday rotation,” (to be published).

**29. **D. R. Smith, J. B. Pendry, and M. C. K. Wiltshire, “Metamaterials and negative refractive index,” Science **305**, 788–792 (2004). [CrossRef] [PubMed]

**30. **C. M. Soukoulis and M. Wegener, “Past achievements and future challenges in the development of three-dimensional photonic metamaterials,” Nature Photon. **5**, 523–530 (2011).

**31. **D. R. Smith, W. J. Padilla, D. C. Vier, S. C. Nemat-Nasser, and S. Schultz, “Composite medium with simultaneously negative permeability and permittivity,” Phys. Rev. Lett. **84**, 4184–4187 (2000). [CrossRef] [PubMed]

**32. **V. M. Shalaev, W. S. Cai, U. K. Chettiar, H. K. Yuan, A. K. Sarychev, V. P. Drachev, and A. V. Kildishev, “Negative index of refraction in optical metamaterials,” Opt. Lett. **30**, 3356–3358 (2005). [CrossRef]

**33. **J. Valentine, S. Zhang, T. Zentgraf, E. Ulin-Avila, D. A. Genov, G. Bartal, and X. Zhang, “Three-dimensional optical metamaterial with a negative refractive index,” Nature **455**, 376–U32 (2008). [CrossRef] [PubMed]

**34. **M. Silveirinha and N. Engheta, “Tunneling of electromagnetic energy through subwavelength channels and bends using epsilon-near-zero materials,” Phys. Rev. Lett. **97**, 157403–4 (2006). [CrossRef] [PubMed]

**35. **A. Alu, M. G. Silveirinha, A. Salandrino, and N. Engheta, “Epsilon-near-zero metamaterials and electromagnetic sources: tailoring the radiation phase pattern,” Phys. Rev. B **75**, 155410–13 (2007). [CrossRef]

**36. **B. Edwards, A. Alu, M. E. Young, M. Silveirinha, and N. Engheta, “Experimental verification of epsilon-near-zero metamaterial coupling and energy squeezing using a microwave waveguide,” Phys. Rev, Lett. **100**, 033903–4 (2008). [CrossRef]

**37. **B. Edwards and N. Engheta, “Experimental verification of displacement-current conduits in metamaterials-inspired optical circuitry,” Phys. Rev. Lett. **108**, 193902–5 (2012). [CrossRef] [PubMed]

**38. **A. Salandrino and N. Engheta, “Far-field subdiffraction optical microscopy using metamaterial crystals: theory and simulations,” Phys. Rev. B **74**, 075103–5 (2006). [CrossRef]

**39. **P. B. Johnson and R. W. Christy, “Optical constants of the noble metals,” Phys. Rev. B **6**, 43704379 (1972). [CrossRef]

**40. **CST Microwave Studio at *www.cst.com*