# Conditions for achieving invisibility of hyperbolic multilayered nanotubes

###### Abstract

Highly-anisotropic plasmonic nanotubes exhibit a dramatic drop of the scattering cross section in the transition regime from hyperbolic to elliptic dispersion. The characterization of a realistic multilayered metamaterial is typically carried out by means of an effective medium approach providing average components of the permittivity tensor and wave fields. Here, the edge effects of the metal-dielectric stratified nanotube for different combinations were thoroughly analyzed. We show how the boundary layers, which in principle remain fully irrelevant in the estimation of the effective permittivity of the nanotube, however play a critical role in the resonant scattering spectra and the near field patterns. A remarkable enhancement of the scattered wave field is unexpectedly experienced at the frequencies of interest when a dielectric layer is chosen to be in contact with the cavity core.

###### keywords:

Hyperbolic Metamaterial, Invisibility, Plasmonics.^{†}

^{†}journal: Optics Communications

## 1 Introduction

Metamaterials with unprecedented optical properties have been recently proposed by several groups for applications in cloaking and invisibility. A vast majority of cloaking nanostructures are engineered in the basis of transformation optics principles Leonhardt06 ; Pendry06 ; Cai07 enabling analytical expressions of the electromagnetic fields and metamaterial properties Forouzeshfard15 ; Yu15 ; Rajput16 . An original concept suggested by Alu and Engheta has also attracted great attention which relies on the use of metamaterial (or metal) coatings to severely drop the scattering efficiency of an object by means of a nonresonant scattering-cancellation approach Alu05 . The strategy of such scattering cancellation takes advantage of the local negative polarizability of metamaterials, and its experimental realization was first demonstrated at microwave frequencies by using an array of metallic fins which are embedded in a high-permittivity environment to generate a metamaterial cloaking shell Edwards09 .

Novel extensions of the previous concepts have more recently made the scene, which are based on the use of double-shell and multilayered plasmonic coatings Tricarico09 ; Zhou16 . For instance, using plasmonic shells with an epsilon-near-zero material enables to reduce substantially the scattering losses and simultaneously providing the shielding of the cloaked region Alu08 ; Filonov12 . As an alternative, a nanotube consisting of a periodic distribution of metal and dielectric alternating layers, where the stratified metamaterial was described as a radial-anisotropic hyperbolic medium, has recently demonstrated narrow-band ultra-low scattering Kim15 . The invisibility spectral band occurs when one of the components of the effective permittivity tensor is near zero. The effective medium theory was adopted to efficiently reproduce the results provided by the analytical Lorenz-Mie scattering method Bussey75 ; Bohren98 .

Here, we study in detail the edge effects of the stratified hyperbolic nanotube employed in Ref. Kim15 for different combinations, that is, when the layer in contact with the environment medium is either the metal or the nonconducting material. We demonstrated that boundary effects play a relevant role in the resulting scattering efficiency of the nanotube, potentially clearing away the characteristic invisibility of the nanotube and boosting additional plasmonic resonances in the visible. Explicitly, when silver is set as the constituent material of the outermost layer, unexpectedly, a significant enhancement of the scattered signal is observed. As a consequence, the effective medium theory enabling a simplified model for the optical characterization of the nanoparticle may apparently lead to fallacious estimations even for metamaterials composed of subwavelength layers with a few-nanometers width.

## 2 The hyperbolic multilayered nanotube

Let us consider a cylindrical shell formed by a multilayered nanostructure as illustrated in Fig. 1. Following Ref. Kim15 , silver and titanium dioxide was used for the stratified medium. A silver nanofilm of width is set by the side of a TiO layer of width , thus forming the unit cell of periodic distribution along the radial coordinate. In our numerical simulations we considered a nanostructured tube of outermost diameter and composed of six layers, giving a total shell thickness of . The relative permittivities of silver and TiO within the visible range of frequencies may be approximated by Kim15

(1) |

and

(2) |

respectively. In the previous equations, the working wavelength is set in micrometers. In this study we examined hollow cylinders immersed in air.

Within a long-wavelength approach, the multilayered metamaterial may be considered as a uniaxial medium whose optic axis is set along the direction of periodicity, that is the radial axis in our particular case. The effective anisotropic medium is then optically characterized by a local permittivity tensor of component Yeh88

(3) |

along the optic axis, and

(4) |

in the perpendicular direction. In the previous equations, the metal filling factor

(5) |

represents the unique geometrical parameter determining the effective permittivities of the metamaterial, regardless the internal one-dimensional distribution of the materials involved in the unit cell.

