International Journal of Antennas and PropagationVolume 2012 (2012), Article ID 782320, 13 pagesdoi:10.1155/2012/782320

Review Article

Analytical Method and Semianalytical Method for Analysis of Scattering by Anisotropic Sphere: A Review

Chao Wan¹ and Hongning Li²

¹Department of Electrical and Computer Engineering, National University of Singapore, Kent Ridge, Singapore, 117583²State Key Laboratory of Integrated Service Networks, Xidian University, Xi'an, Shaanxi 710071, China

Received 14 July 2011; Accepted 3 November 2011

Academic Editor: Joshua Le-Wei Li

Copyright © 2012 Chao Wan and Hongning Li. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

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 Sphere

2.1. Cartesian Anisotropy of Spheres

For Cartesian anisotropic medium, the constitutive relations are given by where the identity dyadic is expressed as . 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 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 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 by The parameters are defined in Cartesian coordinates as (1). We can rewrite the Maxwell's equations in a convenient form [27]: The solution to (7) can be written in the Fourier transform: where the vector wave number is expressed as , the space vector as , and 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 can be seen as a combination of infinite plane wave. Then, for any direction wave, we have . Substituting (8) into (7), the electric field equation is transformed into wherewith In order to make (9) exist all the time, we write the eigenvalue equation as Here, the eigenvalue is a function of the angular , the equation turns to be where and the eigenvectors can also be obtained from (10) at the same time and are given as follows: where q = 1, 2, 3, or 4, here is a plane wave, and with Here, all the angular and eigenvalues should be considered. The E-field in (7) is then given as Only two roots need to be considered, while the periodic function denoting the angular spectrum amplitude can be expressed as where is the associated Legendre function, and means that is from 0 to and is from to . Substituting (20) to (18), we obtain Using the identity to stand for the part , Substituting (22) into (21), we achieve the solution of 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: The other internal field parameters, , , and can be obtained and details are from [27, 33].

Inserting (23) into (21), and integrating with respect to , we have 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 , the incident electromagnetic fields can be expanded by an infinite series of spherical vector wave functions for an isotropic medium as follows: where the scattering fields are expanded as where and are unknown coefficients and stand for the external field information, and , and 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 and we have where The scattering coefficients and , are thus expressed as 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 to rewrite the Maxwell equation [32]: The divergenceless property (31) suggests that be expanded in terms of the vector spherical wave functions and [32, 41, 42]: where is not undetermined. As mentioned above, there are three kinds of VSWF's , and . For this reason, (31) does not involve , 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 withHere, all the coefficients are obtained using the orthogonality of the VSWFs [32, 41, 42].

Therefore, one has where Repeat the procedure and reuse the properties of the VSWFs bellow until (31) can be written in the form of , , and : One gets with Then, we write the equation form: The expression for , , , is in [32]. Let and , denote, respectively, the eigenvalues and the corresponding eigenvectors of eigensystem (40), with representing the index of eigenvalues and corresponding eigenvectors. One can then build a new set of vector in space based on the eigenvectors: Thus, we can express as Use the relation of (6), and fields can be written as 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, For scattered wave, 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 then rewrite the matrix in a clear form: with Afterwards, we can solve the equations and get the coefficients matrix, and the scattering information is in the matrix : where

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 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 . As we can see, when the off-diagonal parameters increase, the backscattering () will also increase. Compared to Figure 1, Figure 2 shows the gyromagnetic influence on the RCS, here , 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 , . When increasing the permittivity elements of the materials, the scattering for larger than area will increase obviously.

Figure 1: Normalized Radar cross-section (RCS) values versus the scattering angle. Different gyroelectric influence is considered. Here, .

Figure 2: Normalized Radar cross-section (RCS) values versus the scattering angle. Different gyromagnetic influence is considered. Here, .

Figure 3: 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 .

3. Electromagnetic Wave Scattering by a Radial Anisotropic Sphere

3.1. Radial Anisotropy of Spheres

