In Situ Planarization of Huygens Metasurfaces by Nanoscale Local

Aug 29, 2017 - We have utilized this approach for the design, fabrication and characterization of amorphous silicon based all-dielectric Huygens metas...
3 downloads 0 Views 8MB Size
Article pubs.acs.org/journal/apchd5

In Situ Planarization of Huygens Metasurfaces by Nanoscale Local Oxidation of Silicon Jonathan Bar-David, Noa Mazurski, and Uriel Levy* Department of Applied Physics, the Benin School of Engineering and Computer Science, the Center for Nanoscienceand Nanotechnology, The Hebrew University of Jerusalem, Jerusalem, 91904, Israel S Supporting Information *

ABSTRACT: Metasurfaces are becoming a flourishing field of research, with diverse applications, such as planar optical components and structural colors. While metallic metasurfaces are typically few tens of nanometers in their thickness, their dielectric counterparts typically span few hundreds of nanometers in thickness variations. This makes the stacking of multilayers a bit challenging. To mitigate this challenge, we have developed a new approach for the realization of dielectric metasurfaces. Our approach is based on the nanoscale local oxidation of silicon (LOCOS), allowing to achieve planar metasurface structures. We have utilized this approach for the design, fabrication and characterization of amorphous silicon based all-dielectric Huygens metasurfaces. These metasurfaces show clear electric and magnetic resonances, which can be structurally tuned. The obtained results are in good agreement with numerical simulations taking into account the unique shape of the nanoantennas. Relying on a robust approach for their realization, and combined with the important feature of in situ planarization, we believe that such planarized metasurfaces will become a viable technology for future applications. KEYWORDS: Huygens metasurface, dielectric metasurface, LOCOS

O

as nanoantennas for a large variety of photonic applications. by properly designing nanoantennas, it is now possible to control the amplitude, phase, and polarization of the transmitted and reflected fields and by fabricating arrays of resonant nanoantennas one may fabricate resonant metasurfaces with tailored optical properties for different applications including, for example, color engineering,11,18,19,21,28,33,34 nonlinear optics,35−38 reflectionless surfaces and perfect absorbers,1,8,39−41 metalenses and other optical elements,10,24,42−44 to name a few. The common to all these demonstrations is the ability to control the polarization, phase, amplitude and frequency of transmitted and reflected light, as well as the near-field of the metasurface, to design arbitrary optical properties. Many metasurfaces use metals and plasmonic phenomena to achieve such effects.17,18,28,29,45,45−51 Yet, metallic metasurfaces suffer significant losses. In contrast, dielectric metasurfaces show only little absorption in most cases and therefore offer high transmission and reflection efficiency reaching the limit of almost no absorption losses and pure control of phase and polarization.1,3,6−9,11−14,39,40,43,52 On the other hand, while metallic metasurface are nearly flat, a typical dielectric metasurface consists of nanoantenna arrays which are situated on the surface of a flat substrate, making the structure

ptical metasurfaces are nanopatterned surfaces exhibiting desired optical characteristics depending on the nano patterns and the structures which are defined on the surface.1−21 This is in contrast to conventional optical elements that rely on the materials’ bulk properties. One may distinguish between two different types of metasurfaces, being resonant or nonresonant. While the latter relies on concepts such as phase accumulation, effective index and geometrical phase to control light,2,3,5−7,9,10,22−27 the optical properties of resonant metasurfaces originate at the optical resonant behavior of the individual scattering elements, which are sometimes defined as optical nanoantennas.1,11,13−18,28,29 The scattering, transmission, and absorption of electromagnetic waves by small particles is an old topic first considered by Mie,30−32 who calculated the inner and outer fields for scattering of plane waves by spheres, concluding they are the sum of the basic multipole fields.31,32 The multipole fields are therefore the electromagnetic eigenmodes of the fields in and around small particles.31,32 These modes and their resonant frequencies are determined mainly by particle size, geometry, aspect ratio, refractive index, and the index contrast between the particle and surrounding medium. Following the rapid progress in nanofabrication techniques, these observations pioneered by Mie are now of great significance to the photonic research, as nowadays researchers are attempting to harness the optical response of small particles and cavities such © 2017 American Chemical Society

Received: June 28, 2017 Published: August 29, 2017 2359

