Monostatic lidar/radar invisibility using coated spheres

doi:10.1364/OE.16.001431

« Show journal navigation

Monostatic lidar/radar invisibility using coated spheres

Peng-Wang Zhai, Yu You, George W. Kattawar, and Ping Yang »View Author Affiliations

Optics Express, Vol. 16, Issue 3, pp. 1431-1439 (2008)
http://dx.doi.org/10.1364/OE.16.001431

View Full Text Article

Acrobat PDF (104 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. Theory and formulas
3. Numerical results
4. Conclusions
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (21)
Cited By
Figures (3)
Metrics

Abstract

The Lorenz-Mie theory is revisited to explicitly include materials whose permeability is different from unity. The expansion coefficients of the scattered field are given for light scattering by both homogeneous and coated spheres. It is shown that the backscatter is exactly zero if the impedance of the spherical particles is equal to the intrinsic impedance of the surrounding medium. If spherical particles are sufficiently large, the zero backscatter can be explained as impedance matching using the asymptotic expression for the radar backscattering cross section. In the case of a coated sphere, the shell can be regarded as a cloak if the product of the thickness and the imaginary part of the refractive index of the outer shell is large.

1. Introduction

Electromagnetic cloaking is an active research topic that has recently drawn a great deal of attention [1

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

, 2

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

, 3

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]

]. A coordinate transformation approach [1

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

, 4

A. J. Ward and J. B. Pendry, “Refraction and geometry in maxwells equations,” J. Mod. Opt. 43, 773–793 (1996). [CrossRef]

] and a conformal mapping approach [2

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

] have been employed to generate a cloak with complex material properties which makes the particle invisible when viewed externally. A two dimensional cylindrical cloak was demonstrated by Schurig et al. [3

]. Schurig et al. simulated this invisible cloak in the geometric optics limit [5

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]

]. Cummer et al. [6

S. A. Cummer, B.-I. Popa, D. Schurig, D. R. Smith, and J. B. Pendry, “Full-wave simulations of electromagnetic cloaking structures,” Phys. Rev. E 74, 036621-1–5 (2006).

] and Zolla et al. [7

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]

] used “COSMOL Multiphysics,” a commercial software package based on a finite element-based electromagnetics solver to study a full-wave 2D cloaking system. Chen et al. [8

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

] reported an analytical solution to the interaction of electromagnetic wave with a 3D metamaterial cloak. Zhang et al. [9

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

] have demonstrated that the incident electromagnetic waves induced magnetic and electric surface currents at the inner boundary of a 2D cylindrical cloak. Other work related to invisible cloaks used plasmonic resonance to achieve transparency [10

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

, 11

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

In Ref. [8

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

], the solution is based on expanding the electromagnetic waves in terms of vector spherical harmonics. The expansion coefficients are determined from the boundary conditions. This approach is essentially a generalization of the Lorenz-Mie theory [12

G. Mie, “Beigrade zur optik truber medien, speziell kolloidaler metallosungen,” Ann. Phys. (Leipzig) 25, 377–455 (1908). [CrossRef]

]. The expansion coefficients of the scattered wave turn out to be zero if the following permittivity and permeability are used [1

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

\overset{=}{ε} = ε_{r} (r) \hat{r} \hat{r} + ε_{t} (\hat{θ} \hat{θ} + \hat{ϕ} \hat{ϕ}),

\overset{=}{μ} = μ_{r} (r) \hat{r} \hat{r} + μ_{t} (\hat{θ} \hat{θ} + \hat{ϕ} \hat{ϕ}),

(1)

where ε_t /ε ₀=R ₂/(R ₂-R ₁), ε_r =ε_t (r-R ₁)²/r ², μ_t /μ ₀=ε_t /ε ₀, and μ_r /μ_t =ε_r /ε_t ; R ₁ and R ₂ are the inner and outer radii of the two concentric spheres, respectively; ε ₀ and μ ₀ are the permittivity and permeability of the surrounding medium; r̂, $\hat{θ}$ , and $\hat{ϕ}$ are the unit vectors of the spherical coordinates. It is worthy noting that the medium defined by Eq. (1) is lossless. If loss exists, the cloak becomes visible and only the exact backscatter is zero. Chen et al. [8

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

] articulated that “This unique property of the spherical cloak indicates the cloaked object can still completely hide from monostatic lidar/radar detection.”

We studied the case carefully and found that zero backscatter is not a unique property of the spherical cloak represented by Eq. (1). In fact, one only needs to require η ₁=η ₀, which is also called the impedance matching condition, for a homogeneous spherical particle to achieve zero backscatter, where

η_{i} = \sqrt{\frac{μ_{i}}{ε_{i}}}, i = 1, 0,

(2)

is the impedance of the spherical particle; subscripts 1 and 0 denote the quantities for the sphere and surrounding medium, respectively. For a coated sphere, the backscatter is also zero if the impedance matching condition η_i =η ₀ is satisfied for each i^th layer of the coated sphere. Furthermore, one can speculate from physical intuition that the backscatter of a coated sphere is still approximately zero even if the impedance matching condition is violated for the inner core, given that its outer layer satisfies the condition and the imaginary part of its refractive index is sufficiently large. This is because the electromagnetic wave cannot penetrate into the inner region. Thus no matter what is inside, as long as it is a passive medium, the scattering pattern is essentially the same. This speculation is proved by numerical study and it can be used to make a cloak invisible for monostatic lidar/radar detection. The condition of this type of cloak is not as strict as that suggested in Ref. [1

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

] if only backscatter is concerned. Materials with impedance matching the surrounding medium can be made by using the Maxwell-Garnett mixing rule [13

J. Maxwell Garnett, “Colours in metal glasses and in metallic films,” Philos. Trans. R. Soc. London 203, 385 (1904). [CrossRef]

] for compounds from known material properties [14

J.R. Liu, M. Itoh, and K.-I Machida, “Frequency dispersion of complex permeability and permittivity on ironbased nanocomposites derived from rare earth-iron intermetallic compounds,” J. Alloys Compd. 408–412, 1396–1399 (2006). [CrossRef]

2. Theory and formulas

The Lorenz-Mie theory is treated in detail in several classical texts [15

H.C. van de Hulst, Light Scattering by Small Particles (Dover, New York, 1981).

, 16

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).

]. For a homogeneous sphere, we define:

x = ka, m = \sqrt{\frac{ε_{1} μ_{1}}{ε_{0} μ_{0}}},

(3)

where x is the size parameter of the sphere; k is the magnitude of the wave vector in the surrounding medium; a is the radius of the sphere; m is the relative refractive index of the sphere; ε is the permittivity; μ is the permeability. Consider the scattering of a x-polarized electromagnetic plane wave by this sphere. The propagation direction of the incident wave is along the z axis. The incident, scattered, and internal fields can be expanded in terms of the vector spherical harmonics. The expansion coefficients are determined by the boundary conditions. Eqs. (4.53) in Ref. [16

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).

] gave the general expressions of the expansion coefficients of the scattered wave. Bohren and Huffman [16

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).

] assumed μ ₁/μ ₀=1 and obtained Eqs. (4.88) in Ref. [16

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).

] for numerical calculations. To derive the corresponding equations of Eqs. (4.88) in Ref. [16

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).

] for the case of μ ₁/μ ₀≠1, the relative impedance η is introduced:

