Direct observation of Anderson localization in plasmonic terahertz devices.

We present the first experimental observation of Anderson localization in the terahertz frequency range using plasmonic structures. To accomplish this goal, we designed THz waveguides consisting of a one-dimensional array of rectangular apertures that were fabricated in a freestanding metal foil. Disorder is introduced into the waveguide by offsetting the position of each aperture by a random distance within a prescribed range. For example, for a waveguide with apertures spaced by 250 μm in a periodic waveguide, 10% disorder would correspond to the apertures being shifted by a random value between ±25 μm along the waveguide axis. We find that for disorder levels below 25%, there is only an increase in the propagation loss along the device. However, for two specific waveguides with 25% disorder, we observe a spatially localized mode that lies just within the stop band of the device and exhibits a double-sided exponential spatial decay away from the maximum.

Introduction

The Anderson localization of waves is an exotic phenomenon in which a propagating wave is spatially localized purely by structural disorder, thereby inhibiting its transport across the system[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19]. Although the generality of this wave phenomenon has been proven by its demonstration in acoustic waves[5], matter waves[6] and electromagnetic waves over different frequency ranges[7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17], it has remained elusive in the terahertz domain. Here, we describe the direct observation of Anderson localization in the THz spectral range using THz plasmonic waveguides. These devices are built using periodically spaced structures fabricated into a metal foil, thus allowing surface-bound guided-wave modes below the Bragg frequency. We find that in the presence of appropriate positional disorder, a localized mode beyond the Bragg frequency becomes evident and exhibits a double-sided exponential spatial decay of the THz electric field. These results dramatically expand the spectral range, types of materials and types of applications in which localization may have an important role.

The experimental demonstration of the mesoscopic electromagnetic phenomenon of Anderson localization is known to be challenging, particularly in three-dimensional systems, wherein a critical degree of disorder is required to manifest localization[17]. In lower dimensions, the fact that all states are localized favors the experimental observation of disorder-induced localization, as long as an appropriate sample size is employed[16]. Even in low dimensions, a ‘terahertz gap’ remains in the literature on localization because dielectric materials, which are commonly employed in such demonstrations, are lossy in this spectral range[20, 21], creating a situation whereby the absorption length may be shorter than the localization length. Metals, by contrast, are characterized by high conductivity and enable the low-loss propagation of surface plasmon polaritons (SPPs)[22, 23, 24]. In this study, we employ metal foils that are periodically structured, which have been shown to support the low-loss propagation of narrowband THz-guided-wave SPP modes[25, 26, 27]. In the simplest embodiment, we have shown that rectangular apertures that fully[25] or partially[27] perforate the metal foil can guide multiple narrowband sets of frequencies. Here, the in-plane length of the aperture orthogonal to the propagation direction determines the frequency of the lowest-order propagation mode, whereas the higher-order modes can be modeled as cavity modes of the aperture, assuming an effective depth[25]. The aperture spacing along the waveguide axis determines the Bragg frequency of the periodic structure, which creates an effective band-edge for the system. Only modes that are characterized by frequencies below the Bragg frequency are observed as propagating modes.

Materials and methods

In this investigation, we examine the guided-wave properties of THz waveguides in which disorder is introduced in the position of the apertures along the propagation direction. The propagating THz electric field is characterized using THz time-domain spectroscopy, which yields both amplitude and phase information. Broadband THz radiation extending from ~0.1 to 1.5 THz is coupled to the waveguide using a single groove etched into the metal foil. The time-domain propagation properties of the coupled THz electric field can be measured at any point on or above the metal surface via electro-optic sampling by simply moving the detection crystal and the optical probe beam[28]. When <110> ZnTe is used as the detection medium, we only measure the z-component (normal to the metal surface) of the electric field, although other probe beam polarizations and crystal orientations can be used to measure the other electric field vector components[28, 29].

