Nonlinear Wavefront Control with All-Dielectric Metasurfaces - Nano

Subscriber access provided by Syracuse University Libraries

Communication

Nonlinear wavefront control with all-dielectric metasurfaces Lei Wang, Sergey S. Kruk, Kirill L. Koshelev, Ivan I. Kravchenko, Barry Luther-Davies, and Yuri S. Kivshar Nano Lett., Just Accepted Manuscript • DOI: 10.1021/acs.nanolett.8b01460 • Publication Date (Web): 11 May 2018 Downloaded from http://pubs.acs.org on May 13, 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 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 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.

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

Nonlinear wavefront control with all-dielectric metasurfaces Lei Wang†, Sergey Kruk†,*,Kirill Koshelev†,§, Ivan Kravchenko‡, Barry Luther-Davies⊥, and Yuri Kivshar †,§ †

Nonlinear Physics Centre, Australian National University, Canberra ACT 2601, Australia

§

ITMO University, St. Petersburg 197101, Russia

‡

Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, TN 37831

⊥Laser

Physics Centre, Australian National University, Canberra ACT 2601, Australia

KEYWORDS Metasurfaces, Nonlinear Optics, Third Harmonic Generation, Mie Resonances, Wavefront Control

ABSTRACT Metasurfaces, two-dimensional lattices of nanoscale resonators, offer unique opportunities for functional flat optics and allow to control transmission, reflection, and polarization of a wavefront of light. Recently, all-dielectric metasurfaces reached remarkable efficiencies, often matching or outperforming conventional optical elements. Exploitation of nonlinear optical response of metasurfaces offers a paradigm shift in nonlinear optics, and dielectric nonlinear metasurfaces are expected to enrich subwavelength photonics by enhancing substantially nonlinear response of natural materials combined with efficient control of the phase

ACS Paragon Plus Environment

1

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

Page 2 of 26

of nonlinear waves. Here we suggest a novel and rather general approach for engineering the wavefront of parametric waves of arbitrary complexity generated by a nonlinear metasurface. We design all-dielectric nonlinear metasurfaces and achieve a highly-efficient wavefront control of a third harmonic field and demonstrate the generation of nonlinear beams at a designed angle and nonlinear focusing of vortex beams. Our nonlinear metasurfaces produce phase gradients over a full 0-2π phase range with 92% diffraction efficiency.

ToC


2

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

In contrast to conventional optics, metasurfaces are composed of subwavelength elements (“metaatoms’) that can manipulate the light–matter interaction by employing a very compact platform. Metasurfaces with spatially varying meta-atoms allow to control polarization, phase, and amplitude of light [1]–[4]. By now, many phenomena have been demonstrated successfully in linear optics, however, nonlinear metasurfaces are expected to meet the growing demand for tailored nonlinearities [5]–[7] and realize nonlinear optical chirality, nonlinear geometric Berry phase, and nonlinear wavefront engineering. Several recent pioneering demonstrations of nonlinear metasurfaces for wavefront control [8]–[15] were based on plasmonic designs. One of the main limitation of plasmonic designs for the nonlinear metasurfaces is the nonlinear conversion efficiency. Different design strategies were implemented to increase nonlinear performance of plasmonic and hybrid metasurfaces [19]–[21] based on resonant coupling of meta-atoms. However, the overall efficiency of the nonlinear frequency conversion in such planar plasmonic nanostructures remained very small, being of the order of ∼10−10. Despite extremely high intrinsic nonlinearities of plasmonic materials, the experimental efficiencies were limited by Ohmic losses, small mode volumes and low laser damage thresholds. In addition, all previous approaches do not provide a direct mechanism to control the forward-to-backward ratio of generated nonlinear signals. All-dielectric nanostructures supporting Mie resonances have recently been suggested as an important pathway to enhance the nonlinear efficiency beyond the limits associated with plasmonics [22]. Indeed, more efficient third-harmonic generation in individual Si and Ge nanoresonators has been recently demonstrated by several groups, showing a huge enhancement of the conversion efficiency by optical pumping in the vicinity of the magnetic dipole Mie mode or composite resonances [23]–[25]. Si and Ge are of a paramount interest for nonlinear photonics


