Abstract

Complex and interesting electromagnetic behavior can be found in spaces with non-flat topology. When considering the properties of an electromagnetic medium under an arbitrary coordinate transformation an alternative interpretation presents itself. The transformed material property tensors may be interpreted as a different set of material properties in a flat, Cartesian space. We describe the calculation of these material properties for coordinate transformations that describe spaces with spherical or cylindrical holes in them. The resulting material properties can then implement invisibility cloaks in flat space. We also describe a method for performing geometric ray tracing in these materials which are both inhomogeneous and anisotropic in their electric permittivity and magnetic permeability.

©2006 Optical Society of America

1. Introduction

Recently the use of coordinate transformations to produce material specifications that control electromagnetic fields in interesting and useful ways has been discussed. We have described such a method in which the transformation properties of Maxwell’s equations and the constitutive relations can yield material descriptions that implement surprising functionality, such as invisibility [1]. Another author described a similar method where the two dimensional Helmholtz equation is transformed to produce similar effects in the geometric limit [2]. We note that these theoretical design methods are of more than academic interest, as the material specifications can be implemented with metamaterial technology [3, 4, 5, 6, 7]. An experimental demonstration of invisibility cloaking using metamaterials has recently been published [8]. There has also been recent work on invisibility cloaking that does not employ the transformation method [9, 10].

In this article we describe how to calculate these material properties directly as Cartesian tensors, though the spherical, and cylindrical cloaks have been previously analyzed in their corresponding coordinate systems [1, 11]. Here we choose to work in Cartesian coordinates for two reasons. First, Cartesian tensors are very straightforward, dimensionally homogenous, and relatively easily applied to any geometry, not just the spherical and cylindrical geometries used as examples here. Second, integration of the ray equations in transformation derived media, which has not yet been described in the literature, is particularly simple to perform in Cartesian coordinates. Obviously, the rays or fields in a transformation medium may be found simply by applying the transform to rays or fields in the simple un-transformed medium, but we believe that integrating the ray equations is an important validation step in the early development of the theory.

Advances in electromagnetic metamaterial technology promise unprecedented flexibility in providing materials with very complex specifications, including: independent control of the permittivity and permeability with both positive and negative values, anisotropy control and designed gradients[12, 8]. These materials should lead to electromagnetic devices spanning a broader scope of capability, but this is quite a bit more flexibility than the optical designer is accustomed to. It is not clear how to begin designing with such complex media. The transformation method provides an intuitive and direct way to design in this enlarged material parameter space. The designer imagines a fictitious space with some topological feature (e.g., a hole) that enacts a desired electromagnetic phenomena (e.g., invisibility), and the transformation method yields, in a direct way, the complex material property specification that implements the behavior.

2. Material properties

We will describe the transformation properties of the electromagnetic material property tensors. We first note that the Minkowski form of Maxwell’s equations [13, 14, 15]

Fαβ,μ+Fβμ,α+Fμα,β=0
G,ααβ=Jβ

is form invariant for general space-time transformations. F αβ is the tensor of electric field and magnetic induction, and G αβ is the tensor density of electric displacement and magnetic field, and J β is the source vector. In component form these tensors are

(Fαβ)=(0E1E2E3E10cB3cB2E2cB30cB1E3cB2cB10)
(Gαβ)=(0cD1cD2cD3cD10H3H2cD2H30H1cD3H2H10)
(Jβ)=(J1J2J3)

where we use SI units (as we will throughout this article), our space-time coordinate vector is (x α)=(ct,x,y,z), and we use the metric signature, +2. All of the information regarding the topology of the space is contained in the constitutive relations

Gαβ=12CαβμνFμν

where C αβµν is the constitutive tensor representing the properties of the medium, including its permittivity, permeability and bianisotropic properties. C αβµν is a tensor density of weight +1, so it transforms as [13]

Cαβμν=det(Λαα)1ΛααΛββΛμμΛννCαβμν

written in terms of the Jacobian transformation matrix

Λαα=xαxα

which is just the derivative of the transformed coordinates with respect to the original coordinates. If we restrict ourselves to transformations that are time invariant, the permittivity and permeability are also tensors individually. Specifically, they are tensor densities of weight +1, which transform as [16, 13]

εij=det(Λii)1ΛiiΛjjεij
μij=det(Λii)1ΛiiΛjjεij

where the roman indices run from1 to 3, for the three spatial coordinates, as is standard practice. Equations (6) are the primary tools for the transformation design method when the base medium does not possess magnetoelectric coupling and the desired device moves or changes shape with speeds much less than that of light, i.e. most devices of practical interest. These equations can be shown to be exactly equivalent to the results derived by Ward and Pendry [17].

If the original medium is isotropic, Eq. (6) can also be written in terms of the metric [16, 11]

εij=det(gij)12gijε
μij=det(gij)12gijμ

where the metric is given by

gij=ΛkiΛljδkl

Maxwell’s equations, Eqs. (1), together with the medium specified by Eq. (4) or Eqs. (6) describe a single electromagnetic behavior, but this behavior can be interpreted in two ways.

One way, the traditional way, is that the material property tensors that appear on the left and right hand sides of Eqs. (6) represent the same material properties, but in different spaces. The components in the transformed space are different form those in the original space, due to the topology of the transformation. We will refer to this view as the topological interpretation.

An alternative interpretation, is that the material property tensors on the left and right hand sides of Eqs. (6) represent different material properties. Both sets of tensor components are interpreted as components in a flat, Cartesian space. The form invariance of Maxwell’s equations insures that both interpretations lead to the same electromagnetic behavior. We will refer to this view as the materials interpretation.

To design something of interest, one imagines a space with some desired property, a hole for example. Then one constructs the coordinate transformation of the space with this desired property. Using Eq. (4) or Eqs. (6) one can then calculate a set of material properties that will implement this interesting property of the imagined space in our own boring, flat, Cartesian space.

