Introduction

Light trapping on metal surfaces has a long history1,2,3,4,5,6 and has been demonstrated in a wide variety of structures for applications ranging from light filtering7,8 to biosensing9 to harmonic generation10,11,12,13. Structures such as lamellar gratings, split-ring resonators, plasmonic coaxial guides and bowties, among others, have been investigated in great detail both theoretically6,14,15,16,17,18 and experimentally19,20,21,22,23,24. A special form of these light trapping nanostructures that can achieve strong broadband absorption and field enhancement is the subwavelength groove nanocavity (NC). However, structures composed of NCs resonant in the visible spectrum have a characteristic dimension6 on the order of 10 nm, which makes realizing the device with an exact resonance wavelength a fabrication challenge. In this work, we present a new method for tuning the plasmon absorption resonance wavelength with a nanometer precision anywhere from the visible (VIS) to the near infrared (NIR) in post-fabrication by deposition of a few nanometers of aluminum oxide inside the grooves via atomic layer deposition (ALD). Furthermore, in order to achieve a high field enhancement within the NCs in the NIR spectrum, the NC must have a higher aspect ratio21 with the width of around 15 nm. In such a case, the ALD-based plasmon resonance tuning becomes a powerful technique for preparing a device that would not otherwise be possible by direct fabrication. Additionally, the ALD coating provides robust protection against surface contamination and can be used to prevent oxidation in such metals as silver and copper.

Light trapping in the visible spectrum via a subwavelength array of rectangular NCs has been proposed theoretically by Le Perchec et al.6 and demonstrated experimentally in our recent work25. The schematic for such a grating is shown in figure 1(a): on resonance, light polarized perpendicular to the orientation of the NCs (denoted p-polarization) excites a standing wave surface plasmon polariton (SPP) mode within the NC, producing the electric field pattern as illustrated in figure 1(b). This SPP mode has its strongest electric field at the mouth of the groove and strongest magnetic field at the bottom. The field enhancement achieved at the mouth of the groove is advantageous for nonlinear optical applications such as harmonic generation and multiphoton photoemission25. In this way, the NCs serve both as an SPP coupling mechanism, as well as resonant cavities with strong field enhancement26 for absorption at the resonant wavelength λres.

Figure 1
figure 1

Light trapping nanocavities on a metallic surface.

A set of grooves arranged in a subwavelength array on a metal surface act collectively as an ensemble to completely absorb light at the resonant wavelength (c). The groove's width is on the order of the metal's skin depth which results in strong on-resonance light confinement in the grooves for p-polarized light (a). However, no SPP coupling occurs for s-polarization and the field profile is that of a bulk material (b).

Figure 1(c) shows the reflectivity spectrum of such a grating, measured at normal incidence; the black solid curve shows the corresponding numerical FDTD calculation. The grating was designed by numerically optimizing the period p, width w and depth h of the NCs to produce a resonance in the VIS-NIR region. In this case, for a grating made on a gold surface: λres = 720 nm, w = 13 nm, h = 42 nm, p = 600 nm. The reflectivity was measured in reference to the s-polarized light, which does not excite the SPP modes inside the NCs and, therefore, has the response of the bulk material. Two reflectivity dips are observed in the reflectivity spectrum shown in figure 1(c). The broader dip, at 720 nm, is due to the resonant SPP modes excited within the NCs, which experimentally yields a 16 times absorption increase compared to the plain gold surface. Theoretically, the reflectivity can be reduced to below 0.01% for NCs with perfectly square edges. In contrast, the classical grating-coupled SPPs are excited on the top surface27 and propagate in between the NCs. The narrow reflectivity dip at 600 nm corresponds to this classical grating-coupled SPP mode. For photovoltaic and ultrafast applications it is desirable to produce large spectral bandwidth light trapping substrates28, while for other applications, such as sensors29, narrower bandwidth is preferred. Both types of response can be achieved by adjusting the pres ratio30.

It has been previously shown that for 150 nm in diameter hole structures, a conformal dielectric coating can be used to effectively tune the transmission resonance as well as to enhance the chemical stability of the plasmonic structure38. Additionally, ALD coating is now emerging as a powerful tool to provide functionalized surfaces of increased sensitivity in bio-sensor applications39,40. In this work for the first time, we show that ALD can be used at a much smaller scale to tune 10 nm sized cavity structures in a completely reversible manner. In the case of a deep subwavelength grating, a small amount–less than a nanometer thick layer–of a dielectric on the cavity surface shifts the resonance position by more than 10 nm. The dielectric layer on the metal surface changes the boundary condition for the SPP dispersion41, shifting the resonance wavelength for the SPP standing wave excited in the NC. The dispersion relations for plasmons propagating in a Metal Insulator Metal (MIM) cavity are well known31 and show that not only does the wave vector of the plasmon travelling in a cavity depend on the metal dielectric function; it also depends on the dielectric function of the insulator that separates the two sides of the cavity. Normally this cavity is empty, so filling it with a dielectric material naturally leads to a mechanism for tuning the plasmon wave vector and hence the position of the absorption maximum. The wavelength position of the absorption maximum can be calculated directly from the normal dispersion relations of MIM cavities. This effect was successfully used in the fabrication process to test for the presence of surface contaminants on the gold grating discussed above.