3

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

Page 4 of 26

as they offer superior nonlinear characteristics while allowing for a potentially low-cost CMOScompatible fabrication [26]. Conversion efficiencies of the order of 10−6 have been achieved experimentally in both individual Si and Ge nanoparticles and metasurfaces [27]. These results clearly illustrate the great potential of all-dielectric resonant nanostructured surfaces for nonlinear nanophotonics and meta-optics. Here we suggest a new concept for embedding any functionality into a nonlinear metasurface, and demonstrate our approach for silicon-based metasurfaces that generate the third-harmonic (TH) field. Constituent elements of our metasurfaces are sets of different elliptical nanopillar resonators, each of them providing similar nonlinear conversion efficiency but accumulating a different phase for the third-harmonic signal covering the full 0-2π range. The resonant phase accumulation of the third harmonic field in the dielectric metasurface occurs in the regime when both electric and magnetic Mie multipoles are excited simultaneously with comparable amplitudes. The simultaneous generation of several Mie multipoles allows to control the directionality of the third-harmonic signal in each individual nanopillar, and in particular to achieve high forward-tobackward ratio of the TH signal. When assembled into a planar structure, these dielectric resonators allow the creation of smooth phase gradients in the generated third-harmonic field. Here we exemplarily present two demonstrations of such metasurfaces. The first metasurface acts as a nonlinear beam deflector that generates a third-harmonic beam at a selected angle with respect to the direction of the pump beam [as visualized in Fig. 1 (a)]. The second metasurface generates a nonlinear focusing vortex beam with orbital angular momentum 1 from a Gaussian pump [as shown in Fig. 1 (b)].


4

Page 5 of 26 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 1. Concept images of functional nonlinear metasurfaces. (a) A nonlinear beam deflector. (b) A nonlinear vortex beam generator. The metasurfaces are assembled from a set of different silicon nanopillars generating third harmonic with different phases. Theoretical approach. We employ a set of silicon nanopillars as building blocks for our nonlinear metasurfaces. We require each nanopillar to support third-harmonic generation with forward directionality. We further require that the different pillars provide similar amplitudes of the third-harmonic signal but have different phases ranging from 0 to 2. This could be achieved via the interplay between the Mie resonances of a single nanopillar at both the pump wavelength and the TH wavelength. Multipolar analysis based on Mie theory proved to be a powerful tool to understand the light scattering from resonant nanoscale objects. In the past it yielded important insights into various fields of nanophotonics. In the past this framework was employed for describes the scattering intensities and radiation pattern of dominant excited multipole modes that yielded important insights into various fields of nanophotonics. [22], [28]-[30]. We extend this approach for the wavefront control in nonlinear regime (see Supporting Information). We consider pillars that at the pump wavelength support electric and magnetic dipole Mie resonances. Correspondingly, at the TH wavelength they support several higher-order multipoles as their


5

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

Page 6 of 26

optical volume increases. This allows an analytical treatment for the intensity of the third harmonic in both the forward (F) and backward (B) directions in terms of multipoles having even or odd parity with respect to the vertical reflection symmetry: (3ω)

(3ω)

(1)

(3ω)

(3ω)

(2)

B = |Eeven − Eodd |2 , F = |Eeven + Eodd |2 , (3ω)

where Eeven(odd) is the TH electric field generated by all even (odd) multipoles. Equation (1) suggests that the necessary condition for forward directionality requires the balance between the even and odd multipoles. (3ω)

(3ω)

Eeven = Eodd = |E|eiφ ,

(3)

Therefore, the total field of the third harmonic in the forward direction E (3ω) is (3ω)

(3ω)

E (3ω) = Eeven + Eodd = 2|E|eiφ .

(4)