2.1. Spherical cloak

The spherical cloak is designed by considering a spherically symmetric coordinate transformation. This transformation compresses all the space in a volume of radius, b, into a spherical shell of inner radius, a, and outer radius, b. Consider a position vector, x. In the original coordinate system (Fig.1(a)) it has components, xi, and in the transformed coordinate system (Fig. 1(b)), xi′. Of course, its magnitude, r, is independent of coordinate system

r=(xixjδij)12=(xixjgij)12

where gi′j′ is the metric of the transformed space. In the materials interpretation, (Fig.1(c)), we consider the components, xi′, to be the components of a Cartesian vector, and its magnitude, which we will call r′, is found using the appropriate flat space metric

r=(xixjδij)12

Perhaps the simplest spherical cloak transformation maps points from a radius, r, to a radius, r′, according to the following linear function

r=babr+a

which we apply over the domain, 0≤rb, (or equivalently, ar′b). Outside this domain we assume the identity transformation. (All equations in the remainder of this article apply only to the transformation domain.) We must always limit the transformation to apply only over a finite region of space if we wish to implement it with materials of finite extent. Note that when r=0 then r′=a, so that the origin is mapped out to a finite radius, opening up a whole in space. Note also that when r=b then r′=b, so that space at the outer boundary of the transformation is undistorted and there is no discontinuity with the space outside the transformation domain.

 figure: Fig. 1.

Fig. 1. The thick blue line shows the path of the same ray in (A) the original Cartesian space, and under two different interpretations of the electromagnetic equations, (B) the topological interpretation and (C) the materials interpretation. The position vector x is shown in both the original and transformed spaces, and the length of the vector where the transformed components are interpreted as Cartesian components is shown in (C).

Download Full Size | PPT Slide | PDF

Now since our transformation is radially symmetric, the unit vectors in materials interpretation and in the original space must be equal.

xir=xirδii

Expressing the components of the position vector in the transformed space in terms of only the components in the original space, using Eq. 12, we obtain.

xi=babxiδii+axirδii

Now that we have this expression, we need not worry about the interpretations of transformed space, we can just proceed in standard fashion to compute the transformation matrix. To take the derivative of this expression we note that

xjxir=xixkδkjr3+1rδji

and obtain the transformation matrix

Λji=rrδjiaxixkδjiδkjr3

The components of this expression written out are

(Λji)=(rrax2r3axyr3axzr3ayxr3rray2r3ayzr3azxr3azyr3rraz2r3)

To calculate the determinant of this matrix we note that we can always rotate into a coordinate system such that

(xi)=(r,0,0)

then the determinant is, by inspection, given by

det(Λji)=rar(rr)2

If we assume that our original medium is free space, then the relative permittivity and permeability will be equal to each other. Working out the algebra, we find that the material properties are then given by

εij=μij=bba[δij2ara2r4xixj]

where we have eliminated any dependence on the components of x in the original space, xi, or the magnitude, r. We can now drop the primes for aesthetic reasons, and we need not make the distinction between vectors and one-forms as we consider this to be a material specification in flat, Cartesian, three-space, where such distinctions are not necessary. Writing this expression in direct notation

ε=μ=bba(I2ara2r4rr)

where rr is the outer product of the position vector with itself, also referred to as a dyad formed from the position vector. We note, for later use, that the determinant can be easily calculated, as above, using an appropriately chosen rotation

det(ε)=det(μ)=(bba)3(rar)2

2.2. Cylindrical cloak

To analyze a cylindrical cloak we will use two projection operators. One which projects onto the cylinder’s axis, (which we will call the third coordinate or z-axis), and one that projects onto the plane normal to the cylinder’s axis.

Zij=δ3iδ3j
Tij=δ1iδ1j+δ2iδ2j

We do not mean to imply that these are tensors. We define these operators to perform these projections onto the third coordinate and the plane normal to the third coordinate in whatever basis (including mixed bases) they are applied to. Thus we will refer to their components with indices up or down, primed or un-primed, at will. We now use the transverse projection operator to define a transverse coordinate.

ρi=Tjixj

The coordinate transformation for the cylindrical case is the same as that of the spherical case in the two dimensions normal to the cylinder’s axis. Along the cylinders axis the transformation is the identity. Thus we have for the transformation matrix.

(Λji)=(ρρax2ρ3axyρ30ayxρ3ρρay2ρ30001)

or written in index form

Λji=ρρTjiaρiρkδiiδkjρ3+Zji

Again, we can easily calculate the determinant by rotating into a coordinate system where

(xi)=(ρi)=(ρ,0,0)

then we find the determinant to be

det(Λji)=ρaρρρ

The material properties in direct notation and dropping the primes are

ε=μ=ρρaT2a2ρ3(ρa)ρρ+(bba)2ρaρZ

Again we note the determinant for later use, which takes the rather simple form

det(ε)=det(μ)=(bba)2ρaρ

3. Hamiltonian and ray equations

The Hamiltonian we will use for generating the ray paths is essentially the plane wave dispersion relation [18]. We derive it here, briefly, to show our choice of dimensionality for the relevant variables. We begin with Maxwell’s curl equations in SI units

×E=Bt×H=Dt

We assume plane wave solutions with slowly varying coefficients, appropriate for the geometric limit

E=E0ei(k0k·xωt)H=1η0H0ei(k0k·xωt)

Here η0=μ0ε0 is the impedance of free space, giving E 0 and H 0 the same units, and k 0=ω/c making k dimensionless. We use constitutive relations with dimensionless tensors ε and µ.

D=ε0εEB=μ0μH

Plugging Eqs. (32) and Eqs. (33) into the curl equations Eqs. (31) we obtain

