Research Article | Open Access

Lei Cao, Yongpin Chen, Kai Kang, "Anisotropic Scattering Characteristics of a Radially Multilayered Gyrotropic Sphere", *International Journal of Antennas and Propagation*, vol. 2017, Article ID 3205198, 8 pages, 2017. https://doi.org/10.1155/2017/3205198

# Anisotropic Scattering Characteristics of a Radially Multilayered Gyrotropic Sphere

**Academic Editor:**Safieddin Safavi-Naeini

#### Abstract

We present a new closed-form solution to the scattering of a monochromatic plane wave by a radially multilayered gyrotropic sphere using the -matrix method. This approach can be utilized to investigate the interactions of a plane wave and a gyrotropic spherical scatterer of multiple layers with each layer characterized by both permittivity and permeability tensors. Based on the completeness and noncoplanar properties of vector spherical wave functions (VSWFs), analytical expressions of the electromagnetic fields in each gyrotropic layer are first derived. The boundary conditions are then applied on each discontinuous interface to obtain the scattering coefficients. Validations are made by first comparing the radar cross section (RCS) values of a 2-layered gyrotropic sphere with that computed from the full-wave finite element method (FEM) simulation and then reducing the general case to uniaxial case to compare the RCS values with the published results computed by Fourier transform combined with VSWFs method; in both cases good agreements are observed. Several specific cases are fully explored to investigate how the RCS are influenced by the parameters of the multilayered spherical structure. The results show that when both electric and magnetic gyrotropy tensors are considered, the RCS of the multilayered spherical scatterer can be suppressed or enhanced, depending on proper configurations of the material parameters.

#### 1. Introduction

With the development of material science and fabrication techniques, the electromagnetic (EM) peculiarities of anisotropic media have attracted continued interests of both physics and engineering communities [1–7]. The effects of anisotropy or gyrotropy of the media must be taken into account in many practical applications, where a number of new devices with unique characteristics are designed by applying these effects [8, 9]. For studying the interactions of EM waves and gyrotropic media, the most important and challenging task is to solve the vector wave equation characterizing EM fields in gyrotropic media. In the last several decades, a few analytical approaches have been developed based on the Cartesian coordinate expression of permittivity and permeability tensors, such as the perturbation expansion method [10, 11], the method of expanding EM fields in gyrotropic media using a complete set of VSWFs [12], the method of adopting VSWFs combined with Fourier transform [13–15], dyadic Green’s functions method [16–18], and the -matrix method [19–23]. To study the spherical anisotropy, Qiu et al. have investigated the characteristics of radial anisotropic uniaxial and gyrotropic sphere, where the permittivity and permeability tensors are expressed in spherical coordinate [24–30]. Due to the completeness and noncoplanar properties of VSWFs, any three-dimensional vector function satisfying the vector Helmholtz equation can be expanded as a linear combinations of VSWFs [31]. In [20], Lin and Chui have developed a complete theory to solve Maxwell equations for gyromagnetic particles. This method has a good solution accuracy since it does not involve lengthy evaluations of complicated integrals. Li and Ong applied this method to study the EM scattering characteristics of a gyroelectric sphere in [21] and then generalized the scattering problem by considering both permittivity and permeability tensors in [22]. In [23], the homogeneous gyroelectric sphere considered in [20, 21] was extended to a multilayered gyroelectric sphere.

In this paper, the method developed in [23] is further extended to the most general case: the gyrotropic sphere under consideration consists of multiple anisotropic layers, each of which is characterized by both permittivity and permeability tensors. A new -matrix formulation is proposed for the spherical multilayered scatterer. It is shown from numerical results that the proposed computation procedure is robust and accurate for a wide range of material parameters. In this method, the EM fields in each layer are first expanded as an infinite summation of VSWFs based on the Mie scattering theory [20, 32]. The expansion coefficients in each layer and the scattering coefficients are then determined recursively by matching the boundary conditions at all interfaces. To validate the derived formulation, we compare two specific RCS results with those obtained by the FEM simulation and the Fourier transform combined with VSWFs method [33] for gyrotropic and uniaxial cases, respectively. Excellent agreements can be observed. The formulation developed in this paper can be used as a reference for other analytical approaches or numerical methods relying on discretizations. Finally, we calculate the RCS for several specific cases to provide some physical insights and possibilities of optimized designs for suppressing or enhancing RCS based on multilayered gyrotropic structures.

