Nonreciprocal Flat Optics with Silicon Metasurfaces - Nano Letters

Jan 25, 2018 - The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.nanolett.7b04646. Numerical ...
2 downloads 21 Views 1MB Size
Subscriber access provided by MT ROYAL COLLEGE

Communication

Nonreciprocal flat optics with silicon metasurfaces Mark Lawrence, David Russell Barton, and Jennifer A. Dionne Nano Lett., Just Accepted Manuscript • DOI: 10.1021/acs.nanolett.7b04646 • Publication Date (Web): 25 Jan 2018 Downloaded from http://pubs.acs.org on January 26, 2018

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.

Nano Letters 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 24 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

Nano Letters

Nonreciprocal flat optics with silicon metasurfaces Mark Lawrence, David R. Barton III, Jennifer A. Dionne Department of Materials Science and Engineering, Stanford University, Stanford, California 94305, USA

ABSTRACT: Metasurfaces enable almost complete control of light through ultra-thin, subwavelength surfaces by locally and abruptly altering the scattered phase. To date, however, all metasurfaces obey time-reversal symmetry, meaning that forward and backward travelling waves will trace identical paths when being reflected, refracted, or diffracted. Here, we use fullfield calculations to design a passive metasurface for nonreciprocal transmission of both direct and anomalously refracted near infrared light over nanoscale optical path lengths. The metasurface consists of a 100-nm-thick, periodically-patterned Si slab. Owing to the highquality-factor resonances of the metasurface and the inherent Kerr nonlinearities of Si, this structure acts as an optical diode for free-space optical signals. This structure also exhibits nonreciprocal anomalous refraction with appropriate patterning to form a phase gradient metasurface. Compared to existing schemes for breaking time-reversal symmetry, our platform enables subwavelength nonreciprocity for arbitrary free-space optical inputs and provides a straightforward path to experimental realization. The concept is also generalizable to other metasurface functions, providing a foundation for one-way lensing and holography.

ACS Paragon Plus Environment

Nano Letters 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

KEYWORDS: nonreciprocal optics, all-dielectric metasurface, nonlinear metasurface, beamsteering The traditional rules of optics have undergone a revolution in recent years, owing in part to the emergence of metasurfaces - two-dimensional arrays of carefully designed nanoantennas1–3. Tasks that previously demanded the path lengths of many complex light rays be accounted for simultaneously, such as lensing4,5, holography6–8 and polarization manipulation9–14, can now be achieved with devices just tens to hundreds of nanometers thick. Such ‘flat’ optics is possible thanks to the ability for the amplitude, polarization and, most importantly, phase of light scattered from each subwavelength antenna to be precisely engineered. Almost every conceivable transfer function has been realized, from diffraction-limited and chromatic aberration-corrected lensing15,16 to retroreflection17 and even invisibility cloaking18, but one important constraint which has yet to be addressed is reciprocity. In other words, for all demonstrations thus far, the measurement gives the same result if the roles of source and detector are reversed. While seemingly innocuous, this symmetry places strict limits on a number of key optical technologies, including telecommunications, photovoltaics and even quantum information19–21. The difficulty of achieving nonreciprocal transmission in the linear regime, i.e. for low light intensities, stems from a need to break time reversal symmetry. To break time reversal symmetry, an effective motion within a device is needed, allowing it to behave differently for an identical wave propagating in forward or backward directions. Traditionally, magneto-optical effects have been employed, in which the application of a DC magnetic field causes molecular Zeeman splitting, leading to an asymmetric permittivity tensor22. A more recent proposal involves a spatiotemporal refractive index modulation used to mimic the presence of a traveling

ACS Paragon Plus Environment

Page 2 of 24

Page 3 of 24 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

Nano Letters