As another type medium, the radial anisotropy materials have parameters given by [15] where the identity dyadic is expressed as . 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] Then, the E-field and H-field can be expressed in the new coordinates: Similarly, the displacement vector is provided: Here, two operators are defined as then the Laplace operator can be written as and Maxwell equations is correspondingly rewritten as 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]: where The field can be expressed in the following form: Also when operate on other components, it has the following properties: we write the (58) in a convenient form: with 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 (), 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: 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, and are not coupled and we encounter the most discussed case [51]: For a plane wave, there are many discussions already; here, we only show the result for the internal and external field: 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 () influence on the backscattering when . Sharp changes can be observed when the radius changes. Compared to Figure 4, Figure 5 shows the positive anisotropy () 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.

Figure 4: Normalized backscattering values versus sphere radius. Different negative anisotropic value influence is considered. Here, we increase the negative values.

Figure 5: Normalized backscattering values versus sphere radius. Different positive anisotropic value influence is considered. Here, we reduce the positive values.

Figure 6: Normalized backscattering values versus sphere radius. Lossy negative anisotropic influence is considered. Here, the imaginary part for the permeability is bigger than 1.

Figure 7: 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 , and . 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.

References

L. Lorentz, “Lysbevaegelsen i og uden for en haf plane lysbolger belyst kulge,” Videnskab Selskab Skrifter, vol. 6, pp. 1–62, 1890.

G. Mie, “Beiträge zur optik trüber medien, speziell kolloidaler metallösungen,” Annalen der Physik, vol. 25, no. 3, pp. 377–445, 1908.

B. S. Luk'yanchuk and V. Ternovsky, “Light scattering by a thin wire with a surface-plasmon resonance: bifurcations of the Poynting vector field,” Physical Review B, vol. 73, no. 23, Article ID 235432, 2006.

B. S. Luk’yanchuk, M. I. Tribelsky, Z. B. Wang et al., “Extraordinary scattering diagram for nanoparticles near plasmon resonance frequencies,” Applied Physics A, vol. 89, no. 2, pp. 259–264, 2007.

M. I. Tribelsky, S. Flach, A. E. Miroshnichenko, A. V. Gorbach, and Y. S. Kivshar, “Light scattering by a finite obstacle and fano resonances,” Physical Review Letters, vol. 100, no. 4, Article ID 043903, 2008.

W. Fang, J. Y. Xu, A. Yamilov et al., “Large enhancement of spontaneous emission rates of InAs quantum dots in GaAs microdisks,” Optics Letters, vol. 27, no. 11, pp. 948–950, 2002.

C. Qiu, L. Gao, J. D. Joannopoulos, and M. Soljačić, “Light scattering from anisotropic particles: propagation, localization, and nonlinearity,” Laser and Photonics Reviews, vol. 4, no. 2, pp. 268–282, 2010.

S. Zhang, Y. Xiong, G. Bartal, X. Yin, and X. Zhang, “Magnetized plasma for reconfigurable subdiffraction imaging,” Physical Review Letters, vol. 106, no. 24, 2011.

J. Yao, X. Yang, X. Yin, G. Bartal, and X. Zhang, “Three-dimensional nanometer-scale optical cavities of indefinite medium,” Proceedings of the National Academy of Sciences of the United States of America, vol. 108, no. 28, pp. 11327–11331, 2011.

S. Liu, W. Lu, Z. Lin, and S. T. Chui, “Molding reflection from metamaterials based on magnetic surface plasmons,” Physical Review B, vol. 84, no. 4, 2011.

S. Liu, J. Du, Z. Lin, R. X. Wu, and S. T. Chui, “Formation of robust and completely tunable resonant photonic band gaps,” Physical Review B, vol. 78, no. 15, 2008.

S. Liu, W. Chen, J. Du, Z. Lin, S. T. Chui, and C. T. Chan, “Manipulating negative-refractive behavior with a magnetic field,” Physical Review Letters, vol. 101, no. 15, Article ID 157407, 2008.

Y. Poo, R.-X. Wu, Z. Lin, Y. Yang, and C. T. Chan, “Experimental realization of self-guiding unidirectional electromagnetic edge states,” Physical Review Letters, vol. 106, no. 9, 2011.

J. Roth and M. J. Dignam, “Scattering and extinction cross sections for a spherical particle coated with an oriented molecular layer,” The Journal of the Optical Society of America, vol. 63, no. 3, pp. 308–311, 1973.

J. C. Monzon, “Three-dimensional field expansion in the most general rotationally symmetric anisotropic material: application to scattering by a sphere,” IEEE Transactions on Antennas and Propagation, vol. 37, no. 6, pp. 728–735, 1989.