In Fig. 2 we represent the real values of and within the visible range of frequencies, for different values of the metal filling factor. We observe that the real part of is near zero around a filling factor , whereas the real part of is near a pole around a filling factor . The hyperbolic regime is determined by the condition , indicated as shaded regions in Fig. 2. The choice corresponds to the so-called Type I hyperbolic metamaterials, whereas the choice denotes a Type II hyperbolic medium Ferrari15 . Note that when occurring at a wavelength nm, both and coincide at a value , denoting that the silver layers and the TiO layers have the same width. In this case, the hyperbolic regime spans the whole spectrum. Let us point out that the extraordinary dispersion features of hyperbolic and epsilon-near-zero metamaterials have inspired us in a plethora of applications such as subwavelength imaging Zapata11e ; Zapata12a , surface-wave engineering Miret10 ; Zapata13b , and double refraction Zapata14b ; Diaz16 , to mention a few.

## 3 Results and discussion

To estimate analytically the scattering efficiency of the multilayered Ag-TiO nanotube we followed the Lorenz-Mie scattering method Bussey75 . Note that the axis of the cylindrical cavity is oriented along the axis. First we assumed that the nanotube is illuminated by a TE-polarized plane wave propagating along the axis, as illustrated in Fig. 1. The magnetic field of the incident plane wave may be set as

(6) |

where is a constant amplitude and is the wavenumber in air. In this case, the scattered magnetic field may be set as

(7) |

where and are the radial and azimuthal cylindrical coordinates, respectively, and is the Hankel function of the first kind and order . By using the Lorenz-Mie theory Bohren98 , it is possible to evaluate analytically the scattering coefficients , which provide the estimation of the scattering efficiency as

(8) |

Plasmonic resonances are determined by the poles of the coefficients . On the other hand, the invisibility condition is established provided that the scattering coefficients arrive simultaneously to a value near zero.

In Fig. 3 we show the scattering efficiency of our Ag-TiO nanotube calculated with Eq. (8), provided that the hyperbolic metamaterial has a metal filling factor and that the hollow particle is illuminated by a TE polarized plane wave. In the numerical simulations we set a layer width and of 10 nm and a cylinder diameter of nm. Equation (8) is accurately estimated by using a finite number of scattering coefficients satisfying . As illustrated in Fig. 3, the calculated exhibits significant discrepancies when the inmost layer is either Ag or TiO, despite and remain the same in both cases. When the interior layer is made of silver, the scattering efficiency of the nanostructure reaches a minimum of at nm, as reported in Ref. Kim15 , where and . In addition, a high peak of is found at nm, where and . However, the number of resonant peaks increases if titanium dioxide is set in the inmost layer. In this case, the highest peak is located at a wavelength of 560 nm where the scattering efficiency yields , and a secondary strong peak is found at nm (). In addition, its minimum in scattering efficiency is shifted at nm where importantly it reaches a value of . We point out that such minimum in scattering efficiency, which is far of being associated with invisibility, is one order of magnitude higher than that estimated previously for an Ag inmost coating. Finally, a dissimilar behavior in scattering is also manifested for TM-polarized incident light (not shown in Fig. 3).

In order to provide a physical insight of such a behavior, we analytically calculated the electromagnetic fields inside the nanotube and in the air core and in the environment medium. In Fig. 4 we represent the modulus of the magnetic field, , at the invisibility wavelength nm for different geometric configurations, calculated by means of the Lorenz-Mie scattering method, provided that the incident wave field is TE polarized. When silver is set in the inmost layer of the hyperbolic metamaterial, we observe a moderate enhancement of the field at three different Ag-TiO interfaces, as shown in Fig. 4(a), which corresponds to a collective excitation of surface plasmon polaritons (SPPs). It is clear that the field exhibits strong irregularities derived from nonlocal effects in the multilayered metal-dielectric nanostructure Elser07 . Though such phenomenon severely puts into question the validity of the effective medium approach followed in Ref. Kim15 , the deviations from an average field inside the hyperbolic nanotube are moderate in this case. Analyzing the case that titanium dioxide is placed in contact with the interior medium, however, we observe even abrupter variations of the field especially at the rear part of the nanotube (in the semi-space ), as shown in Fig. 4(b). Furthermore, now the field distribution is clearly asymmetric with respect to the origin of coordinates, leaving a strong backscattered signal.

Finally, we observe from Fig. 4 that the field amplitude at the center of the nanotube and in its neighborhood is, in average, near the field remaining in the environment medium rather than zero. This fact evidences that the shielding effect that commonly is recognized in epsilon-near-zero core-shell nanocylinders Filonov12 cannot be reproduced in our hyperbolic metamaterial nanotube when the light impinges under TE polarization.

