OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 16, Iss. 9 — Apr. 28, 2008
  • pp: 6134–6145
« Show journal navigation

Invisibility cloaks for irregular particles using coordinate transformations

Yu You, George W. Kattawar, Peng-Wang Zhai, and Ping Yang  »View Author Affiliations


Optics Express, Vol. 16, Issue 9, pp. 6134-6145 (2008)
http://dx.doi.org/10.1364/OE.16.006134


View Full Text Article

Acrobat PDF (1877 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

Invisibility cloaks for ellipsoids, rounded cuboids and rounded cylinders have been studied on the basis of the coordinate transformation approach. The resultant material property tensors for irregular cloaks are more complicated in comparison with those for the spherical invisibility cloak. A generalized Discrete Dipole Approximation (DDA) formalism has been used to simulate the scattered field distribution in the vicinity of the aforementioned irregular cloaks illuminated by an incident plane wave. Simulated scattering efficiencies are on the order of 10-5, and the simulated electric-field distribution outside of a cloak is the same as that of the incident radiation.

© 2008 Optical Society of America

1. Introduction

Controlling electromagnetic fields using material properties given by the coordinate transformation approach [1

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

, 2

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

, 3

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

, 4

4. D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express 14, 9794–9804 (2006). [CrossRef] [PubMed]

, 5

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

, 6

6. U. Leonhardt, “Notes on conformal invisibility devices,” New J. Phys. 8, 118 (2006). [CrossRef]

, 7

7. D. Schurig, J. B. Pendry, and D. R. Smith, “Transformation-designed optical elements,” Opt. Express 15, 14772–14782 (2007). [CrossRef] [PubMed]

] has recently drawn extensive attention in the research community. The coordinate transformation approach has been discussed in a general coordinate system [1

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

, 2

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

, 3

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

] and in Cartesian coordinates [4

4. D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express 14, 9794–9804 (2006). [CrossRef] [PubMed]

]. A coordinate transformation method in the form of optical conformal mapping [5

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

, 6

6. U. Leonhardt, “Notes on conformal invisibility devices,” New J. Phys. 8, 118 (2006). [CrossRef]

] in arbitrary coordinates has also been reported. It has been demonstrated, for example, in Ref. [3

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

], that the transformation properties of Maxwell’s equations under certain coordinate transformations can yield material properties that have interesting functionalities, such as invisibility cloak [2

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

, 3

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

, 4

4. D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express 14, 9794–9804 (2006). [CrossRef] [PubMed]

], perfect lens [3

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

], and magnification [3

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

, 7

7. D. Schurig, J. B. Pendry, and D. R. Smith, “Transformation-designed optical elements,” Opt. Express 15, 14772–14782 (2007). [CrossRef] [PubMed]

], to name a few.

Due to its impressive applications, the invisibility cloak has received a more intensive study. The permittivity and permeability tensors have been given for anisotropic and inhomogeneous spherical [2

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

, 4

4. D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express 14, 9794–9804 (2006). [CrossRef] [PubMed]

] and infinite cylindrical [3

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

, 8

8. S. A. Cummer, B.-I. Popa, D. Schurig, D. R. Smith, and J. B. Pendry, “Full-wave simulations of electromagnetic cloaking structures,” Phys. Rev. E 74, 036621 (2006). [CrossRef]

, 9

9. 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]

] shells that give absolute zero scattering cross section subject to external incident radiation when an arbitrary object is embedded inside. Both analytical and numerical calculations have been performed to verify the cloaking effects. Ray-tracing simulations [2

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

, 4

4. D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express 14, 9794–9804 (2006). [CrossRef] [PubMed]

] in the geometric optics limit have been reported for a cloaked sphere with consistent results. Rigorous solutions to Maxwell’s equations in the spherical case [10

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

] and in the 2-D cylindrical case [11

11. B. Zhang, H. Chen, B.-I. Wu, Y. Luo, L. Ran, and J. A. Kong, “Response of a cylindrical invisibility cloak to electromagnetic waves,” Phys. Rev. B 76, 121101 (2007). [CrossRef]

] have been reported, confirming the cloaking effects on objects with these geometries. Numerical simulations using the full-wave finite-element method were performed to study the effects of cloaking material perturbations on the propagation of the incident waves associated with 2-D invisibility cloaks of various shapes, such as cylinders [8

8. S. A. Cummer, B.-I. Popa, D. Schurig, D. R. Smith, and J. B. Pendry, “Full-wave simulations of electromagnetic cloaking structures,” Phys. Rev. E 74, 036621 (2006). [CrossRef]

], squares [12

12. 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 transformations of Maxwell’s equations,” Photon. Nanostruct.: Fundam. Applic. 6, 87–95 (2008). [CrossRef]

], elliptical cylinders [13

13. H. Ma, S. Qu, Z. Xu, J. Zhang, B. Chen, and J. Wang, “Material parameter equation for elliptical cylindrical cloaks,” Phys. Rev. A 77, 013825 (2008). [CrossRef]

], and eccentric elliptical cylinders [14

14. D.-H. Kwon and D. H. Werner, “Two-dimensional eccentric elliptic electromagnetic cloaks,” Appl. Phys. Lett. 92, 013505 (2008). [CrossRef]

]. Numerical simulations using the Discrete Dipole Approximation method [15

15. Y. You, G. W. Kattawar, P.-W. Zhai, and P. Yang, “Zero-backscatter cloak for aspherical particles using a generalized DDA formalism,” Opt. Express 16, 2068–2079 (2008). [CrossRef] [PubMed]

] were also reported in a 3-D spherical case. In all of these numerical studies, simulated local field distribution in the vicinity of the cloaking object has been found to be the same as that of the incident radiation. With the help of metamaterial technology, the required material specification can be practically implemented. For example, 2-D cylindrical cloaking has been realized using artificially structured metamaterials at microwave frequencies [9

9. 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]

] and at optical frequencies [16

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

].