By shaping the pillar geometry, a full 2π phase coverage with a uniform amplitude can be achieved keeping high forward directionality at the same time. Importantly, as the nanocylinders are optically resonant at both the pump wavelength, and the TH wavelength, the multipolar resonances at the TH are inevitably of higher order multipoles making the decomposition to go beyond the dipole and quadrupole approximations. The balance of multipoles with opposite parities required by Eq. (3) leading to the unidirectional pattern represents the generalized Huygens’ condition [4], [31]-[32] extended to the case of nonlinear harmonic generation. To find a specific set of geometries, we perform full-wave nonlinear simulations of amorphous silicon particles on glass substrate using the finite element method solver of COMSOL


6

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

Multiphysics. We choose a specific excitation wavelength of 1615 nm, and an initial pillar geometry that supports electric and magnetic dipole modes at the excitation wavelength. We note that calculations of third-harmonic are computer-intensive as they require at least a 33 times finer mesh compared to linear calculations. Moreover, our calculations are performed over a 4dimensional parameter space (comprising pillar height, two elliptical axes, and the unit cell size), resulting in large total volume of the 4D parameter space. To minimize the demand for computing power we implemented the gradient descent method in the COMSOL calculations. This method analyses the calculated data sets and learns to find the quickest route to optimization and is commonly used in machine intelligence, in particular, as a subset of the deep-learning method. The optimization was set to find the parameters that provide the exact desired values of TH phase near the maxima of TH amplitude and F/B ratio. Details of the numerical simulations are provided in Methods. The resulting design provides a set of particles with a height of 617 nm in a square array with 550 nm period, and elliptical cross-section with the axes values ranging between 320 nm and 535 nm [see Fig. 2 for the details of the nanopillar geometry]. This range of parameters allows a full 2π phase coverage in the TH to be achieved whilst maintaining similar amplitudes. All the nanopillars have a high forward-to-backward ratio of TH: overall over 90% of the third harmonic is generated in forward direction, and the TH generation in backward direction is suppressed [see Fig. 2 for the details on individual nanopillars]. The optimization is performed for a linearly polarized pump along a-axis of the pillars (as marked in Fig. 2). We note the possibility of independent phase and amplitude control for different polarizations. We numerically estimate the TH generation efficiency of the pillars. For the average pump power of 200mW, a pulse length of 300 fs, and a repetition rate 20MHz, and a diffraction-limited focusing on a single nanopillar, the calculated TH average power is ~300nW. Thus, the overall


7

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

Page 8 of 26

𝑇𝐻 ⁄ 𝑝𝑢𝑚𝑝 conversion efficiency 𝑃𝑎𝑣𝑒 𝑃𝑎𝑣𝑒 = 10−6 which is similar to previous reports on TH efficiency

in Si nanostructures in the vicinity of a magnetic dipole modes and composite resonances [33]3

𝑝𝑢𝑚𝑝 𝑇𝐻 [34]. The peak power-independent conversion efficiency 𝑃𝑎𝑣𝑒 ⁄(𝑃𝑎𝑣𝑒 ) = 10−4 𝑊 −2.

Figure 2. Geometries, nonlinear phases, and TH amplitudes of nanopillar meta-atoms. Shown are sizes of the nanopillars, and corresponding analytical and numerical results for the


8

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

phase, amplitude (in the forward direction), and directionality (the forward-to-backward ratio of intensity) of the third-harmonic field. The results are for the optimal pump wavelength of 1615 nm and linear polarization of the pump along a-axis. We compare our numerical results with the multipolar analysis [see Fig. 2] and find that calculations of phases of the TH are in a perfect agreement. Calculations of the amplitudes and the forward/backward ratios have a small offset that we associate with the effect of the substrate [see details in Supporting Information]. Our multipolar analysis shows that at the pump wavelength the forward-scattering of all the pillars is dominated by the 1st order multipoles: electric and magnetic dipoles, and correspondingly the third-harmonic field is dominated by several low-order multipoles: electric and magnetic dipoles, quadrupoles, and octupoles [see details in Supporting Information]. Importantly, the sum of contributions of all even multipoles to the TH signal is similar to the sum of contributions of all odd multipoles, [in terms of both amplitudes and phases, see details in the Supporting Information]. Thus, our system satisfies Eqs. (3), (4) and achieves the generalized Huygens’ condition in the nonlinear regime, that provides a full 2π phase coverage combined with the similar TH amplitude in the forward direction, and supressed reflection in backward direction for all nanopillars (see Fig. 2). Full-phase coverage combined with nearly the same intensity of forward-generated TH signal allows for an efficient arbitrary complex wavefront shaping in the nonlinear regime. We finally perform the multipolar analysis over the spectral range of pump wavelengths from 1300 nm to 1800 nm with the details provided in Supporting Information.


9

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

Page 10 of 26

Figure 3. Nonlinear beam deflector metasurface. (a) Phase profile of the third-harmonic field encoded into the metasurface. (b) Electron microscope image of the silicon metasurface. (c) Directionality diagram (back-focal plane image) of the forward TH field. 92% of TH is directed into the designed angle θ=5.6o, where kx/k0=-0.098. (d) Cross section of experimental back-focal plane image normalized to the TH maximum. Experimental demonstrations. To demonstrate the capabilities for nonlinear wavefront control, we fabricate two metasurfaces: a nonlinear beam deflector – a metasurface that generates a beam of light at the third-harmonic wavelength propagating at an angle to the k-vector of the pump; and a nonlinear focusing vortex beam generator – a metasurface that generates a focusing TH beam with an orbital angular momentum m=1. The first metasurface acting as a nonlinear beam deflector serves as a canonical example of wavefront control similar to common demonstrations in linear meta-optics [35]–[40]. The second metasurface provides an example of a more complex wavefront control revealing the potential of our general approach. For the beam deflector metasurface we arrange the 10-pillar set from Fig. 2 into a supercell that creates a linear phase gradient at the TH wavelength. The supercell is then arranged into a periodic array [see Fig. 3 (a)]. Such a metasurface generates the TH mainly in the direction determined by


10

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

the condition |kx/k0|=λTH/(10*period)=0.098, where k0 and kx are the full length and in-plane projection of the TH wavevector, and λTH is the TH wavelength.

Figure 4. Spectral performance of nonlinear metasurface. (a) Theoretically calculated and (b) experimentally measured directionality diagrams resolved spectrally. Images show true colours of the TH field for different wavelengths. (c) Diffraction efficiency spectrum defined as the percentage of forward-generated third harmonic directed into the -1st order. Shown are theoretical (solid) and experimental (dots) results. We fabricate the nonlinear beam deflector from amorphous silicon on a glass substrate using electron beam lithography [for details, see Methods]. An electron microscope image of the fabricated structure is shown in Fig. 3 (b). To generate the TH signal, we pump the metasurface at a wavelength of 1615 nm with femtosecond pulses from an optical parametric amplifier pumped


11

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

Page 12 of 26

by a mode-locked Ytterbium laser. To observe the directionality of the TH and to measure the beam deflection we image the back-focal plane onto a camera [see details in Methods]. Nearlyperfect TH beam deflection is observed at an angle of θ=5.6o (where sin[θ]=|kx/k0|) relative to the direction of pump with 92% diffraction efficiency in forward direction. Importantly, virtually zero TH is generated along the direction of pump (i.e. into the zero diffraction order). Considering high forward-to-backward ratio of TH generation, we conclude that over 80% of overall third harmonic energy is generated into the desired channel. The experimentally measured TH conversion 3

𝑝𝑢𝑚𝑝 𝑇𝐻 efficiency 𝑃𝑎𝑣𝑒 ⁄(𝑃𝑎𝑣𝑒 ) = 1.4 × 10−4 𝑊 −2 .