k×E0μH0=0k×H0+εE0=0

Eliminating the magnetic field we find

k×(μ1(k×E0))+εE0=0

Defining the operator, K [15],

Kikεijkkj

the dispersion relation Eq. (35) can be expressed as a single operator on E 0,

(Kμ1K+ε)E0=0

which must be singular for non-zero field solutions. The dispersion relation expresses that this operator must have zero determinate.

det(Kμ1K+ε)=0

Now for material properties derived from transforming free space, ε and µ are the same symmetric tensor, which we will call n. In this case the dispersion relation has an alternate expression

det(Kn1K+n)=1det(n)(knkdet(n))2

This can be proved simply, by evaluating the right and left hand sides in a computer algebra system using arbitrary symmetric components for n, or perhaps some other clever way. The latter expression is clearly fourth order in k, but has only two unique solutions. Thus we discover that media with ε=µ is singly refracting [19], unlike for example, uniaxial dielectrics which exhibit an ordinary and extraordinary ray. This can also be seen by noting that free space is singly refracting. A coordinate transformation cannot separate two degenerate ray paths, so the degenerate ray paths of free space will remain so in the transformed coordinate space and thus also in the equivalent media.

The Hamiltonian then easily factors into two terms that represent degenerate modes. Further it is easy to show (by plugging Eq. (40) into Eqs. (41)) that the Hamiltonian may be multiplied by an arbitrary function of the spatial coordinates without changing the paths obtained from the equations of motion, (only the parameterization is changed), thus we can drop the factor, 1/det (n), and our Hamiltonian is

H=f(x)(knkdet(n))

where f(x) is some arbitrary function of position. The equations of motion are [18]

dxdτ=Hk
dkdτ=Hx

where τ parameterizes the paths. This pair of coupled, first order, ordinary differential equations can be integrated using a standard solver, such as Mathematica’s NDSolve.

4. Refraction

The equations of motion, Eqs. (41), govern the path of the ray throughout the continuously varying inhomogenous media. At discontinuities, such as the outer boundary of a cloak, we must perform boundary matching to determine the discontinuous change in direction of the ray, i.e. refraction. Given k 1 on one side of the boundary we find k 2 on the other side as follows. The transverse component of the wave vector is conserved across the boundary.

(k1k2)×n=0

where here n is the unit normal to the boundary. This vector equation represents just two equations. The third is obtained by requiring the wave vector to satisfy the plane wave dispersion relation of the mode represented by the Hamiltonian.

H(k2)=0

These three equations determine the three unknowns of the vector components of k 2. Since H is quadratic in k, there will be two solutions, one that carries energy into medium 2, the desired solution, and one that carries energy out. The path of the ray, d x/, indicates the direction of energy flow, so the Hamiltonian can be used to determine which is the desired solution. The desired solution satisfies

Hk·n>0

if n is the normal pointing into medium 2. These equations apply equally well to refraction into or out of transformation media. Refracting out into free space is much easier since the Hamiltonian of free space is just, H=k·k-1.

Completely tracing a ray that intersects a cloak involves several steps. The ray is assigned some initial direction and point of origin in the surrounding environment. The ray is traced to its intersection point with the outer boundary of the cloak, (which may be just a straight line intersection if the surrounding environment is homogenous.) Then the ray is refracted into the cloak domain, yielding the initial conditions for integrating the equations of motion, Eqs. (41). The integration is continued until the ray has traversed the cloak and intersects the outer boundary a second time. Then a second refraction yields the direction of the ray as it re-enters the surrounding environment. For the cloaks analyzed in this article, the exiting segment of the ray should be collinear with the initial segment of the ray.

5. Cloak Hamiltonians

We now show specific examples of ray tracing. Below we will choose a specific form for the Hamiltonian, plug in the material properties and display the derivatives of the Hamiltonian, for both the spherical and cylindrical cloak.

5.1. Spherical cloak

For the spherical cloak (Fig.2), the Hamiltonian which yields the simplest equations is

H=12bab(knkdet(n))

Plugging in the material properties from Eq. (21) and Eq. (22) we obtain

H=12k·k122ara2r4(x·k)212[b(ra)r(ba)]2

Taking the derivatives, (which is straight forward particularly in index form), yields

Hk=k2ara2r4(x·k)x
Hx=2ara2r4(x·k)k+3ar2a2r6(x·k)2x(bba)2(ara2r4)x
 figure: Fig. 2.

Fig. 2. Rays traversing a spherical cloak. The transformation media that comprises the cloak lies between the two spheres.

Download Full Size | PPT Slide | PDF

5.2. Cylindrical cloak

For the cylindrical cloak (Fig. 3), the Hamiltonian which yields the simplest equations is

H=12ρaρ(knkdet(n))

Plugging in the material properties form Eq. (29) and Eq. (30) we obtain

H=12kTk122a2ρ4(ρ·k)2+12[b(ρa)ρ(ba)]2(kZk1)

For taking the derivatives we note that the derivative of the transverse position vector with respect to the position vector is the transverse projection operator.

ρx=T
 figure: Fig. 3.

Fig. 3. Rays traversing a cylindrical cloak at an oblique angle. The transformation media that comprises the cloak lies in an annular region between the cylinders.

Download Full Size | PPT Slide | PDF

The derivatives are thus

Hk=Tk2a2ρ4(ρ·k)ρ+[b(ρa)ρ(ba)]2Zk
Hx=32a2ρ6(ρ·k)2ρ2a2ρ4(ρ·k)Tk+(bba)2a2ρ4(kZk1)ρ

6. Conclusion