In Figure 1a, we show a schematic of the plasmonic waveguide, in which the apertures are periodically spaced, along with the experimental excitation and measurement geometry. Disorder is introduced into these waveguides by randomly varying the positions of each aperture by a fractional amount, σ, along the propagation direction. For example, in the periodic waveguides discussed below, the aperture spacing was set to 250 μm. Thus, a waveguide with σ=10% disorder means that each aperture would be shifted by a randomly determined distance between ±25 μm from its original position in the periodic array. A schematic diagram for a waveguide with σ=25% is shown in Figure 1b. Because Anderson localization is a stochastic process, we designed eight different waveguide samples for each of nine levels of disorder, σ, varying between 0 and 25%. All of the waveguides were fabricated via laser ablation in 750-μm-thick stainless steel foils. The measured position tolerance for the apertures was less than 5 μm. The THz time-domain signal was initially measured every 5 mm along each waveguide, with smaller measurement increments for the samples that showed evidence of localization. Detailed information about the fabrication of the waveguides and the measurement scheme is presented in the Supplementary Information.

Figure 1

Schematic diagram of the THz waveguide, including the excitation and detection scheme. (a) Broadband THz radiation is normally incident on a 2-cm-long rectangular cross-section groove, 300 μm wide by 100 μm deep. The groove is used to couple the normally incident freely propagating broadband THz radiation into SPP waves that propagate along the sample surface. The 8-cm-long waveguide consists of periodically spaced (σ=0%) rectangular apertures with the following dimensions: s=500 μm, a=150 μm, h=750 μm and d=250 μm. A <110> ZnTe crystal that can be freely positioned anywhere above the sample surface is used to measure the z-component of the THz electric field via electro-optic sampling. (b) Waveguide with identical aperture dimensions and σ=25%.

Results and discussion

Figure 2 presents our measurement of THz wave propagation through the waveguides with varying degrees of disorder. In Figure 2a, the red profile shows the transmission spectrum measured at the 8 cm terminus from the waveguide input for the periodic sample. The spectrum illustrates three dominant modes, where the lowest-order mode is determined by the fundamental resonance of the subwavelength aperture, ν=c/2s=0.3 THz, and the next two higher-order modes are related to Fabry-Perot resonances within the apertures[24]. Transmission through the waveguide arises from coupling of the THz electric field between subwavelength apertures. The manifestation of other modes associated with higher-order resonances is forbidden by the band-edge. The corresponding Bragg frequency, in our case, is given by c/2d=0.6 THz, beyond which no transmittance is seen to occur, despite the significant power beyond 0.6 THz in the input pulse shown in Supplementary Information. Figure 2a also depicts configurationally averaged transmittance spectra for σ=3% through σ=25%. In these data, the spectrum for σ=25% is averaged over six devices, and the rest are averaged over all eight devices. These additional spectral measurements demonstrate that the introduction of disorder does not alter the spectral properties of the guided-wave modes for almost all of the waveguides and, in particular, does not necessarily remove the frequency constraints imposed by the Bragg frequency.

Figure 2

Propagation properties of the waveguides. (a) Experimentally measured transmission spectra for waveguides with disorder values ranging from σ=0% to σ=25% at a distance of 2 cm from the input groove. The spectra are the averages for all eight waveguides at each disorder level, except σ=25%, where the spectra have been averaged for only six of the eight devices. The dashed vertical line corresponds to the Bragg frequency at 0.6 THz. (b) Amplitude at the peak of the lowest-order resonance (0.28 THz) along the length of the waveguide for σ=0, 10% and 25%. The data are averaged as in a. The lines show the least squares fits to the data. (c) The total decay length for the lowest-order mode as a function of σ.

