Reconfiguring structured light beams using nonlinear metasurfaces
Yun Xu, Jingbo Sun, Jesse Frantz, Mikhail I. Shalaev, Wiktor Walasik, Apra Pandey, Jason D. Myers, Robel Y. Bekele, Alexander Tsukernik, Jasbinder S. Sanghera, Natalia M. Litchinitser
11 Reconfiguring structured light beams using nonlinear metasurfaces
Yun Xu, Jingbo Sun, Jesse Frantz, Mikhail I. Shalaev, Wiktor Walasik, Apra Pandey, Jason D. Myers, Robel Y. Bekele, Alexander Tsukernik, Jasbinder S. Sanghera, Natalia M. Litchinitser Department of Electrical Engineering, the State University of New York, University at Buffalo, Buffalo, NY 14260, USA US Naval Research Laboratory, 4555 Overlook Ave., SW, Washington, DC 20375 University Research Foundation, Greenbelt, MD, USA 20770, USA CST of America, LLC, 3979 Freedom Circle, Suite 750, Santa Clara, CA 95054, USA Toronto Nanofabrication Centre, University of Toronto, Toronto, ON M5S 3G4, Canada * e-mail: [email protected]
Abstract:
Ultra-compact, low-loss, fast, and reconfigurable optical components, enabling manipulation of light by light, could open numerous opportunities for controlling light on the nanoscale. Nanostructured all-dielectric metasurfaces have been shown to enable extensive control of amplitude and phase of light in the linear optical regime. Among other functionalities, they offer unique opportunities for shaping the wave front of light to introduce the orbital angular momentum (OAM) to a beam. Such structured light beams bring a new degree of freedom for applications ranging from spectroscopy and micromanipulation to classical and quantum optical communications. To date, reconfigurability or tuning of the optical properties of all-dielectric metasurfaces have been achieved mechanically, thermally, electrically or optically, using phase-change or nonlinear optical materials. However, a majority of demonstrated tuning approaches are either slow or require high optical powers. Arsenic trisulfide (As S ) chalcogenide glass offering ultra-fast and large nonlinearity as well as a low two-photon absorption coefficient in the near and mid-wave infrared spectral range, could provide a new platform for the realization of fast and relatively low-intensity reconfigurable metasurfaces. Here, we design and experimentally demonstrate an As S chalcogenide glass based metasurface that enables reshaping of a conventional Hermite-Gaussian beam with no OAM into an OAM beam at low-intensity levels, while preserves the original beam’s amplitude and phase characteristics at high-intensity levels. The proposed metasurface could find applications for a new generation of optical communication systems and optical signal processing. KEYWORDS:
Nonlinear optics, structured light, all-dielectric metasurface The discovery of the fact that photons can carry an orbital angular momentum (OAM) opened a new area of optical physics and led to new understanding of a wide range of phenomena [1-4]. The OAM beams possess an azimuthal phase dependence of il exp , where the angle is the azimuthal coordinate and the quantized topological charge is denoted by l ℤ [5]. The OAM beams find important applications, including light-atom interactions [6], manipulation of microscopic objects [7], imaging [8,9], and optical communications [10-13]. Conventionally, OAM beams are generated using spiral phase plates or spatial light modulators [1]. These bulk-optics based devices suit laboratory experiments, but may not be compatible with integrated optics systems that require ultra-compact and flat optical components. Recently, first steps toward the realization of microscale, planar optical components based on liquid crystal technology (q-plates) [14] and optical metasurfaces [16-28] have been made. While many of the first demonstrations of these devices were designed for generating OAM beams with a specific and fixed topological charge, reconfigurability is one of the desired characteristics allowing switching from one charge to another or from the OAM beam to a beam not carrying an OAM. More recently, electrically tunable q-plates [29] as well as mechanically, electrically, thermally and optically tunable metasurfaces have been demonstrated [30-39]. Nevertheless, a majority of the tuning approaches demonstrated to date were not able to simultaneously realize ultra-fast speed and high efficiency. Nonlinear optical interactions offer a promising way for ultra-fast, picosecond scale, and efficient optical switching, tuning and reconfiguration [40-44]. Recently, several material platforms, including silicon [39,40,43] and GaAs, have been used to realize nonlinear light-matter interactions in optical metasurfaces [44-46]. However, while silicon and GaAs have both high linear refractive indices and Kerr nonlinearities, and low linear losses in the telecommunication range, their two-photon absorption (TPA) coefficients are large, resulting in a low figure of merit. On the contrary, As S chalcogenide (ChG) glass displays a very good nonlinear figure of merit [47] in both the near-infrared and the mid-wave infrared spectral bands [48,49]. Moreover, its excellent nano-structuring properties, which enable patterning with resolution superior to polymer photoresists, enable a resist-free approach that simplifies fabrication of the As S -based devices to a single-step process [50]. Therefore, we choose the ChG-glass platform to realize a reconfigurable metasurface. In this work, we design, fabricate, and experimentally demonstrate an As S metasurface capable of reshaping a conventional Hermite-Gaussian beam with no OAM into an OAM beam at low intensity levels, while preserving the original beam’s amplitude and phase characteristics at high intensity levels. Such intensity dependent performance is enabled by the Kerr nonlinearity of ChG glass and carefully designed metasurface that relies on guided resonances of the nanoholes made in the ChG thin film. Figure 1.
Working principle of nonlinear metasurface with reconfigurable output beam.
The color maps show the phase of the input and output beams, and the relative phase introduced by the metasurface. a , At low intensity, upon transmission through the metasurface, the input Hermite-Gaussian beam acquires a non-uniform phase distribution leading to generation of a beam carrying orbital angular momentum. b , At high intensity, the phase difference introduced by the metasurface is uniform; the output beam remains Hermite-Gaussian and does not possess OAM. The basic principle of operation of the proposed metasurface device is illustrated in Fig. 1. The input Hermite-Gaussian beam is transmitted through the metasurface and acquires different phase distribution depending on the input-light intensity. We design our metasurface such that, in the low-intensity regime, the phase acquired in the even quadrants (II and IV) is larger than for the odd quadrants (I and III). The phase distribution directly after the metasurface has a stepwise profile with the values given in row 2 in Table. 1. After propagation, this step-wise phase change smoothens and becomes a continuously varying phase that characterizes OAM-carrying beams. Therefore, in the low-intensity regime, the HG beam, upon transmission through our metasurface, is transformed into a beam that carries OAM. For a high intensity input beam, the phase introduced by the metasurface is uniform, and the input HG beam maintains its phase and intensity distribution and does not acquire the OAM. The output beam reconfigurability is enabled by the design of the metasurface described in detail below, that uses highly nonlinear ChG glass. Switching between low- and high-intensity regimes allows for dynamic introduction of the OAM in the beam. Figure 2. Design of the nonlinear reconfigurable metasurface. a , Metasurface consisting of a square lattice of circular holes perforated in a film of arsenic sulfide (As S ) ChG glass. The film is deposited on a glass substrate and is covered with a layer of the Poly(methyl methacrylate) (PMMA) which also fills the holes. The geometric parameters of the metasurface are: film thickness, h nm; lattice constant, a nm; and the hole diameter is denoted by d . The refractive index of the ChG film is n . The refractive indices of glass and PMMA are PMMA g n n , making the entire structure symmetric. b , c , Spectral dependencies of the transmittance and the phase of the transmitted light as a function of the hole diameter, revealing the presence of guided resonances (see Supplementary Materials for more details). The dashed line indicates the working wavelength nm. The two black crosses indicate the location of the odd and even resonances whose electric- and magnetic-field distributions are shown in panels d and e , respectively. f , g Transmittance and phase of the transmitted light as a function of the hole diameter, for nm and two different values of the refractive index of As S . The two dashed lines at d nm and d nm indicate the diameters used in the final design. The functionality described above is enabled by a ChG-glass-based metasurface illustrated in Fig. 2a. A square array of holes is perforated in a ChG film deposited on a glass substrate. In order to render the structure symmetric, the metasurface is covered with Poly(methyl methacrylate) (PMMA) which also fills the holes. The desired phase difference between the four quadrants of the metasurface is obtained by careful optimization of the geometric parameters. The thickness of the film h nm was a result of the fabrication procedure. The remaining parameters, the lattice constant a and the hole diameter d , were selected based on simulations performed using CST MICROWAVE STUDIO. For the structure to operate in the near-infrared, the lattice constant was chosen to be a nm. Figures d nm, we observe two resonances (odd and even) located at wavelengths o nm and e nm, respectively where the width of the resonance is approximately nm. These resonances originate from the interaction of the incident wave with the modes of the perforated ChG film. For normal incidence, only the doubly degenerate guided resonances can be excited [51], as only for these modes the symmetry of the in-plane electric field components allows for coupling with the incident plane wave [52]. The electric and magnetic field distributions in the metasurface for the light interacting with the odd and even modes are shown in Figs. 2d, 2e. The relative spectral position of the even and odd resonances can be controlled by changing the hole diameter. For d nm, when the resonances are located apart from each other, we observed two dips in transmittance, at wavelengths o and e . In the vicinity of these wavelengths, the phase of the transmitted light changes rapidly by . As the diameter of the hole increases, the central wavelengths of the odd and even guided resonances become closer to each other. When the spectral separation between the resonances is smaller than their width, we observe only one dip in transmittance, and the corresponding phase change spans . For the hole diameter d nm, the resonance is located around the wavelength nm. We choose this wavelength to be the operating wavelength for our nonlinear metasurface. Such combination of the hole diameter d and wavelength allows for effective modification of the phase of the transmitted light while small modification of the refractive index ChG n are introduced. The refractive index of the ChG glasses can be efficiently changed as a function of the input light intensity due to their large nonlinear Kerr response described by ChG n n n n n I . Here, n is the linear refractive index of the ChG glass, n cm /W is the nonlinear coefficient, and n is the refractive index change of the ChG glass corresponding the input beam intensity I . Figures 2f, 2g show the dependence of the transmittance and the phase of the transmitted light on the hole diameter d , at the wavelength , for two different values of the linear refractive index of the ChG glass. Here, the modification of the linear refractive index of the ChG glass is assumed to be uniform in the entire film, and is not a result of nonlinear effects. These simulations are used as guidance for our design, and the validity of this approach is confirmed by the full nonlinear results presented in Fig. 3. As illustrated in Fig. 1, our metasurface design requires the relative phase shift between the odd and even quadrants to be in the low-intensity regime, and in the high-intensity regime. Based on the result shown in Fig. 2g, we observe that these two conditions are fulfilled for the choice of the hole diameter o d nm and e d nm, where o d and e d denote the hole diameters in odd and even quadrants, respectively. For the structure with o d d , the resonance is located far away from the operation wavelength, and the structure is insensitive to small change in refractive index. At the same time, for the chosen parameters the transmittance remains above , as seen in Fig. 2f. Figure 3.
Nonlinear simulations of a plane-wave propagation through a unit cell with light intensities used in experiment. a , Schematic of the unit cell used in the simulations showing the propagation direction of the incident y -polarized light. b , c , Electric field distributions in the x plane (center of the unit cell) for low ( b ) and high ( c ) input intensities and for hole diameters o d (left) and e d (right) normalized to the maximum input beam intensity. Black horizontal line in the top part of the simulation domain helps visualize the relative phase shift of the transmitted light. The color map for resonant structures is saturated. d , e , Distribution of the nonlinear index modification in the z plane (center of the ChG film) corresponding to situations illustrated in b , c . In order to confirm that the predictions based on uniform linear index changes give us a correct design, we performed nonlinear simulations for the hole arrays with the chosen hole sizes o d and e d . The results shown in Fig. 3b confirm that for the low intensity (intensity of the incident plane wave equal to MW/cm , which corresponds to the intensity used in the experiment), the phase difference between the odd and even quadrants is equal to 88°. When the intensity is increased to .2 GW/ cm , the waves transmitted through both structures are in phase, as illustrated in Fig. 3c. The maps of the nonlinear refractive index modification n indicate that in the case of the non-resonant structure ( o d d ), the maximum n , while for the resonant structure ( e d d ), the maximum n . Despite the non-uniform distribution of the refractive index change in the nonlinear case, the phase changes predicted using the results shown in Fig. 2g are in good agreement with the full nonlinear simulations. The nonlinear simulations were performed in COMSOL Multiphysics. Figure 4. Fabrication procedure and the fabricated sample. a , Fabrication process for ChG-based all dielectric metasurfaces involves the following steps: fabrication of a chromium window with size of 150 µm×150 µm; deposition of an As S film; exposure of the As S using electron beam; development of the exposed sample; and spin-coating of a PMMA layer. b , Scanning electron microscopy image of the fabricated metasurface with four quadrants, before the layer of PMMA was spin coated. To realize the functionality described in Fig. 1, a metasurface sample consisting of four quadrants perforated with nano-holes of diameters o d and e d was fabricated. The nano-hole arrays designed according to the guided-resonances theory outlines above were patterned in ChG thin film using electron-beam lithography, as shown in Fig. 4. First, chromium windows were prepared on a glass substrate. Then, a ChG film was deposited on top of the chromium windows with thermal deposition. The linear refractive index of the film was measured using spectroscopic ellipsometry to be n . A solution of diluted MF-319 developer was then used to develop the sample. The developed ChG sample with a thickness of 375 nm is shown in Fig. 4b. Finally, a layer of PMMA was spin coated to fill the holes and cover the fabricated structure in order to obtain a symmetric structure. The fabricated structure contains four quadrants where diagonal quadrants contain holes with the same diameter. The total size of the fabricated metasurface is
132 132 µm . Figure 5.
Metasurface characterization. a , Schematic of the Mach–Zehnder interferometer with a delay line used to characterize the fabricated metasurface. The phase plate and the lens in the dashed rectangle were added in the second step as described in the text. b , Interferogram of the light from the sample plane and the reference beam. c , The y -averaged intensity distributions in the regions enclosed by rectangles shown in b with the colors corresponding to the color of the intensity curve. The relative shift in the peak positions of the red and blue curves allows us to estimate the relative phase shift introduces in the four quadrants, as described in the text. To characterize the fabricated metasurface, we built a Mach-Zehnder interferometer with a delay line as shown in Fig. 5a. A high-power femtosecond laser system (Coherent Inc.) with optical parametric amplifier generating 100-fs-long pulses with the repetition rate of 1 kHz was used as a light source, with the central wavelength set to 1255 nm. To retrieve the phase difference between the different quadrants of the fabricated metasurface, a low-power collimated Gaussian beam from the laser was split into two: one transmitted through the sample and the other used as a reference beam. The interference pattern of the two beams was then recorded by a camera, and the overlap between the two pulses was controlled by the delay line. The interferogram is shown in Fig. 5b and the phase discontinuities are observed between different quadrants and at the sample outline. The y -averaged intensities of the top and bottom parts of the sample are plotted along the x -axis of the interferogram in Fig, 5c to determine the phase difference between the neighboring quadrants. For x , the red curve is shifted toward positive x -direction as compared to the blue curve, indicating that the light transmitted through the quadrant II is delayed with respect to the light transmitted by quadrant III. On the contrary, for x , the red curve is shifted toward negative x -direction as compared to the blue curve, indicating that the light transmitted through the quadrant I is advanced with respect to the light transmitted by quadrant IV. The magnitude of this shift is estimated to be a quarter of the fringe period, which corresponds to a phase shift. This confirms the simulation results presented in Fig. 3b. Figure 6. Experimental and simulation results of a Hermite-Gaussian beam transmitted through the metasurface. a , Main beam with low intensity after propagating through the metasurface. c , Main beam with high intensity after propagating through the metasurface. b and d are interference of Gaussian beam and beams shown in a and c , respectively. e - j, simulation results of a Hermite-Gaussian beam transmitted through the metasurface with the low refractive index and the high refractive index solved by CST time domain solver. e - f , Intensity and phase of the simulated Hermite-Gaussian beam transmitted through the metasurface in linear regime. g , Interference of a Gaussian beam and the beam shown in the e-f . h - i , Intensity and phase of the simulated Hermite-Gaussian beam transmitted through the metasurface in nonlinear regime with higher refractive index. j , Interference of a Gaussian beam and the beam shown in the h-i . Once the 90° phase difference between the odd and even quadrants was confirmed in the low-intensity regime, a phase plate was added to the path of the main beam to produce a HG beam. The beam was next focused on the sample using a lens, as shown in the dashed rectangle in Fig. 5a, yielding a beam size of 60 µm. A near-infrared camera was used to capture the main beam transmitted through the metasurface and the interferograms of the main and reference beams. In the low-intensity regime, when the incident beam has the average power of 5 nW, which is not enough to induce a significant nonlinear index change, odd and even quadrants maintain the 90° phase difference as shown in Figs. 2f, 2g. Upon transmission through such metasurface with two quadrants introducing a 90° phase difference, the HG beam acquires an OAM. The transverse intensity distribution of the main beam transmitted through the metasurface is shown in Fig. 6a, and the pattern resulting from the interference of the main beam and the reference Gaussian beam is shown in Fig 6b. The interference pattern reveals has spiral shape, typical for the OAM beams, due to the helical geometry of the phase front. In the nonlinear regime, when the incident beam has the average power of 2 µW (peak intensity around 1.1 GW/cm2), the 90° phase shift in the even quadrants is removed due to the nonlinear response of the metasurface. Therefore, the phase distribution across the entire metasurface becomes uniform and does not affect the phase of the transmitted HG beam. The transverse intensity distribution of the high-intenity main beam transmitted through the metasurface is shown in Fig. 6c, and the corresponding interference pattern with the reference beam is shown in Fig. 6d. The interferogram features semi-circular shapes characteristic for interference of a HG and a Gaussain beams. This result clearly demonstrates the possibility of dynamic introduction of an OAM into a HG beam that can be controlled by the intensity of the input radiation. The transmittance of the metasurface sample in the low and high intensity regimes were measured to be 50% and 53%, respectively. The simulation results corresponding to the linear and nonlinear regimes are shown in Figs. 6e-g and Figs. 6h-j, respectively. In Fig. 6e-f, we observe a dark beam center and spiral phase characteristic for the beams carrying an OAM. The spiral fringe in Fig. 6g matches the interference experiment shown in Fig. 6b. When the refractive index is increased by 0.04, the beam maintains its HG profile as shown in Fig. 6h-i. The interference pattern in Fig. 6j is corresponding to Fig. 6d. The simulations were performed in CST Microwave Studio with reduced model size of and then propagated in free space for 2.6 mm with the Fast Fourier Transform Beam Propagation Method [53]. The simulations were solved at the wavelength of 1259 nm. Due to the limited computational resources and time, we use a time domain solver and set the mesh to 15 cells per wavelength. With this setting, the 90° phase difference moves 4 nm away from the designed working wavelength. To conclude, we have experimentally demonstrated that the light transmitted through an all-dielectric metasurface can be reshaped between a beam with no OAM and a beam carrying a well-defined OAM, dependent on the intensity of the input light. The reconfigurability is enabled by the strong Kerr effect in the ChG film, leading to an ultra-fast response. The nonlinear response of the film was enhanced by the use of a resonant design of the metasurface pattern. Highly confined energy density inside the resonant chalcogenide structure results in a reconfigurable metasurface operating at low intensity. Moreover, the metasurface is characterized by a high transmission efficiency around 50% in both the low- and high-intensity modes of operation. Finally, the use of the ChG-glass platform allows for simplification of the fabrication procedure as the As S glass itself acts as an e-beam resist. This nonlinear reconfigurable metasurface will find applications in ultra-fast and compact optical communication systems, optical signal processing and all optical switching. Methods Design.
We used CST Microwave Studio Frequency Solver to design the linear metasurface. The refractive indices of materials were measured using a spectroscopic ellipsometer. The nonlinear simulation is performed by Comsol Multiphysics. The nonlinear coefficient is measured by a home-built Z-scan setup.
Sample fabrication.
An array of square chromium windows with size 150 150 m m was fabricated on a glass substrate. Then, an As S film with thickness of 500 nm was deposited on top via thermal evaporation in a Lesker PVD 75 deposition system equipped with a low temperature evaporation source. During deposition, the substrate temperature was maintained at approximately 20° C. Inside each of the Cr windows, the ChG glass was patterned with a square array of holes using electron beam (Vistec EBPG5000+ 100KV) with dosage 10.5 mC/cm . The area of the pattern is 132 132 m m , corresponding to the final sample size. The chromium around the sample plays a double role: (i) it enhances the reflectivity of the substrate enabling automatic sample alignment, and (ii) it increases the sample conductivity allowing to avoid accumulation of charge. After exposure, the sample is immersed in a 1:1 mixture of Microposit MF-319 developer and deionized water for 32 seconds to develop. After development, the exposed area is left and the remaining thickness is 375 nm characterized by atomic force microscopy. Finally, six layers of PMMA was spin-coated on the sample to provide a symmetric refractive index. Experiment.
To characterize the fabricated metasurface, we built a Mach-Zehnder interferometer with a delay line as shown in Figure 5a. To generate a Hermite-Gaussian mode, a phase plate was inserted in the main beam path. The phase plate was a glass substrate spin-coated with a layer of S1813 photoresist to delay the beam path for half wavelength. The photoresist on half of the substrate was removed by photolithography. The edge of the photoresist was placed at the center of the main beam and a spatial light filter was placed after the phase plate to pass the Hermite-Gaussian mode. References Padgett, M. J. Orbital angular momentum 25 years on [Invited].
Opt. Express , 11265-11274 (2017) 2. Padgett, M. Light’s twist.
Proc. Math. Phys. Eng. Sci. , 20140633 (2014). 3.
Yao, M. & Padgett, M. J. Orbital angular momentum: origins, behavior and applications.
Adv. Opt. Photonics , 161–204 (2011). 4. Litchinitser, N. M. Applied physics. Structured light meets structured matter.
Science , 1054–1055 (2012). 5.
Allen, L., Beijersbergen, M. W., Spreeuw, R. J. C. & Woerdman, J. P. Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes . Phys. Rev. A , 8185–8189 (1992). 6. Babiker, M., Power, W. L. & Allen, L. Light-induced torque on moving atoms.
Phys. Rev. Lett. , 1239–1242 (1994). 7. Friese, M. E. J., Nieminen, T. A., Heckenberg, N. R. & Rubinsztein-Dunlop, H. Optical alignment and spinning of laser-trapped microscopic particles.
Nature , 348–350 (1998). 8.
Fürhapter, S., Jesacher, A., Bernet, S. & Ritsch-Marte, M. Spiral phase contrast imaging in microscopy.
Opt. Express , 689–694 (2005). Xu, Y., Sun, J., Walasik, W. & Litchinitser, N. M. Probing metamaterials with structured light.
Opt. Express , 26249-26254 (2016). 10. Tamburini, F. et al. Encoding many channels on the same frequency through radio vorticity: first experimental test.
New J. Phys. , 033001 (2012). 11. Wang, J. et al. Terabit free-space data transmission employing orbital angular momentum multiplexing.
Nat. Photonics , 488–496 (2012). 12. Willner, A. E. et al. Optical communications using orbital angular momentum beams.
Adv. Opt. Photonics , 66 (2015). 13. Bozinovic, N. et al. Terabit-Scale Orbital Angular Momentum Mode Division Multiplexing in Fibers.
Science , 1545 (2013). 14.
Marrucci, L., Manzo, C. & Paparo, D. Optical spin-to-orbital angular momentum conversion in inhomogeneous anisotropic media.
Phys. Rev. Lett.
Bomzon, Z., Biener, G., Kleiner, V. & Hasman, E. Space-variant Pancharatnam-Berry phase optical elements with computer-generated subwavelength gratings.
Opt. Lett.
Yu, N. & Capasso, F. Flat optics with designer metasurfaces.
Nat. Mater . , 139-150 (2014). 17. Yu, N. et al. Light propagation with phase discontinuities: Generalized laws of reflection and refraction.
Science , 333-337 (2011). 18.
Kim, S. W., Yee, K. J., Abashin, M., Pang, L. & Fainman, Y. Composite dielectric metasurfaces for phase control of vector field.
Opt. Lett. , 2453-2456 (2015). 19. Minovich, A. E. et al. Functional and nonlinear optical metasurfaces.
Laser Photon. Rev . , 195-213 (2015). Wu, C. H. et al. Spectrally selective chiral silicon metasurfaces based on infrared Fano resonances.
Nat. Commun . , 1-9 (2014). 21. Ni, X. J., Kildishev, A. V. & Shalaev, V. M. Metasurface holograms for visible light.
Nat. Commun . , 1-6 (2013). 22. Buchnev, O., Podoliak, N., Kaczmarek, M., Zheludev, N. I. & Fedotov, V. A. Electrically controlled nanostructured metasurface loaded with liquid crystal: toward multifunctional photonic switch.
Adv. Opt. Mater . , 674-679 (2015). 23. Sautter, J. et al. Active tuning of all-dielectric metasurfaces.
ACS Nano , 4308-4315 (2015). 24. Pfeiffer, C. et al. Efficient light bending with isotropic metamaterial Huygens' surfaces,
Nano Lett. , 2491-2497 (2014). 25. Holloway, C. L. et al. An overview of the theory and applications of metasurfaces: The two-dimensional equivalents of metamaterials.
IEEE Antennas Propag. Mag. , 10-35 (2012). 26. Lin, D., Fan, P., Hasman, E. & Brongersma, M. L. Dielectric gradient metasurface optical elements.
Science , 298-302 (2014). 27.
Yin, X. B., Ye, Z. L., Rho, J., Wang, Y. & Zhang, X. Photonic spin Hall effect at metasurfaces.
Science , 1405-1407 (2013). 28.
Karimi, E. et al. Generating optical orbital angular momentum at visible wavelengths using a plasmonic metasurface.
Light: Sci. Appl . , e167 (2014). 29. Slussarenko, S. et al. Tunable liquid crystal q-plates with arbitrary topological charge.
Opt. Express , 4085-4090 (2011). Chu, C. H. et al. Active dielectric metasurface based on phase‐change medium.
Laser Photonics Rev . , 986-994 (2016). 31. Li, G., Zhang, S. & Zentgraf, T. Nonlinear photonic metasurfaces.
Nat. Rev. Materials , 1-14 (2017). 32. Ou, J. Y., Plum, E., Jiang, L. & Zheludev, N. I. Reconfigurable photonic metamaterials.
Nano Lett. , 2142–2144 (2011). 33. Kamali, S. M., Arbabi, E., Arbabi, A., Horie, Y. & Faraon, A. Highly tunable elastic dielectric metasurface lenses.
Laser Photon. Rev . , 1002-1008 (2016). 34. Su, Z., Zhao, Q., Song, K., Zhao, X. & Yin, J. Electrically tunable metasurface based on Mie-type dielectric resonators.
Sci. Rep . , 43026 (2017). 35. Rahmani, M. et al. Reversible thermal tuning of all‐dielectric metasurfaces.
Adv. Functional Mat. , 1700580 (2017). 36. Wang, Q. et al. Optically reconfigurable metasurfaces and photonic devices based on phase change materials.
Nat. Photon . , 60-66 (2016). 37. Komar, A. et al. Electrically tunable all-dielectric optical metasurfaces based on liquid crystals.
Appl. Phys. Lett. , 071109 (2017). 38.
Boyd, R. W.
Nonlinear Optics, 3rd ed . (Academic Press, New York, 2008). 39.
Krasnok, A., Tymchenko, M. & Alù, A. Nonlinear metasurfaces: a paradigm shift in nonlinear optics . Mater. Today , 8-21 (2018). 40. Lapine, M., Shadrivov, I. V. & Kivshar, Y. S. Colloquium: Nonlinear metamaterials.
Rev. Mod. Phys. , 1093 (2014). 41. Shadrivov, I. V., Kapitanova, P. V., Maslovski, S. I. & Kivshar, Y. S. Metamaterials controlled with light.
Phys. Rev. Lett. , 083902 (2012). Pandey. A. & Litchinitser, N. M., Nonlinear light concentrators.
Opt. Lett . , 5238-5240 (2012). 43. Shcherbakov, M. R. et al. Ultrafast all-optical switching with magnetic resonances in nonlinear dielectric nanostructures.
Nano Lett. , 6985-6990 (2015). 44. Shcherbakov, M. R. et al. Ultrafast all-optical tuning of direct-gap semiconductor metasurfaces.
Nat. Commun . , 17 (2017). 45. Liu, S. et al. Resonantly enhanced second-harmonic generation using III–V semiconductor all-dielectric metasurfaces.
Nano Lett. , 5426-5432 (2016). 46. Lee, J. et al. Giant nonlinear response from plasmonic metasurfaces coupled to intersubband transitions.
Nature , 65–69 (2014). 47.
Sanghera, J. S. et al. Nonlinear properties of chalcogenide glass fibers.
Int. J. Appl. Glass SCI , 296-308 (2010). 48. Eggleton, B. J., Luther-Davies, B. & Richardson, K. Chalcogenide photonics . Nat. Photon . , 141–148 (2011). 49. Hilton, A. R.
Chalcogenide glasses for infra-red optics (McGraw Hill, 2010). 50.
Vlcek, M. & Jain, H. Nanostructuring of chalcogenide glasses using electron beam lithography.
J. Optoelectron. Adv. M , 2108-2111 (2006). 51. Fan, S. & Joannopoulos, J. D. Analysis of guided resonances in photonic crystal slabs.
Phys. Rev. B , 235112 (2002). 52. Paddon, P. & Young, J. F. Two-dimensional vector-coupled-mode theory for textured planar waveguides.
Phys. Rev. B
61, 2090 (2000). 53.
Feit, M. D. & Fleck, J. A. Light propagation in graded-index optical fibers.
Appl. Opt. , 3990-3998 (1978). Acknowledgments
This work was supported by the U.S. Office of Naval Research Award N000141613020.
Author contributions
Y. X., J. S., and N. M. L. contributed to the initial idea. Y. X., J. S., M. I. S., W. W., and A. P. performed the numerical simulation. J. F., J. D. M, R. Y. B., and J. S. S. prepared the As S ChG film. Y. X., J. S., M. I. S., and A. T. patterned the nanostructures on the samples. Y. X., J. S., M. I. S., W. W., and N. M. L. wrote the manuscript. N. M. L. supervised the work.
Additional information
Supplementary information is available.
Competing financial interests
The authors declare no competing financial interests. Row Beam Phase [ ] in quadrant I II III IV 1 Input 0 0 180 180 2 Output L 0 90 180 270 3 NL 0 0 180 180 Table 1. Phase of the beam in the four quadrants of the Cartesian coordinate frame before and after transmission of the beam through the metasurface for linear (L) and nonlinear (NL) cases, corresponding to low and high intensity of the input beam, respectively. upplementary information In Fig. 2 d , e in the main text, we showed the electric and magnetic field distributions corresponding to two spectrally separated resonances supported by the structure with the hole diameter d nm. To illustrate the origin of the guided resonances and characterize them, in Fig. S1, we show the band structure, transmittance, and phase as functions of frequency for the structure with the same parameters. Figures S1 a , b show the band diagram of the studied structure with the close-up of the spectral range of interest in the vicinity of the center of the first Brillouin zone ( point). At the point, both even and odd the modes of the metasurface can be represented as a superposition of four the modes (with the same parity) of a uniform slab. Normally incident light impinging on the metasurface does not have an in-plane k-vector component and can couple to the modes of the structure located at the point. The coupling efficiency depends on the symmetry of the field distribution of the modes and on its overlap integral with the field of the incident wave. As shown in Ref. [1, 2], normally incident light can only couple to the degenerate modes at the point. These degenerate odd and even modes are indicated in Fig. S1 b , and their spectral position matches with the transmittance dips visible in Fig. S1 c . Figure S1 d reveals two distinct phase jumps at the resonant frequencies/wavelengths with the amplitudes close to . Figure S1. Band structure, transmittance and phase calculated for a square unit cell with periodicity a = 660 nm, thickness h nm, hole diameter d nm and refractive index n . The refractive index of the surrounding medium (substrate, material inside the hole and on the top of the structure) env n . a . Band structure for the wave vectors along the high symmetry direction M and X in reciprocal space ( , ) x y k k . The band structure inside the black rectangle is enlarged in b . The degenerate modes are indicated with the arrows and corresponding frequencies. The two degenerate modes at point couple to the electromagnetic wave incident in the direction normal to the surface. Transmittance ( c ) and phase ( d ) for the waves propagating normally to the photonic crystal slab. Two dips in transmittance and two abrupt phase changes in phase are introduced at the wavelengths corresponding to the degenerated modes shown in b . or a single-unit-cell calculations, the frequency-domain results presented in Fig. 2 in the main text and in Figs. S1c-d are obtained using the CST frequency-domain solver with convergence-based mesh refinement [3]. For these simulations, the mesh was refined in the subsequent steps until the difference in calculated transmittance in the two consecutive steps is smaller than 1%. This refinement process allows us to precisely determine the position of the resonances of the structure. Finally, we found that the size of the simulation domain necessary to model the actual metasurface used in experiments is prohibitively large for the frequency-domain calculations. The simulation results presented in Fig. 6 were obtained using CST time-domain solver. Nevertheless, even in the time-domain simulations, due to the large size of the simulation domain, the limitation of computational resources did not allow us to use an optimum mesh size. The coarser mesh used in the time-domain simulation resulted in small shift in the position of the resonances of the structure. The comparison of the transmittance and phase calculated using the frequency and the time-domain solvers is shown in Fig. S2. It can be seen that, in the time-domain simulation with mesh containing 15 cells per wavelength in the material, the position of the resonance is shifted by 4 nm with respect to the reference frequency-domain result. As a result, the desired phase shift between the quadrants with different hole sizes occurs at a wavelength of 1259 nm. Therefore, the simulations in Figs. 6 e–h are performed at this wavelength. Figure S2. Simulation results of a unit cell with two different diameters d . a , b , Transmittance ( a ) and phase ( b ) obtained using frequency domain solver with mesh refinement described in the text. c and d , Transmittance ( c ) and phase ( d ) with mesh setting 15 cells per effective wavelength solved by time domain solver. References Fan, S. & Joannopoulos. J. D. Analysis of guided resonances in photonic crystal slabs.
Physical Review B , 235112 (2002). . Paddon, P. & Young, J. F. Two-dimensional vector-coupled-mode theory for textured planar waveguides.
Phys. Rev. B , 2090 (2000). 3., 2090 (2000). 3.