## Abstract

We studied the scattering from the simplified cylindrical cloaks analytically at both normal and oblique incidences. We found that these simplified cylindrical cloaks may produce a larger scattering at nonnormal incidences than that from an object without any cloak, making this object more “visible”. Even at normal incidence, the high-order transformation with impedance matched at the outer boundary can produce stronger scattering than the linear simplified one without matched impedance. This is due to the inefficiency of guided waves close to the inner boundary. Therefore, a square root transformation can improve scattering by guiding waves away from the inner boundary.

©2008 Optical Society of America

## 1. Introduction

Recently there is increased interest in studying invisibility cloaks [1–26]. Methods utilizing plasmonic and negative refractive index materials were proposed to achieve invisibility by canceling the dipole moments of small objects to be concealed [1–4]. Inspired by the concept of coordinate transformation, some other methods have also been proposed [5–7]. Especially, the transformation based models proposed in Ref. [6] has been proved to be able to conceal perfectly not only passive objects [12–14], but also active sources [16] from external electromagnetic detection. In studies of the transformation based invisibility cloaks, simplified cylindrical cloaks for normally incident waves with only one polarization are often preferred in experimental demonstration [15], simulations [9, 11] and other practical designs [17–20, 23, 24]. However, the scattering evaluation at only one incident angle is not sufficient to describe the performance of such simplified cloaks. Previously, scattering results for normal incidence were obtained only semi-analytically from commercial simulation tools [21, 22]. Therefore, a fully analytic method to evaluate the performance of simplified cylindrical cloaks at both normal and oblique incidences is needed in designing invisibility cloaks.

The main reasons for scattering from a simplified cloak are the mismatches at the outer and inner boundaries. There was study on the mismatch at the outer boundary that showed that impedance matching at the outer boundary applied by a high-order quadratic transformation in Ref. [18] is able to decrease the scattering efficiently, while there is little study on the mismatch at the inner boundary.

In this paper, the field solutions of the simplified cylindrical cloaks at both normal and oblique incidences are found. The reduced scattering from these simplified cloaks can be achieved only within a limited range of incident angles. The high-order transformation with impedance matched at the outer boundary may produce stronger scattering than the simplified linear one without impedance matched at the outer boundary due to inefficiency of guided waves close to the inner boundary.

## 2. Scattering model and analytic algorithm

The configuration of scattering from a cylindrical cloak with outer radius *R*
_{2} and inner radius *R*
_{1} follows that in Ref. [13] where a time harmonic plane wave
${\overline{E}}_{i}=\left({\hat{\nu}}_{i}{E}_{\nu i}+{\hat{h}}_{i}{E}_{\mathrm{hi}}\right){e}^{i{\overline{k}}_{i}\xb7\overline{r}}$ is incident with *k̄*
_{i}=*x̂*
*kρ*
_{+}
*ẑ*
*k*
_{z}, *k*
^{2}
_{i}=*ω*
^{2}
*µ*
_{0}
*ε*
_{0},
${\hat{h}}_{i}=\frac{\hat{z}\times {\hat{k}}_{i}}{\mid \hat{z}\times {k}_{i}\mid}$ and *v̂*
_{i}=*ĥ*
_{i}×*k̂*
_{i}. In this paper we only consider incident waves with vertical polarization, i.e. *E*
_{hi}=0. Another point here is that the cloak layer between *R*
_{1} and *R*
_{2} with permittivity tensor *ε̿*=*ε*
_{ρ}
*$\widehat{\rho}$$\widehat{\rho}$*+*ε*
_{ϕ}
*$\widehat{\varphi}$$\widehat{\varphi}$*+*ε*
_{z}
*ẑẑ* and permeability tensor *µ̿*=*µ*
_{ρ}
*$\widehat{\rho}$$\widehat{\rho}$*+*µ*
_{ϕ}
*$\widehat{\varphi}$$\widehat{\varphi}$*+*µ*
_{z}
*ẑẑ* is a simplified cylindrical cloak. Without loss of generality, we assume that the concealed region of *ρ*<*R*
_{1} is a perfect electric conductor (PEC). We denote *α*=*π*/2-*θ*
_{i} as the incident angle such that *α*=0 corresponds to the normal incidence.

In the region of *ρ*>*R*
_{2}, we use scalar potentials *ψ*
^{z}
_{TM} and *ψ*
^{z}
_{TE} to describe the TE^{z} and TM^{z} harmonic cylindrical waves [13]. In the following, we will derive the state propagator matrix in the cloak shell using state-variable approach [27]. We denote *Ē*=*Ē*
_{ρ}+*Ē*
_{s} and *H̄*=*H̄*
_{ρ}+*H̄*
_{s}, where *Ē*
_{ρ} and *H̄*
_{ρ} are components parallel to *$\widehat{\rho}$* while *Ē*
_{s} and *H̄*
_{s} are components perpendicular to *$\widehat{\rho}$*, and ∇×=∇_{s}×+∇_{ρ}×, where