We then studied the performance of the metasurface as the pump wavelength is tuned from 1330 nm to 1750 nm, the range available from our optical parametric amplifier. Directionality curves, similar to Fig. 3(b), are experimentally retrieved as the pump wavelength was varied, and then combined into the diagram in Figure 4(b). The colour space (HSL colour-map) is used to represent the TH directionality and wavelength information where the hue (H) shows true colour of TH for given wavelength, saturation (S) is 1, and lightness (L) corresponds to the normalized TH intensity at each wavelength. The complimentary theoretical calculation is shown in Fig. 4(a). The calculation took both the amplitudes and phases of the TH into account across the spectral range. We can see that around the optimal wavelength of 1615 nm the only diffraction order clearly visible is the -1st order. As the pump wavelength deviated from its optimal value, other diffraction orders appear in the directionality diagram. From this we could extract the variation of the diffraction efficiency with wavelength [see Fig. 4(c)]. We observe good correspondence between theoretically calculated and experimentally measured diffraction efficiencies of the wavefront control.


12

Page 13 of 26 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. Nonlinear vortex beam generator. (a) Phase profile of the third-harmonic field encoded into the metasurface. (b) Optical microscope image of the fabricated metasurface. (c) Cross-section of a generated donut-shape vortex beam taken perpendicular to the optical axis at distance z=25 µm from the focus. (d) Cross-section along the optical axis of the focusing and defocusing TH beam. For the vortex metasurface, we design a phase profile shown in Fig. 5 (a). We assemble the metasurface from an increased set of pillars found numerically [see Fig. 5 (b) for the optical image of the fabricated sample]. We pump the metasurface with a 1615 nm wavelength, and scan the formed TH beam along the optical axis. Figure 5 (c) shows the cross-section of a donut-shape vortex TH beam taken perpendicular to the optical axis at 25 µm distance from the focus. Figure 5(c) shows the cross-section along the optical axis of the focusing and de-focusing TH beam. The nonlinear vortex beam generation illustrates the immediate applicability of our approach for engineering complex phase fronts of parametric waves.


13

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

Page 14 of 26

In summary, we have presented a general approach and versatile experimental platform for nonlinear wavefront control with highly-efficient nonlinear dielectric metasurfaces. Our approach is based on the generalized Huygens’ principle extended to nonlinear optics. This allows creating arbitrary phase gradients and wavefronts in nonlinear optics via multipolar nanophotonics, by excitation of electric and magnetic Mie multipoles. The multipolar generation also allows to achieve high forward-to-backward ratio of nonlinear signals. Based on our concept, we have designed and demonstrated experimentally the first nonlinear all-dielectric metasurface that generates the third harmonic signal in forward direction with 92% precision in its wavefront control. METHODS Sample fabrication. The silicon metasurfaces are fabricated on a 4 inch fused silica wafer 500 µm thick. First, a 617 nm polycrystalline silicon layer is deposited onto the substrate by lowpressure chemical vapor deposition (LPCVD). Subsequently, a thin layer of an electron-resist PMMA A4 950 is spin-coated onto the sample, followed by electron-beam lithography (JEOL 100 eV) and development. After this, a 20 nm Cr film is evaporated onto the sample, followed by the lift-off process to generate a Cr mask. Reactive-ion etching (RIE) is used to transfer the Cr mask pattern into the silicon film. Finally, the residual Cr mask is removed via wet Cr etching. Finally, the sample is cleaned with oxygen plasma. Nonlinear optics setup. For nonlinear optical measurements, we pump the metasurfaces with an optical parametric amplifier (MIROPA-fs-M from Hotlight Systems) tunable across the range from 1330-1750 nm and generating 300 fs duration pulses at a repetition rate of 21 MHz and average power at the laser output ranging from 100-300mW. The power is then attenuated to 9


14

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