We have shown how to calculate the material properties associated with a coordinate transformation and use these properties to perform ray tracing. Examples, of spherical and cylindrical cloaks are worked out in some detail. Some of the value in this effort is to provide independent confirmation that the material properties calculated from the transformation do indeed cause electromagnetic waves to behave in the desired and predicted manner. Eventually, the transformation technique will become more accepted and independent confirmation will not be needed. One can see what the waves will do much more easily by applying the transformation to the rays or fields in the original space where the behavior is simpler if not trivial. However, one may still want to perform ray tracing on these media to see the effects of perturbations from the ideal material specification.

We believe that transformation optical design will prove to be a usefulmethodology. Together with advancing metamaterial technology it should lead to the realization of devices that would be difficult, if not impossible, to conceive and fabricate any other way.

Acknowledgments

D. Schurig wishes to acknowledge the IC Postdoctoral Research Fellowship program, and J. B. Pendry wishes to thank the EPSRC for a Senior Fellowship.

References and links

1. J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling Electromagnetic Fields,” Science 312, 1780 (2006). [CrossRef]   [PubMed]  

2. U. Leonhardt, “Optical conformal mapping,” Science 312, 1777 (2006). [CrossRef]   [PubMed]  

3. D. R. Smith, J. B. Pendry, and M. C. K. Wiltshire, “Metamaterials and Negative Refractive Index,” Science 305, 788 (2004). [CrossRef]   [PubMed]  

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

5. T. J. Yen, W. J. Padilla, N. Fang, D. C. Vier, D. R. Smith, J. B. Pendry, D. N. Basov, and X. Zhang, “Terahertz Magnetic Response from Artificial Materials,” Science 303, 1494–1496 (2004). [CrossRef]   [PubMed]  

6. S. Linden, C. Enkrich, M. Wegener, J. Zhou, T. Koschny, and C. M. Soukoulis, “Magnetic Response of Metamaterials at 100 Terahertz,” Science 306, 1351–1353 (2004). [CrossRef]   [PubMed]  

7. D. Schurig, J. J. Mock, and D. R. Smith, “Electric-field-coupled resonators for negative permittivity metamaterials,” Appl. Phys. Lett. 88(4), 041,109- (2006).

8. D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science (2006). In press. [CrossRef]   [PubMed]  

9. A. Alu and N. Engheta, “Achieving transparency with plasmonic and metamaterial coatings,” Phys. Rev. E 72, 016,623 (2005). [CrossRef]  

10. G. W. Milton and N.-A. P. Nicorovici, “On the cloaking effects associated with anomalous localized resonance,” Proc. Roy. Soc. London A 462, 1364 (2006).

11. U. Leonhardt and T. G. Philbin, “General relativity in electrical engineering,” (2006). http://xxx.arxiv.org/abs/cond-mat/0607418.

12. T. Driscoll, D. N. Basov, A. F. Starr, P. M. Rye, S. Nemat-Nasser, D. Schurig, and D. R. Smith, “Free-space microwave focusing by a negative-index gradient lens,” Appl. Phys. Lett. 88, 081,101 (2006). [CrossRef]  

13. E.J. Post, Formal structure of electromagnetics (Wiley, New York, 1962).

14. J. D. Jackson, Classical Electrodynamics (Wiley, New York, 1975).

15. J. A. Kong, Electromagnetic Wave Theory, 2nd ed. (Wiley-Interscience, New York, 1990).

16. D. M. Shyroki, “Exact equivalent-profile formulation for bent optical waveguides,” (2006). Unpublished.

17. A. J. Ward and J. B. Pendry, “Refraction and geometry in maxwell’s equations,” J. Mod. Opt. 43(4), 773–793 (1996). [CrossRef]  

18. Y. A. Kravtsov and Y. I. Orlov, Geometrical optics of inhomogeneous media (Springer-Verlag, Berlin, 1990). [CrossRef]  

19. H. C. Chen, Theory of electromagnetic waves: A coordinate free approach, pp. 216–218 (McGraw-Hill, 1985).

References

  • View by:
  • |
  • |
  • |

  1. J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling Electromagnetic Fields,” Science 312, 1780 (2006).
    [Crossref] [PubMed]
  2. U. Leonhardt, “Optical conformal mapping,” Science 312, 1777 (2006).
    [Crossref] [PubMed]
  3. D. R. Smith, J. B. Pendry, and M. C. K. Wiltshire, “Metamaterials and Negative Refractive Index,” Science 305, 788 (2004).
    [Crossref] [PubMed]
  4. E. Cubukcu, K. Aydin, E. Ozbay, S. Foteinopoulou, and C. M. Soukoulis, “Electromagnetic waves: Negative refraction by photonic crystals,” Nature 423, 604 (2003).
    [Crossref] [PubMed]
  5. T. J. Yen, W. J. Padilla, N. Fang, D. C. Vier, D. R. Smith, J. B. Pendry, D. N. Basov, and X. Zhang, “Terahertz Magnetic Response from Artificial Materials,” Science 303, 1494–1496 (2004).
    [Crossref] [PubMed]
  6. S. Linden, C. Enkrich, M. Wegener, J. Zhou, T. Koschny, and C. M. Soukoulis, “Magnetic Response of Metamaterials at 100 Terahertz,” Science 306, 1351–1353 (2004).
    [Crossref] [PubMed]
  7. D. Schurig, J. J. Mock, and D. R. Smith, “Electric-field-coupled resonators for negative permittivity metamaterials,” Appl. Phys. Lett. 88(4), 041,109- (2006).
  8. D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science (2006). In press.
    [Crossref] [PubMed]
  9. A. Alu and N. Engheta, “Achieving transparency with plasmonic and metamaterial coatings,” Phys. Rev. E 72, 016,623 (2005).
    [Crossref]
  10. G. W. Milton and N.-A. P. Nicorovici, “On the cloaking effects associated with anomalous localized resonance,” Proc. Roy. Soc. London A 462, 1364 (2006).
  11. U. Leonhardt and T. G. Philbin, “General relativity in electrical engineering,” (2006). http://xxx.arxiv.org/abs/cond-mat/0607418.
  12. T. Driscoll, D. N. Basov, A. F. Starr, P. M. Rye, S. Nemat-Nasser, D. Schurig, and D. R. Smith, “Free-space microwave focusing by a negative-index gradient lens,” Appl. Phys. Lett. 88, 081,101 (2006).
    [Crossref]
  13. E.J. Post, Formal structure of electromagnetics (Wiley, New York, 1962).
  14. J. D. Jackson, Classical Electrodynamics (Wiley, New York, 1975).
  15. J. A. Kong, Electromagnetic Wave Theory, 2nd ed. (Wiley-Interscience, New York, 1990).
  16. D. M. Shyroki, “Exact equivalent-profile formulation for bent optical waveguides,” (2006). Unpublished.
  17. A. J. Ward and J. B. Pendry, “Refraction and geometry in maxwell’s equations,” J. Mod. Opt. 43(4), 773–793 (1996).
    [Crossref]
  18. Y. A. Kravtsov and Y. I. Orlov, Geometrical optics of inhomogeneous media (Springer-Verlag, Berlin, 1990).
    [Crossref]
  19. H. C. Chen, Theory of electromagnetic waves: A coordinate free approach, pp. 216–218 (McGraw-Hill, 1985).

2006 (5)

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling Electromagnetic Fields,” Science 312, 1780 (2006).
[Crossref] [PubMed]

U. Leonhardt, “Optical conformal mapping,” Science 312, 1777 (2006).
[Crossref] [PubMed]

D. Schurig, J. J. Mock, and D. R. Smith, “Electric-field-coupled resonators for negative permittivity metamaterials,” Appl. Phys. Lett. 88(4), 041,109- (2006).

G. W. Milton and N.-A. P. Nicorovici, “On the cloaking effects associated with anomalous localized resonance,” Proc. Roy. Soc. London A 462, 1364 (2006).

T. Driscoll, D. N. Basov, A. F. Starr, P. M. Rye, S. Nemat-Nasser, D. Schurig, and D. R. Smith, “Free-space microwave focusing by a negative-index gradient lens,” Appl. Phys. Lett. 88, 081,101 (2006).
[Crossref]

2005 (1)

A. Alu and N. Engheta, “Achieving transparency with plasmonic and metamaterial coatings,” Phys. Rev. E 72, 016,623 (2005).
[Crossref]

2004 (3)

T. J. Yen, W. J. Padilla, N. Fang, D. C. Vier, D. R. Smith, J. B. Pendry, D. N. Basov, and X. Zhang, “Terahertz Magnetic Response from Artificial Materials,” Science 303, 1494–1496 (2004).
[Crossref] [PubMed]

S. Linden, C. Enkrich, M. Wegener, J. Zhou, T. Koschny, and C. M. Soukoulis, “Magnetic Response of Metamaterials at 100 Terahertz,” Science 306, 1351–1353 (2004).
[Crossref] [PubMed]

D. R. Smith, J. B. Pendry, and M. C. K. Wiltshire, “Metamaterials and Negative Refractive Index,” Science 305, 788 (2004).
[Crossref] [PubMed]

2003 (1)

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

1996 (1)

A. J. Ward and J. B. Pendry, “Refraction and geometry in maxwell’s equations,” J. Mod. Opt. 43(4), 773–793 (1996).
[Crossref]

Alu, A.

A. Alu and N. Engheta, “Achieving transparency with plasmonic and metamaterial coatings,” Phys. Rev. E 72, 016,623 (2005).
[Crossref]

Aydin, K.

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

Basov, D. N.

T. Driscoll, D. N. Basov, A. F. Starr, P. M. Rye, S. Nemat-Nasser, D. Schurig, and D. R. Smith, “Free-space microwave focusing by a negative-index gradient lens,” Appl. Phys. Lett. 88, 081,101 (2006).
[Crossref]

T. J. Yen, W. J. Padilla, N. Fang, D. C. Vier, D. R. Smith, J. B. Pendry, D. N. Basov, and X. Zhang, “Terahertz Magnetic Response from Artificial Materials,” Science 303, 1494–1496 (2004).
[Crossref] [PubMed]

Chen, H. C.

H. C. Chen, Theory of electromagnetic waves: A coordinate free approach, pp. 216–218 (McGraw-Hill, 1985).

Cubukcu, E.

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

Cummer, S. A.

D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science (2006). In press.
[Crossref] [PubMed]

Driscoll, T.

T. Driscoll, D. N. Basov, A. F. Starr, P. M. Rye, S. Nemat-Nasser, D. Schurig, and D. R. Smith, “Free-space microwave focusing by a negative-index gradient lens,” Appl. Phys. Lett. 88, 081,101 (2006).
[Crossref]

Engheta, N.

A. Alu and N. Engheta, “Achieving transparency with plasmonic and metamaterial coatings,” Phys. Rev. E 72, 016,623 (2005).
[Crossref]

Enkrich, C.

S. Linden, C. Enkrich, M. Wegener, J. Zhou, T. Koschny, and C. M. Soukoulis, “Magnetic Response of Metamaterials at 100 Terahertz,” Science 306, 1351–1353 (2004).
[Crossref] [PubMed]

Fang, N.

T. J. Yen, W. J. Padilla, N. Fang, D. C. Vier, D. R. Smith, J. B. Pendry, D. N. Basov, and X. Zhang, “Terahertz Magnetic Response from Artificial Materials,” Science 303, 1494–1496 (2004).
[Crossref] [PubMed]

Foteinopoulou, S.

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

Jackson, J. D.

J. D. Jackson, Classical Electrodynamics (Wiley, New York, 1975).