For planar (unperforated) stainless steel metal foils, we have previously found that the 1/e THz electric field decay length for SPPs along the surface of the metal is 21 cm at 0.3 THz (12.6 cm at 0.5 THz and 6.3 cm at 1 THz)[24]. Perforating the metal foil to create a waveguide based on periodically spaced apertures creates an additional loss mechanism. In Figure 2b, we show the magnitude of the THz electric field measured for waveguides with periodically spaced apertures (σ=0%) at the peak of the lowest-order resonance (0.3 THz). The corresponding 1/e propagation length is 8.2 cm, which is still longer than the device length used in this study. Figure 2b also shows the electric field attenuation for two other sets of waveguides with increased disorder. The corresponding 1/e decay lengths for each value of σ are shown in Figure 2c, which appears to follow an approximately exponential dependence for reasons hitherto unknown to us. For each waveguide, we can model the amplitude at any given frequency using an exponential decay of the form , where x is the coordinate along the length of the waveguide and is the 1/e decay length. For σ=25%, for the lowest-order mode corresponds to ~34λ (3.5 cm). Assuming that , where is the decay length in a periodic waveguide (8.2 cm) and is the decay length associated with the disorder, we obtain the disorder-induced decay length =6.1 cm, which is still comparable to the sample length for σ=25%. If we consider the next two higher-order modes for the same disordered waveguides (σ=25%), is ~38λ for both resonances at 0.34 and 0.49 THz.

As noted above, we only showed the spectra and spatial decay properties for six of the eight waveguides with σ=25%. The remaining two waveguides exhibited an additional resonance just beyond the Bragg frequency in the effective bandgap region. In Figure 3a, we show the spectra for one of the two waveguides (because they exhibit very similar spectral properties) as a function of the spatial coordinate. For both of these waveguides, the spectral properties of the three propagating modes are nearly identical to the spectral properties of the other six waveguides with σ=25%, as are the spatial decay properties for the lowest-order mode. However, the decay properties for the additional mode at ~0.62 THz differ dramatically. In Figure 3b, we show the spatial properties of this new mode for both waveguides, which is clearly seen to peak in the middle of the samples away from the edges and is strongly localized over ~100 apertures. We fit the spatial decay properties of this mode in both waveguides at 0.62 THz to a double-sided exponential decay of the form , where xo is the peak position of the mode and ξ is the 1/e decay length. In both cases, we find that ξ is ~8 mm (~16λ at 0.62 THz), which is much smaller than both and .

Figure 3

Spectral and spatial properties for two waveguides with σ=25%. (a) Transmission spectra of one of the two waveguides at distances of 2, 4.5 and 8 cm from the start of the waveguide. (b) Amplitude at the peak of the mode at ~0.62 THz for both waveguides (red and black dots) measured at 2.5 mm intervals about the maximum. The lines represent the least squares fits to the data for the two waveguides.

The introduction of disorder is known to create new states inside the bandgap. For weak disorder, these states tend to be in close proximity to the band-edge, as shown in work involving photonic crystals[11] and slow light[16] waveguides. In the case of photonic crystals, it has been reported that at the band edge, the Bloch wavenumber is given by kBloch=0; we have shown a similar behavior in the case of plasmonic waveguides[27]. Using the analysis in Ref. 19, we find that the product ΔkBloch × ξ~0.83, which is less than 1 and is a signature of Anderson or strong localization. A detailed explanation of this calculation is given in the Supplementary Information. The lowest-frequency resonance inside the bandgap occurs with a resonance peak at ~0.62 THz associated with a TM [201] mode[25], which depends on only the aperture dimensions and not the aperture spacing, which is why the isolated modes in these two waveguides have the same resonance frequency even though they arise from waveguides with different spatial aperture patterns, albeit with identical levels of disorder. By further increasing the level of disorder (σ>25%), it may be possible to realize modes deeper inside the bandgap, corresponding to larger frequencies in the input spectrum.

Conclusion

In summary, we have observed clear signatures of strong localization in the THz spectral range by using THz waveguides created by perforating a metal foil and introducing positional disorder into the perforations. In contrast to work performed at optical or microwave frequencies, for which low-loss materials are readily available, the relative lack of low-loss media in the far-infrared makes such observations challenging. Thus, we only observed localization with a relatively high degree of disorder. We believe that this work will encourage further investigations of localization, both theoretically and experimentally, in regimes where losses are often inherently larger than have been encountered in earlier work (for example, in the microwave and visible-frequency ranges). In such cases, the dissipation length would be an important parameter in the analysis.