Isotropic nonmagnetic flat cloaks degenerated from homogeneous anisotropic trapeziform cloaks

We propose a novel kind of trapeziform cloak requiring only homogeneous anisotropic materials. Large-scale flat cloaks can be degenerated from the general trapeziform cloak with PEC inner boundary, and be realized by isotropic nonmagnetic materials for optical frequencies with controlled index profiles and improved invisibility. With the support of PEC inner boundary, large vehicles and objects of arbitrary shape can be concealed between the PEC and ground, and PEC can be firm by adding pillars in the cloaking space. Full-wave simulations validate the proposed cloaking concept, which is not only based on simple isotropic nonmagnetic materials but also realizable in practice. ©2010 Optical Society of America OCIS codes: (160.2710) Inhomogeneous optical materials: Anisotropic optical materials; (290.5839) Scattering, invisibility. References and links 1. J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science 312(5781), 1780–1782 (2006). 2. W. X. Jiang, T. J. Cui, X. M. Yang, Q. Cheng, R. P. Liu, and D. R. Smith, “Invisibility cloak without singularity,” Appl. Phys. Lett. 93(19), 194102 (2008). 3. M. Schmiele, C. Rockstuhl, and F. Lederer, “Strategy for cloaking of twisted domains,” Phys. Rev. A 79(5), 053854 (2009), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-3-9-332. 4. T. C. Han, X. H. Tang, and F. Xiao, “The petal-shape cloak,” J. Electromagn. Waves Appl. 23(14), 2055–2062 (2009). 5. Y. You, G. W. Kattawar, and P. Yang, “Invisibility cloaks for toroids,” Opt. Express 17(8), 6591–6599 (2009), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-17-8-6591. 6. J. Hu, X. Zhou, and G. Hu, “Design method for electromagnetic cloak with arbitrary shapes based on Laplace’s equation,” Opt. Express 17(3), 1308–1320 (2009), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-1715-13070. 7. H. S. Chen, B. I. Wu, B. Zhang, and J. A. Kong, “Electromagnetic wave interactions with a metamaterial cloak,” Phys. Rev. Lett. 99(6), 063903 (2007). 8. C.-W. Qiu, A. Novitsky, H. Ma, and S. Qu, “Electromagnetic interaction of arbitrary radial-dependent anisotropic spheres and improved invisibility for nonlinear-transformation-based cloaks,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 80(1), 016604 (2009). 9. A. Novitsky, C.-W. Qiu, and S. Zouhdi, “Transformation-based spherical cloaks designed by an implicit transformation-independent method: Theory and optimization,” N. J. Phys. 11(11), 113001 (2009). 10. M. Rahm, S. A. Cummer, D. Schurig, J. B. Pendry, and D. R. Smith, “Optical design of reflectionless complex media by finite embedded coordinate transformations,” Phys. Rev. Lett. 100(6), 063903 (2008). 11. H. Y. Chen, and C. T. Chan, “Electromagnetic wave manipulation by layered systems using the transformation media concept,” Phys. Rev. B 78(5), 054204 (2008). 12. Y. Luo, J. Zhang, H. Chen, L. Ran, B.-I. Wu, and J. A. Kong, “A rigorous analysis of plane-transformed invisibility cloaks,” IEEE Trans. Antenn. Propag. 57(12), 3926–3933 (2009). 13. J. Li, and J. B. Pendry, “Hiding under the carpet: a new strategy for cloaking,” Phys. Rev. Lett. 101(20), 203901 (2008). 14. R. Liu, C. Ji, J. J. Mock, J. Y. Chin, T. J. Cui, and D. R. Smith, “Broadband ground-plane cloak,” Science 323(5912), 366–369 (2009). 15. C.-W. Qiu, L. Hu, X. Xu, and Y. Feng, “Spherical cloaking with homogeneous isotropic multilayered structures,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 79(4), 047602 (2009). #127030 $15.00 USD Received 15 Apr 2010; revised 29 May 2010; accepted 29 May 2010; published 2 Jun 2010 (C) 2010 OSA 7 June 2010 / Vol. 18, No. 12 / OPTICS EXPRESS 13038 16. C.-W. Qiu, L. Hu, B. Zhang, B. I. Wu, S. G. Johnson, and J. D. Joannopoulos, “Spherical cloaking using nonlinear transformations for improved segmentation into concentric isotropic coatings,” Opt. Express 17(16), 13467–13478 (2009), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-17-16-13467.


Introduction
In the past few years, invisibility cloaks have attracted great attention due to their exciting properties of deflecting incoming electromagnetic (EM) waves and directing them to propagate around the cloaked objects without any scatterings.Coordinate transformation was proposed by Pendry et al [1] to deform an empty space into a shell region, which makes the object inside the shell invisible to incident EM waves.According to this concept, many different kinds of cloaks have been proposed, including elliptical cloak without singularity [2], cloaks of twisted domains [3], petal-shaped cloak [4], toroids cloak [5], and arbitrary shaped cloaks [6].Analytical methods based on Maxwell's equations have been used to reveal the nature of invisibility cloaks and characterize the wave interaction [7,8].More recently, the inverse mechanism of cloaking designs without specifying the needed coordinate transformation beforehand has been developed [9].Different from these aforementioned approaches, in which the coordinate transformation is operated along the radial direction, onedirectional transformation is proposed [10] to manipulate EM wave by layered systems [11].Recently, a conical cloak is investigated in detail using plane transformation [12].Different from the complete cloaks [1][2][3][4][5][6][7][8][9][10][11][12], a "carpet cloak" has been proposed to mimic a flat ground plane [13] and has been experimentally verified [14].Such a carpet cloak could avoid both material and geometry singularities.However, the material for the carpet cloak is still inhomogeneous and needs to be carefully designed using quasi-conformal transformation.
In this paper, we propose a general trapeziform cloak that can be realized by isotropic nonmagnetic materials.The cloak could degenerate into a triangular cloak as well as flat cloak by particularly setting the geometry of the trapeziform cloak.Note that the degenerated flat cloak only requires homogeneous uniaxial medium, when the size of the cloaked object is much smaller than the cloak or the trapeziform's slopes tend to zero.The flat cloak can be a large-scale device and easily realized due to the simplicity in the parameters of the medium above the PEC.Meanwhile, the position of the PEC cover can be adjusted based on the practical sizes of military targets which can be perfectly concealed and protected from the changes in external environments.Moreover, since the space below the PEC cover is blind, one can make the PEC cover thicker and use hard pillars between the PEC and ground as firm supports.The proposed isotropic trapeziform cloaks and its degenerated cloaks suggest insightful scenarios for exciting applications, and this concept may nurture further degrees of freedom in the design of large-scale cloaks.

Design scheme and analysis
A schematic diagram illustrating a trapeziform cloak is shown in Fig. 1.The region within the dashed trapezoid in the virtual space [Fig.1(a)] is compressed into two triangles and one rectangle regions in the physical space [Fig.1(b)] with inner PEC boundaries, resulting in a trapeziform region blind to external illumination.Thus the spatial deformation is only along y axis in a uniformly linear manner within the cloak area.
Due to the form-invariance property of Maxwell's equations throughout the transformation, the permittivity and permeability tensors of the medium in the physical space can be expressed [1] T det( ) where ( , , ) / ( , , ) is the Jacobian transformation matrix between the virtual space ( , , ) x y z ′ ′ ′ and physical space ( , , ) x y z .It should be noted that the permittivity and permeability are relative parameters throughout.
We can find that the trapeziform cloak is composed of three regions in Fig.
Similarly, the transformation equations for region II become Based on Eq. ( 4), the parameters of region II can be expressed From Eq. ( 6), the parameters of region III become It is interesting to note that the constitutive parameters in Eqs. ( 3), ( 5) and ( 7) are spatially invariant, which facilitates the experimental verification by periodic metamaterial structures and could also be easily realized through an alternating layered system of two homogeneous isotropic materials.

Simulation and representative results
In order to verify the designed formulae, we make full-wave simulations based on the finite element method (FEM).The working frequency is set to be 2 GHz under TE polarization.We create a trapeziform cloak with 1    The flat cloak in Fig. 2(d) composed of uniaxial materials is less demanding on the requirement of parameters and design complexity compared to [14] and it could be easily realized through an alternating layered isotropic mediums [11,15,16].The permittivity and permeability of the flat cloak could be expressed as III III diag[ , , ] ; However, an improved set of effective mediums is based on the index matching, to further lower down the total cross section of spherical cloaks.The new set can be expressed as Using old and new sets of material parameters respectively, the flat cloak with homogeneous uniaxial materials could be easily imitated by multilayered homogeneous isotropic materials.The electric field distributions of the flat cloak with the old (left column) and new (right column) sets of isotropic materials are shown in Fig. 3, respectively.It is obvious that the cloaking effect is quite pronounced when the discretization increases.8)] and new [Eq.( 9)] effective medium theories with M = 5.
We simulated the average power flow at the observation circle with the radius of 1.8m whose origin is at (2m, 0).In Fig. 3(c), from the angles 0 to 35 , the dashed circle is located below and inside the flat cloak, so the observation angle starts from 35 .In Fig. 4(a), it is obvious that the flat cloak with the new set of isotropic parameters behaves much better than that with the old set of isotropic material parameters.The reflective beam by the flat cloak with the new set is very close to the ideal case, while the old set still leads to measurable deviation even though the layer number M is up to 60 [Fig.4(b)].It can be attributed to the fact that the relative index of the new set of isotropic materials equal to 1 resulting in the improved invisibility, which is less oscillative and varied more smoothly than the old set of isotropic materials with step-index profile, as shown in Fig. 5.Although the impedance mismatch of the new set will lead to slightly increased backscattering, the index matching will bring down the forward scattering more significantly, resulting in improved overall scattering reduction.This feature in turn provides us more degrees of freedom to make the trade-offs between perfect invisibility and fabrication cost-effectiveness.Then near-perfect flat cloak with restored invisibility performance based on the new set of isotropic parameters can be an economic alternative, especially when the coating process of layered isotropic dielectrics has to be terminated at a certain medium value of M (e.g., M = 40 or even smaller).The examples discussed above are with complete EM parameters, which can be realized by metamaterial.For application at optical wavelength, we should reduce the material parameters to a nonmagnetic form.By keeping the products of x z ε µ and y z ε µ unchanged, the permittivity and permeability of the flat cloak can be simplified as , respectively.Figure 6 shows the magnetic field distributions for four different cases.Figure 6(a) corresponds to the magnetic field of a flat reflective mirror.When a triangular PEC is placed on the flat boundary, the magnetic field scattered by the obstacle is quiet irregular shown in Fig. 6(b).In order to make the obstacle invisible, a nonmagnetic flat cloak with 0.2 m h = and 0.8 m b = is placed above the obstacle, as shown in Fig. 6(c).Clearly, the reflective magnetic field is confined highly in the specular direction, which is in consistence with the case of Fig. 6(a).Based on the effective medium theory, the nonmagnetic flat cloak could be realized through an alternating layered isotropic mediums with A 0.625 ε = and B 2.5 ε = , as shown in Fig. 6(d).It can be seen that the isotropic nonmagnetic flat cloak perfectly mimics the case of a flat reflective mirror [Fig.6(a)].In order to design a realizable isotropic nonmagnetic flat cloak, assuming the cloak is embedded in a composite background medium with permittivity of 1.6, the alternating isotropic dielectrics proportionally become air ( A 1 ε = ) and silicon ( B 4 ε = ), respectively.

Conclusion
In this paper, a trapeziform cloak that only requires homogeneous anisotropic material is proposed, which can degenerate into flat or triangular.The flat cloak is constructed by isotropic mediums based on new set of effective medium theory, which furthermore improves the performance of actual flat cloak.We also design an isotropic nonmagnetic flat cloak that is able to work at optical frequency and could be easily realized through alternating layered dielectrics.This general cloaking concept may pave a realizable way to the practical applications in large scales and in higher frequency.

Fig. 1 . 1 k
Fig. 1.The scheme of coordinate transformation from the virtual space (a) to physical space (b). 1 k represents the slope of the inner PEC boundary in Region I (identical to that in Region 1 k = and 0.5 m d b h = = = illuminated by a Gaussian beam incident from the right at 45 as in Fig. 2. Figure 2(a) demonstrates the electric field distribution of a perfectly reflective flat boundary corresponding to Fig. 1(a), and Fig. 2(b) shows the electric field distribution of a trapeziform cloak with PEC inner boundary #127030 -$15.00USD Received 15 Apr 2010; revised 29 May 2010; accepted 29 May 2010; published 2 Jun 2010 (C) 2010 OSA corresponding to Fig. 1(b).Comparing Figs.2(a) to 2(b), the fields remain unperturbed in the presence of the cloak, as if no scatterers were present above the PEC ground.Assuming 0 d → , region III in Fig. 1(b) vanishes and the trapeziform cloak then degenerates into a triangular cloak as Fig. 2(c) illustrates.The distribution of the reflected beam by the triangular cloak demonstrates the perfect cloaking performance.

Fig. 2 .
Fig. 2. Snapshot of the electric field for a Gaussian beam incident from the right at an oblique direction of 45 .(a) A perfectly reflective flat boundary.(b) A trapeziform cloak with inner PEC boundaries as in Fig. 1(b).(c) The triangular cloak degenerated from a trapeziform cloak when Region III vanishes.(d) The degenerate flat cloak when b h = .Assuming 1 0 k → while h and b being equal and invariant, the trapeziform cloak degenerates into a flat cloak only composed of region III III III diag[2, 1/ 2, 2]

Fig. 3 .
Fig. 3. Electric field distributions of the flat cloak with M-layered isotropic materials based on two sets of effective medium theories.(a-c) correspond to the old set [Eq. (8)] while (d-f) correspond to the new set [Eq. (9)].(a, d) M = 20.(b, e) M = 40.(c, f) M = 60.The dashed semi-circle is the observation radius for measuring the average power flow for the reflective beam.
effective medium theory, the flat cloak is equally discretized into M layers, and each layer #127030 -$15.00USD Received 15 Apr 2010; revised 29 May 2010; accepted 29 May 2010; published 2 Jun 2010 (C) 2010 OSA consists of medium-A and medium-B.The old set of material parameters of isotropic medium-A and medium-B is represented by

Fig. 6 .
Fig. 6. (Color online) Magnetic field distributions for a Gaussian beam incident from the right at 45 on (a) a perfectly reflective flat boundary, (b) a bumped reflective boundary, (c) a nonmagnetic flat cloak placed above the bumped boundary, (d) an isotropic nonmagnetic flat cloak placed above the bumped boundary with M = 40.