η = \frac{η_{1}}{η_{0}} = \frac{\sqrt{\frac{μ_{1}}{ε_{1}}}}{\sqrt{\frac{μ_{0}}{ε_{0}}}} .

(4)

Following steps similar to those in Ref. [16

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).

], the expansion coefficients of scattered field can be written as:

a_{n} = \frac{[D_{n} (m x) η + \frac{n}{x}] ψ_{n} (x) - ψ_{n - 1} (x)}{[D_{n} (m x) η + \frac{n}{x}] ξ_{n} (x) - ξ_{n - 1} (x)},

b_{n} = \frac{[\frac{D_{n} (m x)}{η} + \frac{n}{x}] ψ_{n} (x) - ψ_{n - 1} (x)}{[\frac{D_{n} (m x)}{η} + \frac{n}{x}] ξ_{n} (x) - ξ_{n - 1} (x)},

(5)

where the standard nomenclature for the special functions involved in the Lorenz-Mie theory is adopted [16

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).

]. Specifically, ψ_n (ρ)=ρj_n (ρ) and ξ_n (ρ)=ρh ⁽¹⁾ _n (ρ) are the Riccati-Bessel functions; j_n and h ⁽¹⁾ _n =j_n +iy_n are the spherical Bessel functions of the first and third kinds, respectively; y_n is the spherical Bessel function of the second kind. D_n (ρ)=dlnψ_n (ρ)/dρ is the logarithmic derivative of the Riccati-Bessel function ψ_n .