Results

The sensitivity of the plasmonic grating to the presence of a dielectric coating was effectively used to tune the grating absorption resonance wavelength in post-fabrication. A well-controlled thickness of the dielectric coating was conformally deposited onto the structure by an ALD growth process, yielding a high-quality conformal coating32,33,34.

Due to the narrow opening of the NCs the alumina layer growth was non uniform inside the NCs with the middle of the NC filled with a less dense material then at the edges of the NC. The transmission electron microscopy (TEM) image showing the cross section view of the completely filled NC is shown in figure 2(a). While the top-surface alumina coating has uniform density and appears of constant shade, the alumina filling of the NC varies greatly in shade corresponding to the density variation of the alumina filling: as the dimensions of the cavity approach values less than 10 nm in width, the alumina growth is no longer uniform and instead forms a less dense material.

Figure 2
figure 2

Alumina composition variation within the nanocavity.

The TEM image of the filled NC shows a variation in the alumina density from the NC wall to the middle of the groove. This variation was used as a basis for the FDTD model: plotting the average pixel intensity–integrated over a region indicated by a rectangle in (a)–yields an exponential decay curve (b). The alumina index of refraction is therefore modeled as an exponential decay function from the nominal alumina refractive index of 1.6 to a lower value at the center of the NC (c). The decrease in the value of the refractive index is attributed to the alumina density decrease away from the cavity walls.

The optical properties of the alumina filling were modeled by approximating the non-uniformity of the alumina growth inside the NCs as an exponential decay in the index of refraction of the alumina filling corresponding to the average intensity variation inside the NC as indicated in figure 2(a) and plotted in figure 2(b). The refractive index is modeled to vary from the nominal value of nox = 1.6 at the NC edge to the reduced refractive index at the center of the NC, nc, which is treated as a fit parameter. The mathematical model for the alumina index of refraction is then given by:

where d is the distance away from the edge of the groove and t is the exponential decay length; figure 2(c) shows the plot of the modeled index of refraction inside the NCs. Two gold gratings of closely matched dimensions were prepared to measure the dielectric-based tuning: sample A with w = 13 nm, h = 42 nm, p = 600 nm and sample B with w = 15 nm, h = 42 nm, p = 600 nm. Grating A was designed for resonant absorption at around 705 nm wavelength with the reflectivity spectrum shown in figure 3(a). The structure was then coated with 20 cycles of alumina (Al2O3) ALD growth resulting in the resonance shift past 815 nm. The alumina was then removed by an hydrofluoric acid (HF) sample cleaning, which returned the grating to its original state–see figure 3(a). The same process was repeated for 50 cycles of ALD growth producing a resonance shift of 160 nm at 865 nm starting from the original resonance at 705 nm.

Figure 3
figure 3

Plasmon resonance tuning.

Depositing a layer of alumina inside the NCs strongly affects the resonance wavelength. The original response, however, can be fully recovered by etching out the alumina; (a) shows four plots: the original reflectivity spectrum, the spectrum after 20 cycles of alumina ALD compared to 50 cycles and a plot of the recovered original response after HF cleaning. This effect allows tuning of the plasmon resonance in post fabrication (b). Two gold gratings, samples A and B, were prepared: for sample A, the alumina growth was incremental, with the sample exposed to air in between each successive deposition; for sample B, the alumina was grown in one deposition from 0 nm to the target thickness, then alumina completely removed returning the resonance to the original wavelength. And the entire process repeated. Both samples are tuned from the resonance position up to ≈ 860 nm. This effect saturates once the NC is completely filled–this critical alumina thickness is marked by the vertical dashed lines in (b).

Using the process described above the resonance position was measured as a function of the alumina thickness as shown in figure 3(b). For sample A the alumina was first grown via 2, 4, or 8 ALD cycles in a deposition; then the sample was taken to air to measure the reflectivity spectrum. A new layer of alumina was then incrementally grown. Blue circles in figure 3(a) correspond to each measured resonance position; the solid line is the FDTD model fit based on the dielectric constant as described above (the steps in the theory curve are due to the finite size of the mesh used). When the sample was taken out to air, a reacted layer would form that resulted in the overall alumina growth to be less dense yielding nc = 1.2 and t = 2.5 nm. The growth rate on the top surface was measured to be 1.1 Å/cycle based on the TEM picture shown in figure 2(a). The growth inside the NC is modeled to be twice that amount (2.2 Å/cycle) as discussed above in good agreement with the TEM cross section showing the light region of the alumina in the center of the NC.