A. D. Kiselev, V. Y. Reshetnyak, and T. J. Sluckin, “Light scattering by optically anisotropic scatterers: T-matrix theory for radial and uniform anisotropies,” Physical Review E, vol. E65, no. 5, Article ID 056609, pp. 1–16, 2002.

M. F. Yanik, H. Cinar, H. N. Cinar, A. D. Chisholm, Y. Jin, and A. Ben-Yakar, “Functional regeneration after laser axotomy,” Nature, vol. 432, no. 7019, p. 822, 2004.

R. D. Weglein, “A model for predicting acoustic material signatures,” Applied Physics Letters, vol. 34, no. 3, pp. 179–181, 1979.

J. Kushibiki, A. Ohkubo, and N. Chubachi, “Anisotropy detection in sapphire by acoustic microscope using line-focus beam,” Electronics Letters, vol. 17, no. 15, pp. 534–536, 1981.

A. C. Guo, J. R. Petrella, J. Kurtzberg, and J. M. Provenzale, “Evaluation of white matter anisotropy in Krabbe disease with diffusion tensor MR imaging: initial experience,” Radiology, vol. 218, no. 3, pp. 809–815, 2001.

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science, vol. 312, no. 5781, pp. 1780–1782, 2006.

X. Q. Sheng and Z. Peng, “Analysis of scattering by large objects with off-diagonally anisotropic material using finite element-boundary integral-multilevel fast multipole algorithm,” IET Microwaves, Antennas and Propagation, vol. 4, pp. 492–500, 2010.

W. Ren, “Contributions to the electromagnetic wave theory of bounded homogeneous anisotropic media,” Physical Review E, vol. 47, no. 1, pp. 664–673, 1993.

Y. Geng, X. Wu, and L.-W. Li, “Analysis of electromagnetic scattering by a plasma anisotropic sphere,” Radio Science, vol. 38, no. 6, pp. 12/1–12/12, 2003.

Y. L. Geng and X. B. Wu, “A Plane Electromagnetic Wave Scattering By a Ferrite Sphere,” Journal of Electromagnetic Waves and Applications, vol. 18, no. 2, pp. 161–179, 2004.

Y.-L. Geng, X.-B. Wu, L.-W. Li, and B.-R. Guan, “Mie scattering by a uniaxial anisotropic sphere,” Physical Review E, vol. 70, no. 5, article 056609, pp. 1–8, 2004.

Y. L. Geng and C. W. Qiu, “Analytical spectral-domain scattering theory of a general gyrotropic sphere,” Materials Science. In press, http://arxiv.org/abs/1102.4057

Y.-L. Geng, “Scattering of a plane wave by an anisotropic ferrite-coated conducting sphere,” IET Microwaves, Antennas and Propagation, vol. 2, no. 2, pp. 158–162, 2008.

B. Stout, M. Nevière, and E. Popov, “Mie scattering by an anisotropic object. Part I. Homogeneous sphere,” Journal of the Optical Society of America A, vol. 23, no. 5, pp. 1111–1123, 2006.

B. Stout, M. Nevière, and E. Popov, “Mie scattering by an anisotropic object. Part II. Arbitrary-shaped object: differential theory,” Journal of the Optical Society of America A, vol. 23, no. 5, pp. 1124–1134, 2006.

J. Roth and M. J. Dignam, “Scattering and extinction cross sections for a spherical particle coated with an oriented molecular layer,” Journal of the Optical Society of America A, vol. 63, no. 3, pp. 308–311, 1973.

W. L. Ong, “New solutions to EM scattering by a gyrotropic sphere,” Bachelor Thesis, National University of Singapore, 2010

Y. L. Geng, “Application of Spherical Vector Wave Function to Electromagnetic Scattering by Anisotropic Media,” PhD. Thesis, Xi’an, China, Xidian University, 2006.

C. T. Tai, Dyadic Green’s Functions in Electromagnetic Theory, IEEE Press, New York, NY, USA, 1993.

N. K. Uzunoglu, P. G. Cottis, and J. G. Fikioris, “Excitation of electromagnetic waves in a gyroelectric cylinder,” IEEE Transactions on Antennas and Propagation, vol. 33, no. 1, pp. 90–99, 1985.

