Abstract
The phase matching condition relating the real transverse wave vectors across a periodic boundary has been generalized to the case of complex transverse wave vectors. Based on this generalization, we describe diffraction of a complex Bloch wave propagating within a composite prism, and show that the resulting light in free space is an inhomogeneous plane wave in the presence of losses within the prism.
©2010 Optical Society of America
1. Introduction
Negative refraction in lossless dielectric photonic crystals [1,2] has been well understood using the equi-frequency contour (EFC) approach, which is based on a phase matching condition across an interface of periodicity d, [3]. Here, is the magnitude of the interface reciprocal latice vector, m is an integer for denoting different diffraction orders, and and are real wave vectors of a Bloch wave and a plane wave, respectively. This condition is a corollary of Bloch's Theorem [3], and connects the transverse components (denoted by subscript t) of and . It has also been used for describing the diffraction of quantum mechanical Bloch waves into free space in angle resolved photoemission spectroscopy experiments [4] that probe the electronic band structure of a solid. In the more general case of lossy periodic media, however, Bloch waves acquire an imaginary component of the wave vector [5]. To our knowledge, diffraction of complex Bloch waves by periodic photonic structures has not been generally treated.
Negative refraction has also been observed in negative index metamaterials (NIMs) [6–8], which are typically periodic structures with subwavelength unit cells containing lossy metal-dielectric composites. Although complex band structure of metamaterials has been routinely calculated, experiments [6–8] that demonstrate negative refraction by composite NIM prisms are generally interpreted as the behavior of a homogeneous NIM with effective permittivity,ε, and permeability, μ.
Recently, we designed a near-infrared NIM with subwavelength unit cells [9]. In that work, assuming zero loss and using the phase matching condition above, we considered diffraction of a real Bloch wave by the periodic hypotenuse of the metamaterial prism, and derived a general relationship between the band structure of a NIM unit cell and the effective index. We applied this relation approximately to the full complex band structure of the NIM unit cell, and achieved reasonable agreement with full-wave simulations of prism refraction [9]. However, substantial differences persist. In addition, it was assumed [9] that the Bloch wave vector is either parallel or anti-parallel to the group velocity, which is not necessarily true in the more general case of unit cells with anisotropy [2]. Moreover, it is not clear whether the neglect of loss is always justified. In fact, the issue of how lossy materials might influence experimental interpretation has been controversial [10] since the first demonstration of negative refraction by NIMs at microwave frequencies [6]. In a subsequent experiment [7], numerical simulations assuming a homogeneous NIM with varying degree of losses have been performed to show independence of the direction of the refracted beam on loss. The effective parameters of the homogeneous NIM [7] have been obtained from the experimental unit cells with the scattering (S) parameter method [11], a homogenization approach that assigns effective ε and μ to inhomogeneous composites. However, these simulations neglect the inherent structure of the unit cell and are therefore somewhat problematic, since the sizes of most NIM unit cells fabricated to date remain a considerable fraction of wavelength (λ). Additionally, the S-parameter method is known to sometimes yield anomalies such as unphysical negative imaginary parts of ε or μ [12].
In this work, we generalize previous treatment of NIM refraction [9] by including the effect of losses and possible anisotropy without using homogenized effective parameters. Based on a generalized phase matching condition with complex transverse wave vectors for periodic media, we describe the diffraction of a complex Bloch wave propagating within a composite prism, and show that the detected light is an inhomogeneous plane wave due to losses in the prism.
2. Diffraction of a complex Bloch wave
Figure 1 shows a schematic of a generalized NIM negative refraction experiment [6–8] that employs bulk prisms made of cubic unit cells of size a. A plane wave is normally incident on the bottom of the prism along and excites Bloch modes within the composite, where is normal to the plane of incidence. We assume a sufficiently large but finite prism, such that multiple reflections within the prism and diffraction by its corners can be neglected. The wave vector of the Bloch mode in the first Brillouin zone must be along , due to conservation of its transverse component. The complex dispersion relation of this Bloch mode, , can be calculated [5]. Here, ω is the (real) angular frequency, , and . This choice of real frequency and complex wave vector is consistent with typical experiments that are conducted in the frequency domain with near monochromatic illumination [7]. The electric and magnetic field of this Bloch mode are and , respectively. Here, , and and are both functions periodic on the cubic lattice. The component of the time averaged power flow of this Bloch mode per unit cell is given by the surface integral of the Poynting vector (proportional to momentum density), i.e.,
where is the cross section of a single unit cell at in the plane. Therefore, , and has no or dependence, due to the periodicity of and . Since power must flow away from the surface, , corresponding to decay of along due to dissipation within the composite, although in general the localPoynting vector and wave vector may not be in the same direction [13]. Similar considerations of and components of the power flow per unit cell show that they are periodic, but do not decay, along their respective axes. These components are not necessarily zero, due to the possible anisotropy of the unit cell [2]. Therefore, we do not assume a total power flow along , although typical unit cells [6–8] are, by design, approximately symmetric with respect to the plane.For a sufficiently large prism, only the complex Bloch mode with the lowest loss (corresponding to the smallest positive ) contributes to wave propagation within the bulk of the composite [9,14,15]. Modes with correspond to an anti-parallel phase velocity and power flow along . The figure of merit (FOM) used for gauging the loss of the unit cell is the inverse of the loss tangent, . Typical low-loss NIM unit cells at microwave frequencies [7] are characterized by a high FOM~100 and a narrow negative index band in the first Brillouin zone. This band has opposite phase and group velocities [16] (i.e. ), which is consistent with the Bloch mode in Fig. 1 with , since the direction of the component of group velocity, , is also along that of the power flow () in the limit of zero loss [17].
Due to the finite size of the unit cell, the hypotenuse of the prism is stepped [18], with prism angle . Here, is the interface periodicity, and there are one and l unit cells along and per step, respectively. The Bloch wave is impingent on the hypotenuse directed along at an apparent incidence angle, θ, and excites reflected Bloch waves and transmitted plane waves. Obtaining transmission and reflection coefficients requires matching of these waves along the stepped hypotenuse with boundary conditions. In previous work on negative refraction in lossless dielectric photonic crystals [1,2], details of the photonic crystal surface have not been considered, and as a result transmission and reflection coefficients are not available [6–8]. Independent of details of a period on the interface, Bloch’s Theorem dictates [3] the phase matching condition
where , and is the complex transverse wave vector of the m th transmitted diffraction order and reflected Bloch wave in the extended Brillouin zone scheme, as shown in Fig. 1. Since Bloch’s Theorem does not require to be real [19], this equation applies even when , i.e. in the presence of loss within the prism. In scalar form, this generalized phase matching condition can be written as:where the subscript y denotes the scalar component along . Both and are along also as a result of this phase matching condition. Since Bloch’s Theorem derives from the discrete translational symmetry of the physical structure, Eq. (3) is independent of the possible nonzero power flow of the Bloch mode transverse to due to anisotropy [2]. Equation (3) is also independent of the types of waves involved and may apply to other similarly periodic media.In the far field, the propagating wave is a solution to the Helmholtz Equation with translational symmetry associated with phase factor, . This eigenmode in a lossless medium must be an inhomogeneous plane wave [20] with orthogonal planes of constant phase and amplitude. Its wave vector has both real and imaginary components along the normal to these two planes, as indicated by and , respectively, in Fig. 1. The inhomogeneous plane wave amplitude diverges at , but in physical systems the wave is always bounded by a finite aperture. Declercq et al. [21] found that a single inhomogeneous plane wave component dominates the behavior of the bounded wave. Indeed, the time-averaged power flow of a single inhomogeneous plane wave is predicted to be along the real wave (ray) vector [20]. Such predictions, including those based on Eq. (3), have been experimentally verified at acoustical frequencies extensively [22–24].
The dispersion relation of an inhomogeneous plane wave is [20]
where k0 = 2π/λ is the vacuum wave number. For m = 0, solutions of to Eqs. (3) and (4) are:In the lossless case, , and Eq. (5) reduces to or . The former solution corresponds to a homogeneous plane wave () when , while the latter corresponds to an evanescent wave propagating along the surface in the presence of total internal reflection when . When , both and are between the two solutions in Eq. (5), and the smaller (negative) solution is discarded, since Eq. (4) requires that . Higher order waves corresponding to m≠0 can also be excited. Given that typically [6] d~λ/2 and, therefore, , Eq. (4) shows for m≠0. The imaginary component of Eq. (3) therefore implies that the corresponding rays will be at near-grazing angles and bound to the interface, as indicated by and in Fig. 1.3. Comparison with NIM refraction experiments
Equation (3) implies that the real part of the effective index follows Snell’s Law, viz:
where is the refraction angle for m = 0. This index corresponds to the far-field power flow direction and is used for explaining experimental results. In the lossless case, Eq. (6) gives:where the sgn function ensures that the group velocity is along , consistent with previous results [9]. Assuming parameters (θ, λ/a, and FOM) typical of experimental unit cells [6,7,9,25,26], we calculate nr as a function of , with results shown in Fig. 2 . FOM = 3 typical of unit cells in the optical domain [9,26] leads to a deviation of nr from Eq. (7) that is dependent on both θ and , particularly toward the edge of the Brillouin zone. These observations are consistent with steady state full-wave simulations of and prisms consisting of unit cells of a near-infrared frequency NIM structure [9].In contrast, the lossless case is nearly identical to the case for FOM = 100, which is often observed at microwave frequencies [7]. This suggests that loss can be neglected for typical microwave NIM structures in refraction experiments. We calculated band structures of the unit cells in Refs. [7,8] and derived the effective indices using Eq. (7). The results are plotted together with the experimental data in Fig. 3 . Band structures [16] calculated for the first reported structure are also consistent with the experimentally measured indices [6] using , and are characterized by a dispersion relation within the first Brillouin zone. This “backward” Bloch band also accounts for negative refraction in lossless dielectric photonic crystals [1].
4. Discussion
The above general treatment explicitly accounts for the inhomogeneity of the unit cell through a complex Bloch wave, but requires no assumptions about ε and μ. Previous experiments [6–8] have generally been interpreted by assuming a homogeneous NIM, whose optical properties are due to a superposition of the negative ε of metal wires [27] and the negative μ of the split ring resonators [28] that comprise the prism. Although both interpretations can describe the experiments, their predictions may differ in some occasions. For example, Eq. (2) suggests that more than one far-field beam may be excited by the Bloch wave when the unit cell size is sufficiently large but still subwavelength, as is the case for some NIM unit cells. In previous NIM refraction experiments [18,25], such phenomena have been observed and explained as the diffraction by the stepped interface grating between air and a homogeneous NIM, which is equivalent to interpreting the Bloch wave vector in Eq. (2) to be that of a plane wave. If the stepped interface is linearized by using partial unit cells, this homogeneous NIM interpretation predicts the transmission of a plane wave through the planar interface between two homogeneous media, where the additional far-field beam disappears. However, the interpretation in terms of Bloch waves suggests that this beam is still present, since Eq. (2) derives from the translational symmetry of the periodic interface only, independent of the details of a period. The observation of two propagating orders from a photonic crystal prism with a linear hypotenuse has been previously reported [29], where the angular positions of the far-field beams are consistent with Eq. (2). This difference is evidence for spatial dispersion, the deviation of the behavior of the NIM from its approximate local homogeneous model.
The results presented here may also suggest the existence of a minimum unit cell size of optical NIMs. The maximum ratio of wavelength to periodicity currently known for optical NIMs is about 7 [9,30]. Tsukerman [31] has recently predicted a minimum unit cell size imposed by a backward Bloch band for photonic crystals consisting of non-dispersive dielectrics. Since negative refraction by metamaterials is shown to be governed by the same Bloch band in the zero-loss limit, whether the minimum periodicity for NIMs containing dispersive metals at optical frequencies is restricted by a similar fundamental limit (apart from technological constraints) remains an open question.
5. Conclusion
In conclusion, based on a generalized phase matching condition, we have shown that the negative refractive behavior of lossless dielectric photonic crystals and lossy metal-dielectric periodic NIMs can be given a unified explanation. We have also discussed its implications for the local homogeneous model of NIMs as well as the possible existence of a minimum unit cell size of optical NIMs.
Acknowledgements
The authors thank the AFOSR MURI contract No. FA 9550-06-01-0279 for financial support. We also thank Prof. Paul Berman for reading the manuscript and useful discussions.
References and links
1. M. Notomi, “Theory of light propagation in strongly modulated photonic crystals: Refraction like behavior in the vicinity of the photonic band gap,” Phys. Rev. B 62(16), 10696–10705 (2000). [CrossRef]
2. C. Luo, S. G. Johnson, J. D. Joannopoulos, and J. B. Pendry, “All-angle negative refraction without negative effective index,” Phys. Rev. B 65(20), 201104 (2002). [CrossRef]
3. J. D. Joannopoulos, R. D. Meade, and J. N. Winn, Photonic Crystals: Molding the Flow of Light (Princeton University Press, Princeton and Oxford, 2008), pp. 221–222.
4. A. Damascelli, Z. Hussain, and Z. X. Shen, “Angle-resolved photoemission studies of the cuprate superconductors,” Rev. Mod. Phys. 75(2), 473–541 (2003). [CrossRef]
5. J. B. Pendry, “Photonic Band Structures,” J. Mod. Opt. 41(2), 209–229 (1994). [CrossRef]
6. R. A. Shelby, D. R. Smith, and S. Schultz, “Experimental verification of a negative index of refraction,” Science 292(5514), 77–79 (2001). [CrossRef] [PubMed]
7. C. G. Parazzoli, R. B. Greegor, K. Li, B. E. C. Koltenbah, and M. Tanielian, “Experimental verification and simulation of negative index of refraction using Snell’s law,” Phys. Rev. Lett. 90(10), 107401 (2003). [CrossRef] [PubMed]
8. R. B. Greegor, C. G. Parazzoli, K. Li, B. E. C. Koltenbah, and M. Tanielian, “Experimental determination and numerical simulation of the properties of negative index of refraction materials,” Opt. Express 11(7), 688–695 (2003). [CrossRef] [PubMed]
9. X. H. Zhang, M. Davanço, Y. Urzhumov, G. Shvets, and S. R. Forrest, “From scattering parameters to Snell’s law: a subwavelength near-infrared negative-index metamaterial,” Phys. Rev. Lett. 101(26), 267401 (2008). [CrossRef] [PubMed]
10. N. Garcia and M. Nieto-Vesperinas, “Is there an experimental verification of a negative index of refraction yet?” Opt. Lett. 27(11), 885–887 (2002). [CrossRef]
11. D. R. Smith, S. Schultz, P. Markos, and C. M. Soukoulis, “Determination of effective permittivity and permeability of metamaterials from reflection and transmission coefficients,” Phys. Rev. B 65(19), 195104 (2002). [CrossRef]
12. T. Koschny, P. Markos, D. R. Smith, and C. M. Soukoulis, “Resonant and antiresonant frequency dependence of the effective parameters of metamaterials,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 68(6 Pt 2), 065602 (2003). [CrossRef]
13. D. de Ceglia, M. A. Vincenti, M. G. Cappeddu, M. Centini, N. Akozbek, A. D'Orazio, J. W. Haus, M. J. Bloemer, and M. Scalora, “Tailoring metallodielectric structures for superresolution and superguiding applications in the visible and near-ir ranges,” Phys. Rev. A 77(3), 033848 (2008). [CrossRef]
14. D. Seetharamdoo, R. Sauleau, K. Mahdjoubi, and A. C. Tarot, “Effective parameters of resonant negative refractive index metamaterials: Interpretation and validity,” J. Appl. Phys. 98(6), 063505 (2005). [CrossRef]
15. C. Rockstuhl, C. Menzel, T. Paul, T. Pertsch, and F. Lederer, “Light propagation in a fishnet metamaterial,” Phys. Rev. B 78(15), 155102 (2008). [CrossRef]
16. R. A. Shelby, D. R. Smith, S. C. Nemat-Nasser, and S. Schultz, “Microwave transmission through a two-dimensional, isotropic, left-handed metamaterial,” Appl. Phys. Lett. 78(4), 489–491 (2001). [CrossRef]
17. K. Sakoda, Optical Properties of Photonic Crystals (Springer, Berlin, 2001), pp. 30–32.
18. D. R. Smith, P. M. Rye, J. J. Mock, D. C. Vier, and A. F. Starr, “Enhanced diffraction from a grating on the surface of a negative-index metamaterial,” Phys. Rev. Lett. 93(13), 137405 (2004). [CrossRef] [PubMed]
19. N. W. Ashcroft, and N. D. Mermin, Solid State Physics, Solid State Physics (Holt, Rinehart and Winston, New York, 1976), pp. 368–369.
20. P. C. Clemmow, The Plane Wave Spectrum Representation of Electromagnetic Fields (Pergamon Press, Oxford and New York, 1966), pp. 13–14.
21. N. F. Declercq, J. Degrieck, and O. Leroy, “The Laplace transform to describe bounded inhomogeneous waves,” J. Acoust. Soc. Am. 116(1), 51–60 (2004). [CrossRef]
22. W. Huang, R. Briers, S. I. Rokhlin, and O. Leroy, “Experimental-Study of Inhomogeneous Wave Reflection from a Solid-Air Periodically Rough Boundary Using Leaky Rayleigh-Waves,” J. Acoust. Soc. Am. 96(1), 363–369 (1994). [CrossRef]
23. R. Briers, O. Leroy, O. Poncelet, and M. Deschamps, “Experimental verification of the calculated diffraction field generated by inhomogeneous waves obliquely incident on a periodically rough liquid-solid boundary,” J. Acoust. Soc. Am. 106(2), 682–687 (1999). [CrossRef]
24. N. F. Declercq, R. Briers, J. Degrieck, and O. Leroy, “The history and properties of ultrasonic inhomogeneous waves,” IEEE Trans. Ultrason. Ferroelectr. Freq. Control 52(5), 776–791 (2005). [CrossRef] [PubMed]
25. A. A. Houck, J. B. Brock, and I. L. Chuang, “Experimental observations of a left-handed material that obeys Snell’s law,” Phys. Rev. Lett. 90(13), 137401 (2003). [CrossRef] [PubMed]
26. 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(7211), 376–379 (2008). [CrossRef] [PubMed]
27. J. B. Pendry, A. J. Holden, W. J. Stewart, and I. Youngs I, “Extremely low frequency plasmons in metallic mesostructures,” Phys. Rev. Lett. 76(25), 4773–4776 (1996). [CrossRef] [PubMed]
28. J. B. Pendry, A. J. Holden, D. J. Robbins, and W. J. Stewart, “Magnetism from conductors and enhanced nonlinear phenomena,” IEEE Trans. Microw. Theory Tech. 47(11), 2075–2084 (1999). [CrossRef]
29. S. Foteinopoulou and C. M. Soukoulis, “Negative refraction and left-handed behavior in two-dimensional photonic crystals,” Phys. Rev. B 67(23), 235107 (2003). [CrossRef]
30. X. H. Zhang, M. Davanco, Y. Urzhumov, G. Shvets, and S. R. Forrest, “A Subwavelength Near-Infrared Negative-Index Material,” Appl. Phys. Lett. 94(13), 131107 (2009). [CrossRef]
31. I. Tsukerman, “Negative refraction and the minimum lattice cell size,” J. Opt. Soc. Am. B 25(6), 927–936 (2008). [CrossRef]