Revolution analysis of three-dimensional arbitrary cloaks

doi:10.1364/OE.17.022603

Revolution analysis of three-dimensional arbitrary cloaks

Guillaume Dupont, Sébastien Guenneau, Stefan Enoch, Guillaume Demesy, André Nicolet, Frédéric Zolla, and André Diatta »View Author Affiliations

Optics Express, Vol. 17, Issue 25, pp. 22603-22608 (2009)
http://dx.doi.org/10.1364/OE.17.022603

View Full Text Article

Acrobat PDF (495 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. Change of coordinates and pullbacks from optical to physical space
2.2. Singularity analysis of the transformation matrix
3. Finite elements computations
4. Conclusion
– References and links

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

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (17)
Cited By
Figures (4)
Metrics

Abstract

We extend the design of radially symmetric three-dimensional invisibility cloaks through transformation optics [1] to cloaks with a surface of revolution. We derive the expression of the transformation matrix and show that one of its eigenvalues vanishes on the inner boundary of the cloaks, while the other two remain strictly positive and bounded. The validity of our approach is confirmed by finite edge-elements computations for a non-convex cloak of varying thickness.

1. Introduction

It was recently found that transformation optics open new avenues in electromagnetic cloaking, either through their heterogeneous anisotropic effective material parameters (transformation optics, [1

J.B. Pendry, D. Shurig, 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]

]) or through low index materials [3

A. Alu and N. Engheta, “Achieving Transparency with Plasmonic and Metamaterial Coatings,” Phys. Rev. E 95 016623 ( 2005). [CrossRef]

] or negative refractive index materials [4

G. Milton and N. A. Nicorovici, “On the cloaking effects associated with anomalous localized resonance,” Proc. R. Soc. London A 462 3027 ( 2006). [CrossRef]

]. Interestingly, the invisibility is preserved in the case of an intense near field [5

F. Zolla, S. Guenneau, A. Nicolet, and J.B. Pendry, “Electromagnetic analysis of cylindrical invisibility cloaks and the mirage effect,” Opt. Lett. 32, 1069 ( 2007). [CrossRef] [PubMed]

], when the ray optics picture breaks down. The mathematics behind the scene have been known from researchers working in the area of inverse conductivity problems [6

A. Greenleaf, M. Lassas, and G. Uhlmann, “On nonuniqueness for Calderons inverse problem,” Math. Res. Lett. 10, 685–693 ( 2003).

]. The first experimental realization of an invisibility cloak, chiefly achieved in the microwave regime [7

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

], suggests that cloaking will be limited to a very narrow range of frequencies. However, it will not be perfect since the cloak is necessarily dissipative and dispersive, and some of its tensor components are singular on the inner boundary. The latter drawback can be overcome by considering a generalized transform [8

P. Zhang, Y. Jin, and S. He, “Obtaining a nonsingular two-dimensional cloak of complex shape from a perfect three-dimensional cloak,” Appl. Phys. Lett. 93, 243502 ( 2008). [CrossRef]

, 9

U. Leonhardt and T. Tyc, “Broadband invisibility by non-euclidean cloaking,” Science 323, 110 ( 2009). [CrossRef]

] in an upper dimensional space and then projecting the resulting metric on the physical space, and this leads to non-singular tensors of permittivity and permeability. Alternatively, one can design an approximate structured cloak via homogenization [10

W. Cai, U. K. Chettiar, A. V. Kildiev, and V. M. Shalaev, “Optical Cloaking with metamaterials,” Nat. Photon. 1, 224–227 ( 2007). [CrossRef]

, 11

M. Farhat, S. Guenneau, A. B. Movchan, and S. Enoch, “Achieving invisibility over a finite range of frequencies,” Opt. Express 16, 5656–5661 ( 2008). [CrossRef] [PubMed]

, 12

C.W. Qiu, L. Hu, X. Xu, and Y. Feng, “Spherical cloaking with homogeneous isotropic multilayered structures,” Phys. Rev. E 79, 047602 ( 2009). [CrossRef]

], although the rapid growth of the field, fueled by a keen interest of the optics community, promises a large panel of new technological applications.

