Cylindrical optimized nonmagnetic concentrator with minimized scattering

By using optimized transformation function, we research on a minimized scattering nonmagnetic concentrator, which can realize impedance matching at the inner and the outer boundaries. It has been demonstrated that the optimized transformation function method can improve the concentrating performance remarkably. The cylindrical anisotropic shell can be mimicked by radial symmetrical sectors which alternate in composition between two profiles of isotropic dielectrics, and the permittivity in each sector can be properly determined by the effective medium theory. The nonmagnetic concentrator has been validated by fullwave finite element simulations. We can believe that this work will improve the flexibilities for the EM concentrator design. ©2013 Optical Society of America OCIS codes: (290.5839) Scattering, invisibility; (220.1770) Concentrators; (050.2065) Effective medium theory. References and links 1. J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science 312(5781), 1780–1782 (2006). 2. U. Leonhardt, “Optical conformal mapping,” Science 312(5781), 1777–1780 (2006). 3. D. R. Smith, W. J. Padilla, D. C. Vier, S. C. Nemat-Nasser, and S. Schultz, “Composite medium with simultaneously negative permeability and permittivity,” Phys. Rev. Lett. 84(18), 4184–4187 (2000). 4. D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science 314(5801), 977–980 (2006). 5. D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express 14(21), 9794–9804 (2006). 6. 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,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 80(1), 016604 (2009). 7. Y. Huang, Y. Feng, and T. Jiang, “Electromagnetic cloaking by layered structure of homogeneous isotropic materials,” Opt. Express 15(18), 11133–11141 (2007). 8. M. Farhat, S. Guenneau, and S. Enoch, “Ultrabroadband elastic cloaking in thin plates,” Phys. Rev. Lett. 103(2), 024301 (2009). 9. C. W. Qiu, L. Hu, X. Xu, and Y. Feng, “Spherical cloaking with homogeneous isotropic multilayered structures,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 79(4), 047602 (2009). 10. C. W. Qiu, L. Hu, and S. Zouhdi, “Isotropic non-ideal cloaks providing improved invisibility by adaptive segmentation and optimal refractive index profile from ordering isotropic materials,” Opt. Express 18(14), 14950–14959 (2010). 11. W. Cai, U. K. Chettiar, A. V. Kildishev, and V. M. Shalaev, “Optical cloaking with metamaterials,” Nat. Photonics 1(4), 224–227 (2007). 12. W. Yan, M. Yan, and M. Qiu, “Non-magnetic simplified cylindrical cloak with suppressed zeroth order scattering,” Appl. Phys. Lett. 93(2), 021909 (2008). 13. X. F. Xu, Y. J. Feng, L. Zhao, T. Jiang, C. H. Lu, and Z. Z. Xu, “Designing the coordinate transformation function for non-magnetic invisibility cloaking,” J. Phys. D Appl. Phys. 41(21), 215504 (2008). 14. J. Zhang, Y. Luo, and N. A. Mortensen, “Minimizing the scattering of a nonmagnetic cloak,” Appl. Phys. Lett. 96(11), 113511 (2010). 15. L. Peng, L. Ran, and N. A. Mortensen, “The scattering of a cylindrical invisibility cloak: reduced parameters optimization,” J. Phys. D Appl. Phys. 44(13), 135101 (2011). #176873 $15.00 USD Received 26 Sep 2012; revised 11 Jan 2013; accepted 11 Jan 2013; published 29 Jan 2013 (C) 2012 OSA 11 March 2013 / Vol. 21, No. S2 / OPTICS EXPRESS A231 16. L. J. Huang, D. M. Zhou, J. Wang, Z. F. Li, X. S. Chen, and W. Lu, “Generalized transformation for nonmagnetic invisibility cloak with minimized scattering,” J. Opt. Soc. Am. B 28(4), 922–928 (2011). 17. L. J. Huang, D. M. Zhou, J. Wang, Z. F. Li, X. S. Chen, and W. Lu, “Scattering characteristics of nonmagnetic invisibility cloak with minimized scattering,” Opt. Commun. 284(24), 5523–5530 (2011). 18. M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, D. R. Smith, and J. B. Pendry, “Design of electromagnetic cloaks and concentrators using form-invarant coordinate transformations of Maxwell’s equations,” Photon. Nanostruct. Fundam. Appl. 6(1), 87–95 (2008). 19. J. J. Yang, M. Huang, C. F. Yang, Z. Xiao, and J. H. Peng, “Metamaterial electromagnetic concentrators with arbitrary geometries,” Opt. Express 17(22), 19656–19661 (2009). 20. W. Wang, L. Lin, J. X. Ma, C. T. Wang, J. H. Cui, C. L. Du, and X. G. Luo, “Electromagnetic concentrators with reduced material parameters based on coordinate transformation,” Opt. Express 16(15), 11431–11437 (2008). 21. W. X. Jiang, T. J. Cui, Q. Cheng, J. Y. Chin, X. M. Yang, R. P. Liu, and D. R. Smith, “Design of arbitrarily shaped concentrators based on conformally optical transformation of nonuniform rational B-spline surfaces,” Appl. Phys. Lett. 92(26), 264101 (2008). 22. K. Zhang, Q. Wu, J. H. Fu, and L. W. Li, “Cylindrical electromagnetic concentrator with only axial constitutive parameter spatially variant,” J. Opt. Soc. Am. B 28(6), 1573–1577 (2011). 23. S. Xi, H. S. Chen, B. Zhang, B. I. Wu, and J. A. Kong, “Route to low-scattering cylindrical cloaks with finite permittivity and permeability,” Phys. Rev. B 79(15), 155122 (2009). 24. B. Wood, J. B. Pendry, and D. P. Tsai, “Directed subwavelength imaging using a layered metal-dielectric system,” Phys. Rev. B 74(11), 115116 (2006). 25. W. Cai, U. K. Chettiar, A. V. Kildishev, and V. M. Shalaev, “Designs for optical cloaking with high-order transformations,” Opt. Express 16(8), 5444–5452 (2008).


Introduction
In the past few years, the technique of transformation optic [1,2] is a powerful and systematic approach in manipulating electromagnetic waves and producing new functionalities with artificial engineered structures [3,4]. A series of applications have been reported based on this theoretical framework. Among all these applications, free-space invisibility cloaks [5][6][7][8][9][10] have definitely received most scientific interests. To achieve invisibility cloak at optical frequencies, Cai et al. proposed a nonmagnetic cloak with simplified material parameters for transverse magnetic (TM) waves [11]. Compared to the perfect set, such a nonmagnetic cloak has a drawback that significant scattering outside the cloak exists due to a mismatched impedance at the outer boundary. In order to match the wave impedance at the outer boundary, various coordinate transformation functions [12][13][14][15], such as quadratic, power, and piece power functions, have been applied to design the nonmagnetic cloak. Huang et al. proposed the generalized transformation functions and a kind of transformation functions with two parameters for the nonmagnetic invisibility cloak [16,17]. The purpose of using these optimized functions is to avoid the reflection at the outer boundary of the nonmagnetic cloak. Recently, the optical transformation theory has been applied to the electromagnetic (EM) concentrator, and some theoretical analyses, numerical simulations, and parameter designs have been devoted to EM concentrating devices [18][19][20][21][22]. It is expected that the electromagnetic concentrator design possesses potential applications in optical microscopes, solar cells, and devices in which high field intensities are needed [22]. However, nonmagnetic concentrator with simplified material parameters also has significant scattering outside the structure, which is similar to that of nonmagnetic cloak. Thus, how to minimize the scattering of the nonmagnetic concentrator and improve the concentrating performance is highly urgent. In the process of design of optimized nonmagnetic concentrator, in order to minimize the scattering outside, the impedance of the inner and the outer boundaries must be considered simultaneously. The increasing of limited conditions will induce more complicated optimization process. Although people have achieved success in the minimizing the scattering of invisibility cloak theory, how to reduce the scattering of the electromagnetic concentrator has not been revealed.
In this paper, we research on how to further improve the impedance-matched nonmagnetic concentrator by properly designing the coordinate transformation function. The appropriate transformation function is given with the matched impedance at the inner and the outer boundaries of the nonmagnetic concentrator. With the aid of far field scattering patterns, we show the full-wave solution of the scattering by nonmagnetic concentrator. We can find that this kind of concentrator has the perfect performance. Moreover, we find an approach to realize the practical fabrication, in which radial symmetrical layered structures of isotropic materials instead of using the artificial metamaterials with subwavelength structured inclusions. Just by properly choosing the permittivity in each isotropic dielectric layer, excellent concentrating performance can be achieved, which will improve the concentrator design.

Principle of transformation
For simplification purposes, we consider a bi-dimensional case of the nonmagnetic concentrator in a cylindrical coordinate system as schematically illustrated in Fig. 1 , namely, the circular region. Here, r and r represent the radii of the physical space and the virtual space, respectively.
(1) According to the standard procedure of transformation optics, the relative permittivity and the permeability after the transformation can be written in terms of () fr as Here, in order to remove the magnetic response of the material, which is especially important for making devices at optical frequency, the constitutive parameters of the proposed nonmagnetic concentrator for TM illumination can be further simplified as 2 2 () , Compared to the ideal properties as depicted in Eq. (2), the reduced parameters can provide the similar wave trajectory inside the cylindrical region. According to Eq. (3), the wave impedance at the boundary becomes which depends on the transformation function.
According to the theory in Ref.
[18], the transformation function can be easily selected as the linear function in core region Through the transformation function of Eq. (4), we observe that the enhancing ratio can be expressed as the ratio of 2 R to 1 R , and enhancement theoretically diverges to infinity as 1 R tends to zero. In circular region, the linear transformation function can be selected as Then the reduced parameters of the core and circular region can be derived as Parameters presented in Eq. (6) minimize the most challenging part of experimental design. But this set of parameter specification leads to impedance mismatch at the inner and outer boundary. The reduced parameters in core region have the impedance of in circular region at the two boundaries. So the design approach using linear transformation function will not only induce some unfavorable scattering at the inner and the outer boundaries due to the impedance mismatch, but also lead to energy scattering. Meanwhile, the concentrating performance is degenerated based on this design approach.
To improve the concentrating efficiency, we need to derive another set of simplified material tensors, which also have simple tensor components and make the inner and the outer boundaries perfectly matched. In order to achieve the nonmagnetic concentrator with minimized scattering, the transformation functions are required to satisfy following conditions:

Simulation and discussion
We notice that the choice for the transformation function which satisfy these restrictions is not unique, which indicates that the optimization might be accomplished by using some numerical methods like genetic algorithm [15,23]. If the variations in the distribution of the material parameter values along the radius are gently, the scattering of the proposed concentrator is smaller. We find that the high polynomial functions which satisfy Eq. (7) where, 33  2  1  1  3  1  3  1 3  1  3  2  12  3  3  3  1  3  1  3  1  3 3 (8) and (9) into Eq. (3), we can obtain   Figure 2 compares the results of two transformation functions: the linear and the polynomial function. Here the geometry parameters are selected as 3 We can find that the polynomial function curve is even smoother than the linear function curve, especially, on the inner boundary which can realize the impedance matched on inner boundary perfectly. Based on the above results, the constitutive parameters of the nonmagnetic concentrator can be calculated through Eqs. (10) and (11), as illustrated in Fig. 3. It can be seen that the permittivity values of different regions on the inner boundary ( 1 rR  ) are identical. We can expect that the concentrator with the new set of parameters would induce smaller scattering. Moreover, all material parameters have finite values. It is indicated that the materials are relatively easily realized with the metamaterials.  To compare the performance of the nonmagnetic concentrators with different transformation functions, we calculate the EM field scattering using the finite-elementmethod (FEM) based full-wave EM simulation with the commercial COMSOL Mutiphysics package. We consider the TM incident wave with the magnetic field polarized along the z axis. The TM wave propagates along the x direction. The geometry parameters are also selected as 3 2 1 2 3 0.6 , R R R m    and the incident wave frequency is set at 1.5GHz. The results of the spatial distributions of the magnetic field and the corresponding power flow distribution are shown in Fig. 4. Figure 4(a) shows the magnetic field distribution of concentrator with linear transformation. It is obvious that part of magnetic field is disturbed by mismatched impedance when crossing the concentrator, so the scattering field obviously emerges in the outside region. Figure 4(c) shows the magnetic field distribution with polynomial transformation function, it can be seen that the magnetic field is smoothly concentrated into the inner core region, and the field outside is undisturbed. Furthermore, the corresponding power flow distribution with different transformations are calculated and shown in Fig. 4(b) and Fig. 4(d). Comparing the case of linear function with the case of polynomial function, the difference between the energy distribution of the two cases are conspicuous. For the case of linear transformation, part of the energy is scattered by the concentrator, which results in the disturbance of the energy in the whole space. As for the case of polynomial transformation, although the distribution of the energy in the core region is somewhat inhomogeneous, most of energy is concentrated.
In order to quantify the perfect concentrating performance of the proposed concentrator, we calculate the intensity of power flow on the centre line of core region with linear and polynomial transformation function, which is shown in Fig. 5. Comparing the two cases with different transformation functions, we can find that the concentrating performance with polynomial transformation is much better than that of linear transformation. To be specific, the normalized power flow of optimized concentrator is close to 1, but the maximal normalized power flow value of linear transformed concentrator just reaches 0.8 and what is more, the distribution of the energy is highly non-uniform. Comparing curves in Fig. 5, we can find that the method using polynomial transformation function can enhance the concentrating performance by 25%. In order to make a further comparison, we plot the normalized far field scattering patterns of the nonmagnetic concentrator with different transformation functions, as shown in Fig. 6. By comparing the two cases with different transformation functions, we find that the backscattering reduction property can be distinctly improved by choosing the optimized polynomial transformation function other than the conventional linear transformation function. Based on the effective medium theory [24], a series of cylindrical or spherical anisotropic invisible cloaks have been mimicked by concentric layered structures consisting of alternating layers of two different homogenous isotropic dielectric materials in subwavelengh scale. Next, we restrict the practical fabrication of the nonmagnetic concentrator device. It can be seen from Eqs. (10) and (11) that the constitutive materials are all inhomogeneous and anisotropic in core and circular regions, which should be approximated by the structure of homogeneous and isotropic materials based on effective media theory.
Furthermore, from the parameters derived in Eqs. (10) and (11), we can find that the radial component is always larger than the angular components ( r    ) which is different from all the previous cases of invisibility cloaks [12][13][14][15]. Such required cylindrical anisotropy shell can be mimicked by utilizing radial symmetrical sectors alternate in composition between two profiles of isotropic dielectrics instead of using the concentric alternating layers.
Supposing the tangential width of the dielectric profile-A and profile-B is identical, we find a new mapping of parameters between the cylindrical anisotropic shell and isotropic layered shell in interlaced profile-A and profile-B as  for the two dielectric profiles. Figure 7 illustrates the permittivity patterns of the two interlaced dielectric profiles obtained by radial discretizing the circular region into 2 160 M  constituents in the core region and the circular region. It can be seen that the value of permittivity for both profiles is positive, indicating the desirable cylindrical anisotropy can be effectively demonstrated by nonmagnetic right-handed material. In order to validate such an isotropic layered concentrator device, the simulation results illuminated with a plane wave are shown in Fig. 8 and Fig. 9. Figure 8 shows that the transverse magnetic field of ( 2 80 M  ) layered nonmagnetic concentrator with the linear and the polynomial transformation. It is found that the scattering of the nonmagnetic concentrator using optimized function is much smaller than that of linear transformation.   Figure 9 shows the comparison of different coordinate transformation functions of ( 2 160 M  ) layered nonmagnetic concentrator. Following the comparison between Fig. 8(b) and Fig. 9(b), we can see that the concentrating performance of the nonmagnetic layered shell is quite pronounced when the discretizing number is large enough. Furthermore, from Fig.  9(b), the performance of the sufficient discretizing ( 2 160 M  ) layered nonmagnetic concentrator using optimized polynomial function is remarkably close to the performance of ideal optimized nonmagnetic concentrator which is shown in Fig. 4(c). Therefore, we can conclude that the method of using optimized transformation function can generate the concentrating performance remarkably.

Conclusion
In summary, we have proposed an optimized transformation function for nonmagnetic concentrator with minimized scattering. Such transformation function offers us a great deal of flexibility in designing concentrators of high performance. Based on the effective medium theory, we present a radial symmetrical layered system in composition of two profiles as nonmagnetic concentrator device. Finite element method has been carried out to demonstrate the performance of the nonmagnetic isotropic layered concentrator, and excellent performance has been observed, respectively.