The history of methods for the electromagnetic scattering by an anisotropic sphere has been reviewed. Two main methods, angular expansion method and T-matrix method, which are widely used for the anisotropic sphere, are expressed in Cartesian coordinate firstly. The comparison of those and the further exploration on the scattering field are illustrated afterwards. Based on the most general form concluded by variable separation method, the coupled electric field and magnetic field of radial anisotropic sphere can be derived. By simplifying the condition, simpler case of uniaxial anisotropic media is expressed with confirmed coefficients for the internal and external field. Details of significant phenomenon are presented.
1. Introduction
Since the properties of isotropic homogeneous sphere is formulated by Lorentz [1] and Mie [2], due to plenty of corresponding significant phenomena [3–6] and increasing applications in atmospheric optics, remote sensing, computational electromagnetic and photonics, electromagnetic wave propagation and scattering for a sphere has been a subject of great interest in recent years. The amazing properties of solid anisotropic materials in nature and classes of artificial metamaterials of anisotropy would possibly produce plentiful technological applications. In the meantime, fabrication technology development has greatly enlarged the area of the anisotropic materials.
Cartesian anisotropic medium and Radial anisotropic medium are considered as two classic kinds of anisotropic materials [7]. Some novel properties of Cartesian anisotropy medium have been studied: magnetic plasma realizing the subdiffraction imaging [8], optical cavity made by indefinite medium of anisotropy [9], the sharply asymmetric reflection (SAR) effect [10], formation of tunable resonant photonic band gaps [11], negative-refractive behavior manipulation [12], and self-guiding unidirectional electromagnetics edge [13]. Based on the general anisotropic materials, those amazing properties can greatly enlarge the possibility to control light and wave in special structure and also open a new research area in the interaction between the electromagnetic wave and materials. For radial anisotropic materials, the spherical particle is of great interest [14–16]. The radial anisotropic particle scattering provides an insight into the interaction between the microparticles and the microwave or optical illumination by external devices [17–19], which could help to detect the position for some abnormal proteins [7, 20]. The classic spherical cloaks have been achieved by radial anisotropic materials [21]. The surface plasma resonance invisibility and extraordinary scattering of the RA particle provide a good opportunity to study the physical phenomena of microparticles.
For most of these published works, many properties of anisotropic material have been explored and applied. However, little work has been done so far on summaries and analysis for the methods to deal with problems of such anisotropic material. Scientists and engineers have to check the anisotropic class first and then choose a convenient numerical method or some commercial software. Either way will encounter redundant difficulties: for numerical method like FDTD, FEM, MOM, and other improved ones [22], prediction of the properties of medium is quite a difficulty, let alone the new structure design; for commercial software, the methods embedded in them are commonly standard, which are not able to deal with the possible novel materials with abnormal parameters. Therefore, modeling of those novel materials becomes extremely difficult. Fortunately, the analytical method and semianalytical method such as angular expansion method and T-matrix method will be very helpful in concluding good predictions for the anisotropic problems while keeping the accuracy [23–33].
It is of major purpose to summarize the analytical and semianalytical methods and their material parameter's effects on the scattering patterns, which will be of great use in electromagnetic wave properties for Cartesian and the Radial anisotropic spheres. In Section 2, besides two general analytical and semianalytical algorithms of anisotropy in Cartesian are reviewed, advantages and disadvantages of these two methods are also illustrated. In Section 3, the variables separation method for a single radial anisotropic sphere is demonstrated and numerical results obtained. A short conclusion is provided in Section 4.
2. Electromagnetic Wave Scattering by a Cartesian Anisotropic Sphere2.1. Cartesian Anisotropy of Spheres
For Cartesian anisotropic medium, the constitutive relations are given byε̅=εs(εr-iεκ0iεκεr0001)μ̅=μs(μr-iμκ0iμκμr0001),
where the identity dyadic is expressed as Î=x̂x̂+ŷŷ+ẑẑ. By extending Mie theory in isotropic cases, many methods such as FDTD, FEM, and MOM have been established to analyzed wave interaction issue in the medium. However, they cannot provide an explicit explanation for the physical phenomenon. Being able to deal with multilayered problems, the dyadic green's function is considered as one analytical approach [34]. However, the researchers have to encounter the singularities while they choose the integration path, then the expression turns to be rather complicated. In [35], Uzunoglu first employed the Fourier form to state the field, followed by a lot of work on this direction for its clear computation procedure [36, 37]. Using the angular expansion to present electromagnetic field, Ren successfully transferred the three-dimensional problem into a two-dimension integration [23]. The field can be treated as combination of infinite plane wave, while only four kinds of Bessel function need to be considered according to the wave direction. After that, Dr. Sarkar gives an explicit expression of the plane wave using vector spherical wave functions, where the field in the media is a form of that infinite plane wave [38, 39]. Based on these clear field and plane wave expression, Dr. Geng proposed a transpicuous procedure and made a fundamental computation [24–28]. The same result has been achieved by other researcher when compared with Geng's work [29]. Besides, the method firstly dealing with the photonic band problem can be seen as another way [40]. Instead of giving an explicit expression and using the integral form on the surface, Lin used the transfer matrix method to gather the information in the medium [41]. By using the complete matrix form, they express the wave number k by a matrix vector instead of an expression [32, 42], and this approach obtains great success [10–13]. Herein, we mainly focus on the review and analysis of these two analytical and semianalytical methods.
2.2. Angular Expansion Method
Inside the sourceless and homogeneous anisotropic sphere, the constitutive relationships as used in the Maxwell equations are∇×E=iωB,∇×H=-iωD,∇⋅D=0,∇⋅B=0.
The relations between the electric displacement vector D, the magnetic induction B, the electric field E, and the magnetic field H inside the medium are given byD=ε̅⋅E,B=μ̅⋅H.
The parameters are defined in Cartesian coordinates as (1). We can rewrite the Maxwell's equations in a convenient form [27]:∇×[μ̅-1⋅∇×E]-ω2ε̅⋅E=0.
The solution to (7) can be written in the Fourier transform:E(r)=∫-∞+∞dkx∫-∞+∞dky∫-∞+∞E(k)eik⋅rdkz,
where the vector wave number is expressed as k=kxx̂+kyŷ+kzẑ, the space vector as r=xx̂+yŷ+zẑ, and x̂,ŷ,ẑ being the unit vectors in Cartesian coordinates. Here, the problem is transferred from spatial domain to spectrum domain. For each point in the medium, the electric field E can be seen as a combination of infinite plane wave. Then, for any direction wave, we have ∇=ik. Substituting (8) into (7), the electric field equation is transformed into∫-∞+∞dkx∫-∞+∞dky∫-∞+∞A̅(k)⋅E(k)eik⋅rdkz=0,
whereA̅(k)=(-b1kz2-b3kz2+a1b3kxky-ib2kz2-ia2b1kxkz+ib2kykzb3kxky-ib2kz2-ia2-b1kz2-b3kz2+a1b1kykz-ib2kxkzb1kxkz-ib2kykzb1kykz+ib2kxkz-b1(kx2+ky2)+a3)witha1=ω2εsεra2=ω2εsεκa3=ω2εs,b1=μrμs(μr2-μκ2)b2=μκμs(μr2-μκ2)b3=1μs.
In order to make (9) exist all the time, we write the eigenvalue equation asDet[A̅(k)]=0.
Here, the eigenvalue k is a function of the angular (θk,ϕk), the equation turns to beA(θk,ϕk)k4-B(θk,ϕk)k2+C=0,
whereA(θk,ϕk)=[b1b3sin2θk+(b12-b22)cos2θk]×[a1sin2θk+a3cos2θk],B(θk,ϕk)=[b1(a12-a22)+b3a1a3]sin2θk+2a3(b1a1+b2a2)cos2θk,C=a3(a12-a22),k2=kx2+ky2+kz2,θk=tan-1(kx2+ky2kz),ϕk=tan-1(kykx),k1,32=B+B2-4AC2A,k2,42=B-B2-4AC2A,
and the eigenvectors can also be obtained from (10) at the same time and are given as follows:Eq=Fqefq(θk,ϕk)=[Fqxe(θk,ϕk)x̂+Fqye(θk,ϕk)ŷ+Fqze(θk,ϕk)ẑ]fq(θk,ϕk),
where q = 1, 2, 3, or 4, here Eq is a plane wave, andFqxe=-Δ1Δsinϕk+Δ2Δcosϕk,Fqye=Δ1Δcosϕk+Δ2Δsinϕk,Fqze=1
withΔ1=i(b1a2+b2a1)kq2sinθkcosθk,Δ2=[b1b3kq2sin2θk+(b12-b22)kq2cos2θk]×kq2sinθkcosθk-(b1a1+b2a2)Δ=-(b2kq2cos2θk+a2)2+(b1kq2cos2θk-a1)×(b1kq2cos2θk+b3kq2sin2θk-a1).
Here, all the angular and eigenvalues should be considered. The E-field in (7) is then given asE(r)=∑q=12∫0π∫02πFqe(θk,ϕk)fq(θk,ϕk)eikq⋅rkq2sinθkdθkdϕk,kq=kqsinθkcosϕkx̂+kqsinθksinϕkŷ+kqcosθkẑ.
Only two roots need to be considered, while the periodic function fq(θk,ϕk) denoting the angular spectrum amplitude can be expressed asfq(θk,ϕk)=∑m′,n′Gm′n′qPn′m′(cosθk)eim′ϕk,
where Pnm(x) is the associated Legendre function, and ∑m′,n′ means that n′ is from 0 to +∞ and m′ is from -n′ to n′. Substituting (20) to (18), we obtainE(r)=∑q=12∑m′,n′Gm′n′q∫0π∫02πFqe(θk,ϕk)Pn′m′×(cosθk)eim′ϕkeikq⋅rkq2sinθkdθkdϕk.
Using the identity to stand for the part eikq·r,eik⋅r=∑n=0∞in(2n+1)jn(kr)×[∑m=0n(n-m)!(n+m)!Pnm(cosθk)Pnm(cosθ)eim(ϕ-ϕk)+∑m=0n(n-m)!(n+m)!Pnm(cosθk)Pnm(cosθ)e-im(ϕ-ϕk)].
Substituting (22) into (21), we achieve the solution of E(r) for a general homogeneous gyrotropic anisotropic media. Then, the vector spherical wave functions is employed to express the eigenvector, which is one key step. After this, the internal field can be in the form of VSWFs:Fqe(θ,ϕ)eikq⋅r=∑m,n[Amnqe(θk)Mmn(1)(r,kq)+Bmnqe(θk)Nmn(1)(r,kq)+Cmnqe(θk)Lmn(1)(r,kq)]e-imϕk.
The other internal field parameters, Amnqe(θk), Bmnqe(θk), and Cmnqe(θk) can be obtained and details are from [27, 33].
Inserting (23) into (21), and integrating with respect to ϕk, we haveE(r)=∑q=12∑m,n∑n′2πGmn′q×∫0π[Amnqe(θk)Mmn(1)(r,kq)+Bmnqe(θk)Nmn(1)(r,kq)+Cmnqe(θk)Lmn(1)(r,kq)]×Pn′m′(cosθk)eim′ϕkkq2sinθkdθkdϕk,
which is an eigenfunction representation of the E-field in gyrotropic anisotropic media. We can get the H-field in a similar form [27, 33].
From the result in (24), it is found that we can use the VSWFs to present the field, and this procedure is also suitable for the external area. Afterwards, we take the plane wave as an example. Assuming that an incident plane wave is E=x̂E0eik0z, the incident electromagnetic fields can be expanded by an infinite series of spherical vector wave functions for an isotropic medium as follows:Einc=E0∑n,m[δm,1+δm,-1]×[amnxMmn(1)(k0,r)+bmnxNmn(1)(k0,r)],Hinc=k0iωμ0E0∑n,m[δm,1+δm,-1]×[amnxMmn(1)(k0,r)+bmnxNmn(1)(k0,r)],
whereamnx={in+12n+12n(n+1),m=1,in+12n+12,m=-1,bmnx={in+12n+12n(n+1),m=1,-in+12n+12,m=-1,δs,l={1,s=l,0,s≠l,
the scattering fields are expanded asEs=∑n,m[AmnsMmn(3)(k0,r)+BmnxNmn(3)(k0,r)],Hs=k0iωμ0∑n,m[AmnsMmn(3)(k0,r)+BmnxNmn(3)(k0,r)],
where Amns and Bmns are unknown coefficients and stand for the external field information, and k0=ωμ0ε0, μ0 and ε0 denote the free space, wave number, permeability, and permittivity. Until here, this method is being discussed in an unbounded material. We should match the boundary condition on the surface when we limit it to an anisotropic sphere. The tangential components of the electromagnetic field continues at r=a and we have∑q=12∑n′∞2πGmn′q∫0πQmnqPn′m′(cosθk)kq2sinθkdθk=E0[δm,1+δm,-1]amnx⋅i(k0a)2,∑q=12∑n′∞2πGmn′q∫0πRmnqPn′m′(cosθk)kq2sinθkdθk=E0[δm,1+δm,-1]bmnx⋅i(k0a)2,
whereQmnq={Amnqe1k0rddr[rhn(1)(k0r)]jn(kqr)-iωμ0k0[Bmnqh1kqrddr[rjn(kqr)]+Cmnqhjn(kqr)r]⋅hn(1)(k0r)Amnqe1k0r}r=a,Rmnq={iωμ0k0Amnqh1k0rddr[rhn(1)(k0r)]jn(kqr)-[Bmnqe1kqrddr[rjn(kqr)]+Cmnqejn(kqr)r]⋅hn(1)(k0r)iωμ0k0Amnqh}r=a.
The scattering coefficients Amns and Bmns, are thus expressed asAmns=1hn(1)(k0a)[∑n′=0∞∑q=122πGmnq∫0πAmnqejn(kqa)×Pn′′m′(cosθk)kq2sinθkdθk-E0[δm,1+δm,1]amnxjn(k0a)∫],Bmns=1hn(1)(k0a)[iωμ0k0∑n′=0∞∑q=122πGmnq∫0πAmnqhjn(kqa)×Pn′m′(cosθk)kq2sinθkdθk-E0[δm,1+δm,1]bmnxjn(k0a)iωμ0k0∑n′=0∞∑q=122πGmnq∫0πAmnqhjn(kqa)iωμ0k0].
Then the internal and external field can be expressed in the form of the coefficients.
2.3. T-Matrix Method
Different from the angular expansion method, here we use magnetic induction B to rewrite the Maxwell equation [32]:∇×[εsε̅-1⋅(∇×μsμ̅-1⋅B)]-ks2B=0.
The divergenceless property (31) suggests that B be expanded in terms of the vector spherical wave functions Mmn(1)(k,r) and Nmn(1)(k,r) [32, 41, 42]:B=∑n,mE̅mn[dmnMmn(1)(k,r)+cmnNmn(1)(k,r)],
where k is not undetermined. As mentioned above, there are three kinds of VSWF's Mmn(1)(k,r), Nmn(1)(k,r) and Lmn(1)(k,r). For this reason, (31) does not involve Lmn(1)(k,r), because its curl is 0, for details of the parameters see [32, 41, 42]. Since the VSWFs is a complete frame and can stand for all vectors in the space, we can use them to present the elements in (31), written asεsε̅-1⋅Mmn=∑v=0+∞∑u=-vv[õuvmnMuv+p̃uvmnNuv+q̃uvmnLuv],εsε̅-1⋅Nmn=∑v=0+∞∑u=-vv[o̅uvmnMuv+p̅uvmnNuv+q̅uvmnLuv],μsμ̅-1⋅Mmn=∑v=0+∞∑u=-vv[g̃uvmnMuv+ẽuvmnNuv+f̃uvmnLuv],μsμ̅-1⋅Nmn=∑v=0+∞∑u=-vv[g̅uvmnMuv+e̅uvmnNuv+f̅uvmnLuv],
withõuvmn=δnvδmu+[(n2+n-m2)ε̅r′+mεκ′]δnvδmun(n+1),p̃uvmn=i(n+m)[mε̅r′-(n+1)εκ′]δn-1,vδmun(2n+1)+i(n-m+1)[mε̅r′+nεκ′]δn+11,vδmu(n+1)(2n+1),q̃uvmn=-i(n+m)[mε̅r′-(n+1)εκ′]δn-1,vδmu(2n+1)+i(n-m+1)[mε̅r′+nεκ′]δn+11,vδmu(2n+1),o̅uvmn=-i(n+m)(n+1)[mε̅r′+(n-1)εκ′]δn-1,vδmun(n-1)(2n+1)-i(n-m+1)n[mε̅r′-(n+2)εκ′]δn+11,vδmu(n+1)(n+2)(2n+1),p̅uvmn=δnvδmu+{[(2n2+2n+3)m2+(2n2+2n-3)n(n+1)]ε̅r′+(4n2+4n-3)mεκ′}δnvδmun(n+1)(2n-1)(2n+3)-(n+1)(n+m-1)(n+m)ε̅r′δn-2,vδmu(n-1)(2n-1)(2n+1)-n(n-m+1)(n-m+2)ε̅r′δn+2,vδmu(n+2)(2n+1)(2n+3),q̅uvmn=-[(n2+n-3m2)ε̅r′-m(2n-1)(2n+3)εκ′]δnvδmu(2n-1)(2n+3)+(n+1)(n+m-1)(n+m)ε̅r′δn-2,vδmu(2n-1)(2n+1)-n(n-m+1)(n-m+2)ε̅r′δn+2,vδmu(2n+1)(2n+3),g̃uvmn=δnvδmu+[(n2+n-m2)μ̅r′+mμκ′]δnvδmun(n+1),ẽuvmn=i(n+m)[mμ̅r′-(n+1)μκ′]δn-1,vδmun(2n+1)+i(n-m+1)[mμ̅r′+nμκ′]δn+11,vδmu(n+1)(2n+1),f̃uvmn=-i(n+m)[mμ̅r′-(n+1)μκ′]δn-1,vδmu(2n+1)+i(n-m+1)[mμ̅r′+nμκ′]δn+11,vδmu(2n+1),g̅uvmn=-i(n+m)(n+1)[mμ̅r′+(n-1)μκ′]δn-1,vδmun(n-1)(2n+1)-i(n-m+1)n[mμ̅r′-(n+2)μκ′]δn+11,vδmu(n+1)(n+2)(2n+1),e̅uvmn=δnvδmu+{[(2n2+2n+3)m2+(2n2+2n-3)n(n+1)]μ̅r′+(4n2+4n-3)mμκ′}δnvδmun(n+1)(2n-1)(2n+3)-(n+1)(n+m-1)(n+m)μ̅r′δn-2,vδmu(n-1)(2n-1)(2n+1)-n(n-m+1)(n-m+2)μ̅r′δn+2,vδmu(n+2)(2n+1)(2n+3),f̅uvmn=-[(n2+n-3m2)μ̅r′-m(2n-1)(2n+3)μκ′]δnvδmu(2n-1)(2n+3)+(n+1)(n+m-1)(n+m)μ̅r′δn-2,vδmu(2n-1)(2n+1)-n(n-m+1)(n-m+2)μ̅r′δn+2,vδmu(2n+1)(2n+3).Here, all the coefficients are obtained using the orthogonality of the VSWFs [32, 41, 42].
Therefore, one hasμsμ̅-1⋅B=∑n,mE̅mn[d̅mnMmn(1)(k,r)+c̅mnNmn(1)(k,r)+w̅mnLmn(1)(k,r)]+w̅00L00(1)(k,r),
where
d̅mn=∑v,uE̅uvE̅mn[g̃mnuvduv+g̅mnuvcuv],c̅mn=∑v,uE̅uvE̅mn[ẽmnuvduv+e̅mnuvcuv],w̅mn=∑v,uE̅uvE̅mn[f̃mnuvduv+f̅mnuvcuv],w̅00=-[23μκ′d01+215μ̅r′c02]E0.
Repeat the procedure and reuse the properties of the VSWFs bellow until (31) can be written in the form of Mmn(1)(k,r), Nmn(1)(k,r), and Lmn(1)(k,r):∇×Mmn(j)-kNmn(j)=0,∇×Nmn(j)-k2Mmn(j)=0,∇×Lmn(j)=0.
One gets∑n,mE̅mn[d≈mnMmn(1)(k,r)+c≈mnNmn(1)(k,r)]=0
withd≈mn=k2∑v,uE̅uvE̅mn[p̅mnuvd̅uv+p̃mnuvc̅uv]-ks2dmn,c≈mn=k2∑v,uE̅uvE̅mn[o̅mnuvd̅uv+õmnuvc̅uv]-ks2cmn.
Then, we write the equation form:(P̃P̅ÕO̅)(dc)=λ(dc).
The expression for P̃, P̅, Õ, O̅ is in [32]. Let λl and (dmn,l,cmn,l)T, denote, respectively, the eigenvalues and the corresponding eigenvectors of eigensystem (40), with l representing the index of eigenvalues and corresponding eigenvectors. One can then build a new set of vector Vl in space based on the eigenvectors:Vl=-klω∑n,mE̅mn[dmn,lMmn(1)(kl,r)+cmn,lNmn(1)(kl,r)].
Thus, we can express B asB=∑lαlVl=-∑lαlklω∑n,mE̅mn[dmn,lMmn(1)(kl,r)+cmn,lNmn(1)(kl,r)].
Use the relation of (6), H and E fields can be written asH=μ̅-1⋅B=-∑n,mE̅mn∑lklωμsαl×[d̅mn,lMmn(1)(kl,r)+c̅mn,lNmn(1)(kl,r)+wmn,lLmn(1)(kl,r)]-∑lklωμsαl[w00,lL00(1)(k,r)],E=-i∑n,mE̅mn∑lαl[wmn,lλlcmn,lMmn(1)(kl,r)+dmn,lNmn(1)(kl,r)+wmn,lλlLmn(1)(kl,r)]-i∑lαl[w̅̅00,lλlL00(1)(k,r)].
As we have discussed in the angular expansion method, here, we also give the form of the incident wave and scatted wave.
For incident wave,Einc=-∑n,miE̅mn[pmnNmn(1)(k,r)+qmnMmn(1)(k,r)],Hinc=-k0ωμ0∑n,mE̅mn[qmnNmn(1)(k,r)+pmnMmn(1)(k,r)].
For scattered wave,Es=∑n,miE̅mn[amnNmn(3)(k0,r)+bmnMmn(3)(k0,r)],Hs=k0ωμ0∑n,mE̅mn[bmnNmn(3)(k0,r)+amnMmn(3)(k0,r)].
The relevant coefficients are in [41]. Then, match boundary conditions on the sphere surface, we can get a system of equations for the coefficients, which is[ξn′(x)ψn′(x)]amn+∑l[1msk̅lψn′(k̅lmsx)ψn′(x)dmn,l]αl+∑l[μ0μsjn(k̅lmsx)ψn′(x)ωmn,l]αl=pmn,[ξn(x)ψn(x)]bmn+∑l[1msk̅lψn(k̅lmsx)ψn(x)cmn,l]αl=qmn,[ξn(x)ψn(x)]amn+∑l[μ0λlμsψn(k̅lmsx)ψn(x)dmn,l]αl=pmn,[ξn′(x)ψn′(x)]bmn+∑l[μ0λlμsψn′(k̅lmsx)ψn′(x)cmn,l]αl+∑l[μ0μsjn(k̅lmsx)ψn′(x)ωmn,l]αl=qmn,
then rewrite the matrix in a clear form:(Λ00Λ̅)(ab)+(VV̅)α̃=(pq),(Λ̅00Λ)(ab)+(U̅U)α̃=(pq)
withΛmn,uv=Sn(x)δnvδmu,Λ̅mn,uv=S̅n(x)δnvδmu,Umn,l=1msk̅lT̅n(x,msk̅lx)cmn,l,U̅mn,l=1msk̅lT̅n(x,msk̅lx)dmn,l,Vmn,l=μ0λlμsTn(x,msk̅lx)dmn,l,V̅mn,l=μ0λlμsT̅n(x,msk̅lx)cmn,l+Wmn,l,Wmn,l=μ0μs1msk̅lxT̅n(x,msk̅lx)Dn(1)(msk̅lx)wmn,l,Sn(x)=ξn(x)ψn(x)S̅n(x)=ξn′(x)ψn′(x),Tn(x,z)=ψn(z)ψn(x)T̅n(x,z)=ψn′(z)ψn′(x).
Afterwards, we can solve the equations and get the coefficients matrix, and the scattering information is in the matrix (ab):α̃=R(pq),(ab)=S(pq),
whereZ=(Y00-Y)-1(V-U̅V̅-U),S=ZR,Y=Λ̅-Λ.
2.4. Numerical Results and Comparison
For the two methods mentioned above, difference between them will be analyzed as follows. In angular expansion method, the plane wave expansion is employed and the wave in such complex medium is a combination of plane wave of different direction. The whole procedure is clear and easy to understand. At the same time, the wave number k has its physical meaning, which can be expressed. Theoretically, this method can deal with any complex medium. Different from angular expansion method, the T-matrix method does not give an expression for the wave number, the properties of spherical vector wave function are used instead and all the wave information is gathered from the matrix. In computation, the T-matrix method has more advantage over angular expansion method, which does not need to consider grid and can be applied to calculate larger problem.
As many discussions and comparisons published already, here we mainly show the strength of the two methods [32]. Figure 1 shows the gyroelectric influence on the Radar cross-section, here μκ=0,μsμr=4μ0μs=2μ0,εsεr=4ε0,εs=2ε0. As we can see, when the off-diagonal parameters increase, the backscattering (180∘) will also increase. Compared to Figure 1, Figure 2 shows the gyromagnetic influence on the RCS, here εsεr=2.4ε0,εs=2ε0, the less scattering areas is reduced from three to one. Figure 3 gives a general case, the influences of permittivity and permeability are both considered; and this is also a lossy sphere, where εs=(2+1.0i)ε0, μs=(2+1.0i)μ0. When increasing the permittivity elements of the materials, the scattering for larger than 80∘ area will increase obviously.
Normalized Radar cross-section (RCS) values versus the scattering angle. Different gyroelectric influence is considered. Here, k0r=π.
Normalized Radar cross-section (RCS) values versus the scattering angle. Different gyromagnetic influence is considered. Here, k0r=0.75π.
Normalized Radar Cross Section (RCS) values versus the scattering angle. Here the permeability and the permittivity are complex number. The Loss gyrotropic influence is considered. Here k0r=0.75π.
3. Electromagnetic Wave Scattering by a Radial Anisotropic Sphere3.1. Radial Anisotropy of Spheres
As another type medium, the radial anisotropy materials have parameters given by [15]ε̅=(εrr000τσ0-στ),μ=(μrr000γξ0-ξγ),
where the identity dyadic is expressed as Î=r̂r̂+θ̂θ̂+ϕ̂ϕ̂. Different from the Cartesian anisotropic (CA), less discussion have been published, because the radial anisotropic (RA) material is difficult to be fabricated [7]. The research attracts much attention recently due to the properties of materials with novel parameters. How RA particles interact with waves is an essential topic which will provide us a more physical insight into the invisibility phenomenon [43–47], enhanced surface plasmon resonance [48, 49], and Fano resonance [50]. Although Monzon developed the variables separation method and demonstrated a perfect explanation to the solution form [15], researchers prefer the Debye potential or other potentials similar [51–53]. Recently, Novitsky has stated a general form for the material, and then extend the T-matrix method to a multilayer case [54]. When a Hamilton operator is used, the Maxwell equations will be greatly simplified. By employing the impedance matrix, the field area is both simple and clear. However, this solution is just another form of variables separation method and only suitable to the uniaxial anisotropic case, which does not satisfy the general situation [54]. Herein, we explain Monzon's method and conclude an explanation to the solution.
3.2. Variables Separation Method
Firstly, we define new vectors to be [15]ξ̂±=θ̂±jϕ̂2sinθ.
Then, the E-field and H-field can be expressed in the new coordinates:E(r)=Err̂+e-ξ̂++e+ξ̅-e±=sinθ2(Eθ±jEϕ),H(r)=Hrr̂+h-ξ̂++h+ξ̅-h±=sinθ2(Hθ±jHϕ),
Similarly, the displacement vector is provided:ε̅⋅E=εrrErr̂++α+e-ξ̅++α-e+ξ̅-,α±=τ±jσ,μ̅⋅H=μrrHrr̂++β+h-ξ̅++β-h+ξ̅-,β±=γ±jζ.
Here, two operators are defined asD±=12r(sinθ∂∂θ±j∂∂ϕ),Dr=1r∂∂r⋅r,
then the Laplace operator can be written as∇=r∂∂r+ξ̂+D-+ξ̂-D+,
and Maxwell equations is correspondingly rewritten as∇×E=jr̂sin2θ{D+e--D-e+}+jξ̅+{D-Er-Dre-}+jξ̅-{Dre+-D+Er},∇×H=jr̂sin2θ{D+h--D-h+}+jξ̅+{D-Hr-Drh-}+jξ̅-{Drh+-D+Hr},
we insert (54) and (57) into (2) and (3). After some calculations including repeated operations on Maxwell’s equations, a coupled set of differential equations involving only the radial field components can be obtained [15]:{μrrγrDr2(rHr)+2sin2θD+D-Hr+ω2τμrr(τ2+σ2)Hr}+jωεrrγτr(σγ+τζ)Dr(rEr)=0,{εrrτrDr2(rEr)+2sin2θD+D-Er+ω2τεrr(γ2+ζ2)Er}-jωμrrγτr(σγ+τζ)Dr(rHr)=0,
whereD+D-=12r2[∂2∂ϕ2+sinθ∂∂θsinθ∂∂θ].
The field can be expressed in the following form:Er=∑n,mEmn(r)Pnm(cosθ)ejmϕ,Hr=∑n,mHmn(r)Pnm(cosθ)ejmϕ.
Also when D+D- operate on other components, it has the following properties:D+D-=-n(n+1)sin2θ2r2,
we write the (58) in a convenient form:d2f+α1(ρ)f-δdg=0d2g+α2(ρ)g-δdf=0
withd=ddρρ=κrκ=ωεrrμrrδ=(σγ+τζ)τγεrrμrr,α1(ρ)=A1-B1ρ2α2(ρ)=A2-B2ρ2,A1=τγ(γ2+ζ2)εrrμrr,A2=γτ(τ2+σ2)εrrμrr,B1=τεrrn(n+1),B2=γμrrn(n+1).
For (62), there is an analytical solution to the differential equations only when the coupled component is 0 (the isotropic or uniaxial anisotropic case). For a general case (δ≠0), there is no analytical solution for this kind of differential equations in mathematics [55]. Therefore, we have to use the series to stand for the equation roots and the procedure is as follows:f=∑n=-∞∞ρλ+nKn,g=∑n=-∞∞ρλ+nQn[(λ+n)(λ+n-1)-B1]Kn+A1Kn-2=δ(λ+n-1)Qn-1,[(λ+n)(λ+n-1)-B2]Qn+A2Qn-2=δ(λ+n-1)Kn-1,Q1=δ(τ1+1/2)[2τ1+1+(B1-B2)],K2=δ(τ1+3/2)Q1-A14(τ1+1).
When the original condition is given, the coefficients can be calculated step by step and then an expression for the radial component as we needed from (66) and (67) can be achieved. If the off-diagonal parameters are zero, Er and Hr are not coupled and we encounter the most discussed case [51]:Er=ωikt2(∂2∂r2+kt2)ΨTM,Eθ=ωikt21r∂2ΨTM∂r∂θ-1ε0εt1rsinθ∂ΨTE∂ϕ,Eϕ=ωikt21rsinθ∂2ΨTM∂r∂ϕ+1ε0εt1r∂ΨTE∂θ,Hr=ωikt2(∂2∂r2+kt2)ΨTE,Hθ=1μ0μt1rsinθ∂ΨTM∂ϕ+ωikt21r∂2ΨTE∂r∂θ,Hϕ=-1μ0μt1r∂ΨTM∂θ+ωikt21rsinθ∂2ΨTE∂r∂ϕ.
For a plane wave, there are many discussions already; here, we only show the result for the internal and external field:an=μt/εtjn(k0a)jv1′(kta)-jn′(k0a)jv1(kta)hn(2)′(k0a)jv1(kta)-μt/εthn(2)(k0a)jv1′(kta)Tn,bn=μt/εtjn′(k0a)jv2(kta)-jn(k0a)jv2′(kta)hn(2)(k0a)jv2′(kta)-μt/εthn(2)′(k0a)jv2(kta)Tn,cn=iμt/εthn(2)(k0a)jv1′(kta)-hn(2)′(k0a)jv1(kta)Tn,dn=iμt/εthn(2)(k0a)jv2′(kta)-μt/εthn(2)′(k0a)jv2(kta)Tn,Tn=i-n(2n+1)n(n+1).
Using these coefficients, there are many interesting results [7], which will not be discussed here.
3.3. Numerical Results
Results of the method using the above coefficients are presented in the following. Figure 4 states the negative anisotropy (μrr<γ) influence on the backscattering when εrr=τ=ε0. Sharp changes can be observed when the radius changes. Compared to Figure 4, Figure 5 shows the positive anisotropy (μrr>γ) influence on the backscattering, from which the oscillations are more regular. Figures 6 and 7 give a general case for the absorbing sphere. Taking the influence of negative absorbing into accounting, we can see from Figure 6 that the backscattering will approach to a certain value as the radius increases. There will be more backscattering as the sphere has much loss. The positive absorbing properties are discussed in Figure 7, which behaves a similar phenomenon with Figure 6.
Normalized backscattering values versus sphere radius. Different negative anisotropic value influence is considered. Here, we increase the negative values.
Normalized backscattering values versus sphere radius. Different positive anisotropic value influence is considered. Here, we reduce the positive values.
Normalized backscattering values versus sphere radius. Lossy negative anisotropic influence is considered. Here, the imaginary part for the permeability is bigger than 1.
Normalized backscattering values versus sphere radius. Lossy positive anisotropic influence is considered. Here, the imaginary part for the permeability is smaller than 1.
4. Discussion and Conclusion
In this review, two different types of anisotropy have been discussed in Cartesian and spherical coordinates, respectively. The role of anisotropy in scattering properties is characterized analytically. Two main methods are introduced in Cartesian coordinates. As the first one, the angular expansion method can be used in the infinitely large area and present an obvious explanation for the properties of the medium. Different from the first method, the T-matrix method can avoid the complex function by using the orthogonal properties for Mmn(l),Nmn(l), and Lmn(l). There is no need for researchers to consider the mesh before the calculation. Since these two methods are analytical and semianalytical, respectively, the parameter’s influence on every part is very clear and thus gives us an insight into anisotropy materials.
For a spherical coordinate anisotropic material, no matter which operators are used, the variables separation remains the main method to settle the problem. Though the Debye potentials can be widely used based on the properties of the Maxwell equations, it is still not suitable to solve the coupled gyrotropic materials. Monzon’s work demonstrated the reason why and how the electric field and magnetic field coupled. When a specific situation is considered, for example, a uniaxial anisotropic media, the electric field and the magnetic field can be separated and be both expressed in the radial function form. This review provides the convenience for researchers who want to explore new properties of anisotropic materials and, therefore, design new structures.
Acknowledgments
The authors wish to thank Professor Li Lewei, and Ms. Ong Weeling for helpful discussions. This work is supported by CSC scholarship and Key Program of NSFC-Guangdong Union Foundation under Grant no. U0835004.
LorentzL.Lysbevaegelsen i og uden for en haf plane lysbolger belyst kulge18906162MieG.Beiträge zur optik trüber medien, speziell kolloidaler metallösungen1908253377445Luk'yanchukB. S.TernovskyV.Light scattering by a thin wire with a surface-plasmon resonance: bifurcations of the Poynting vector field200673232-s2.0-3374565327610.1103/PhysRevB.73.235432235432Luk’yanchukB. S.TribelskyM. I.WangZ. B.ZhouY.HongM. H.ShiL. P.ChongT. C.Extraordinary scattering diagram for nanoparticles near plasmon resonance frequencies20078922592642-s2.0-3454825926310.1007/s00339-007-4099-1TribelskyM. I.FlachS.MiroshnichenkoA. E.GorbachA. V.KivsharY. S.Light scattering by a finite obstacle and fano resonances200810042-s2.0-3884915292710.1103/PhysRevLett.100.043903043903FangW.XuJ. Y.YamilovA.CaoH.MaY.HoS. T.SolomonG. S.Large enhancement of spontaneous emission rates of InAs quantum dots in GaAs microdisks200227119489502-s2.0-0036602053QiuC.GaoL.JoannopoulosJ. D.SoljačićM.Light scattering from anisotropic particles: propagation, localization, and nonlinearity2010422682822-s2.0-7704909559910.1002/lpor.200810078ZhangS.XiongY.BartalG.YinX.ZhangX.xiang@berkeley.eduMagnetized plasma for reconfigurable subdiffraction imaging20111062410.1103/PhysRevLett.106.243901YaoJ.YangX.YinX.BartalG.ZhangX.xiang@berkeley.eduThree-dimensional nanometer-scale optical cavities of indefinite medium201110828113271133110.1073/pnas.1104418108LiuS.LuW.LinZ.ChuiS. T.Molding reflection from metamaterials based on magnetic surface plasmons201184410.1103/PhysRevB.84.045425LiuS.DuJ.LinZ.WuR. X.ChuiS. T.Formation of robust and completely tunable resonant photonic band gaps2008781510.1103/PhysRevB.78.155101LiuS.ChenW.DuJ.LinZ.ChuiS. T.ChanC. T.Manipulating negative-refractive behavior with a magnetic field2008101152-s2.0-5484943630410.1103/PhysRevLett.101.157407157407PooY.WuR.-X.rxwu@nju.edu.cnLinZ.YangY.ChanC. T.Experimental realization of self-guiding unidirectional electromagnetic edge states2011106910.1103/PhysRevLett.106.093903RothJ.DignamM. J.Scattering and extinction cross sections for a spherical particle coated with an oriented molecular layer19736333083112-s2.0-0015593148MonzonJ. C.Three-dimensional field expansion in the most general rotationally symmetric anisotropic material: application to scattering by a sphere19893767287352-s2.0-002468161310.1109/8.29359KiselevA. D.ReshetnyakV. Y.SluckinT. J.Light scattering by optically anisotropic scatterers: T-matrix theory for radial and uniform anisotropies2002E6551162-s2.0-4134910427410.1103/PhysRevE.65.056609056609YanikM. F.CinarH.CinarH. N.ChisholmA. D.JinY.Ben-YakarA.Functional regeneration after laser axotomy200443270198222-s2.0-1654436319710.1038/432822aWegleinR. D.A model for predicting acoustic material signatures19793431791812-s2.0-000058691410.1063/1.90741KushibikiJ.OhkuboA.ChubachiN.Anisotropy detection in sapphire by acoustic microscope using line-focus beam198117155345362-s2.0-0019596041GuoA. C.PetrellaJ. R.KurtzbergJ.ProvenzaleJ. M.Evaluation of white matter anisotropy in Krabbe disease with diffusion tensor MR imaging: initial experience200121838098152-s2.0-0035111103PendryJ. B.SchurigD.SmithD. R.Controlling electromagnetic fields20063125781178017822-s2.0-3374548713810.1126/science.1125907ShengX. Q.PengZ.Analysis of scattering by large objects with off-diagonally anisotropic material using finite element-boundary integral-multilevel fast multipole algorithm20104492500RenW.Contributions to the electromagnetic wave theory of bounded homogeneous anisotropic media199347166467310.1103/PhysRevE.47.664GengY.gengyoulin@yahoo.com.cnWuX.wuxb@50.sh.cnLiL.-W.LWLi@nus.edu.sgAnalysis of electromagnetic scattering by a plasma anisotropic sphere200338612/112/12GengY. L.WuX. B.A Plane Electromagnetic Wave Scattering By a Ferrite Sphere200418216117910.1163/156939304323062022GengY.-L.WuX.-B.LiL.-W.GuanB.-R.Mie scattering by a uniaxial anisotropic sphere2004705, article 0566091810.1103/PhysRevE.70.056609GengY. L.QiuC. W.Analytical spectral-domain scattering theory of a general gyrotropic sphereMaterials Science. In press, http://arxiv.org/abs/1102.4057GengY.-L.gengyoulin@yahoo.com.cnScattering of a plane wave by an anisotropic ferrite-coated conducting sphere20082215816210.1049/iet-map:20070075StoutB.NevièreM.PopovE.Mie scattering by an anisotropic object. Part I. Homogeneous sphere2006235111111232-s2.0-3374523556510.1364/JOSAA.23.001111StoutB.NevièreM.PopovE.Mie scattering by an anisotropic object. Part II. Arbitrary-shaped object: differential theory2006235112411342-s2.0-3374524189310.1364/JOSAA.23.001124RothJ.DignamM. J.Scattering and extinction cross sections for a spherical particle coated with an oriented molecular layer19736333083112-s2.0-0015593148OngW. L.New solutions to EM scattering by a gyrotropic sphereBachelor Thesis, National University of Singapore, 2010GengY. L.Application of Spherical Vector Wave Function to Electromagnetic Scattering by Anisotropic Media2006Xidian UniversityXi’an, ChinaTaiC. T.1993New York, NY, USAIEEE PressUzunogluN. K.CottisP. G.FikiorisJ. G.Excitation of electromagnetic waves in a gyroelectric cylinder198533190992-s2.0-0021786698PapadakisS. N.UzunogluN. K.ChristosN.Scattering of a plane wave by a general anisotropic dielectric ellipsoid19907699199710.1364/JOSAA.7.000991ChewW. C.WangY. M.Fast algorithm for solution of a scattering problem using a recursive aggregate T-matrix method199035164169SarkarD.HalasN. J.General vector basis function solution of Maxwell's equations1997561110211122-s2.0-0000246382SarkarD.1996Houston, Tex, USARice UniversityLiuS.Electromagnetic properties for three dimensional anisotropic photonic crystal2007Shanghai, ChinaFudan UniversityLinZ.ChuiS. T.Electromagnetic scattering by optically anisotropic magnetic particle2004695, article 056614114LiuS.041019018@fudan.edu.cnLinZ.phlin@fudan.ac.cnOpening up complete photonic bandgaps in three-dimensional photonic crystals consisting of biaxial dielectric spheres200673610.1103/PhysRevE.73.066609GaoL.leigao@suda.edu.cnFungT. H.YuK. W.kwyu@phy.cuhk.edu.hkQiuC. W.chengwei.qiu@nus.edu.sgElectromagnetic transparency by coated spheres with radial anisotropy2008784, article 04660910.1103/PhysRevE.78.046609QiuC.-W.cwq@mit.eduNovitskyA.MaH.cwq@mit.eduQuS.Electromagnetic interaction of arbitrary radial-dependent anisotropic spheres and improved invisibility for nonlinear-transformation-based cloaks200980110.1103/PhysRevE.80.016604QiuC. W.HuL.ZhangB.WuB. I.JohnsonS. G.JoannopoulosJ. D.Spherical cloaking using nonlinear transformations for improved segmentation into concentric isotropic coatings2009171613467134782-s2.0-6834915444710.1364/OE.17.013467QiuC.-W.eleqc@nus.edu.sgHuL.XuX.FengY.Spherical cloaking with homogeneous isotropic multilayered structures2009794, article 04760210.1103/PhysRevE.79.047602NovitskyA.QiuC. W.ZouhdiS.Transformation-based spherical cloaks designed by an implicit transformation-independent method: theory and optimization2009112-s2.0-7204909014210.1088/1367-2630/11/11/113001113001Luk'yanchukB. S.Boris_L@dsi.a-star.edu.sgQiuC.-W.eleqc@nus.edu.sgEnhanced scattering efficiencies in spherical particles with weakly dissipating anisotropic materials200892477377610.1007/s00339-008-4572-5NiY. X.GaoD. L.SangZ. F.GaoL.leigao@suda.edu.cnQiuC. W.eleqc@nus.edu.sgInfluence of spherical anisotropy on the optical properties of plasmon resonant metallic nanoparticles2011102367367910.1007/s00339-010-6090-5ChenW.-J.LeeJ. C. W.DongJ.-W.dongjwen@mail.sysu.edu.cnQiuC.-W.WangH.-Z.stswhz@mail.sysu.edu.cnFano resonance of three-dimensional spiral photonic crystals: paradoxical transmission and polarization gap2011988, article 08111610.1063/1.3560338QiuC.-W.LiL.-W.lwli@nus.edu.sgYeoT.-S.ZouhdiS.Scattering by rotationally symmetric anisotropic spheres: potential formulation and parametric studies200775210.1103/PhysRevE.75.026609BohrenC. F.HuffmanD. R.1998New York, NY, USAJohn Wiley & SonsQiuC.-W.LiL.-W.elelilw@nus.edu.sgWuQ.qwu@hit.edu.cnYeoT.-S.eleyeots@nus.edu.sgField representations in general gyrotropic media in spherical coordinates20054146747010.1109/LAWP.2005.862055NovitskyA.BarkovskyL.Matrix approach for light scattering from a multilayered rotationally symmetric bianisotropic sphere20087732-s2.0-4154916337210.1103/PhysRevA.77.033849033849ArfkenG.WeberH. J.20056thElsevier Academic Press