The following analyzes the effects of increasing the number of layers in our nanotube maintaining the thickness nm and the metal filling factor of the shell, implying that the layer width and of the metamaterial is reduced. Considering 12 alternating layers of silver and titanium dioxide, the layer width (and the period of the multilayered nanostructure) is reduced to one half with respect to the cases analyzed above and illustrated in Figs. 3 and 4. The scattering efficiency of such additional nanotube is shown in Fig. 5, considering the examples where the inmost layer is either Ag or TiO; the efficiency is again evaluated by means of the Lorenz-Mie scattering method. When the layer in contact with the core is silver, we find a minimum of equal to at a wavelength nm, representing an exiguous variation in comparison with the hyperbolic nanotube of 6 layers. In addition, the maximum of scattering efficiency reaches a value of at nm. If now we consider a nanotube with the interior layer made of TiO, its scattering efficiency at the invisibility wavelength nm reaches a value of , representing a factor of 25 with respect to the case previously analyzed. In fact, the minimum of the scattering efficiency is , which is found at nm. The overall shape of , however, becomes relatively similar in both cases.

The profile of the magnetic field at the new invisibility wavelength nm again exhibits abrupt changes at the metal-dielectric interfaces, as shown in Fig. 6(a) for a nanotube with inmost layer made of Ag, but the envelope enabling the estimation of an average value presents a slow variation inside the metamaterial. Moreover, a minimum of intensity is found near the origin of coordinates, however such intensity is far of zero disabling the gain of a shielding effect. In Fig. 6(b) we represent the magnetic field of the nanotube at the invisibility wavelength of nm, considering a nanotube of interior layer made of TiO. The resonant excitation of collective SPPs at the rear part of the nanocylinder, developing a chief (inhomogeneous) Bloch wave in the multilayered metamaterial Torrent09 , seems to be responsible of the significant increase in the scattering efficiency and subsequently the predominant visibility of the nanotube. Now this is clear that the appearance of a predominant radially-evanescent Bloch wave also occurs for a 6-layer nanostructure, as shown in Fig. 4(b).

Obviously, by including a higher number of layers and, at the same time, maintaining the thickness of the nanotube and the filling factor of the metamaterial, both the scattering efficiency spectrum and the field pattern will continuously approach for the Ag and TiO ending layers, leading to the response of a purely hyperbolic nanotube Kim15 . Figure 7(a) shows the scattering efficiency of a hollow nanocylinder of thickness nm and total diameter nm, which is composed of a hyperbolic metamaterial with components of its permittivity tensor as given by Eqs. (3) and (4), considering that the filling factor . The calculations were carried out by using the Lorenz-Mie scattering method described in Refs. Chen12b ; Chen13 ; Chen15 . The minimum efficiency () is found at nm establishing the invisibility wavelength, where and . The main efficiency peak reaching is localized at nm. The field distribution in the nanotube is also plotted in Fig. 7(b) at the invisibility wavelength nm. We observe the formation of surface waves in the two air-hyperbolic metamaterial interfaces. In average, this field pattern is closer to the field distribution shown in Fig. 6(a) for an Ag inmost layer than that exhibited for a nanostructured tube of Ag ending coating depicted in Fig. 6(b). At the rear part of the nanotube, the field oscillations inside the multilayered metamaterial becomes stronger near the core; however, a surface wave at the metamaterial-core interface is only produced when a silver layer remains in contact with the core, as shown in Fig. 6(a) but not in Fig. 6(b).

## 4 Conclusions

In summary, we investigated numerically the scattering efficiency, with emphasis in the invisibility regime, of multilayered Ag-TiO nanotubes with configurations of a different ending layer. We demonstrated that when the nonconducting layer remains in contact with the core of the hollow cylinder, the characteristic invisibility of the hyperbolic nanotube is cleared away, even considering elementary layers of a few nanometers. In this case there is an additional boosting of plasmonic-Bloch resonances markedly observed on the wave fields localized at the rear part of the multilayered nanocylinder. Furthermore, the effective medium theory enabling a simplified model for the optical characterization of the nanoparticle will lead to fallacious estimations. Such effect becomes weaker for metamaterials composed of subwavelength layers with less-than-5-nanometers width.

## Acknowledgments

This work was supported by the Spanish Ministry of Economy and Competitiveness (MINECO) (TEC2014-53727-C2-1-R).

## References