Justice, B. J.

D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science (2006). In press.
[Crossref] [PubMed]

Kong, J. A.

J. A. Kong, Electromagnetic Wave Theory, 2nd ed. (Wiley-Interscience, New York, 1990).

Koschny, T.

S. Linden, C. Enkrich, M. Wegener, J. Zhou, T. Koschny, and C. M. Soukoulis, “Magnetic Response of Metamaterials at 100 Terahertz,” Science 306, 1351–1353 (2004).
[Crossref] [PubMed]

Kravtsov, Y. A.

Y. A. Kravtsov and Y. I. Orlov, Geometrical optics of inhomogeneous media (Springer-Verlag, Berlin, 1990).
[Crossref]

Leonhardt, U.

U. Leonhardt, “Optical conformal mapping,” Science 312, 1777 (2006).
[Crossref] [PubMed]

U. Leonhardt and T. G. Philbin, “General relativity in electrical engineering,” (2006). http://xxx.arxiv.org/abs/cond-mat/0607418.

Linden, S.

S. Linden, C. Enkrich, M. Wegener, J. Zhou, T. Koschny, and C. M. Soukoulis, “Magnetic Response of Metamaterials at 100 Terahertz,” Science 306, 1351–1353 (2004).
[Crossref] [PubMed]

Milton, G. W.

G. W. Milton and N.-A. P. Nicorovici, “On the cloaking effects associated with anomalous localized resonance,” Proc. Roy. Soc. London A 462, 1364 (2006).

Mock, J. J.

D. Schurig, J. J. Mock, and D. R. Smith, “Electric-field-coupled resonators for negative permittivity metamaterials,” Appl. Phys. Lett. 88(4), 041,109- (2006).

D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science (2006). In press.
[Crossref] [PubMed]

Nemat-Nasser, S.

T. Driscoll, D. N. Basov, A. F. Starr, P. M. Rye, S. Nemat-Nasser, D. Schurig, and D. R. Smith, “Free-space microwave focusing by a negative-index gradient lens,” Appl. Phys. Lett. 88, 081,101 (2006).
[Crossref]

Nicorovici, N.-A. P.

G. W. Milton and N.-A. P. Nicorovici, “On the cloaking effects associated with anomalous localized resonance,” Proc. Roy. Soc. London A 462, 1364 (2006).

Orlov, Y. I.

Y. A. Kravtsov and Y. I. Orlov, Geometrical optics of inhomogeneous media (Springer-Verlag, Berlin, 1990).
[Crossref]

Ozbay, E.

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

Padilla, W. J.

T. J. Yen, W. J. Padilla, N. Fang, D. C. Vier, D. R. Smith, J. B. Pendry, D. N. Basov, and X. Zhang, “Terahertz Magnetic Response from Artificial Materials,” Science 303, 1494–1496 (2004).
[Crossref] [PubMed]

Pendry, J. B.

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling Electromagnetic Fields,” Science 312, 1780 (2006).
[Crossref] [PubMed]

D. R. Smith, J. B. Pendry, and M. C. K. Wiltshire, “Metamaterials and Negative Refractive Index,” Science 305, 788 (2004).
[Crossref] [PubMed]

T. J. Yen, W. J. Padilla, N. Fang, D. C. Vier, D. R. Smith, J. B. Pendry, D. N. Basov, and X. Zhang, “Terahertz Magnetic Response from Artificial Materials,” Science 303, 1494–1496 (2004).
[Crossref] [PubMed]

A. J. Ward and J. B. Pendry, “Refraction and geometry in maxwell’s equations,” J. Mod. Opt. 43(4), 773–793 (1996).
[Crossref]

D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science (2006). In press.
[Crossref] [PubMed]

Philbin, T. G.

U. Leonhardt and T. G. Philbin, “General relativity in electrical engineering,” (2006). http://xxx.arxiv.org/abs/cond-mat/0607418.

Post, E.J.

E.J. Post, Formal structure of electromagnetics (Wiley, New York, 1962).

Rye, P. M.

T. Driscoll, D. N. Basov, A. F. Starr, P. M. Rye, S. Nemat-Nasser, D. Schurig, and D. R. Smith, “Free-space microwave focusing by a negative-index gradient lens,” Appl. Phys. Lett. 88, 081,101 (2006).
[Crossref]

Schurig, D.

T. Driscoll, D. N. Basov, A. F. Starr, P. M. Rye, S. Nemat-Nasser, D. Schurig, and D. R. Smith, “Free-space microwave focusing by a negative-index gradient lens,” Appl. Phys. Lett. 88, 081,101 (2006).
[Crossref]

D. Schurig, J. J. Mock, and D. R. Smith, “Electric-field-coupled resonators for negative permittivity metamaterials,” Appl. Phys. Lett. 88(4), 041,109- (2006).

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling Electromagnetic Fields,” Science 312, 1780 (2006).
[Crossref] [PubMed]

D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science (2006). In press.
[Crossref] [PubMed]

Shyroki, D. M.

D. M. Shyroki, “Exact equivalent-profile formulation for bent optical waveguides,” (2006). Unpublished.

Smith, D. R.

T. Driscoll, D. N. Basov, A. F. Starr, P. M. Rye, S. Nemat-Nasser, D. Schurig, and D. R. Smith, “Free-space microwave focusing by a negative-index gradient lens,” Appl. Phys. Lett. 88, 081,101 (2006).
[Crossref]

D. Schurig, J. J. Mock, and D. R. Smith, “Electric-field-coupled resonators for negative permittivity metamaterials,” Appl. Phys. Lett. 88(4), 041,109- (2006).

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling Electromagnetic Fields,” Science 312, 1780 (2006).
[Crossref] [PubMed]