The interaction between a spherical invisibility cloak and incident radiation from internal sources has also been investigated [17

17. B. Zhang, H. Chen, B.-I. Wu, and J. A. Kong, “Extraordinary surface voltage effect in the invisibility cloak with an active device inside,” Phys. Rev. Lett. 100, 063904 (2008). [CrossRef] [PubMed]

]. The electromagnetic field solution for a spherical cloak with an active device inside shows that extraordinary surface voltages induced at the inner boundary of the cloak prevent electromagnetic waves from going out.

The Discrete Dipole Approximation (DDA) [18

18. E. M. Purcell and C. R. Pennypacker, “Scattering and absorption of light by nonspherical dielectric grains,” Astrophys. J. 186, 705–714 (1973). [CrossRef]

, 19

19. B. T. Draine, “The discrete-dipole approximation and its application to interstellar graphite grains,” Astrophys. J. 333, 848–872 (1988). [CrossRef]

, 20

20. B. T. Draine and P. J. Flatau, “Discrete-dipole approximation for scattering calculations,” J. Opt. Soc. Am. A 11, 1491–1499 (1994). [CrossRef]

] method is a robust numerical technique for the computation of elastic scattering and absorption by dielectric particles with arbitrary shapes. The conventional DDA formalism has been generalized to computations pertaining to particles with permeabilities μ≠1 [15

15. Y. You, G. W. Kattawar, P.-W. Zhai, and P. Yang, “Zero-backscatter cloak for aspherical particles using a generalized DDA formalism,” Opt. Express 16, 2068–2079 (2008). [CrossRef] [PubMed]

] by representing the particle by a cubic array of both electric and magnetic oscillating dipoles. This formalism has been proven to be quite useful for simulating the field distributions and scattering cross sections of cloaked spheres.

To date, studies on invisibility cloaks have been limited to spheres and several 2-D irregular particles. To the authors’ best knowledge, the application of invisibility cloak to 3-D particles with irregular shapes has not been reported in the literature.

In this study, we first apply the coordinate transformation approach to ellipsoids, 3-D rounded cuboids, and 3-D rounded cylinders, and give the material properties required for the cloaks with these geometries. Furthermore, we simulate the electric field distributions and the scattering cross sections associated with cloaked objects using the DDA method.

2. Transformation equations

In Cartesian coordinates, a point in space is described by its Cartesian components xi, with i=1,2,3. Under a coordinate transformation that maps point xi to point xi=xi(xi) , a rank 2 contravariant tensor Tij is transformed as follows:

Tij=det(Λii)1ΛiiΛjjTij,
(1)

where

Λii=xixi
(2)

is the Jacobian transformation matrix. As rank 2 contravariant tensors that describe material properties, the permittivity tensor εij and the permeability tensor μij are transformed as

εij=det(Λii)1ΛiiΛjjεij,
(3a)
μij=det(Λii)1ΛiiΛjjμij.
(3b)

