## Abstract

We investigate the influence of metamaterials on the scaling laws of the transmission on multilayered structures composed of random sequences of ordinary dielectric and metamaterial layers. The spectrally averaged transmission in a frequency range around the fully transparent resonant mode is shown to decay with the total number of layers as 1/*N*. Such thickness dependence is faster than the 1/*N*
^{1/2} decay recently reported to take place in random sequences of ordinary dielectric slabs. The interplay of strong localization and the emergence of resonant modes within the gap leads to a non-monotonous disorder dependence of the transmission that reaches a minimum at an intermediate disorder strength.

©2008 Optical Society of America

## 1. Introduction

Metamaterials are optical structures designed to present both negative permissivity and negative permeability. This characteristic leads to unusual electromagnetic phenomena, for instance, **E**×**H** lies along the direction of -**k** for propagating plane waves. A theoretical study of metamaterials was firstly suggested a half-century ago by Veselago [1]. However, the more recent proposal of Pendry [2] for the development of such materials stimulated a wealthy of theoretical and experimental researches on this subject [3, 4, 5, 6, 7].

Novel and fascinating phenomena have been shown to emerge in structures that include metamaterials. Pendry comproved that a slab of negative refractive index material has the power to focus all Fourier components of a 2*D* image. He named this device “superlens”[8]. This structure also provides a negative lateral displacement of an obliquely incident Gaussian beam[9].

Multilayer structures have been attracting a widespread interest because of their deep implications on both fundamental and technological aspects. A periodically modulated dielectric function is the simplest way to construct a one-dimensional photonic crystal (PC) [10]. For a PC containing metamaterials, some unusual phenomena such as spurious modes with complex frequencies, discrete modes and photon tunneling modes are observed in the band structure [11]. The mixing of metamaterial and ordinary dielectric slabs creates new features such as the rising of a non-Bragg band-gap called *zero-n̄ gap* [12, 13], which emerges naturally when the volume average of the effective refractive index equals to zero. Such gap has a special characteristic: it remains unchanged when the lattice constant varies. Experimental verifications of *zero-n̄ gap* were made using *double-S shaped* metamaterials [14].

In 1*D* disordered multilayered structures, multiple incoherent superposition of waves reflected by the layers interfaces give rises to the exponential localization of most of the electromagnetic modes within the structure. A few resonant necklace modes may remain extended[15, 16, 19]. This phenomenon is similar to the electronic eigenfunction localization in disordered atomic chains and is usually termed as the Anderson localization of light. In a recent paper [18], we investigated a 1D multilayered structure composed of a random sequence of *N* ordinary dielectric slabs with refraction indices satisfying the Bragg condition. An experimental observation of resonant transport in similar structures was reported by Bertolotti *et al* [19]. We showed that the half-wavelength mode is insensitive to disorder and fully transparent. The spectrally averaged transmission in a frequency range around this mode decays as 1/*N*
^{1/2} and the localization length diverges quadratically as this resonant mode is approached. In the vicinity of the quarter-wavelength mode, the localization length diverges logarithmically and the frequency averaged transmission exhibits an stretched exponential dependence on the total thickness. These scaling laws were discussed in the context of the Anderson localization of electrons in chains with off-diagonal and correlated disorder. It was recently demonstrated theoretically that correlations may indeed strongly modify the transmission properties of electronic wires and waveguides in small energy windows [20].

New features have been also reported when metamaterials are introduced in disordered optical system. Asatryan *et. al.* [21] studied wave propagation in mixed 1*D* disordered stacks of alternating right- and left-handed layers and revealed that the introduction of metamaterials substantially suppresses the Anderson localization. At long wavelengths, the localization length in mixed stacks is several orders of magnitude larger than in normal structures.

In this work, we study the influence of metamaterials on the scaling laws for the transmission of an electromagnetic plane wave inciding normally on a multilayer random structure. We will show that, for layers do not satisfying the Bragg condition, the transmission peak around the fully transmitting mode is quite sensitive to disorder. The spectrally averaged transmission develops a faster decay as a function of the total number of layers when compared with random structures composed just of ordinary dielectric slabs. The scaling behavior of the localization length remains unchanged. The interplay between the faster finite-size scaling and the emergence of resonance necklace modes within the gap leads to a non-monotonous dependence of the transmission as a function of the disorder strength.

## 2. Model and formalism

For a monochromatic plane wave of angular frequency *ω* propagating along the *z* axis direction, the electric and magnetic fields at the interfaces of a dielectric slab of thickness *d* can be numerically obtained by a transfer matrix formalism. Assuming a normally incident wave, linearly polarized in such a way that the electric field amplitude can be written as *E⃗*(*z*)=*E*(*z*)*x̂*, the relation between the electric and magnetic field at the interface placed at *z*=*z*
_{1} and the fields at the interface placed at *z*=*z*
_{1}+*d* can be expressed in transfer matrix form as

where the phase change *δ*=±|*n*|*ωd*/*c* [the choice of sign corresponding to right-handed (RH) and left-handed (LH) slabs, respectively], |*n*| is the absolute refractive index of the medium, and
$p=\sqrt{\frac{\epsilon}{\mu}}$
, where *ε* and *µ* are the dielectric constant and magnetic permeability, respectively. For a stratified medium consisting of a sequence of *N* dielectric layers, the boundary conditions across the interface require the continuity of the parallel components of the fields. This relation can be written as a product of individual transfer matrices as