- (1) U. Leonhardt, Optical conformal mapping, Science 312 (2006) 1777–1780.
- (2) J. B. Pendry, D. Schurig, D. R. Smith, Controlling electromagnetic fields, Science 312 (2006) 1780–1782.
- (3) W. Cai, U. K. Chettiar, A. V. Kildishev, V. M. Shalaev, Optical cloaking with metamaterials, Nat. Photon. 1 (2007) 224–227.
- (4) M. R. Forouzeshfard, M. H. Farzad, Electromagnetic wave propagation through two coaxial transformation-based cylindrical media, Plasmonics 10 (2015) 1345–1357.
- (5) G. X. Yu, R. Cao, M. Luo, The property of energy flow velocity in the cloak shell, Optik 126 (2015) 1990–1993.
- (6) A. Rajput, K. V. Srivastava, Design of thin simplified cloak with finite and small dynamic range constitutive tensors, Appl. Phys. A 122 (2016) 230.
- (7) A. Alu, N. Engheta, Achieving transparency with plasmonic and metamaterial coatings, Phys. Rev. E 72 (2005) 016623.
- (8) B. Edwards, A. Alu, M. G. Silveirinha, N. Engheta, Experimental verification of plasmonic cloaking at microwave frequencies with metamaterials, Phys. Rev. Lett. 103 (2009) 153901.
- (9) S. Tricarico, F. Bilotti, L. Vegni, Scattering cancellation by metamaterial cylindrical multilayers, J. Eur. Opt. Soc, Rapid Publ. 4 (2009) 09021.
- (10) X. Zhou, Y. Chen, Y. Zhou, Z. Luo, The effect of coordinate transformation function on scattering characteristics of cylindrical cloaks with a quantity of discrete layers, Opt. Commun. 364 (2016) 67–71.
- (11) A. Alu, N. Engheta, Multifrequency optical invisibility cloak with layered plasmonic shells, Phys. Rev. Lett. 100 (2008) 113901.
- (12) D. S. Filonov, A. P. Slobozhanyuk, P. A. Belov, Y. S. Kivshar, Double-shell metamaterial coatings for plasmonic cloaking, Phys. Status Solidi RRL 6 (2012) 46–48.
- (13) K.-H. Kim, Y.-S. No, S. Chang, J.-H. Choi, H.-G. Park, Invisible hyperbolic metamaterial nanotube at visible frequency, Sci. Rep. 5 (2015) 16027. doi:10.1038/srep16027.
- (14) H. E. Bussey, J. H. Richmond, Scattering by a lossy dielectric circular cylindrical multilayer, numerical values, IEEE Trans. Antennas Propag. 23 (1975) 723–725.
- (15) C. F. Bohren, D. R. Huffman, Absorption and scattering of light by small particles, Wiley, 1998.
- (16) P. Yeh, Optical waves in layered media, Wiley, New York, 1988.
- (17) L. Ferrari, C. Wu, D. Lepage, X. Zhang, Z. Liu, Hyperbolic metamaterials and their applications, Prog. Quant. Electron. 40 (2015) 1–40.
- (18) C. J. Zapata-Rodríguez, D. Pastor, V. Camps, M. T. Caballero, J. J. Miret, Three-dimensional point spread function of multilayered flat lenses and its application to extreme subwavelength resolution, J. Nanophoton. 5 (2011) 051807.
- (19) C. J. Zapata-Rodríguez, D. Pastor, M. T. Caballero, J. J. Miret, Diffraction-managed superlensing using plasmonic lattices, Opt. Commun. 285 (2012) 3358–3362.
- (20) J. J. Miret, C. J. Zapata-Rodríguez, Diffraction-free propagation of subwavelength light beams in layered media, J. Opt. Soc. Am. B 27 (7) (2010) 1435–1445.
- (21) C. J. Zapata-Rodríguez, J. J. Miret, S. Vuković, M. R. Belić, Engineered surface waves in hyperbolic metamaterials, Opt. Express 21 (2013) 19113–19127.
- (22) C. J. Zapata-Rodríguez, D. Pastor, J. J. Miret, S. Vuković, Uniaxial epsilon-near-zero metamaterials: from superlensing to double refraction, J. Nanophoton. 8 (2014) 083895.
- (23) C. Díaz-Aviñó, D. Pastor, C. J. Zapata-Rodríguez, M. Naserpour, R. Kotyński, J. J. Miret, Some considerations on the transmissivity of trirefringent metamaterials, J. Opt. Soc. Am. B 33 (2016) 116–125.
- (24) J. Elser, V. A. Podolskiy, I. Salakhutdinov, I. Avrutsky, Nonlocal effects in effective-medium response of nanolayered metamaterials, Appl. Phys. Lett. 90 (2007) 191109.
- (25) D. Torrent, J. Sánchez-Dehesa, Radial wave crystals: Radially periodic structures from anisotropic metamaterials for engineering acoustic or electromagnetic waves, Phys. Rev. Lett. 103 (2009) 064301.
- (26) H. L. Chen, L. Gao, Anomalous electromagnetic scattering from radially anisotropic nanowires, Phys. Rev. A 86 (2012) 033825.
- (27) H. L. Chen, L. Gao, Tunablity of the unconventional Fano resonances in coated nanowires with radial anisotropy, Opt. Express 21 (2013) 23619–23630.
- (28) H. L. Chen, L. Gao, Anomalous optical forces on radially anisotropic nanowires, Appl. Phys. A 121 (2015) 1053–1056.