Pendry et al. [2

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

] and Schurig et al. [4

4. D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express 14, 9794–9804 (2006). [CrossRef] [PubMed]

] suggested that the left hand sides of Eq.(3a) and Eq.(3b) be interpreted as either the properties of the same material in the transformed coordinate system (the topological interpretation), or the permittivity and permeability of an actual material or metamaterial in a flat Cartesian coordinate system (the material interpretation).

As has been discussed in Ref. [4

4. D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express 14, 9794–9804 (2006). [CrossRef] [PubMed]

], the material properties of an invisibility cloak can be determined as illustrated in Fig. 1. Consider a closed domain that is to be transformed (region (I) in Fig. 1(a)). The internal and external regions of this domain are assumed to be vacuum. Next, consider a coordinate transformation that maps this region into a region (region (I′) in Fig. 1(b)) that has the same outer boundary, but contains a hole (region (III)) bounded by an inner boundary. Then, determine the material properties εij and μij in region (I′), which no longer represent a vacuum, on the basis of Eqs.(3a) and (3b). Meanwhile, the external region (region (II)) remains undistorted, and the material properties εij and μij in this region still represent a vacuum. In this manner, we have defined a cloaking region. Subject to radiation from outside sources, this object compresses the radiation field in region (I) into region (I′), the cloaking region, and leaves region (III), the cloaked region, radiation-free. This implies that there is no interaction between anything in the cloaked region and sources in the outside domain. Outside of the cloaking region, the radiation fields remain unchanged, as if neither the cloaking material nor any cloaked object exists. In principle, the regions (I), (I′), and (III) can be of arbitrary shapes, although the transformation matrix Eq.(2) associated with irregular shapes may be quite complicated.

Fig. 1. Illustration of the coordinate transformation used to determine the material properties of a cloak.

3. Coordinate transformations for ellipsoidal cloaks

The coordinate transformation used in the studies of spherical cloaks [2

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

, 4

4. D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express 14, 9794–9804 (2006). [CrossRef] [PubMed]

] can be readily generalized to ellipsoids. For an ellipsoid with its semi-axes lying along x-, y-, and z-directions with lengths α 1 b,α 2 b, and α 3 b, respectively, the outer boundary is described by

(x1α1b)2+(x2α2b)2+(x3α3b)2=1,
(4)

where α 1, α 2, and α 3 are positive numbers describing the aspect ratio of the ellipsoid, and do not transform as a vector. Inside this boundary, we start with the coordinate transformation determined by Eq.(14) in Ref. [4

4. D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express 14, 9794–9804 (2006). [CrossRef] [PubMed]

],

xi=babxiδii+axirδii,(a<b)
(5)

where δii is the Kronecker tensor, r=(x1)2+(x2)2+(x3)2 is the distance between the point and the origin. Furthermore, we re-define r as a scaled distance

r=[i=13(xiαi)2]12.
(6)

This definition of r guarantees that r=b for any point xi on the outer boundary. Therefore, xi=xi . This means that the outer boundary is not changed by the transformation.

By writing the coordinate transformation in the fashion of Eq.(5), we have a singularity at the origin. Therefore, to investigate how the origin transforms, we need to avoid the singularity and study an infinitesimal ellipsoid bounded by the surface described by

(x1α1b)2+(x2α2b)2+(x3α3b)2=ε2,
(7)

where ε is a small positive number. In the limit when ε approaches 0, this surface becomes the origin. A little algebra will show that the transformed surface becomes

(x1α1a)2+(x2α2a)2+(x3α3a)2=1,
(8)

which is the boundary of a smaller concentric ellipsoid. Therefore, the coordinate transformation given by Eq.(5) and Eq.(6) transforms a closed ellipsoid with lengths of the outer semi-axes α 1 b,α 2 b, and α 3 b into an ellipsoidal shell with the same lengths of the outer semi-axes and lengths of the inner semi-axes α 1 a,α 2 a, and α 3 a. Outside this region, we assume the identity transformation. In the following discussions, all equations apply only to the internal region.

In Fig. 2 we show the application of the transformation Eq.(5) to a closed ellipsoidal region. Here we characterize the size of an irregular region in terms of the size parameter of an equal-volume sphere x=2πa eff/λ where a eff=(3V/4π)1/3 with V the volume of the irregular region, and λ is the wavelength of some incident radiation. The ellipsoid shown in Fig. 2 has a size parameter of x=8, with α 1=α 2=1, and α 3=2, or b=21/38λ/(2π). This ellipsoid is practically a spheroid with its axis of symmetry lying along the z direction.