where *M _{i}* is the transfer matrix of the

*i*th layer,

*E*

_{0}and

*B*

_{0}are the electric and magnetic field amplitudes at the first interface, and

*E*and

_{N}*B*the field amplitudes at the last interface. Assuming that the incident wave is coming from the left with electric field amplitude normalized to a unitary value and also that one has only the outgoing transmitted wave at the right of the sample, the complex transmission coefficient of such sample can be written as

_{N}where *m _{ij}*’s are the elements of the total transfer matrix

*M*.

*p*and

_{i}*p*are related to the input and output media. All quantities computed here are related with the ratio between the amplitudes of the outgoing and incoming waves, which given the transmission $T\left(\omega \right)=\frac{{p}_{o}}{{p}_{i}}{\mid t\left(\omega \right)\mid}^{2}$ .

_{o}## 3. Results

For the numerical analysis, we assume non-absorbing dispersionless layers. The multilayer structure is composed of a random sequence of two distinct dielectric slabs, namely, an ordinary dielectric *A* and a left-handed dielectric *B*. The refraction indices and thicknesses were taken in such a way that the optical lengths on the slabs *A* and *B* were
$\mid {n}_{a}{d}_{a}\mid =\frac{5}{3}\mid {n}_{b}{d}_{b}\mid $
. If the slabs were chosen to satisfy the Bragg condition |*n _{a}d_{a}*|=|

*n*|, we would have the emergence of a

_{b}d_{b}*zero-n̄ gap*. This kind of gap is insensitive to disorder [

**?**]. This condition taken on a dispersionless multilayered structure creates a large photonic bang-gap that covers the whole spectrum, except for singular modes where the phase change

*δ*=

*ωnd*/

*c*=

*iπ*(

*i*=0,1,2…). This scenario does not support transmitting bands in finite-size structures. The disorderwas implemented assuming that the i-

*th*layer of the sequence has the same probability of being type

*A*or

*B*. In order to be able to explicitly show the direct influence of the metamaterial slabs, all quantities were obtained both for a mixed structure (consisting of matematerials and ordinary dielectric layers) as well as for an ordinary structure (consisting of only ordinary dielectric layers). In both cases, we assumed the same absolute values of refraction indices and thicknesses.