D. R. Smith, J. B. Pendry, and M. C. K. Wiltshire, “Metamaterials and Negative Refractive Index,” Science 305, 788 (2004).
[Crossref] [PubMed]

T. J. Yen, W. J. Padilla, N. Fang, D. C. Vier, D. R. Smith, J. B. Pendry, D. N. Basov, and X. Zhang, “Terahertz Magnetic Response from Artificial Materials,” Science 303, 1494–1496 (2004).
[Crossref] [PubMed]

D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science (2006). In press.
[Crossref] [PubMed]

Soukoulis, C. M.

S. Linden, C. Enkrich, M. Wegener, J. Zhou, T. Koschny, and C. M. Soukoulis, “Magnetic Response of Metamaterials at 100 Terahertz,” Science 306, 1351–1353 (2004).
[Crossref] [PubMed]

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

Starr, A. F.

T. Driscoll, D. N. Basov, A. F. Starr, P. M. Rye, S. Nemat-Nasser, D. Schurig, and D. R. Smith, “Free-space microwave focusing by a negative-index gradient lens,” Appl. Phys. Lett. 88, 081,101 (2006).
[Crossref]

D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science (2006). In press.
[Crossref] [PubMed]

Vier, D. C.

T. J. Yen, W. J. Padilla, N. Fang, D. C. Vier, D. R. Smith, J. B. Pendry, D. N. Basov, and X. Zhang, “Terahertz Magnetic Response from Artificial Materials,” Science 303, 1494–1496 (2004).
[Crossref] [PubMed]

Ward, A. J.

A. J. Ward and J. B. Pendry, “Refraction and geometry in maxwell’s equations,” J. Mod. Opt. 43(4), 773–793 (1996).
[Crossref]

Wegener, M.

S. Linden, C. Enkrich, M. Wegener, J. Zhou, T. Koschny, and C. M. Soukoulis, “Magnetic Response of Metamaterials at 100 Terahertz,” Science 306, 1351–1353 (2004).
[Crossref] [PubMed]

Wiltshire, M. C. K.

D. R. Smith, J. B. Pendry, and M. C. K. Wiltshire, “Metamaterials and Negative Refractive Index,” Science 305, 788 (2004).
[Crossref] [PubMed]

Yen, T. J.

T. J. Yen, W. J. Padilla, N. Fang, D. C. Vier, D. R. Smith, J. B. Pendry, D. N. Basov, and X. Zhang, “Terahertz Magnetic Response from Artificial Materials,” Science 303, 1494–1496 (2004).
[Crossref] [PubMed]

Zhang, X.

T. J. Yen, W. J. Padilla, N. Fang, D. C. Vier, D. R. Smith, J. B. Pendry, D. N. Basov, and X. Zhang, “Terahertz Magnetic Response from Artificial Materials,” Science 303, 1494–1496 (2004).
[Crossref] [PubMed]

Zhou, J.

S. Linden, C. Enkrich, M. Wegener, J. Zhou, T. Koschny, and C. M. Soukoulis, “Magnetic Response of Metamaterials at 100 Terahertz,” Science 306, 1351–1353 (2004).
[Crossref] [PubMed]

Appl. Phys. Lett. (2)

D. Schurig, J. J. Mock, and D. R. Smith, “Electric-field-coupled resonators for negative permittivity metamaterials,” Appl. Phys. Lett. 88(4), 041,109- (2006).

T. Driscoll, D. N. Basov, A. F. Starr, P. M. Rye, S. Nemat-Nasser, D. Schurig, and D. R. Smith, “Free-space microwave focusing by a negative-index gradient lens,” Appl. Phys. Lett. 88, 081,101 (2006).
[Crossref]

J. Mod. Opt. (1)

A. J. Ward and J. B. Pendry, “Refraction and geometry in maxwell’s equations,” J. Mod. Opt. 43(4), 773–793 (1996).
[Crossref]

Nature (1)

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

Phys. Rev. E (1)

A. Alu and N. Engheta, “Achieving transparency with plasmonic and metamaterial coatings,” Phys. Rev. E 72, 016,623 (2005).
[Crossref]

Proc. Roy. Soc. London A (1)

G. W. Milton and N.-A. P. Nicorovici, “On the cloaking effects associated with anomalous localized resonance,” Proc. Roy. Soc. London A 462, 1364 (2006).

Science (5)

T. J. Yen, W. J. Padilla, N. Fang, D. C. Vier, D. R. Smith, J. B. Pendry, D. N. Basov, and X. Zhang, “Terahertz Magnetic Response from Artificial Materials,” Science 303, 1494–1496 (2004).
[Crossref] [PubMed]

S. Linden, C. Enkrich, M. Wegener, J. Zhou, T. Koschny, and C. M. Soukoulis, “Magnetic Response of Metamaterials at 100 Terahertz,” Science 306, 1351–1353 (2004).
[Crossref] [PubMed]

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling Electromagnetic Fields,” Science 312, 1780 (2006).
[Crossref] [PubMed]

U. Leonhardt, “Optical conformal mapping,” Science 312, 1777 (2006).
[Crossref] [PubMed]

D. R. Smith, J. B. Pendry, and M. C. K. Wiltshire, “Metamaterials and Negative Refractive Index,” Science 305, 788 (2004).
[Crossref] [PubMed]

Other (8)

U. Leonhardt and T. G. Philbin, “General relativity in electrical engineering,” (2006). http://xxx.arxiv.org/abs/cond-mat/0607418.

D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science (2006). In press.
[Crossref] [PubMed]

Y. A. Kravtsov and Y. I. Orlov, Geometrical optics of inhomogeneous media (Springer-Verlag, Berlin, 1990).
[Crossref]

H. C. Chen, Theory of electromagnetic waves: A coordinate free approach, pp. 216–218 (McGraw-Hill, 1985).

