In Situ Planarization of Huygens Metasurfaces by ... - ACS Publications

Aug 29, 2017 - Department of Applied Physics, the Benin School of Engineering and Computer Science, the Center for Nanoscienceand. Nanotechnology, The...
1 downloads 9 Views 5MB Size
Subscriber access provided by UNIVERSITY OF ADELAIDE LIBRARIES

Article

In-Situ planarization of Huygens metasurfaces by nanoscale local oxidation of silicon Jonathan Bar-David, Noa Mazurski, and Uriel Levy ACS Photonics, Just Accepted Manuscript • DOI: 10.1021/acsphotonics.7b00688 • Publication Date (Web): 29 Aug 2017 Downloaded from http://pubs.acs.org on August 29, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

ACS Photonics is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

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. *Corresponding author: [email protected] Keywords: Huygens Metasurface, Dielectric Metasurface, LOCOS 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 alldielectric Huygens metasurfaces. This 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,

ACS Paragon Plus Environment

ACS Photonics

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

we believe that such planarized metasurfaces will become a viable technology for future applications. INTRODUCTION. Optical metasurfaces are nano-patterned surfaces exhibiting desired optical characteristics depending on the nano patterns and the structures which are defined on the surface1–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 non-resonant. While the latter relies on concepts such as phase accumulation, effective index and geometrical phase to control light2,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 nanoantennas1,11,13–18,28,29. The scattering, transmission and absorption of electromagnetic waves by small particles is an old topic, first considered by Mie30–32, who calculated the inner and outer fields for scattering of plane waves by spheres, concluding they are the sum of the basic multipole fields31,32. The multipole fields are therefore the electromagnetic eigenmodes of the fields in and around small particles31,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 nano-fabrication techniques, these observations pioneered by Mie are now of great significance to the photonic research, as nowadays researchers as attempting to harness the optical response of small particles and cavities such as nanoantennas for 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,

ACS Paragon Plus Environment

Page 2 of 21

Page 3 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

including for example color engineering11,18,19,21,28,33,34, nonlinear optics35–38, reflectionless surfaces and perfect absorbers1,8,39–41, metalenses and other optical elements10,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 effects17,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 polarization1,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 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 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

ACS Paragon Plus Environment

. In this technique, a protective

ACS Photonics

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

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 obtain quasi-planar structures without jeopardizing the smoothness of the obtained structures55. 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.

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

ACS Paragon Plus Environment

Page 4 of 21

Page 5 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

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 quasiplanar metasurface.

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 (ebeam) 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 inter-antenna 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 ±25nm. The relation between antenna layer thickness and final surface

ACS Paragon Plus Environment

ACS Photonics

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

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 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.

ACS Paragon Plus Environment

Page 6 of 21

Page 7 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

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 ~420nm, 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.

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.

ACS Paragon Plus Environment

ACS Photonics

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

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

ACS Paragon Plus Environment

Page 8 of 21

Page 9 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

nanometers in the simulated antenna geometry, e.g. 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 950 to 1500 nm. The gap was kept constant at 200 nm, i.e. the period is increasing with the radius. As may be seen in figure 5(A), 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 e.g.8,31,40,57.

ACS Paragon Plus Environment

ACS Photonics

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 10 of 21

Page 11 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