#### 2. Basic Formulations

Throughout the paper, a time dependence is assumed but always suppressed. The model of a radially multilayered gyrotropic sphere is depicted in Figure 1, where the radii of the spherical layers are denoted by . Region-0 is the free space characterized by scalar permittivity and permeability . Region- is homogeneous and characterized by both permittivity and permeability tensors, which are explicitly expressed in Cartesian coordinate aswhere .

##### 2.1. EM Fields in Multilayered Gyrotropic Sphere

In general, the EM fields in Region- of a multilayered gyrotropic sphere can be expressed as [22, 23] When , the summation subscript must take both 1 and 3 to obtain the complete solution of the vector wave functions, since the first-kind spherical Bessel function and the third-kind spherical Bessel function (or sometimes referred to as the first-kind spherical Hankel function) are regular in these regions; however, when , the summation subscript only takes 1 due to the singularity of the second-kind spherical Bessel function (and hence the third-kind spherical Bessel function) at the origin. A detailed derivation for obtaining the expansion coefficients and in (2a) and (2b) can be found in [22, 23] and will not be repeated here for simplicity.

##### 2.2. Matching Boundary Conditions

The EM boundary conditions on the interfaces are for , andfor . In the above equations, is the unit vector of the radial component in spherical coordinate, and are the scattered and incident fields expanded by VSWFs (see (26) and (28) in [23] for details). Based on the boundary conditions, we can obtain the following recurrent formulation:whereThe matrix elements are defined asAfter a careful comparison of () and () in [23], we observe that, by adding the gyromagnetic tensors, additional terms are introduced into the matrix elements of and . In (7), and are the Riccati-Bessel functions given bywith and being the spherical Bessel functions of the first and the third kind, respectively. According to the recurrence relations in (5), the expansion coefficients in Region-1 can be related to in Region- aswhere . Then we apply the EM boundary condition in (4a) and (4b) at to solve a scattering problem asIn (10a) and (10b), and are the scattering and incident coefficients (see detailed expressions in [23]). The matrices and are expressed as [20]where . By substituting (9) into (10a) and (10b), the scattering coefficients and and the expansion coefficients for Region- can all be solved numerically. After obtaining , the field expansion coefficients for internal layers (from Region-() to Region-1) can be readily calculated by applying the inward to outward recurrence procedure.

When , the asymptotic forms of the third-kind VSWFs are expressed as [34]where and are the unit vectors of the and components in the spherical coordinate, respectively, andConsequently, the RCS of the multilayered gyrotropic scatterer can be calculated by using the obtained scattering coefficients as

#### 3. Numerical Results and Discussions

The incident plane wave is assumed to be propagating along the positive direction, where the electric field is polarized in the direction with unit amplitude. The size parameters of the multilayered sphere are represented as , which are essentially the electrical dimensions with respect to the incident wavelength . In the practical computations, the summation index in (2a) and (2b) should be truncated at . As suggested in Appendix in [32] and () in [35], can be chosen to bewhere denotes the size parameter of the spherical scatterer. In this paper, is uniformly chosen as to make the matrices in (5) have the same dimensions. The linear equations (10a) and (10b) are solved independently of each due to the Kronecker symbols appearing in (11a) and (11b). In the following discussions, the numerical results of RCS are normalized by the square of the incident wavelength ().

The convergence criteria (15) are commonly used for homogeneous spherical scatterers [35, 36]. As for the gyrotropic case considered in this work, the validation of truncation number has been checked by computing the relative error of truncated series of the scattering efficiency factor [20, 36]. Explicitly, the relative error is defined as , where is the truncated series and is the converged value (calculated for a sufficiently large ). Four cases with different electrical dimensions are depicted in Figure 4. It can be seen that the truncated series converge well when the summation index is truncated at the value given by (15), and the convergence accuracy is at least .