For sample B, the alumina was grown to the target thickness all in one deposition starting from a completely cleaned surface each time. After measuring the resonance position at the target alumina thickness, the sample was cleaned using a combination of the HF, acetone and isopropyl alcohol (IPA) and the resonance position of the cleaned sample was consistently confirmed to be before growing the new target thickness of the alumina layer (no appreciable change in the surface quality of the nanopattern was evident in the SEM images). This way, a total target thickness of alumina was deposited onto the clean structure without any contamination. This process was repeated nine times and the result is plotted in figure 3(a) (black squares) with the corresponding FDTD model. Alumina grown this way yielded a more dense coating with the FDTD model fit parameters of nc = 1.4 and t = 4 nm. Unlike sample A, sample B was not exposed to air in between the successive layers, therefore the alumina growth rate is slightly higher, 1.25 Å/cycle on the top surface, or 2.5 Å/cycle inside the NC, which is consistent with growth rates reported for alumina ALD at room temperate under similar conditions34,35,36.

Both samples exhibit the same tuning behavior, in good agreement with the FDTD: in the region where the NC is not completely filled, sub nanometer size increase in the thickness of the dielectric layer partially filling the NC produces a large shift in the resonance position. In this regime, the resonance can be tuned anywhere up to the NIR region. Once the NC is filled, the tuning saturates. The thickness of the alumina layer at this point is half the NC width and this critical thickness is marked with vertical dashed lines in figure 3(a). Any continued growth is equivalent to a small decrease in the periodicity of the subwavelength grating yielding only a minor further shift in the resonance wavelength.

Discussion

ALD deposition of a conformal dielectric coating inside the metallic NCs allows a wide spectral bandwidth plasmon resonance tuning in a completely reversible manner. This process can be repeated many times–reverting to the original grating response by HF cleaning–in order to achieve the optimal response. This flexibility offered by the dielectricbased tuning is also useful beyond post-fabrication resonance adjustment. Preparing a plasmonic grating resonant in the NIR with large field enhancement requires extreme aspect ratio, which is a major fabrication challenge37. However, NIR resonance can be achieved by designing a grating for resonance towards the red part of the visible spectrum and later converting it into the NIR grating by controllable growth of a thin dielectric coating. This technique is completely reversible: the alumina is readily removed from the gold grating by an HF etch completely recovering the initial optical response. This effect also opens the possibility of using such gratings as sensors for probing sub-nanometer thick layers of a dielectric.

In summary, we have presented a design and post-fabrication tuning of a subwavelength plasmonic metallic grating for complete control over the resonance wavelength, the absorption strength, as well as spectral and angular bandwidths. Exquisite control over the resonance position was demonstrated in this study of gold gratings that were tuned in a completely reversible manner from a wavelength of 720 nm up to 860 nm by ALD growth of a few nanometer thick layer of alumina. This coating also serves as a protective layer for the metal structure, allowing fabrication of highly controlled resonant light trapping metal surface with resistance to environmental contaminants. Both types exhibit optical characteristics in good agreement with the FDTD model presented. Combining the design of the grating with the tunability offered by the dielectric coating allows production of a robust, completely absorbing substrate with a tunable resonance wavelength.

Methods

The plasmonic grating fabrication is based on the evaporation of gold onto a template prepared by electron beam lithography (see ref. 37 for detailed fabrication description). The evaporated metal is then attached to a support plate and peeled-off from the template, with the residual electron beam resist etched out using HF. A properly cleaned structure–rinsed in acetone and IPA–had a resonance at λres = 705 nm. If, however, a thin residue layer formed during the final cleaning, the resonance position could shift by up to 15 nm to λres = 720 nm. Due to the nanometer dimensions of the NCs, even a few atomic layers of a dielectric result in an appreciable shift in the resonance wavelength.

For this study we employed a room temperature alumina plasma activated ALD deposition (25 C to avoid deterioration of the NC shape on the gold surface) using an Oxford Instruments FlexAL reactor. The process consisted of the following steps: (1) 30 ms exposure to trimethylaluminium (TMA) dose at 15 mTorr without any oxygen flow; (2) followed by 1 s TMA purge and (3) oxygen inflow at 60 sccm for another 1.5 s; (4) with the final 2 s 400 W oxygen plasma oxidation at 15 mTorr with oxygen flow maintained at 60 sccm. Throughout the entire process argon gas flow was set at 100 sccm.