## Abstract

Electromagnetic external cloak is an important device, which can make an object outside its domain invisible, meanwhile the object can exchange information with the outer region. Based on optical transformation method, we design a simplified cylindrical electromagnetic external cloak with only axial material parameter spatially variant in this paper. The general expressions of material parameters are derived, and then the performance of the external cloak is simulated using the full wave simulations. The advantage of this external cloak is that transverse material parameters are constants, which makes it easier to realize with two-dimensional metamaterials. Besides, the effects of loss and perturbations of parameters on the performance of the cloak are also investigated. This work provides a feasible way for the fabrication of the metamaterial-assisted external cloak.

©2011 Optical Society of America

## 1. Introduction

Control of electromagnetic wave with metamaterials is of great topical interest, and is fuelled by rapid progress in electromagnetic cloaks [1–6]. The approach to achieve invisibility of an object generally includes optical transformation [1,2], scattering cancellation [7], and transmission-line technique [8,9]. A review of these works has been made by Alitalo and Tretyakov [10]. Among various approaches, the optical transformation method plays an important role in the design of metamaterial devices [11–15]. Through a specific transformation, Pendry et al. [1] firstly proposed the cylindrical and spherical cloak that is capable of guiding electromagnetic wave around the cloaking region without scattering. Later, this theoretical prediction was confirmed by full-wave simulations [3] and realized in experiments at microwave and optical frequencies [4,5] with help of the metamaterials. Inspired by these pioneering works, many further investigations on cloak have been conducted, including elliptical-cylindrical cloak [16], arbitrary shaped cloak [17], line-transformed cloak [18], area-transformed cloak [19], and cloaks with a twin cavity [20]. The above cloaks can conceal an object by steering waves around an enclosed domain so that any object located inside the domain is hidden from observation. However, the objects hidden by these cloaks are “blind”, since no outside electromagnetic waves can reach into the cloaking region. It would thus be desirable to design new invisibility cloaks in which the hidden object can see the outside world.

Recently, Lai et al. [21] proposed a new recipe for invisibility cloak based on the concepts of complementary media and optical transformation. This cloak is composed of a dielectric core and an “anti-object” embedded inside a negative index shell, and it can cloak an object with a pre-specified shape and size outside the cloaking shell. Later, the generalized material parameter equations for the external cloak with arbitrary cross section were developed by our group [22]. In the foregoing investigations, however, the practical fabrication of the external cloak has so far not been achievable, for all components of material parameters are the functions of radius. Towards the practical and flexible realizations of the external cloak, a simplified cylindrical electromagnetic external cloak is designed based on optical transformation method. The radial and tangential (transverse i.e. $r$and $\phi $direction) material parameters of the designed external cloak are constants, and only axial ($z$direction) parameter is a function of radius. Thus, it is possible for the external cloak to be constructed with two-dimensional (2D) metamaterials. Full wave simulations based on the finite element method verified the performance of the designed cloak. The simplicity of the material parameters would move the 2D metamaterial external cloak a step further towards the practical realization.

## 2. Method and simulation model

The schematic diagram for the design of the external cloak is shown in Fig. 1
, where three circles with radius of${R}_{1}$,${R}_{2}$,${R}_{3}$divide the original space into three regions, i.e., $r\text{'}<{R}_{1}$, ${R}_{1}<r\text{'}<{R}_{2}$and${R}_{2}<r\text{'}<{R}_{3}$. To design the cloak, we first fold the region${R}_{2}<r\text{'}<{R}_{3}$intoregion${R}_{1}<r<{R}_{2}$, and then compress the region${R}_{1}<r\text{'}<{R}_{3}$into region$r<{R}_{1}$. Here *r* and *r’* represent the radius of the transformed space and the original space. According to the coordinate transformation method and the form invariance of the Maxwell’s equations, the relative permittivity and permeability in the transformed space can be given as [23]:

As for the material parameters in the core material layer ($r\in (0,{R}_{1})$), it can be obtained through the linear coordinate transformation of${f}_{cor}(r)=r{R}_{3}/{R}_{1}$. Then according to Eq. (1), the material parameters can be derived as:

where the subscript*cor*represents the core material. Obviously, all three components of material parameters in the core material layer are constants, which make it feasible to be realized by the ordinary anisotropic material.

Suppose a blue object with permittivity$\epsilon $and permeability$\mu $is located in the outer air layer. In order to make it invisible, we need to add a red “anti-object” with parameters$\epsilon \text{'}=\epsilon {\epsilon}_{com}$and $\mu \text{'}=\mu {\mu}_{com}$into the complementary media layer, as shown in Fig. 1(b). It should be noted that the “anti-object” is mapped into complementary media layer according to Eq. (3). The cloak is composed of the modified complementary layer embedded with the “anti-object” and a core material. Its working principle can be summarized as follows: First, the scattering of the object and surrounding space is optically cancelled by the modified complementary media layer embedded with the “anti-object”. Then, the optical path in the cancelled space is restored by the dielectric core material. In the next section, we will make full-wave simulation based on the finite element software COMSOL Multiphysics to demonstrate the designed material parameters and the performance of the external cloak.

## 3. Simulation results and discussions

In the simulation, we only consider the case of TE polarization, of which a plane wave or a cylindrical wave is impinging onto the cloak from a specific direction. As for TM polarization case, the simulation can be done in the same way. It is not included herein for brevity. The whole computational domain is surrounded by a perfectly matched layer that absorbs waves propagating outward from the bounded domain. For the 2D external cloak, the geometry parameters are chosen as${R}_{1}=0.5m$,${R}_{2}=\text{1m}$, and${R}_{3}=\text{2m}$. Based on these geometry parameters, the material parameters for the core material layer and complementary media layer of the cloak can be easily calculated through Eqs. (4) and (5), and shown in Figs. 2(a)
and 2(b), respectively. It is clear from Fig. 2(a) that all components of material parameters for the core material layer are homogeneous and greater than or equal to one. From Fig. 2(b), we can observe that only *z* component of material parameters for complementary media layer is space dependent. Although the values of material parameters are less than zero, they are accessible from the metamaterials. In what follows, the performance of the external cloak will be investigated from the following aspects.

We first demonstrate the scheme shown in Fig. 1(a), i.e. a system composed of a core material layer ($r<0.5m$) and the complementary media layer ($0.5m<r<1m$) is optically equal to a large circle of air ($r\text{'}<2m$). Figure 3(a)
shows the electric field distribution in the vicinity of this system when a plane wave with wavelength$\lambda =1$m propagates from left to right. The absence of scattered waves clearly verifies the invisibility of the whole system. Then, we demonstrate the scheme shown in Fig. 1(b), i.e. the cloaking of an object by placing its “anti-object” in the complementary layer. An circular dielectric object with radius *r* = 0.3m and parameters$\epsilon =1.5$,$\mu =1$is centered at (−1m, −1m), as shown in Fig. 3(b), which also show its scattering pattern under plane wave irradiation. In order to make it invisible, we include an “anti-object” with parameters of ${\epsilon}_{z}\text{'}=1.5{\epsilon}_{comz}\text{'}$and$\mu \text{'}={\mu}_{com}\text{'}$ into the complementary media layer, as shown in Fig. 3(c). As can be seen from Fig. 3(c), although the incident plane waves are distorted in the transformation region, they restore original wave fronts when passing through the cloak. Hence, the perfect cloaking effect is verified. It is worth noticing that the object to be cloaked is placed outside the cloaking shell, and the cloaking effect comes from its “anti-object” embedded in the complementary media layer. The white flecks in the figure represent overvalued fields which are caused by the surface mode resonance. In addition, the cloaking of a circular dielectric object with linearly changing permittivity$\epsilon =1-r/5$under cylindrical wave irradiation is also simulated, and shown in Fig. 3(d). A line source with a current of 1A/m is located at $(-\text{2}.\text{25m},-\text{2}.\text{25m})$to generate the cylindrical wave. In this case, the wavelength of cylindrical wave is set as$\lambda =0.5m$, thus the size of object to be cloaked is larger than the wavelength. It can be clearly seen that the wave fronts of the cylindrical wave are perfectly restored when the wave exits the cloak. That is to say, whether the object being cloaked is larger than the wavelength or not, the object can be perfectly cloaked when it fit into the region bounded by${R}_{2}$and${R}_{3}$, and its “anti-object” is located in the complementary layer. From above simulation results, we can conclude that arbitrary object with homogeneous or anisotropic material parameters can be cloaked when its anti-object is “custom-made” according to Eq. (3), and the performance of the cloak is independent on the excitation source and incident direction.

Next, we will illustrate the cloaking scheme with another two examples. Figure 4(a) displays the scattering pattern of two circular dielectric objects. The one on the top-left has a linearly changing permittivity of $\epsilon =1-r/5$, while the permittivity of another one on the downside is$\epsilon =1.5$. In such a case, two “anti-objects” with material parameters of ${\epsilon}_{z}\text{'}=(1-r/5){\epsilon}_{comz}\text{'}$, $\mu \text{'}={\mu}_{com}\text{'}$and${\epsilon}_{z}\text{'}=1.5{\epsilon}_{comz}\text{'}$,$\mu \text{'}={\mu}_{com}\text{'}$are embedded at the corresponding “image” positions in the complementary media layer, respectively. The perfect recovered wave fronts of the impinging plane wave shown in Fig. 4(b) demonstrate the effectiveness of the external cloak for the two objects. Actually, more objects located outside the cloak can also be cloaked, and there is no geometrical or material constraint on the object to be cloaked, as long as it fits into the region bounded by${R}_{2}$and${R}_{3}$, and its “anti-object” is located in the complementary layer. In Fig. 4(c), we show the scattering pattern of a dielectric shell with anisotropic permeability of ${\mu}_{r}=0.3$ and${\mu}_{\phi}=-0.6$. The shell is bounded between the circles of $r=1.5m$and $r=1.8m$. Under these circumstances, the “anti-object” is a complementary “image” shell with ${\mu}_{r}\text{'}=0.3{\mu}_{comr}\text{'}$, ${\mu}_{\phi}\text{'}=-0.6{\mu}_{com\phi}\text{'}$, and ${\epsilon}_{z}\text{'}={\epsilon}_{comz}\text{'}$. The electric field distribution in the vicinity of the cloak is shown in Fig. 4(d). The perfect plane wave pattern manifests the excellent cloaking effect.

Since metamaterials are always lossy in real applications, it is very necessary to study the effect of loss on the performance of the cloak. The electric field distributions in the vicinity of the external cloak with electric and magnetic loss tangents (tgδ) of 0.001, 0.01, 0.02, and 0.03 are displayed in panels (a), (b), (c) and (d) of Fig. 5 . Here, the permittivity of circular dielectric object to be cloaked is$\epsilon =1.5$. It can be seen in Figs. 5(a)–5(c) that the electric field distributions are basically undisturbed when loss tangents of 0.001, 0.01 and 0.02 are added into both permittivity and permeability of the material. But when the loss tangent of metamaterials is 0.03 or more than, it deteriorates the performance of the external cloak in the forward scattering region, as shown in Fig. 5(d). In order to give a more intuitive understanding of loss effect, electric field distributions along x axis of the external cloak with different loss tangents are shown in Fig. 6 . We can clearly observe that when loss tangents are less than or equal to 0.02, electric field distribution generally overlaps with the lossless case. Therefore, loss tangents of metamaterials being less than or equal to 0.02 is acceptable.

Finally, the robustness of external cloak with non-ideal parameters is studied. We multiply transverse material parameters for both the complementary media layer and core material by coefficient ${\delta}_{1}$, and axial material parameter by coefficient${\delta}_{2}$. It represents the slight changes in the transverse and axial parameters. When there is a slight change in the axial parameter, research on the robustness of the external cloak is carried out under three cases: (i) the transverse parameters are kept invariant, (ii) the impedance ($Z=\sqrt{{\mu}_{\phi}/{\epsilon}_{z}}$and$\sqrt{{\mu}_{r}/{\epsilon}_{z}}$) is kept invariant, (iii) the refractive index ($n=\sqrt{{\mu}_{\phi}{\epsilon}_{z}}$and$\sqrt{{\mu}_{r}{\epsilon}_{z}}$) is kept invariant. In the three cases, the corresponding changes in the transverse parameters can be easily obtained. In the simulation, the negative perturbation (${\delta}_{2}=0.9$) and positive perturbation (${\delta}_{2}=1.1$) in the axial parameter are taken into consideration, and the parameters of circular dielectric object to be cloaked is $\epsilon =1.5$and$\mu =1$. Figures 7(b) –7(d) and Figs. 8(b) –8(d) show the corresponding electric field distributions in the vicinity of the external cloak. Simulation results for the cloak without perturbation are also simulated for comparisons, as shown in Figs. 7(a) and 8(a). It is clear that whether the transverse parameters change or not, the performance of the external cloak will be affected by the positive or negative perturbation in the axial parameters. In practice, if we want to keep the cloaking effects of the external cloak while at the same time minimizing its scattering field, it is the best choice to keep the refractive index invariant.

Far field intensities of the external cloak with negative and positive perturbation in axial parameter are shown in Figs. 9(a) and 9(b), respectively. It is seen that once the refraction index is kept invariant, negative or positive perturbation in axial parameter will have little influence on the performance of the external cloak. This is in good agreement with the conclusions obtained from the near field distribution shown in Figs. 7 and 8.

## 4. Conclusion

In conclusion, a cylindrical electromagnetic external cloak with simplified material parameters is designed by means of optical transformation method. We show that axial permittivity and permeability of the external cloak is dependent on the radius, while the transverse material parameters are constants. The cloaking effect is basically undisturbed when loss tangents of metamaterials is less than or equal to 0.02. Moreover, the external cloak is robust against the perturbation in the axial material parameters when the refractive index is kept invariant. We believe that this external cloak can be constructed with 2D metamaterials, and the design method can also be extended to other transformation devices.

## Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant No. 60861002), Training Program of Yunnan Province for Middle-aged and Young Leaders of Disciplines in Science and Technology (Grant No. 2008PY031), the Scientific Research Foundation of Yunnan University (Grant No. 2010YB025), and NSFC-YN (Grant No. U1037603).

## References and links

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

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

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

**4. **D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science **314**(5801), 977–980 (2006). [CrossRef] [PubMed]

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

**6. **U. Leonhardt and T. Tyc, “Broadband invisibility by non-Euclidean cloaking,” Science **323**(5910), 110–112 (2009). [CrossRef] [PubMed]

**7. **A. Alù and N. Engheta, “Achieving transparency with plasmonic and metamaterial coatings,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. **72**(1), 016623 (2005). [CrossRef] [PubMed]

**8. **P. Alitalo, S. Ranvier, J. Vehmas, and S. Tretyakov, “A microwave transmission-line network guiding electromagnetic fields through a dense array of metallic objects,” Metamaterials (Amst.) **2**(4), 206–212 (2008). [CrossRef]

**9. **R. B. Hwang and C. Y. Chin, “Broadband cloaking using composite dielectrics,” AIP Advances **1**(1), 012112 (2011). [CrossRef]

**10. **P. Alitalo and S. Tretyakov, “Electromagnetic cloaking with metamaterials,” Mater. Today **12**(3), 22–29 (2009). [CrossRef]

**11. **W. X. Jiang, J. Y. Chin, and T. J. Cui, “Anisotropic metamaterial devices,” Mater. Today **12**(12), 26–33 (2009). [CrossRef]

**12. **D. H. Kwon and D. H. Werner, “Transformation electromagnetics: an overview of the theory and applications,” IEEE Antennas Propag. Mag. **52**(1), 24–46 (2010). [CrossRef]

**13. **J. J. Yang, M. Huang, C. F. Yang, Z. Xiao, and J. H. Peng, “Metamaterial electromagnetic concentrators with arbitrary geometries,” Opt. Express **17**(22), 19656–19661 (2009). [CrossRef] [PubMed]

**14. **M. Huang and J. J. Yang, *Wave Propagation* (Intech Press, 2011), Chap. 2.

**15. **J. J. Yang, T. H. Li, M. Huang, and M. Cheng, “Transparent device with homogeneous material parameters,” Appl. Phys., A Mater. Sci. Process. , doi:. [CrossRef]

**16. **W. X. Jiang, T. J. Cui, G. X. Yu, X. Q. Lin, Q. Cheng, and J. Y. Chin, “Arbitrarily elliptical–cylindrical invisible cloaking,” J. Phys. D Appl. Phys. **41**(8), 085504 (2008). [CrossRef]

**17. **C. Li and F. Li, “Two-dimensional electromagnetic cloaks with arbitrary geometries,” Opt. Express **16**(17), 13414–13420 (2008). [CrossRef] [PubMed]

**18. **W. X. Jiang, H. F. Ma, Q. Cheng, and T. J. Cui, “A class of line-transformed cloaks with easily realizable constitutive parameters,” J. Appl. Phys. **107**(3), 034911 (2010). [CrossRef]

**19. **X. H. Wang, S. B. Qu, X. Wu, J. F. Wang, Z. Xu, and H. Ma, “Area-transformation method for designing invisible cloaks,” J. Appl. Phys. **108**(7), 073108 (2010). [CrossRef]

**20. **T. Y. Chen and C. N. Weng, “Invisibility cloak with a twin cavity,” Opt. Express **17**(10), 8614–8620 (2009). [CrossRef] [PubMed]

**21. **Y. Lai, H. Y. Chen, Z. Q. Zhang, and C. T. Chan, “Complementary media invisibility cloak that cloaks objects at a distance outside the cloaking shell,” Phys. Rev. Lett. **102**(9), 093901 (2009). [CrossRef] [PubMed]

**22. **C. F. Yang, J. Yang, M. Huang, Z. Xiao, and J. H. Peng, “An external cloak with arbitrary cross section based on complementary medium and coordinate transformation,” Opt. Express **19**(2), 1147–1157 (2011). [CrossRef] [PubMed]

**23. **W. S. Cai, U. K. Chettiar, A. V. Kildishev, V. M. Shalaev, and G. W. Milton, “Nonmagnetic cloak with minimized scattering,” Appl. Phys. Lett. **91**(11), 111105 (2007). [CrossRef]