wave which can impart a fixed amount of linear or angular momentum to an incident photon, similar to the Doppler effect for sound23,24. In both cases, the induced change in the material properties is generally very small at optical frequencies, so long optical path lengths are required to observe reasonable levels of nonreciprocity. Though a great many schemes for miniaturization have been proposed, involving plasmonics25–28, photonic crystals29,30, micro-ring resonators24,31, experimental demonstrations remain limited to the scale of a few hundred microns32–34. Moreover, generating the control bias, i.e. creating a magnetic field with a solenoid or a nonuniform voltage bias with electrical contacts, adds further to the complexity and footprint of such devices. Here, we propose an approach to nonreciprocally manipulate free-space near-infrared light with ultrathin metasurfaces. By exploiting High Quality factor (high Q) resonances in periodic Silicon nanostructures, we theoretically demonstrate one-way transmission through devices as thin as 100nm. By engineering the distribution of these Si nanostructures, as to generate a spatial variation in the transmitted phase, we also show that nonreciprocal beam steering can be achieved. Our approach relies on the nonlinear Kerr effect: With sufficiently high powerdensities, photon-photon interactions mediated by anharmonic oscillations of the electron cloud cause the refractive index of a medium to change in proportion to the local electric field intensity. Our high Q resonators are designed to have highly-asymmetric internal fields in the linear regime, so that the incorporation of a nonlinear refractive index projects the local asymmetry to the far-field, generating nonreciprocal transmission. This phenomenon provides a means to break reciprocity without the need for any sort of external biasing mechanism. Moreover, this platform works for arbitrary optical mode inputs, promising the potential for nonreciprocal beam steering, one-way flat lensing and efficient holography.

ACS Paragon Plus Environment

Nano Letters 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

The nonlinear Kerr effect can be described by including a term in the permittivity of the dielectric material which is proportional to the square of the local electric field strength, () =   + ( ) | ()| . Trapping light over many optical cycles causes the electric field inside a resonator to become greatly amplified compared to the incident wave, leading to far higher index changes compared to travelling wave platforms. Moreover, the scattering efficiency varies dramatically with wavelength across a resonance and, therefore, small shifts in the central wavelength, coming from a change in refractive index, can produce large differences in transmission. Just how pronounced both of these effects are for a given system depends on the associated Q factor. As ( ) of most materials is very small, large Q is typically required to reduce the modulation power threshold to reasonable levels. Thankfully, some of the largest nonlinear indices are found in semiconductors such as Silicon, which, as well as being easy to pattern, also have negligible absorption coefficients at infrared frequencies35. By utilizing the nonlinear Kerr effect, high-Q cavities of Si (and InP) have exhibited optical diode action for waveguide modes36,37. Nonreciprocal transmission has also been realized in similar arrangements with active materials, based on gain saturation38,39, optomechanics40–43 and parametric wave-mixing44. However, these systems achieved low radiative loss rates through physical isolation from other scattering structures; i.e., they relied on evanescent coupling between a nearby waveguide and the cavity. Passive nonlinear nonreciprocity which relies on the precise placement of a subwavelength scatterer within a Fabry-Perot cavity has also been proposed45. Although useful in the Microwave regime, the much weaker nonlinear coefficients and smaller wavelengths make this approach less practical at optical frequencies. When dealing with free-space optical radiation modes, the cavity’s position is no longer a useful degree of

ACS Paragon Plus Environment

Page 4 of 24

Page 5 of 24 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

Nano Letters

freedom. Instead, the symmetries of a structure must be exploited, suppressing the radiative loss via destructive interference. Figure 1a illustrates our chosen geometry to controllably reduce free-space radiative coupling and realize a non-reciprocal high Q cavity. It consists of shallow periodic grooves cut into both sides of a 100 nm thick poly-crystalline Si slab with a refractive index of 3.746. Without grooves, the plane-wave illumination of such a slab, suspended in vacuum, reveals only broad Fabry-Perot (FP) fringes in the transmittance spectrum, as seen in the solid black curve in figure 1b. However, the system also supports guided waves, which are completely isolated from the vacuum solutions because of the need for momentum conservation (grey curve in figure 1b). Superpositions of these guided waves can form bound states, an example of which is denoted by the pink star in figure 1b; this state is calculated by folding the continuous slab dispersion about the fictitious period, p=650 nm, represented by the dashed grey curve in figure 1b. Inserting grooves with a period of 650 nm causes the bound state to leak, appearing as a resonance in transmission (see figure 1c). By keeping the depth of the grooves small compared to the slab thickness, coupling to free-space remains weak, leading to a huge amplification of the internal electric field. The local field is seen to reach 130 times that of the incident plane wave for 10 nm deep grooves in figure 1c. Importantly, Silicon can be considered lossless in the infrared and so the enhancement associated with these so-called guided mode resonances (GMR) can be made arbitrarily large, as indicated by figure 1d which shows the Q factor diverging as the groove depth tends to zero. As such, they have been shown to make an excellent platform for sensing47, modulation48,49 and lasing applications50–52. While the in-plane geometry of a patterned slab is responsible for the overall GMR lifetimes, the relative contribution from leakage above and below the slab can be controlled by altering the