The first quantity analyzed was the transmission spectrum. Fig. 1 represents an average over 10^{3} realizations of disorder in a structure with *N*=200 layers. The plots are related to the ordinary structure Fig.1(a) and the mixed structure containing a metamaterial Fig.1(b). We observe the same pattern of transmission peaks irrespective to the presence or not of the metamaterial. These resonant modes occur when multiple integers of the half-wavelength equal the optical length of one of the constituent layers. Under this condition, the electromagnetic wave is not scattered by such layer, independently of the refraction index sign. We labeled the resonance frequencies as *ω _{i,j}*, where

*i*and

*j*are associated with the number of half-wavelength fitting in the optical length of each layer. The index

*i*is related to the larger optical length layer and

*j*is related to the smaller one. If

*i*and

*j*are both integers, the mode becomes fully transparent. Formally, it reflects the fact that all matrices in the product (2) turn out to be ±

*I*, where

*I*is the identity operator. Other resonant modes occur when just one of the indices is integer. In this case, only the matrices associated with the layers that fit integer numbers of half-wavelength become ±

*I*. In this case, the system behaves like a single slab composed of the other material.

An important feature is observed in the transmission spectrum around the fully transparent modes *ω*
_{0,0} and *ω*
_{5,3}. In these frequency regions, the transmission band width is much narrower for the mixed structure containing the metamaterial. In order to investigate more accurately this fact, we compute the spectral average of the transmission in the frequency range around these peaks, as shown in Fig. 2(a–b). The result represents an average over 10^{2} realizations of disorder. We compare the finite-size scaling behavior for the ordinary [Fig. 2(a)] and mixed [Fig. 2(b)] multilayer structures. The spectrally averaged transmission for the purely ordinary case decays as 1/*N*
^{1/2}. The same exponent was found in Ref. [18] where the layer parameters satisfied the Bragg condition. We found a new scaling law when one of the layers is composed of a metamaterial. In this case, the averaged transmission has a stronger dependence on the total number of layers, decaying asymptotically as 1/*N*. Another important quantity to characterize the Anderson localization of light in disordered dielectric structures is the localization length (*ξ*), considered as the inverse of the Lyapunov exponent [*ξ*=1/Λ=-lim_{N→∞}(*N*/ln*T*)]. Fig. 2(c) shows an average over 10 realizations of disorder, computed in the vicinity of *ω*
_{5,3} mode. We observe the localization length diverging quadratically as one approaches the resonance. It reproduces the same scaling law previously reported for the ordinary case [18].

The interplay of the strong finite-size decay of the transmission and the stabilization of resonance necklace modes by disorder[15, 16, 19] suggests that the influence of disorder on the transmission properties may show new features in the presence of metamaterials. In order to explore this point, we investigate the spectral average of the transmission as a function of the disorder strength. Starting with an initial periodic structure of alternating *A* and *B* layers, the disorder strength was controlled by introducing the probability *q* of the i*th* layer to be replaced by the other species. For *q*=0, one retains the periodic sequence, while the *q*=1/2 limit recovers the uncorrelated fully disordered sequence. Our results are shown in Fig. 3 for *N*=10^{2} layers, averaged over 5×10^{3} disorder configurations and within a spectral range around the *ω*
_{5,3} mode. “O” labels the ordinary case and “M” corresponds to the mixed case.

