Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Design of invisibility cloaks with an open tunnel

Open Access Open Access

Abstract

In this paper we apply the methodology of transformation optics for design of a novel invisibility cloak which can possess an open tunnel. Such a cloak facilitates the insertion (retrieval) of matter into (from) the cloak’s interior without significantly affecting the cloak’s performance, overcoming the matter exchange bottleneck inherent to most previously proposed cloak designs. We achieve this by applying a transformation which expands a point at the origin in electromagnetic space to a finite area in physical space in a highly anisotropic manner. The invisibility performance of the proposed cloak is verified by using full-wave finite-element simulations.

© 2010 Optical Society of America

1. Introduction

Transformation optics is an elegant and straightforward formalism used to design optical devices that can control light in an unprecedented manner. In recent years, transformation optics has been applied to the design of various electromagnetic devices such as invisibility cloaks [1, 2], concentrators [3], cylindrical superlenses [4] and field-rotating cloaks [5] etc. Invisibility cloaks are arguably the most interesting optical devices which are designed using this design methodology. To design an invisibility cloak, one first chooses a geometrical transformation which maps a point, line or sheet to a finite region of space (void space), this is followed by an interpretation of the geometrical transformation as a set of material parameters of an effective medium; the resultant cloaks correspond to a point-transformed, line-transformed and surface-transformed (or carpet) cloaks respectively [6]. The material obtained would guide light around the void created by the aforementioned transformation. Any object placed in this void would not scatter incident electromagnetic waves.

There have been several different types of cloaks put forward in the literature. In a 2D scenario, most cloaks have a circular annular cross-section [1, 2, 5], or an elliptical cross-section [7,8]. These cloaking designs have either been verified by geometrical ray tracing [9] or numerical full-wave simulations [10]. For all these cloaking designs, it is not possible to insert (retrieve) objects into (from) the cloaking domain due to their closed annular structures; we call this the matter-exchange bottleneck. In many practical cloaking applications it is desirable for objects to be able to freely enter or leave the cloaking device. Reference [11] reports the design of an open cloak, which allows for material exchange between the cloak’s exterior and interior. The open cloak was designed by positioning the null point before transformation very close to the cloak’s outer boundary, such that there exists only a trivial (close to an identity) coordinate transformation along a certain angular direction. Such a cloak was realized by breaking it down into quadrilaterals; therefore although the approach can easily generate cloaks with different shapes, its design procedure is tedious and the final cloak lacks of closed-form descriptions. In this paper, we propose an open cloak design with a less complicated procedure. Its geometry and in turn the material parameters have analytical expressions. Such a global view of the cloak is especially advantageous for optimization purposes. The open cloak is realized with a transformation function which deforms space in a highly anisotropic manner. We verify the designed cloak, both with and without an opening, using full-wave finite-element simulations.

2. Cardioid cloak

The anisotropic expansion of space is achieved by expanding a point, for our case, at the origin in cylindrical coordinates (r, θ) in electromagnetic (EM) space to a cardioid in cylindrical coordinates (r’, θ’) in physical space. Suppose we have a mapping between the EM and physical spaces in cylindrical coordinates specified in the form r = r(r′,θ′), θ = θ′ and that in EM space we have ɛ = μ = I, where I is the identity matrix and ɛ and μ are the dielectric permittivity and magnetic permeability tensors respectively. Our choice of mapping function between EM and physical space is given by,

r=ra(1+cosθ)ba(1+cosθ)b,
where b denotes the radius of the cloak’s outer boundary and a is a constant determining the size of the cloak’s inner boundary. This mapping function can equivalently be written as
r=ba(1+cosθ)br+a(1+cosθ).
Effectively the transformation maps the origin in EM space r = 0 to a cardioid r′ = a(1 +cosθ) with its cusp located at the origin in physical space. The Jacobian transformation matrix from EM to physical space is given by,
Λ=[rrrrθ00rrθθ0001].
The derivatives in Eq. (3) above are given by,
rr=ba(1+cosθ)b,rθ=a(rb)bsinθ.
The permeability and permittivity tensors in physical space in cylindrical coordinates, following the formula
ɛ=μ=ΛΛT/det(Λ),
are derived as
ɛrr=μrr=(rr)2+(rrθ)2rrrr,
ɛθr=ɛrθ=μθr=μrθ=rrθrr,
ɛθθ=μθθ=rrrr,
ɛzz=μzz=1rrrr.
The above permeability and permittivity tensors in physical space can be re-cast in Cartesian coordinates by using the rotation matrix (R) from cylindrical coordinates to Cartesian coordinates [7], as
ɛcart=μcart=RɛRT/det(R)=RμRT/det(R).

It is interesting to notice that along the angular direction θ′ = π, the mapping function in Eq. (1) reduces to r = r′. This means spaces along this direction undergo an identity transformation; therefore no difference in material parameters would occur there after the transformation. The above statement is also approximately valid for θ′ values close to π. This is indeed observed in the plots of the material parameters for a cardioid cloak with a = 1m and b = 2.5m, shown in Fig. 1. From the plot, we also notice that the zz-component of the material tensors have a finite value, while the remaining components exhibit an infinity at the inner boundary of the cloak. In the following we use COMSOL Multiphysics to verify that our designed cloak can indeed function as desired. For simplicity, we only deal with the 2D in-plane propagation scenario particularly with the transverse-electric (TE) polarization (with Ez, Hx, Hy field components)). Figure 2 shows the scattering pattern of the cardioid cloak (as characterized in Fig. 1) when a plane wave at 0.5 GHz is incident on it, both from the left and from the bottom. The interior surface of the cloak is coated with a perfect electric conductor (PEC). It is noticed that the plane wave passes through the cloak experiencing little perturbation.

 figure: Fig. 1

Fig. 1 Material parameters for the cardioid cloak; a = 1m and b = 2.5m. It should be noticed that the xy- and yx-components of the tensors are equal. In the plot, we have only shown the material parameter values up to 20. Higher values are denoted with a saturated dark-red color.

Download Full Size | PDF

 figure: Fig. 2

Fig. 2 Ez field distribution for a plane wave interacting with a cardioid cloak. (a) left-to-right propagation; (b) bottom-to-top propagation.

Download Full Size | PDF

3. Open cardioid cloak

Next, based on the perfect cardioid cloak, we construct the open cardioid cloak by removing a section of its annular body along the θ′ = π direction. As such a tunnel connecting the interior and the exterior of the cloak can be formed for exchanging objects. We refer to such a cloak as an open cardioid cloak. As a specific example, an open cardioid cloak with a tunnel 2 m and 0.1 m wide at its outer and inner boundaries respectively is discussed in the following. The open cardioid cloak is constructed in COMSOL with its interior surface (except at the opening) coated with PEC. Figure 3 shows the results of the full-wave simulation for this open cloak, again operating at 0.5 GHz. By comparing the field plots in Fig. 3 with those in Fig. 2, one notices that the presence of the tunnel does not significantly affect the cloak’s invisibility performance. Depending on application requirements, tunnels with other opening sizes can be used.

 figure: Fig. 3

Fig. 3 Ez field distribution for a plane wave interacting with an open cardioid cloak. (a) left-to-right propagation; (b) bottom-to-top propagation.

Download Full Size | PDF

To further assess the invisibility performance of the open cardioid cloak, we calculated its far field scattering patterns, as shown in Fig. 4 for plane waves incident from the left (a) and the bottom (b) of the cloak. In the same figure, we also superimposed the scattering patterns for a bare cardioid shaped PEC, and those for the full cardioid cloak studied in section 2. The full cardioid cloak can greatly reduce the scattering of the bare PEC cardioid shaped domain in all directions; the reduction in scattering is on average about 15dB, with a maximum reduction of ∼30dB for the forward scattering direction. Due to the presence of a tunnel, the open cardioid cloak has a slightly larger scattering cross-section than the cardioid cloak. For left-to-right plane wave incidence, the forward scattering is noticeably larger; however compared to the bare PEC, a reduction of ∼20 dB is still achieved. For bottom-to-top plane wave incidence, the open cardioid cloak exhibits slightly larger scattering along top-right and bottom-left directions; but still reductions around ∼10dB are achieved in those directions compared to the bare PEC scatterer case.

 figure: Fig. 4

Fig. 4 Far field scattering patterns of a bare PEC cardioid, a PEC cardioid coated with a cloak, and a PEC cardioid (with an openning) coated with an open cloak, for plane waves incident from (a) left-to- right and (b) bottom-to-top. The maximum far field caused by the bare cardioid PEC is taken as reference, i.e. 0 dB.

Download Full Size | PDF

4. Reduced open cloak

We note that the ideal cloak as well as the open cloak designed above require infinitely large permeability parameters (Fig. 1), which can’t be achieved experimentally. Here we examine the performance of such a cloak at a reduced material parameter space. That is, instead of using unrestricted permeability values, we limit the parameters to a finite range of values. Such a study in part reveals the practicality of the designed cloaks. We refer to cloaks with a reduced material parameter space as reduced cloaks. The emphasis of our study is on the cloak design with an open channel.

Figure 5 shows the far-field scattering patterns of the same open cloak studied in the previous section, but with the maximum permeability set to various μmax values. The simulations are carried out for the left-to-right propagation case. One sees that at μmax = 10 the reduced open cloak still reduces the overall scattering cross-section compared to the bare-PEC case by 5–10dB in most directions. The scattering becomes slightly more inhomogeneous with respect to certain directions, especially around +90° and −90° directions. The scattering cross-section further reduces when one sets μmax = 20. The reduction in forward scattering is as large as 15dB. When μmax = 30, the scattering pattern is approaching that of the open cloak studied in the previous section, i.e. with μmax = ∞, especially along the forward direction. Similar findings are observed with bottom-to-top incidence case (not shown). We therefore conclude that even with a limited material parameter space, one can achieve a reasonable level of reduction in scattering with the designed open cardioid cloak.

 figure: Fig. 5

Fig. 5 Far-field scattering patterns of a bare PEC cardioid, an ideal open cardioid cloak, and reduced open cardioid cloaks designed with various μmax values, for plane waves incident from left-to-right. The maximum far field caused by the bare PEC cardioid is taken as reference i.e. 0 dB.

Download Full Size | PDF

Though challenging at the optical wavelength regime, realization of such a cloaking device is possible at microwave frequencies with current metamaterial fabrication technology. At the microwave regime, the relatively large permeability value can be readily realized with metamaterials based on split-ring resonators [12]. Alternatively for the transverse-magnetic polarization (with Hz, Ex and Ey field components) scenario, the electric-circuit metamaterial proposed in [13] can be deployed to obtain an effective material with a maximum permittivity value around 20.

5. Conclusion

In conclusion, with transformation optics we have designed an invisibility cloak using a highly anisotropic transformation function. Such a cloak can have close-to-background material parameters along a certain angular direction, such that a section of the annular structure can be cut away without affecting its cloaking performance. This effectively opens a tunnel that allows objects to pass into or out of the cloak’s interior freely. Despite the fact that the open cardioid cloak exhibits a slightly deteriorated invisibility performance compared to the cardioid cloak in certain scattering directions, it can still significantly reduce the scattering of, e.g. a bare PEC object.

Acknowledgments

This work is supported by the Swedish Foundation for Strategic Research (SSF), the Swedish Research Council (VR) and the European Commission through an Erasmus Mundus Fellowship.

References and links

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

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

3. 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 the Maxwell’s equations,” Photonics Nanostruct. Fundam. Appl. 6, 87 (2008). [CrossRef]  

4. M. Yan, W. Yan, and M. Qiu, “Cylindrical superlens by a coordinate transformation,” Phys. Rev. B 78, 125113 (2008). [CrossRef]  

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

6. W. Yan, M. Yan, Z. Ruan, and M. Qiu, “Coordinate transformations make perfect invisibility cloaks with arbitrary shape,” N. J. Phys. 10, 043040 (2008). [CrossRef]  

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

8. D. Kwon and D. H. Werner, “Two-dimensional electromagnetic cloak having a uniform thickness for elliptic cylindrical regions,” App. Phys. Lett. 92, 113502 (2008). [CrossRef]  

9. D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express14, 9794 (2006), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-14-21-9794. [CrossRef]   [PubMed]  

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

11. H. Ma, S. Qu, Z. Xu, and J. Wang, “The open cloak,” Appl. Phys. Lett. 94, 103501 (2009). [CrossRef]  

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

13. D. Schurig, J. J. Mock, and D. R. Smith, “Electric field-coupled resonators for negative permittivity metamaterials,” Appl. Phys. Lett. 88, 041109 (2006). [CrossRef]  

Cited By

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

Alert me when this article is cited.


Figures (5)

Fig. 1
Fig. 1 Material parameters for the cardioid cloak; a = 1m and b = 2.5m. It should be noticed that the xy- and yx-components of the tensors are equal. In the plot, we have only shown the material parameter values up to 20. Higher values are denoted with a saturated dark-red color.
Fig. 2
Fig. 2 Ez field distribution for a plane wave interacting with a cardioid cloak. (a) left-to-right propagation; (b) bottom-to-top propagation.
Fig. 3
Fig. 3 Ez field distribution for a plane wave interacting with an open cardioid cloak. (a) left-to-right propagation; (b) bottom-to-top propagation.
Fig. 4
Fig. 4 Far field scattering patterns of a bare PEC cardioid, a PEC cardioid coated with a cloak, and a PEC cardioid (with an openning) coated with an open cloak, for plane waves incident from (a) left-to- right and (b) bottom-to-top. The maximum far field caused by the bare cardioid PEC is taken as reference, i.e. 0 dB.
Fig. 5
Fig. 5 Far-field scattering patterns of a bare PEC cardioid, an ideal open cardioid cloak, and reduced open cardioid cloaks designed with various μmax values, for plane waves incident from left-to-right. The maximum far field caused by the bare PEC cardioid is taken as reference i.e. 0 dB.

Equations (10)

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

r = r a ( 1 + cos θ ) b a ( 1 + cos θ ) b ,
r = b a ( 1 + cos θ ) b r + a ( 1 + cos θ ) .
Λ = [ r r r r θ 0 0 r r θ θ 0 0 0 1 ] .
r r = b a ( 1 + cos θ ) b , r θ = a ( r b ) b sin θ .
ɛ = μ = Λ Λ T / det ( Λ ) ,
ɛ r r = μ r r = ( r r ) 2 + ( r r θ ) 2 r r r r ,
ɛ θ r = ɛ r θ = μ θ r = μ r θ = r r θ r r ,
ɛ θ θ = μ θ θ = r r r r ,
ɛ z z = μ z z = 1 r r r r .
ɛ cart = μ cart = R ɛ R T / det ( R ) = R μ R T / det ( R ) .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.