mW. The metasurface is mounted on a three-dimensional stage, and is placed at normal incidence in the focal region of two confocal Olympus Plan N objective lenses. A low-NA lens is used to focus the pump (NA=0.1, 4x). The TH is collected using a lens with NA= 0.4. The beam is focused from the front side of the sample. To detect the TH radiation, we use a cooled CCD camera (Starlight Xpress Ltd, SXVR-H9). The TH signal was separated from the pump using shortpass filter Thorlabs FESH0650 for spectra measurements. A pair of confocal lenses are used to transfer the back-focal plane image of the TH radiation onto the camera. The spectrum of the pump beam is cleared up by a dichroic mirror and long-pass infrared filters (Thorlabs FEL 1000 and FELH1150) before the first objective. Numerical calculations. For numerical simulations of the nonlinear amplitudes and phases of the forward-scattered TH wave, we use the finite-element-method solver in COMSOL Multiphysics in the frequency domain. All calculations are realized for a single nanopillar of a specific size placed on a semi-infinite substrate with periodic boundary conditions mimicking a square array. Thus, the simulation takes into account the effect of coupling between identical neighbour particles but neglects the difference in coupling strengths between non-identical near-neighbours. We employ the approach based on an undepleted pump approximation, using two steps to calculate the radiated nonlinear emission. First, we simulate the linear scattering at the fundamental wavelength, and then obtain the nonlinear polarization induced inside the nanopillar. Then, we employ this as a source for the electromagnetic simulation at the harmonic wavelength to obtain the generated TH field. The nonlinear susceptibility tensor χ(3) was considered as a constant scalar value 2.45 × 10−19 m2 /V 2 . ASSOCIATED CONTENT


15

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

Page 16 of 26

Supporting Information. The following files are available free of charge. Details of the analytical approach. (PDF)

AUTHOR INFORMATION Corresponding Author *e-mail: [email protected] Author Contributions L.W. and S.K. conceived the research. L.W and K.K. performed numerical simulations and theoretical calculations. S.K. and I.K. fabricated the samples. L.W., S.K., and B.L.D. provided experimental data. S.K., K.K. and Y.K. wrote the manuscript based on the input from all authors. Y.K. supervised the project. All authors contributed to the editing of the manuscript. Funding Sources The authors acknowledge a financial support from the Ministry of Education and Science of the Russian Federation (3.1668.2017), the Australian Research Council and the Australian National University. A part of this research was conducted at the Center for Nanophase Materials Sciences which is a DOE Office of Science User Facility. K.K. acknowledges FASIE (10864GU/2016) for the valuable support. Notes The authors declare no competing financial interest. ACKNOWLEDGMENT


16

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

The authors acknowledge useful discussions with H. Atwater, M. Brongersma, D. Christodoulides, A. Kuznetsov, W. Liu, S. Zhang, T. Zentgraf, A. Zayats, and also thank D. Smirnova for a help with the numerical code and advices for numerical simulations.


17

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

Page 18 of 26

REFERENCES [1]

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

[2]

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

[3]

Arbabi, A.; Horie, Y.; Bagheri, M.; Faraon, A. Nat. Nanotechnol. 2015, 10, 937–943.

[4]

Kruk, S.; Hopkins, B.; Kravchenko, I. I.; Miroshnichenko, A.; Neshev, D. N.; Kivshar, Yu.

S. APL Photonics 2016, 1, 030801. [5]

Zheludev, N. I.; Kivshar, Yu. S. Nat. Mater. 2012, 11, 917–924.

[6]

Li, G.; Zhang, S.; Zentgraf, T. Nat. Rev. Mater. 2017, 2, 17010.

[7]

Krasnok, A.; Tymchenko, M.; Alù, A. Mater. Today 2017, 21, 8-21.

[8]

Ellenbogen, T.; Voloch-Bloch, N.; Ganany-Padowicz, A.; Arie, A. Nat. Photonics 2009,

3, 395-398. [9]

Li, G.; Chen, S.; Pholchai, N.; Reineke, B.; Wong, P. W. H.; Pun, E. Y. B., Cheah, K. W.;

Zentgraf, T.; Zhang, S., Nat. Mater. 2015, 14, 607–612. [10]