Faraday’s law and Ampere’s law can be written as

$${\nabla}_{s}\times {\overline{E}}_{s}=i\omega {\mu}_{\rho}{\overline{H}}_{\rho};$$

$${\nabla}_{s}\times {\overline{H}}_{s}=-i\omega {\epsilon}_{\rho}{\overline{E}}_{\rho}.$$

Noting that for harmonic cylindrical waves of order *n*, the field has *e*
^{inϕ} and
${e}^{{\mathrm{ik}}_{z}z}$ dependencies, hence the state transfer equation can be deduced after eliminating *E*
_{ρ} and *H*
_{ρ} in Eqs. (2) and (3) as follows:

We let the state function be the four-dimensional vector of *V̄*=[*E*
_{z}
*E*
_{ϕ}
*H*
_{z}
*H*
_{ϕ}]^{T} and denote the matrix in Eq. (4) as *T̿*. We cut the cloak shell between *ρ*=*R*
_{1} and *ρ*=*R*
_{2} into *N* layers, each of which has thickness of Δ*ρ*. It follows from Eq. (4) that

Subsequently,

By choosing an sufficiently large N and matching the boundary conditions at *ρ*=*R*
_{1} and *ρ*=*R*
_{2}, we are able to solve the field solution for various cylindrical cloaks with continuously varying parameters.

## 3. Numerical results and analysis

Now we consider different transformations from virtual space *ρ*
^{′} to physical space *ρ* as follows: linear transformation as *ρ*=*R*
_{1}+(*R*
_{2}-*R*
_{1})*ρ*
^{′}/*R*
^{2} [6, 9]; quadratic transformation as *ρ*=[1-*R*
_{1}/*R*
_{2} +(*ρ*
^{′}-*R*
_{2})*R*
_{1}/*R*
^{2}
_{2}]*ρ*
^{′}+*R*
_{1} [18, 23]; and square root transformation as
$\rho ={R}_{1}+\left({R}_{2}-{R}_{1}\right)\sqrt{\frac{\rho \prime}{{R}_{2}}}$. The linear transformation allows the formation of virtual space within physical space, which is compressed linearly in the radial direction during the process. In the quadratic transformation, such compression is denser near the inner boundary. In the square root transformation, such compression is denser near the outer boundary. They can form different simplified cylindrical cloaks: (a) Linear simplified cloak [9, 15, 21] with *ε*
_{z}=*ε*
_{0}[*R*
_{2}/(*R*
_{2}-*R*
_{1}]^{2}, *µ*
_{ρ}=*µ*
_{0}[(*ρ*-*R*
_{1})/*ρ*]^{2} and *µ*
_{ϕ}=*µ*
_{0;} (b) Impedance matched linear simplified cloak [22] with *ε*
_{z}=*ε*
_{0}
*R*
_{2}/(*R*
_{2}-*R*
_{1}), *µ*
_{ρ}=*µ*
_{0}
*R*
_{2}/(*R*
_{2}-*R*
_{1})[(*ρ*-*R*
_{1})/*ρ*]^{2} and *µ*
_{ϕ}=*µ*
_{0}
*R*
_{2}/(*R*
_{2}-*R*
_{1}); (c) Impedance matched quadratic simplified cloak [18, 23] with *ε*
_{z}=*ε*
_{0}/[2*R*
_{1}
*ρ*
^{′}/*R*
^{2}
_{2}-2*R*
_{1}/*R*
_{2}+1]^{2}, *µ*
_{ρ}=*µ*
_{0}(*ρ*
^{′}/*ρ*)^{2}[2*R*
_{1}
*ρ*
^{′}/*R*
^{2}
_{2}-2*R*1/*R*
_{2}+1]^{2} and *µ*
_{ϕ}=*µ*
_{0}; and (d) Impedance matched square root simplified cloak with *ε*
_{z}=*ε*
_{0}2*R*
_{2}(*ρ*-*R*
_{1})^{2}/(*R*
_{2}-*R*
_{1})^{3}, *µ*
^{ρ}=*µ*
_{0}
*R*
_{2}(*ρ*-*R*
_{1})^{2}/[2(*R*
_{2}-*R*
_{1})*ρ*
^{2}] and *µ*
_{ϕ}=*µ*
_{02}
*R*
_{2}/(*R*
_{2}-*R*
_{1}). Note that (a) and (b) are both from linear transformation and (b), (c) and (d) have matched impedance at the outer boundary [18, 22, 23].

Figure 1 shows the electric field distribution in xy plane of different simplified cylindrical cloaks illuminated by a normally incident plane wave with vertical polarization. In region *ρ*>*R*
_{2}, only the scattered field is plotted. In Fig. 1(a), the linear simplified cloak has some intrinsic scattering [21], especially in the backward direction due to the impedance mismatch. From Fig. 1(a) to Fig. 1(b), the field pattern inside the cloak shell (the bending effect of waves) does not change much, since these cases are both from linear transformation. However, the scattered field (normalized radar cross section (RCS) *Q*
_{scat} [13]) is reduced significantly because the impedance at the outer boundary has been matched. The ring-like scattered field pattern is mainly from the zeroth order scattering [21], which is impossible to be completely eliminated by simplified cloaks due to the lack of surface magnetic current in this case [13, 26]. When the waves reach the inner boundary, instead of inducing magnetic current along *$\widehat{\varphi}$* direction in the case of ideal cylindrical cloak [13, 26], they induce electric current along *ẑ* direction on the surface of PEC core which reradiates and contributes to scattering. From Fig. 1(b) to Fig. 1(c), the impedance at the outer boundary stays matched, but the quadratic transformation in Fig. 1(c) forces most waves to be bent close to the inner boundary. Thus, more energy of the incident waves is blocked by the PEC core in Fig. 1(c) than in Fig. 1(b), forming a shadow behind the cloak characterized by a strong forward scattering as shown in Fig. 1(c). It is worth noting that the normalized RCS in Fig. 1(c) is even larger than that in Fig. 1(a). This result does not contradict with that in Ref. [18], since in the case studied in Ref. [18], the surface electric current on the PEC core happens to be helpful for complete cloaking [13, 26]. In Fig. 1(d), the square root transformation forces most waves to be bent away from the inner boundary and thus produces a very small RCS. Figure 2 shows the far-field differential normalized RCS [13] of cases (a) to (d) in Fig. 1. It can be seen that all the simplified cloaks having matched impedance at the outer boundary are able to suppress the backward scattering efficiently. But the impedance matched quadratic simplified cloak has so large forward scattering which leads to an even larger total RCS, while the impedance matched square root simplified cloak is able to suppress both backward and forward scattering, which therefore results in the smallest total RCS.

We can further calculate the scattering at nonnormal incidences. For each simplified cloak, since only *ε*
_{z}, *µ*
_{ρ} and *µ*
_{ϕ} are specified, we require that *µ*
_{z}/*µ*
_{0}=*ε*
_{z}/*ε*
_{0}, *ε*
_{ρ}/*ε*
_{0}=*µ*
_{ρ}/*µ*
_{0} and *ε*
_{ϕ}/*ε*
_{0}=*µ*
_{ϕ}/*µ*
_{0}. As shown in Fig. 3, the reduced RCS can only be achieved within a limited range of incident angles, beyond which the scattering is even larger than a bare PEC without any cloak. This result is reasonable since the simplified cloaks are all designed for only normal incidence while the scattering performance at other incident angles has not been considered. Therefore it can be expected that as the incident angle increases, the scattering will increase. But an important question is what is the critical incident angle beyond which a simplified invisibility cloak no longer makes an object “invisible” but makes it even more “visible”. Obviously, this critical incident angle is of great importance in practical applications of simplified invisibility cloaks.

## 4. Conclusion

In summary, we calculated the scattering from the simplified cylindrical cloaks at both normal and oblique incidences. We found that simplified cylindrical cloaks may produce a larger scattering at nonnormal incidences than that from an object without any cloak. The impedance matched high-order transformation forces waves to be bent close to the inner boundary and thus may produce stronger scattering than the simplified linear one without impedance matched at the outer boundary.

## Acknowledgments

This work is sponsored by ONR under Contract N00014-01-1-0713, the Department of the Air Force under Air Force Contract F19628-00-C-0002, and the Chinese NSF under Grant No. 60531020.

## References and links

**1. **A. Alu and N. Engheta, “Achiveing transparency with plasmonic and metamaterial coatings,” Phys. Rev. E **72**, 016623 (2005). [CrossRef]

**2. **G. W. Milton and N. A. P. Nicorovici, “On the cloaking effects associated with anomalous localized resonance,” Proc. R. Soc. London A **462**, 3027 (2006). [CrossRef]

**3. **A. Sihvola, “Peculiarities in the dielectric response of negative-permittivity scatterers,” Prog. Electromagn. **Res. pier-66**, 191–198 (2006). [CrossRef]

**4. **N. A. P. Nicorovici, G. W. Milton, R. C. McPhedran, and L. C. Botten, “Quasistatic cloaking of two-dimensional polarizable discrete systems by anomalous resonance,” Opt. Express **15**, 6314–6323 (2006). [CrossRef]

**5. **A. Greenleaf, M. Lassas, and G. Uhlmann, “Anisotropic conductivities that cannot be detected by EIT,” Physiol. Meas. **24**, 413–419 (2003). [CrossRef] [PubMed]

**6. **J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling Electromagnetic Fields,” Science **312**, 1780–1782 (2006). [CrossRef] [PubMed]

**7. **U. Leonhardt, “Optical Conformal Mapping,” Science **312**, 1777–1780 (2006). [CrossRef] [PubMed]

**8. **D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express **14**, 9794–9804 (2006). [CrossRef] [PubMed]

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

**10. **F. Zolla, S. Guenneau, A. Nicolet, and J. B. Pendry , “Electromagnetic analysis of cylindrical invisibility cloaks and the mirage effect,” Opt. Lett. **32**, 1069–1071 (2007). [CrossRef] [PubMed]

**11. **C. Blanchard, J. Porti, B. I. Wu, J. A. Morente, A. Salinas, and J. A. Kong, “Time domain simulation of electromagnetic cloaking structures with TLM method,” Opt. Express **16**, 6461–6470 (2008). [CrossRef] [PubMed]

**12. **H. Chen, B. I. Wu, B. Zhang, and J. A. Kong, “Electromagnetic Wave Interactions with a Metamaterial Cloak,” Phys. Rev. Lett. **99**, 063903 (2007). [CrossRef] [PubMed]

**13. **B. Zhang, H. Chen, B. I. Wu, Y. Luo, L. Ran, and J. A. Kong, “Response of a cylindrical invisibility cloak to electromagnetic waves,” Phys. Rev. B **76**, 121101(R) (2007). [CrossRef]

**14. **Z. Ruan, M. Yan, C. W. Neff, and M. Qiu, “Ideal Cylindrical Cloak: Perfect but Sensitive to Tiny Perturbations,” Phys. Rev. Lett. **99**, 113903 (2007). [CrossRef] [PubMed]

**15. **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**, 977–980 (2006). [CrossRef] [PubMed]

**16. **B. Zhang, H. Chen, B. I. Wu, and J. A. Kong, “Extraordinary Surface Voltage Effect in the Invisibility Cloak with an Active Device Inside,” Phys. Rev. Lett. **100**, 063904 (2008). [CrossRef] [PubMed]

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

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

**19. **Y. Huang, Y. Feng, and T. Jiang, “Electromagnetic cloaking by layered structure of homogeneous isotropic materials,” Opt. Express **15**, 11133–11141 (2007). [CrossRef] [PubMed]

**20. **H. Chen, Z. Liang, P. Yao, X. Jiang, H. Ma, and C. T. Chan, “Extending the bandwidth of electromagnetic cloaks,” Phys. Rev. B **76**, 241104(R) (2007). [CrossRef]

**21. **M. Yan, Z. Chao, and M. Qiu, “Cylindrical Invisibility Cloak with Simplified Material Parameters is Inherently Visible,” Phys. Rev. Lett. **99**, 233901 (2007). [CrossRef]

**22. **M. Yan, Z. Chao, and M. Qiu, “Scattering characteristics of simplified cylindrical invisibility cloaks,” Opt. Express **15**, 17772–17782 (2007). [CrossRef] [PubMed]

**23. **W. Cai, U. K. Chettiar, A. V. Kildishev, and V. M. Shalaev, “Designs for optical cloaking with high-order transformations,” Opt. Express **16**, 5444–5452 (2008). [CrossRef] [PubMed]

**24. **M. Farhat, S. Guenneau, A. B. Movchan, and S. Enoch, “Achieving invisibility over a finite range of frequencies,” Opt. Express **16**, 5656–5661 (2008). [CrossRef] [PubMed]

**25. **S. Xi, H. Chen, B. I. Wu, B. Zhang, J. Huangfu, D. Wang, and J. A. Kong, “Effects of different transformations on the performance of cylindrical cloaks,” J. Electrom. Waves and Appl. **22**, 1489–1497 (2008). [CrossRef]

**26. **A. Greenleaf, Y. Kurylev, M. Lassas, and G. Uhlmann, “Improvement of cylindrical cloaking with the SHS lining,” Opt. Express **15**, 12717–12734 (2007). [CrossRef] [PubMed]

**27. ** W. C. Chew, *Waves and Fields in inhomogeneous Media*, 2nd ed, (IEEE Press, 1995).