Fig. 2. Application of the coordinate transformation to a spheroidal region subject to a plane-wave radiation propagating in the z-direction. Shown in the figures is the situation in the x-z plane, with the x- and z-axes as depicted in Fig. (a). The scales are in unit of λ. Figure (a) shows the closed region (green region) in Cartesian coordinates, with vacuum both inside and outside. The light rays (dark red lines) and wave fronts (dark blue lines) of the incident radiation are also shown. Figure (b) shows the transformed region according to the transformation Eq.(5) with a=0.5b, with the material properties in the cloaking region determined by Eq.(3a) and Eq.(3b). The light rays and wave fronts transform accordingly.

Figure 2(a) shows the ellipsoidal region, along with the Cartesian coordinate grid, before the coordinate transformation. As we assume a vacuum everywhere, the coordinate grid lines can also be interpreted as light rays and wave fronts (loci of points having the same phase) if the incident radiation propagates in the z-direction. A coordinate transformation Eq.(5) with a=0.5b transforms this closed region into the cloaking region shown in Fig. 2(b). The grid lines in the region transform accordingly, and the transformed grid lines, which are geodesics in the transformed coordinates, can be interpreted as light rays and wave fronts in the presence of the cloak. As expected, the radiation field never penetrates into the cloaked region. Evidently, the light rays and wave fronts beyond the outer boundary remain the same, as an identity transformation is applied in this region.

Note that the coordinate transformation does not depend on the orientation of the incident radiation beam with respect to the cloak. Therefore, it also works for an incident radiance making an arbitrary angle θ to the z-direction. We just need to do the coordinate transformation in the same Cartesian coordinates spanned by the three semi-axes of the ellipsoid as can be seen in Fig. 3, where we show the transformation applied to the same ellipsoid, but subject to an incident radiance in the x-z plane making an angle of θ=30° to the z-axis.

Here, a rotated coordinate grid represents the light rays and the wave fronts. Again, we notice from Fig. 3(b) that the light rays and wave fronts inside the boundary are transformed such that the cloaked region is never reached.

Fig. 3. Same as Fig. 2, but for an incident radiation in the x-z plane making an angle of θ=30° to the z-axis.

The material properties εij and μij in the cloaking region can be determined by substituting Eqs.(5) and (6) into Eq.(2), leading to a transformation matrix as follows:

Λji=rrδjiaxixkδiiδkjαj2r3
(9)

in Cartesian coordinates, where the Einstein summation convention has been assumed for indices i and k. For the vacuum, εij=μij=δij. Therefore, a combination of Eq.(1) and Eq.(9) gives the material properties in the cloaking region as follows:

εij=μij=bba[δijxixj(ar3(αi2+αj2)a2r~2r4k=13(x~k)2αk4)].
(10)

As we take the material interpretation, we have dropped all primes and relabeled the variables as follows: xii, r, xixi , and r′→r, with r given by Eq.(6), and i given by

x~i=xibab+ar=andr~=raba·b.
(11)

Therefore, we have obtained the permittivity and permeability tensors of an ellipsoidal cloak. In the special case of spheres, αi=1, Eq.(10) reduces to the simple results described by the Eq.(20) in Ref. [4

4. D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express 14, 9794–9804 (2006). [CrossRef] [PubMed]

].

4. Coordinate transformations for rounded-cuboids

Rahm et al. [12

12. 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 transformations of Maxwell’s equations,” Photon. Nanostruct.: Fundam. Applic. 6, 87–95 (2008). [CrossRef]

] have studied the application of the coordinate transformation approach to 2-D squares. The transformation they used can be easily generalized to 3-D cases. This transformation, however, leads to discontinuities in the variance of material properties of the cloak. We will present material properties of cloaks for rounded-cuboids with arbitrary aspect ratio, determined by an alternative coordinate transformation. In our approach, discontinuities are not involved.

We seek the possibility of representing a cuboid approximately by a superellipsoid [22

22. I. D. Faux and M. J. Pratt, Computational geometry for design and manufacture, Wiley, Chichester (1979).

, 23

23. T. Wriedt, “Using the T-Matrix method for light scattering computations by non-axisymmetric particles: superellipsoids and realistically shapes particles,” Part. Part. Syst. Charact. 19, 256–268 (2002). [CrossRef]

] bounded by the surface

