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

Isotropic nonmagnetic flat cloaks degenerated from homogeneous anisotropic trapeziform cloaks

Open Access Open Access

Abstract

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

1. 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, one-directional 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 [112], 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.

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

 figure: Fig. 1

Fig. 1 The scheme of coordinate transformation from the virtual space (a) to physical space (b). k1 represents the slope of the inner PEC boundary in Region I (identical to that in Region II).

Download Full Size | PDF

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]

ε=μ=AAT/det(A) 
where A=(x,y,z)/(x,y,z) 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. 1(b). The transformation equations for region I are

x=x',  y=k1(x+d)+kry+h,  z=z'
where kr=b/(b+h). Then we can obtain the parameters of region I from Eqs. (1-2)
εI=μI=1kr[1k10k1k12+kr20001] .
Similarly, the transformation equations for region II become
x=x,  y=k1(dx)+kry+h,  z=z' .
Based on Eq. (4), the parameters of region II can be expressed

εII=μII=1kr[1k10k1k12+kr20001] .

The coordinate transformations of region I and II compress two triangular regions into another two triangular regions, and the coordinate transformation of region III deforms one rectangular region into another in accordance with the deformations of the triangles on both sides, which can be expressed as

x=x',  y=kry+h,  z=z'.
From Eq. (6), the parameters of region III become

εIII=μIII=diag[b+hb,bb+h,b+hb] .

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.

3. 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 k1=1 and d=b=h=0.5  m 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 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.

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

Download Full Size | PDF

Assuming d0, 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.

Assuming k10 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]. Figure 2(d) presents the electric field distribution of the flat cloak. Compared Fig. 2(a) to Fig. 2(d), it can be seen that the directions of reflected waves are identical. Obviously, the flat cloak maintains good performance as long as the width of the cloak is large enough. Hence, the flat cloak with uniaxial materials could be extensively utilized to cloak groups of static objects in the cloaked space from y=0 to y=h beneath the PEC.

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[εx,εy,εx]. According to effective medium theory, the flat cloak is equally discretized into M layers, and each layer consists of medium-A and medium-B. The old set of material parameters of isotropic medium-A and medium-B is represented by

εA=μA=εx+εx2εxεy;  εB=μB=εxεx2εxεy .
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

εA=μB=εx+εx2εxεy;  εB=μA=εxεx2εxεy.

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.

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

Download Full Size | PDF

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

 figure: Fig. 4

Fig. 4 Comparison of the average power flow of the isotropic flat cloak in Fig. 3(c) and Fig. 3(f), mimicked by old [Eq. (8)] and new [Eq. (9)] effective medium theories. (a) M = 40, (b) M = 60. The same incidence is used as Fig. 2.

Download Full Size | PDF

 figure: Fig. 5

Fig. 5 The relative index profiles of the isotropic flat cloak in Fig. 2(d) mimicked by old [Eq. (8)] and new [Eq. (9)] effective medium theories with M = 5.

Download Full Size | PDF

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μzand εyμz unchanged, the permittivity and permeability of the flat cloak can be simplified as εIII=diag[(b+hb)2,1,1] and μIII=diag[1,1,1], 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 h=0.2  m and b=0.8  mis 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.

 figure: Fig. 6

Fig. 6 (Color online) Magnetic field distributions for a Gaussian beam incident from the right at 45on (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.

Download Full Size | PDF

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

Acknowledgement

The authors are grateful for the support from National University of Singapore under the Grant No. R-263-000-574-133.

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. 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). [CrossRef]  

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. [CrossRef]  

4. T. C. Han, X. H. Tang, and F. Xiao, “The petal-shape cloak,” J. Electromagn. Waves Appl. 23(14), 2055–2062 (2009). [CrossRef]  

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. [CrossRef]   [PubMed]  

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-17-15-13070. [CrossRef]   [PubMed]  

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). [CrossRef]   [PubMed]  

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). [CrossRef]   [PubMed]  

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). [CrossRef]  

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). [CrossRef]   [PubMed]  

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). [CrossRef]  

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). [CrossRef]  

13. J. Li and J. B. Pendry, “Hiding under the carpet: a new strategy for cloaking,” Phys. Rev. Lett. 101(20), 203901 (2008). [CrossRef]   [PubMed]  

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). [CrossRef]   [PubMed]  

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). [CrossRef]   [PubMed]  

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. [CrossRef]   [PubMed]  

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 (6)

Fig. 1
Fig. 1 The scheme of coordinate transformation from the virtual space (a) to physical space (b). k 1 represents the slope of the inner PEC boundary in Region I (identical to that in Region II).
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 .
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.
Fig. 4
Fig. 4 Comparison of the average power flow of the isotropic flat cloak in Fig. 3(c) and Fig. 3(f), mimicked by old [Eq. (8)] and new [Eq. (9)] effective medium theories. (a) M = 40, (b) M = 60. The same incidence is used as Fig. 2.
Fig. 5
Fig. 5 The relative index profiles of the isotropic flat cloak in Fig. 2(d) mimicked by old [Eq. (8)] and new [Eq. (9)] effective medium theories with M = 5.
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.

Equations (9)

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

ε = μ = A A T / det ( A )  
x = x ' ,    y = k 1 ( x + d ) + k r y + h ,    z = z '
ε I = μ I = 1 k r [ 1 k 1 0 k 1 k 1 2 + k r 2 0 0 0 1 ]
x = x ,    y = k 1 ( d x ) + k r y + h ,    z = z '
ε II = μ II = 1 k r [ 1 k 1 0 k 1 k 1 2 + k r 2 0 0 0 1 ]
x = x ' ,    y = k r y + h ,    z = z '
ε III = μ III = diag [ b + h b , b b + h , b + h b ]
ε A = μ A = ε x + ε x 2 ε x ε y ;    ε B = μ B = ε x ε x 2 ε x ε y
ε A = μ B = ε x + ε x 2 ε x ε y ;    ε B = μ A = ε x ε x 2 ε x ε y
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.