A non trivial question to ask is whether one can design cloaks of non-spherical shapes. A parameterization of the cloak’s boundaries was proposed by three of us to design cylindrical cloaks of an arbitrary cross-section [17

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

]. Its extension to the general three-dimensional case now requires to parameterize the inner and outer boundaries of the cloak as some surfaces with varying radii

ρ (θ, ϕ) = a_{0, 0} + \underset{(m, n) \in ℕ^{2} ∖ {(0, 0)}}{Σ} {a_{m, n} \cos (m θ + n ϕ) + b_{m, n} \sin (m θ + n ϕ)},

(1)

which is somewhat out of reach running the COMSOL multiphysics package on a 256 Gb RAM computer: Compared to a spherical cloak, the tetrahedral mesh needs be further refined due to the complexity of surfaces involved. However, to demonstrate the versatility of our Fourier-based approach, it is enough to analyse arbitrary cloaks built with surfaces of revolution i.e. generated by rotating a curve about an axis.

In the present paper, we discuss the design of axially invariant three-dimensional cloaks with an arbitrary cross-section described by two functions R ₁(ϕ) and R ₂(ϕ) giving an angle dependent distance from the origin. These functions correspond respectively to the interior and exterior boundary of the cloak. We shall only assume that these two boundaries can be represented by a differentiable function. Their finite Fourier expansions are thought in the form R_j (ϕ)=a ^j _0,0 +∑ ^p _n=1 a ^j _0,n cos(nϕ), j=1,2, where p can be a small integer, and a ^j _0,n=0 for n>p for computational easiness. Note that our approach encompasses the case of toroidal cloaks [13

Y. You, G. W. Kattawar, and P. Yang, “Invisibility cloaks for toroids,” Opt. Express 17, 6591 ( 2009). [CrossRef] [PubMed]

] but cloaks with irregular boundaries such as cubes, fall beyond the scope of this study. Nevertheless, contrarily to previous three dimensional numerical studies [13

Y. You, G. W. Kattawar, and P. Yang, “Invisibility cloaks for toroids,” Opt. Express 17, 6591 ( 2009). [CrossRef] [PubMed]

, 12

C.W. Qiu, L. Hu, X. Xu, and Y. Feng, “Spherical cloaking with homogeneous isotropic multilayered structures,” Phys. Rev. E 79, 047602 ( 2009). [CrossRef]

], our scheme can be used in the design of cloaks with non-convex boundaries.

To illustrate our methodology, we compute the electromagnetic field diffracted by a cloak with rotational symmetry about the z-axis. We perform full-wave finite element simulations in the commercial package COMSOL, when the cloak is illuminated by an approximate plane wave (generated by a constant electric field on the upper surface of the computational domain).

2. Change of coordinates and pullbacks from optical to physical space

2.1. Material properties of the heterogeneous anisotropic cloak

The geometric transformation which maps the field within the full domain ρ≤R ₂(ϕ) onto the annular domain R ₁(ϕ)≤ρ′≤R ²(ϕ) can be expressed as:

{\begin{matrix} ρ' (ρ, ϕ) = R_{1} (ϕ) + ρ \frac{R_{2} (ϕ) - R_{1} (ϕ)}{R_{2} (ϕ)}, \\ θ' = θ, 0 < θ \leq 2 π, \\ ϕ' = ϕ, - π ⁄ 2 \leq π ⁄ 2 \end{matrix}

(2)

where 0≤ρ≤R ₂(ϕ). Note that the transformation maps the field for ρ>R ₂(ϕ) onto itself through the identity transformation.

This change of co-ordinates is characterized by the transformation of the differentials through the Jacobian:

J' (ρ', ϕ') = \frac{\partial (ρ (ρ', ϕ'), θ, ϕ)}{\partial (ρ', θ', ϕ')} .