Tymchenko, M.; Gomez-Diaz, J. S.; Lee, J.; Nookala, N.; Belkin, M. A.; Alù, A. Phys.

Rev. Lett. 2015, 115, 207403. [11]

Wolf, O.; Campione, S.; Benz, A.; Ravikumar, A. P.; Liu, S.; Luk, T. S.; Kadlec, E. A.;

Shaner, E. A.; Klem, J. F.; Sinclair, M. B.; Brener, I. Nat. Commun. 2015, 6, 7667. [12]

Almeida, E.; Bitton, O.; Prior, Y.; Nat. Commun. 2016, 7, p. 12533.

[13]

Almeida, E.; Shalem, G.; Prior, Y.; Nat. Commun. 2016, 7, 10367.

[14]

Keren-Zur, S.; Avayu, O.; Michaeli, L.; Ellenbogen, T. ACS Photonics 2016, 3, 117–123.

[15]

Ye, W.; Zeuner, F.; Li, X.; Reineke, B.; He, S.; Qiu, C.-W.; Liu, J.; Wang, Y.; Zhang, S.;

Zentgraf, T. Nat. Commun. 2016, 7, 11930. [16]

Walter, F.; Li, G.; Meier, C.; Zhang, S.; Zentgraf, T; Nano Lett. 2017, 17, 3171-3175.


18

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

[17]

Schlickriede, C.; Waterman, N.; Reineke, B.; Georgi, G.; Zhang, S.; Zentgraf, T. Adv. Opt.

Mater. 2018, 1703843. [18]

Schlickriede, C.; Waterman, N.; Reineke, B.; Georgi, P.; Li, G.; Zhang, S.; Zentgraf, T.

Adv. Mater. 2018, 30, 1703843. [19]

Kauranen, M.; Zayats, A. V. Nat. Photonics 2012, 6, 737–748.

[20]

Celebrano, M.; Wu, X.; Baselli, M.; Großmann, S.; Biagioni, P.; Locatelli, A.; De Angelis,

C.; Cerullo, G.; Osellame, R.; Hecht, B; Duò, L.; Ciccacci, F.; Finazzi, M. Nat. Nanotechnol. 2015, 10, 412–417. [21]

Lee, J.; Tymchenko, M.; Argyropoulos, C.; Chen, P.-Y.; Lu, F.; Demmerle, F.; Boehm, G.;

Amann, M.-C.; Alù, A.; Belkin, M. A. Nature 2014, 511, 65–69. [22]

Kruk, S.; Kivshar, Yu. ACS Photonics 2017, 4, 2638-2649.

[23]

Smirnova, D; Kivshar, Yu. S. Optica 2016, 3, 1241–1255.

[24]

Grinblat, G.; Li, Y.; Nielsen, M. P.; Oulton, R. F.; Maier, S. A. Nano Lett. 2016, 16, 4635–

4640. [25]

Wang, L.; Kruk, S.; Xu, L.; Rahmani, M.; Smirnova, D.; Solntsev, A.; Kravchenko, I.;

Neshev. D.; Kivshar, Yu. Nanoscale 2017, 9, 2201-2206. [26]

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

[27]

Yang, Y.; Wang, W.; Boulesbaa, A.; Kravchenko, I. I.; Briggs, D. P.; Puretzky, A.;

Geohegan, D.; Valentine, J. Nano Lett. 2015, 15, 7388–7393. [28]

Liu, W.; Zhang, J.; Lei, B.; Ma, H.; Xie, W.; Hu, H. Opt. Express 2014, 22, 16178-16187.

[29]

Kuznetsov, A. I.; Miroshnichenko, A. E.; Brongersma, M. L.; Kivshar, Y. S.; Luk’yanchuk,

B. Science 2016, 354, aag2472. [30]

Kivshar, Y.; Miroshnichenko, A. Opt. Photonics News 2017, 28, 24-31.


19

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

[31]

Page 20 of 26