The most prominent disorder effect is to promote the Anderson localization and, consequently, to reduce the transmission band. This phenomenon is due to the multiple incoherent superposition of waves scattered at the layers interfaces. This feature is clearly observed for the purely ordinary structure (see Fig. 3) on which the spectral average of the transmission decreases monotonically as the disorder strength increases. A similar trend is observed for the mixed structure in the weak disorder limit. However, for such mixed structure, the spectrally averaged transmission presents a minimum value around *q*=0.15. This unexpected behavior is related to a secondary disorder effect, which becomes more prominent in the mixed structure. Disorder also promotes the emergence of a few resonance necklace modes that may lead to a small transmission in frequency gaps. Such secondary effect is hidden by the slow development of the gap in ordinary structures, as shown in Fig. 4(a). In the mixed structure containing a metamaterial, the gap develops rather quickly due to the stronger finite-size scaling dependence of the average transmission. This leads to the strong decay of the transmission in the weak disorder regime reported in Fig. 4(b). At weak disorder, where the Anderson localization is predominant, the band structure of the periodic case (*q*=0) gives place to a wide low-transmission frequency range. The contribution of the resonant necklace states to the transmission thus becomes predominant in the high-disorder regime. Once these modes become more frequent at strong disorder, increasing disorder promotes an enhancement of the light transmission in the spectral range around the fully transparent mode, as seen in Fig. 4(c).