(3)

This change of coordinates amounts to replacing a homogeneous isotropic medium with scalar permittivity and permeability ε and µ, by a metamaterial described by anisotropic heterogeneous matrices of permittivity and permeability given by [5

F. Zolla, S. Guenneau, A. Nicolet, and J.B. Pendry, “Electromagnetic analysis of cylindrical invisibility cloaks and the mirage effect,” Opt. Lett. 32, 1069 ( 2007). [CrossRef] [PubMed]

, 14

A. Nicolet, J.F. Remacle, B. Meys, A. Genon, and W. Legros, “Transformation methods in computational electromagnetics,” J. Appl. Phys. 75 6036 ( 1994). [CrossRef]

]

\underset{=}{ε'} = ε T^{- 1}, and \underset{=}{μ'} = μ T^{- 1},

(4)

where T=J ^T J/det(J) is a representation of the metric tensor in the so called stretched radial coordinates. Note that there is no change in the impedance of the media since the permittivity and permeability undergo the same transformation.

After some elementary algebra, we find that

T^{- 1} = (\begin{matrix} \frac{{c_{13}^{2} + ρ' (ρ', ϕ)}^{2}}{c_{11} {ρ'}^{2}} & 0 & - \frac{c_{13}}{ρ'} \\ 0 & c_{11} & 0 \\ - \frac{c_{13}}{ρ'} & 0 & c_{11} \end{matrix}),

(5)

where

c_{11} (ϕ') = \frac{R_{2} (ϕ')}{R_{2} (ϕ') - R_{1} (ϕ')},

(6)

and

c_{13} (ϕ') = R_{2} (ϕ') \frac{ρ' - R_{2} (ϕ')}{{(R_{2} (ϕ') - R_{1} (ϕ'))}^{2}} \frac{d R_{1} (ϕ')}{d ϕ'} + R_{1} (ϕ') \frac{R_{1} (ϕ') - ρ'}{{(R_{2} (ϕ') - R_{1} (ϕ'))}^{2}} \frac{d R_{2} (ϕ')}{d ϕ'},

(7)

for R ₁(ϕ′)≤ρ′≤R ₂(ϕ′). Elsewhere, T ^-1 reduces to the identity matrix (c ₁₁=1, c ₁₃=0 and ρ=ρ′ for ρ′>R ₂(ϕ′)).

2.2. Singularity analysis of the transformation matrix

To exemplify the symmetric nature of the coefficients of T ^-1, and due to the role played by its first entry in the following singularity analysis, we show the variation of (T ^-1)11 in Fig. 1. We note that this coefficient varies between zero and three. The fine mesh on the inner boundary of the cloak is also apparent from Fig. 1, as are the non-convex and symmetric features of the cloak.

In the case of a spherical cloak, c ₁₃ vanishes and T ^-1 reduces to

Diag (\frac{ρ^{2}}{c_{11} {ρ'}^{2}}, c_{11}, c_{11})

. We note that the first eigenvalue vanishes on the inner boundary and the other two remain constant, which is consistent with the singularity analysis led in [15

R. V. Kohn, H. Shen, M. S. Vogelius, and M. I. Weinstein, “Cloaking via change of variables in electric impedance tomography,” Inverse Probl. 24, 015016 ( 2008). [CrossRef]

]. This is unlike the circular cylindrical case whereby one eigenvalue goes to zero while the other one goes to infinity on the inner boundary [5

F. Zolla, S. Guenneau, A. Nicolet, and J.B. Pendry, “Electromagnetic analysis of cylindrical invisibility cloaks and the mirage effect,” Opt. Lett. 32, 1069 ( 2007). [CrossRef] [PubMed]

Fig. 1. 3D plot of (T ^-1)₁₁, as given by Eq. 5, Eq. 6 and Eq. 7 within the cloak with boundaries given by Eq. 10. The symmetries are noted.

However, in our case, the cloak is of an arbitrary shape, and it is therefore illuminating to look at the behaviour of its permittivity and permeability tensors’s eigenvalues. The eigenvalues of (5) are found to be

λ_{j} = \frac{c_{13}^{2} + ρ^{2} + c_{11}^{2} {ρ'}^{2}}{2 c_{11} {ρ'}^{2}} + \frac{{(- 1)}^{j}}{2} \sqrt{{(\frac{c_{13}^{2} + ρ^{2} + c_{11}^{2} {ρ'}^{2}}{c_{11} {ρ'}^{2}})}^{2} - 4 \frac{ρ^{2}}{{ρ'}^{2}}}, j = 1, 2, and λ_{3} = c_{11} .

(8)

We can therefore see that λ _j,j=1,2 are spatially varying functions of ρ′(ρ,ϕ), such that λ ₁=0 when ρ=0 i.e. at the inner boundary, while λ ₂>0. Importantly, λ ₃ is a strictly positive constant for a given angle ϕ′(i.e. a function independent upon ρ). Last, we checked that when there is no longer a symmetry of revolution about one axis, all three eigenvalues are also spatially varying with θ′, but λ ₃ remains independent upon ρ. This reflects the fact that the geodesics for light within the arbitrary cloak follow more and more complex trajectories when we perturb the geometry away from the spherical cloak.

3. Finite elements computations

We would now like to further investigate the electromagnetic response of the cloak to an incident plane wave from above, see Fig. 2. For this, we choose the electric field E as the unknown:

\nabla \times ({\underset{=}{μ'}}^{- 1} \nabla \times E) - k^{2} \underset{=}{ε'} E = 0

(9)

where

k = ω \sqrt{μ_{0} ε_{0}} = ω ⁄ c

is the wavenumber, c being the speed of light in vacuum, and ε̳′ and µ̳′ are defined by Eqs. (4). Also, E=E _i +E _d , where E _i is the incident field and E _d is the diffracted field which satisfies the usual outgoing wave conditions (to ensure existence and uniqueness of the solution). The weak formulation associated with Eq. 9 is discretised using second order finite edge elements (or Whitney forms) which behave nicely under geometric transforms (pull-back properties) [14

A. Nicolet, J.F. Remacle, B. Meys, A. Genon, and W. Legros, “Transformation methods in computational electromagnetics,” J. Appl. Phys. 75 6036 ( 1994). [CrossRef]

Fig. 2. 3D plot of the magnitude

\underset{=}{ε'}

of the total electric field for a plane wave of wavenumber k=2π/0.3 incident from above on an non-convex invisibility cloak.

For the sake of illustration, cf. Fig. 2, let us consider a cloak with inner and outer boundaries expressed as

R_{1} (ϕ) = 0.2 + 0.02 \cos (4 ϕ), R_{2} (ϕ) = 0.4 + 0.02 \cos (8 ϕ) .

(10)

Fig. 3. 2D plot of

\sqrt{E_{1}^{2} + E_{2}^{2} + E_{3}^{2}}

generated by a slice of Fig. 2 in the xz-plane for y=0.

We report the computations for the magnitude of the total electric field in Fig. 2 (3D plot), Fig. 3 (2D plot in the xz-plane) and Fig. 4 (2D plot in the yz-plane) for a plane wave incident from above at wavenumber k=2π/0.3 (units are in inverse of a length, say µm-1 for nearly visible light (UV)).

Around 4.10⁵ elements were used in this computation, which corresponds to about 2.7 10⁶ degrees of freedom. While the convergence of the numerical scheme has been checked by considering different types of meshes, the large size of the system means we were not able to further refine the mesh for the computational resources at hand.

Fig. 4. 2D plot of

\sqrt{E_{1}^{2} + E_{2}^{2} + E_{3}^{2}}

generated by a slice of Fig. 2 in the yz-plane for x=0.

4. Conclusion

In conclusion, we have proposed a design of an arbitrarily shaped cloak using a Fourier approach. We only assumed that the cloak displays a symmetry of revolution about one axis in order to reduce the computational complexity of the problem (extension to completely arbitrarily shaped cloaks is a straightforward matter). Cloaking has been confirmed numerically for an incident plane wave in resonance with the concealed region and an analysis of the cloak’s singularity has been carried out.

A. Diatta and S. Guenneau acknowledge funding from EPSRC grant EP/F027125/1.

References and links

1.	J.B. Pendry, D. Shurig, 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.	A. Alu and N. Engheta, “Achieving Transparency with Plasmonic and Metamaterial Coatings,” Phys. Rev. E 95 016623 ( 2005). [CrossRef]
4.	G. Milton and N. A. Nicorovici, “On the cloaking effects associated with anomalous localized resonance,” Proc. R. Soc. London A 462 3027 ( 2006). [CrossRef]
5.	F. Zolla, S. Guenneau, A. Nicolet, and J.B. Pendry, “Electromagnetic analysis of cylindrical invisibility cloaks and the mirage effect,” Opt. Lett. 32, 1069 ( 2007). [CrossRef] [PubMed]
6.	A. Greenleaf, M. Lassas, and G. Uhlmann, “On nonuniqueness for Calderons inverse problem,” Math. Res. Lett. 10, 685–693 ( 2003).
7.	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, 977 ( 2006). [CrossRef] [PubMed]
8.	P. Zhang, Y. Jin, and S. He, “Obtaining a nonsingular two-dimensional cloak of complex shape from a perfect three-dimensional cloak,” Appl. Phys. Lett. 93, 243502 ( 2008). [CrossRef]
9.	U. Leonhardt and T. Tyc, “Broadband invisibility by non-euclidean cloaking,” Science 323, 110 ( 2009). [CrossRef]
10.	W. Cai, U. K. Chettiar, A. V. Kildiev, and V. M. Shalaev, “Optical Cloaking with metamaterials,” Nat. Photon. 1, 224–227 ( 2007). [CrossRef]
11.	M. Farhat, S. Guenneau, A. B. Movchan, and S. Enoch, “Achieving invisibility over a finite range of frequencies,” Opt. Express 16, 5656–5661 ( 2008). [CrossRef] [PubMed]
12.	C.W. Qiu, L. Hu, X. Xu, and Y. Feng, “Spherical cloaking with homogeneous isotropic multilayered structures,” Phys. Rev. E 79, 047602 ( 2009). [CrossRef]
13.	Y. You, G. W. Kattawar, and P. Yang, “Invisibility cloaks for toroids,” Opt. Express 17, 6591 ( 2009). [CrossRef] [PubMed]
14.	A. Nicolet, J.F. Remacle, B. Meys, A. Genon, and W. Legros, “Transformation methods in computational electromagnetics,” J. Appl. Phys. 75 6036 ( 1994). [CrossRef]
15.	R. V. Kohn, H. Shen, M. S. Vogelius, and M. I. Weinstein, “Cloaking via change of variables in electric impedance tomography,” Inverse Probl. 24, 015016 ( 2008). [CrossRef]
16.	W. X. Jiang, J. Y. Chin, Z. Li, Q. Cheng, R. Liu, and T. J. Cui, “Analytical design of conformally invisible cloaks for arbitrarily shaped objects,” Phys. Rev. E 77, 066607 ( 2008) [CrossRef]
17.	A. Nicolet, F. Zolla, and S. Guenneau, “Electromagnetic analysis of cylindrical cloaks of an arbitrary cross section,” Opt. Lett. 33, 1584–1586 ( 2008). [CrossRef] [PubMed]

OCIS Codes
(000.3860) General : Mathematical methods in physics
(160.1190) Materials : Anisotropic optical materials
(260.2110) Physical optics : Electromagnetic optics
(160.3918) Materials : Metamaterials

ToC Category:
Physical Optics

History
Original Manuscript: October 12, 2009
Revised Manuscript: November 12, 2009
Manuscript Accepted: November 14, 2009
Published: November 24, 2009

Citation
Guillaume Dupont, Sébastien Guenneau, Stefan Enoch, Guillaume Demesy, André Nicolet, Frédéric Zolla, and André Diatta, "Revolution analysis of three-dimensional arbitrary cloaks," Opt. Express 17, 22603-22608 (2009)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-17-25-22603

Sort: Author | Year | Journal | Reset

References

J.B. Pendry, D. Shurig and D.R. Smith, "Controlling electromagnetic fields," Science 312, 1780 (2006). [CrossRef] [PubMed]
U. Leonhardt, "Optical conformal mapping," Science 312,1777 (2006). [CrossRef] [PubMed]
A. Alu and N. Engheta, "Achieving Transparency with Plasmonic and Metamaterial Coatings," Phys. Rev. E 95,016623 (2005). [CrossRef]
G. Milton and N. A. Nicorovici, "On the cloaking effects associated with anomalous localized resonance," Proc. R. Soc. London A 462,3027 (2006). [CrossRef]
F. Zolla, S. Guenneau, A. Nicolet, and J.B. Pendry, "Electromagnetic analysis of cylindrical invisibility cloaks and the mirage effect," Opt. Lett. 32, 1069 (2007). [CrossRef] [PubMed]
A. Greenleaf, M. Lassas and G. Uhlmann, "On nonuniqueness for Calderons inverse problem," Math. Res. Lett. 10, 685-693 (2003).
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, 977 (2006). [CrossRef] [PubMed]
P. Zhang, Y. Jin, and S. He, "Obtaining a nonsingular two-dimensional cloak of complex shape from a perfect three-dimensional cloak," Appl. Phys. Lett. 93, 243502 (2008). [CrossRef]
U. Leonhardt and T. Tyc, "Broadband invisibility by non-euclidean cloaking," Science 323, 110 (2009). [CrossRef]
W. Cai, U. K. Chettiar, A. V. Kildiev, and V. M. Shalaev, "Optical Cloaking with metamaterials," Nat. Photon. 1, 224-227 (2007). [CrossRef]
M. Farhat, S. Guenneau, A. B. Movchan, and S. Enoch, "Achieving invisibility over a finite range of frequencies," Opt. Express 16, 5656-5661 (2008). [CrossRef] [PubMed]
C.W. Qiu, L. Hu, X. Xu, and Y. Feng, "Spherical cloaking with homogeneous isotropic multilayered structures," Phys. Rev. E 79, 047602 (2009). [CrossRef]
Y. You, G. W. Kattawar, and P. Yang, "Invisibility cloaks for toroids," Opt. Express 17, 6591 (2009). [CrossRef] [PubMed]
A. Nicolet, J.F. Remacle, B. Meys, A. Genon and W. Legros, "Transformation methods in computational electromagnetics," J. Appl. Phys. 756036 (1994). [CrossRef]
R. V. Kohn, H. Shen,M. S. Vogelius, andM. I.Weinstein, "Cloaking via change of variables in electric impedance tomography," Inverse Probl. 24, 015016 (2008). [CrossRef]
W. X. Jiang, J. Y. Chin, Z. Li, Q. Cheng, R. Liu, and T. J. Cui, "Analytical design of conformally invisible cloaks for arbitrarily shaped objects," Phys. Rev. E 77, 066607 (2008). [CrossRef]
A. Nicolet, F. Zolla, and S. Guenneau, "Electromagnetic analysis of cylindrical cloaks of an arbitrary cross section," Opt. Lett. 33, 1584-1586 (2008). [CrossRef] [PubMed]

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