E.J. Post, Formal structure of electromagnetics (Wiley, New York, 1962).

J. D. Jackson, Classical Electrodynamics (Wiley, New York, 1975).

J. A. Kong, Electromagnetic Wave Theory, 2nd ed. (Wiley-Interscience, New York, 1990).

D. M. Shyroki, “Exact equivalent-profile formulation for bent optical waveguides,” (2006). Unpublished.

Cited By

OSA participates in Crossref's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (3)

Fig. 1.
Fig. 1. The thick blue line shows the path of the same ray in (A) the original Cartesian space, and under two different interpretations of the electromagnetic equations, (B) the topological interpretation and (C) the materials interpretation. The position vector x is shown in both the original and transformed spaces, and the length of the vector where the transformed components are interpreted as Cartesian components is shown in (C).
Fig. 2.
Fig. 2. Rays traversing a spherical cloak. The transformation media that comprises the cloak lies between the two spheres.
Fig. 3.
Fig. 3. Rays traversing a cylindrical cloak at an oblique angle. The transformation media that comprises the cloak lies in an annular region between the cylinders.

Equations (58)

Equations on this page are rendered with MathJax. Learn more.

F α β , μ + F β μ , α + F μ α , β = 0
G , α α β = J β
( F α β ) = ( 0 E 1 E 2 E 3 E 1 0 c B 3 c B 2 E 2 c B 3 0 c B 1 E 3 c B 2 c B 1 0 )
( G α β ) = ( 0 c D 1 c D 2 c D 3 c D 1 0 H 3 H 2 c D 2 H 3 0 H 1 c D 3 H 2 H 1 0 )
( J β ) = ( J 1 J 2 J 3 )
G α β = 1 2 C α β μ ν F μ ν
C α β μ ν = det ( Λ α α ) 1 Λ α α Λ β β Λ μ μ Λ ν ν C α β μ ν
Λ α α = x α x α
ε i j = det ( Λ i i ) 1 Λ i i Λ j j ε ij
μ i j = det ( Λ i i ) 1 Λ i i Λ j j ε ij
ε i j = det ( g i j ) 1 2 g i j ε
μ i j = det ( g i j ) 1 2 g i j μ
g i j = Λ k i Λ l j δ kl
r = ( x i x j δ ij ) 1 2 = ( x i x j g i j ) 1 2
r = ( x i x j δ i j ) 1 2
r = b a b r + a
x i r = x i r δ i i
x i = b a b x i δ i i + a x i r δ i i
x j x i r = x i x k δ kj r 3 + 1 r δ j i
Λ j i = r r δ j i a x i x k δ j i δ kj r 3
( Λ j i ) = ( r r a x 2 r 3 axy r 3 axz r 3 ayx r 3 r r ay 2 r 3 ayz r 3 azx r 3 azy r 3 r r az 2 r 3 )
( x i ) = ( r , 0 , 0 )
det ( Λ j i ) = r a r ( r r ) 2
ε i j = μ i j = b b a [ δ i j 2 ar a 2 r 4 x i x j ]
ε = μ = b b a ( I 2 ar a 2 r 4 r r )
det ( ε ) = det ( μ ) = ( b b a ) 3 ( r a r ) 2
Z ij = δ 3 i δ 3 j
T ij = δ 1 i δ 1 j + δ 2 i δ 2 j
ρ i = T j i x j
( Λ j i ) = ( ρ ρ ax 2 ρ 3 axy ρ 3 0 ayx ρ 3 ρ ρ ay 2 ρ 3 0 0 0 1 )
Λ j i = ρ ρ T j i a ρ i ρ k δ i i δ kj ρ 3 + Z j i
( x i ) = ( ρ i ) = ( ρ , 0 , 0 )
det ( Λ j i ) = ρ a ρ ρ ρ
ε = μ = ρ ρ a T 2 a 2 ρ 3 ( ρ a ) ρ ρ + ( b b a ) 2 ρ a ρ Z
det ( ε ) = det ( μ ) = ( b b a ) 2 ρ a ρ
× E = B t × H = D t
E = E 0 e i ( k 0 k · x ωt ) H = 1 η 0 H 0 e i ( k 0 k · x ω t )
D = ε 0 ε E B = μ 0 μ H
k × E 0 μ H 0 = 0 k × H 0 + ε E 0 = 0
k × ( μ 1 ( k × E 0 ) ) + ε E 0 = 0
K ik ε ijk k j
( K μ 1 K + ε ) E 0 = 0
det ( K μ 1 K + ε ) = 0
det ( K n 1 K + n ) = 1 det ( n ) ( knk det ( n ) ) 2
H = f ( x ) ( knk det ( n ) )
d x d τ = H k
d k d τ = H x
( k 1 k 2 ) × n = 0
H ( k 2 ) = 0
H k · n > 0
H = 1 2 b a b ( knk det ( n ) )
H k = k 2 ar a 2 r 4 ( x · k ) x
H x = 2 ar a 2 r 4 ( x · k ) k + 3 ar 2 a 2 r 6 ( x · k ) 2 x ( b b a ) 2 ( ar a 2 r 4 ) x
H = 1 2 ρ a ρ ( knk det ( n ) )
H = 1 2 kTk 1 2 2 a 2 ρ 4 ( ρ · k ) 2 + 1 2 [ b ( ρ a ) ρ ( b a ) ] 2 ( kZk 1 )
ρ x = T
H k = Tk 2 a 2 ρ 4 ( ρ · k ) ρ + [ b ( ρ a ) ρ ( b a ) ] 2 Zk
H x = 3 2 a 2 ρ 6 ( ρ · k ) 2 ρ 2 a 2 ρ 4 ( ρ · k ) Tk + ( b b a ) 2 a 2 ρ 4 ( kZk 1 ) ρ

Metrics