Figure 5. FDTD simulation of antenna transmission. (A) Transmission spectrum of an antenna array for 600-1700 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 condition58 - a decrease and narrowing of the resonant reflection, is met by antennas of radius ~ 150nm at a wavelength of about 900 nm. The dashed line marks the antenna radius for which we present the spectra in sub-figure (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 inwards. These resonances are actually the two first Mie resonances, which "switch order" - i.e. 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 height of 275 nm and radii of 330 nm. Figure 6 shows the calculated electric and magnetic fields for the two modes indicated in figure 5(B). It may be seen that under illumination with a x-polarized plane wave, the MD modal field at a wavelength of 1306 nmhas the expected strong magneticfield in the yො direction and the 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 theasymmetric 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.

ACS Paragon Plus Environment

ACS Photonics

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 21

Figure 6. Electric and magnetic fields at resonance, for antenna with a 330 nm radius.Light is propagates towards the antenna along the z axis. The black lines denote the antenna crosssection. 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. 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, 50x, NA 0.45) and directed into a spectrometer (Horiba, microHR \

ACS Paragon Plus Environment

Page 13 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

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 circa 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 increasing 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.

Figure 7. Measured transmission spectrum for antenna with radii of 70 to 300 nm. With the reduction of aspect ratio, the resonant modes begin to separate, and for nanoantennas with radius

ACS Paragon Plus Environment

ACS Photonics

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 21

>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.

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 towards 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.

ACS Paragon Plus Environment

Page 15 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

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 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 towards 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 towards longer wavelengths as the antenna length increases.

ACS Paragon Plus Environment

ACS Photonics

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 21

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 inwards.

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-

ACS Paragon Plus Environment

Page 17 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

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’s compatibility with current CMOS techniques. Following the results reported here, one may use this technique for the implementation of planar metasurface devices, e.g. 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. Author Contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. Funding Sources 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. REFERENCES (1)

(2)

(3) (4) (5) (6)

Kuznetsov, A. I.; Miroshnichenko, A. E.; Brongersma, M. L.; Kivshar, Y. S.; Luk’yanchuk, B. Optically Resonant Dielectric Nanostructures. Science 2016, 354 (6314), aag2472. Bomzon, Z.; Biener, G.; Kleiner, V.; Hasman, E. Space-Variant Pancharatnam–Berry Phase Optical Elements with Computer-Generated Subwavelength Gratings. Opt. Lett. 2002, 27 (13), 1141–1143. Desiatov, B.; Mazurski, N.; Fainman, Y.; Levy, U. Polarization Selective Beam Shaping Using Nanoscale Dielectric Metasurfaces. Opt. Express 2015, 23 (17), 22611–22618. Yariv, A. Photonics Optical Electronics in Modern Communications, 6th ed.; Oxford University Press: New York ;;Oxford, 2007. Lin, D.; Fan, P.; Hasman, E.; Brongersma, M. L. Dielectric Gradient Metasurface Optical Elements. Science 2014, 345 (6194), 298–302. 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.

ACS Paragon Plus Environment

ACS Photonics

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(7)

(8) (9)

(10) (11) (12) (13)

(14)

(15) (16) (17)

(18) (19) (20)

(21) (22)

(23)

(24) (25)

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. Iyer, P. P.; Butakov, N. A.; Schuller, J. A. Reconfigurable Semiconductor Phased-Array Metasurfaces. ACS Photonics 2015. 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. Yu, N.; Capasso, F. Flat Optics with Designer Metasurfaces. Nat. Mater. 2014, 13 (2), 139–150. Proust, J.; Bedu, F.; Gallas, B.; Ozerov, I.; Bonod, N. All-Dielectric Colored Metasurfaces with Silicon Mie Resonators. ACS Nano 2016, 10 (8), 7761–7767. Jahani, S.; Jacob, Z. All-Dielectric Metamaterials. Nat. Nanotechnol. 2016, 11 (1), 23–36. 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. 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. Kuznetsov, A. I.; Miroshnichenko, A. E.; Fu, Y. H.; Zhang, J.; Luk’yanchuk, B. Magnetic Light. Sci. Rep. 2012, 2, 492. 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. 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. Hedayati, M. K.; Elbahri, M. Review of Metasurface Plasmonic Structural Color. Plasmonics 2016, 1–17. 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. 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. Huang, Y.-W.; Chen, W. T.; Tsai, W.-Y.; Wu, P. C.; Wang, C.-M.; Sun, G.; Tsai, D. P. Aluminum Plasmonic Multicolor Meta-Hologram. Nano Lett. 2015, 15 (5), 3122–3127. Devlin, R. C.; Ambrosio, A.; Wintz, D.; Oscurato, S. L.; Zhu, A. Y.; Khorasaninejad, M.; Oh, J.; Maddalena, P.; Capasso, F. Spin-to-Orbital Angular Momentum Conversion in Dielectric Metasurfaces. ArXiv160503899 Phys. 2016. 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. Lalanne, P.; Chavel, P. Metalenses at Visible Wavelengths: Past, Present, Perspectives. ArXiv161002507 Phys. 2016. Lalanne, P.; Hazart, J.; Chavel, P.; Cambril, E.; Launois, H. A Transmission Polarizing Beam Splitter Grating. J. Opt. Pure Appl. Opt. 1999, 1 (2), 215.

ACS Paragon Plus Environment

Page 18 of 21

Page 19 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

(26) Lalanne, P.; Hugonin, J. P.; Chavel, P. Optical Properties of Deep Lamellar Gratings: A Coupled Bloch-Mode Insight. J. Light. Technol. 2006, 24 (6), 2442. (27) Lalanne, P. Waveguiding in Blazed-Binary Diffractive Elements. JOSA 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 & 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. (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 Patch-Dipole Nanoantenna Reflectarrays. Nano Lett. 2014. (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.

ACS Paragon Plus Environment

ACS Photonics

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(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. (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. (48) Chen, P.-Y.; Argyropoulos, C.; D’Aguanno, G.; Alù, A. Enhanced Second-Harmonic Generation by Metasurface Nanomixer and Nanocavity. ACS Photonics 2015. (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. (51) Wan, W.; Gao, J.; Yang, X. Full-Color Plasmonic Metasurface Holograms. ACS Nano 2016. (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 AllDielectric 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 2005, 72 (15), 155412. (58) Kerker, M.; Wang, D.-S.; Giles, C. L. Electromagnetic Scattering by Magnetic Spheres. JOSA 1983, 73 (6), 765–767.

ACS Paragon Plus Environment

Page 20 of 21

Page 21 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

TOC GRAPHIC

ACS Paragon Plus Environment