# Mie Scattering by a Conducting Sphere Coated Uniaxial Single-Negative Medium

Authors: You-Lin Geng
Journal: International Journal of Antennas and Propagation
 (2012)

---

## Abstract

We propose an accurate analytical method to compute the electromagnetic scattering from three-dimensional (3D) conducting sphere coated uniaxial anisotropic single-negative (SNG) medium. Based on the spherical vector wave functions (SVWFs) in uniaxial anisotropic medium, the electromagnetic field in homogeneous uniaxial SNG medium and free space can be expressed by the SVWFs in uniaxial SNG medium and free space. The continued boundary conditions of electromagnetic fields between the uniaxial SNG medium and free space are applied, and the tangential electrical field is vanished in the surface of conducting sphere, the coefficients of scattering fields in free space can be derived, and then the character of scattering of conducting sphere coated homogeneous uniaxial SNG medium can be obtained. Some numericals are given in the end.


---

## Body

## 1. Introduction

Recently, the development of metamaterials, which was first introduced by Veselago in 1968 [1], has attracted more and more attention in the physical and engineering communities. Metamaterials, which are a kind of new artificial media, have been studied intensively since the concept of perfect lenses was proposed by Pendry in 2000 [2] and the negative index of refraction was experimentally verified by Smith et al. in 2001 [3]. Another breakthrough in metamaterials happened in 2006, when a method for controlling electromagnetic (EM) fields was proposed using inhomogeneous and anisotropic metamaterials [4] and the reduced invisible cloak was successfully realized using metamaterials at microwave frequencies [5]. Since then, more and more researchers have investigated the properties and applications of anisotropic metamaterials.The EM scattering problem for penetrable objects is very important to study material properties and arbitrary dielectric lenses. There are many classical methods for analysing traditional two-dimensional (2D) and three-dimensional (3D) dielectric scattering problems [6–11]. In recent years, there have also been many theoretical investigations on the study of EM characteristics for double-negative (DNG) media [12–15], where both the permittivity and permeability are negative.However, “single-negative” (SNG) materials in which only one of the material parameters, not both, has a negative real value may also possess interesting properties when they are paired in a conjugate manner. These media include the epsilon-negative (ENG) media, in which the real part of permittivity is negative, but the real permeability is positive, and the mu-negative (MNG) media, in which the real part of permeability is negative, but the real permittivity is positive. For instance, the idea of constructing an effective metamaterials by having layers of SNG media has been explored by Fredkin and Ron in [16], and some works have been done by Alù and Engheta [17–19]. It is well known that the negative medium is not isotropic; it is anisotropy; some works have been done between acoustic or electromagnetic waves and anisotropic metamaterials [20–22], and this paper discusses the scattering by a conducting sphere coated uniaxial anisotropic SNG medium.Based on the spherical vector wave functions in uniaxial sphere [11] and uniaxial sphere of left-handed material [13], electromagnetic fields in homogeneous uniaxial anisotropic single-negative medium can be expressed as the addition of the first and second SVWFs in uniaxial SNG medium. Applying the continued boundary conditions of electromagnetic fields, and the tangential electric field vanishing in the interface of the conducting sphere, the scattering coefficients of scattering fields in free space can be derived. The numerical results between the present method in this paper and Mie theory of scattering by a conducting sphere coated isotropic SNG medium are given, a good agreement is obtained as we expected, and some new numerical results are given in this paper.In the subsequent analysis, anexp(-iωt) time dependence is assumed for electromagnetic field quantities but suppressed throughout the treatment.


## 2. Formulas