Decker, M.; Staude, I.; Falkner, M.; Dominguez, J.; Neshev, D. N.; Brener, I.; Pertsch, T.;

Kivshar, Yu. S. Adv. Opt. Mater. 2015, 3, 813–820. [32]

Wang, L.; Kruk, S.; Tang, H.; Li, T.; Kravchenko, I.; Neshev, D. N.; Kivshar, Yu. S.

Optica 2016, 3, 1504. [33]

Shcherbakov, M. R.; Neshev, D. N.; Hopkins, B.; Shorokhov, A. S.; Staude, I.; Melik-

Gaykazyan, E. V.; Decker, M.; Ezhov, A. A.; Miroshnichenko, A. E.; Brener, I.; Fedyanin, A. A. Nano Lett. 2014, 14, 6488-6492. [34]

Yang, Y.; Wang, W.; Boulesbaa, A.; Kravchenko, I. I.; Briggs, D. P.; Puretzky, A.;

Geohegan, D.; Valentine, J. Nano Lett. 2015, 15, 7388-7393. [35]

Yu, Y. F.; Zhu, A. Y.; Paniagua‐Domínguez, R.; Fu, Y. H.; Luk’yanchuk, B.; Kuznetsov,

A. I. Laser Photon. Rev. 2015, 9, 412–418. [36]

Lin, D.; Melli, M.; Poliakov, E.; Hilaire, P. S.; Dhuey, S.; Peroz, C.; Cabrini, S.;

Brongersma, M.; Klug, M. Sci. Rep. 2017, 7, 2286. [37]

Zhou, Z.; Li, J.; Su, R.; Yao, B.; Fang, H.; Li, K.; Zhou, L.; Liu, J.; Stellinga, D.; Reardon,

C. P.; Krauss, T. F.; Wang, X. ACS Photonics 2017, 4, 544–551. [38]

Chalabi, H.; Ra’di, Y.; Sounas, D. L.; Alù, A.; Phys. Rev. B 2017, 96, 75432.

[39]

Khaidarov, E.; Hao, H.; Paniagua-Domínguez, R.; Yu, Y. F.; Fu, Y. H.; Valuckas, V.; Yap,

S. L. K.; Toh, Y. T.; Ng, J. S. K.; Kuznetsov, A. I. Nano Lett. 2017, 17, 6267–6272. [40]

Sell, D.; Yang, J.; Doshay, S.; Yang, R.; Fan, J. A. Nano Lett. 2017, 17, 3752–3757.


20

-1st

Page 21 of 26 1 2 3 4 5 6 7 8 9 10 11 12

th Nano Letters

0

3ω


ω

1st

a

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

Nano Letters

b

3ω

ω


Page 22 of 26

Geometry

b 26 Page 23 of

1 2 3 a (nm) 4 5 b (nm) 6 7 2π 8 9 10 11 12 π 13 14 15 16 17 0 18 19 20ave 21ave 22 23 24 25 26 27 28 0 29 100 30 31 32 33 34 35 10 36 37 38 39 40 1

495

410

450

350

340

340

320

525

450

535

390

505

380

505

470

400

450

460

535

360

numerical simulations multipolar analysis


Forw./Backw. Ratio

TH Amplitude (a.u.)

THG Phase (rad)

Nano Letters

a



a

c

-3rd

-2nd

st

-1Letters Nano

1st

0th

2nd

3rd

max

Page 24 of 26

10μm

0

b

5µm

0

THG Intensity (a.u.)

1 2 3 4 5 6 7 8 9 10 11 12 13

2π

d

1

0.5

0 -0.4

-0.3

ACS Paragon Plus Environment -0.2

-0.1

0

kx/ko

10μm

0.1

0.2

0.3

0.4

Page 25 of 26 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


a 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

b

Theory

c

Bright eld

Nano Letters

10μm

10μm

y, µm

d


z, µm

TH Beam ProlePage 26 of 26

10μm

Nonlinear Wavefront Control with All-Dielectric Metasurfaces - Nano

Recommend Documents