ACS Paragon Plus Environment

Nano Letters 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

mirror symmetry in the out-of-plane direction. As the period has been chosen to align the first TE polarized GMR with the first quarter-wave minimum of the FP background, a completely symmetric perturbation,  = 0, can be seen to produce ideal Fano interference. This highly asymmetric line shape is a signature of a simple superposition between resonant and nonresonant scattering53, illustrated schematically in the left panel of figure2a. Starting with the configuration investigated in figure 1c, figure 2b shows the evolution of the transmittance as the notched region is shifted by  with respect to the symmetric configuration. The transmittance spectrum exhibits both a blue-shift and a reduced maximum. This change can be explained by a direct coupling pathway emerging between the FP mode and the GMR. Importantly, unlike directly exciting it with a plane wave, the FP-mediated excitation of the GMR depends on the incident direction. Consequently, destructive and constructive interference between the various pathways leads to a marked difference in the field enhancement when illuminating the slab from above (i.e. “forward” direction) and below (i.e. “reverse” direction), respectively, as depicted schematically in the right panel of figure 2a. Indeed, figure 2c, which tracks the peak volume average of the electric field within the silicon, shows the ratio of forward to reverse field enhancement growing monotonically with , until  = 14 . In this configuration, the direct and FP mediated GMR excitation strengths are matched, producing maximal asymmetry. While the field asymmetry peaks at  = 14 , the maximum enhancement in the forward direction occurs at  = 4 , after which it begins to fall. This result is consistent with the decrease in Q factor, shown in the inset of figure 2c, which can be compensated for by revisiting the in-plane design of the grooves. The maximum transmission drops inversely with increasing field asymmetry, which has previously been shown to be a result of time reversal symmetry54,55. Each of these features must

ACS Paragon Plus Environment

Page 6 of 24

Page 7 of 24 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

Nano Letters

be considered when building a nonreciprocal device in order to achieve the desired threshold power, operational power range and insertion loss. To realize nonreciprocity over nanoscale optical paths, we include the nonlinear behavior of the slabs. We incorporate the third order susceptibility of Poly-Silicon, ( ) = 2.79 × 10  /



, into our frequency domain

simulations and search for the steady state solutions with increasing incident power46. Figure 3a shows the transmittance spectra through an asymmetric corrugated membrane with  = 10  when illuminated from each side; we consider an incident intensity of 10 !"/# . Comparing the dip in forward transmission with the linear case, the GMR is seen to red-shift by roughly one full-width at half maximum at this intensity. However, the reverse transmission spectrum remains essentially unchanged. Consequently, the device becomes highly nonreciprocal, behaving as a diode for $~1472.69 . Details regarding the discontinuity visible in the forward spectrum are provided in the supplementary information. We next consider the change in transmittance as a function of input power. Analogous to a p-n junction diode, the dashed curves in figure 3b represent effective optical I-V curves for light with a wavelength of $ = 1472.69 . In the linear regime, the slab is highly reflective from both sides. As the incident intensity is increased beyond 2 !"/# the transmittance begins to climb in the forward direction, while remaining very low in the reverse direction. One-way transparency persists up until around 20 !"/# when the membrane also starts to transmit in the reverse direction. This behavior, which is functionally equivalent to the response of an electrical diode at the breakdown voltage, results from an illumination intensity which is high enough to shift the GMR across the wavelength of interest, even in the reverse direction. This structure with  = 10  exhibits highly nonreciprocal behavior over a large range of incident powers, completely blocking light from one side while allowing it to pass through from

ACS Paragon Plus Environment

Nano Letters 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

the other. However, the transmittance in the forward direction is only around 25%. While comparable to previous demonstrations of non-reciprocal transmission, it still represents a fairly significant loss for the incoming signal. To increase the peak transmittance, the field asymmetry must be reduced (see Figure 2). Nevertheless, small amounts of asymmetry can still be sufficient to achieve diode action. Figures 3b-c shows the results for a slab with  = 4  . In this case, above a threshold of about 3 !"/# , high transmittance of around 75% is observed in the forward direction. At the same time, the range of incident powers for which the transmission is suppressed in the reverse direction is reduced compared with  = 10 , (~6 !"/ # ,compared to ~20 !"/# ). This result makes sense, given that under reverse illumination a higher intensity is reached within the Silicon for the  = 4  membrane. Accordingly, modification of the structural asymmetry in this metasurface provides continuous tunability of the non-reciprocal response, from exhibiting diode action over a wide range of illuminating powers to exhibiting low insertion loss over a reduced power range. We now show how this platform can be used to realize a nonreciprocal metasurface. Instead of a continuous Silicon sheet, we consider a periodic structure with period ( ≈ √2$, where $ is the freespace wavelength. Here, two new solutions emerge alongside the zeroth order transmission, with diffraction angles + ≈ ∓45°. The scattering into one of these off normal diffraction orders becomes optimal for a transmitted field with uniform intensity and a linearly varying phase; by using nanoantenna designs that approximate a linear phase function, many groups have observed this phenomenon, which is known as anomalous refraction, or beam steering1,5,56,57. Here, we employ the structure depicted in figure 4a, consisting of a simple Si bar with a height of 700 nm and a variable width. Figure 4c plots the intensities of the various diffraction orders when combining bars with widths of 210 nm, 240 nm and 400 nm, taken from figure 4b, into a single

ACS Paragon Plus Environment

Page 8 of 24

Page 9 of 24 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

Nano Letters

unit cell with a period of 2121 nm. As expected, around the design wavelength of 1500 nm, a normally incident plane wave is redirected to the left (+1st order) with very little energy diffracted into other orders. This behavior is also verified in the field map shown in the right panel of figure 4d. Importantly, as for the continuous films considered in Figure 1, the phase gradient metasurface in Figure 4d also supports bound states made up of guided modes which propagate orthogonally to the plane of diffraction. Next, we introduce a high Q mode into the diffraction spectrum. Following the above prescription, we add periodic grooves, 70 nm in width and 25 nm deep, to each side of the 400 nm wide silicon bar, as illustrated in figure 5a. Subtly perturbing the translationally invariant grating gives rise once again to a sharp GMR which perturbs the linear phase gradient, scattering light in to other diffraction orders. Because the overall change to the geometry is very small, away from the sharp spectral feature the diffraction spectrum remains unaffected. As a result, we are free to adjust the period of the grooves in order to center the GMR close to the operating wavelength of 1500 nm. From the black curve in figure 5b it can be seen that a period of 570 nm produces a resonant wavelength of around 1509 nm. Unlike the slab investigated in the first part of this paper which, without grooves, is invariant under a transformation converting between forward and reverse illumination, the metasurface designed in figure 4 is highly asymmetric with respect to the normal and oblique incident waves being considered. Because of this inherent asymmetry, different electric field strengths build up inside the Si for forward and reverse illumination, defined in figure 5c, even with symmetric grooves, i.e. equal groove depths on each side. Details regarding the consequences of tuning the geometric asymmetry can be found in the supporting information.

ACS Paragon Plus Environment

Nano Letters 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

The nonlinear, nonreciprocal response of this metasurface is shown in figures 5b and d. With an illumination intensity of 8.3 !"/# , forward illumination (normal incidence) generates a pronounced red-shift of the GMR. For reciprocal illumination (i.e., oblique incidence matching the angle of the +1st diffracted order), a much smaller red shift is observed. Therefore, the Kerr modulation of the GMR can be seen to produce nonreciprocal beam steering for $~1509.8 . Specifically, highly efficient beam steering occurs in the forward direction, while diffraction into both the 0th and -1st orders is found in the reverse direction. Figure 5d plots the power dependence of this system, showing highly nonreciprocal behavior at $ = 1509.81  for 6 !"/# 01 12 !"/# . In summary, dielectric metasurfaces represent an exciting platform for breaking the reciprocity of free-space electromagnetic signals. Here, we designed a subwavelength grating which, thanks to the nonlinear Kerr effect, acts as an ultrathin (100-nm thick) optical diode for illumination intensities of a few !"/# . Then, using the same platform but starting with a three-element phase gradient metasurface, we achieved large angle (~45o) nonreciprocal beam steering. While our calculations focused on nonreciprocal beam steering as a proof of principle, this scheme could be easily extended to more complicated transfer functions such as nonreciprocal lensing and holography. Furthermore, the proposed scheme of embedding a high Q resonance within a phase gradient metasuarface is generalizable to other nonlinear phenomena, such as parametric amplification and lasing. And while high Q structures require precise fabrication, extremely high Q GMRs, at and beyond those found in the current manuscript, have been experimentally observed52,58–60. So, with the simplicity of the structures presented and the maturity of Silicon processing, the relevant tolerances for our scheme are well within current technological capabilities. As a final note, we have only considered systems which rely on translational

ACS Paragon Plus Environment

Page 10 of 24

Page 11 of 24 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

Nano Letters

symmetry in at least one dimension, with the excitation being planar in that dimension. However, a number of recent studies have shown that trapped states can be found in structures confined in all three directions61–63. By carefully tailoring a nonlinear resonator to a particular, non-planar incident wave, nonreciprocal freespace devices could be miniaturized even further. ASSOCIATED CONTENT Supporting Information The supporting information is available free of charge. Numerical specifications, simulations confirming robustness of scheme to variations in fabrication and characterization conditions, effects of two photon absorption, structural symmetry breaking with phase gradient metasurface AUTHOR INFORMATION Corresponding Author *E-mail: [email protected] *E-mail: [email protected] Author Contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. Notes The authors declare no competing financial interest

Acknowledgements

ACS Paragon Plus Environment

Nano Letters 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 acknowledge the generous funding from the AFOSR PECASE grant (FA9550-15-1-0006) and NSF EFRI grant (1641109), as well as funding from Northrop Grumman. DB acknowledges support from a Stanford Graduate Fellowship.

References (1)

Yu, N.; Genevet, P.; Kats, M. A.; Aieta, F.; Tetienne, J.-P.; Capasso, F.; Gaburro, Z. Science 2011, 334 (6054), 333–337.

(2)

Yu, N.; Capasso, F. Nat. Mater. 2014, 13 (2), 139–150.

(3)

Zhao, Y.; Liu, X.-X.; Alù, A. J. Opt. 2014, 16 (12), 123001.

(4)

Chen, X.; Huang, L.; Mühlenbernd, H.; Li, G.; Bai, B.; Tan, Q.; Jin, G.; Qiu, C.-W.; Zhang, S.; Zentgraf, T. Nat. Commun. 2012, 3, 1198.

(5)

Lin, D.; Fan, P.; Hasman, E.; Brongersma, M. L. Science 2014, 345 (6194), 298–302.

(6)

Wen, D.; Yue, F.; Li, G.; Zheng, G.; Chan, K.; Chen, S.; Chen, M.; Li, K. F.; Wong, P. W. H.; Cheah, K. W.; Yue Bun Pun, E.; Zhang, S.; Chen, X. Nat. Commun. 2015, 6, 8241.

(7)

Zheng, G.; Mühlenbernd, H.; Kenney, M.; Li, G.; Zentgraf, T.; Zhang, S. Nat. Nanotechnol. 2015, 10 (4), 308–312.

(8)

Ni, X.; Kildishev, A. V.; Shalaev, V. M. Nat. Commun. 2013, 4, 2807.

(9)

Zhao, Y.; Alù, A. Nano Lett. 2013, 13 (3), 1086–1091.

(10)

Yu, N.; Aieta, F.; Genevet, P.; Kats, M. A.; Gaburro, Z.; Capasso, F. Nano Lett. 2012, 12 (12), 6328–6333.

(11)

Ding, F.; Wang, Z.; He, S.; Shalaev, V. M.; Kildishev, A. V. ACS Nano 2015, 9 (4),

ACS Paragon Plus Environment

Page 12 of 24

Page 13 of 24 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

Nano Letters

4111–4119. (12)

Wu, C.; Arju, N.; Kelp, G.; Fan, J. A.; Dominguez, J.; Gonzales, E.; Tutuc, E.; Brener, I.; Shvets, G. Nat. Commun. 2014, 5, ncomms4892.

(13)

Li, Z.; Gokkavas, M.; Ozbay, E. Adv. Opt. Mater. 2013, 1 (7), 482–488.

(14)

Menzel, C.; Helgert, C.; Rockstuhl, C.; Kley, E.-B.; Tünnermann, A.; Pertsch, T.; Lederer, F. Phys. Rev. Lett. 2010, 104 (25), 253902.

(15)

Khorasaninejad, M.; Chen, W. T.; Devlin, R. C.; Oh, J.; Zhu, A. Y.; Capasso, F. Science 2016, 352 (6290), 1190–1194.

(16)

Aieta, F.; Kats, M. A.; Genevet, P.; Capasso, F. Science 2015, 347 (6228), 1342–1345.

(17)

Arbabi, A.; Arbabi, E.; Horie, Y.; Kamali, S. M.; Faraon, A. Nat. Photonics 2017, 11 (7), 415–420.

(18)

Ni, X.; Wong, Z. J.; Mrejen, M.; Wang, Y.; Zhang, X. Science 2015, 349 (6254), 1310– 1314.

(19)

Hadad, Y.; Soric, J. C.; Alu, A. .

(20)

Green, M. A. Nano Lett. 2012, 12 (11), 5985–5988.

(21)

Stannigel, K.; Rabl, P.; Zoller, P. New J. Phys. 2012, 14 (6), 63014.

(22)

Dionne, G. F. J. Appl. Phys. 2005, 97 (10), 10F103.

(23)

Bhandare, S.; Ibrahim, S. K.; Sandel, D.; Hongbin Zhang; Wust, F.; Noe, R. IEEE J. Sel. Top. Quantum Electron. 2005, 11 (2), 417–421.

ACS Paragon Plus Environment

Nano Letters 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

(24)

Sounas, D. L.; Caloz, C.; Alù, A. Nat. Commun. 2013, 4, 2407.

(25)

Chin, J. Y.; Steinle, T.; Wehlus, T.; Dregely, D.; Weiss, T.; Belotelov, V. I.; Stritzker, B.; Giessen, H. Nat. Commun. 2013, 4, 1599.

(26)

Floess, D.; Hentschel, M.; Weiss, T.; Habermeier, H.-U.; Jiao, J.; Tikhodeev, S. G.; Giessen, H. Phys. Rev. X 2017, 7 (2), 21048.

(27)

Davoyan, A. R.; Engheta, N. New J. Phys. 2013, 15 (8), 83054.

(28)

Mousavi, S. H.; Khanikaev, A. B.; Allen, J.; Allen, M.; Shvets, G. Phys. Rev. Lett. 2014, 112 (11), 117402.

(29)

Tao, K.; Xiao, J.-J.; Yin, X. Appl. Phys. Lett. 2014, 105 (21), 211105.

(30)

He, C.; Sun, X.-C.; Zhang, Z.; Yuan, C.-S.; Lu, M.-H.; Chen, Y.-F.; Sun, C. Opt. Express 2013, 21 (23), 28933.

(31)

Yu, Z.; Fan, S.; Fan, S. Nat. Photonics 2009, 3 (2), 91–94.

(32)

Bi, L.; Hu, J.; Jiang, P.; Kim, D. H.; Dionne, G. F.; Kimerling, L. C.; Ross, C. A.; Ross, C. A. Nat. Photonics 2011, 5 (12), 758–762.

(33)

Tzuang, L. D.; Fang, K.; Nussenzveig, P.; Fan, S.; Lipson, M. Nat. Photonics 2014, 8 (9), 701–705.

(34)

Lira, H.; Yu, Z.; Fan, S.; Lipson, M. Phys. Rev. Lett. 2012, 109 (3), 33901.

(35)

Leuthold, J.; Koos, C.; Freude, W. Nat. Photonics 2010, 4 (8), 535–544.

(36)

Fan, L.; Wang, J.; Varghese, L. T.; Shen, H.; Niu, B.; Xuan, Y.; Weiner, A. M.; Qi, M.

ACS Paragon Plus Environment

Page 14 of 24

Page 15 of 24 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

Nano Letters

Science 2012, 335 (6067), 447–450. (37)

Yu, Y.; Chen, Y.; Hu, H.; Xue, W.; Yvind, K.; Mork, J. Laser Photon. Rev. 2015, 9 (2), 241–247.

(38)

Peng, B.; Ozdemir, S. K.; Lei, F.; Monifi, F.; Gianfreda, M.; Long, G. L.; Fan, S.; Nori, F.; Bender, C. M.; Yang, L. Nat. Phys. 2013, 10.

(39)

Chang, L.; Jiang, X.; Hua, S.; Yang, C.; Wen, J.; Jiang, L.; Li, G.; Wang, G.; Xiao, M. Nat. Photonics 2014, 8 (7), 524–529.

(40)

Shen, Z.; Zhang, Y.-L.; Chen, Y.; Zou, C.-L.; Xiao, Y.-F.; Zou, X.-B.; Sun, F.-W.; Guo, G.-C.; Dong, C.-H. Nat. Photonics 2016, 10 (10), 657–661.

(41)

Dong, C.-H.; Shen, Z.; Zou, C.-L.; Zhang, Y.-L.; Fu, W.; Guo, G.-C. Nat. Commun. 2015, 6, 6193.

(42)

Fang, K.; Luo, J.; Metelmann, A.; Matheny, M. H.; Marquardt, F.; Clerk, A. A.; Painter, O. Nat. Phys. 2017, 13 (5), 465–471.

(43)

Ruesink, F.; Miri, M.-A.; Alù, A.; Verhagen, E. Nat. Commun. 2016, 7, 13662.

(44)

Guo, X.; Zou, C.-L.; Jung, H.; Tang, H. X. Phys. Rev. Lett. 2016, 117 (12), 123902.

(45)

Mahmoud, A. M.; Davoyan, A. R.; Engheta, N. Nat. Commun. 2015, 6, 8359.

(46)

Yang, Y.; Wang, W.; Boulesbaa, A.; Kravchenko, I. I.; Briggs, D. P.; Puretzky, A.; Geohegan, D.; Valentine, J. Nano Lett. 2015, 15 (11), 7388–7393.

(47)

Liu, Y.; Zhou, W.; Sun, Y. Sensors 2017, Vol. 17, Page 1861 2017, 17 (8), 1861.

ACS Paragon Plus Environment

Nano Letters 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

(48)

Kim, S. J.; Brongersma, M. L. Opt. Lett. 2017, 42 (1), 5.

(49)

Xiao, G.; Zhu, Q.; Shen, Y.; Li, K.; Liu, M.; Zhuang, Q.; Jin, C. Nanoscale 2015, 7 (8), 3429–3434.

(50)

Kodigala, A.; Lepetit, T.; Gu, Q.; Bahari, B.; Fainman, Y.; Kanté, B. Nature 2017, 541 (7636), 196–199.

(51)

Hsu, C. W.; Zhen, B.; Stone, A. D.; Joannopoulos, J. D.; Soljačić, M. Nat. Rev. Mater. 2016, 1 (9), 16048.

(52)

Wei Hsu, C.; Zhen, B.; Lee, J.; Chua, S.-L.; Johnson, S. G.; Joannopoulos, J. D.; Soljači, M. Nature 2013, 499.

(53)

Fan, S.; Suh, W.; Joannopoulos, J. D. J. Opt. Soc. Am. A 2003, 20 (3), 569.

(54)

Sounas, D. L.; Alù, A. Phys. Rev. Lett. 2017, 118 (15), 154302.

(55)

Wang, K. X.; Yu, Z.; Sandhu, S.; Fan, S. Opt. Lett. 2013, 38 (2), 100.

(56)

Huang, L.; Chen, X.; Mühlenbernd, H.; Li, G.; Bai, B.; Tan, Q.; Jin, G.; Zentgraf, T.; Zhang, S. Nano Lett. 2012, 12 (11), 5750–5755.

(57)

Qin, F.; Ding, L.; Zhang, L.; Monticone, F.; Chum, C. C.; Deng, J.; Mei, S.; Li, Y.; Teng, J.; Hong, M.; Zhang, S.; Alu, A.; Qiu, C.-W. Sci. Adv. 2016, 2 (1), e1501168–e1501168.

(58)

Nicolaou, C.; Lau, W. T.; Gad, R.; Akhavan, H.; Schilling, R.; Levi, O. Opt. Express 2013, 21 (25), 31698.

(59)

Lee, J.; Zhen, B.; Chua, S.-L.; Qiu, W.; Joannopoulos, J. D.; Soljačić, M.; Shapira, O.

ACS Paragon Plus Environment

Page 16 of 24

Page 17 of 24 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

Nano Letters

Phys. Rev. Lett. 2012, 109 (6), 67401. (60)

Liu, Y.; Wang, S.; Zhao, D.; Zhou, W.; Sun, Y. Opt. Express 2017, 25 (9), 10536.

(61)

Taghizadeh, A.; Chung, I.-S. Appl. Phys. Lett. 2017, 111 (3), 31114.

(62)

Grepstad, J. O.; Greve, M. M.; Holst, B.; Johansen, I.-R.; Solgaard, O.; Sudbø, A. Opt. Express 2013, 21 (20), 23640.

(63)

Barton, J. H.; Rumpf, R. C.; Smith, R. W.; Kozikowski, C.; Zellner, P. Prog. Electromagn. Res. B 2012, 41, 269–283.

ACS Paragon Plus Environment

Nano Letters 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 1: Guided mode resonance. a) schematic of proposed setup, periodically patterned silicon slab, t=100nm, w=40nm and period p=650nm, suspended in vacuum and illuminated with plane waves. b) Transmittance spectrum and guided wave dispersion for ∆=t and δ=0. Dotted lines indicate projections of the dispersion into the first Brillouin zone, based on p=650nm, and pink

ACS Paragon Plus Environment

Page 18 of 24

Page 19 of 24 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

Nano Letters

star marks the first bound state with field plotted in inset. c) Transmittance spectrum for t∆=20nm. Inset maps the peak electric field normalized to incident plane wave, scale bar is 100nm. d) Quality factor of resonance with increasing groove depth.

ACS Paragon Plus Environment

Nano Letters 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 2: Guided mode resonance in asymmetric structure. a) schematic of modal interactions. For a symmetric structure, the guided mode and Fabry-Perot resonances simply interfere in the farfield, whereas, for asymmetric structures, the two resonances couple in the nearfield. b) Transmittance spectrum evolution with increasing δ, in 2nm steps. Inset maps the electric field normalized to incident plane wave in forward and reverse directions, bordered with red and blue respectively, for δ=10nm (dashed curve). scale bar is 100nm. c) Red dots plot the volume average of the electric field amplitude normalized to incident wave. Black stars plot the ratio of the volume averages of the electric field for forward vs reverse illumination. Inset plots Q factor of resonance for same δ range. ∆=80nm throughout. Field plots and volume averages are plotted at central resonant wavelength, as defined in supporting information.

ACS Paragon Plus Environment

Page 20 of 24

Page 21 of 24 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

Nano Letters

Figure 3: Nonreciprocity from nonlinear Kerr effect: a and c) Transmittance spectra for incident plane waves with 10kW/cm2 and δ=10nm (a), 5kW/cm2 and δ=4nm (c). Black curves represent linear spectra for reference. b) Transmittance as a function of input power at wavelength of 1.47269µm and δ=10nm (dashed curves), 1.47315µm and δ=4nm (solid curves). Red and blue curves represent forward and reverse illumination, respectively, as defined in figure 1a.

ACS Paragon Plus Environment

Nano Letters 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: Phase gradient metasurface: a) Periodic array of silicon structures, used to design metasurface. b) Transmitted intensity (grey) and phase (black) through geometry shown in figure 4a, as a function of width ‘w’. Stars denote parameters chosen for metasurface elements, w1=210nm, w2=240nm and w3=400nm. c) Diffraction transmittance spectrum. d) schematic of phase gradient metasurface, left panel, and H field map verifying anomalous refraction. Element colours correspond to stars in Figure 4b.

ACS Paragon Plus Environment

Page 22 of 24

Page 23 of 24 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

Nano Letters

Figure 5: Nonreciprocal anomalous refraction: a) Phase gradient metasurface with periodic grooves of period p2=570nm, 70 nm width and 25 nm depth. b) Metasurface diffraction intensity for incident power of 8.3kW/cm2. black curve represents the linear forward +1st intensity for reference, while the colours of the other curves corresponds to arrows shown in c. c) Schematic of diffraction setup through metasurface, with incident angle for reverse illumination matching +1st diffraction angle for forward illumination θ=sin-1λ/p. d) Diffraction intensity at wavelength of 1.50981µm as a function of input power, curve designation same as c.

ACS Paragon Plus Environment

Nano Letters 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

TOC 102x47mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 24 of 24