DOI: 10.1021/acsphotonics.7b00688 ACS Photonics 2017, 4, 2359−2366

ACS Photonics

Article

Figure 1. LOCOS fabrication scheme. (A) a-Si, Si3N4 and e-beam resist layers are deposited on a Quartz substrate. (B) ZEP etch mask is defined by electron beam lithography. (C) RIE etch is performed through the Si3N4 and some of the a-Si layer, leaving a thin a-Si layer close to substrate. (D) Silicon oxidation followed by removal of Si3N4 layer. The end result is a quasi-planar metasurface.

Figure 2. SEM micrographs of aSi LOCOS nanoantennas. (Left) cylindrical nanoantennas with radius of ∼130 nm and gap of ∼150 nm. (Right) rectangular nanoantennas with length of ∼420 nm, width of ∼160 nm, and gaps of ∼220 nm. In both antenna types the a-Si thickness is ∼270 nm. The a-Si antennas are seen as dark gray areas, and the lighter “frame” around them is the oxide layer.

in an oxygen reach furnace, with the end result of having silicon dioxide replacing the silicon in desired location across the device. More recently, this approach has been used for the purpose of defining smooth, low-loss optical waveguides on silicon substrate.53−56 In this technique, a protective layer is deposited above a silicon (Si) substrate. Next, photo- or electron beam-lithography patterns are transferred to the protective layer by etching. Finally, these patterns are transferred to the silicon substrate by oxidation, rather than by further etch steps. The oxidation process allows for smooth structures to be formed, significantly reducing roughness and consequently undesired scattering of light. The original approach results in a highly nonplanar structure due to the higher volume of the silicon dioxide layer that is formed instead of the silicon. However, Recent publications also show that by partial etching of the Si substrate prior to oxidation, one may

nonplanar. As a result, special efforts needs to be carried out in order to planarize the surface if the stacking of additional layers is needed. Furthermore, the need for etching of the dielectric metasurface structure leads to unavoidable surface roughness and consequently scattering loss. In this work, we demonstrate a novel approach which allows to fabricate smooth and planar dielectric metasurfaces with optical resonances at the telecom regime. By doing so, we address some of the major issues of current metasurface technology and provide a path for even brighter future for metasurfaces.



FABRICATION Our approach relies on the local oxidation of silicon (LOCOS) fabrication technique. LOCOS is a fabrication method originally used for the isolation of adjacent MOS transistors from each other. In short, it is based on the oxidation of silicon 2360

DOI: 10.1021/acsphotonics.7b00688 ACS Photonics 2017, 4, 2359−2366

ACS Photonics

Article

Figure 3. AFM measurement. (Left) AFM scan of nanoantenna array (same as in Figure 2, left). Blue line marks the cross section. (Right) AFM cross-section reveals a ±25 nm height variation due to oxidation swelling.

obtain quasi-planar structures without jeopardizing the smoothness of the obtained structures.55 An illustration of our fabrication process is shown in Figure 1. Unlike previous demonstrations, in this work, we chose to utilize the LOCOS technique to fabricate quasi-planar nanoantenna arrays in an amorphous silicon (a-Si) layer. The fabrication process begins by depositing a 275 nm thick hydrogenated a-Si layer on a quartz substrate by PECVD (Oxford). Atop this layer, we deposited the oxidation-shield layer, a 160 nm thick silicon nitride (Si3N4) by PECVD. Sample was then spin coated with electron beam (e-beam) resist (ZEP 520), and nanoantenna patterns were created by e-beam lithography (Raith eline). Note that ZEP is a positive resist, hence, e-beam patterns are the interantenna gaps. After e-beam lithography, sample was developed and the pattern was transferred to the silicon nitride layer by a reactive ion etch (RIE, Corial), where the ZEP e-beam resist is used as the etch mask. In this step we etched the full height of the silicon nitride layer and exposed the a-Si layer in the bottom of the created trenches.Next, a second etch step was performed, partially etching the a-Si layer, to compensate for oxidation swelling. Finally, the e-beam resist was removed and the sample was oxidized by heating the sample to ∼1000 °C in an H2O rich atmosphere. During oxidation, the silicon nitride layer acts as an oxidation shield defining the oxidized regions. Finally, the silicon nitride layer is removed, and the planarized structure of amorphous-silicon nanoantennas, embedded in SiO2 is formed. A scanning electron microscope (SEM) image showing the top view of some of the fabricated samples is shown in Figure 2. Atomic force microscope (AFM) scans of the fabricated samples are shown in Figure 3. An SEM cross section revealing the detailed subterranean structure of the nanoantennas is shown in Figure 4. These scans clearly show the oxidized ring surrounding each of the Si nanoantennas. The surface height variation is of the order of ±25 nm. The relation between antenna layer thickness and final surface variation was not investigated directly. However, based on current results we assume that the antenna thickness does not have a significant effect on the final flatness, which is governed mostly by the accuracy of the etch and the oxidation steps. We have chosen to fabricate arrays of identical antennas with radii and period varying between arrays in order to better understand and control the fabrication process and the optical properties of the antenna structures.

Figure 4. Scanning electron microscope (SEM) micrograph image (colored for clarity) showing a cross section in fabricated nanoantenna. The nanoantennas are seen as rounded-side trapezoids in middle part of the image. Due to the resonant nature of these nanoantennas, the exact fabricated shape has a major influence on the transmission spectrum of these structures. Interface between the amorphous silicon layer and quartz substrate is clearly seen.

The cross section shown in Figure 4 further reveals the antenna sidewall geometry which is curved inward. This geometry is an outcome of the oxidation process, the curved structure is a result of oxygen diffusion into the silicon and the in- and outward swelling of the formed silicon dioxide. This unique structure was considered in the numerical simulations, as will be discussed later on.



NUMERICAL SIMULATIONS The optical transmission spectrum of the fabricated structure was simulated by 3D FDTD software (Lumerical FDTD solutions). The transmission of nanoantenna arrays was considered by applying periodic boundary conditions and plane-wave illumination. Due to the resonant nature of these nanoantennas, even a slight change of few nanometers in the simulated antenna geometry, for example, its height, radius, or sidewall profile, results in a noticeable variation in optical performance and shift in the resonance frequency. Figure 5 presents the results of such simulations, calculating the transmitted amplitude for cylindrical nanoantennas with height of 275 nm and varying radii of 60 to 400 nm, for a broad SWIR spectrum with wavelength ranging from 600 to 1600 nm. The gap was kept constant at 200 nm, that is, the period is 2361

DOI: 10.1021/acsphotonics.7b00688 ACS Photonics 2017, 4, 2359−2366

ACS Photonics

Article

Figure 5. FDTD simulation of antenna transmission. (A) Transmission spectrum of an antenna array for 600−1600 nm. The MD-ED mode splitting is clearly seen at a wavelength of ∼1200 nm for antenna radii of 250−300 nm. Diffraction effects are seen in the upper left corner, for short wavelengths and (relatively) large antennas, and the Kerker condition,58 a decrease and narrowing of the resonant reflection, is met by antennas of radius ∼150 nm at a wavelength of about 900 nm. The dashed line marks the antenna radius for which we present the spectra in subfigure (B). (B) Transmission spectrum of a single antenna array with antenna radius of 330 nm. MD, “magnetic dipole” mode; ED, “electric dipole” mode. The ED mode is at longer wavelength due to antenna shape and aspect ratio. The curved sidewall has been approximated by a parabolic shape with a curvature of about 30 nm inward.

circular vortex-shaped electric field in the xz plane, as depicted in Figure 6 (top row), while the ED modal fields at the wavelength of 1446 nm have a strong electric field at the x̂ direction and a vortex-shaped magnetic field loop at the yz plane, as depicted in Figure 6 (bottom row). It may also be noticed that the ED modal fields are actually concentrated close to the nanoantenna boundary due to the asymmetric antenna structure and (relatively) low index contrast between the antenna and the surrounding medium (here, the substrate), which increase the effective volume of the nanoantenna.

increasing with the radius. As may be seen in Figure 5A, the transmission properties of these nanoantennas change with the antenna radius. The change is mostly observable as a wavelength shift and separation of the two distinct transmission minima. These two resonances can also be observed in Figure 5b, which is a line scan of Figure 5a at a constant cylinder radius of 329 nm. Indeed, the two distinct minima, labeled as “electric dipole” (ED) and “magnetic dipole” (MD) resonances are seen, as also reported in, for example, refs 8, 31, 40, and 57. These resonances are actually the two first Mie resonances, which “switch order”, that is, the ED mode is found at a longer wavelength (or lower frequency) than the MD mode, due to the small aspect ratio of the nanoantennas, with a height of 275 nm and a radii of 330 nm. Figure 6 shows the calculated electric and magnetic fields for the two modes indicated in Figure 5B. It may be seen that under illumination with a x-polarized plane wave, the MD modal field at a wavelength of 1306 nm has the expected strong magnetic field in the ŷ direction and the



EXPERIMENTAL CHARACTERIZATION Following the design and fabrication of the LOCOS based metasurface, we have experimentally characterized its optical response. Transmission measurements were done by illuminating the sample with a white light source (tungsten-halogen lamp) through a microscope condenser lens. The transmitted light is collected by an objective lens (Nikon, 50×, NA 0.45) 2362

DOI: 10.1021/acsphotonics.7b00688 ACS Photonics 2017, 4, 2359−2366

ACS Photonics

Article

Figure 6. Electric and magnetic fields at resonance, for antenna with a 330 nm radius. Light propagates toward the antenna along the z axis. The black lines denote the antenna cross-section. As described above, sidewall curvature is approximated by a parabola. Black arrows are added as guides to the eye to indicate the direction of electric and magnetic fields. (Top row) The magnetic dipole mode at λ = 1306 nm exhibits an “electric field loop” at xz plane and strong magnetic field at the y direction. (Bottom row) the electric dipole mode at λ = 1446 nm exhibits a “magnetic field loop” at the yz plane and a strong electric field at the x direction.

Figure 7. Measured transmission spectrum for antenna with radii of 70−300 nm. With the reduction of aspect ratio, the resonant modes begin to separate, and for nanoantennas with radius >250 nm the modes are easy to distinguish. On the other hand, for smaller antenna radius of about 150 nm, one may observe spectral overlap between the MD and the ED modes, followed by enhanced transmission around 800 nm wavelength. This is known as the Kerker condition. The transmission is limited by the band to band absorption in the a-Si layer.

ing radius. This comparison reveals a good agreement between the two. Furthermore, the gradual spectral separation between the electric and magnetic modes can be clearly observed. In addition to the characterization of cylindrical nanoantennas, we have also fabricated and characterized rectangular nanoantenna arrays (such as depicted in Figure 2). The transmission through such arrays was measured by the same method described above, with light polarized along the nanoantennas’ long axis. Due to the geometry of these nanoantennas, with length of 300−400 nm and width of ∼150 nm, the resonant frequencies of these antennas are closer

and directed into a spectrometer (Horiba, microHR\Ocean Optics, Flame). All measurements are later normalized to the transmission through the system. Figure 7 presents the measured transmission of our fabricated cylindrical nanoantenna arrays. As may be seen, the main features predicted by simulation are clearly seen in our measurements, namely the two resonant modes and their overlap for nanoantennas of radius about 150 nm. To further compare our measured results to the simulated prediction, we present in Figure 8 a comparison of the numerical simulations (left) to the measurements (right) of nanoantennas with gradually increas2363

DOI: 10.1021/acsphotonics.7b00688 ACS Photonics 2017, 4, 2359−2366

ACS Photonics

Article

Figure 8. Transmission spectrum of cylindrical nanoantennas with radii of ∼90 to ∼250 nm. (Left) full wave numerical simulation, (Right) measured results. The shifting of modes toward longer wavelengths and the splitting of the MD and ED modes are clearly seen in both simulation and measured data. The Kerker condition, in which MD and ED resonances overlap is denoted by the red box. Full transmission is not obtained as the wavelength is below the material band gap.

Figure 9. Transmission spectrum of rectangular nanoantennas, with light linearly polarized along the antennas’ long axis: (left) simulation, (right) measurement. Notice the gradual change in resonant frequencies as the antenna length is increased from 300 nm at the bottom row to 400 nm at the top. Spectrum complexity arises primarily from the intrinsic absorption of a-Si for frequencies above the material band gap. The curved sidewall has been approximated by a segment shape with a curvature of about 30 nm inward.

2364

DOI: 10.1021/acsphotonics.7b00688 ACS Photonics 2017, 4, 2359−2366

ACS Photonics



to the visible spectrum (compared to the cylindrical nanoantennas shown previously). As such, the demonstration of LOCOS based nanoantennas allows for an additional degree of freedom in mitigating specific applications toward the visible band, while maintaining dimensions which are not too difficult to fabricate. Figure 9 shows the simulated (left) and the measured (right) transmission of nanoantenna arrays with length varying from 300 to 400 nm. Here, the transmission spectra are more complex due to the intrinsic silicon absorption at these frequencies resulting from interband transition. Nevertheless, we are able to identify the major trends in both simulation and measured data, namely, the general spectral location of the antenna resonances and their gradual shift toward longer wavelengths as the antenna length increases.

SUMMARY We have demonstrated the viability of the LOCOS technique for the fabrication of planar, optically resonant dielectric metasurfaces in a-Si layers. Measurements of our fabricated metasurfaces show good agreement with simulation predictions, indicating that these metasurfaces can operate in the visible and NIR spectral regimes, above and below the Si band gap frequency. We have designed and fabricated metasurfaces based on nanoantennas of different sizes and geometries and by this we have shown the fabrication robustness provided by this approach, as well as it is compatibility with current CMOS techniques. Following the results reported here, one may use this technique for the implementation of planar metasurface devices, for example, metalenses, metafilters, metasplitters, metaholograms, and more. Finally, the LOCOS technique and the use of PECVD deposited materials on commercially available substrates pave the way for future implementation in optical devices and in few-layer stack structures. ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsphotonics.7b00688. Figure S1 (PDF).



REFERENCES

(1) Kuznetsov, A. I.; Miroshnichenko, A. E.; Brongersma, M. L.; Kivshar, Y. S.; Luk’yanchuk, B. Optically Resonant Dielectric Nanostructures. Science 2016, 354 (6314), aag2472. (2) Bomzon, Z.; Biener, G.; Kleiner, V.; Hasman, E. Space-Variant Pancharatnam−Berry Phase Optical Elements with ComputerGenerated Subwavelength Gratings. Opt. Lett. 2002, 27 (13), 1141− 1143. (3) Desiatov, B.; Mazurski, N.; Fainman, Y.; Levy, U. Polarization Selective Beam Shaping Using Nanoscale Dielectric Metasurfaces. Opt. Express 2015, 23 (17), 22611−22618. (4) Yariv, A. Photonics Optical Electronics in Modern Communications, 6th ed.; Oxford University Press: New York; Oxford, 2007. (5) Lin, D.; Fan, P.; Hasman, E.; Brongersma, M. L. Dielectric Gradient Metasurface Optical Elements. Science 2014, 345 (6194), 298−302. (6) Bar-David, J.; Stern, L.; Levy, U. Dynamic Control over the Optical Transmission of Nanoscale Dielectric Metasurface by Alkali Vapors. Nano Lett. 2017, 17 (2), 1127−1131. (7) Levy, U.; Kim, H.-C.; Tsai, C.-H.; Fainman, Y. Near-Infrared Demonstration of Computer-Generated Holograms Implemented by Using Subwavelength Gratings with Space-Variant Orientation. Opt. Lett. 2005, 30 (16), 2089−2091. (8) Iyer, P. P.; Butakov, N. A.; Schuller, J. A. Reconfigurable Semiconductor Phased-Array Metasurfaces. ACS Photonics 2015.2107710.1021/acsphotonics.5b00132 (9) Arbabi, A.; Horie, Y.; Bagheri, M.; Faraon, A. Dielectric Metasurfaces for Complete Control of Phase and Polarization with Subwavelength Spatial Resolution and High Transmission. Nat. Nanotechnol. 2015, 10 (11), 937−943. (10) Yu, N.; Capasso, F. Flat Optics with Designer Metasurfaces. Nat. Mater. 2014, 13 (2), 139−150. (11) Proust, J.; Bedu, F.; Gallas, B.; Ozerov, I.; Bonod, N. AllDielectric Colored Metasurfaces with Silicon Mie Resonators. ACS Nano 2016, 10 (8), 7761−7767. (12) Jahani, S.; Jacob, Z. All-Dielectric Metamaterials. Nat. Nanotechnol. 2016, 11 (1), 23−36. (13) Evlyukhin, A. B.; Novikov, S. M.; Zywietz, U.; Eriksen, R. L.; Reinhardt, C.; Bozhevolnyi, S. I.; Chichkov, B. N. Demonstration of Magnetic Dipole Resonances of Dielectric Nanospheres in the Visible Region. Nano Lett. 2012, 12 (7), 3749−3755. (14) Fu, Y. H.; Kuznetsov, A. I.; Miroshnichenko, A. E.; Yu, Y. F.; Luk’yanchuk, B. Directional Visible Light Scattering by Silicon Nanoparticles. Nat. Commun. 2013, 4, 1527. (15) Kuznetsov, A. I.; Miroshnichenko, A. E.; Fu, Y. H.; Zhang, J.; Luk’yanchuk, B. Magnetic Light. Sci. Rep. 2012, 2, 492. (16) Glybovski, S. B.; Tretyakov, S. A.; Belov, P. A.; Kivshar, Y. S.; Simovski, C. R. Metasurfaces: From Microwaves to Visible. Phys. Rep. 2016, 634, 1−72. (17) Guo, R.; Rusak, E.; Staude, I.; Dominguez, J.; Decker, M.; Rockstuhl, C.; Brener, I.; Neshev, D. N.; Kivshar, Y. S. Multipolar Coupling in Hybrid Metal−Dielectric Metasurfaces. ACS Photonics 2016, 3 (3), 349−353. (18) Hedayati, M. K.; Elbahri, M. Review of Metasurface Plasmonic Structural Color. Plasmonics 2016, 1−17. (19) Kristensen, A.; Yang, J. K. W.; Bozhevolnyi, S. I.; Link, S.; Nordlander, P.; Halas, N. J.; Mortensen, N. A. Plasmonic Colour Generation. Nat. Rev. Mater. 2016, 2, 16088. (20) High, A. A.; Devlin, R. C.; Dibos, A.; Polking, M.; Wild, D. S.; Perczel, J.; de Leon, N. P.; Lukin, M. D.; Park, H. Visible-Frequency Hyperbolic Metasurface. Nature 2015, 522 (7555), 192−196. (21) Huang, Y.-W.; Chen, W. T.; Tsai, W.-Y.; Wu, P. C.; Wang, C.M.; Sun, G.; Tsai, D. P. Aluminum Plasmonic Multicolor MetaHologram. Nano Lett. 2015, 15 (5), 3122−3127. (22) Devlin, R. C.; Ambrosio, A.; Wintz, D.; Oscurato, S. L.; Zhu, A. Y.; Khorasaninejad, M.; Oh, J.; Maddalena, P.; Capasso, F. Spin-toOrbital Angular Momentum Conversion in Dielectric Metasurfaces. Opt. Express 2017, 25, 4239.





Article

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Jonathan Bar-David: 0000-0002-4464-636X Uriel Levy: 0000-0002-5918-1876 Author Contributions

The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. Funding

This research was supported by the ActiPlAnt program of Einstein foundation, Berlin, Germany, and the three-dimensional printing project from the Israeli Ministry of Science. Notes

The authors declare no competing financial interest. 2365

DOI: 10.1021/acsphotonics.7b00688 ACS Photonics 2017, 4, 2359−2366

ACS Photonics

Article

(23) Holloway, C. L.; Kuester, E. F.; Gordon, J. A.; O’Hara, J.; Booth, J.; Smith, D. R. An Overview of the Theory and Applications of Metasurfaces: The Two-Dimensional Equivalents of Metamaterials. IEEE Antennas Propag. Mag. 2012, 54 (2), 10−35. (24) Lalanne, P.; Chavel, P. Metalenses at Visible Wavelengths: Past, Present, Perspectives. Laser Photonics Reviews 2017, 11, 1600295. (25) Lalanne, P.; Hazart, J.; Chavel, P.; Cambril, E.; Launois, H. A Transmission Polarizing Beam Splitter Grating. J. Opt. A: Pure Appl. Opt. 1999, 1 (2), 215. (26) Lalanne, P.; Hugonin, J. P.; Chavel, P. Optical Properties of Deep Lamellar Gratings: A Coupled Bloch-Mode Insight. J. Lightwave Technol. 2006, 24 (6), 2442. (27) Lalanne, P. Waveguiding in Blazed-Binary Diffractive Elements. J. Opt. Soc. Am. A 1999, 16 (10), 2517−2520. (28) Clausen, J. S.; Højlund-Nielsen, E.; Christiansen, A. B.; Yazdi, S.; Grajower, M.; Taha, H.; Levy, U.; Kristensen, A.; Mortensen, N. A. Plasmonic Metasurfaces for Coloration of Plastic Consumer Products. Nano Lett. 2014, 14 (8), 4499−4504. (29) Huang, L.; Chen, X.; Mühlenbernd, H.; Zhang, H.; Chen, S.; Bai, B.; Tan, Q.; Jin, G.; Cheah, K.-W.; Qiu, C.-W.; Li, J.; Zentgraf, T.; Zhang, S. Three-Dimensional Optical Holography Using a Plasmonic Metasurface. Nat. Commun. 2013, 4, 2808. (30) Mie, G. Beiträge Zur Optik Trüber Medien, Speziell Kolloidaler Metallösungen. Ann. Phys. 1908, 330 (3), 377−445. (31) Bohren, C. F.; Huffman, D. R. Absorption and Scattering of Light by Small Particles; John Wiley and Sons, 2008. (32) Jackson, J. D. Classical Electrodynamics, 3rd ed.; Wiley: New York, NY, 1999. (33) Mirshafieyan, S. S.; Guo, J. Silicon Colors: Spectral Selective Perfect Light Absorption in Single Layer Silicon Films on Aluminum Surface and Its Thermal Tunability. Opt. Express 2014, 22 (25), 31545−31554. (34) Kumar, K.; Duan, H.; Hegde, R. S.; Koh, S. C. W.; Wei, J. N.; Yang, J. K. W. Printing Colour at the Optical Diffraction Limit. Nat. Nanotechnol. 2012, 7 (9), 557−561. (35) Salomon, A.; Zielinski, M.; Kolkowski, R.; Zyss, J.; Prior, Y.; Size. and Shape Resonances in Second Harmonic Generation from Silver Nanocavities. J. Phys. Chem. C 2013, 117, 22377. (36) Almeida, E.; Bitton, O.; Prior, Y. Nonlinear Metamaterials for Holography. Nat. Commun. 2016, 7, 12533. (37) Cai, W.; Vasudev, A. P.; Brongersma, M. L. Electrically Controlled Nonlinear Generation of Light with Plasmonics. Science 2011, 333 (6050), 1720−1723. (38) Segal, N.; Keren-Zur, S.; Hendler, N.; Ellenbogen, T. Controlling Light with Metamaterial-Based Nonlinear Photonic Crystals. Nat. Photonics 2015, 9 (3), 180−184. (39) Staude, I.; Miroshnichenko, A. E.; Decker, M.; Fofang, N. T.; Liu, S.; Gonzales, E.; Dominguez, J.; Luk, T. S.; Neshev, D. N.; Brener, I.; Kivshar, Y. Tailoring Directional Scattering through Magnetic and Electric Resonances in Subwavelength Silicon Nanodisks. ACS Nano 2013, 7 (9), 7824−7832. (40) Decker, M.; Staude, I.; Falkner, M.; Dominguez, J.; Neshev, D. N.; Brener, I.; Pertsch, T.; Kivshar, Y. S. High-Efficiency Dielectric Huygens’ Surfaces. Adv. Opt. Mater. 2015, 3 (6), 813−820. (41) Yu, Y. F.; Zhu, A. Y.; Paniagua-Domínguez, R.; Fu, Y. H.; Luk’yanchuk, B.; Kuznetsov, A. I. High-Transmission Dielectric Metasurface with 2π Phase Control at Visible Wavelengths. Laser Photonics Rev. 2015, 9 (4), 412−418. (42) Yifat, Y.; Eitan, M.; Iluz, Z.; Hanein, Y.; Boag, A.; Scheuer, J. Highly Efficient and Broadband Wide-Angle Holography Using PatchDipole Nanoantenna Reflectarrays. Nano Lett. 2014, 14, 2485. (43) Yifat, Y.; Iluz, Z.; Eitan, M.; Friedler, I.; Hanein, Y.; Boag, A.; Scheuer, J. Quantifying the Radiation Efficiency of Nano Antennas. Appl. Phys. Lett. 2012, 100 (11), 111113. (44) Khorasaninejad, M.; Chen, W. T.; Devlin, R. C.; Oh, J.; Zhu, A. Y.; Capasso, F. Planar Lenses at Visible Wavelengths. ArXiv160502248 Phys. 2016, na.

(45) Almeida, E.; Shalem, G.; Prior, Y. Subwavelength Nonlinear Phase Control and Anomalous Phase Matching in Plasmonic Metasurfaces. Nat. Commun. 2016, 7, 10367. (46) Arbabi, A.; Horie, Y.; Ball, A. J.; Bagheri, M.; Faraon, A. Subwavelength-Thick Lenses with High Numerical Apertures and Large Efficiency Based on High-Contrast Transmitarrays. Nat. Commun. 2015, 6, 7069. (47) Black, L.-J.; Wang, Y.; de Groot, C. H.; Arbouet, A.; Muskens, O. L. Optimal Polarization Conversion in Coupled Dimer Plasmonic Nanoantennas for Metasurfaces. ACS Nano 2014, 8, 6390. (48) Chen, P.-Y.; Argyropoulos, C.; D’Aguanno, G.; Alù, A. Enhanced Second-Harmonic Generation by Metasurface Nanomixer and Nanocavity. ACS Photonics 2015, 2, 1000. (49) Karimi, E.; Schulz, S. A.; De Leon, I.; Qassim, H.; Upham, J.; Boyd, R. W. Generating Optical Orbital Angular Momentum at Visible Wavelengths Using a Plasmonic Metasurface. Light: Sci. Appl. 2014, 3 (5), e167. (50) Pors, A.; Albrektsen, O.; Radko, I. P.; Bozhevolnyi, S. I. Gap Plasmon-Based Metasurfaces for Total Control of Reflected Light. Sci. Rep. 2013, 3, 02155. (51) Wan, W.; Gao, J.; Yang, X. Full-Color Plasmonic Metasurface Holograms. ACS Nano 2016, 10, 10671. (52) Chong, K. E.; Hopkins, B.; Staude, I.; Miroshnichenko, A. E.; Dominguez, J.; Decker, M.; Neshev, D. N.; Brener, I.; Kivshar, Y. S. Observation of Fano Resonances in All-Dielectric Nanoparticle Oligomers. Small 2014, 10 (10), 1985−1990. (53) Cardenas, J.; Poitras, C. B.; Robinson, J. T.; Preston, K.; Chen, L.; Lipson, M. Low Loss Etchless Silicon Photonic Waveguides. Opt. Express 2009, 17 (6), 4752−4757. (54) Nezhad, M. P.; Bondarenko, O.; Khajavikhan, M.; Simic, A.; Fainman, Y. Etch-Free Low Loss Silicon Waveguides Using Hydrogen Silsesquioxane Oxidation Masks. Opt. Express 2011, 19 (20), 18827− 18832. (55) Naiman, A.; Desiatov, B.; Stern, L.; Mazurski, N.; Shappir, J.; Levy, U. Ultrahigh-Q Silicon Resonators in a Planarized Local Oxidation of Silicon Platform. Opt. Lett. 2015, 40 (9), 1892−1895. (56) Goykhman, I.; Desiatov, B.; Khurgin, J.; Shappir, J.; Levy, U. Locally Oxidized Silicon Surface-Plasmon Schottky Detector for Telecom Regime. Nano Lett. 2011, 11 (6), 2219−2224. (57) Mayergoyz, I. D.; Fredkin, D. R.; Zhang, Z. Electrostatic (Plasmon) Resonances in Nanoparticles. Phys. Rev. B: Condens. Matter Mater. Phys. 2005, 72 (15), 155412. (58) Kerker, M.; Wang, D.-S.; Giles, C. L. Electromagnetic Scattering by Magnetic Spheres. J. Opt. Soc. Am. 1983, 73 (6), 765−767.

2366

DOI: 10.1021/acsphotonics.7b00688 ACS Photonics 2017, 4, 2359−2366