S. N. Papadakis, N. K. Uzunoglu, and N. Christos, “Scattering of a plane wave by a general anisotropic dielectric ellipsoid,” Journal of the Optical Society of America A, vol. 7, no. 6, pp. 991–997, 1990.

W. C. Chew and Y. M. Wang, “Fast algorithm for solution of a scattering problem using a recursive aggregate T-matrix method,” Microwave and Optical Technology Letters, vol. 3, no. 5, pp. 164–169, 1990.

D. Sarkar and N. J. Halas, “General vector basis function solution of Maxwell's equations,” Physical Review E, vol. 56, no. 1, pp. 1102–1112, 1997.

D. Sarkar, Vector basis function solution of Maxwell's equations, Ph.D. thesis, Rice University, Houston, Tex, USA, 1996.

S. Liu, “Electromagnetic properties for three dimensional anisotropic photonic crystal,” PhD. Thesis, Fudan University, Shanghai, China, 2007.

Z. Lin and S. T. Chui, “Electromagnetic scattering by optically anisotropic magnetic particle,” Physical Review E, vol. 69, no. 5, article 056614, pp. 1–14, 2004.

S. Liu and Z. Lin, “Opening up complete photonic bandgaps in three-dimensional photonic crystals consisting of biaxial dielectric spheres,” Physical Review E, vol. 73, no. 6, 2006.

L. Gao, T. H. Fung, K. W. Yu, and C. W. Qiu, “Electromagnetic transparency by coated spheres with radial anisotropy,” Physical Review E, vol. 78, no. 4, article 046609, 2008.

C.-W. Qiu, A. Novitsky, H. Ma, and S. Qu, “Electromagnetic interaction of arbitrary radial-dependent anisotropic spheres and improved invisibility for nonlinear-transformation-based cloaks,” Physical Review E, vol. 80, no. 1, 2009.

C. W. Qiu, L. Hu, B. Zhang, B. I. Wu, S. G. Johnson, and J. D. Joannopoulos, “Spherical cloaking using nonlinear transformations for improved segmentation into concentric isotropic coatings,” Optics Express, vol. 17, no. 16, pp. 13467–13478, 2009.

C.-W. Qiu, L. Hu, X. Xu, and Y. Feng, “Spherical cloaking with homogeneous isotropic multilayered structures,” Physical Review E, vol. 79, no. 4, article 047602, 2009.

A. Novitsky, C. W. Qiu, and S. Zouhdi, “Transformation-based spherical cloaks designed by an implicit transformation-independent method: theory and optimization,” New Journal of Physics, vol. 11, Article ID 113001, 2009.

B. S. Luk'yanchuk and C.-W. Qiu, “Enhanced scattering efficiencies in spherical particles with weakly dissipating anisotropic materials,” Applied Physics A, vol. 92, no. 4, pp. 773–776, 2008.

Y. X. Ni, D. L. Gao, Z. F. Sang, L. Gao, and C. W. Qiu, “Influence of spherical anisotropy on the optical properties of plasmon resonant metallic nanoparticles,” Applied Physics A, vol. 102, no. 3, pp. 673–679, 2011.

W.-J. Chen, J. C. W. Lee, J.-W. Dong, C.-W. Qiu, and H.-Z. Wang, “Fano resonance of three-dimensional spiral photonic crystals: paradoxical transmission and polarization gap,” Applied Physics Letters, vol. 98, no. 8, article 081116, 2011.

C.-W. Qiu, L.-W. Li, T.-S. Yeo, and S. Zouhdi, “Scattering by rotationally symmetric anisotropic spheres: potential formulation and parametric studies,” Physical Review E, vol. 75, no. 2, 2007.

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles, John Wiley & Sons, New York, NY, USA, 1998.

C.-W. Qiu, L.-W. Li, Q. Wu, and T.-S. Yeo, “Field representations in general gyrotropic media in spherical coordinates,” IEEE Antennas and Wireless Propagation Letters, vol. 4, no. 1, pp. 467–470, 2005.

A. Novitsky and L. Barkovsky, “Matrix approach for light scattering from a multilayered rotationally symmetric bianisotropic sphere,” Physical Review A, vol. 77, no. 3, Article ID 033849, 2008.

G. Arfken and H. J. Weber, Mathematical Method for Physicists, Elsevier Academic Press, 6th edition, 2005.