To verify the correctness of the formulations derived for the multilayered gyrotropic sphere in this paper, two numerical experiments are conducted. First, the normalized RCS values of a 2-layered gyrotropic sphere are compared with that computed from the full-wave finite element method (FEM) simulation, as depicted in Figure 2. The normalized root mean square deviations (NRMSD) for -plane () and - plane () are 0.0041 and 0.0052, respectively. In the FEM simulation, a mesh with 53294 tetrahedral elements has been used and the basis functions are constructed from polynomials of order 3. Second, the gyrotropic case is reduced to the uniaxial case i.e., and , for , so that a direct comparison can be made with the published results (Figure in [33]), as shown in Figure 3. The NRMSD for - and -planes are 0.007 and 0.0099, respectively. In both cases, excellent agreements of the RCS values are achieved in both -plane and -plane, confirming the correctness of the newly derived formulations.

In the following, three specific cases will be presented to further investigate the scattering characteristic of the multilayered gyrotropic sphere.

In Figure 5, the bistatic RCS of a 2-layered lossless gyrotropic sphere with small-size parameters is presented, where two cases of different EM parameter settings are studied comparatively. The RCS results are shown in Figures 5(a) and 5(b). We observe that if the permittivity and permeability tensors in each layer satisfy the following conditions,the RCS in -plane and -plane coincide with each other. More importantly, a sharp RCS reduction in the backscattering direction can be obviously observed. Therefore, this property may have potential applications in stealth or cloaking design.

**(a) Case**

**(b) Case**

To further understand this scattering property, another 4-layered lossy gyrotropic sphere with medium-size parameters is considered and the bistatic RCS values are depicted in Figure 6. In this case, the parameters are chosen based on (16), and the property presented in Figure 5(b) still holds for this case. It can be possibly concluded that if the parameters in each layer satisfy (16), the RCS in the backward direction can be made too weak to be detected, no matter for lossless or lossy anisotropic materials.

In Figure 7, the bistatic RCS of a 2-layered lossless gyrotropic sphere with large-size parameters is presented. The EM parameters are configured specifically to give another insight of the scattering property. From Figure 7, it can be observed that the RCS becomes more oscillatory than those depicted in Figures 5 and 6. When the scattering angle increases from 0° to 180°, the RCS decreases first and then increases to a high level in both -plane and -plane.

#### 4. Conclusions

The analytical solutions of the EM scattering by a radially multilayered gyrotropic sphere are studied in this paper. First, by applying the completeness and noncoplanar properties of the VSWFs, theoretical formulations of this scattering problem have been derived. Second, we compare the numerical results of RCS values with those obtained from the full-wave simulation or from the literature to verify the validity and correctness of the proposed method. Third, specific cases are investigated to reveal some interesting physical insights of the scattering property of the multilayered gyrotropic sphere. It is found that if the gyrotropic tensors are configured properly, the RCS in the backward scattering direction can be suppressed or enhanced significantly.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

The authors are grateful to the useful discussions with Mr. Chao Wan at the National University of Singapore. L. Cao acknowledges the financial support from Special Talent Program at University of Electronic Science and Technology of China (UESTC). This is in memory of Professor Joshua Le-Wei Li, and taking this opportunity, the authors wish to express their genuine gratitude to him for his kind and patient guidance.

#### References

- J. Luo, W. Lu, Z. Hang et al., “Arbitrary control of electromagnetic flux in inhomogeneous anisotropic media with near-zero index,”
*Physical Review Letters*, vol. 112, no. 7, Article ID 073903, 2014. View at: Publisher Site | Google Scholar - Z.-J. Li, H.-Y. Li, and Z.-S. Wu, “Analysis of electromagnetic scattering by uniaxial anisotropic bispheres,”
*Journal of the Optical Society of America A: Optics and Image Science, and Vision*, vol. 28, no. 2, pp. 118–125, 2011. View at: Publisher Site | Google Scholar - Y. Petrov, “Anisotropic spherical head model and its application to imaging electric activity of the brain,”
*Physical Review E - Statistical, Nonlinear, and Soft Matter Physics*, vol. 86, no. 1, Article ID 011917, 2012. View at: Publisher Site | Google Scholar - K. Sainath and F. L. Teixeira, “Tensor Green's function evaluation in arbitrarily anisotropic, layered media using complex-plane Gauss-Laguerre quadrature,”
*Physical Review E - Statistical, Nonlinear, and Soft Matter Physics*, vol. 89, no. 5, Article ID 053303, 2014. View at: Publisher Site | Google Scholar - T.-J. Guo, T.-F. Li, M. Yang et al., “Nonreciprocal optical diffraction by a single layer of gyromagnetic cylinders,”
*Optics Express*, vol. 22, no. 1, pp. 537–546, 2014. View at: Publisher Site | Google Scholar - R. C. Haussman and M. Deserno, “Effective field theory of thermal Casimir interactions between anisotropic particles,”
*Physical Review E - Statistical, Nonlinear, and Soft Matter Physics*, vol. 89, no. 6, Article ID 062102, 2014. View at: Publisher Site | Google Scholar - J. Christensen and F. J. G. de Abajo, “Anisotropic metamaterials for full control of acoustic waves,”
*Physical Review Letters*, vol. 108, no. 12, Article ID 124301, 2012. View at: Google Scholar - Z. Wu, M. Levy, V. J. Fratello, and A. M. Merzlikin, “Gyrotropic photonic crystal waveguide switches,”
*Applied Physics Letters*, vol. 96, no. 5, Article ID 051125, 2010. View at: Publisher Site | Google Scholar - N. K. Dissanayake, M. Levy, A. A. Jalali, and V. J. Fratello, “Gyrotropic band gap optical sensors,”
*Applied Physics Letters*, vol. 96, no. 18, Article ID 181105, 2010. View at: Publisher Site | Google Scholar - G. W. Ford, J. K. Furdyna, and S. A. Werner, “Electromagnetic excitations of a small gyrotropic sphere,”
*Physical Review B*, vol. 12, no. 4, pp. 1452–1465, 1975. View at: Publisher Site | Google Scholar - R.-J. Tarento, K.-H. Bennemann, P. Joyes, and J. Van de Walle, “Mie scattering of magnetic spheres,”
*Physical Review E - Statistical, Nonlinear, and Soft Matter Physics*, vol. 69, no. 2, Article ID 026606, pp. 1–26606, 2004. View at: Publisher Site | Google Scholar - G. W. Ford and S. A. Werner, “Scattering and absorption of electromagnetic waves by a gyrotropic sphere,”
*Physical Review B*, vol. 18, no. 12, pp. 6752–6769, 1978. View at: Publisher Site | Google Scholar - 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 ID 056609, 2004. View at: Google Scholar - Y.-L. Geng, X.-B. Wu, and L.-W. Li, “Characterization of electromagnetic scattering by a plasma anisotropic spherical shell,”
*IEEE Antennas and Wireless Propagation Letters*, vol. 3, no. 1, pp. 100–103, 2004. View at: Publisher Site | Google Scholar - Y.-L. Geng, C.-W. Qiu, and N. Yuan, “Exact solution to electromagnetic scattering by an impedance sphere coated with a uniaxial anisotropic layer,”
*Institute of Electrical and Electronics Engineers. Transactions on Antennas and Propagation*, vol. 57, no. 2, pp. 572–576, 2009. View at: Publisher Site | Google Scholar | MathSciNet - A. Eroglu and J. K. Lee, “Dyadic green’s functions for an electrically gyrotropic medium,”
*Progress In Electromagnetics Research*, vol. 58, pp. 223–241, 2006. View at: Google Scholar - A. Eroglu and J. K. Lee, “Simplified formulation of dyadic green’s functions and their duality relations for general anisotropic media,”
*Progress In Electromagnetics Research*, vol. 77, pp. 391–408, 2007. View at: Publisher Site | Google Scholar - L.-W. Li, N.-H. Lim, W.-Y. Yin, and J.-A. Kong, “Eigenfunctional expansion of dyadic Green's functions in gyrotropic media using cylindrical vector wave functions,”
*Progress in Electromagnetics Research*, vol. 43, pp. 101–121, 2003. View at: Publisher Site | Google Scholar - B. Stout, M. Nevière, and E. Popov, “T matrix of the homogeneous anisotropic sphere: applications to orientation-averaged resonant scattering,”
*J. Opt. Soc. Am. A*, vol. 24, no. 4, pp. 1120–1130, 2007. View at: Google Scholar - Z. Lin and S. T. Chui, “Electromagnetic scattering by optically anisotropic magnetic particle,”
*Physical Review E*, vol. 69, no. 5, Article ID 056614, 2004. View at: Google Scholar - J. L. Li and W.-L. Ong, “A new solution for characterizing electromagnetic scattering by a gyroelectric sphere,”
*Institute of Electrical and Electronics Engineers. Transactions on Antennas and Propagation*, vol. 59, no. 9, pp. 3370–3378, 2011. View at: Publisher Site | Google Scholar | MathSciNet - J. L.-W. Li, W.-L. Ong, and K. H. R. Zheng, “Anisotropic scattering effects of a gyrotropic sphere characterized using the T-matrix method,”
*Physical Review E - Statistical, Nonlinear, and Soft Matter Physics*, vol. 85, no. 3, Article ID 036601, 2012. View at: Publisher Site | Google Scholar - L. Cao, J. L. Li, and J. Hu, “A new T-matrix formulation for electromagnetic scattering by a radially multilayered gyroelectric sphere,”
*Institute of Electrical and Electronics Engineers. Transactions on Antennas and Propagation*, vol. 61, no. 2, pp. 836–842, 2013. View at: Publisher Site | Google Scholar | MathSciNet - 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. View at: Publisher Site | Google Scholar - 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 - Statistical, Nonlinear, and Soft Matter Physics*, vol. 75, no. 2, Article ID 026609, 2007. View at: Publisher Site | Google Scholar - C.-W. Qiu, S. d. Zouhdi, and A. Razek, “Modified spherical wave functions with anisotropy ratio: application to the analysis of scattering by multilayered anisotropic shells,”
*Institute of Electrical and Electronics Engineers. Transactions on Antennas and Propagation*, vol. 55, no. 12, pp. 3515–3523, 2007. View at: Publisher Site | Google Scholar | MathSciNet - C.-W. Qiu and B. Luk'yanchuk, “Peculiarities in light scattering by spherical particles with radial anisotropy,”
*Journal of the Optical Society of America A*, vol. 25, no. 7, pp. 1623–1628, 2008. View at: Publisher Site | Google Scholar - C.-W. Qiu, L. Hu, X. Xu, and Y. Feng, “Spherical cloaking with homogeneous isotropic multilayered structures,”
*Physical Review E - Statistical, Nonlinear, and Soft Matter Physics*, vol. 79, no. 4, Article ID 047602, 2009. View at: Publisher Site | Google Scholar - 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 - Statistical, Nonlinear, and Soft Matter Physics*, vol. 80, no. 1, Article ID 016604, 2009. View at: Publisher Site | Google Scholar - C.-W. Qiu, L. Gao, J. D. Joannopoulos, and M. Soljacic, “Light scattering from anisotropic particles: propagation, localization, and nonlinearity,”
*Laser & Photonics Reviews*, 2010. View at: Google Scholar - D. Sarkar and N. J. Halas, “General vector basis function solution of Maxwell's equations,”
*Physical Review E. Statistical, Nonlinear, and Soft Matter Physics*, vol. 56, no. 1, part B, pp. 1102–1112, 1997. View at: Publisher Site | Google Scholar | MathSciNet - C. F. Bohren and D. R. Huffman,
*Absorption and Scattering of Light by Small Particles*, John Wiley & Sons, NY, USA, 1998. - Y.-L. Geng, S.-X. Guo, and L.-W. Li, “EM scattering by radially multilayered uniaxial anisotropic spheres,” in
*Proceedings of the Asia Pacific Microwave Conference 2009, APMC 2009*, pp. 669–672, sgp, December 2009. View at: Publisher Site | Google Scholar - Y.-L. Xu, “Electromagnetic scattering by an aggregate of spheres: far field,”
*Applied Optics*, vol. 36, no. 36, pp. 9496–9508, 1997. View at: Publisher Site | Google Scholar - W. J. Wiscombe, “Improved mie scattering algorithms,”
*Applied Optics*, vol. 19, no. 9, pp. 1505–1509, 1980. View at: Publisher Site | Google Scholar - J. R. Allardice and E. C. Le Ru, “Convergence of mie theory series: criteria for far-field and near-field properties,”
*Applied Optics*, vol. 53, no. 31, pp. 7224–7229, 2014. View at: Google Scholar

#### Copyright

Copyright © 2017 Lei Cao et al. 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.