Let us consider a conducting sphere coated homogeneous uniaxial SNG medium illuminated by an incident plane wave. As illuminated in Figure1, the coated sphere with outer radius a1 and inner radius a2 is located at the coordinate origin. On the surface of the inner conducting sphere, the uniaxial SNG medium with permittivity tensor (ϵ¯) and permeability tensor (μ¯) is coated with thickness d(=a1-a2). It is assumed that the incident wave propagates in the +z^ direction, the incident electric field has unity of amplitude, and is polarized in the +x^ direction.Figure 1
Geometry of a plane wave scattering by a plasma anisotropic spherical shell.The electric field vector wave equation in such a source-free uniaxial anisotropic SNG medium can be written in the following form [6, 11]:
(1)∇×[μ--1·∇×E(r)]-ω2ϵ-·E(r)=0,
where E denotes the electric field, while ϵ- and μ¯ represent the permittivity tensor and the permeability tensor of uniaxial SNG medium; the expression are [11, 13]
(2)ϵ-=[ϵt000ϵt000ϵz],μ=[μt000μt000μz].Some works were done about this negative medium in the past years [21, 22], and in this paper, either of the permittivity or permeability is negative; that is, the permittivity ϵ- is negative while the permeability μ¯ is positive, or on the contrary, the permeability μ¯ is negative while the permittivity ϵ- is positive. Using Fourier transform [6, 11], the expansion of plane wave factors in terms of spherical vector wave functions in isotropic medium [23], and the properties of spherical Bessel functions [24], the electromagnetic fields in the uniaxial anisotropic SNG medium can be obtained as follows:
(3)E=∑l=12∑q=12∑mnn′Fmn′q(l)∫0π[Amnqe(θk)Mmn(l)(r,kq)+Bmnqe(θk)Nmn(l)(r,kq)+Cmnqe(θk)Lmn(l)(r,kq)]×Pn′m(cosθk)kq2sinθkdθk,(4)H=∑l=12∑q=12∑mnn′Fmn′q(l)∫0π[Amnqh(θk)Mmn(l)(r,kq)+Bmnqh(θk)Nmn(l)(r,kq)+Cmnqh(θk)Lmn(l)(r,kq)]×Pn′m(cosθk)kq2sinθkdθk,
where n′ and n are summed up both from 0 to +∞ while m is summed up from -n to n, and r is pointing to the (θ,ϕ)-direction in the spherical coordinates. The coefficients, Fmnq(l), are unknown, as in [11, 13]. Amnqp(θk), Bmnqp(θk), Cmnqp(θk) (where p=e or h), and kq are functions of θk and they have been derived in [11, 13]. The vector wave functions, Mmn(l), Nmn(l), and Lmn(l) are spherical vector wave functions and they are also shown in [11, 23].The incident electromagnetic fields (designated by the superscriptinc) can be expanded in an infinite series in isotropic spherical vector wave functions [11, 23]:
(5)Einc=∑mn[amnxMmn(1)(r,k0)+bmnxNmn(1)(r,k0)]×[δm,1+δm,-1],(6)Hinc=k0iωμ0∑mn[amnxNmn(1)(r,k0)+bmnxMmn(1)(r,k0)]×[δm,1+δm,-1],
where the expansion coefficients are defined as:
(7)amnx={in+12n+12n(n+1),m=1,in+12n+12,m=-1;(8)bmnx={in+12n+12n(n+1),m=1,-in+12n+12,m=-1;(9)δs,l={1,s=l,0,s≠l.According to the radiation condition of an outgoing wave (attenuating to zero at infinity) and the asymptotic behavior of spherical Bessel functions, onlyhn(1) should be retained in the radial functions; therefore, the expansion of scattered fields (designated by the superscript s) are
(10)Es=∑mn[AmnsMmn(3)(r,k0)+BmnsNmn(3)(r,k0)],(11)Hs=k0iωμ0∑mn[AmnsNmn(3)(r,k0)+BmnsMmn(3)(r,k0)],
where the coefficients, Amns and Bmns (n varies from 0 to +∞ while m changes from -n to n), are unknowns to be determined, Mmn(l)(r,k0) and Nmn(l)(r,k0) denote the spherical vector wave functions defined in [11, 23], and k0=ω(ϵ0μ0)1/2 identifies the wave number of free space, respectively.Applying the boundary conditions at the surface of uniaxial anisotropic SNG medium, for example, whenr=a2, the expansion coefficients of electromagnetic fields in uniaxial anisotropic medium can be obtained by the following equations:
(12)∑l=12∑q=12∑n′=0∞Fmn'q(l)∫0πAmnqezn(l)(kqa2)Fmn'q(l)Fmn'q(l)×Pn′m(cosθk)kq2sinθkdθk=0,(13)∑l=12∑q=12∑n′=0∞Fmn′q(l)∫0π{zn(l)(kqr)rBmnqe1kqrddr(rzn(l)(kqr))+Cmnqezn(l)(kqr)r}r=a2∑l=12∑q=12∑n′=0∞Fm(l)×Pn′m(cosθk)kq2sinθkdθk=0
and r=a1 and can be obtained the following expression:
(14)∑l=12∑q=12∑n′=0∞Fmn′q(l)∫0πQmnq(l)Pn′m(cosθk)kq2sinθkdθk=[δm,1+δm,-1]amnxi(k0a1)2,(15)∑l=12∑q=12∑n′=0∞Fmn′q(l)∫0πRmnq(l)Pn′m(cosθk)kq2sinθkdθk=[δm,1+δm,-1]bmnxi(k0a1)2,
where expansion coefficients amnx and bmnx can be expressed in (7) and (8), and zn(l) (where l=1, 2, 3, and 4) denotes an appropriate kind of spherical Bessel functions, jn, yn, hn(1), and hn(2), respectively. Qmnq(l) and Rmnq(l) have the following expression:
(16)Qmnq(l)={zn(l)(kqr)rAmnqe1k0rddr(rhn(1)(k0r))zn(l)(kqr)-iωμ0k0[zn(l)(kqr)rBmnqh1kqrddr(rzn(l)(kqr))+Cmnqhzn(l)(kqr)r]·hn(1)(k0r)}r=a1,(17)Rmnq(l)={zn(l)(kqr)r1kqriωμ0k0Amnqh1k0rddr(rhn(1)(k0r))zn(l)(kqr)-[Bmnqe1kqrddr(rzn(l)(kqr))+Cmnqezn(l)(kqr)r]·hn(1)(k0r)iωμ0k0zn(l)(kqr)r}r=a1.The integral in (12) to (15) can be easily calculated by Gauss integral [24], and from (12) to (15), it shows that(i)
firstly, the unknown coefficients (Fmn'q(l)) of electromagnetic fields in the uniaxial anisotropic SNG medium can be obtained;(ii)
secondly, the coefficients (Amns, Bmns) of scattered fields in region 0 (free space) are calculated; and(iii)
lastly, the far scattering field of electromagnetic fields from a conductor sphere coated uniaxial anisotropic SNG medium by a plane wave and the radar cross-section are thus obtained.


## 3. Numerical Results and Discussion

In the last section, we have presented the necessary theoretical formulation of the electromagnetic fields of a plane wave scattered by a conducting sphere coated uniaxial anisotropic SNG medium. To gain more physics insight into the problem, we will provide in this section some numerical solutions to the problem of electromagnetic scattering by a conducting sphere coated uniaxial anisotropic SNG medium.Numerical computations have been performed by applying the theoretical formulae derived earlier in the previous sections. Since there are no numerical results of scattering by a conducting sphere coated uniaxial anisotropic SNG medium, in order to check the accuracy of the newly obtained numerical results, we performed one trial; that is, we calculated the radar cross-sections using the present method and the Mie theory (scattering by a conducting sphere coated homogeneous isotropic SNG medium). The results are shown in Figure2, where electric dimensions of outer and inner spherical surfaces are k0a1=1.1π and k0a2=π, while the permittivity and permeability tensor elements are ϵt=ϵz=2ϵ0, μt=μz=-2μ0, respectively, (where and subsequently, ϵ0 and μ0 stand for the free space permittivity and permeability, resp.).Figure 2
Radar cross sections (RCSs) versus scattering angleθ (in degrees): results of this paper (solid curve) and of Mie theory (black round).From Figure2, it is seen apparently that the radar cross sections calculated by using the two methods (i.e., the present method in this paper and Mie theory) are in very good agreement in both the E- and H-planes, where the maximum number of n' used in (12) to (15) is only 6 to achieve the convergence. It partially verifies the correctness and applicability of our theory as well as the Fortran program codes.After this, we obtain some new results unavailable elsewhere in the literature. Three examples are considered herein, and their radar cross sections are plotted in Figures3, 4, and 5.Figure 3
Radar cross-sections (RCSs) versus scattering angleθ (in degrees) in the E-plane (solid curve) and in the H-plane (short-dashed curve).Figure 4
Radar cross sections (RCSs) versus scattering angleθ (in degrees) in the E-plane (solid curve) and in the H-plane (short-dashed curve).Figure 5
Radar cross sections (RCSs) versus the thickness of the uniaxial SNG medium in the forward (solid curve) and in the back (dashed curve) directions, when inner radius is fixed ata2=0.75λ. The parameters of the uniaxial SNG coating of the conducting sphere are assumed to be ϵt=-1.5ϵ0, ϵt=-3ϵ0, μt=3μ0, and μz=1.5μ0, respectively.Figure3 represents radar cross sections of a conducting sphere coated uniaxial anisotropic SNG medium of more general uniaxial medium, where the permittivity and permeability tensor elements are characterized by ϵt=-2ϵ0, ϵz=-3ϵ0, and μt=3μ0, μz=2μ0, and the electric size of the uniaxial anisotropic SNG spherical shell is chosen as k0a1=2.1π and k0a2=2π. The maximum number n' in (12) to (15) to achieve a good convergence is found to be 16.Figure4 gives a numerical result of scattering by a conducting sphere coated loss uniaxial SNG medium. The electric size of the uniaxial anisotropic SNG spherical shell is chosen as k0a1=3π and k0a2=2π under the illumination by an incident plane wave. The permittivity and permeability tensor parameters used for this case are ϵt=(1+0.001i)ϵ0, ϵz=(2+0.001i)ϵ0, and μt=(-2+0.001i)μ0, μz=(-1+0.001i)μ0. As the electric dimension of the sphere is increased, the maximum number of n' used in (12) to (15) must be significantly increased to 22 to achieve the convergence.Later in this section, monostatic RCSs are given; for example, the forward RCS (the scattering angleθ=0) and the back RCS (the scattering angle θ=π) are given in Figure 5. The parameters of the coated uniaxial SNG medium are ϵt=-1.5ϵ0, ϵz=-3ϵ0, and μt=3μ0,μz=1.5μ0, and the electric size of conducting sphere is k0a2=1.5π, and the change of coated uniaxial SNG medium is from 0.1λ to 1.25λ (λ is the wavelength of incident wave in free space) at the interval of 0.025. From this Figure, the forward RCS of the core-shell system almost increases as the thickness of the uniaxial SNG coating increases. However, it is interesting to note that the back scattering will be oscillating. The RCS in the forward direction is always larger than that in the back direction.


## 4. Conclusions

The spherical vector wave function expansion solution to the plane wave scattering by a conducting sphere coated uniaxial anisotropic SNG medium is obtained analytically in this paper. The solution has only one-dimensional integral which can be calculated by Gauss integral [24] easily. Numerical results are obtained using the present method and compared with Mie theory and a fairly good agreement is observed. It is shown that the obtained solution is stable even for almost isotropic scatterers, since the proposed solution is an analytical one of a conducting sphere coated uniaxial SNG medium, and the result of the Mie theory is a special case of the present method. The general numerical results, including the conducting sphere coated lossy uniaxial SNG medium and resonance region, are given and are found reducible to those of spacial cases. The present analysis is believed to be useful in antenna and wave propagation.



---

*Source: 856476-2012-12-12.xml*