[(xα1b)2e+(yα2b)2e]em+(zα3b)2m=1,
(12)

which is well known in computer graphics and used to model a wide range of shapes, including rounded-cubes and rounded-cylinders [22

22. I. D. Faux and M. J. Pratt, Computational geometry for design and manufacture, Wiley, Chichester (1979).

]. Wriedt [23

23. T. Wriedt, “Using the T-Matrix method for light scattering computations by non-axisymmetric particles: superellipsoids and realistically shapes particles,” Part. Part. Syst. Charact. 19, 256–268 (2002). [CrossRef]

] introduced this shape to the light scattering community in his study of the T-Matrix method.

To approximate a cuboid, we let e=m=2/n. Therefore, Eq.(12) becomes

(x1α1b)n+(x2α2b)n+(x3α3b)n=1.
(13)

Eq.(13) is a direct generalization of Eq.(4) by replacing the exponent 2 with an integer n. We call this shape an order-n-cuboid for convenience. As can be seen in Fig. 4(a), an order-n-cuboid approaches a cube as the exponent n increases. In Fig. 4(b) we show the scattering patterns associated with an order-10-cuboidal scattering particle with α 1=α 2=α 3 and with a cubical scattering particle. As can be seen in this figure, an order-n-cuboid with n=10 is already a good approximation of a cuboid as far as light scattering is concerned.

Fig. 4. (a) Approximations of a cube by order-n-cuboids with various n values, and α 1=α 1=α 3. (b) Comparison of phase functions pertaining to an order-10-cuboidal scattering particle with α 1=α 2=1, α 3=2 and pertaining to a cubical scattering particle with the same aspect ratio. Both particles are homogeneous and have the same size parameter x=4 and the same permittivity ε=1.44.

To determine the material property tensors of an order-n-cuboidal cloak, we can still use the coordinate transformation Eq.(5), and redefine the scaled magnitude r again as

r=[i=13(xiαi)n]1n.
(14)

It can be easily verified that the boundary transforms to itself and the origin transforms to an inner boundary with b replaced by a.

Substituting Eq.(5) and Eq.(14) into Eq.(2), we obtain a transformation matrix as follows:

Λji=rrδjiaxi(xk)n1δiiδkjαjnrn+1,
(15)

and a combination of Eq.(1) and Eq.(15) gives the permittivity and permeability tensors for an order-n-cuboidal cloak

εij=μij
=bba[δijxixj×
(ar3r~n2((x~i)n2αin+(x~j)n2αjn)a2r~2(n1)r4k=13(x~k)2n2αk2n)],
(16)

with r given by Eq.(14), and i given by Eq.(11).

In Fig. 5, we show a rounded cuboidal cloak, approximated by an order-10-cuboid with α 1=α 2=1, α 3=2, and size parameter x=8. The light rays and wave fronts in the vicinity of the cloak as predicted by the coordinate transformation are shown. Two particle orientations were considered. Again, the light rays and wave fronts deviate from the cloaked region. Therefore, any object in the cloaked region is hidden from outside observers.

Fig. 5. Same as Fig. 2(b), but for a rounded cuboidal cloak approximated by an order- 10-cuboid with x=8, a=0.5b, α 1=α 2=1, and α 3=2, subject to incident radiation propagating in various directions. (a) The incident radiation is along the z-axis; (b) The incident radiation is in the x-z plane, making an angle of θ=30° to the z-axis.

5. Coordinate transformations for rounded cylinders

To approximate an elliptic cylinder of finite height, we let e=1 and m=2/n. Therefore, Eq.(12) becomes

[(x1α1b)2+(x2α2b)2]n2+(x3α3b)n=1,
(17)

which defines a shape we will call an order-n-cylinder. As Fig. 6(a) shows, an order-n-cylinder approaches a cylinder as the exponent n increases. Figure 6(b) implies that an order-n-cylinder with n=10 gives a good approximation to a cylinder as far as light scattering is concerned.

We can then use the same coordinate transformation Eq.(5) to determine material property tensors of an order-n-cylindrical cloak by redefining the scaled magnitude r as

r=([(x1α1)2+(x2α2)2]n2+(x3α3)n)1n.
(18)

Again, the boundary transforms to itself and the origin transforms to an inner boundary with

Fig. 6. (a) Approximations of a cylinder by order-n-cylinders with various n values, and α 1=α 1=α 3. (b) Same as Fig. 4(b), but for an order-10-cylindrical scattering particle of the same size parameter and aspect ratio, and a cylindrical scattering particle.

b replaced by a. The corresponding transformation matrix will then be

Λji=rrδjiarn+1xiδii[xkδkjρn2αj2(δ1j+δ2j)+(xk)n1δkjαj2δ3j],
(19)

and the permittivity/permeability tensor for an order-n-cylindrical cloak is a hybrid of the corresponding tensors for an ellipsoidal cloak and for an order-n-cuboidal cloak, given by

εij=μij
=bba[δijxixj×
(ar3r~n2[ρn2(1αi2(δ1i+δ2i)+1αj2(δ1j+δ2j))+
(x~i)n2αinδ3i+(x~j)n2αjnδ3j]
a2r~2(n1)r4[ρ2n4((x~1)2α14+(x~2)2α24)+(x~3)2n2α32n])],
(20)

with r given by Eq.(18), and i given by Eq.(11).

Note that when z=0, Eq.(17) and Eq.(4) are identical, and Eq.(18) and Eq.(6) are identical; When x=0 or y=0, Eq.(17) and Eq.(13) are identical, and Eq.(18) and Eq.(14) are identical. Namely, for an order-n-cylindrical cloak, the wave fronts and light rays behave the same as those associated with an ellipsoidal cloak in the x-y plane (similar to the situation shown in Fig. 2(b) and Fig. 3(b)), and the same as those associated with an order-n-cuboidal cloak in the x-z plane and the y-z plane (same as the situation shown in Fig. 5).

To conclude this section, in Fig. 7 we present 3-D views of light rays and wave fronts associated with the three irregular invisibility cloaks we discussed.

Fig. 7. 3-D views of light rays and wave fronts associated with the three irregular invisibility cloaks: (a) ellipsoid, (b) rounded-cuboid, and (c) rounded cylinder.

6. DDA simulations

To check if these coordinate transformations really work, we simulated the light scattering pertaining to spheroidal cloaks, rounded cuboidal cloaks, and rounded cylindrical cloaks with the permittivity and permeability tensors given by Eq.(10), Eq.(16), and Eq.(20), respectively, using the DDA formalism discussed in Ref. [15

15. Y. You, G. W. Kattawar, P.-W. Zhai, and P. Yang, “Zero-backscatter cloak for aspherical particles using a generalized DDA formalism,” Opt. Express 16, 2068–2079 (2008). [CrossRef] [PubMed]

]. The situations shown in Fig. 2(b), Fig. 3(b), and Fig. 5 are considered. Figure 8 shows the simulated electric-field distributions in the vicinity of an ellipsoidal cloak. DDA simulations for electric field distributions in the vicinity of a rounded cuboidal cloak are presented in Fig. 9. Both figures show the situation in the x-z plane. The DDA simulations in the x-z plane for a rounded cylindrical cloak look similar to Fig. 9. It can be noticed that the plane-wave feature outside of the cloak is perfectly kept. The field in the cloaking region is compressed with patterns that are consistent with the predictions of the coordinate transformation approach as can be seen in Fig. 2(b), Fig. 3(b), and Fig. 5. The field in the cloaked region is close to 0 with a leakage of about 10% of the radiating field, or about 1% of the radiative energy, into the cloaked region due to the discretization of material properties in the DDA calculations.

As a numerical method, the DDA is not expected to give exact zero scattering. In our calculations for the three cloaks with size parameter x=8, cloak parameter a=0.5b, and two particle orientations, all simulated scattering efficiencies Q sca=C sca/(πa 2 eff) are on the order of 10-3, which is 3 orders lower than that for a regular dielectric particle. The simulated scattering efficiencies are on the order of 10-5 for smaller cloaks (x=5) applied to smaller cloaked regions (a=0.3b).

Fig. 8. The DDA simulations of the electric-field distribution in the vicinity of an ellipsoidal cloak in the x-z plane.
Fig. 9. The DDA simulations of the electric-field distribution in the vicinity of an order-10- cuboidal cloak in the x-z plane.

7. Conclusions

Invisibility cloaks for particles with three irregular geometries have been studied on the basis of the coordinate transformation approach. The permittivity and permeability tensors have been given for ellipsoidal cloaks, 3-D rounded-cuboidal cloaks approximated by order-n-cuboidal cloaks, and 3-D rounded-cylindrical cloaks approximated by order-n-cylindrical cloaks. Numerical calculations using a generalized Discrete Dipole Approximation method were carried out to simulate the light scattering of plane-wave incident radiation associated with such cloaking objects. The simulated electric-field distributions in the vicinity of the cloak and scattering efficiencies suggest that these cloaks do not change the electric field outside at all. Therefore, outside observers cannot detect any object embedded within the cloaked region. Moreover, simulated scattering efficiencies of as low as 10-5 have been observed.

Unlike spherical cloaks, for which the permittivity and permeability tensors are rotationally uniaxial, cloaks for irregular particles have more complicated material properties. This may cause extra difficulties in the fabrication of these cloaks using metamaterial or other materials.

Acknowledgments

References and links

1.

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

2.

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

3.

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

4.

D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express 14, 9794–9804 (2006). [CrossRef] [PubMed]

5.

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

6.

U. Leonhardt, “Notes on conformal invisibility devices,” New J. Phys. 8, 118 (2006). [CrossRef]

7.

D. Schurig, J. B. Pendry, and D. R. Smith, “Transformation-designed optical elements,” Opt. Express 15, 14772–14782 (2007). [CrossRef] [PubMed]

8.

S. A. Cummer, B.-I. Popa, D. Schurig, D. R. Smith, and J. B. Pendry, “Full-wave simulations of electromagnetic cloaking structures,” Phys. Rev. E 74, 036621 (2006). [CrossRef]

9.

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]

10.

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

11.

B. Zhang, H. Chen, B.-I. Wu, Y. Luo, L. Ran, and J. A. Kong, “Response of a cylindrical invisibility cloak to electromagnetic waves,” Phys. Rev. B 76, 121101 (2007). [CrossRef]

12.

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 transformations of Maxwell’s equations,” Photon. Nanostruct.: Fundam. Applic. 6, 87–95 (2008). [CrossRef]

13.

H. Ma, S. Qu, Z. Xu, J. Zhang, B. Chen, and J. Wang, “Material parameter equation for elliptical cylindrical cloaks,” Phys. Rev. A 77, 013825 (2008). [CrossRef]

14.

D.-H. Kwon and D. H. Werner, “Two-dimensional eccentric elliptic electromagnetic cloaks,” Appl. Phys. Lett. 92, 013505 (2008). [CrossRef]

15.

Y. You, G. W. Kattawar, P.-W. Zhai, and P. Yang, “Zero-backscatter cloak for aspherical particles using a generalized DDA formalism,” Opt. Express 16, 2068–2079 (2008). [CrossRef] [PubMed]

16.

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

17.

B. Zhang, H. Chen, B.-I. Wu, and J. A. Kong, “Extraordinary surface voltage effect in the invisibility cloak with an active device inside,” Phys. Rev. Lett. 100, 063904 (2008). [CrossRef] [PubMed]

18.

E. M. Purcell and C. R. Pennypacker, “Scattering and absorption of light by nonspherical dielectric grains,” Astrophys. J. 186, 705–714 (1973). [CrossRef]

19.

B. T. Draine, “The discrete-dipole approximation and its application to interstellar graphite grains,” Astrophys. J. 333, 848–872 (1988). [CrossRef]

20.

B. T. Draine and P. J. Flatau, “Discrete-dipole approximation for scattering calculations,” J. Opt. Soc. Am. A 11, 1491–1499 (1994). [CrossRef]

21.

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

22.

I. D. Faux and M. J. Pratt, Computational geometry for design and manufacture, Wiley, Chichester (1979).

23.

T. Wriedt, “Using the T-Matrix method for light scattering computations by non-axisymmetric particles: superellipsoids and realistically shapes particles,” Part. Part. Syst. Charact. 19, 256–268 (2002). [CrossRef]

OCIS Codes
(160.1190) Materials : Anisotropic optical materials
(290.0290) Scattering : Scattering
(160.3918) Materials : Metamaterials
(280.1350) Remote sensing and sensors : Backscattering
(260.2710) Physical optics : Inhomogeneous optical media

ToC Category:
Metamaterials

History
Original Manuscript: February 20, 2008
Revised Manuscript: April 4, 2008
Manuscript Accepted: April 4, 2008
Published: April 16, 2008

Citation
Yu You, George W. Kattawar, Peng-Wang Zhai, and Ping Yang, "Invisibility cloaks for irregular particles using coordinate transformations," Opt. Express 16, 6134-6145 (2008)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-9-6134


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. A. J. Ward and J. B. Pendry, "Refraction and geometry in Maxwell??s equations," J. Mod. Opt. 43, 773-793 (1996). [CrossRef]
  2. J. B. Pendry, D. Schurig, and D. R. Smith, "Controlling electromagnetic fields," Science 312, 1780 (2006). [CrossRef] [PubMed]
  3. U. Leonhardt and T. G. Philbin, "General relativity in electrical engineering," New J. Phys. 8, 247 (2006). [CrossRef]
  4. D. Schurig, J. B. Pendry, and D. R. Smith, "Calculation of material properties and ray tracing in transformation media," Opt. Express 14, 9794-9804 (2006). [CrossRef] [PubMed]
  5. U. Leonhardt, "Optical conformal mapping," Science 312, 1777-1780 (2006). [CrossRef] [PubMed]
  6. U. Leonhardt, "Notes on conformal invisibility devices," New J. Phys. 8, 118 (2006). [CrossRef]
  7. D. Schurig, J. B. Pendry, and D. R. Smith, "Transformation-designed optical elements," Opt. Express 15, 14772-14782 (2007). [CrossRef] [PubMed]
  8. S. A. Cummer, B.-I. Popa, D. Schurig, D. R. Smith, and J. B. Pendry, "Full-wave simulations of electromagnetic cloaking structures," Phys. Rev. E 74, 036621 (2006). [CrossRef]
  9. 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]
  10. H. Chen, B.-I. Wu, B. Zhang, and J. A. Kong, "Electromagnetic wave interaction with a metamaterial cloak," Phys. Rev. Lett. 99, 063903 (2007). [CrossRef] [PubMed]
  11. B. Zhang, H. Chen, B.-I. Wu, Y. Luo, L. Ran, and J. A. Kong, "Response of a cylindrical invisibility cloak to electromagnetic waves," Phys. Rev. B 76, 121101 (2007). [CrossRef]
  12. 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 transformations of Maxwell??s equations," Photon. Nanostruct. Fundam. Appl. 6, 87-95 (2008). [CrossRef]
  13. H. Ma, S. Qu, Z. Xu, J. Zhang, B. Chen, and J. Wang, "Material parameter equation for elliptical cylindrical cloaks," Phys. Rev. A 77, 013825 (2008). [CrossRef]
  14. D.-H. Kwon and D. H. Werner, "Two-dimensional eccentric elliptic electromagnetic cloaks," Appl. Phys. Lett. 92, 013505 (2008). [CrossRef]
  15. Y. You, G. W. Kattawar, P.-W. Zhai, and P. Yang, "Zero-backscatter cloak for aspherical particles using a generalized DDA formalism," Opt. Express 16, 2068-2079 (2008). [CrossRef] [PubMed]
  16. W. Cai, U. K. Chettiar, A. V. Kildishev, and V. M. Shavlaev, "Optical cloaking with metamaterials," Nat. Photon. 1, 224-227 (2007). [CrossRef]
  17. B. Zhang, H. Chen, B.-I. Wu, and J. A. Kong, "Extraordinary surface voltage effect in the invisibility cloak with an active device inside," Phys. Rev. Lett. 100, 063904 (2008). [CrossRef] [PubMed]
  18. E. M. Purcell and C. R. Pennypacker, "Scattering and absorption of light by nonspherical dielectric grains," Astrophys. J. 186, 705-714 (1973). [CrossRef]
  19. B. T. Draine, "The discrete-dipole approximation and its application to interstellar graphite grains," Astrophys. J. 333, 848-872 (1988). [CrossRef]
  20. B. T. Draine and P. J. Flatau, "Discrete-dipole approximation for scattering calculations," J. Opt. Soc. Am. A 11, 1491-1499 (1994). [CrossRef]
  21. J. D. Jackson, Classical Electrodynamics (John Wiley & Sons, New York, 1975).
  22. I. D. Faux and M. J. Pratt, Computational geometry for design and manufacture (Wiley, Chichester 1979).
  23. T. Wriedt, "Using the T-Matrix method for light scattering computations by non-axisymmetric particles: superellipsoids and realistically shapes particles," Part. Part. Syst. Charact. 19, 256-268 (2002). [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.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited