Gain-assisted transformation optics

Loss severely degrades the cloaking effect of the device designed by traditional transformation. In this letter, we propose gain-assisted transformation optics to overcome the loss problem by introducing gain media into a spherical cloak. The gain media, which can amplify the electromagnetic fields, is controlled precisely to compensate the inherent loss in experimental realization of cloaks. We discuss the significance of controlling embedded gain materials in the context of the inverse design mechanism, which allows us to wisely select realizable materials with constant gain and loss along the radius. For practical realizations, isotropic spherical gain-assisted cloak is designed. Full-wave simulations validate the proposed design concept, which can be utilized to alleviate the inevitable loss problem in transformational optical devices. ©2011 Optical Society of America OCIS codes: (160.1190) Anisotropic optical materials; (160.3918) Metamaterials; (230.0230) Optical devices; (260.2065) Effective medium theory. References and links 1. U. Leonhardt, “Optical conformal mapping,” Science 312(5781), 1777–1780 (2006). 2. J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science 312(5781), 1780–1782 (2006). 3. J. Li, and J. B. Pendry, “Hiding under the carpet: a new strategy for cloaking,” Phys. Rev. Lett. 101(20), 203901 (2008). 4. Y. Lai, H. Y. Chen, Z. Q. Zhang, and C. T. Chan, “Complementary media invisibility cloak that cloaks objects at a distance outside the cloaking shell,” Phys. Rev. Lett. 102(9), 093901 (2009). 5. 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). 6. R. Liu, C. Ji, J. J. Mock, J. Y. Chin, T. J. Cui, and D. R. Smith, “Broadband ground-plane cloak,” Science 323(5912), 366–369 (2009). 7. J. Valentine, J. Li, T. Zentgraf, G. Bartal, and X. Zhang, “An optical cloak made of dielectrics,” Nat. Mater. 8(7), 568–571 (2009). 8. B. Zhang, H. Chen, and B. I. Wu, “Practical limitations of an invisibility cloak,” Prog. Electromagn. Res. 97, 407–416 (2009). 9. H. Hashemi, B. L. Zhang, J. D. Joannopoulos, and S. G. Johnson, “Delay-bandwidth and delay-loss limitations for cloaking of large objects,” Phys. Rev. Lett. 104(25), 253903 (2010). 10. B. Zhang, Y. Luo, X. Liu, and G. Barbastathis, “Macroscopic invisibility cloak for visible light,” Phys. Rev. Lett. 106(3), 033901 (2011). 11. X. Chen, Y. Luo, J. Zhang, K. Jiang, J. B. Pendry, and S. Zhang, “Macroscopic invisibility cloaking of visible light,” Nat Commun 2, 176 (2011). 12. A. A. Govyadinov, V. A. Podolskiy, and M. A. Noginov, “Active metamaterials: sign of refractive index and gain-assisted dispersion management,” Appl. Phys. Lett. 91(19), 191103 (2007). 13. S. M. Xiao, V. P. Drachev, A. V. Kildishev, X. J. Ni, U. K. Chettiar, H. K. Yuan, and V. M. Shalaev, “Loss-free and active optical negative-index metamaterial,” Nat. Lett. 466(7307), 735–738 (2010). 14. S. Wuestner, A. Pusch, K. L. Tsakmakidis, J. M. Hamm, and O. Hess, “Overcoming losses with gain in a negative refractive index metamaterial,” Phys. Rev. Lett. 105(12), 127401 (2010). 15. A. Novitsky, C. W. Qiu, and S. Zouhdi, “Transformation-based spherical cloaks designed by an implicit transformation-independent method: theory and optimization,” N. J. Phys. 11(11), 113001 (2009). 16. C. W. Qiu, L. Hu, X. F. Xu, and Y. J. Feng, “Spherical cloaking with homogeneous isotropic multilayered structures,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 79(4 Pt 2), 047602 (2009). 17. C. A. Balanis, Advanced Engineering Electromagnetics, (John Wiley & Sons, 1989). #143723 $15.00 USD Received 7 Mar 2011; accepted 28 Mar 2011; published 18 Apr 2011 (C) 2011 OSA 25 April 2011 / Vol. 19, No. 9 / OPTICS EXPRESS 8610 18. A. Boltasseva, and H. A. Atwater, “Materials science. Low-loss plasmonic metamaterials,” Science 331(6015), 290–291 (2011).

However, the inevitable loss will appear in the process of constructing artificial materials.Generally the losses are orders of magnitude too large for the proposed applications, and the reduction of losses with optimized designs seems to be out of reach.More recently, it is reported that the existence of loss will significantly frustrate the performance of invisibility cloaks [8,9].Namely, the loss problem plagues the invisibility cloaks and is one of the major restrictions preventing them from wider industrial applications.One solution is to use low-loss dielectric uniaxial materials as macroscopic cloaks experimentally verified in [10,11].However, the transformation applies to a particular geometry resulting in homogeneous uniaxial materials and also uses high-permittivity background medium to scale up the cloak's parameters (all bigger than, so as to use available calcites [10,11]).When inhomogeneities or parameters smaller than unity are inevitable, another possible solution to this bottleneck is the use of active materials, which can amplify the waves and compensate the loss [12][13][14].Unfortunately, the traditional transformation optics usually involve only ideal lossless passive media and cannot overcome the loss problem effectively.
In this work, gain-assisted transformation optics is proposed to overcome the loss problem of 3D spherical invisibility cloaks.Different from traditional transformation that usually involves only ideal lossless passive media, the proposed transformation involves lossy media with negative imaginary part and gain media with positive imaginary part.The gain-assisted transformation has the advantage of traditional transformation of controlling the wave transmission, but compensates its deficiency in dealing with loss problem.Once the loss in cloak is determined, the needed active media can be obtained according to the theory.By compensating the attenuation caused by the lossy parameters, it is, a priori, possible to realize an overall ideal cloak.

Theoretical analysis
Gain-assisted transformation optics is introduced with the assistance of inverse design method [15], in which the constitutive parameters can be independently obtained without any knowledge of the corresponding coordinate transformation.We adopt this inverse design because it allows us to effectively spot a simple gain material that is homogenous.The simulations are carried out in MathematicaTM.We know that a spherical cloak is obtained by compressing a spherical region ( in real space ( , , )  r   .The well-known expression for the constitutive parameters in the real space transformed from vacuum is where ( , , ) / ( , , ) J r r is the Jacobian matrix.The transformation is operated along the radial direction (   and  )           , and thus the Jacobian matrix is diagonal, though we still have no information of what the coordinate transformation is.The constitutive parameters can be expressed as Different from the previous study [15] that regarded constitutive parameters to be real, here the parameters are complex.By setting and substituting it into Eq.( 3), one can obtain Then we introduce ( ) Then the transverse parameter can be represented , the radial parameter can be determined Nevertheless, the unknown coordinate transformation can be found To validate the proposed method, we select a specific generating function ( ) g r = 1 and , where C is a constant.Its corresponding coordinate transformation can be expressed as Then the transverse and radial parameters can be derived Obviously, If C = 0, Eq. ( 9) degenerates into a traditional transformation, and the complex constitutive parameters in Eq. ( 10) degenerate into real parameters [1].

Simulation results and discussion
When a = 0.1 m and b = 0.2 m in Eq. ( 10), Fig. 1 shows the complex permittivity composed of real part and imaginary part.The imaginary part can be negative and positive constants corresponding to lossy or active materials respectively.Obviously, attenuation and amplification of the EM field will coexist in the realized cloak.Finally, the gain-assisted cloak would behave effectively like a lossless one, though the constitutive parameters are non-ideal.From Fig. 1, it is found that the selection of the value of C significantly affects the material parameters, especially the imaginary part of the parameters.Increasing C may lead to bigger imaginary parts of the parameters, and if C approaches zero, the constitutive parameters will approach the values obtained by the traditional transformation.More importantly, the imaginary part of the transverse parameter keeps homogeneous in Fig. 1(d), which represents the maximum loss or gain according to Eq. ( 10)a).In what follows, we show the full-wave simulations to validate the proposed method.The incident plane wave is along z-axis and the frequency is 2 GHz.A dielectric spherical core is placed in cloaking region with refractive index n = 1.45.Figure 2 shows the total electric field distributions of the gain-assisted cloak with C = 0.5 in x-z plane (Fig. 2(a)) and y-z plane (Fig. 2(b)).As predicted, the EM waves smoothly flow around the cloaked region and then restore their original propagation paths.However, compared to the traditional cloak [16], it is observed that the field distribution is asymmetric in the shell region.This is because the gainassisted cloak is composed of loss materials and gain materials, hence the phenomena of loss (corresponding to the field attenuation) and compensation (corresponding to the field amplifying) in the shell region can be found.In order to examine the compensation effect quantitatively, we calculate the cross section of the gain-assisted cloak in Fig. 2 with different amounts of gain embedded, as shown in Fig. 3(a).Obviously, the fully compensated cloak (corresponding to 100% gain) performs much more perfectly than the non-compensated cloak (corresponding to 0% gain), though both of them have very small backscattering which accords very well with the theoretical results.When the gain is introduced insufficiently (corresponding to 50% gain) or overly (corresponding to 150% gain), the backward scattering is very small though evident scattering is induced in the forward direction.To give more insights, the scattering cross sections (σ) normalized by that of the non-compensated cloak (σL) are plotted versus the percentage of the gain (from 0 to 200%) in Fig. 3(b).Diverging from the ideal compensation point (100% gain), σ always increases regardless of that the gain is either decreasing or increasing.Nevertheless, the insufficient compensation makes σ increase much faster than the over compensation.Approximately above 40% gain, the total scattering cross section can decrease more than 10 dB compared to the non-compensated case.It should be noted that if it is difficult to introduce exactly 100% gain, it is better to introduce slightly more gain.The spherical gain-assisted cloak could be easily realized through alternating layered isotropic mediums [16].According to effective medium theory, the shell region is equally discretized into N layers, and each layer consists of medium-A and medium-B with equal thickness.Here we introduce two kinds of parameter profiles including old set and new set.The old set of material parameters for isotropic medium-A and medium-B is respectively.It is obvious that the old set corresponds to impedance matching and the new set corresponds to refractive index matching.Figure 4(a) illustrates the real part and imaginary part of medium-A and medium-B.It is found that the parameters of medium-A and media-B keep nearly invariant along the radius, and this property will make the fabrication easier.We also calculate the differential cross section of the isotropic cloak in Fig. 4(b).Clearly, the forward scattering of the old set is large, which can be improved much more by considering the new set.Although there is a small sacrifice in backscattering of the new set compared with the backscattering of the old set, it is evidently worthwhile since the backscattering of the new set is still negligible.It should be noted that the performance is quite pronounced when the discretization increases.

Conclusion
To conclude, we propose a general transformation that can degenerate into traditional transformation, namely gain-assisted transformation optics, which provides a powerful theoretical tool to overcome the loss problem in invisibility cloaks.Although the current gain materials are not enough to compensate large metallic loss [18], certain heavily doped semiconductors shed the bright light to this problem owing to the low loss.Nevertheless, this paper theoretically considers material's inevitable loss (high or low) into the design process and compensates it by introducing the gain of active materials properly, though the model is still magnetic now.The current method is not limited to the design of invisibility cloak, but can also be widely used in designing other transformational devices, such as perfect lenses, rotators, and illusion devices.It paves an alternative path of alleviating the inevitable loss problem in practical transformational devices.

Fig. 1 .
Fig. 1.Permittivity along the radial direction with different C. (a) Real part of radial parameter.(b) Imaginary part of radial parameter.(c) Real part of transverse parameter.(d) Imaginary part of transverse parameter.The permeability curve will be exactly the same since r r    and

Fig. 2 .
Fig. 2. Snapshots of the total electric field for the spherical gain-assisted cloak with 0.5 C  .(a) x-z plane, (b) y-z plane.

Fig. 3 .
Fig. 3. (a) Differential cross section (exactly the same definition in Ref [15]) for the gainassisted cloak with different compensation when C = 0.5.(b) The normalized total scattering cross sections σ/σL (normalized by σL of the loss cloak with non-compensation) as a function of the gain introduced.The scattering cross section is defined exactly the same as that in Ref [17].

Fig. 4 .
Fig. 4. (a) Constitutive parameters of media-A and media-B along radial direction for isotropic cloak with C = 0.5.(b) Differential cross section for isotropic cloaks designed by the old set and new set parameter profiles with different layer number (N = 20 and N = 40).
r to generate the real part of the transverse parameter, which is