OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 17, Iss. 23 — Nov. 9, 2009
  • pp: 21326–21335
« Show journal navigation

Equivalent configurations of optical 
transformation media

Tungyang Chen, Shen-Wen Cheng, and Chung-Ning Weng  »View Author Affiliations


Optics Express, Vol. 17, Issue 23, pp. 21326-21335 (2009)
http://dx.doi.org/10.1364/OE.17.021326


View Full Text Article

Acrobat PDF (235 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

We demonstrate that a medium consisting of two adjoining distinct layers of transformation materials, corresponding respectively to two linear coordinate transformations, can behave effectively as that of the same region transformed by another linear transformation. The equivalence means that, irrespective of the direction of incident wave, the fields of the medium exterior to the transformed regions of the two configurations are exactly the same. This property can also apply to a domain that is transformed by a piecewise linear transformation function, and to a medium that is mapped by a general curved function. This proof is shown analytically based on a rigorous Fourier-Bessel analysis. The equivalence suggests that, for a given transformed domain, one can find an infinite number of complementary media that altogether can give a desired effective response of certain transformation path.

© 2009 OSA

1. Introduction

The subject of transformation optics is fascinating with far-reaching implications. Transformation optics shows that the space for light can be controlled in an almost arbitrary way, in which the idea is to exploit coordinate transformations to yield an “exotic” optical path. The key feature to achieve this effect is that, with appropriately designed material parameters, the governing Maxwell’s equations remain the same under curvilinear coordinate transformation. Various concepts were proposed recently which offer new possibilities to design and engineer optical devices through the use of metamaterials. Among the latest developments, Pendry et al. [1

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science 312(5781), 1780–1782 ( 2006). [CrossRef] [PubMed]

] and Leonhardt [2

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

] showed in 2006 the possibility of designing an invisibility cloak, which could bend light around an object and return to their original propagation direction. Among relevant developments, we mention the theoretical work [3

G. W. Milton, M. Briane, and J. R. Willis, “On cloaking for elasticity and physical equations with a transformation invariant form,” N. J. Phys. 8(10), 248 ( 2006). [CrossRef]

] which provided a rigorous proof on the possibility of form invariance of various physical phenomena, and the paper [4

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 314(5801), 977–980 ( 2006). [CrossRef] [PubMed]

] for experimental demonstrations of cloak operating at microwave frequency. Other related works include the study of cloaks for different physical phenomena and various geometric shapes [5

S. A. Cummer and D. Schurig, “One path to acoustic cloaking,” N. J. Phys. 9(3), 45 ( 2007). [CrossRef]

12

T. Chen, C. N. Weng, and J. S. Chen, “Cloak for curvilinearly anisotropic media in conduction,” Appl. Phys. Lett. 93(11), 114103 ( 2008). [CrossRef]

], see also the reviews [13

U. Leonhardt and T. G. Philbin, “Transformation optics and the geometry of light,” Prog. Opt. 53, 70–152 ( 2009).

,14

V. M. Shalaev, “Physics. Transforming light,” Science 322(5900), 384–386 ( 2008). [CrossRef] [PubMed]

] for a general scope. In addition to cloaking, which expands a point into a finite domain, we mention the work [15

H. Chen and C. T. Chan, “Transformation media that rotate electromagnetic fields,” Appl. Phys. Lett. 90(24), 241105 ( 2007). [CrossRef]

] which demonstrates how to manipulate the direction of electromagnetic waves by introducing a rotation mapping of coordinates. This was later referred to as field rotators. Another recent progress [16

Z. L. Mei and T. J. Cui, “Arbitrary bending of electromagnetic waves using isotropic materials,” J. Appl. Phys. 105(10), 104913 ( 2009). [CrossRef]

] is to employ the theory of geometrical optics in demonstrating bending of electromagnetic waves. They showed that the designed wave bending structures could be isotropic. In addition, we mention an important notion of folded geometry [17

G. W. Milton, N.-A. P. Nicorovici, R. C. McPhedran, K. Cherednichenko, and Z. Jacob, “Solutions in folded geometries, and associated cloaking due to anomalous resonance,” N. J. Phys. 10(11), 115021 ( 2008). [CrossRef]

,18

U. Leonhardt and T. G. Philbin, “General relativity in electrical engineering,” N. J. Phys. 8(10), 247 ( 2006). [CrossRef]

], which will lead to materials with negatively refracting index [19

V. G. Veselago, “The electrodynamics of substances with simultaneously negative values of ε and μ,” Sov. Phys. Usp. 10 509–14 (Engl. Transl.) (1968).

]. The use of negatively refracting materials enables applications in the design of perfect electric conductor (PEC) reshaper [20

H. Chen, X. Zhang, X. Luo, H. Ma, and C. T. Chan, “Reshaping the perfect electrical conductor cylinder arbitrarily,” N. J. Phys. 10(11), 113016 ( 2008). [CrossRef]

], superscatterers [21

T. Yang, H. Chen, X. Luo, and H. Ma, “Superscatterer: enhancement of scattering with complementary media,” Opt. Express 16(22), 18545–18550 ( 2008). [CrossRef] [PubMed]

] and light concentrators [22

M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, D. R. Smith, and J. B. Pendry, “Design of electromagnetic cloaks and concentrators using form-invariant coordinate transfromations of Maxwell’s equations,” Photon. Nanostruct.: Fundam. Appl. 6(1), 87–95 ( 2008). [CrossRef]

,23

Y. Luo, H. Chen, J. Zhang, L. Ran, and J. A. Kong, “Design and analytical full-wave validation of the invisibility cloaks, concentrators, and field rotators created with a general class of transformations,” Phys. Rev. B 77(12), 125127 ( 2008). [CrossRef]

]. Anti-cloaks [24

H. Chen, X. Luo, H. Ma, and C. T. Chan, “The anti-cloak,” Opt. Express 16(19), 14603–14608 ( 2008). [CrossRef] [PubMed]

,25

G. Castaldi, I. Gallina, V. Galdi, A. Alù, and N. Engheta, “Cloak/anti-cloak interactions,” Opt. Express 17(5), 3101–3114 ( 2009). [CrossRef] [PubMed]

], which also invoke the concept of folded-geometries, can annihilate some effects of invisibility. Lai et al. [26

Y. Lai, J. Ng, H. Chen, D. Han, J. Xiao, Z. Q. Zhang, and C. T. Chan, “Illusion optics: The optical transformation of an object into another object,” Phys. Rev. Lett. 102(25), 253902 ( 2009). [CrossRef] [PubMed]

] recently outlined a mathematical formalism which showed that an object of arbitrary shape may appear to be like another object of other shape, see also Pendry [27

J. Pendry, “All smoke and metamaterials,” Nature 460(7255), 579–580 ( 2009). [CrossRef]

].

In this work, we show that two adjoining distinct pieces of metamaterials, transformed respectively from two linear coordinate transformations, can behave effectively as one single piece of metamaterial corresponding to another linear transformation. The equivalence means that, irrespective of the direction of incident wave, the electric and magnetic fields outside the transformed domains are exactly the same for the two configurations. The procedure, based on the method of transformation optics, is to first set up a linear functional relationship between the virtual space and the physical space of the two layers. This will give the prescription of the material parameters in the physical domain. As in invisibility cloaks, the constituent material parameters of the two regions are anisotropic and position varying. Next, we consider the same domain in the physical space mapped from the virtual space by one single linear transformation. In the analysis, we consider that the domains outside the transformed region are free space. We restrict our analysis to two-dimensional cylindrical configuration for simplicity, and employ the cylindrical wave expansion method to show analytically the equivalence of the two configurations. We mention that the proof can be extended to a domain that is constituted by multiple layers corresponding to a piecewise linear transformation function, and also to a medium that is mapped from a general curved function. In addition we show that the equivalence property remains valid when the core region is a PEC.

This finding implies that, for a given path of linear transformed domain, there are infinite pairs of linear coordinate transformations that will give the same field outside the transformed region. One may consider that the two transformation media as one pair of regions in which one medium is complementary to the other. The fact of equivalence, offering us new insight of transformation optics, could be of theoretical interest, relating to neutral inclusions and effective media of heterogeneous and composite materials [28

G. W. Milton, The Theory of Composites (Cambridge University Press, Cambridge, 2002).

]. Also, identifications of equivalent paths of different geometric configurations may find practical applications, for example in establishing the effective response of interactions between cloak and anti-cloak and other novel optical devices.

2. Two layers by linear coordinate transformations

We consider two adjoining layers in the physical space x that are mapped from a virtual space x via coordinate transformations. Specifically, we consider one linear transformation [Fig. 1(a) ] that maps the circular annulus with the inner radius r=a onto r= a, and the outer radius r= r0 onto r= r0; and also another linear transformation that maps the layer with inner radius r= r0 onto r= r0, and outer radius r=b onto r= b. This mapping, illustrated in Fig. 1(a), can be written as
r= α1r+ β1= r0 a r0ar+ a r0a r0 r0a, ar r0, θ=θ, z=z,
(1)
r= α2r+ β2= b r0 b r0r+ b r0 b r0 b r0, r0rb, θ=θ, z=z.
(2)
The mapping function associated with the vertical axis, the physical r space, needs to be single-valued and thus a r0 b. While that associated with the horizontal axis, the virtual space, could be multiple-valued (folded). In other words, albeit b>a, the value of r0 could be either less or greater than that of a. Specifically, the linear functions, (1) and (2), transforms the circular layer arbin the virtual space onto the circular domain ar b in the physical space. Milton et al. [3

G. W. Milton, M. Briane, and J. R. Willis, “On cloaking for elasticity and physical equations with a transformation invariant form,” N. J. Phys. 8(10), 248 ( 2006). [CrossRef]

] and earlier work [29

E. J. Post, Formal Structure of Electromagnetics: General Covariance and Electromagnetics, (North-Holland, Amsterdam, 1962).

] showed that Maxwell’s equations under coordinate transformation are invariant in their forms provided that the basic material parameters follow ε= Aε AT | det( A)| 1, where A= ( r, θ, z)/ ( r,θ,z). The transformed electric permittivity in the physical space x can be derived as [3

G. W. Milton, M. Briane, and J. R. Willis, “On cloaking for elasticity and physical equations with a transformation invariant form,” N. J. Phys. 8(10), 248 ( 2006). [CrossRef]

]
ε 1r= r β1 r ε0, ε 1θ= r r β1 ε0, ε 1z= ( 1 α1)2 r β1 r ε0, for   a r r0,
(3)
ε 2r= r β2 r ε0, ε 2θ= r r β2 ε0, ε 2z= ( 1 α2)2 r β2 r ε0, for   r0 r b.
(4)
The transformed magnetic permeability takes the same expressions as in Eqs. (3) and (4) except that ε0 is now replaced by μ0. Here ε0 and μ0 are vacuum permittivity and permeability. We mention that the transformed material parameters are anisotropic and position varying. Further one can verify that the impedance is perfectly matched at r= r0. For convenience, in the following exposition, we shall denote the regions r a, a r r0, r0 r b, r b, respectively, as regions I, II, III and IV. We note that when the point P is inside the region ABCD [Fig. 1(a)], then the transformed material parameters of the two regions II and III are positive. In the case where the point P locates in the region EFBC, the transformed material parameters for the region II will be negative (equivalent to negatively refracting material), while those corresponding to the region III are positive. This can be seen geometrically from Fig. 1(a). We note from Eqs. (1) and (2) that the value α1 represents the slope of the line segment CP¯, while the value of β1 is the ordinate of line CP¯ joining the vertical axis. From Eq. (3), the sign of the transformed parameters in region II, in which a r r0, is the same as that of ( r β1). Thus, in the case of point P lying in the region EFBC, the transformed parameters in II corresponding to PC¯are always negative. Similar observations hold for the line segment AP¯. When the point P locates inside the region GADH, the transformed parameters associated with region III are negative. We conclude that the material parameters of the two transformed domains (regions II and III) can be composed either of two positively refracting materials, or of one positively and one negatively refracting materials.

Fig. 1 (a) An illustration of the coordinate transformation of Eqs. (1) and (2) in which ( r0, r0)= ( 0.12,0.20), a=0.07, a=0.1, b=0.25, b=0.3. (b) the value of εr/ ε0 versus r. (c) the value of εθ/ ε0 versus r. (d) the value of εz/ ε0 versus r.

We first consider that the regions inside r a and outside r b are free space and thus their corresponding material parameters are simply ε0and μ0. A transverse-electric (TE) polarized EM incident field together with harmonic time dependence exp (iωt) is considered. The coordinate system is the cylindrical coordinate coinciding with the cylinder we consider. The wave equation of Ez in each region is governed by
1 εz 1r r ( r μθ Ez r)+ 1 εz 1 r2 ( 1 μr 2 Ez θ2)+ k02 Ez=0,
(5)
where k0= ω/ c0=ω μ0 ε0is the wave number of the EM wave in vacuum. The electric fields in each domain can be expressed as [30

Z. Ruan, M. Yan, C. W. Neff, and M. Qiu, “Ideal cylindrical cloak: perfect but sensitive to tiny perturbations,” Phys. Rev. Lett. 99(11), 113903 ( 2007). [CrossRef] [PubMed]

]
EzI= n dnI Jn ( k0 r)exp ( inθ), for   r a
(6)
Ez II= n { dn II Jn ( k1( r β1))+ fn II Hn (1) ( k1( r β1))}exp ( inθ), for   a r r0,
(7)
Ez III= n { dn III Jn ( k2( r β2))+ fn III Hn (1) ( k2( r β2))}exp ( inθ), for   r0 r b,
(8)
Ez IV= n { dn in Jn ( k0 r)+ fn sc Hn (1) ( k0 r)}exp ( inθ), for   r b
(9)
where k1= k0 ( r0a)/ ( r0 a) and k2= k0 ( b r0)/ ( b r0). Jnand Hn (1) are the n-order Bessel and Hankel functions of the 1st kind. The coefficients dn in correspond to the prescribed incident field and dnI, dn II, dn III, fnI, fn II, fn III and the scattering coefficient fn sc are unknown coefficients to be determined from the interface conditions at the interfaces. The interface conditions suggest that the tangential fields Ez and Hθ= 1 iω μθ Ez r be continuous at r= b, r= r0 and r= a. By virtue of orthogonality of exp ( inθ), we find that the matching conditions at r=bare
dn in Jn ( k0 b)+ fn sc Hn (1) ( k0 b)= dn III Jn ( k2 ( b β2))+ fn III Hn (1) ( k2 ( b β2)),
(10)
k0 dn in Jn ( k0 b)+ k0 fn sc Hn (1) ( k0 b) = k2 μ 2θ ( b) dn III Jn ( k2 ( b β2))+ k2 μ 2θ ( b) fn III Hn (1) ( k2 ( b β2)).
(11)
where the prime denotes differentiation with respect to the entire argument. Since
k2( b β2)= b r0 b r0 k0 bb r0b b r0= k0b,
(12)
k2 μ 2θ( b)= b r0 b r0 k0 b β2 b= b r0 b r0 k0 b bb r0b b r0= k0 bb,
(13)
Equations (10) and (11) simplify to
dn in Jn ( k0 b)+ fn sc Hn (1) ( k0 b)= dn III Jn ( k0b)+ fn III Hn (1) ( k0b),
(14)
b dn in Jn ( k0 b)+ b fn sc Hn (1) ( k0 b)=b dn III Jn ( k0b)+b fn III Hn (1) ( k0b).
(15)
At the interface of r= r0, we have
dn III Jn ( k2 ( r0 β2))+ fn III Hn (1) ( k2 ( r0 β2)) = dn II Jn ( k1 ( r0 β1))+ fn II Hn (1) ( k1 ( r0 β1)),
(16)
k2 μ 2θ( r0) dn III Jn ( k2 ( r0 β2))+ k2 μ 2θ( r0) fn III Hn (1) ( k2 ( r0 β2)) = k1 μ 1θ ( r0) dn II Jn ( k1 ( r0 β1))+ k1 μ 1θ ( r0) fn II Hn (1) ( k1 ( r0 β1)).
(17)
Again since
k2( r0 β2)= b r0 b r0 k0 b r0 r0 r0 b r0= k0 r0,
(18)
k1( r0 β1)= r0a r0 a k0 r0 r0 a r0 r0a= k0 r0,
(19)
and also
k2 μ 2θ( r0)= b r0 b r0 k0 r0 β2 r0= b r0 b r0 k0 r0 b r0 r0 r0 b r0= k0 r0 r0,
(20)
k1 μ 1θ( r0)= r0a r0 a k0 r0 β1 r0= r0a r0 a k0 r0 r0 r0 a r0 r0a= k0 r0 r0,
(21)
Equations (16) and (17) are reduced to the simple forms
( dn III dn II) Jn ( k0 r0)+ ( fn III fn II) Hn (1) ( k0 r0)=0,
(22)
( dn III dn II) Jn ( k0 r0)+ ( fn III fn II) Hn (1) ( k0 r0)=0,
(23)
which suggest that
dn II= dn III, fn II= fn III.
(24)
Lastly, at the interface r= a, we have

dn II Jn ( k0a)+ fn II Hn (1) ( k0a)= dnI Jn ( k0 a),
(25)
a dn II Jn ( k0a)+a fn II Hn (1) ( k0a)= a dnI Jn ( k0 a).
(26)

Now for a given incident field, Eqs. (14), (15), (24), (25), (26) altogether, allow us to determine the unknown coefficients dnI, dn II, dn III, fn II, fn III and fn sc, and thus the electric field as well as the magnetic field in each region can be determined readily. We shall not proceed to solve the linear algebraic system. Rather we will demonstrate that a different configuration will give the same fields in the regions I and IV (Fig. 2 ).

Fig. 2 A schematic illustration of equivalence of two configurations obtained from linear coordinate transformations. The equivalence means the electromagnetic fields in regions I and IV are identical for both configurations.

Let us now consider another linear coordinate transformation in which the circular annulus arb in the virtual space is mapped onto the region a r b. The domains in the regions I and IV are again considered as free space. This linear transformation is written as [red line of Fig. 1(a)]
r= α3r+ β3= b a bar+ ab ba ba, θ=θ, z=z,
(27)
and the corresponding material parameters for this transformed domain are
ε 3r= r β3 r ε0, ε 3θ= r r β3 ε0, ε 3z= ( 1 α3)2 r β3 r ε0.
(28)
Here we use a subscript ‘3′ to distinguish the notations in Eqs. (3) and (4). Numerical illustrations of the material parameters for Eqs. (3), (4) and (28) for a selected point of ( r0, r0) are exemplified in Fig. 1(b)1(d). Analogous to (6)-(9), the electric fields in each domain can be written as
E^zI= n d^nI Jn ( k0 r)exp ( inθ), for   r a,
(29)
E^z II+III= n { d^n II+III Jn ( k3( r β3))+ f^n II+III Hn (1) ( k3( r β3))}exp ( inθ), for   a r b,
(30)
E^z IV= n { dn in Jn ( k0 r)+ f^n sc Hn (1) ( k0 r)}exp ( inθ), for   r b,
(31)
where k3= k0 ( ba)/ ( b a) and the ‘hat’ quantities are used to distinguish the fields with those in Eqs. (6)-(9). The objective now is to demonstrate that, under the same prescribed incident field dn in, the electric field Eqs. (29)-(31), and also the magnetic field, inside the regions I and IV will be the same with those of the configuration outlined in Eqs. (6)-(9) (Fig. 2). To show this, we again demand the interface continuity conditions at r= a and r= b. This will give
dn in Jn ( k0 b)+ f^n sc Hn (1) ( k0 b)= d^n II+III Jn ( k0b)+ f^n II+III Hn (1) ( k0b),
(32)
b dn in Jn ( k0 b)+ b f^n sc Hn (1) ( k0 b)=b d^n II+III Jn ( k0b)+b f^n II+III Hn (1) ( k0b),
(33)
and
d^n II+III Jn ( k0a)+ f^n II+III Hn (1) ( k0a)= d^nI Jn ( k0 a),
(34)
a d^n II+III Jn ( k0a)+a f^n II+IIII Hn (1) ( k0a)= a d^nI Jn ( k0 a).
(35)
Now in comparison of the interface conditions of the two configurations, that is Eqs. (14), (15), (24), (25), (26) and Eqs. (32)-(35), we find interestingly that
fn sc= f^n sc, dnI= d^nI, dn II= dn III= d^n II+III, fn II= fn III= f^n II+III.
(36)
This suggests that, in the innermost region I ( r a)and also the outermost region IV ( r b), the electric and the magnetic fields of the two configurations are identical
EzI= E^zI, Ez IV= E^z IV, HθI= H^θI, Hθ IV= H^θ IV.
(37)
This observation means that the transformed media described by the two linear coordinate transformations [black lines, Fig. 1(a)] are effectively equivalent to that by one linear transformed domain [red line, Fig. 1(a)], irrespective of the direction of incident wave. The two adjoining transformed media (regions II and III) can be regarded as two-phase domains in which one medium is complementary to the other. The equivalence of both configurations is akin to the effect of neutral inclusion and effective medium of heterogeneous and composite media [28

G. W. Milton, The Theory of Composites (Cambridge University Press, Cambridge, 2002).

].

3. Some simplifications and generalizations

The proof of equivalence can be extended to transformed media that are described by a multiple piecewise linear transformation with continued end points [Fig. 3(a) ]. This can be perceived geometrically. First, one can replace any two linear transformations with one linear transformation. Repeating the process, the transformation medium can be eventually equivalent to that of a single linear transformation that links the end points. Note that in Fig. 3, the piecewise linear transformations could also be composed of a few line segments that may cross the equivalent line segment (red line). To see this, we can identify the effect of the first two line segments with that of the dashed line, and then the dashed line together with the third line segment will lead to the net effect of the red line. In fact, for a general transformation characterized by a curved function rather than a straight linear function (Fig. 3, right panel), the equivalence is still valid. Note that, however, the mapping function along with the r axis needs to be single-valued since rcorresponds to the physical space. The effective transformation of a curved function is then equivalent to one linear transformation function linked by the two ends of the curve. The reason is simple. We can discretize the curve into many line segments with sufficiently small length and thus, in principle, it is possible to utilize one single linear transformation to mimic a transformed region corresponding to a curved one.

Fig. 3 (left panel) A schematic illustration of equivalence of a multiple piecewise linear transformation to a single linear coordinate transformation (red line). (right panel) equivalence of a general curved transformation function to a single linear coordinate transformation.

Suppose now that the region I of Fig. 2 is a PEC rather than free space. In this case, Eqs. (25) and (26) simplify to
dn II Jn ( k0a)+ fn II Hn (1) ( k0a)=0.
(38)
Thus, Eqs. (14), (15), (24) and (38) can be used to determine the 5 unknown coefficients dn II, dn III, fn II, fn IIIand fn sc. Let us now consider the configuration of the right panel of Fig. 2. Again, the region I is a PEC. Instead of (34) and (35), we now have
d^n II+III Jn ( k0a)+ f^n II+III Hn (1) ( k0a)=0.
(39)
Hence the electric field of this configuration can be resolved from (32), (33) and (39). Now in comparison of these two linear algebraic systems, we find again that
fn sc= f^n sc, dn II= dn III= d^n II+III, fn II= fn III= f^n II+III.
(40)
This suggests that when the region I is a PEC, the electromagnetic fields in regions I and IV in both configurations remains the same.

Lastly we consider the simple situation in which a=aand b= b. In this case, the interface conditions of Eqs. (14) and (15) simplify to
( dn in dn III) Jn ( k0b)+ ( fn sc fn III) Hn (1) ( k0b)=0,
(41)
( dn in dn III) Jn ( k0b)+ ( fn sc dn III) Hn (1) ( k0b)=0,
(42)
at r= b, and those of Eqs. (25) and (26) to
( dn II dnI) Jn ( k0a)+ ( fn II) Hn (1) ( k0a)=0,
(43)
( dn II dnI) Jn ( k0a)+ ( fn II) Hn (1) ( k0a)=0,
(44)
at r= a. Note that Eq. (24) remains unchanged. All these together allow us to solve the unknown coefficients as
dn in= dn II= dn III, fn sc= fn II= fn III=0.
(45)
This implies that there is no scattering field and that the electric fields in regions I and IV are exactly identical. We mention, however, that inside the region II and III, the field solutions are not the same.

4. Discussions and implications

In summary, we have shown analytically in this work, by resolving the electromagnetic field solutions of the two different configurations, that a medium consisting of two adjoining distinct layers of transformation materials, corresponding respectively to two linear coordinate transformations, can behave effectively as that of the same region transformed by another linear transformation. We note that as the r axis could be folded, the point P( r0, r0) may lie on any point inside the semi-infinite strip of EFGH. Figure 4 depicts a schematic illustration of possible equivalent transformation paths that give to the same effect. In fact, different transformation paths may correspond to different optical behavior. For example, the path indicated by the blue dashed line (Fig. 4) could simulate the effects of anti-cloaks [24

H. Chen, X. Luo, H. Ma, and C. T. Chan, “The anti-cloak,” Opt. Express 16(19), 14603–14608 ( 2008). [CrossRef] [PubMed]

], while others (green and black dashed lines) can be modeled as concentrators.

Fig. 4 Possible transformation paths that are equivalent to the linear one (red).

It should be noted that the transformed material properties (3) and (4) are derived based on cylindrical configurations, which in fact involve the geometric information of the adjoining layers. We conjecture that the equivalence of the configurations could also be derivable from Fermat’s principle based on extremal optical path [31

V. G. Veselago, “About the wording of Fermat’s principle for light propagation in media with negative refraction index,” (2002) http://arxiv.org/abs/cond-mat/0203451.

,32

M. Born, and E. Wolf, Principles of Optics, 7th (Expanded) Edition, (Cambridge University Press, New York, 2002).

]. In a separate study we have performed a numerical check on the equivalence property based on finite element calculations. The equivalence between two configurations offers a new insight at the fundamentals of transformation optics. For instance, given a transformed domain, one can always find an infinite number of complementary media that altogether can give an effective response of certain transformed path. Identifications of equivalent paths of different geometric configurations may find useful in some practical applications, for example in establishing the effective response of combined transformed domains. For potential applications, when a certain desired optical ray-path outcome would require a medium which is realistically unfeasible one could look for a series of adjacent media which are realistically feasible and lead to the same net effect. Particularly, with the recent advances of fabrication methods of metamaterials, one can design synthetic structures with a desired distribution of electric permittivity and magnetic permeability for practical implementations.

In a future study, as in the design of cloaks of arbitrary cross-sections [9

A. Nicolet, F. Zolla, and S. Guenneau, “Electromagnetic analysis of cylindrical cloaks of an arbitrary cross section,” Opt. Express 33, 1584–1586 ( 2008).

], one can explore whether the equivalence property applies to other geometrical configurations. To analytically prove the equivalence, one may first think of extending the proof to three-dimensional spherical regions and also to regions with elliptical geometry. For spherical regions, one can derive the full wave solution based on Mie scattering [33

H. Chen, B. I. Wu, B. Zhang, and J. A. Kong, “Electromagnetic wave interactions with a metamaterial cloak,” Phys. Rev. Lett. 99(6), 063903 ( 2007). [CrossRef] [PubMed]

]. Similar electromagnetic scattering of elliptical layers, in terms of series expansions of Mathieu functions [34

S. Caorsi, M. Pastorino, and M. Raffetto, “Electromagnetic scattering by a multilayer elliptic cylinder under transverse-magnetic illumination: series solution in terms of Mathieu functions,” IEEE Trans. Antenn. Propag. 45(6), 926–935 ( 1997). [CrossRef]

], can be envisaged based on elliptical coordinates. More generally, one may wonder if the equivalence property holds for any configurations that can be characterized by orthogonal coordinates [10

H. Chen, “Transformation optics in orthogonal coordinates,” J. Opt. A, Pure Appl. Opt. 11(7), 075102 ( 2009). [CrossRef]

], in which a full wave series solution could be analyzed in closed forms.

Acknowledgment

This work was supported by the National Science Council, Taiwan, under contract NSC97-2211-E-006-117-MY3.

References and links

1.

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science 312(5781), 1780–1782 ( 2006). [CrossRef] [PubMed]

2.

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

3.

G. W. Milton, M. Briane, and J. R. Willis, “On cloaking for elasticity and physical equations with a transformation invariant form,” N. J. Phys. 8(10), 248 ( 2006). [CrossRef]

4.

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 314(5801), 977–980 ( 2006). [CrossRef] [PubMed]

5.

S. A. Cummer and D. Schurig, “One path to acoustic cloaking,” N. J. Phys. 9(3), 45 ( 2007). [CrossRef]

6.

H. Chen and C. T. Chan, “Acoustic cloaking in three dimensions using acoustic metamaterials,” Appl. Phys. Lett. 91(18), 183518 ( 2007). [CrossRef]

7.

G. X. Yu, W. X. Jiang, and T. J. Cui, “Invisible slab cloaks via embedded optical transformation,” Appl. Phys. Lett. 94(4), 041904 ( 2009). [CrossRef]

8.

J. Hu, X. Zhou, and G. Hu, “Design method for electromagnetic cloak with arbitrary shapes based on Laplace’s equation,” Opt. Express 17(3), 1308–1320 ( 2009). [CrossRef] [PubMed]

9.

A. Nicolet, F. Zolla, and S. Guenneau, “Electromagnetic analysis of cylindrical cloaks of an arbitrary cross section,” Opt. Express 33, 1584–1586 ( 2008).

10.

H. Chen, “Transformation optics in orthogonal coordinates,” J. Opt. A, Pure Appl. Opt. 11(7), 075102 ( 2009). [CrossRef]

11.

N. Kundtz, D. A. Roberts, J. Allen, S. Cummer, and D. R. Smith, “Optical source transformations,” Opt. Express 16(26), 21215–21222 ( 2008). [CrossRef] [PubMed]

12.

T. Chen, C. N. Weng, and J. S. Chen, “Cloak for curvilinearly anisotropic media in conduction,” Appl. Phys. Lett. 93(11), 114103 ( 2008). [CrossRef]

13.

U. Leonhardt and T. G. Philbin, “Transformation optics and the geometry of light,” Prog. Opt. 53, 70–152 ( 2009).

14.

V. M. Shalaev, “Physics. Transforming light,” Science 322(5900), 384–386 ( 2008). [CrossRef] [PubMed]

15.

H. Chen and C. T. Chan, “Transformation media that rotate electromagnetic fields,” Appl. Phys. Lett. 90(24), 241105 ( 2007). [CrossRef]

16.

Z. L. Mei and T. J. Cui, “Arbitrary bending of electromagnetic waves using isotropic materials,” J. Appl. Phys. 105(10), 104913 ( 2009). [CrossRef]

17.

G. W. Milton, N.-A. P. Nicorovici, R. C. McPhedran, K. Cherednichenko, and Z. Jacob, “Solutions in folded geometries, and associated cloaking due to anomalous resonance,” N. J. Phys. 10(11), 115021 ( 2008). [CrossRef]

18.

U. Leonhardt and T. G. Philbin, “General relativity in electrical engineering,” N. J. Phys. 8(10), 247 ( 2006). [CrossRef]

19.

V. G. Veselago, “The electrodynamics of substances with simultaneously negative values of ε and μ,” Sov. Phys. Usp. 10 509–14 (Engl. Transl.) (1968).

20.

H. Chen, X. Zhang, X. Luo, H. Ma, and C. T. Chan, “Reshaping the perfect electrical conductor cylinder arbitrarily,” N. J. Phys. 10(11), 113016 ( 2008). [CrossRef]

21.

T. Yang, H. Chen, X. Luo, and H. Ma, “Superscatterer: enhancement of scattering with complementary media,” Opt. Express 16(22), 18545–18550 ( 2008). [CrossRef] [PubMed]

22.

M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, D. R. Smith, and J. B. Pendry, “Design of electromagnetic cloaks and concentrators using form-invariant coordinate transfromations of Maxwell’s equations,” Photon. Nanostruct.: Fundam. Appl. 6(1), 87–95 ( 2008). [CrossRef]

23.

Y. Luo, H. Chen, J. Zhang, L. Ran, and J. A. Kong, “Design and analytical full-wave validation of the invisibility cloaks, concentrators, and field rotators created with a general class of transformations,” Phys. Rev. B 77(12), 125127 ( 2008). [CrossRef]

24.

H. Chen, X. Luo, H. Ma, and C. T. Chan, “The anti-cloak,” Opt. Express 16(19), 14603–14608 ( 2008). [CrossRef] [PubMed]

25.

G. Castaldi, I. Gallina, V. Galdi, A. Alù, and N. Engheta, “Cloak/anti-cloak interactions,” Opt. Express 17(5), 3101–3114 ( 2009). [CrossRef] [PubMed]

26.

Y. Lai, J. Ng, H. Chen, D. Han, J. Xiao, Z. Q. Zhang, and C. T. Chan, “Illusion optics: The optical transformation of an object into another object,” Phys. Rev. Lett. 102(25), 253902 ( 2009). [CrossRef] [PubMed]

27.

J. Pendry, “All smoke and metamaterials,” Nature 460(7255), 579–580 ( 2009). [CrossRef]

28.

G. W. Milton, The Theory of Composites (Cambridge University Press, Cambridge, 2002).

29.

E. J. Post, Formal Structure of Electromagnetics: General Covariance and Electromagnetics, (North-Holland, Amsterdam, 1962).

30.

Z. Ruan, M. Yan, C. W. Neff, and M. Qiu, “Ideal cylindrical cloak: perfect but sensitive to tiny perturbations,” Phys. Rev. Lett. 99(11), 113903 ( 2007). [CrossRef] [PubMed]

31.

V. G. Veselago, “About the wording of Fermat’s principle for light propagation in media with negative refraction index,” (2002) http://arxiv.org/abs/cond-mat/0203451.

32.

M. Born, and E. Wolf, Principles of Optics, 7th (Expanded) Edition, (Cambridge University Press, New York, 2002).

33.

H. Chen, B. I. Wu, B. Zhang, and J. A. Kong, “Electromagnetic wave interactions with a metamaterial cloak,” Phys. Rev. Lett. 99(6), 063903 ( 2007). [CrossRef] [PubMed]

34.

S. Caorsi, M. Pastorino, and M. Raffetto, “Electromagnetic scattering by a multilayer elliptic cylinder under transverse-magnetic illumination: series solution in terms of Mathieu functions,” IEEE Trans. Antenn. Propag. 45(6), 926–935 ( 1997). [CrossRef]

OCIS Codes
(160.1190) Materials : Anisotropic optical materials
(260.2110) Physical optics : Electromagnetic optics
(160.2710) Materials : Inhomogeneous optical media

ToC Category:
Physical Optics

History
Original Manuscript: August 25, 2009
Revised Manuscript: October 21, 2009
Manuscript Accepted: November 2, 2009
Published: November 6, 2009

Citation
Tungyang Chen, Shen-Wen Cheng, and Chung-Ning Weng, "Equivalent configurations of optical 
transformation media," Opt. Express 17, 21326-21335 (2009)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-17-23-21326


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science 312(5781), 1780–1782 (2006). [CrossRef] [PubMed]
  2. U. Leonhardt, “Optical conformal mapping,” Science 312(5781), 1777–1780 (2006). [CrossRef] [PubMed]
  3. G. W. Milton, M. Briane, and J. R. Willis, “On cloaking for elasticity and physical equations with a transformation invariant form,” N. J. Phys. 8(10), 248 (2006). [CrossRef]
  4. 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 314(5801), 977–980 (2006). [CrossRef] [PubMed]
  5. S. A. Cummer and D. Schurig, “One path to acoustic cloaking,” N. J. Phys. 9(3), 45 (2007). [CrossRef]
  6. H. Chen and C. T. Chan, “Acoustic cloaking in three dimensions using acoustic metamaterials,” Appl. Phys. Lett. 91(18), 183518 (2007). [CrossRef]
  7. G. X. Yu, W. X. Jiang, and T. J. Cui, “Invisible slab cloaks via embedded optical transformation,” Appl. Phys. Lett. 94(4), 041904 (2009). [CrossRef]
  8. J. Hu, X. Zhou, and G. Hu, “Design method for electromagnetic cloak with arbitrary shapes based on Laplace’s equation,” Opt. Express 17(3), 1308–1320 (2009). [CrossRef] [PubMed]
  9. A. Nicolet, F. Zolla, and S. Guenneau, “Electromagnetic analysis of cylindrical cloaks of an arbitrary cross section,” Opt. Express 33, 1584–1586 (2008).
  10. H. Chen, “Transformation optics in orthogonal coordinates,” J. Opt. A, Pure Appl. Opt. 11(7), 075102 (2009). [CrossRef]
  11. N. Kundtz, D. A. Roberts, J. Allen, S. Cummer, and D. R. Smith, “Optical source transformations,” Opt. Express 16(26), 21215–21222 (2008). [CrossRef] [PubMed]
  12. T. Chen, C. N. Weng, and J. S. Chen, “Cloak for curvilinearly anisotropic media in conduction,” Appl. Phys. Lett. 93(11), 114103 (2008). [CrossRef]
  13. U. Leonhardt and T. G. Philbin, “Transformation optics and the geometry of light,” Prog. Opt. 53, 70–152 (2009).
  14. V. M. Shalaev, “Physics. Transforming light,” Science 322(5900), 384–386 (2008). [CrossRef] [PubMed]
  15. H. Chen and C. T. Chan, “Transformation media that rotate electromagnetic fields,” Appl. Phys. Lett. 90(24), 241105 (2007). [CrossRef]
  16. Z. L. Mei and T. J. Cui, “Arbitrary bending of electromagnetic waves using isotropic materials,” J. Appl. Phys. 105(10), 104913 (2009). [CrossRef]
  17. G. W. Milton, N.-A. P. Nicorovici, R. C. McPhedran, K. Cherednichenko, and Z. Jacob, “Solutions in folded geometries, and associated cloaking due to anomalous resonance,” N. J. Phys. 10(11), 115021 (2008). [CrossRef]
  18. U. Leonhardt and T. G. Philbin, “General relativity in electrical engineering,” N. J. Phys. 8(10), 247 (2006). [CrossRef]
  19. V. G. Veselago, “The electrodynamics of substances with simultaneously negative values of ε and μ,” Sov. Phys. Usp. 10 509–14 (Engl. Transl.) (1968).
  20. H. Chen, X. Zhang, X. Luo, H. Ma, and C. T. Chan, “Reshaping the perfect electrical conductor cylinder arbitrarily,” N. J. Phys. 10(11), 113016 (2008). [CrossRef]
  21. T. Yang, H. Chen, X. Luo, and H. Ma, “Superscatterer: enhancement of scattering with complementary media,” Opt. Express 16(22), 18545–18550 (2008). [CrossRef] [PubMed]
  22. M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, D. R. Smith, and J. B. Pendry, “Design of electromagnetic cloaks and concentrators using form-invariant coordinate transfromations of Maxwell’s equations,” Photon. Nanostruct.: Fundam. Appl. 6(1), 87–95 (2008). [CrossRef]
  23. Y. Luo, H. Chen, J. Zhang, L. Ran, and J. A. Kong, “Design and analytical full-wave validation of the invisibility cloaks, concentrators, and field rotators created with a general class of transformations,” Phys. Rev. B 77(12), 125127 (2008). [CrossRef]
  24. H. Chen, X. Luo, H. Ma, and C. T. Chan, “The anti-cloak,” Opt. Express 16(19), 14603–14608 (2008). [CrossRef] [PubMed]
  25. G. Castaldi, I. Gallina, V. Galdi, A. Alù, and N. Engheta, “Cloak/anti-cloak interactions,” Opt. Express 17(5), 3101–3114 (2009). [CrossRef] [PubMed]
  26. Y. Lai, J. Ng, H. Chen, D. Han, J. Xiao, Z. Q. Zhang, and C. T. Chan, “Illusion optics: The optical transformation of an object into another object,” Phys. Rev. Lett. 102(25), 253902 (2009). [CrossRef] [PubMed]
  27. J. Pendry, “All smoke and metamaterials,” Nature 460(7255), 579–580 (2009). [CrossRef]
  28. G. W. Milton, The Theory of Composites (Cambridge University Press, Cambridge, 2002).
  29. E. J. Post, Formal Structure of Electromagnetics: General Covariance and Electromagnetics, (North-Holland, Amsterdam, 1962).
  30. Z. Ruan, M. Yan, C. W. Neff, and M. Qiu, “Ideal cylindrical cloak: perfect but sensitive to tiny perturbations,” Phys. Rev. Lett. 99(11), 113903 (2007). [CrossRef] [PubMed]
  31. V. G. Veselago, “About the wording of Fermat’s principle for light propagation in media with negative refraction index,” (2002) http://arxiv.org/abs/cond-mat/0203451 .
  32. M. Born, and E. Wolf, Principles of Optics, 7th (Expanded) Edition, (Cambridge University Press, New York, 2002).
  33. H. Chen, B. I. Wu, B. Zhang, and J. A. Kong, “Electromagnetic wave interactions with a metamaterial cloak,” Phys. Rev. Lett. 99(6), 063903 (2007). [CrossRef] [PubMed]
  34. S. Caorsi, M. Pastorino, and M. Raffetto, “Electromagnetic scattering by a multilayer elliptic cylinder under transverse-magnetic illumination: series solution in terms of Mathieu functions,” IEEE Trans. Antenn. Propag. 45(6), 926–935 (1997). [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.

Figures

Fig. 1 Fig. 2 Fig. 3
 
Fig. 4
 

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited