Plasmonic metasurfaces, which can be considered as the two-dimensional analog of metal-based metamaterials, have attracted progressively increasing attention in recent years because of the ease of fabrication and unprecedented control over the reflected or transmitted light while featuring relatively low losses even at optical wavelengths. Among all the different design approaches, gap-surface plasmon metasurfaces – a specific branch of plasmonic metasurfaces – which consist of a subwavelength thin dielectric spacer sandwiched between an optically thick metal film and arrays of metal subwavelength elements arranged in a strictly or quasi-periodic fashion, have gained awareness from researchers working at practically any frequency regime as its realization only requires a single lithographic step, yet with the possibility to fully control the amplitude, phase, and polarization of the reflected light. In this paper, we review the fundamentals, recent developments, and opportunities of gap-surface plasmon metasurfaces. Starting with introducing the concept of gap-surface plasmon metasurfaces, we present three typical gap-surface plasmon resonators, introduce generalized Snell’s law, and explain the concept of Pancharatnam-Berry phase. We then overview the main applications of gap-surface plasmon metasurfaces, including beam-steerers, flat lenses, holograms, absorbers, color printing, polarization control, surface wave couplers, and dynamically reconfigurable metasurfaces. The review is ended with a short summary and outlook on possible future developments.
The capability of manipulating light at will is tantalizingly attractive, promising numerous applications. Conventional methods to mold the flow of light are typically relying on gradually accumulated amplitude, phase, and polarization changes during light propagation , with the resulting devices featuring certain curved surfaces and complex shapes. These bulky configurations do not comply with current trends aiming at very dense integration and miniaturization in photonics/plasmonics.
In recent years, metasurfaces, i.e. thin planar arrays of resonant subwavelength elements arranged in a periodic or aperiodic manner, have attracted progressively increasing attention and become a rapidly growing field of research because of their unprecedented control over electromagnetic (EM) fields by modifying boundary conditions for impinging waves , . Generally, metasurfaces significantly benefit from their simplified fabrication process compared with their three-dimensional (3D) analog – metamaterials that are usually made up of complicated 3D artificial nanostructures processing exceptional EM properties , , . As such, metasurfaces are promising for on-chip integration in the area of photonics and plasmonics owing to their planar profiles. Among all the metasurfaces, gap-surface plasmon metasurfaces (GSPMs) – a specific branch of plasmonic metasurfaces – which consist of a subwavelength thin dielectric spacer sandwiched between an optically thick metal film and arrays of metal subwavelength elements arranged in a strictly or quasi-periodic fashion, have gained awareness from researchers working at practically any frequency regime as its realization only requires a single lithographic step, yet with the possibility to fully control the amplitude, phase, and polarization of the reflected light.
As metasurfaces comprise a rapidly growing field of research, several excellent review articles can be found, some presenting a comprehensive review of the filed , , , , , , , , , , , while others are mainly focused on a certain type of metasurfaces , , , ,  or specific branches of applications , , , , . In this review, we review the recent progress of GSPMs during the past few years and attempt to provide our prospective on this specific branch of metasurfaces. The rest of this paper is organized in the following sections. In Section 2, we briefly discuss and classify the typical gap-surface plasmon (GSP) resonators and then introduce the concept of generalized Snell’s law and Pancharatnam-Berry (PB) phase. Section 3 is devoted to the applications of GSPMs, including beam-steerers, flat lenses, holograms, absorbers, color printing, polarization control, surface wave couplers, and dynamically reconfigurable metasurfaces. Finally, we summarize and provide perspective for future developments in Section 4.
Surface plasmon polaritons (SPPs) are EM excitations propagating along and tightly bound to metal-dielectric interfaces, which arise from the coupling of EM field in dielectric and collective electron oscillations in metal . Owing to the unique properties of subwavelength confinement and large field enhancement, SPPs have motivated numerous potential applications, ranging from biochemical sensing , quantum optics , and information communication technologies  to sustained energy sources . However, in the long-wavelength part of the visible spectrum, SPPs are quite close to the light line and unwieldy to be utilized to achieve strongly enhanced local field. To realize strong field enhancement at any wavelength, one effective method is to place two metal-dielectric interfaces close to each other, which introduces the near-field coupling (NFC) between the SPPs of individual interface . Among the multilayer structures, the metal-insulator-metal (MIM) configurations involving GSPs are superior to insulator-metal-insulator resonators in terms of high-quality factor and strong mode confinement , . Here we will briefly discuss three typical GSP resonators with MIM configuration and present their scattering properties.
We begin the discussion of GSP resonator by revisiting the retardation-based resonant nanostructure realized by truncating the infinitely long MIM configuration, which consists of two gold strips separated by a silicon dioxide (SiO2) spacer (left panel of Figure 1A). Because of the terminations of the MIM structure, multiple reflections of GSP modes can be achieved, resulting in lateral standing-wave GSP modes, whose resonance positions can be described by a simple Fabry-Perot (FP) resonator formula , :
where λ0 is the wavelength in free space, w is the strip width, nGSP is the real part of the effective mode index of GSP, m is an integer defining the order of the GSP mode, and ϕ is an additional phase shift acquired upon reflection of the GSPs at the boundaries of the MIM configuration. By properly tuning the dimensions according to equation (1), GSP resonance at ~800 nm can be achieved with a gap size of d=50 nm, which shows a dominant scattering process, making this structure a good candidate for reflective metasurfaces (middle panel of Figure 1A) . Additionally, the GSP mode at 800 nm shows a characteristic of magnetic dipole (MD) resonance with out-of-phase currents in the two metal strips and a maximum of the magnetic field in the gap (corresponding to a minimum of the electric field; right panel of Figure 1A) , , . It is worth pointing that when the gap size is small, the GSP resonator is mainly an absorbing element because of the strong NFC between two metallic layers , , . Once the gap size goes too large, the NFC becomes too week to support GSP resonance, and the resonance undergoes a transition from a lateral MD resonance to a vertical Fabry-Perot resonance .
Besides the highly symmetric configuration, the GSP resonator can be constructed with metal strip placed close to an infinitely extended metal surface, which shows similar scattering resonances related to the GSP mode . From Figure 1B, one can clearly see that the first-order GSP mode (i.e. MD resonance) has a red-shift of ~25 nm, which is ascribed to the less sharp termination of the resonator, allowing the GSP resonator to have a slight larger effective width . In addition, such resonator shows stronger scattering with a two-fold increase in the scattering cross section. One should note that such resonator can couple part of the scattered light into SPPs propagating along the metal-air interface.
As a third kind of GSP resonator, nanoscale truncation is carried out only for the top metal layer in a MIM structure, which can be termed as continuous-layer GSP resonator  (Figure 1C). In addition to the attractive properties of usual GSP resonators, this continuous-layer GSP resonator can be easily achieved by depositing gold in the resist profile with a single lithography step, thus relaxing the lift-off process. Compared with the second configuration, the GSP mode shows a strong red-shift and is excited at λ=915 nm, resulted from the removal of horizontal refractive index contrast between air and SiO2 spacer . Because of the simplified fabrication process, most of the GSMPs are constructed by the third kind of GSP resonators, in which the upper metal layer usually consists of periodic arranged antennas.
As a final comment, it ought to be mentioned that the GSP resonator could be well described by the coupled-mode theory , , . Base on coupled-mode theory, Zhou’s group has established a generic phase diagram where the optical properties of GSP resonator are fully controlled by the intrinsic and radiative losses, thereby greatly facilitating the design of appropriate GSMPs with tailored functionalities .
The abrupt phase shift, denoted as phase discontinuity, at the interface between two media was first introduced by Capasso’s group in 2011, which allows us to revisit the laws of reflection and refraction by applying Fermat’s principle, thereby realizing generalized Snell’s laws of reflection and refraction . Figure 2 schematically illustrates the one-dimensional (1D) system used to derive the generalized Snell’s laws, where the interface between two media is consisting of artificially metasurface in order to introduce a position-dependent phase shift Φ(x) in the light path. Considering an incident plane wave impinging at angle of θi relative to the surface normal of the metasurface-decorated interface, the generalized laws of reflection and refraction can be written as
where ni and nt are the refractive indexes of the two media, λ0 is the wavelength in free space, and θr and θt are the reflected and refracted angles, respectively. From equations (2) and (3), it is evident that the reflection/refracted beam can have an arbitrary direction, provided that a suitable, constant, and nonzero phase gradient along the interface (dΦ/dx) is introduced. In the case of dΦ/dx=0, we recover the usual laws of reflection of refraction that imply continuity of the in-plane wave vector.
To realize a metasurface with interface discontinuities, a V-shaped antenna, supporting symmetrical and antisymmetrical electric dipole resonance, was first proposed and experimentally demonstrated to control the reflection/refraction of linearly polarized (LP) light at middle-infrared range governed by the generalized Snell’s law . Soon after, the idea was extended to the near-infrared regime, and the functionality was found to be wideband . It should be noted that while great successes have already been achieved, there are two issues unsolved in the previous metasurfaces based on a V-shaped antenna , . In these cases, the 2π phase control is achieved with reflected/transmitted light polarized orthogonally to the incident wave as the V-shaped antenna only supports electric dipole resonance and thus allows one to control the phase in a limited range, from 0 to π, because of its Lorentz-like polarizabilities , . Such polarization conversion will cause inconvenience for some applications. Moreover, the theoretical upper bound of the coupling efficiency between the two polarizations for such single nonmagnetic metasurfaces is 25% in the limit of negligible absorption , . To solve these two issues, one can place the metal antennas in close proximity to a back metal reflector separated by an optically thin dielectric spacer, making the GSPMs, which ensure an almost full control of the phase space together with a high efficiency reaching 100% for the lossless system for small angles .
One final point ought to be mentioned is that the generalized Snell’s law is not enough for realizing perfectly performing GSPMs because of significant parasitic reflections in undesired directions, especially for large incident angles . To perfectly control reflection, one may use properly engineered bi-anisotropic meta-atoms ,  or strong nonlocal response .
In the previous examples, an abrupt phase discontinuity is introduced by varying meta-atoms’ geometric parameters; another completely different technique, known as PB phase or geometric phase, achieves a full phase control of the cross-polarized light by using anisotropic meta-atoms with identical geometry but spatially varying orientations for circularly polarized (CP) light , . Here it should be emphasized that the recent development of metasurfaces based on PB phase has been largely following the early work on space-variant subwavelength gratings by Hasman and co-workers , , , . The concept behind PB phase is simple, which can be easily revealed by Jones matrix . In general, the Jones matrix of an anisotropic GSP meta-atom that rotated an angle θ from the x-axis can be written as
where Mxx and Myy are the reflection or transmission coefficients for LP light along the two primary axes of the anisotropic meta-atom, respectively, and R(θ) is the rotation matrix. When the incident light is CP, the reflected or transmitted light can be written as
where represents the incident left-handed CP (LCP) and right-handed CP (RCP) light, respectively. The first term in equation (5) represents CP light with the same helicity as the incident wave, while the second term stands for CP light with the opposite helicity that gains an additional PB phase of ±2θ, whose sign depends on the helicity state. Thereby, a 2π phase coverage can be achieved if the meta-atom is rotated from 0 to 180°. In this sense, the PB phase can be regarded as dispersionless. However, the conversion efficiency may be strongly wavelength-dependent, as it relies on the ability of individual meta-atom to tune Mxx and Myy. It is clear that once the condition Myy=−Mxx is satisfied, the co-polarized light will vanish and the entire light is converted to the cross-polarization. As such, GSPMs play an important role in realizing highly efficient PB phase metasurface with the conversion efficiency up to 100% in the limit of negligible Ohmic loss .
GSPMs have been utilized in a wide range of applications, ranging from flat and compact versions of conventional optical elements with control over the amplitude, phase, and polarization of incident wave to dynamically reconfigurable GSPMs. Here we try to give an overview of the many, yet very different, applications.
As described in Section 2.2, beam steering is a fascinating functionality that can be realized by metasurfaces with constant phase gradient along the surface, leading to generalized Snell’s law of reflection and refraction that allows light to propagate in anomalous directions. Previous attempts to demonstrate beam steering typically relied on single-layered metasurfaces, which suffers from low efficiency , . In contrast, GSPMs enabling full control of phase with high efficiency are an excellent candidate to realize highly efficient beam steering for both the LP light and CP light.
The first experimental demonstration of GSPM for beam steering is shown in Figure 3A, where H-shaped antennas are on top of a dielectric layer coated metal film, designed to anomalously reflect x-polarized incident wave in the microwave regime . Through engineering the top H-shaped meta-atoms properly, a constant phase gradient is introduced along the x-axis, thereby resulting in diffraction into +1 order, verified by the angle-resolved far-field measurements (right panel of Figure 3A). As there is no polarization converted to the cross-polarized component and Ohmic loss is negligible, this GSPM shows close to 100% efficiency. Here it is worth noting that this type of GSPMs can be easily extended to near-infrared range with slightly reduced efficiency of ~80% because of increased Ohmic loss , .
Besides the high efficiency in controlling the reflection phase, GSPMs show even greater control of the reflected light by tailoring the reflection coefficients for two polarizations independently, thereby allowing for the design of birefringent GSPMs that achieve beam steering and splitting of individual LP light simultaneously , , , , . As an example, Figure 3B shows a MIM unit cell for a birefringent GSPM working at 800 nm . By properly selecting nanobricks with different widths, the incident light experiences an incremental (decremental) phase step within the supercell for the x-polarized (y-polarized) light, resulting in anomalous diffraction into +1 (−1) diffraction order (right panel of Figure 3B). The proof-of-concept experiment shows a reduced efficiency of ~50% compared with the theoretical efficiency around 80%, which is ascribed to the fabrication imperfections and uncertainties in the material properties. In principle, the polarization splitter functionality can be easily extended to visible by scaling the dimensions of the unit cell, which however requires small sizes and high accuracy of rectangular antennas. To relax these fabrication constraints, high-order GSP resonance was introduced to realize polarization beam splitter with excellent quality in terms of angular-independence and high polarization contrast .
Although the high-order GSPM could relax the difficulty in fabrication, it is rather challenging to experimentally realize metasurfaces with high-performance that operate in the visible range because of difficulties in fabricating several meta-atoms of subwavelength size with high uniformity. As an alternative, continuous rather than discrete phase delay using single spatial-variant gradient meta-atom is proposed to achieve beam steering , , . Here we highlight a broadband and highly efficient anomalous meta-reflectarray across the visible frequency range by using trapezoid-shaped nanorod as the building block, which is quite easy to fabricate with high accuracy (Figure 3C) . The broadband anomalous reflection and spectrum splitting have been experimentally confirmed at visible and near-infrared frequencies, shown in the right panel of Figure 3C. Specifically, most of the incoming power is reflected toward the anomalous direction (+1 order diffraction), while other diffraction modes are significantly suppressed. The average power ratio of +1 order diffraction to other diffractions for visible bandwidth from 500 to 800 nm is on the order of 103.
The resonant phase realized by tailoring the geometry dimensions of meta-atoms is intrinsically dispersive, and thus, the operating bandwidth is restricted. Distinct from the resonant phase, the PB phase is dispersionless in the sense that the induced phase on the cross-polarized light is based on identical meta-atom with spatially varying orientations as discussed in Section 2.3. Within the concept of GSPMs, Figure 3D depicts a highly efficient PB phase metasurface working at near-infrared regime, which consists of identical split rings with different opening orientations on the surface of a SiO2-silver bilayer . The right panel of Figure 3D shows the measured intensity of LCP and RCP components of two anomalous beams when the LP light impinges normally on the surface, confirming that the anomalous beam +1 is RCP and the anomalous beam −1 is LCP. Because of the high polarization conversion ratio and broad bandwidth realized by the rotating split ring, this structure possesses significantly strong anomalous reflection beams (over 70% of incident light intensity) over a broad frequency range (1100–1750 nm). We note that a similar GSPM design has been realized in the microwave range with a record-high conversion efficiency of >90% .
In conventional optical lenses, the wavefront engineering is typically based on the propagation phase gradually accumulated along the optical path. In this way, the dielectric lenses or parabolic reflectors with particular surface topography are bulky and curved components, thereby preventing the trend of very dense integration and miniaturization in photonics. Metasurfaces, on the other hand, enable an unprecedented opportunity to engineer the wavefronts at will, because of their planar profiles and relative ease of fabrication. To realize flat lenses, metasurfaces should impose the following parabolic phase profile to convert the incident planar wavefronts into spherical ones in transmission or reflection:
where the metasurface is assumed to lie in the xy-plane, f is the focal length, and λ0 is the wavelength in free space. Here we review some of planar metalenses based on GSPMs, while also addressing the important topic of achromatic lenses. Starting with polarization-sensitive GSP metalens for 1D focusing , , Figure 4A conveniently illustrates the scanning electron microscope (SEM) image of the fabricated metalens, consisting of an array of gold nanobricks atop a thin dielectric spacer and continuous gold film. The dimensions of the meta-atoms are varied along the horizontal direction so as to achieve a parabolic phase profile. Here the continuously varying phase profile is approximated by a step-like function, where the phase takes on a finite number of values, thus making the metasurfaces consist of a series of aperiodic meta-atoms. The fabricated metalens can focus an LP (along the direction of nanobrick size variation) incident beam in the plane of its polarization with the focal length, which changes from ~15 to 11 μm when tuning the light wavelength from 750 to 950 nm, respectively, verifying the broadband focusing ability. As shown in the lower panel of Figure 4A, the highest efficiency (~27%) was measured at the wavelength of 800 nm. Even though the measured efficiency is considerably smaller than the calculated one because of various fabrication-related circumstances, it is substantially larger than the upper limit of ~10% estimated for V-antenna metasurfaces , . We emphasize that polarization-independent focusing can be realized using symmetrical meta-atoms.
As stated in the previous section, birefringent metasurfaces have even greater capability of controlling the reflected light with different polarizations independently by engineering the reflection amplitudes and phases for two orthogonal polarizations simultaneously , , , , . Therefore, by designing proper birefringent metasurfaces and imposing different parabolic phase profiles, a multifunctional GSP-based metalens (PSFML) has been experimentally demonstrated (Figure 4B), which can split orthogonal linear light polarizations and focus them into different focal spots . The functionalities of polarization splitting and focusing are visualized in the right lower panel of Figure 4B, where it is readily seen that reflected light is focused in different positions for different polarizations. As can be seen, the PSFML deflects TM polarization to the right-hand side and TE polarization to the left-hand side from the origin, while two spots are symmetrically generated when two polarizations are involved. The proof-of-concept PSFML features high efficiency (up to ~65%) and polarization extinction ratio (up to ~30 dB), exhibiting broadband response in the near-infrared ranging from 750 to 950 nm.
Aside from the excellent focusing function, metalenses can provide scanning ability with controllable focus . Figure 4C schematically shows the operating principle of the frequency-controlled focus scanning reflective metasurfaces, where the focal position is well controlled by varying the operating frequency for the plane waves incident at an angle of 45°. Specifically, the proposed flat metalens is successfully implemented by carefully controlling the dispersions of the elements with different physical dimensions, which are consisting ring and I-shaped resonators on top of a ground dielectric substrate. The frequency-controlled focus scanning ability has been numerically and experimentally confirmed for several wavelengths, with the lower panel of Figure 4C displaying the power intensity distribution along the x-axis at z=40 mm and y=0 mm. It can be seen that the measured results agree well with simulation results, which clearly demonstrate the focal-scanning ability of the flat metalens.
As conventional lenses suffer from severe chromatic aberrations because of the dispersion of the phase accumulated during light propagation, metalenses also face the same problem because of the intrinsic spectral dispersion of the meta-atoms and the fact that phase gradient metasurfaces follow the diffraction theory . As such, metalenses imply a change in focal length proportional to λ−1, thus degrading the performance of an imaging system. Recently, several efforts have been dedicated to eliminate chromatic aberration of metalens by either dispersive phase compensation using multiple meta-atoms ,  or birefringent unit cells with equal focal length for orthogonally polarized light at two different wavelengths . However, the aforementioned metalenses can only compensate spectral dispersion at several discrete wavelengths, meaning that wavelengths of light deviating from the designed ones feature varying focal lengths. To achieve broadband achromatic metasurface devices, Tsai’s group employed the integrated resonant unit element in metasurfaces with smooth and linear phase dispersion combined with the PB phase . The proof-of-concept broadband achromatic metalens is capable of focusing light with arbitrary wavelength at the same focal plane with the efficiency on the order of ~12% over a continuous wavelength region from 1200 to 1680 nm for CP incidences in a reflection scheme (Figure 4D). Additionally, this approach is not limited to achromatic lens and thus can be easily used to realize various flat achromatic devices, such as achromatic gradient metasurfaces. Freshly, a polarization-insensitive achromatic metalens was demonstrated by Capasso’s group to suppress the chromatic effect over a 60-nm bandwidth in the visible for LP incidences, which consists of titanium dioxide (TiO2) nanopillar array on a silica-coated silver substrate. However, it is worth pointing out that the phase gradient is mainly due to the waveguide mode supported by the top TiO2 nanopillar. In this way, the structure cannot be regarded as GSPM .
Conventional holograms, based on the phase accumulation along the optical path, suffer from wavelength-scale thickness and large pixel size, which is usually tens of wavelengths for visible light, resulting in severe issues of high-order diffraction and twin image. In contrast, arbitrary amplitude, phase, and polarization can be locally engineered by introducing resonant meta-atoms and tailoring their lateral dimensions or orientations while the thickness is at the subwavelength scale. As such, metasurfaces are a good platform to realize all types of computer-generated holograms . Figure 5 shows several typical examples of GSPM-based holograms with excellent performance.
To reconstruct a holographic image, a 2π-phase modulation with quasi-uniform amplitude distribution is usually required, which can be achieved using resonant phase through deliberately optimizing the structure geometry for LP light at the design wavelength. Figure 5A schematically illustrates the first demonstration of GSPM-based hologram that generates reflective hologram at visible wavelengths . Particularly, by utilizing the polarization-sensitive response of the gold cross nanoantennas, two holograms with different patterns are simultaneously encoded on the same GSPM, and each image is reconstructed with an LP light with orthogonal polarizations, which works as a metasurface-based polarization-multiplexing technique. The performance of polarization controlled dual images was experimentally verified, where the projected patterns transitioned from NTU to RCAS once the input polarization is changed from x- to y-direction. The measured efficiency at λ=780 nm reaches 18% at 15° incidence and drops down to 8.5% when the angle is increased to 45°, which is greater than the existing metaholograms with fishnet  and complementary V-shaped structures . We note that the efficiency of GSPM-based hologram can be greatly increased to around 50% at a near-infrared wavelength of 1550 nm, resulting from the decreased metal loss .
Besides polarization-multiplexing, wavelength-multiplexed metahologram operating under LP light illumination was also demonstrated, where different holographic images were constructed by addressing the incident wavelength . With proper design of the aluminum nanorods with different sizes, resonances with narrow bandwidths are realized, thereby allowing for implementation of the multicolor metahologram with two-phase modulation at three primary colors (Figure 5B). However, it should be noted that the efficiency of wavelength-multiplexing metaholograms is inevitably low, less than 1%, ascribed to the trade-off between wavelength and phase control to minimize the crosstalk between different color images by spectrally separating the multiple resonances .
In the previous GSPM-based holograms, phase modulation is achieved by geometry-dependent resonance, thereby limiting the operating bandwidth because of the intrinsic dispersive resonance of meta-atoms. As an alternative approach to metaholograms, PB phase holograms gain the required phase on the cross-polarized light with identical geometric parameters but spatially varying orientations, thereby releasing the process of design and fabrication. However, it is worth noting that although PB phase modulation is dispersionless, the working bandwidth of PB phase metaholograms is still limited, which is in line with Section 2.3. Because of the linear relationship between the PB phase and the orientated angle of the individual meta-atom, the design of multilevel-phase holograms becomes much simpler. As an example of GSPM holograms based on PB phase, Figure 5C schematically illustrates the holographic image reconstruction using a GSPM that introduced 16-level geometric phases , where the MIM structure was optimized as a broadband half-wave plate maintaining high polarization conversion. Remarkably, this metahologram has a diffraction efficiency of 80% at 825 nm resulting from the antenna-orientation controlled PB profile combined with the GSPM for high polarization conversion. Additionally, the dispersionless nature of the PB phase GSPM results in broadband operation between 630 and 1050 nm with a high window efficiency larger than 50%.
Similar to the resonant phase holograms, holography multiplexing, such as spin-multiplexing, can be easily accomplished with PB phase metasurfaces. Figure 5D illustrates a good example of helicity multiplexing with two sets of PB phase GSPMs, which functions as synthetic holograms for normal incident CP light in reflection . Two symmetrically distributed off-axis images are simultaneously reconstructed in one identical hologram by controlling the helicity of the incident light. Moreover, the positions of the symmetrically reconstructed holographic images are exchanged once the spin of the incident light is altered. The experimentally measured conversion efficiency is higher than 40% over a broad wavelength band ranging from 620 to 1020 nm, with the maximum conversion efficiency reaching 59.2% at λ=860 nm. Such high image quality and wide bandwidth are ascribed to the simple structure geometry and the geometric nature of the PB phase, which is more robust against fabrication tolerances and variation of material properties.
Besides the aforementioned functionalities, enhancing absorption of light in a controllable manner is also of central importance to numerous applications such as photovoltaics , photodetectors , , thermal emitters , and biomedicine . Whereas bulky plasmonic materials, for example, noble metals, already exhibit substantial absorption, coexisted strong reflection and poorly-defined spectra features severely limit their further applicability and practicality. To resolve this issue, a particular branch of metasurfaces, perfect absorbers with near-unity absorption and designer spectra, has emerged and garnered considerable interest in recent years. In this context, GSP resonators with intriguing characteristics provide novel possibilities to enhance light absorption processes. To achieve perfect absorption, a system needs to prohibit all transmission and reflection channels at a certain frequency range. In view of this, optically thick continuous films are generally adopted as bottom back reflectors to block all transmission, making GSPMs as single-port reflective systems. By rationally designing topmost metallic nanostructures and carefully choosing spacer thickness, reflection can be eliminated via the generation of strongly confined GSP resonances. This exotic effect has triggered intense research efforts to realize GSP absorbers with various functions, aiming at many different uses in a diverse set of scenarios , , . Here we divide GSP absorbers into two main categories based on their different operating bandwidth, i.e. narrowband absorbers and broadband absorbers.
Figure 6A illustrates a representative narrowband perfect absorber using GSP resonators as basic building blocks . The top-layer gold disks were laid out with subwavelength periodicity, in a similar manner to other GSPMs. However, unlike reflective GSPMs, GSP absorbers usually require a much thinner dielectric spacer to ensure strong interactions between the top nanoparticles and the bottom reflector, forming a closely coupled antisymmetric current distribution , . Such a field distribution resembles an electric dipole near a mirror, giving rise to a mirrored dipole approximately in antiphase and thereby canceling out far-field radiation of the system, which can be also regarded as MD resonance. A maximum storage of near-field energy is thus generated, and 99% absorption can be achieved in the near-infrared range at normal incidence, as shown in the right panel of Figure 6A. Because of the rotational symmetry of the design and its ultrathin thickness, such a high absorption is polarization independent and robust against a wide range of incident angles. Moreover, when the GSP absorber is covered with different dielectric surroundings, the spectra show evident shifts or noticeable changes in the reflectivity/absorptivity at a fixed frequency. These interesting features further provide the feasibility of GSP absorbers working as plasmonic sensors.
After absorbing light, the confined EM field inside GSP resonators then dissipates into thermal energy, resulting in localized heating response and associated photothermal effects. With nearly full absorption, GSP absorbers manifest a strikingly enhanced photothermal process that is not reachable with natural materials or conventional plasmonic nanoparticles. Even with a moderate incident power (i.e. 2.3 mW), GSP absorbers can generate an ultrahigh temperature elevation up to 900 K . As a consequence, the top metallic nanoparticles would undergo melting and reshaping by light irradiation, as displayed in Figure 6B. Further transient analysis demonstrated that such a remarkable temperature increase occurs in just a few nanoseconds . The strong photothermal effects and associated fast response make GSP absorbers a promising candidate for many temperature-related applications including biotechnology , modulation , , , laser-induced particle transfer , and the rapidly growing field of plasmonic color laser printing .
So far, the foregoing absorbers featuring single GSP resonance all possess narrowband response, which is desirable for certain practical scenes such as plasmonic sensing , . However, for many other applications, especially photovoltaics and energy harvesting, broadband absorption is usually highly in demand. The most intuitive way to realize broadband GSP absorbers is to utilize GSP resonators with multiresonant properties. For example, by combining distinct yet overlapped GSP resonances within one design, a broadband response can be obtained, as shown in Figure 6C. When four differently sized GSP resonators and nine identical nanoparticles constitute a single unit cell, average absorption of ~90% for unpolarized light can be achieved in the entire visible spectrum (400–750 nm) . This simple, yet powerful design principle can also be leveraged to other configurations and other material platforms . Indeed, metallic nanostructures on the top of GSP resonators can be further tailored to introduce additional functionalities such as polarization dependence and dual- or multi-band absorption , , , , , . We also note that chemically synthesized nanoparticles can serve as the top metallic nanostructures , , . By doing this, GSP absorbers can overcome the spatial limitation of top-down lithography techniques, realizing perfect absorption over truly macroscopic areas. In addition, by controlling the size distribution and stacking layers of self-assembled nanoparticles, spectral response of GSP absorbers can be further engineered toward different practical applications.
Figure 6D illustrate another route to attain broadband absorption, i.e. the exploitation of highly lossy metals instead of commonly used coinage metals . Periodic titanium nanodisks were placed at the top of a SiO2 spacer layer and a continuous gold film. The pronounced inherent loss in titanium decreases the quality factor of GSP resonances and consequently broadens the absorption bandwidth. The accompanied propagating SPP along the SiO2-gold interface provides an extra channel for light dissipation. Absorption exceeding 90% thus was experimentally observed over a wide spectrum ranging from 900 to 1825 nm. Other highly lossy materials such as tungsten and nickel also hold the promise for broadband absorption , , .
The printing techniques in industry, such as inkjet and laserjet methods, use micro-sized ink spots and result in limited dots per inch (DPI), which is usually less than 10,000. Worst still, the positioning of multiple color pigments is inaccurate. To increase the resolution of color printing, significant attention has been drawn to designing proper micro/nanostructures for structural colors, inspired by the natural colors. In this sense, metasurfaces are a good platform to generate vivid colors through a strong interaction between light and subwavelength nanostructures , , , . In particular, GSPMs, which enable deep-subwavelength light confinement , can print subwavelength color elements beating the diffraction limit as MIM structures can be tailored for efficient light absorption by properly changing the optical constant or the thickness of the middle spacer layer that separates the two metal films, as discussed in the previous section. Additionally, the operating wavelength of MIM structures can be easily tuned with the different dimensions, leading to flexible spatial and spectral localization of light. In this section, we will discuss recent advances in GSPM-based high-resolution color printing, and several kinds of representative GSPMs that have excellent quality or special functionalities of structured colors are presented.
Figure 7A shows an example of polarization-insensitive GSP color printing with MIM structure containing gold nanodisk arrays hovering on a continuous gold film with a thin dielectric layer in between . When the size of the top nanodisks and the thickness of the dielectric layer are changed, the reflected spectra shift accordingly, resulting in tunable color. With properly selected disk sizes, an optical image of the University of Southern Denmark logo was produced, revealing a high-quality color print in terms of high contrast, high color uniformity, and high reproduction fidelity. Remarkably, this GSPM color printing is insensitive to the incident angle, superior to other structured color that is angular-dependent. Additionally, a cover of transparent dielectric overlay on top the fabricated GSPM would not affect the color too much, as shown in the right panel of Figure 7A. One can clearly see that the overall color image remains except for a slight contrast variation, ascribed to the nature of GSP resonance that entirely concentrates the light within the MIM structure itself. Therefore, this GSPM color device can be protected with a transparent dielectric overlay without destroying the colors. Similarly, aluminum-based GSPMs have also been considered for full-color subwavelength printing .
To tackle the challenge of highly saturated and bright color generation with plasmonic structures and maintain the pixel design functionality, tandem nanodisk arrays were demonstrated to achieve full colors  (Figure 7B). Here the MIM sandwich nanodisks are employed, which support enhanced in-phase electric dipole mode showing blue shift with respect to a single metal disk. By combining this dipole mode with the Wood’s anomaly in periodic structures, a high narrow peak in reflection and a deep valley in transmission is produced, thus providing vivid red-green-blue colors in reflection and cyan-magenta-yellow colors in transmission. The right panel of Figure 7B displays the experimental reflection and transmission colors for fabricated capital letters in serif font Georgia captured by optical microscope, respectively, where the colors are kept complementary in reflection and transmission and the full details of all letters are clearly reproduced.
The previous GSPM-based color patterns are relying on electron-beam lithography, which is expensive and not suitable for post-processing customization. Worst still, the color patterns need to be redesigned and thus are not scalable when the electron-beam lithography is conducted. More recently, printing new colors on a prefabricated GSPM color substrate using laser post-writing was successfully demonstrated (Figure 7C) . Laser pulses induce transient local heat generation that results in melting and reshaping of the nanoimprinted gap-plasmon nanostructures , . Depending on the laser pulse energy density, different surface morphologies can be created, corresponding to different gap-plasmon resonances and consequently different color modulation. Based on this technique, all primary colors can be printed with a speed of 1 ns per pixel, resolution up to 127,000 DPI, and power consumption down to 0.3 nJ per pixel. As a final comment, it should be mentioned that although laser-induced morphology changes are irreversible, the structures can be repeatedly rewritten to exhibit new color resonances associated with morphologies of lower surface energy.
Polarization is an intrinsic property of EM waves. The capability of controlling the polarization states is of vital importance in a variety of practical applications. Typically, the conversion of polarization states is achieved by utilizing birefringent crystals, where the phase retardation between two orthogonally polarized wave components is gradually accumulated during light propagation. As a result, the polarization control components are usually bulky and suffer from specific thickness limitations, going against the general trend of integration and miniaturization in photonics, which in turn boosts the development of metasurface-based waveplates operating at different wavelengths. In this section, we review some examples of the GSPM waveplates with excellent performance (Figure 8).
Generally, GSPMs featuring anisotropic optical responses for two orthogonally polarized wave components, i.e. Ex and Ey, can be applied to the realization of ultrathin metasurfaces based quarter-wave plate , , ,  or half-wave plate (HWP) , , , , , , ,  when the reflection amplitudes Rxx and Ryy are equal, and the phase difference ΔΦ=Φxx−Φyy is π/2 or π, respectively. For instance, Anders et al. demonstrated that the MIM configurations in which the top metal layer consists of a periodic arrangement of nanobricks whose main axes function as orthogonal detuned scatterers, thus supporting detuned GSP resonances, can be designed to function as reflective broadband half-wave plates  (Figure 8A). By properly tuning the dimensions of nanobricks, two resonant dips can be found under the illumination of horizontally and vertically LP light, corresponding to the GSP resonances along the long- and short-axes, respectively. Therefore, the reflection efficiencies between the x- and y-polarized light exhibited nearly equal reflection amplitude ratio (Rxx/Ryy) accompanied by a phase difference close to 180° in the off-resonance regime between two GSP resonances, leading to broadband linear-to-linear polarization conversion, as shown in the left lower panel of Figure 8A. From the measured reflectivity at different analyzer angles (right lower panel of Figure 8A), it is clearly demonstrated that the reflected light is rotated by 90° with the maximum reflection reaching ~50% at the analyzer angle 130° in the wavelength range of 740–880 nm. The broadband polarization conversion can also be interpreted in the general mechanism of dispersion engineering by integrating a metallic metasurface possessing a strong yet dispersive interaction with light, with a dielectric interlayer, which has accumulative yet broadband interaction with light . In this way, the intrinsic dispersion of the metallic structures is perfectly canceled out by the thickness-dependent dispersion of the dielectric spacing layer, resulting in dispersion-free metasurface and thus wideband linear-to-linear or linear-to-circular polarization conversion.
Because of increased Ohmic loss, the efficiencies are slightly reduced when GSPMs are extended to near-infrared and visible regimes. To increase the efficiency, a hybrid GSPM was demonstrated to create broadband HWP . The metasurface in Figure 8B consists of an array of 45°-rotated silicon bricks, 200-nm PMMA spacer, and a sliver back-reflector, thereby considerably reducing the level of absorption within the metasurface in near-infrared range. The experimental measurements demonstrate that the cross-polarized reflection |rcr|2 remains above 97% while the co-polarized reflection |rco|2 is greatly suppressed from 1420 to 1620 nm, in good agreement with the simulations. Additionally, the measured polarization conversion ratio (PCR), defined as |rcr|2/(|rco|2+|rcr|2), remains above 98% over this bandwidth across a 200-nm bandwidth. Furthermore, by varying the meta-atom dimensions, the 2π phase coverage of the reflected cross-polarized light has been achieved while maintaining high conversion efficiency, opening the door to ultracompact reflective phase plates for vortex beam generation.
Though the GSPM-based wave plates can achieve excellent polarization conversion over a wide wavelength spectrum, the co-polarized light cannot be completely eliminated. As such, there will always be co-polarized reflected light remaining in space, merging together with the converted cross-polarized light and affecting the performance. To separate the co-polarized and cross-polarized reflection spatially, a background-free GSPM-based HWP was demonstrated by integrating several HWPs into a supercell to create a linear phase gradient for the cross-polarized reflected wavefront, resulting in anomalous reflection and routing the cross-polarized reflected light away (Figure 8C) . Specifically, four antennas are introduced by numerical simulations to form nearly 2π phase shift coverage for the cross-polarized reflected light with a phase increment of about π/2 between adjacent elements, and each antenna can serve as a HWP with high PCR. As shown in the right panel of Figure 8C, the cross-polarized reflection power (red curve) coincides with the total reflection power (blue curve), and the co-polarized reflection power (black curve) is approaching zero, confirming the required performance. In addition, the ratio between the power of the desired and orthogonal polarizations is above 20, and the integrated PCR exceeds 95%.
In contrast to the aforementioned ultrathin wave plates based on resonant phases, the PB phase can be used to realize the broadband wave plate , , . Here we would like to highlight a GSPM polarization generator, which is capable of producing reflected light beams of any polarizations from an LP light source (left panel of Figure 8D) . The metasurface chip is composed of six types of supercells, in which each supercell is made up of two subunits with different offsets. Each subunit features a constant phase gradient by rotating the orientations of aluminum nanoantennas gradually with a constant step size of θ from 0 to π. Upon the incidence of LP light, two CP reflection beams with opposite handedness are redirected into the ±1 diffraction orders. If two supercells with different offsets and opposite nanoantenna orientations are introduced, different LP states could be produced. From the right panel of Figure 8D, it is clearly seen that six light beams with different polarizations are achieved in different reflection angles simultaneously, including four different linear polarizations and two circular polarizations.
The rapid development of metasurface-based polarization conversion has motivated the implementation of polarimeters, which enables direct detection of polarization states. As the polarization state of light is uncorrected with its intensity and wavelength, it is intrinsically difficult to be measured. Typically, polarization states are probed by using a set of discrete polarizers and wave plates consecutively placed in the light path in front of a detector. By measuring the transmitted intensity through these components for polarization control, the Stokes parameters that uniquely define the polarization state can thus be determined. In this way, polarimeters based on conventional optical components are bulky and complicated, not compatible with the integrated photonic system. In recent years, polarization detection using GSPMs has gained much attention because of their ease of fabrication, compact footprint, and design flexibility , , . For example, photonic spin Hall effect has been utilized to implement chiroptical spectroscopy based on PB phase GSPM, which can determine certain parts of the polarization state, i.e. the degree of circular polarization . Such chiroptical spectroscope can spatially separate the LCP and RCP light, and the scattering intensity ratio of these two CP states can be used to retrieve the ellipticity and handedness of the incident light. However, this type of chiroptical spectroscopy has only two channels and cannot fully characterize the polarization state.
Recently, on-chip GSPM polarimeter has been proposed to simultaneously determine all polarization states by a metagrating . As shown in Figure 9A, the metagrating consists of three interweaved GSPMs in which each functions as a polarization beam-splitter for a certain polarization basis. Once an arbitrary polarized incident beam impinges on the metasurface, the beam will be diffracted spatially and routed in six predesigned directions. By conducting parallel measurements of the six diffraction intensities, the Stokes parameters can be quickly retrieved, there allowing one to analyze an arbitrary polarization state. At the design wavelength of 800 nm, the experimental diffraction contrast obtained by averaging three successive measurements agrees well with the input normalized Stokes parameters, repeating well the Poincaré sphere with the two-norm deviation between Stokes parameters and measured diffraction contrast being around ~0.1 (right panel of Figure 9A). Additionally, the designed GSPM could operate well in the wavelength range of 750–850 nm. Following this concept, a waveguide-based polarimeter was proposed based on GSPM, where the polarization state can be directly determined by calculating the coupling efficiencies of launched waveguide modes in six different directions .
To further extend the functionality of polarimeter, segmented and interleaved GSPMs have been proposed and demonstrated to realize spectropolarimeters, which enable simultaneous characterization of the polarization state and spectrum of incident light , , . Starting with segmented GSPM spectropolarimeter, Figure 9B schematically illustrates the principle of a spectropolarimeter metadevice composed of six GSPMs arranged in a 2×3 array corresponding to horizontal (0°), vertical (90°), ±45°, RCP and LCP analyzers, which steers different polarization and spectrum components into distinct positions . Once a probe beam with a certain polarization state is incident on the metadevice, six intensity peaks are generated in the far-field. After careful calibration, the polarization detection at a certain wavelength is determined. Additionally, the measured angular dispersion for the LCP and RCP channels are 0.053°/nm and 0.024°/nm, respectively (right panel of Figure 9B), proving the potential of spectral measurement. Later on, an interleaved GSPM spectropolarimeter was demonstrated by randomly interspersing three linear phase profiles associated with different meta-atom subarrays (Figure 9C) . As shown in the middle panel of Figure 9C, the measured Stokes parameters are in good correspondence with the calculated ones on the Poincaré sphere for an arbitrary incident beam at the design wavelength of 760 nm, revealing the excellent capability of polarization probe. Furthermore, the interleaved spectropolarimeter has good spectral resolving power (λ/Δλ), which was measured to be 13 when the diameter is only 50 μm.
Although these two spectropolarimeters show good performance, they are very sensitive to the size of incident beam, and the calibration is always needed. To increase the robustness of polarization and spectrum detection, a self-calibrating segmented GSPM spectropolarimeter, featuring the beam-size-invariant nature, has been demonstrated . Figure 9D shows the GSPM spectropolarimeter with center-symmetrical configuration, which consists of three GSPMs that occupy 120° circular sectors each. Upon the excitation of a normally incident beam with arbitrary polarization state, six diffraction spots are generated in the predesigned directions, whose polar angles are proportional to the wavelength while the contrasts in the corresponding diffraction intensities provide a direct measure of the incident polarization state. Here it should be pointed out that uneven illumination of the three GSPMs caused by any misalignment will not affect the final result and no calibration is needed at all as the retrieved Stokes parameters are only associated with the relative diffraction contrasts for three polarization bases. The proof-of-concept 96-μm-diameter spectropolarimeter operating in the wavelength range of 750–950 nm exhibits excellent polarization sensitivity with beam-size-invariant property. Additionally, the experimentally measured spectral resolving power is estimated to be ~15.2 (right panel of Figure 9D).
In addition to manipulating the propagating waves (PWs) in free space, another exciting feature of GSPMs is the capability of efficiently launching surface waves (SWs), such as SPPs and spoof SPPs in low-frequency range, in which arbitrary in-plane wave vector can be designed to match with the wave vector of SWs, distinct from the prism or grating couplers based on resonant coupling between PW and SW . In this way, a nearly perfect PW-SW unidirectional conversion for any incidence angle larger than a critical value can be achieved , . However, the generated SW is not the eigenmode of the GSPM and more likely to scatter out because of the spatial inhomogeneity of the GSPM. Additionally, the conversion efficiency decreases significantly with the increased size of incident beam, arising from the significant scattering caused by inter-supercell discontinuities . To solve these issues, a new SPP meta-coupler has been demonstrated, which consists of a transparent gradient metasurface placed at a certain distance above the target plasmonic metal . As shown in Figure 10A, the incident wave is first converted into a driven SW bound on the metasurface and then resonantly coupled to the eigenmode (i.e. SPP wave in optical frequencies) on the plasmonic metal. Based on this new configuration, a theoretical efficiency up to 94% has been predicted by model calculations, eliminating the nonnegligible issues that severely affect the coupling efficiency . As a practical realization, a realistic device operating in the microwave regime has been fabricated, which exhibits a spoof-SPP conversion efficiency of ~73% in both near-field and far-field experiments (right panel of Figure 10A). The measured conversion efficiency is much higher than those of all other available devices operating in this frequency domain.
Though the aforementioned GSPM couplers show good performance, the direction of launched SWs is predefined without any tunability. Very recently, unidirectional polarization-controlled SPP coupling has been demonstrated using GSPMs, so that the direction of SPP excitation was switched by the polarization of an LP light  or the handedness of a CP incident beam , . Here we would like to highlight an efficient unidirectional polarization-controlled SPP coupler by using GSPMs operating at telecom wavelengths . Specifically, arrays of GSP resonators that would produce two independent orthogonal reflection phase gradients in two respective linear polarizations have been proposed, where basic GSP unit cell consists of a gold nanobrick atop a 50-nm-thick dielectric spacer layers covering 80-nm-thick gold films supported by glass substrates (Figure 10B). Therefore, the incident radiation (with arbitrary polarization) can efficiently (up to 40%) be converted into SPPs propagating in orthogonal directions dictated by the phase gradients. The right panel of Figure 10B displays the recorded images from leakage radiation microscopy for a free-space wavelength of 1500 nm, featuring the coupling efficiency C of ~25% for either of two linear polarizations and the directivity D of SPP excitation exceeding 100.
To further extend the functionalities of GSPM-based SW coupler, a bifunctional GSPM for visible light has been proposed, which is capable of efficiently converting x-polarized radiation into SPPs and anomalously steering the reflected light for y-polarization simultaneously (Figure 10C) . Here two different linear phase gradients, ζx and ζy, are introduced along the x-direction in respective linear polarizations, distinct from the previous polarization-controlled GSPM coupler processing orthogonal phase gradients . Specifically, the reflection phase gradient for the x-polarization is equal to the wave-vector of excited SPPs bounded on the air-dielectric-metal interface (i.e. ζx=kSPP), while the phase-gradient is smaller than the wave-vector of propagating light in free space when the incident light is y-polarized (i.e. ζy<k0). Using an x-polarized Gaussian beam optimally positioned to maximize the coupling efficiency of the right-propagating SPPs, unidirectional excitation of SPPs over a wide spectrum range was observed. Middle panel of Figure 10C shows the measured coupling efficiencies over the wavelength range from 600 to 650 nm, where an averaged coupling efficiency Cr of above 25% has been achieved, and the measured extinction ratio is larger than 20 dB. Apart from the unidirectional SPP excitation under x-polarization, the designed metasurface can function as a broadband beam steerer for the y-polarized light. From the right panel of Figure 10C, one can clearly see that most of the reflected light is contained within the +1 diffraction order in a broad wavelength regime, verifying the broadband steering for the y-polarized light, albeit with some discrepancies compared with the simulation results regarding the efficiencies, especially at short-wavelength range.
The aforementioned unprecedented phenomena and applications of GSPMs are mainly relying on their passive properties. An as-yet largely unexplored milestone in the field of metasurfaces is to achieve a dynamically reconfigurable metasurfaces, which can extend their exotic passive properties and enables dynamic wavefront control over an ultrathin surface, such as active spatial light modulation, dynamic beam steering, and tunable ultrathin lenses , .
Reconfigurable GSPMs can be realized by altering the shape of individual meta-atoms, or by manipulating the near-field interactions between them, which can be accomplished through microelectromechanical systems , . For example, a voltage-tunable Fano metamaterial absorber, which simultaneously supports thermomechanically coupled optical and mechanical resonances, has been utilized to control mechanical damping with light (Figure 11A) . Based on the plasmomechanical parametric gain mechanism, optically pumped coherent mechanical oscillations with a very low threshold have been experimentally demonstrated over a ~4-THz bandwidth. In addition, optical damping of mechanical resonance was also achieved through the reverse effect.
Besides the structure aspect, an alternative strategy to realize reconfigurable GSPMs is to incorporate active functional materials into metasurfaces that provide tunable optical constants through voltage bias, thermal heating, optical pump, and magnetic field. In particular, owing to the excellent mechanical and optoelectronic properties and the largely tunable carrier density, graphene processes largely tunable optical conductivity in the mid-infrared and THz frequency ranges once the doping is adjusted through changing the bias voltage, resulting in the ultimate reconfigurability and multiple functionalities with graphene GSPMs , , , . Figure 11B schematically shows the wide-range phase modulator, which is consisting of magnetic resonators and gate-controlled graphene . A variation in the graphene resistance, associated with the gate, drastically changes the resonance behavior of the proposed metasurface, namely, driving an underdamped resonator to an overdamped resonator, leading to a phase modulation of ±180°. It should be noted that the maximum phase modulation achieved using the present scheme is not limited to 180°. A much larger maximum phase modulation can be achieved using two slightly different graphene metasurfaces gated independently (right panel of Figure 11B). The large phase modulation is ascribed to the fact that gate doping modulates the phase in the range from approximately −180° to 0° at frequencies immediately above the resonance frequency in device A, while simultaneously suppressing the phase in the range from approximately 180° to 0° at frequencies slightly below the resonance in device B. Consequently, an extremely large phase modulation range of 243° at 0.48 THz has been achieved. The combination of graphene and GSPMs generates a versatile new platform for applications in the mid-infrared and THz frequency regime. Here it should be emphasized that highly doped graphene sheet or graphene nanostructures are needed if one wants to realize tunable graphene metasurfaces in the near-infrared and even the visible range .
As alternative functional materials with excellent reconfigurability, phase-change materials has been extensively utilized in optical data-storage systems and photonic devices because of their outstanding switchable optical properties by thermal, laser, or electrical current pulses with controlled duration and intensity , , , , , , , , . On the transformation from the amorphous to crystalline state, their optical properties change drastically. For instance, germanium-antimony-tellurium (GST) chalcogenide glass has recently been used to demonstrate all-optical, non-volatile GSPM switch . However, GST needs very high temperature to return back to amorphous state, which may limit the reversibility of GST GSPMs. As such, vanadium dioxide (VO2), a classical transition metal oxide that behaves as an insulator and metal at temperature below and above the insulator-to-metal transition temperature (~340 K in bulk crystals), respectively, is a good choice to facilitate reconfigurable GSPMs , , . Very recently, Werner’s group demonstrated a hybrid metasurface platform that achieves electrically triggered multifunctional control in the mid-infrared region by integrating a VO2 film of nanoscale thickness into an optical metasurface absorber, where the active VO2 layer is enabled as a function of the electrical current flowing through the nanoengineered metal layer supporting the optical resonances and the electrical control simultaneously (Figure 11C) . When the electrical current (I) increases to 1.20 A, absolute tunings of 80% and 75%, which correspond to 75-fold and 5-fold relative modulations relative to the static reflectance when I=0 A, are observed at wavelengths of 3.05 and 3.85 μm, respectively. In addition to the drastic tuning, a continuous spectrum tuning is observed for in-between current values before saturation. The proposed electrically active GSPM was further used to tune the spatial dependence of infrared signals, verifying a more sophisticated spatial modulation scheme in the infrared range.
In addition to graphene and phase-change materials, a variety of other functional materials has been successfully used to realize active GSPMs, for instance, transparent conducting oxides , , . Among the transparent conducting oxide materials, indium tin oxide (ITO) has been widely used as significant voltage-tunable changes in the refractive index can be achieved through accumulation of carriers through the electric field effect when applying an electrical bias between metal and ITO. In particular, the electric field can change the real part of the permittivity of ITO in the accumulation layer from positive to negative, lending to large electric field enhancement in the epsilon-near-zero region. Here we would like to highlight a gate-tunable GSPM that enables dynamic electrical control of the phase and amplitude of the reflected plane wave . As shown in Figure 11D, the metasurface is composed of a continuous gold film and a thin ITO film followed by a thin alumina film on which the connected gold stripe nanoantenna array is patterned. Because of the field-effect modulation of the complex refractive index of ITO layers incorporated into metasurface antenna elements, a phase shift of 180° and ~30% change in the reflectance was experimentally measured by applying a 2.5-V gate bias. Furthermore, by electrically controlling subgroups of metasurface elements properly, electrical switching of ± first order diffracted beams has been achieved. When the applied voltage is increased above ~1.5 V, the ±1 order diffracted beams appear and become dominant while the zero-order diffracted beam intensity reduces, in good agreement with the simulation results (right panel of Figure 11D).
As a final comment, it should be noted that while here we are mainly focused on the reconfigurable GSPMs in high frequencies, similar reconfigurability can be achieved by varying the voltages applied on the varactor diodes incorporated in the GSP meta-atoms at low frequencies , , .
As a rapidly developing research area, GSPMs are expanding very fast every day, and novel functionalities continue to be explored benefitting from their unprecedented control of EM waves and recent advances in nanofabrication technologies. Recently, compact analog computing based on metasurfaces has emerged as a promising candidate for real-time and parallel continuous data processing , , , . Motivated by the theoretical work , the first metasurface that allows for mathematical operations on visible light was demonstrated using GSPM working in a reflection setup (left panel of Figure 12A) . Proof-of-concept experiments have been conducted to realize a high-quality poor-man’s integrator operating at the wavelength of 800 nm, as shown in the right panel of Figure 12A.
The cloaking device that can cause objects to be partially or wholly invisible to the incoming waves has attracted a great deal of interest, becoming one of the most important applications of GSPMs. It was numerically proved that a carpet cloak can be achieved to mimic the reflection from a planar ground by compensating the local phase with gradient metasurfaces, thereby hiding arbitrary objects in a predesigned region , . Lately, an ultrathin invisibility skin cloak for visible light was experimentally demonstrated by MIM resonators with spatially distributed phase shifts (Figure 12B) . With the cloak off, a strong contrast between the cloaked region and the surrounding reflective surface revealed the object, but with the cloak on, the contrast became indistinguishable, rendering the object completely invisible. In contrast to bulky cloaks realized via 3D metamaterials, the GSPM cloak was only 80 nm thick. Because of the scalable feature of the metasurface, the GSPM-based cloak can be extended to microwave for full polarization , .
Metasurface-based vortex-beam generators have recently gained tremendous interest because of potential applications in optical trapping, high-resolution microscopy, and quantum communication technology , . A variety of GSPMs has been demonstrated to generate vortex-beams with arbitrary orbital angular momentum (OAM) values , , , . Very recently, Chen’s group proposed and experimentally demonstrated an efficient approach to generate multiple OAM states with different topological charges in four separate channels and arbitrarily control their superpositions using PB phase GSPM (Figure 12C) . By changing the helicity of incident light from RCP to LCP, the positions of all OAM beams are changed to their centrosymmetric positions with respect to the axis of incident light. Meanwhile, the signs of all topological charges are changed from plus to minus.
In addition to unprecedented applications relying on linear light-matter interaction, metasurfaces offer great capability of designing nonlinear photonic devices with remarkable performance, which is ascribed to the unique ability to localize EM fields at nanoscale volumes, together with significantly relaxed phase matching condition because of subwavelength thickness , . In this sense, the combination of strong field enhancement within the gap area and intrinsic nonlinearities of metals make GSPM an excellent platform to realize effective nonlinear optical processes , , , , . To further enhance the nonlinear response, metal/multiple-quantum-well (MQW) hybrid metasurfaces were experimentally demonstrated by coupling EM mode in GSPM with intersubband transitions of MQW semiconductor heterostructure , , . The GSPMs not only enhance both the fundamental and second harmonic frequencies but also tailor the near-field polarizations, converting giant nonlinear susceptibility of MQW heterostructures and thus resulting in giant nonlinearity (Figure 12D) .
The research filed for GSPMs has been rapidly growing as a platform for novel EM properties and functionalities for the reflected light, because of its simplified realization that only requires a single step of lithography. In this paper, we have briefly reviewed the fundamental physics and the latest developments of representative applications in the research area of GSPMs.
Owing to the unprecedented control over reflected light with surface-confined planar components, GSPMs are expected to have an important impact in other promising areas converting both fundamental research and real applications, beyond what we have discussed in this review article. For example, GSPMs provide an ideal platform to create multifunctional metadevices, which efficiently integrate multiple functionalities into a single ultrathin device with compact footprint , . Additionally, the so-called quantum metasurfaces, which combine existing quantum systems with metasurfaces, is an emerging area of research , , .
This work was funded by the European Research Council (the PLAQNAP project, grant 341054) and the University of Southern Denmark (SDU2020 funding).
 Born M, Wolf E. Principles of optics: electromagnetic theory of propagation, interference and diffraction of light. New York, NY, Elsevier, 2013. Search in Google Scholar
 Yu N, Genevet P, Kats MA, et al. Light propagation with phase discontinuities: generalized laws of reflection and refraction. Science 2011;334:333–7.10.1126/science.121071321885733 Search in Google Scholar
 Cai W, Shalaev V. Optical metamaterials: fundamentals and applications. New York, NY, Springer, 2010. Search in Google Scholar
 Luo X, Pu M, Ma X, Li X. Taming the electromagnetic boundaries via metasurfaces: from theory and fabrication to functional devices. Int J Antennas Propag 2015;2015:204127. Search in Google Scholar
 Minovich AE, Miroshnichenko AE, Bykov AY, Murzina TV, Neshev DN, Kivshar YS. Functional and nonlinear optical metasurfaces. Laser Photon Rev 2015;9:195–213.10.1002/lpor.201400402 Search in Google Scholar
 Yu N, Genevet P, Aieta F, et al. Flat optics: controlling wavefronts with optical antenna metasurfaces. IEEE J Sel Top Quant Electron 2013;19:4700423.10.1109/JSTQE.2013.2241399 Search in Google Scholar
 Hum SV, Perruisseau-Carrier J. Reconfigurable reflectarrays and array lenses for dynamic antenna beam control: a review. IEEE Trans Antennas Propag 2014;62:183–98.10.1109/TAP.2013.2287296 Search in Google Scholar
 Rui G, Zhan Q. Tailoring optical complex fields with nano-metallic surfaces. Nanophotonics 2015;4:2–25. Search in Google Scholar
 Bozhevolnyi SI, Søndergaard T. General properties of slow-plasmon resonant nanostructures: nano-antennas and resonators. Opt Express 2007;15:10869–77.10.1364/OE.15.01086919547444 Search in Google Scholar
 Moreau A, Cirací C, Mock JJ, et al. Controlled-reflectance surfaces with film-coupled colloidal nanoantennas. Nature 2012;493:86–9. Search in Google Scholar
 Ma S, Xiao S, Zhou L. Resonant modes in metal/insulator/metal metamaterials: an analytical study on near-field couplings. Phys Rev B 2016;93:045305.10.1103/PhysRevB.93.045305 Search in Google Scholar
 Jung J, Søndergaard T, Bozhevolnyi SI. Gap plasmon-polariton nanoresonators: scattering enhancement and launching of surface plasmon polaritons. Phys Rev B 2009;79:035401.10.1103/PhysRevB.79.035401 Search in Google Scholar
 Qu C, Ma S, Hao J, et al. Tailor the functionalities of metasurfaces based on a complete phase diagram. Phys Rev Lett 2015;115:235503.10.1103/PhysRevLett.115.23550326684124 Search in Google Scholar
 Miao Z, Wu Q, Li X, et al. Widely tunable terahertz phase modulation with gate-controlled graphene metasurfaces. Phys Rev X 2015;5:041027. Search in Google Scholar
 Pors A, Nielsen MG, Valle GD, Willatzen M, Albrektsen O, Bozhevolnyi SI. Plasmonic metamaterial wave retarders in reflection by orthogonally oriented detuned electrical dipoles. Opt Lett 2011;36:1626–8.10.1364/OL.36.00162621540949 Search in Google Scholar
 Monticone F, Estakhri NM, Alù A. Full control of nanoscale optical transmission with a composite metascreen. Phys Rev Lett 2013;110:203903.2516741110.1103/PhysRevLett.110.203903 Search in Google Scholar
 Ding X, Monticone F, Zhang K, et al. Ultrathin Pancharatnam-Berry metasurface with maximal cross-polarization efficiency. Adv Mater 2015;27:1195–200.10.1002/adma.20140504725545285 Search in Google Scholar
 Sun S, He Q, Xiao S, Xu Q, Li X, Zhou L. Gradient-index meta-surfaces as a bridge linking propagating waves and surface waves. Nat Mater 2012;11:426–31.10.1038/nmat329222466746 Search in Google Scholar
 Asadchy VS, Albooyeh M, Tcvetkova SN, Daz-Rubio A, Ra’di Y, Tretyakov SA. Perfect control of reflection and refraction using spatially dispersive metasurfaces. Phys Rev B 2016;94:075142.10.1103/PhysRevB.94.075142 Search in Google Scholar
 Ra’di Y, Sounas DL, Alù A. Metagratings: beyond the limits of graded metasurfaces for wave front control. Phys Rev Lett 2017;119:067404.10.1103/PhysRevLett.119.06740428949646 Search in Google Scholar
 Wong AMH, Eleftheriades GV. Perfect anomalous reflection with a bipartite Huygens’ metasurface. Phys Rev X 2018;8:011036. Search in Google Scholar
 Daz-Rubio A, Asadchy VS, Elsakka A, Tretyakov SA. From the generalized reflection law to the realization of perfect anomalous reflectors. Sci Adv 2017;3:e1602714.2881964210.1126/sciadv.1602714 Search in Google Scholar
 Pancharatnam S. Generalized theory of interference, and its applications. Proc Indian Acad Sci Sect 1956;44:247–62. Search in Google Scholar
 Bomzon Z, Kleiner V, Hasman E. Pancharatnam–Berry phase in space-variant polarization-state manipulations with subwavelength gratings. Opt Lett 2001;26:1424–6.1804962610.1364/OL.26.001424 Search in Google Scholar
 Bomzon Z, Biener G, Kleiner V, Hasman E. Space-variant Pancharatnam–Berry phase optical elements with computer-generated subwavelength gratings. Opt Lett 2002;27:1141–3.1802638710.1364/OL.27.001141 Search in Google Scholar
 Biener G, Niv A, Kleiner V, Hasman E. Formation of helical beams by use of Pancharatnam–Berry phase optical elements. Opt Lett 2002;27:1875–7.1803338710.1364/OL.27.001875 Search in Google Scholar
 Hasman E, Kleiner V, Biener G, Niv A. Polarization dependent focusing lens by use of quantized Pancharatnam-Berry phase diffractive optics. Appl Phys Lett 2003;82:328–30.10.1063/1.1539300 Search in Google Scholar
 Li Z, Palacios E, Butun S, Aydin K. Visible-frequency metasurfaces for broadband anomalous reflection and high-efficiency spectrum splitting. Nano Lett 2015;15:1615–21.2566481510.1021/nl5041572 Search in Google Scholar
 Jiang S-C, Xiong X, Hu Y-S, et al. High-efficiency generation of circularly polarized light via symmetry-induced anomalous reflection. Phys Rev B 2015;91:125421.10.1103/PhysRevB.91.125421 Search in Google Scholar
 Huang Y, Zhao Q, Kalyoncu SK, et al. Phase-gradient gap-plasmon metasurface based blazed grating for real time dispersive imaging. Appl Phys Lett 2014;104:161106.10.1063/1.4872170 Search in Google Scholar
 Zhang L, Hao J, Qiu M, et al. Anomalous behavior of nearly-entire visible band manipulated with degenerated image dipole array. Nanoscale 2014;6:12303–9.10.1039/C4NR03163F25163776 Search in Google Scholar
 Gao S, Yue W, Park C-S, Lee S-S, Kim E-S, Choi D-Y. Aluminum plasmonic metasurface enabling a wavelength-insensitive phase gradient for linearly polarized visible light. ACS Photonics 2017;4:322–8.10.1021/acsphotonics.6b00783 Search in Google Scholar
 Pors A, Nielsen MG, Eriksen RL, Bozhevolnyi SI. Broadband focusing flat mirrors based on plasmonic gradient metasurfaces. Nano Lett 2013;13:829–34.10.1021/nl304761m23343380 Search in Google Scholar
 Boroviks S, Deshpande RA, Mortensen NA, Bozhevolnyi SI. Multifunctional meta-mirror: polarization splitting and focusing. ACS Photonics 2017 (in press). doi: 10.1021/acsphotonics.7b01091. Search in Google Scholar
 Yi H, Qu S-W, Chen BJ, Bai X, Ng KB, Chan CH. Flat terahertz reflective focusing metasurface with scanning ability. Sci Rep 2017;7:3478.10.1038/s41598-017-03752-328615684 Search in Google Scholar
 Aieta F, Genevet P, Kats MA, et al. Aberration-free ultrathin flat lenses and axicons at telecom wavelengths based on plasmonic metasurfaces. Nano Lett 2012;12:4932–6.2289454210.1021/nl302516v Search in Google Scholar
 Khorasaninejad M, Aieta F, Kanhaiya P, et al. Achromatic metasurface lens at telecommunication wavelengths. Nano Lett 2015;15:5358–62.2616832910.1021/acs.nanolett.5b01727 Search in Google Scholar
 Khorasaninejad M, Shi Z, Zhu AY, et al. Achromatic metalens over 60 nm bandwidth in the visible and metalens with reverse chromatic dispersion. Nano Lett 2017;17:1819–24.10.1021/acs.nanolett.6b0513728125234 Search in Google Scholar
 Chen WT, Yang KY, Wang CM, et al. High-efficiency broadband meta-hologram with polarization-controlled dual images. Nano Lett 2013;14:225–30.24329425 Search in Google Scholar
 Yifat Y, Eitan M, Iluz Z, Hanein Y, Boag A, Scheuer J. Highly efficient and broadband wide-angle holography using patch-dipole nanoantenna reflectarrays. Nano Lett 2014;14:2485–90.2464605710.1021/nl5001696 Search in Google Scholar
 Schuller JA, Barnard ES, Cai W, Jun YC, White JS, Brongersma ML. Plasmonics for extreme light concentration and manipulation. Nat Mater 2010;9:193–204.10.1038/nmat263020168343 Search in Google Scholar
 Liu X, Tyler T, Starr T, Starr AF, Jokerst NM, Padilla WJ. Taming the blackbody with infrared metamaterials as selective thermal emitters. Phys Rev Lett 2011;107:045901.10.1103/PhysRevLett.107.04590121867022 Search in Google Scholar
 Ra’di Y, Simovski CR, Tretyakov SA. Thin perfect absorbers for electromagnetic waves: theory, design, and realizations. Phys Rev Appl 2015;3:037001.10.1103/PhysRevApplied.3.037001 Search in Google Scholar
 Nielsen MG, Pors A, Albrektsen O, Bozhevolnyi SI. Efficient absorption of visible radiation by gap plasmon resonators. Opt Express 2012;20:13311–9.2271435910.1364/OE.20.013311 Search in Google Scholar
 Ding F, Dai J, Chen Y, Zhu J, Jin Y, Bozhevolnyi SI. Broadband near-infrared metamaterial absorbers utilizing highly lossy metals. Sci Rep 2016;6:39445.2800071810.1038/srep39445 Search in Google Scholar
 Jain PK, Huang X, El-Sayed IH, El-Sayed MA. Noble metals on the nanoscale: optical and photothermal properties and some applications in imaging, sensing, biology, and medicine. Accounts Chem Res 2008;41:1578–86.10.1021/ar7002804 Search in Google Scholar
 Shi Y, Chen X, Lou F, et al. All-optical switching of silicon disk resonator based on photothermal effect in metal–insulator–metal absorber. Opt Lett 2014;39:4431–4.2507819510.1364/OL.39.004431 Search in Google Scholar
 Gong H, Chen X, Qu Y, Li Q, Yan M, Qiu M. Photothermal switching based on silicon Mach–Zehnder interferometer integrated with light absorber. IEEE Photonics J 2016;8:1–10. Search in Google Scholar
 Gong H, Yang Y, Chen X, et al. Gold nanoparticle transfer through photothermal effects in a metamaterial absorber by nanosecond laser. Sci Rep 2014;4:6080. Search in Google Scholar
 Meng L, Zhao D, Ruan Z, Li Q, Yang Y, Qiu M. Optimized grating as an ultra-narrow band absorber or plasmonic sensor. Opt Lett 2014;39:1137–40.2469069010.1364/OL.39.001137 Search in Google Scholar
 Aydin K, Ferry VE, Briggs RM, Atwater HA. Broadband polarization-independent resonant light absorption using ultrathin plasmonic super absorbers. Nat Commun 2011;2:517.10.1038/ncomms152822044996 Search in Google Scholar
 Zhang B, Zhao Y, Hao Q, et al. Polarization-independent dual-band infrared perfect absorber based on a metal-dielectric-metal elliptical nanodisk array. Opt Express 2011;19:15221–8.2193488510.1364/OE.19.015221 Search in Google Scholar
 Ding F, Jin Y, Li B, Cheng H, Mo L, He S. Ultrabroadband strong light absorption based on thin multilayered metamaterials. Laser Photon Rev 2014;8:946–53.10.1002/lpor.201400157 Search in Google Scholar
 Hedayati MK, Javaherirahim M, Mozooni B, et al. Design of a perfect black absorber at visible frequencies using plasmonic metamaterials. Adv Mater 2011;23:5410–4.2199737810.1002/adma.201102646 Search in Google Scholar
 Chirumamilla M, Roberts AS, Ding F, et al. Multilayer tungsten-alumina-based broadband light absorbers for high-temperature applications. Opt Mater Express 2016;6:2704–14.10.1364/OME.6.002704 Search in Google Scholar
 Wang W, Qu Y, Du K, et al. Broadband optical absorption based on single-sized metal-dielectric-metal plasmonic nanostructures with high-ε″ metals. Appl Phys Lett 2017;110:101101.10.1063/1.4977860 Search in Google Scholar
 Kristensen A, Yang JKW, Bozhevolnyi SI, et al. Plasmonic colour generation. Nat Rev Mater 2016;2:16088. Search in Google Scholar
 Yang Y, Wang W, Moitra P, Kravchenko II, Briggs DP, Valentine J. Dielectric meta-reflectarray for broadband linear polarization conversion and optical vortex generation. Nano Lett 2014;14:1394–9.10.1021/nl404448224547692 Search in Google Scholar
 Ding F, Wang Z, He S, Shalaev VM, Kildishev AV. Broadband high-efficiency half-wave plate: a supercell-based plasmonic metasurface approach. ACS Nano 2015;9:4111–9.10.1021/acsnano.5b00218 Search in Google Scholar
 Markovich DL, Andryieuski A, Zalkovskij M, Malureanu R, Lavrinenko AV. Metamaterial polarization converter analysis: limits of performance. Appl Phys B 2013;112:143–52.10.1007/s00340-013-5383-8 Search in Google Scholar
 Chen Z, Wang C, Xu F, Lou Y, Cao B, Li X. Reflective plasmonic waveplates based on metal-insulator-metal subwavelength rectangular annular arrays. Photonics Nanostruct Fundam Appl 2014;12:189–98.10.1016/j.photonics.2014.01.001 Search in Google Scholar
 Jiang S-C, Xiong X, Hu Y-S, et al. Controlling the polarization state of light with a dispersion-free metastructure. Phys Rev X 2014;4:021026. Search in Google Scholar
 Hao J, Yuan Y, Ran L, et al. Manipulating electromagnetic wave polarizations by anisotropic metamaterials. Phys Rev Lett 2007;99:063908.10.1103/PhysRevLett.99.06390817930829 Search in Google Scholar
 Dai Y, Ren W, Cai H, Ding H, Pan N, Wang X. Realizing full visible spectrum metamaterial half-wave plates with patterned metal nanoarray/insulator/metal film structure. Opt Express 2014;22:7465–72.2471812110.1364/OE.22.007465 Search in Google Scholar
 Lévesque Q, Makhsiyan M, Bouchon P, et al. Plasmonic planar antenna for wideband and efficient linear polarization conversion. Appl Phys Lett 2014;104:111105.10.1063/1.4869127 Search in Google Scholar
 Ribaudo T, Taylor A, Nguyen B-M, Bethke D, Shaner EA. High efficiency reflective waveplates in the midwave infrared. Opt Express 2014;22:2821–9.2466357310.1364/OE.22.002821 Search in Google Scholar
 Guo Y, Wang Y, Pu M, et al. Dispersion management of anisotropic metamirror for super-octave bandwidth polarization conversion. Sci Rep 2015;5:8434.10.1038/srep0843425678280 Search in Google Scholar
 Yu N, Aieta F, Genevet P, Kats MA, Gaburro Z, Capasso F. A broadband, background-free quarter-wave plate based on plasmonic metasurfaces. Nano Lett 2012;12:6328–33.2313097910.1021/nl303445u Search in Google Scholar
 Shaltout A, Liu J, Kildishev A, Shalaev V. Photonic spin hall effect in gap plasmon metasurfaces for on-chip chiroptical spectroscopy. Optica 2015;2:860–3.10.1364/OPTICA.2.000860 Search in Google Scholar
 Maguid E, Yulevich I, Veksler D, Kleiner V, Brongersma ML, Hasman E. Photonic spin-controlled multifunctional shared-aperture antenna array. Science 2016;352:1202–6.10.1126/science.aaf341727103668 Search in Google Scholar
 Ding F, Pors A, Chen Y, Zenin VA, Bozhevolnyi SI. Beam-size-invariant spectropolarimeters using gap-plasmon metasurfaces. ACS Photonics 2017;4:943–9.10.1021/acsphotonics.6b01046 Search in Google Scholar
 Maier SA. Plasmonics: fundamentals and applications. Berlin, Springer Science & Business Media, 2007. Search in Google Scholar
 Qu C, Xiao S, Sun S, He Q, Zhou L. A theoretical study on the conversion efficiencies of gradient meta-surfaces. Europhys Lett 2013;101:54002.10.1209/0295-5075/101/54002 Search in Google Scholar
 Pors A, Nielsen MG, Bernardin T, Weeber JC, Bozhevolnyi SI. Efficient unidirectional polarization-controlled excitation of surface plasmon polaritons. Light Sci Appl 2014;3:e197.10.1038/lsa.2014.78 Search in Google Scholar
 Ding F, Deshpande R, Bozhevolnyi SI. Bifunctional gap-plasmon metasurfaces for visible light: polarization-controlled unidirectional surface plasmon excitation and beam steering at normal incidence. Light Sci Appl 2018;7:e17178.10.1038/lsa.2017.178 Search in Google Scholar
 Mühlenbernd H, Georgi P, Pholchai N, et al. Amplitude- and phase-controlled surface plasmon polariton excitation with metasurfaces. ACS Photonics 2016;3:124–9.10.1021/acsphotonics.5b00536 Search in Google Scholar
 Abad PC, Ou J-Y, Plum E, Zheludev N, et al. Electro-mechanical light modulator based on controlling the interaction of light with a metasurface. Sci Rep 2017;7:5405.10.1038/s41598-017-05906-928710432 Search in Google Scholar
 Liu L, Kang L, Mayer TS, Werner DH. Hybrid metamaterials for electrically triggered multifunctional control. Nat Commun 2016;7:13236.27807342 Search in Google Scholar
 Yao Y, Shankar R, Kats MA, et al. Electrically tunable metasurface perfect absorbers for ultrathin mid-infrared optical modulators. Nano Lett 2014;14:6526–32.10.1021/nl503104n25310847 Search in Google Scholar
 Emani NK, Kildishev AV, Shalaev VM, Boltasseva A. Graphene: a dynamic platform for electrical control of plasmonic resonance. Nanophotonics 2015;4:214–23. Search in Google Scholar
 Sherrott MC, Hon PWC, Fountaine KT, et al. Experimental demonstration of >230° phase modulation in gate-tunable graphene–gold reconfigurable mid-infrared metasurfaces. Nano Lett 2017;17:3027–34.10.1021/acs.nanolett.7b0035928445068 Search in Google Scholar
 Wang Z, Li T, Almdal K, Mortensen NA, Xiao S, Ndoni S. Experimental demonstration of graphene plasmons working close to the near-infrared window. Opt Lett 2016;41:5345–8.10.1364/OL.41.00534527842143 Search in Google Scholar
 Tittl A, Michel A-KU, Schäferling M, et al. A switchable mid-infrared plasmonic perfect absorber with multispectral thermal imaging capability. Adv Mater 2015;27:4597–603.10.1002/adma.20150202326173394 Search in Google Scholar
 Wang Q, Rogers ETF, Gholipour B, et al. Optically reconfigurable metasurfaces and photonic devices based on phase change materials. Nat Photon 2016;10:60–5.10.1038/nphoton.2015.247 Search in Google Scholar
 Kats MA, Blanchard R, Zhang S, et al. Vanadium dioxide as a natural disordered metamaterial: perfect thermal emission and large broadband negative differential thermal emittance. Phys Rev X 2013;3:041004. Search in Google Scholar
 Zhu Z, Evans PG, Haglund Jr RF, Valentine JG. Dynamically reconfigurable metadevice employing nanostructured phase-change materials. Nano Lett 2017;17:4881–5.10.1021/acs.nanolett.7b0176728731722 Search in Google Scholar
 Ding F, Zhong S, Bozhevolnyi SI. Vanadium dioxide integrated metasurfaces with switchable functionalities at terahertz frequencies. Adv Opt Mater 2018;1701204. https://onlinelibrary.wiley.com/doi/abs/10.1002/adom.201701204. Search in Google Scholar
 Park J, Kang J-H, Kim SJ, Liu X, Brongersma ML. Dynamic reflection phase and polarization control in metasurfaces. Nano Lett 2017;17:407–13.2793678410.1021/acs.nanolett.6b04378 Search in Google Scholar
 Xu H-X, Tang S, Ma S, et al. Tunable microwave metasurfaces for high-performance operations: dispersion compensation and dynamical switch. Sci Rep 2016;6:38255.10.1038/srep3825527901088 Search in Google Scholar
 Silva A, Monticone F, Castaldi G, Galdi V, Alù A, Engheta N. Performing mathematical operations with metamaterials. Science 2014;343:160–3.10.1126/science.124281824408430 Search in Google Scholar
 Chizari A, Abdollahramezani S, Jamali MV, Salehi JA. Analog optical computing based on a dielectric meta-reflect array. Opt lett 2016;41:3451–4.2747259110.1364/OL.41.003451 Search in Google Scholar
 Lee J, Tymchenko M, Argyropoulos C, et al. Giant nonlinear response from plasmonic metasurfaces coupled to intersubband transitions. Nature 2014;511:65–9.10.1038/nature1345524990746 Search in Google Scholar
 Estakhri NM, Alù A. Ultra-thin unidirectional carpet cloak and wavefront reconstruction with graded metasurfaces. IEEE Antennas Wireless Propag Lett 2014;13:1775–8.10.1109/LAWP.2014.2371894 Search in Google Scholar
 Orazbayev B, Estakhri NM, Alù A, Beruete M. Experimental demonstration of metasurface-based ultrathin carpet cloaks for millimeter waves. Adv Opt Mater 2017;5:1600606.10.1002/adom.201600606 Search in Google Scholar
 Yue F, Wen D, Xin J, Gerardot BD, Li J, Chen X. Vector vortex beam generation with a single plasmonic metasurface. ACS Photonics 2016;3:1558–63.10.1021/acsphotonics.6b00392 Search in Google Scholar
 Shen S, Meng L, Zhang Y, et al. Plasmon-enhanced second-harmonic generation nanorulers with ultrahigh sensitivities. Nano Lett 2015;15:6716–21.2637242510.1021/acs.nanolett.5b02569 Search in Google Scholar
 Butet J, Brevet P-F, Martin OJF. Optical second harmonic generation in plasmonic nanostructures: from fundamental principles to advanced applications. ACS Nano 2015;9:10545–62.10.1021/acsnano.5b0437326474346 Search in Google Scholar
 Wang F, Martinson ABF, Harutyunyan H. Efficient nonlinear metasurface based on nonplanar plasmonic nanocavities. ACS Photonics 2017;4:1188–94.10.1021/acsphotonics.7b00094 Search in Google Scholar
 Lee J, Nookala N, Gomez-Diaz JS, et al. Ultrathin second-harmonic metasurfaces with record-high nonlinear optical response. Adv Opt Mater 2016;4:664–70.10.1002/adom.201500723 Search in Google Scholar
 Cai T, Tang SW, Wang GM, et al. High-performance bifunctional metasurfaces in transmission and reflection geometries. Adv Opt Mater 2017;5:1600506.10.1002/adom.201600506 Search in Google Scholar