Calculation of effective permittivity, permeability, and surface impedance of negative-refraction photonic crystals

doi:10.1364/OE.15.008340

« Show journal navigation

Calculation of effective permittivity, permeability, and surface impedance of negative-refraction photonic crystals

Zhaolin Lu and Dennis W. Prather »View Author Affiliations

Optics Express, Vol. 15, Issue 13, pp. 8340-8345 (2007)
http://dx.doi.org/10.1364/OE.15.008340

View Full Text Article

Acrobat PDF (374 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. Calculation of effective permittivity and permeability
3. Calculation of effective surface impedance
4. Summary
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Multimedia
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (12)
Cited By
Figures (4)
Multimedia (2)
Metrics

Abstract

We consider the eigen-fields of a two-dimensional negative-refraction photonic crystal and obtain negative effective permittivity and negative effective permeability. Effective permittivity, permeability, and surface impedance are calculated by averaging the eigen-fields. The value of the surface impedance is shown to be location-dependent and is validated by finite-difference time-domain simulations. The unique power propagation mechanism in the photonic crystal is demonstrated through time-evolution of eigen-fields.

1. Introduction

Recently, the work on metamaterial [1–4

J. B. Pendry, A. J. Holden, W. J. Stewart, and I. Youngs, “Extremely low frequency plasmons in metallic mesostructures,” Phys. Rev. Lett. 76, 4773(1996). [CrossRef] [PubMed]

] and “perfect lenses” [5

J. B. Pendry, “Negative refraction makes a perfect lens,” Phys. Rev. Lett. 85, 3966 (2000). [CrossRef] [PubMed]

] revived Veselago’s lefthanded materials [6

V. G. Veselago, “The electrodynamics of substances with simultaneously negative values of permittivity and permeability,” Sov. Phys. Usp. 10, 509 (1968). [CrossRef]

] and triggered intensive discussion on negative refraction. Meanwhile, negative refraction [7

E. Cubukcu, K. Aydin, E. Ozbay, S. Foteinopoulou, and C. M. Soukoulis, “Electromagnetic wave: negative refraction by photonic crystals,” Nature 423, 604 (2003). [CrossRef] [PubMed]

] and imaging by negative refraction [8

P. V. Parimi, W. T. Lu, P. Vodo, J. Sokoloff, J. S. Derov, and S. Sridhar, “Negative refraction and left-handed electromagnetism in microwave photonic crystals,” Phys. Rev. Lett. 92, 127401 (2004). [CrossRef] [PubMed]

Z. Lu, C. Chen, C. A. Schuetz, S. Shi, J. A. Murakowski, G. J. Schneider, and D. W. Prather, “Sub-wavelength imaging by a flat cylindrical lens using optimized negative refraction,” Appl. Phys. Lett. 87, 091907 (2005). [CrossRef]

] were also demonstrated using 2D photonic crystals (PhCs) by dispersion-engineering. Subsequently, we experimentally achieved subwavelength imaging by a 3D PhC flat lens exhibiting full 3D negative refraction [10

Z. Lu, J. A. Murakowski, C. A. Schuetz, S. Shi, G. J. Schneider, and D. W. Prather, “Three-dimensional subwavelength imaging by a photonic-crystal flat lens using negative refraction at microwave frequencies,” Phys. Rev. Lett. 95, 153901 (2005). [CrossRef] [PubMed]

]. Although significant experimental progress has been made in observing and employing this novel phenomenon, several basic questions are still unclear — for example, what is the relation between metamaterial and negative-refraction PhCs, and specifically whether negative-refraction PhCs have negative effective permittivity and negative permeability. More recently, Decoopman, et al., proposed to calculate the impedance of a negative-refraction PhC by observing the reflection into its matched medium, and numerically demonstrated that the value of the impedance of a PhC strongly depends on its boundary — especially the way of truncation of the elementary cells located on this boundary [11

T. Decoopman, G. Tayeb, S. Enoch, D. Maystre, and B. Gralak, “Photonic Crystal Lens: From Negative Refraction and Negative Index to Negative Permittivity and Permeability,” Phys. Rev. Lett. 97, 073905 (2006). [CrossRef] [PubMed]

]. In their work, they assumed that the relations among effective permittivity, permeability, impedance, and refractive index,

η_{eff} = \sqrt{\frac{μ_{eff}}{ε_{eff}}}

and

n_{eff} = \sqrt{μ_{eff}} \sqrt{ε_{eff}}

, are still valid. Their work provides a very novel approach to calculate effective impedance, but several issues are still open to discuss. First, effective permittivity and permeability as bulk parameters for material should not depend on the boundary. If they only describe the surface properties of the PhC, then the effective parameters inside the PhC are still unknown. Second, negative-refraction PhCs work by the unique photonic dispersion, which however does not provide any clue for negative permittivity and negative permeability. In fact, these two parameters in a PhC are known and accurately given as periodic functions, i.e. μ_r(r) = μ_r(r + R _lmn and ε_r = ε_r (r + R), where R _lmn represents a unit lattice translation vector. Moreover, in their work the determination of the effective parameters relies on numerically simulated results and each simulation provides the parameters for only one boundary. To address these issues, we present a method based on averaging eigen-fields to calculate the effective permittivity, permeability, and surface impedance on all possible boundaries of negative-refraction PhCs.

2. Calculation of effective permittivity and permeability

We consider a 2D triangular-lattice PhC formed by patterning air holes in a low loss material with dielectric constant ε_r = 20. The lattice constant is a and the diameter of the air holes is 2r = 0.7a, as shown in Fig. 1(a). Subwavelength imaging by this PhC was reported in our recent work [9

]. For transverse-magnetic-like (TM-like, electric field is polarized along the z-axis) modes, the electromagnetic fields have three independent components, E_z , H_x , and H_y . The coordinate system as illustrated in Fig. 1(a) will be used throughout this paper. Using μ_r(r) = 1 in Maxwell equations, we obtain the governing function for E_z

\frac{1}{ε_{r}} (\frac{\partial^{2} E_{z}}{{\partial x}^{2}} + \frac{\partial^{2} E_{z}}{{\partial y}^{2}}) + {(\frac{ω}{c})}^{2} E_{z} = 0 .

(1)

According to Bloch’s Theorem,

E_{z} (k, r, t) = \exp (j k ∙ r - jωt) E_{z}^{(k)} (r),

(2)

where E ^(k) _z(r) = E ^(k) _z(r + R _lmn) is a function with the same periodicity as the PhC. Plugging ε_r(r) = ε_r (r + R _lmn) into the governing function and using the plane-wave method [12

S. G. Johnson and J. D. Joannopoulos, “Block-iterative frequency-dormain methods for Maxwell’s equations in a planewave basis,” Opt. Express 8, 173 (2001). [CrossRef] [PubMed]

], we calculate the corresponding frequencies for a series of wave vectors and depict them as the dispersion surfaces [Fig. 1(b)], where the wave vectors are normalized to k _a/2π, and the frequency to ωa/2πc (c is the speed of light in vacuum). To avoid clutter, we only depict the first two bands and the light cone. As can be seen from Fig. 1(b), the second band is curved downward. As a result, the group velocity, v _g = ∇_k ω(k), is opposite to the phase velocity and negative refraction ensues. In this work, we focused on the intersection of the second band and the light cone because it was suggested that the corresponding effective index n _eff = -1. In order to see that region clearly, we project the second band in the vicinity of the intersection into a series of equi-frequency contours. As shown in Fig. 1(c), these contours are inward-growing and centered at the origin (k_x = 0, k_y = 0). In particular, we can see that the intersection occurs at ω = 0.236. Since the contour of ω = 0.236 is approximately circular, we can treat it as an isotropic homogenous medium. Using the plane wave method, we can also calculate the eigen-field, E ^(k) _z(r), for a given wave vector. Note that the Bloch wave, E_z (k,r,t) = exp(j k ∙ r - jωt)E(k)_z, can be understood as the eigen-field E(k)_z(r) packaged in a plane wave, exp(j k ∙ r - jωt). The corresponding eigen-fields H_x and H_y can be calculated through

{\begin{matrix} H_{x} (k, r, t) = \frac{1}{jωμ} \frac{{\partial E}_{z} (k, r, t)}{\partial y} \\ H_{y} (k, r, t) = \frac{1}{- jωμ} \frac{{\partial E}_{z} (k, r, t)}{\partial x} \end{matrix}

(3)

Fig. 1. (a) The illustration of the PhC and coordinate system for the study. (b) The photonic dispersion diagram for the PhC and the light in this frequency range. (c) The equi-frequency dispersion contours for the second band and the light cone.

For simplicity, we consider the case where light propagates along the x-axis. We first tried (k_x = 0.236, k_y = 0) and found that the power propagates opposite to the x-axis, which is predicted by the dispersion of the PhC. So we switched the wave vector to (k_x = -0.236, k_y = 0) and repeated the computation. Figures 2 (a, b) show the amplitude and phase distributions of the corresponding eigen-field, E(k)_z(r), respectively: the amplitude distribution has the same periodicity as ε_r(r) = ε_r (r + R _lmn) and the phase of E(k)_z(r) increases along the x-axis. For comparison, the plane wave, exp(j k ∙ r - jωt), has longer period along the x-axis than that of the PhC, and its phase decreases along the x-axis. See Fig. 2(c). In other words, the phase growth of the package and the carrier waves is along opposite directions, i.e. the carrier wave propagates backward while the package wave propagates forward. The phase distribution of the package wave determines the power flow direction because S(av)_x(r) = -0.5Re{E_z (k,r,t)[H_y (k,r,t)]^*} = -0.5Re{E^(k) _z (r)[H^(k) _y (r]^*}. As a result, the power flows forward. Compared to H ^(k) _y(r),H ^(k) _x(r) is very small and can be neglected. When they evolve with time, the Bloch wave moves forwards, while the carrier wave moves backwards. The supplemental movies illustrate Re{exp(j k ∙ r - jωt)E ^(k) _z(r)} and Re{exp(j k ∙ r - jωt)}, respectively. Different phase growth directions of the carrier and package explain why the group velocity is opposite to the phase velocity.

To obtain the effective permittivity and effective permeability, we need to treat the PhC as a homogenous medium. In a homogenous medium, the solution of Maxwell equations for a single wave vector is plane waves. Therefore, we need to approximate the Bloch waves as plane waves, or equivalently approximate the package waves as constants. A direct approach is to take the average value of the package.

{\begin{matrix} E_{z}^{(2 Dav)} = \frac{\int \int_{Unit Cell} E_{z}^{(k)} (r) dxdy}{Area} \\ H_{y}^{(2 Dav)} = \frac{\int \int_{Unit Cell} H_{y}^{(k)} (r) dxdy}{Area} \end{matrix},

(4)

where the averaging region is chosen as unit cells. In our case, we chose the hexagonal cell located in - √3a/4 < x < √3a/4 and -0.5a < y < 0.5a. To avoid positive and negative values cancelling each other, like the case in sin-function, as well as to keep the “handness” of the electromagnetic wave, the points r are grouped into region I with Re{H ^(k) _y(r)} > 0 and region II with Re{H ^(k) _y(r)} < 0, and the averaging is applied in these two regions separately. In region I, we obtained E ^(2Dav) _z = 0.0202 + 3×10^-5, j, H ^(2Dav) _y = 0.0861 + 6×10^-5; in region II, we obtained E ^(2Dav) _z = -0.0241 - 8×10^-5, H ^(2Dav) _y = -0.0822 - 2×10^-4. In both regions, the imaginary part of average values is negligible. If the electric and magnetic fields of a plane wave are both known, we can in turn calculate the effective parameters

{\begin{matrix} μ_{eff} = \frac{1}{jω H_{y}} \frac{{\partial E}_{z}}{\partial x} = \frac{k}{ω} \frac{E_{z}^{(2 Dav)}}{H_{y}^{(2 Dav)}} \\ ε_{eff} = \frac{1}{jω E_{z}} \frac{{\partial H}_{y}}{\partial x} = \frac{k}{ω} \frac{H_{y}^{(2 Dav)}}{E_{z}^{(2 Dav)}} \end{matrix} .

(5)

In region I, we obtain ε^(I) _eff = - 4.270 and μ^(I) _eff = -0.235. In region II, we obtain ε^(II) _eff = -3.405 and μ^(II) _eff = -0.294. The final effective parameters can be estimated to be

ε_{eff} = \sqrt{ε_{eff}^{(I)}} \sqrt{ε_{eff}^{II}} = - 3.813

and

μ_{eff} = \sqrt{μ_{eff}^{(I)}} \sqrt{μ_{eff}^{II}} = - 0.263

(If we divide the region by Re{E ^(k) _z(r)} > 0 and Re{E ^(k) _z(r)} < 0, very close result can be obtained because the signs of electric and magnetic fields ensure S ^(av) _x(r) > 0 at every point). Consequently, the PhC simultaneously exhibits negative effective permittivity and negative effective permeability. Note that

n_{eff} = \sqrt{ε_{eff}} \sqrt{μ_{eff}} = - 1

is valid and ε_eff μ_eff = (k/ω)² is always satisfied. In addition, these effective parameters are derived from the eigen-fields using effective medium theory [2

J. B. Pendry, A. J. Holden, D. J. Robbins, and W. J. Stewart, “Magnetism from conductors and enhanced nonlinear phenomena,” IEEE Trans. Microwave Theor Tech 47, 2075 (1999). [CrossRef]

]. In this approach,

n_{eff} = \sqrt{ε_{eff}} \sqrt{μ_{eff}} = - 1

is a derived result, rather than a presumption as used in Ref. [11

Fig. 2. (a) The amplitude distribution of the eigen-field, E ^(k) _z(r). (b) The phase distribution of the eigen-field,E ^(k) _z(r). [Media 1] (c) The phase distribution of the plane wave, exp(j k ∙ r - jωt). [Media 2] (d) The impedance with respect to the location of PhC surface.

3. Calculation of effective surface impedance

To obtain the effective impedance, we average the field distribution on the surface according to

{\begin{matrix} E_{z}^{(1 Dav)} = \frac{1}{L_{y}} \int_{Unit Cell} E_{z}^{(k)} (r) dy \\ H_{y}^{(1 Dav)} = \frac{1}{L_{y}} \int_{Unit Cell} H_{y}^{(k)} (r) dy \end{matrix} .

(6)

Note that in the spatial frequency domain,

E_{z}^{(k)} (x_{0}; k_{y}) = \frac{1}{L_{y}} \int_{Unit Cell} E_{z}^{(k)} (x_{0}, y) \exp (- j k_{y} y) dy

and

E_{z}^{(k)} (x_{0}; k_{y} = 0) = \frac{1}{L_{y}} \int_{Unit Cell} E_{z}^{(k)} (x_{0}, y) dy = E_{z}^{(1 Dav)}

. The convolution with the Fourier transform of the plane wave, exp(j k ∙ r), will shift spatial frequency to k. As a result, the averaging process results in the component with the spatial frequency matching that of the plane wave the best. In this fashion, the effective impedance can be calculated using

η_{eff} = - \frac{E_{z}^{(1 Dav)}}{H_{z}^{(1 Dav)}} .

Fig. 3. The simulations (amplitudes) of propagation of light from (a, c) the matched medium and (b, d) the air into the PhC. In (a, b), the PhC surface is located at x = 0; in (c, d), the PhC surface is located at x = √3a/4.

We shifted the termination surface along the x-axis and obtained different values, which mean the impedance is location-dependent. So we refer to the effective impedance as effective surface impedance. The imaginary part of surface impedance is negligible. As can be seen from Fig. 2(d), the surface has the maximum impedance when it goes through the center of holes (x = 0), while the surface has the minimum impedance when it is located at the half way between two neighboring hole columns (x = √3a/4). The effective impedance is not constant, so

η_{eff} = \sqrt{\frac{μ_{eff}}{ε_{eff}}}

is not valid. The reason lies in: 1D-averaging is taken to calculate the surface impedance, while 2D-averaging is used to calculate the effective permittivity and permeability.

To validate the results, we apply the effective impedance into 2D finite-difference time-domain simulations. In each simulation, a collimated light beam (with width 6λ) is normally incident into the PhC from a matched material and the mesh size is 0.004λ (The smaller the mesh size is, the more accurate result can be obtained). The PhC has a specific boundary and the matched material has the corresponding impedance calculated from field-averaging. The edges of the simulated region are paved by perfectly-matched layers (PML). Good transmission requires impedance match, but the refractive index of matched material can be arbitrary — or “ ε_m = ±μ_eff and μ_m = ±μ_eff” are sufficient but not necessary. There are many possibilities to match the surface impedance of the photonic crystal only if

\sqrt{\frac{μ_{m}}{ε_{m}}} = η_{PhC}

, but for convenience we assume the permeability of the matched medium, μ_m = 1. The selection of permeability is different from that presented in Ref. [12

S. G. Johnson and J. D. Joannopoulos, “Block-iterative frequency-dormain methods for Maxwell’s equations in a planewave basis,” Opt. Express 8, 173 (2001). [CrossRef] [PubMed]

], where the constraint,

\sqrt{ε_{m}} \sqrt{μ_{m}} = - 1

, was applied.

In particular, when the boundary goes through the center of holes (x = 0), we apply the calculated effective impedance 0.283. As can be seen in Fig. 3(a), very high transmission is achieved. The back-refraction is less than 0.1%. In contrast, the back-reflection increases to 31.2% if the matched medium is switched to the air, as shown in Fig. 3(b). When the boundary goes along the middle line between two hole columns (x = √3a/4), we apply the calculated effective impedance 0.221. As can be seen in Fig. 3(c), very high transmission is achieved. The back-refraction is less than 0.1%. In contrast, the back-reflection increases to 40.7% if the matched medium is switched to the air, as shown in Fig. 3(d).

In addition, different wave vectors have different eigen-fields and these effective parameters depend on eigen-fields. As a result, these parameters are also different when light is incident at different angles. We calculated the eigen-fields for all wave vectors along the equi-frequency dispersion contour corresponding to ω = 0.236 as shown in Fig. 1(c). Figures 4(a,b) show the amplitude and phase distributions, respectively, when wave vector is at 45° to the x-axis. The phase distribution of E ^(k) _z(r) is also opposite to that of the carrier as shown in Fig. 4(c). The surface impedance strongly depends on incident angles. Figure 4(d) plots the effective impedance for the truncation surfaces located at x = 0 and x = √3a/4 changing with respect to incident angle, respectively.

Fig. 4. (a) The amplitude distribution of the eigen-field, E ^(k) _z(r)(ω_n = 0.236 and θ = 45° ). (b) The phase distribution of the eigen-field, E ^(k) _z(r)(ω_n = 0.236 and θ = 45° ). (c) The phase distribution of the plane wave, exp(j k ∙ r - jωt). (d) The impedance with respect to the location of PhC surface (ω_n = 0.236 and θ = 45°).

4. Summary

To summarize, a plane wave method was employed to solve the eigen-fields of a negative refraction PhC. By averaging the eigen-fields, negative effective permittivity and permeability were obtained and calculated. Effective surface impedance was shown to change with truncation surfaces and verified by the finite-difference time-domain simulations. The power propagation mechanism was also investigated by observing the time-evolution of eigen-fields. The application of this method can be extended to analyze the impedance of PhCs for other applications, such as self-collimation, super-prism, photonic bandgap waveguide, and slow light.

References and links

1.	J. B. Pendry, A. J. Holden, W. J. Stewart, and I. Youngs, “Extremely low frequency plasmons in metallic mesostructures,” Phys. Rev. Lett. 76, 4773(1996). [CrossRef] [PubMed]
2.	J. B. Pendry, A. J. Holden, D. J. Robbins, and W. J. Stewart, “Magnetism from conductors and enhanced nonlinear phenomena,” IEEE Trans. Microwave Theor Tech 47, 2075 (1999). [CrossRef]
3.	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 (2000). [CrossRef] [PubMed]
4.	R. A. Shelby, D. R. Smith, and S. Schultz, “Experimental verification of a negative index of refraction,” Science 292, 77 (2001). [CrossRef] [PubMed]
5.	J. B. Pendry, “Negative refraction makes a perfect lens,” Phys. Rev. Lett. 85, 3966 (2000). [CrossRef] [PubMed]
6.	V. G. Veselago, “The electrodynamics of substances with simultaneously negative values of permittivity and permeability,” Sov. Phys. Usp. 10, 509 (1968). [CrossRef]
7.	E. Cubukcu, K. Aydin, E. Ozbay, S. Foteinopoulou, and C. M. Soukoulis, “Electromagnetic wave: negative refraction by photonic crystals,” Nature 423, 604 (2003). [CrossRef] [PubMed]
8.	P. V. Parimi, W. T. Lu, P. Vodo, J. Sokoloff, J. S. Derov, and S. Sridhar, “Negative refraction and left-handed electromagnetism in microwave photonic crystals,” Phys. Rev. Lett. 92, 127401 (2004). [CrossRef] [PubMed]
9.	Z. Lu, C. Chen, C. A. Schuetz, S. Shi, J. A. Murakowski, G. J. Schneider, and D. W. Prather, “Sub-wavelength imaging by a flat cylindrical lens using optimized negative refraction,” Appl. Phys. Lett. 87, 091907 (2005). [CrossRef]
10.	Z. Lu, J. A. Murakowski, C. A. Schuetz, S. Shi, G. J. Schneider, and D. W. Prather, “Three-dimensional subwavelength imaging by a photonic-crystal flat lens using negative refraction at microwave frequencies,” Phys. Rev. Lett. 95, 153901 (2005). [CrossRef] [PubMed]
11.	T. Decoopman, G. Tayeb, S. Enoch, D. Maystre, and B. Gralak, “Photonic Crystal Lens: From Negative Refraction and Negative Index to Negative Permittivity and Permeability,” Phys. Rev. Lett. 97, 073905 (2006). [CrossRef] [PubMed]
12.	S. G. Johnson and J. D. Joannopoulos, “Block-iterative frequency-dormain methods for Maxwell’s equations in a planewave basis,” Opt. Express 8, 173 (2001). [CrossRef] [PubMed]

OCIS Codes
(130.3130) Integrated optics : Integrated optics materials
(160.4760) Materials : Optical properties

ToC Category:
Photonic Crystals

History
Original Manuscript: March 2, 2007
Revised Manuscript: May 18, 2007
Manuscript Accepted: May 21, 2007
Published: June 19, 2007

Citation
Zhaolin Lu and Dennis W. Prather, "Calculation of effective permittivity, permeability, and surface impedance of negative-refraction photonic crystals," Opt. Express 15, 8340-8345 (2007)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-15-13-8340

Sort: Year | Journal | Reset

References

J. B. Pendry, A. J. Holden, W. J. Stewart, and I. Youngs, "Extremely low frequency plasmons in metallic mesostructures," Phys. Rev. Lett. 76, 4773 (1996). [CrossRef] [PubMed]
J. B. Pendry, A. J. Holden, D. J. Robbins, and W. J. Stewart, "Magnetism from conductors and enhanced nonlinear phenomena," IEEE Trans. Microwave Theory Tech. 47, 2075 (1999). [CrossRef]
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 (2000). [CrossRef] [PubMed]
R. A. Shelby, D. R. Smith, and S. Schultz, "Experimental verification of a negative index of refraction," Science 292, 77 (2001). [CrossRef] [PubMed]
J. B. Pendry, "Negative refraction makes a perfect lens," Phys. Rev. Lett. 85, 3966 (2000). [CrossRef] [PubMed]
V. G. Veselago, "The electrodynamics of substances with simultaneously negative values of permittivity and permeability," Sov. Phys. Usp. 10, 509 (1968). [CrossRef]
E. Cubukcu, K. Aydin, E. Ozbay, S. Foteinopoulou, and C. M. Soukoulis, "Electromagnetic wave: negative refraction by photonic crystals," Nature 423, 604 (2003). [CrossRef] [PubMed]
P. V. Parimi, W. T. Lu, P. Vodo, J. Sokoloff, J. S. Derov, and S. Sridhar, "Negative refraction and left-handed electromagnetism in microwave photonic crystals," Phys. Rev. Lett. 92, 127401 (2004). [CrossRef] [PubMed]
Z. Lu, C. Chen, C. A. Schuetz, S. Shi, J. A. Murakowski, G. J. Schneider, and D. W. Prather, "Sub-wavelength imaging by a flat cylindrical lens using optimized negative refraction," Appl. Phys. Lett. 87, 091907 (2005). [CrossRef]
Z. Lu, J. A. Murakowski, C. A. Schuetz, S. Shi, G. J. Schneider, and D. W. Prather, "Three-dimensional subwavelength imaging by a photonic-crystal flat lens using negative refraction at microwave frequencies," Phys. Rev. Lett. 95, 153901 (2005). [CrossRef] [PubMed]
T. Decoopman, G. Tayeb, S. Enoch, D. Maystre, and B. Gralak, "Photonic Crystal Lens: from Negative Refraction and Negative Index to Negative Permittivity and Permeability," Phys. Rev. Lett. 97, 073905 (2006). [CrossRef] [PubMed]
S. G. Johnson and J. D. Joannopoulos, "Block-iterative frequency-dormain methods for Maxwell’s equations in a planewave basis," Opt. Express 8, 173 (2001). [CrossRef] [PubMed]

Appl. Phys. Lett.

Z. Lu, C. Chen, C. A. Schuetz, S. Shi, J. A. Murakowski, G. J. Schneider, and D. W. Prather, "Sub-wavelength imaging by a flat cylindrical lens using optimized negative refraction," Appl. Phys. Lett. 87, 091907 (2005). [CrossRef]

IEEE Trans. Microwave Theory Tech.

J. B. Pendry, A. J. Holden, D. J. Robbins, and W. J. Stewart, "Magnetism from conductors and enhanced nonlinear phenomena," IEEE Trans. Microwave Theory Tech. 47, 2075 (1999). [CrossRef]

Nature

E. Cubukcu, K. Aydin, E. Ozbay, S. Foteinopoulou, and C. M. Soukoulis, "Electromagnetic wave: negative refraction by photonic crystals," Nature 423, 604 (2003). [CrossRef] [PubMed]

Opt. Express

S. G. Johnson and J. D. Joannopoulos, "Block-iterative frequency-dormain methods for Maxwell’s equations in a planewave basis," Opt. Express 8, 173 (2001). [CrossRef] [PubMed]

Phys. Rev. Lett.

J. B. Pendry, A. J. Holden, W. J. Stewart, and I. Youngs, "Extremely low frequency plasmons in metallic mesostructures," Phys. Rev. Lett. 76, 4773 (1996). [CrossRef] [PubMed]
P. V. Parimi, W. T. Lu, P. Vodo, J. Sokoloff, J. S. Derov, and S. Sridhar, "Negative refraction and left-handed electromagnetism in microwave photonic crystals," Phys. Rev. Lett. 92, 127401 (2004). [CrossRef] [PubMed]
Z. Lu, J. A. Murakowski, C. A. Schuetz, S. Shi, G. J. Schneider, and D. W. Prather, "Three-dimensional subwavelength imaging by a photonic-crystal flat lens using negative refraction at microwave frequencies," Phys. Rev. Lett. 95, 153901 (2005). [CrossRef] [PubMed]
T. Decoopman, G. Tayeb, S. Enoch, D. Maystre, and B. Gralak, "Photonic Crystal Lens: from Negative Refraction and Negative Index to Negative Permittivity and Permeability," Phys. Rev. Lett. 97, 073905 (2006). [CrossRef] [PubMed]
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 (2000). [CrossRef] [PubMed]
J. B. Pendry, "Negative refraction makes a perfect lens," Phys. Rev. Lett. 85, 3966 (2000). [CrossRef] [PubMed]

Science

R. A. Shelby, D. R. Smith, and S. Schultz, "Experimental verification of a negative index of refraction," Science 292, 77 (2001). [CrossRef] [PubMed]

Sov. Phys. Usp.

V. G. Veselago, "The electrodynamics of substances with simultaneously negative values of permittivity and permeability," Sov. Phys. Usp. 10, 509 (1968). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper