Design of invisibility cloaks with an open tunnel

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 OCIS codes: (160.1190) Anisotropic optical materials; (160.3918) Metamaterials; (230.3205) Invisibility cloaks. References and links 1. J.B. Pendry, D. Schurig and D.R. Smith, “Controlling electromagnetic fields,” Science 312, 1780 (2006). 2. U. Leonhardt, “Optical conformal mapping,” Science 312, 1777 (2006). 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,” Photon. Nanostruct. : Fundam. Appl. 6, 87 (2008). 4. M. Yan, W. Yan and M. Qiu, “Cylindrical superlens by a coordinate transformation,” Phys. Rev. B 78, 125113 (2008). 5. H. Chen and C.T. Chan, “Transformation media that rotate electromagnetic fields,” Appl. Phys. Lett. 90, 241105 (2007). 6. W. Yan, M. Yan, Z. Ruan and M. Qiu, “Coordinate transformations make perfect invisibility cloaks with arbitrary shape,” New J. Phys. 10, 043040 (2008). 7. H. Chen, “Transformation optics in orthogonal coordinates,” J. Opt. A: Pure Appl. Opt. 11, 075102 (2009). 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). 9. D. Schurig, J.B. Pendry, and D.R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express 14, 9794 (2006). http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-14-21-9794 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). 11. H. Ma, S. Qu, Z. Xu and J. Wang, “The open cloak,” App. Phys. Lett. 94, 103501 (2009). 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). 13. D. Schurig, J.J. Mock and D.R. Smith, “Electric field-coupled resonators for negative permittivity metamaterials,” Appl. Phys. Lett. 88, 041109 (2006).


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 surfacetransformed (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.Ref [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.

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, 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 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, The derivatives in Eq. 3 above are given by, The permeability and permittivity tensors in physical space in cylindrical coordinates, following the formula are derived as 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 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 E z , H x , H y 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.

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

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.
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 H z , E x and E y 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.

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.

Fig. 1 .
Fig. 1.Material parameters for the cardioid cloak; a = 1m and b = 2.5m.It should be noticed that the xyand 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. E z field distribution for a plane wave interacting with a cardioid cloak.(a) left-to-right propagation; (b) bottom-to-top propagation.

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

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.