An immediate observation of Eq. (5) is that a_n =b_n if the relative impedance η=1. From Eqs. (4.74) in Ref. [16

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).

], the complex amplitude scattering matrix elements can be written as:

S_{1} = S_{2} = \sum_{n} \frac{2 n + 1}{n (n + 1)} a_{n} (π_{n} + τ_{n}),

(6)

which are equal to zero at θ=180° since π_n (180°)+τ_n (180°)=0 for all n. The radar backscattering cross section σ_b is defined as the area which is 4π times the differential scattering cross section for scattering into a unit solid angle around the backscattering direction [16

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).

]. It is then evident that σ_b is equal to zero for η=1. For large absorbing spherical particles, the asymptotic expression of the radar back scattering efficiency Q_b =σ_b /π a ² has a value of (Eq. 4.83 in Ref. [16

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).

])

\lim_{x \to \infty} Q_{b} = R (0^{o}),

(7)

where R(0°) is the reflectance at normal incidence at the surface of the sphere. R(0°) vanishes if the impedance matching is satisfied, namely, η=1 (see pp. 306, Ref. [17

J. D. Jackson, Classical Electrodynamics, 3rd Edition (Wiley-VCH, 1998),

]).

For spherical particles, the scattering matrix can be written in terms of the amplitude scattering matrix as follows (see Eqs. (4.77) in Ref. [16

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).

]):

\begin{matrix} S_{11} = \frac{{∣ S_{2} ∣}^{2} + {∣ S_{1} ∣}^{2}}{2}, & S_{12} = \frac{{∣ S_{2} ∣}^{2} - {∣ S_{1} ∣}^{2}}{2}, \\ S_{33} = \frac{S_{2}^{*} S_{1} + S_{2} S_{1}^{*}}{2}, & S_{34} = \frac{{i (S_{1} S_{2}^{*} - S_{2} S_{1}^{*})}^{}}{2} . \end{matrix}

(8)

Eqs. (8) shows S ₁₂=S ₃₄=0 and S 11=S ₃₃ for particles with η=1. Only diagonal elements survive for this case.

Light scattering by a coated sphere is also treated in Ref. [16

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).

]. Eqs. (8.1) in Ref. [16

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).

] give the relations for the eight expansion coefficients of the fields. The permeability μ=1 is then assumed for the particle and surrounding medium to obtain a_n and b_n . To reinstate the dependence of μ in a_n and b_n , the following definitions are convenient:

x = ka, m_{1} = \sqrt{\frac{ε_{1} μ_{1}}{ε_{3} μ_{3}}}, η_{1} = \frac{\sqrt{\frac{μ_{1}}{ε_{1}}}}{\sqrt{\frac{μ_{3}}{ε_{3}}}},

y = kb, m_{2} = \sqrt{\frac{ε_{2} μ_{2}}{ε_{3} μ_{3}}}, η_{2} = \frac{\sqrt{\frac{μ_{2}}{ε_{2}}}}{\sqrt{\frac{μ_{3}}{ε_{3}}}},

(9)