Finally, we would like to stress that the narrow transmission peak in the close vicinity of the fully transparent mode is quite insensitive to disorder in the mixed structure, as can be seen in Fig.4(b–c). Actually, this peak resembles the singular frequency points of mixed structures with *zero-n̄* [12]. In finite structures, the singularity is replaced by a narrow peak whose width is inversely proportional to the total number of layers, which is consistent with our previous finding of a 1/*N* scaling law for the spectrally averaged transmission.

## 4. Summary and conclusions

In this work, we reported that the inclusion of metamaterials on the composition of a random binary multilayered structure changes the finite-size scaling law for the transmission in the frequency range around the fully transparent mode. The spectrally averaged transmission decays as 1/*N* which is faster than the scaling-law obtained for ordinary structures. In spite of this, the localization length obeys the ordinary quadratic scaling-law as the frequency approaches the resonance mode. The faster finite-size scaling and the emergence of a few transmission necklace modes inside the band-gap at the strong disorder regime promote a non-monotonous disorder dependence of the transmission in the frequency range around the fully transparent mode. This new feature is in contrast with the usual Anderson localization phenomenon that implies in a decreasing transmission as a function of the disorder strength. In structures containing metamaterials there is a characteristic disorder strength for which the system shows a minimum on the spectrally averaged transmission. Although we limited our investigation to non-absorbing structures, the main role played by absorption shall be to promote an overall reduction of the transmission in the vicinity of the resonant modes, with no direct impact on the finite-size scaling law or on the non-monotonicity of the spectrally averaged transmission. Our findings indicate that the strength of disorder shall be carefully chosen in the development of devices, such as band filters, that exploit the joint properties of binary structures, light localization and left-handed materials.

## Acknowledgments

We would like to thank the partial financial support from CAPES, CNPq, CNPq-Rede Nanobioestruturas and FINEP (Brazilian Research Agencies), as well as from FAPEAL (Alagoas State Research Agency).

## References and links

**1. **V. G. Veselago, “The electrodynamics of substances with simultaneously negative values of permissivity and permeability,” Sov. Phys. Usp. **10**, 509–514 (1968). [CrossRef]

**2. **J. B. Pendry, A. J. Holden, W. J. Stewart, and I. Youngs, “Extremely low frequency plasmons in metallic mesostructures,” Phys. Rev. Lett. **76**, 4773 (1996). [CrossRef] [PubMed]

**3. **P. Kolinko and D. R. Smith, “Numerical study of electromagnetic waves interacting with negative index materials,” Opt. Express **11**, 640 (2003). [CrossRef] [PubMed]

**4. **K. Guven, M. D. Caliskan, and E. Ozbay, “Experimental observation of left-handed transmission in a bilayer metamaterial under normal-to-plane propagation,” Opt. Express **30**, 8685 (2006). [CrossRef]

**5. **Z. Tang, R. Peng, D. Fan, S. Wen, H. Zhang, and L. Qian, “Absolute left-handed behaviors in a triangular elliptical-rod photonic crystal,” Opt. Express **13**, 9796 (2005). [CrossRef] [PubMed]

**6. **D. R. Smith and N. Kroll, “Negative refractive index in left-handed materials,” Phys. Rev. Lett. **84**, 4184 (2000). [CrossRef] [PubMed]

**7. **P. Markǒs and C. M. Soukoulis, “Transmission properties and effective electromagnetic parameters of double negative metamaterials,” Opt. Express **11**, 649 (2003). [CrossRef] [PubMed]

**8. **J. B. Pendry, “Negative refraction makes a perfect lens,” Phys. Rev. Lett. **85**, 3966 (2000). [CrossRef] [PubMed]

**9. **J. A. Kong, B.-I. Wu, and Y. Zhang, “A unique lateral displacement of a gaussian beam transmitted through a slab with a negative permittiviy and permeability,” Microwave Opt. Technol. Lett. **33**, 136–139 (2002). [CrossRef]

**10. **E. Yablonovitch, “Inhibited spontaneous emission in solid-state physics and electronics,” Phys. Rev. Lett. **58**, 2059 (1987) [CrossRef] [PubMed]

**11. **Liang Wu, Sailing He, and Linfang Shen, “Band structure for a one-dimensional photonic crystal containing left-handed materials,” Phys. Rev. B **67**, 235103 (2003) [CrossRef]

**12. **Jensen Li, Lei Zhou, C. T. Chan, and P. Sheng, “Photonic band gap from a stack of positive and negative index materials,” Phys. Rev. Lett. **90**, 083901 (2003) [CrossRef] [PubMed]

**13. **J. A. Monsoriu, R. A. Depine, M. L. Martínez-Ricci, and E. Silvestre, “Interaction between non-Bragg band gaps in 1D metamaterials photonic crystals,” Opt. Express **14**, 12958 (2006). [CrossRef] [PubMed]

**14. **Y. Yuan, L. Ran, J. Huangfu, H. Chen, L. Shen, and J. Au Kong, “Experimental verification of zero order bandgap in a layered stack of left-handed and right-handed materials,” Opt. Express **14**, 2220 (2006) [CrossRef] [PubMed]

**15. **Mher Ghulinyan, “Formation of optimal-order necklace modes in one-dimensional random photonic superlattices,” Phys. Rev. A **76**, 013822 (2007) [CrossRef]

**16. **Mher Ghulinyan, “Pediodic oscillations in transmission decay of Anderson localized one-dimensional dieletric system,” Phys. Rev. Lett. **99**, 063905 (2007) [CrossRef] [PubMed]

**17. **J. Bertolotti, S. Gottardo, D. S. Wiersma, M. Ghulinyan, and L. Pavesi, “Optical necklace states in Anderson localized 1D systems,” Phys. Rev. Lett. **94**, 113903 (2005) [CrossRef] [PubMed]

**18. **E. M. Nascimento, F. A. B. F. de Moura, and M. L. Lyra, “Scaling laws for the transmission of random binary dielectric multilayered structures,” Phys. Rev. B **76**, 115120 (2007) [CrossRef]

**19. **J. Bertolotti, M. Galli, R. Sapienza, M. Ghulinyan, S. Gottardo, L. C. Andreani, L. Pavesi, and D. S. Wiersma, “Wave transport in random systems: Multiple resonance character of necklace modes and their statistical behavior,” Phys. Rev. E **74**035602(R) (2006) [CrossRef]

**20. **M. Titov and H. Schomerus, “Nonuniversality of Anderson Localization in Short-Range Correlated Disorder,” Phys. Rev. Lett. **95**126602 (2005) [CrossRef] [PubMed]

**21. **A. A. Asatryan, L. C. Botten, M. A. Byrne, V. D. Freilikher, S. A. Gredeskul, I. V. Shadrivov, R. C. McPhedran, and Y. S. Kivsharet, “Suppression of Anderson localization in disordered metamaterials,” Phys. Rev. Lett. **99**193902 (2007) [CrossRef]