where a and b are the inner and outer radii of the coated sphere, respectively; m_i, i=1,2 are the relative refractive indices in the region i; the subscripts i=1,2 denote that the quantities are for the regions of 0<r<a and a<r<b, respectively; ε3 and μ3 are the permittivity and permeability of the surrounding medium, respectively. The expansion coefficients of the scattered field can be solved from Eqs. (8.1) in Ref. [16

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).

]; given by:

a_{n} = \frac{[{η_{2} \tilde{D}}_{n} + \frac{n}{y}] ψ_{n} (y) - ψ_{n - 1} (y)}{[η_{2} {\tilde{D}}_{n} + \frac{n}{y}] ξ_{n} (y) - ξ_{n - 1} (y)},

b_{n} = \frac{[\frac{{\tilde{G}}_{n}}{η_{2}} + \frac{n}{y}] ψ_{n} (y) - ψ_{n - 1} (y)}{[\frac{{\tilde{G}}_{n}}{η_{2}} + \frac{n}{y}] ξ_{n} (y) - ξ_{n - 1} (y)},

(10)

where

{\tilde{D}}_{n} = \frac{\frac{D_{n} (m_{2} y) - A_{n} χ_{n}^{'} (m_{2} y)}{ψ_{n} (m_{2} y)}}{\frac{1 - A_{n} χ_{n} (m_{2} y)}{ψ_{n} (m_{2} y)}},

{\tilde{G}}_{n} = \frac{\frac{D_{n} (m_{2} y) - B_{n} χ_{n}^{'} (m_{2} y)}{ψ_{n} (m_{2} y)}}{\frac{1 - B_{n} χ_{n} (m_{2} y)}{ψ_{n} (m_{2} y)}},

A_{n} = ψ_{n} (m_{2} x) \frac{η_{1} D_{n} (m_{1} x) - η_{2} D_{n} (m_{2} x)}{η_{1} D_{n} (m_{1} x) χ_{n} (m_{2} x) - η_{2} χ_{n}^{'} (m_{2} x)},

B_{n} = ψ_{n} (m_{2} x) \frac{η_{1} D_{n} (m_{2} x) - η_{2} D_{n} (m_{1} x)}{η_{1} χ_{n}^{'} (m_{2} x) - η_{2} D_{n} (m_{1} x) χ_{n} (m_{2} x)} .

(11)

χ_n (ρ)=-ρy_n (ρ) is used in Eqs. (11). If η ₁=η ₂=1, it is obvious that A_n =B_n , D̃_n =G̃_n , and a_n =b_n . According to Eqs. (8), it is the same as the homogeneous sphere since S ₁₁(θ=180°)=0. Furthermore, S ₁₂=S ₃₄=0 and S ₁₁=S ₃₃ also hold for coated spheres if η ₁=η ₂.

3. Numerical results

Fig. 1. The phase function as a function of the scattering angle for homogeneous spherical particles. η=1 and the size parameter x=10 are used. The refractive indices plotted are m=1.1,1.3,1.7,1.1+0.1i,1.3+0.1i,1.7+0.1i

Table 1. Extinction and scattering efficiencies for the cases in Fig. 1

	m=1.1	m=1.3	m=1.7	m=1.1+0.1i	m=1.3+0.1i	m=1.7+0.1i
Q_ext	1.7129	2.6575	1.7716	1.9358	2.4718	2.3655
Q_sca	1.7129	2.6575	1.7716	0.8926	1.2506	1.1172

Two codes were written on the basis of the theory and formulas in Sec. 2 to calculate the scattering properties of homogeneous and coated spheres whose relative permeabilities are not unity. The various numerical techniques suggested by Bohren and Huffman [16

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).

], Wiscombe [18

W. J. Wiscombe, “Mie scattering calculation,” NCAR Tech. Note TN-140+STR (National Center for Atmospheric Research, Boulder, Colo., 1979).

, 19

W.J. Wiscombe, “Improved Mie scattering algorithms,” Appl. Opt. 19, 1505–1509 (1980). [CrossRef] [PubMed]

], and Toon and Ackerman [20

O. B. Toon and T. P. Ackerman, “Algorithms for the calculation of scattering by starified spheres,” Appl. Opt. 20, 3657–3660 (1981). [CrossRef] [PubMed]

] are employed in this study. For homogeneous spheres with η=1, it would be interesting to show how fast S ₁₁ approaches zero around θ=180° for different values of m. Figure 1 plots the phase function, the normalized scattering matrix element S ₁₁, for m=1.1,1.3,1.7,1.1+0.1i,1.3+0.1i, and 1.7+0.1i. In this computation, η=1 and the size parameter x=10 are assumed. Evidently, the phase function approaches zero faster in the region of 170°<θ<180° as the relative refractive index decreases. For m=1.3 and 1.7, the rates of decrease of the phase functions are almost the same. The rate of decrease of the phase functions stops if the refractive index reaches a certain value. In addition, an imaginary part of the refractive index causes the phase function to decrease faster. Table 1 shows the extinction efficiencies Q_ext and scattering efficiencies Q_sca for the cases in Fig. 1. The cases with zero imaginary part of the refractive indices, the extinction and scattering efficiencies are equal. It is interesting that the extinction efficiency for m=1.3+0.1i is smaller than that of m=1.3, while the extinction efficiencies of m=1.1+0.1i and m=1.7+0.1i are larger than those of m=1.1 and m=1.7, respectively.

Fig. 2. The phase function as a function of the scattering angle for homogeneous spherical particles. The size parameter is x=10. The refractive indices is m=1.3. η=0.9999,0.999,0.99,0.91 are plotted.

Table 2. Extinction/scattering efficiencies for cases in Fig. 2

η=0.9999	η=0.999	η=0.99	η=0.91
2.6575	2.6575	2.6575	2.6594

To understand how the sensitivity of the phase function depends on η, Fig. 2 shows the phase functions for spherical particles with η=0.9999,0.999,0.99,0.91for scattering angles ranging from 160° to 180°. A size parameter of x=10 and a relative refractive index of m=1.3 are used. It is shown that the backscatter is small if η remains close to 1. As η deviates from 1, the backscatter also gradually deviates from 0. Table 2 shows the extinction/scattering efficiencies for cases in Fig. 2. The extinction and scattering efficiencies are equal for the four cases because the imaginary part of the refractive indices are zero. The extinction efficiencies do not have significant changes until η reaches 0.91.

For coated spheres, the phase function is also zero in the backscattering direction if η _1,2=1 is satisfied for each layer, which means that it can be completely hidden from monostatic lidar/radar detection. An important question is to find under what conditions a coated sphere can still hide from monostatic lidar/radar detection even if its core region has η ₁≠1. If these conditions are found, the shell of a coated sphere can be regarded as a cloak. From physical intuition, it is speculated that these conditions may include a thick outer shell with a large imaginary refractive index. Ref. [21

P. -W. Zhai, Y. -K. Lee, G. W. Kattawar, and P. Yang, “Implementing the Near- to Far-Field Transformation in the Finite-Difference Time-Domain Method,” Appl. Opt. 43, 3738–3746 (2004). [CrossRef] [PubMed]

] has shown the distribution of the internal electric field of a homogeneous sphere for different values of refractive index. It was shown that the electric field is mainly distributed around the boundary of the particle for a large refractive index of m=7.1499+2.914i. The physical reason is that the large absorption caused by the large imaginary part prevents the penetration of electromagnetic waves into the particle. It can be concluded that the electromagnetic waves cannot “see” the inner region of a particle provided it has a large imaginary refractive index. Applying this rule to a coated sphere, we have the following conditions under which a coated sphere can still hide from a monostatic lidar/radar even if the inner region has η ₁≠1:

η_{2} = 1, τ = k (b - a) Im (m_{2}) = (y - x) Im (m_{2}) ≫ 1,

(12)

where k is the wave vector in the surrounding medium; x=ka and y=kb are the size parameters of the inner and outer spheres, respectively; Im(m₂) is the imaginary part of the relative refractive index of the outer layer.

Fig. 3. The phase function as a function of the scattering angle for coated spherical particles. The size of the inner sphere is x=5. and m₁=1.5, η ₁=2/3. Two refractive indices for the outer sphere are chosen: m₂=2+0.5i and m₂=2+i. η ₂=1. For each of the refractive indices, a proper value of y is used to produce τ=1 and τ=4 according to τ=(y-x)Im(m₂).

Table 3. Extinction and scattering efficiencies for cases in Fig. 3. Also shown are the ratios of the extinction and scattering cross sections of the coated particles and those of the “Not coated” case.

	Not coated	Im(m₂)=0.5		Im(m₂)=1
	Not coated	τ=1	τ=4	τ=1	τ=4
Q_ext	3.9278	2.5314	2.3368	2.3572	2.4121
Q_sca	3.9278	1.1852	1.1507	0.9987	1.1544
C_ext/C^Nc _ext	-	1.2632	4.0218	0.8642	1.9897
C_sca/C^Nc _sca	-	0.5914	1.9804	0.3661	0.9522

Figure 3 shows the phase functions for four coated spheres. The refractive indices for the outer shell are taken as m₂=2+0.5i and m₂=2+i. η ₂=1 is used to insure the outer layer produces a zero backscatter. The size parameters for the inner sphere is x=5. A refractive index of m₁=1.5 is assumed for the inner core to test its sensitivity. The impedance for the inner sphere is η ₁=2/3. Thus the permeability of the inner core is taken to be 1. We noticed that τ=(y-x) Im(m₂) is the critical factor to achieve a zero backscatter. In Fig. 3, two values of τ=1 and 4 are chosen to test the sensitivity. For each value of τ and Im(m₂), the corresponding value of y can be determined by τ=(y-x)Im(m₂). The phase function for the inner sphere without the coated sphere is also shown for comparison. It is evident from Fig. 3 that the coated sphere decreases the phase function at the backscatter. The phase functions for τ=4 have negligible values while the phase functions for τ=1 still have significant values at the backscatter. This feature confirms that τ is the dominant factor that determines the magnitude of the backscatter. This is a remarkable result because τ=k(b-a)Im(m₂)=4 is just a fairly moderate value and the backscatter has already been reduced by several orders of magnitude. Table 3 shows the extinction and scattering efficiencies for cases in Fig. 3. The particles shown in Fig. 3 all have different sizes, which makes it hard to interpret the efficiencies. To overcome this problem, the ratio of the extinction/scattering cross sections for the coated particle and the “Not coated” particle is introduced.

\frac{C_{ext, sca}}{C_{ext, sca}^{Nc}} = \frac{Q_{ext, sca} \cdot y^{2}}{C_{ext, sca}^{Nc} \cdot x^{2}},

(13)

where C^Nc _ext,sca is the extinction/scattering cross section for the “Not coated” case; C_ext,sca is the extinction/scattering cross section for the coated particle. The proportionality between the size parameter and the radius of a spherical particle is used to get the right hand side of Eq. (13). The size parameter x of the “Not coated” case is equal to that of the cores of the coated particles. Equation (13) represents the ratio of the power eliminated/scattered from the incident light for the coated and “Not coated” particles. Table 3 shows the extinction cross sections for the coated particles generally larger than that of the “Not coated” particle. The only exception is in the case of τ=1 and Im(m₂)=0.5. It means that the negligible backscatters are achieved at the cost of larger extinction cross sections.

4. Conclusions

The expansion coefficients of the scattered field are given for light scattering by both homogeneous and coated spheres whose relative permeability values are not assumed to be 1. Two Fortran codes were written to calculate the scattering matrix for these particles. An interesting feature is that a zero backscatter for homogeneous and coated spheres can be achieved if the impedance for each region of the particle is equal to that of the surrounding medium. If the inner core of a coated sphere does not satisfy this condition, the back scattering can still be approximately zero under the conditions that the impedance η ₂=1 holds for its outer layer and τ=(y-x)Im(m₂) is sufficiently large. In other words, this coated sphere can act as an invisible cloak for monostatic lidar/radar detection. Furthermore, the extinction cross section for the coated sphere is generally larger than the “Not coated” particle.

Acknowledgments

The authors thanks Dr. Craig F. Bohren, professor in Department of Meteorology in the Pennsylvania State University, Dr. Warren J. Wiscombe, the Chief Scientist of the Atmospheric Radiation Measurement (ARM) Program, and the anonymous reviewers of this paper, for their insightful comments and suggestions on impedance matching. This research was supported by the Office of Naval Research under contracts N00014-02-1-0478 and N00014-06-1-0069. For

References and links

1.	J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling Electromagnetic Fields,” Science 312, 1780–1782 (2006). [CrossRef] [PubMed]
2.	U. Leonhardt, “Optical Conformal Mapping,” Science 312, 1777–1780 (2006). [CrossRef] [PubMed]
3.	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]
4.	A. J. Ward and J. B. Pendry, “Refraction and geometry in maxwells equations,” J. Mod. Opt. 43, 773–793 (1996). [CrossRef]
5.	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]
6.	S. A. Cummer, B.-I. Popa, D. Schurig, D. R. Smith, and J. B. Pendry, “Full-wave simulations of electromagnetic cloaking structures,” Phys. Rev. E 74, 036621-1–5 (2006).
7.	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]
8.	H. Chen, B. I. Wu, B. Zhang, and J. A. Kong, “Electromagnetic Wave Interactions with a Metamaterial Cloak,” Phys. Rev. Lett. 99, 063903-1–4 (2007). [CrossRef] [PubMed]
9.	B. Zhang, H. Chen, B. Wu, Y. Luo, L. Ran, and J. A. Kong, “Response of a cylindrical invisibility cloak to electromagnetic waves,” Phys. Rev. B 76, 121101-1–4 (R) (2007) [CrossRef]
10.	A. Alu and N. Engheta, “Achieving transparency with plasmonic and metamaterial coatings,” Phys. Rev. E 72, 016623-1–9 (2005). [CrossRef]
11.	G. W. Milton and N.-A. P. Nicorovici, “On the cloaking effects associated with anomalous localized resonance,” Proc. R. Soc. London, Ser. A 462, 3027–3059 (2006). [CrossRef]
12.	G. Mie, “Beigrade zur optik truber medien, speziell kolloidaler metallosungen,” Ann. Phys. (Leipzig) 25, 377–455 (1908). [CrossRef]
13.	J. Maxwell Garnett, “Colours in metal glasses and in metallic films,” Philos. Trans. R. Soc. London 203, 385 (1904). [CrossRef]
14.	J.R. Liu, M. Itoh, and K.-I Machida, “Frequency dispersion of complex permeability and permittivity on ironbased nanocomposites derived from rare earth-iron intermetallic compounds,” J. Alloys Compd. 408–412, 1396–1399 (2006). [CrossRef]
15.	H.C. van de Hulst, Light Scattering by Small Particles (Dover, New York, 1981).
16.	C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).
17.	J. D. Jackson, Classical Electrodynamics, 3rd Edition (Wiley-VCH, 1998),
18.	W. J. Wiscombe, “Mie scattering calculation,” NCAR Tech. Note TN-140+STR (National Center for Atmospheric Research, Boulder, Colo., 1979).
19.	W.J. Wiscombe, “Improved Mie scattering algorithms,” Appl. Opt. 19, 1505–1509 (1980). [CrossRef] [PubMed]
20.	O. B. Toon and T. P. Ackerman, “Algorithms for the calculation of scattering by starified spheres,” Appl. Opt. 20, 3657–3660 (1981). [CrossRef] [PubMed]
21.	P. -W. Zhai, Y. -K. Lee, G. W. Kattawar, and P. Yang, “Implementing the Near- to Far-Field Transformation in the Finite-Difference Time-Domain Method,” Appl. Opt. 43, 3738–3746 (2004). [CrossRef] [PubMed]

OCIS Codes
(290.5850) Scattering : Scattering, particles
(280.1350) Remote sensing and sensors : Backscattering

ToC Category:
Scattering

History
Original Manuscript: October 17, 2007
Revised Manuscript: January 15, 2008
Manuscript Accepted: January 17, 2008
Published: January 18, 2008

Citation
Peng-Wang Zhai, Yu You, George W. Kattawar, and Ping Yang, "Monostatic lidar/radar invisibility using coated spheres," Opt. Express 16, 1431-1439 (2008)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-3-1431

Sort: Year | Journal | Reset

References

J. B. Pendry, D. Schurig, and D. R. Smith, "Controlling Electromagnetic Fields," Science 312, 1780-1782 (2006). [CrossRef] [PubMed]
U. Leonhardt, "Optical Conformal Mapping," Science 312, 1777-1780 (2006). [CrossRef] [PubMed]
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]
A. J. Ward and J. B. Pendry,"Refraction and geometry in maxwells equations," J. Mod. Opt. 43, 773-793 (1996). [CrossRef]
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]
S. A. Cummer, B.-I. Popa, D. Schurig, D. R. Smith, and J. B. Pendry, "Full-wave simulations of electromagnetic cloaking structures," Phys. Rev. E 74, 036621-1-5 (2006).
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]
H. Chen, B. I. Wu, B. Zhang, and J. A. Kong, "Electromagnetic Wave Interactions with a Metamaterial Cloak," Phys. Rev. Lett. 99, 063903-1-4 (2007). [CrossRef] [PubMed]
B. Zhang, H. Chen, B. Wu, Y. Luo, L. Ran, and J. A. Kong, "Response of a cylindrical invisibility cloak to electromagnetic waves," Phys. Rev. B 76, 121101-1-4 (R) (2007) [CrossRef]
A. Alu and N. Engheta, "Achieving transparency with plasmonic and metamaterial coatings," Phys. Rev. E 72, 016623-1-9 (2005). [CrossRef]
G. W. Milton and N.-A. P. Nicorovici, "On the cloaking effects associated with anomalous localized resonance," Proc. R. Soc. London, Ser. A 462, 3027-3059 (2006). [CrossRef]
G. Mie, "Beigrade zur optik truber medien, speziell kolloidaler metallosungen," Ann. Phys. (Leipzig) 25, 377- 455 (1908). [CrossRef]
J. Maxwell Garnett, "Colours in metal glasses and in metallic films," Philos. Trans. R. Soc. London 203, 385 (1904). [CrossRef]
J.R. Liu, M. Itoh, and K.-I Machida, "Frequency dispersion of complex permeability and permittivity on ironbased nanocomposites derived from rare earth-iron intermetallic compounds," J. Alloys Compd. 408-412, 1396- 1399 (2006). [CrossRef]
H.C. van de Hulst, Light Scattering by Small Particles (Dover, New York, 1981).
C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).
J. D. Jackson, Classical Electrodynamics, 3rd Edition (Wiley-VCH, 1998),
W. J. Wiscombe, "Mie scattering calculation," NCAR Tech. Note TN-140+STR (National Center for Atmospheric Research, Boulder, Colo., 1979).
W.J. Wiscombe, "Improved Mie scattering algorithms," Appl. Opt. 19, 1505-1509 (1980). [CrossRef] [PubMed]
O. B. Toon and T. P. Ackerman, "Algorithms for the calculation of scattering by starified spheres," Appl. Opt. 20, 3657-3660 (1981). [CrossRef] [PubMed]
P. -W. Zhai, Y. -K. Lee, G. W. Kattawar, and P. Yang, "Implementing the Near- to Far-Field Transformation in the Finite-Difference Time-Domain Method," Appl. Opt. 43, 3738-3746 (2004) [CrossRef] [PubMed]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper