A Review of Graphene Plasmons and its Combination with Metasurface

Article information

J. Korean Ceram. Soc.. 2017;54(5):349-365
Publication date (electronic) : 2017 September 29
doi : https://doi.org/10.4191/kcers.2017.54.5.10
*Key Laboratory of Environmental Fracture (Ministry of Education), University of Science and Technology Beijing, Beijing 100083, China
**State Key Laboratory of New Ceramics and Fine Processing, Tsinghua University, Beijing 100084, China
***State Kay Laboratory of Tribology, Department of Mechanical Engineering, Tsinghua University, Beijing 100084, China
Corresponding author : Yang Bai, E-mail : baiy@mater.ustb.edu.cn, Tel : +86-10-6233-4493, Fax : +86-10-6233-4493
Received 2017 July 23; Revised 2017 September 5; Accepted 2017 September 6.

Abstract

Graphene has attracted a lot of attentions due to the unique electrical and optical properties. Compared with the noble metal plasmons in the visible and near-infrared frequencies, graphene can support surface plasmons in the lower frequencies of terahertz and mid-infrared and it demonstrates an extremely large confinement at the surface because of the particular electronic band structures. Especially, the surface conductivity of graphene can be tuned by either chemical doping or electrostatic gating. These features make graphene a promising candidate for plasmonics, biosensing and transformation optics. Furthermore, the combination of graphene and metasurfaces presents a powerful tunability for exotic electromagnetic properties, where the metasurfaces with the highly-localized fields offer a platform to enhance the interaction between the incident light and graphene and facilitate a deep modulation. In this paper, we provide an overview of the key properties of graphene, such as the surface conductivity, the propagating surface plasmon polaritons, and the localized surface plasmons, and the hybrid graphene/metasurfaces, either metallic and dielectric metasurfaces, from terahertz to near-infrared frequencies. Finally, there is a discussion for the current challenges and future goals.

1. Introduction

Graphene, a single layer of carbon atoms arranged in honeycomb structures, has attracted extensive attentions since the advent of free-standing graphene sheets in 2004.1) Because of the linear dispersion relation of Dirac electronic states, graphene demonstrates the unique electrical and optical properties,2,3) such as high carrier mobility, broadband operation, and a tunable band gap, which make it a fantastic material for photonics and optoelectronics.46) The optical response of graphene can be described by the surface conductivity which is originated from the contributions of inter-band and intra-band transitions, related to the chemical potential or Fermi level.7,8) In the near-infrared and visible frequencies, the graphene manifests a strong absorption caused by the inter-band transitions.7,9) In the terahertz and mid-infrared frequencies, the graphene behaves as an ultrathin metal film, supporting the surface plasmon polaritons (SPPs) with the greatly enhanced local fields at the graphene surface.1013) Compared with the conventional metal palsmons primarily working in the visible frequencies, the SPPs in graphene demonstrate a stronger confinement and lower loss, that is, an extremely short wavelength and a long propagation distance.12) Most important, SPPs in graphene can be dynamically tuned by the chemical doping14) or electrostatic gating.15) As a result, graphene has promising applications for high performance SPP-based nano-devices.16)

Due to the large mismatch between the wave vector of SPPs and the incident light, many methods, such as nanoemitters/tips scattering,12,13,1723) gratings2427) and prism coupling, 28,29) have been proposed. However, the coupling efficiency in the above methods is very low. Alternatively, the structured graphene (sub-wavelength nanoparticles) can demonstrate localized surface plasmons (LSPs) resonance with a strongly confined local field at the surface of nanoparticles. 3037) Compared with the propagating SPPs, LSPs resonance is a non-propagating regime and can be directly excited by a light with appropriate polarization and frequency. Apart from the strong light-matter interaction for sensing applications, the structured graphene is capable of modulating phase and amplitude used in polarization converter, 38,39) optical activity40) and wavefront reshaping applications. 4144)

The optical response for the monolayer graphene usually has a limited modulation depth due to the fact of single atomic thickness. To obtain a more powerful capacity of modulation, the hybrid graphene/metasurfaces have been proposed as a possible solution.4548) Metasurfaces are composed of sub-wavelength scatters which can interact with the electric or magnetic field components of incident light to produce the desired responses, including the phase shift, amplitude and polarization. By arranging the spatial distribution of the designed scatters treated as the tiny secondary sources, the metasurface can arbitrarily manipulate the electromagnetic wavefront, leading to myriad applications, such as holograms,4951) meta-lens5255) and abnormal refraction. 5658) Compared with the conventional optical components or bulky metamaterials59,60) depending on the electromagnetic wave propagating in the bulky media to gradually accumulate the changes of wavefront, the metasurfaces can realize the abrupt changes in a sub-wavelength-thin thickness, endowing the superior advantages of miniaturization, lightening and thinning. A drawback of conventional metasurfaces is non-tunability after they were fabricated. The versatile functionality of metasurfaces depends on the unit cells which are implemented by metallic patterns or dielectric cubes based on the plasmonic or Mie resonances.56,6163) It opens up an avenue to dynamically control light when the active graphene integrates with metasurfaces, which demonstrates more extensive application prospects, such as fast optical modulators,46,6467) switches45,47) and other optical devices.68,69) In addition, graphene is compatible with the nanofabrication technology and easy to be integrated into the nano/micro-sized metasurfaces.

This paper is organized as described below. In section 2, we overview the surface conductivity of graphene from terahertz to visible frequencies by the random-phase approximation (RPA) theory. In section 3, we introduce the propagating SPPs in graphene, including the characteristics, SPPs excitation, and SPPs hybridization. In section 4, we introduce LSPs in the structured graphene, modulation depth of absorption and electromagnetic applications, such as polarizer and meta-lens. In section 5, we overview the hybrid graphene/metallic metasurface working from terahertz to near-infrared frequencies, and the hybrid graphene/dielectric metasurface mainly working in the near-infrared frequencies. In the last section, we provide some concluding remarks and present an outlook on future research directions.

2. Surface Conductivity of Graphene

2.1. Surface conductivity at terahertz and infrared frequencies

The optical response of graphene is derived from the dynamic surface conductivity closely related to chemical potential or Fermi level. It can be derived with the random-phase approximation (RPA) theory8) and originate from the contribution of intra-band and inter-band transitions as below,

(1) σs(ω)=i·2e2kBTπ2(ω+iτ-1)ln[2cosh(EF2kBT)]+e24[12+1π+arctan(ω-EF2kBT)-i2πln(ω+EF)2(ω-EF)2+(2kBT)2]

where e is the electron charge, kB is the Boltzmann’s constant, ħ is the reduced Planck’s constant, T is the temperature, and ω is the angular frequency of incident light. τ is the relaxation time and estimated by expression τ=uEF/evf2, where vf = 106 m/s is the Fermi velocity, u is the carrier mobility of graphene and ranges from 1000 cm2/Vs in chemical vapor deposition (CVD)70,71) to 100 m2/Vs in suspended graphene.72) EF is the Fermi level which translates to the carrier concentration of n=EF2/π2vf2. In most cases, the Fermi energy EF is limited range of −1.0 eV to 1.0 eV as the consideration of device requirements to guarantee against the breakdown electric field of the gate dielectric.73) The first and second terms are attributed to the intra- and inter-band transitions, respectively.

Figure 1 is the calculated results based on Eq. (1). In the terahertz and mid-infrared frequencies where the photonic energy is far below the double Fermi level (ħω << 2|EF|), the surface conductivity of graphene are dominated by the intra-band transitions (Fig. 1(a)). The inter-band transitions are diminished due to the presence of electrons and holes near the band edges (Pauli blocking). The surface conductivity can be reduced to

Fig. 1

(a) Real parts of intra-band (solid) and inter-band (dash) surface conductivities of graphene with EF = 0.2 eV (black) and EF = 0.4 eV (red). The carrier mobility is assumed to be μ = 1000 cm2/Vs. The inset shows the band structure of graphene. (b) The optical response of moderately doped graphene in different spectral ranges. In the terahertz, near-infrared or midinfrared frequencies, graphene exhibits a lossy (|Im(σs)/Re(σs) < 1|) or plasmonic behavior (Im(σs)/Re(σs) > 1), respectively.

(2) σs(ω)=i·2e2kBTπ2(ω+iτ-1)ln[2cosh(EF2kBT)].

When the Fermi level is EF >> kBT, the surface conductivity is further simplified to

(3) σs(ω)=i·e2EFπ2(ω+iτ-1).

From Eqs. (2) and (3), the graphene just behaves as a conductive film and can be described by a classical Drude model. As a result, the graphene can support the transverse magnetic (TM) polarized SPPs. The surface conductivity presents a linear relationship with the Fermi level EF. Note that a considerable loss caused by the intra-band transitions is inevitable in the terahertz frequencies (Fig. 1(b)) which can be utilized to modulate amplitude by the absorption in graphene instead.

However, in the mid-infrared frequencies, the surface conductivity approaches the minimum due to the Pauli blocking and a relatively large incident frequency, degenerating the inter-band and intra-band transitions (Fig. 1(a)). The imaginary part of conductivity becomes dominant over the real part, facilitating to support the transverse magnetic SPPs with an extremely low loss (Fig. 1(b)). In the terahertz and mid-infrared frequencies, the conductivity of graphene has a largely tunable range which is facilitated to large-scalely tune the amplitude or frequency of the plasmon resonance when the graphene critically couples with the meta-surface of enhanced local field. Besides, the very thin “metal” layer can be further tailored to different shapes to obtain the desired resonances and reshape wavefront.

For the numerical simulations, graphene is often modeled as an anisotropic medium with in-plane effective permittivity as

(4) ɛeff,t=1+iσsɛ0ωtg=1-ωp2ω2

where ɛ0 is the vacuum permittivity, tg is the thickness of graphene, ωp=(2e2kBTπ2ɛ0tgln[2cosh (EF2kBT)])1/2 is the plasm frequency. The out-plane effective permittivity is equal to the vacuum permittivity of ɛeff,n = 1.0.

2.2. Surface conductivity at near-infrared and visible frequencies

On the contrary, the inter-band transitions become a decisive factor (ħω > 2|EF|) in the near-infrared and frequencies, as shown in Fig. 1(a). The free carrier scattering rate τ1 can be treated as zero due to the small quantity compared with the incident frequency.74) The real part of the surface conductivity is changed from 0 to e2/4ħ (universal optical conductivity) by electrically tuning Fermi level, which results in a modulation of absorption from 0 to 2.3% under the normal incidence. Obviously, the limited range of surface conductivity also limits the modulation depth of an isolated graphene, no matter for the transmission, reflection and absorption.75) To get a large modulation, the metasurface which can capture light from free space and offer a local filed enhancement in plane (two to three orders of magnitude larger than the incident field) is critical coupling with the tunable graphene. Besides, the imaginary part of the graphene permittivity becomes dominant over the real part so that the surface plasmons cannot be sustained at the graphene surface.10) In the visible frequencies, graphene can be visualized despite of only a single atom thick, due to a significant absorption πα = 2.3%, where α = 1/137 is the fine structure constant that describes coupling between light and electrons.

It is noted that the surface conductivity can also be tuned by the photonic excitation,76) magnetic field77) and strain mechanism,78) in addition to the common electrostatic gating and chemical doping.

3. Propagating Surface Plasmon Polaritons in Graphene

3.1. Properties of SPPs

The SPP is a collective oscillation of conduction electrons or holes along the boundary surface of a metal or semiconductor under the drive of electric field component of an incident light. It can highly confine the plasmonic wave and strongly enhance the electric field in the vicinity of surface, promising applications for the light-mater interactions, such as chemical, biological sensing and surface-enhanced Raman scattering.79) As stated above, the graphene behaves as a metal film in the terahertz and infrared frequencies, i.e. it is capable of supporting transverse magnetic polarized SPPs. For the simplest dielectric/continuous graphene sheet (the size larger than SPPs wavelength)/dielectric prototype, the dispersion relationship of SPPs is given by10,80)

(5) β(ω)=i(ɛr1+ɛr2)ɛ0ωσs(ω)=π2ɛ0(ɛr1+ɛr2)e2EFω(ω+iτ-1),

where β is the wave-vector of SPPs in graphene, ɛr1 and ɛr2 are the permittivity of dielectric above and below the graphene. Note that the precondition for the establishment of Eq. (5) is k0 << β << kF (in non-retarded regime and long wavelength limit), where k0 and kF are the vacuum and Fermi wave-vectors. For a finite β compared with kF, the dispersion of SPPs will deviate the Eq. (5) and present a linear relationship or damping regime.81) In addition, the dispersion of SPPs is also influenced by the polar medium, where the surface optical phonons couple with the SPPs via a long-range Frohlich coupling.82)

Compared with the conventional noble metal (silver and gold) plasmons, the graphene supported SPPs have many prominent properties. (1) Huge modal index or extensive confinement. For clarity, the Eq. (5) is simplified to

(6) β(ω)(ɛr1+ɛr2)ω4αEFk0,

where the scattering rate τ1 is assumed to zero, α = e2/4π ħɛ0c = 1/137 is the fine structure constant, the weight of ħω/EF is order of 1. Hence, the wave-vector of SPPs in graphene is two orders of magnitude bigger than the vacuum wave-vector.13,19,20) The wave-vector of SPPs in the noble metal, by contrast, is slightly larger than the vacuum wave-vector.83) (2) Relatively low losses or longer propagation. In the mid-infrared frequencies, the surface conductivity reaches the minimum, as shown in Fig. 1. The real part of permittivity in the graphene dominate the imaginary part, resulting in a relatively low damping loss of SPPs and a longer propagation length about hundreds of wavelength of SPPs, comparing to a range of tens of optical cycle for SPPs on metals. (3) Dynamic tunability. From Eq. (5), the propagation of SPPs in graphene is related to the Fermi level which can be tuned by electric field, magnetic field, electrostatic gating and chemical doping in real time. (4) Broadband operation. For the SPPs at the metal/dielectric surface which only works for the frequencies close to the plasma frequency of metal. At terahertz and infrared frequencies, the metal behaves as a perfect conductor with delocalized SPPs. On the contrary, the SPPs can be supported by the graphene from terahertz to near-infrared frequencies by controlling the Fermi level. These unique features make the graphene a promising platform for the active SPP-based nano-optoelectronic devices, plasmonics and transformation optical devices by spatially varying the graphene Fermi level in the terahertz and infrared frequencies. 84)

3.2. SPPs excitation

The major challenge is to efficiently excite SPPs in graphene with a free space incident light due to the large wave-vector mismatch between two waves. Conventional plasmonics approaches, such as the diffractive grating,2427) shifts the wave-vector of the incident light to compensate the mismatch due to Bragg scattering following the expression of β = sinθk0 + 2nπ/Λ, where n is the diffraction order, Λ is the periodicity, θ is the incident angle. Fig. 2(a) and (b) show the silicon-diffractive grating-assisted graphene plasmon excitation by guided-wave resonance.24) The Fourier transform infrared spectroscopy (FTIR) measurement result of T+/T//, where T+ and T// are the transmission spectrum under the normal incidence polarized perpendicular and parallel to the gratings lines, demonstrates a resonant dip due to the absorption of SPPs in the graphene (Fig. 2(c)). It is noted that, to efficiently excite the SPPs, a large coupling length is required (in the order of incident wavelength), which goes against the compactness of the nanodevices and bring a challenge to the nanofabrication.

Fig. 2

(a) Schematic of silicon-diffractive grating-assisted graphene plasmon excitation by guided-wave resonance (GWR). (b) Scanning electron microscopy (SEM) image of transferred graphene on a silicon grating, where the dark regions are the cracks. (c) Electrical tuning of the fabricated GWR graphene plasmonic device where arrows indicate the position of the resonant dips. (d) Schematic of the proposed coupling device, consisting of an Otto prism (couple the incident light to SPs), a tapered slab (compress the wavelength of SPs) and a thin slab (couple the SPs to SPPs in graphene). (e, f) Snapshots of the real part of the vertical component of the electric field Re(Ez) for SiC and n-GaAs waveguides. (a) – (f) used with permission from24,29) copyright by 2013, 2014 by American Chemical Society.

Another conventional method to excite the SPPs is using the prism coupling,28,29) such as the Otto or Kretschmann configurations,85,86) where the evanescent wave propagates along the metal-dielectric interface, due to the total internal reflection under the condition of the incident angle (θ) larger than the critical angle (θc). However, the huge difference of wave-vector between the evanescent wave and the incident light requires a high refractive index of the prism (np≈10–100), according to the equation β = k0npsinθ, which is not available practically. Nikitin et al.29) proposed a tactful configuration (Fig. 2(d)) to excite the SPPs in graphene based on the compression of surface polaritons (SPs) on the tapered slabs of a polaritonic material supporting either surface phonon polaritons (SPhPs) or SPPs, as shown in Fig. 2(e) and (f). The SPhPs or SPPs generating by the prism coupler propagate along the tapered slab and gradually yield a compression of wavelength to match the SPPs in graphene finally.

Besides the above methods, the excitation of SPPs is also proposed by coupling the near field of the nano-emitters with graphene.12,13,1723) Typically, the metallized tip of an atomic force microscope (AFM) illuminated by a focused infrared laser beam generates strong evanescent electric fields with a wide range of in-plane momenta which can effectively excite the graphene plasmons. In combination with the scattering-type scanning near-field optical microscopy (s-SNOM),13,19,20,23) the SPPs reflected or scattered by the external and internal edges, atomic defects, and boundaries of single-layer or multilayer graphene are visualized including amplitude and phase. Fig. 3(a) shows the illustration of the tip-induced SPPs excitation.20) Under the illumination of a CO2 laser with a wavelength of λ = 10.8 μm, the s-SNOM amplitude and phase exhibits the standing waves for SPPs due to the existence of external (transition from single and bilayer graphene) and internal (folds and grain boundaries) defects. The SPPs wavelength (λspp) in single and bilayer graphene are (260 ± 10) nm and (190 ± 10) nm, with a reflection coefficient of = −1, respectively (Fig. 3(b)). The SPPs in both sides of the line defects have different wavelength of 240 nm and 260 nm, as well as different reflection coefficients of = 0.45 and 0.55 (Fig. 3(c)). The ratio of λspp/λ0 is in the order of the fine structure constant. Besides, the change of Fermi level can effectively tune the SPPs wavelength. The full amplitude and phase spatially resolved imaging provides nano-scale insight into the local electronic structure and dynamics in graphene.

Fig. 3

(a) Illustration of metallic tip-induced SPP excitation and subsequent interference due to the emission of scattered and reflected SPP waves. (b, c) AFM topography (left), near-field amplitude (middle), and near-field phase (right) of a graphene wedge and a region of single-layer graphene with a high concentration of grain boundaries and folds. (d) Schematic of an antenna-based graphene SPPs launcher. (e, f) Experimental near-field images Re(Es,p) of a convex and a concave antenna extremity, respectively (λ0 = 11.06 μm). (g) Experimental near-field image Re(Es,p) of a bilayer graphene prism on a monolayer graphene (λ0 = 10.20 μm). (a) – (c) used with permission from20) copyright 2014 by American Physical Society, (d) – (g) used with permission from87) copyright 2014 by AAAS.

Figure 3(d) demonstrates the propagating SPPs in graphene excited by the dipole nano-antenna under the incident light (λ0 = 10.20 μm or 11.06 μm) polarized parallel to the long axis and imaged by s-SNOM with a dielectric Si-tip which passively scatter the near field of the graphene plasmon and gather the cross-polarization with a polarizer in front of the detector.87) The study shows the efficiency of resonance is enhanced by a factor of 28 compared with the nonresonant antennas. Besides, the wavefront of SPPs can be engineered by the shape of resonant antennas, such as concave and convex extremities (Fig. 3(e) and (f)). Further, the author demonstrated the propagating direction of SPPs can be controlled by a 2D graphene bilayer prism based on the conductivity patterns84) and revealed that the wavenumber in graphene is also influenced by the number of graphene layers (λspp ~ N1/2, N is the number of layers) which qualitatively followed the Snell’ law (Fig. 3(g)).37)

Furthermore, the SPPs in graphene can be visualized by the electron energy-loss (EEL) spectrum imaging22) at the atomic scale which offers a platform to study the effect of single atom on monolayer graphene sheets and brings the possibilities for designing atomic-scale optoelectronic and plasmonic devices.

3.3. SPPs hybridization

So far, we mainly review the excitation of SPPs in monolayer graphene sheet with a large size compared with the plasmon wavelength. Actually, the interaction of spatially separated graphene sheets is also an important subfield of graphene plasmon, because it is not only a fundamental issue but also provides a wider range of application prospects. Fig. 4(a) shows a graphene waveguide splitter.88) The propagating SPPs in the high doped graphene (ħω << EF) are completely transferred to the adjacent graphene by the weak SPPs coupling in the coupling region (top, Fig. 4(a)). By decreasing the Fermi level of the graphene in the upper or low coupling region to the condition of 2|EF| < ħω, where the graphene no longer supports SPPs (middle and bottom, Fig. 4(a)), the configuration realizes a tunable optical splitter or switch. Then, the author studied the strong coupling of SPPs in the monolayer graphene sheet arrays (MGSA) (Fig. 4(b)),89) which are physically different from the above case with weak coupling and demonstrate interesting properties (Fig. 4(c)). Beyond that, the SPPs in the monolayer graphene ribbon and the plasmonic hybridization of the spatially separated graphene ribbons are also intensively studied.90,91) The graphene ribbon supports an edge SPPs (E-SPPs) with a strongly localized field at the ribbon edges except the waveguide SPPs with the field concentration within the ribbon area. Especially, the graphene ribbon pairs demonstrate a more intensive enhancement of field.91)

Fig. 4

(a) Electric field (Ex, perpendicular to graphene surface) distribution of SPPs with bias on the output arm by a gate voltage: (top) no bias as splitter, upper arm (middle) or lower arm (bottom) biased for digital switch. (b) Schematic of a monolayer graphene sheet array (MGSA). (c) Intensity distributions of SPPs in the MGSA for different initial phase differences (ϕ0) of SPPs between neighboring graphene. (a) used with permission from88) copyright 2012 by AIP Publishing, (b) and (c) used with permission from89) copyright 2012 by American Physical Society.

Benefiting from the fascinating properties of SPPs in graphene in the terahertz and mid-infrared frequencies, the graphene opens up an exciting opportunity for the active sub-wavelength optics and the extremely compact optoelectronic devices.

4. Localized Surface Plasmons in Structured Graphene

4.1. Properties of LSPs

Different from the propagating SPPs, the LSPs are the non-propagating excitation of conductions of sub-wavelength metallic particles coupled with the incident light. It can be resonantly excited with a light of appropriate frequency and polarization, irrespective of the wave-vector match required in the propagating SPPs. At the plasmon resonant frequency, the LSPs are strongly confined to the surface of sub-wavelength particles and result in a significant electric field enhancement, facilitating to capture the incident light in solar cell and light-mater interaction for sensing applications. In the terahertz and mid-infrared frequencies, the noble metals are too metallic so that the metallic micro-structures suffer from a relatively poor field confinement, going against the non-destructive sensing. Alternatively, the graphene, as an excellent plasmonic material with an adjustable plasma frequency, can be structured into different shapes such as micro/nano-ribbons,3034) nano-tubes35) and disks36,37) which are capable of exploring the tunable LSPs in the terahertz and infrared frequencies.

Figure 5(a) and (b) are the graphene micro-ribbons32) and nano-ribbons34) tuned by the electrostatic electrolyte-gating and back-gating configurations working in the terahertz and mid-infrared frequencies, respectively. The extinction spectral in transmission is defined as 1−T/TR, where T is the transmission through the graphene ribbons with an electric field perpendicular to the ribbons, TR is the referenced transmission with different types, such as transmission through substrate without graphene ribbons,34) transmission through graphene ribbons at the charge neutral point (CNP)32) and transmission through graphene ribbons with electric fields parallel to the ribbons.33) The extinction spectral indicates a strong absorption and exhibits a prominent peak at the LSPs resonance. The frequency and strength of the maximum absorption can be tuned by the bias voltage, which clearly demonstrates a plasmon tunability through electrical gating, as shown in Fig. 5(c). In particular, when the duty ratio is below 0.5, the coupling of the graphene ribbons is negligible that the resonant wavelength of LSPs is identical to the case of isolated ribbon.92) The resonant wavelength of LSPs approximately satisfied

Fig. 5

(a) AFM and (b) SEM images of a graphene micro-ribbon array and a nano-ribbon array. (c) Control of terahertz resonance of plasmon excitations through electrical gating. Terahertz radiation was polarized perpendicular to the graphene ribbons. The inset shows extinction spectra for terahertz radiation polarized parallel to the ribbons. (d) Plasmon resonance frequency ωp as a function of |EF| or |n|1/2 for micro-ribbon arrays of different widths. (e) Normalized plasmon resonance frequency ωpw1/2 as a function of |EF| or |n|1/2 for micro-ribbon arrays of different widths. (f) SEM image (false color) of a stacked graphene/insulator micro-disk array arranged in a triangular lattice. (g) Extinction spectra for stacked plasmonic devices with one, two and five graphene layers under normal incidence. (h) Normalized plasmonic resonance frequency versus normalized total carrier density for single-layer graphene and graphene/insulator graphene stack. (i) Extinction spectra in single-layer graphene plasmonic devices with identical d (4.4 μm) but different a (4.8 and 9 μm, respectively). Inset: schematic of the graphene disk-disk interaction. (a) and (c) – (e) used with permission from32) copyright 2011 by Nature Publishing, (b) used with permission from34) copyright 2015 by AAAS, (f) – (i) used with permission from37) copyright 2012 by Nature Publishing.

(7) λp=2π/Re(β)2ω/(2m+1),

where w is the width of ribbon, m is the number of nodes with zero fields. Substituting Eq. (6) into Eq. (7), we can conclude the resonant wavelength is related to the ribbon width and the carrier concentration following the scaling behavior of λpw1/2|σs|1/2w1/2|EF|1/2w1/2n1/4 (Fig. 5(d) and 5(e)), which is different from the conventional two-dimensional electron gas (2DEG), where λpw1/2n1/2. Note that the LSPs can be further influenced by the shape of particles, the distance between particles and the number of layers.

However, in the monolayer graphene, LSPs resonance depending on the carrier concentration is usually very weak, which leads to a limited modulation depth of frequency and strength (ΔA ~ 10%) by varying the Fermi level of graphene. Recently, the multilayer graphene/insulator structures,37) shown in Fig. 5(f), demonstrates a larger tunability. The absorption strength is greatly enhanced (ΔA ~ 50%) and the resonant frequency is blue-shift with the increasing number of graphene disk (Fig. 5(g)). These phenomena are similar to the case of the monolayer graphene with the rise of carrier concentration or Fermi level. For the multilayer graphene/insulator stacks, the resonant wavelength of LSPs needs to be taken the layer number into consideration and scaled as λpN1/2 (Fig. 5(h)), leading to a relatively large tunable range compared with the monolayer graphene. In addition, the resonant wavelength of LSPs is significantly influenced by the Coulomb interactions among graphene disks, as shown in Fig. 5(i).

Although the multilayer graphene structure can realize a relatively large tunable range, the complicated preparation method limits the development of particle applications. Therefore, we will focus on the applications using the monolayer structured graphene in the following of this section.

4.2. Modulation depth of graphene absorption

A simple situation is first considered where the incident wave at an angle θ1 moves from a medium of refractive index n1 to a second medium of refractive index n2. A monolayer graphene with the surface conductivity of σs at the interface between two dielectrics is treated as a zero thickness conductive film. The transmission and reflection coefficients for TE and TM polarized waves can be expressed by the modified Fresnel equations:93)

(8) tTE=2cosθ1n1cosθ1n1+cosθ2n2+Z0σs
(9) rTE=cosθ1n1-cosθ2n2-Z0σscosθ1n1+cosθ2n2+Z0σs
(10) tTM=2cosθ1n1cosθ2n1+cosθ1n2+cosθ1cosθ2Z0σs
(11) rTM=cosθ2n1-cosθ1n2-cosθ1cosθ2Z0σscosθ2n1+cosθ1n2+cosθ1cosθ2Z0σs

where Z0 = 120π is the wave impedance of vacuum, θ2 = asin(sinθ1n1/n2) is the refractive angle. As a result, we get the power transmittance ( T=n2cosθ2n1cosθ1t2), reflectance (R = |r|2) and absorbance (A = 1−TR). The absorbance A is a function of surface conductivity (Re(σs) and Im(σs)). Based on the conditions of the extreme value for two variable functions, the maximum absorbances for TE and TM polarizations are94)

(12) ATEmax=11+Re(f),
(13) ATMmax=(ɛ1/ɛ2)f2Re(f)+(ɛ1/ɛ2)f2.

where f=(ɛ2/ɛ1)-sinθ12/cosθ1 and ɛ1,2=n1,22 is the relative permittivity. In the case of the normal incidence, Eqs. (12) and (13) are the same and equal to

(14) Amax=1/(1+n2/n1),

when the imaginary part of conductivity Im(σs) = 0 and the real part Re(σs) = (n1 + n2)/Z0 are satisfied simultaneously. From Eqs. (12) and (13), it is possible to achieve a perfect absorption, when the incident angle is larger than the critical angle (θ1 > asin(n2/n1), n1 > n2), where the transmission wave is suppressed and converted into the evanescent wave absorbed by the graphene. Note that the maximum of absorbance is just a theoretical upper limit which does not contains any clues about how to structure the materials or whether some material is possible to achieve it. For example, there is an upper limit absorbance of 50% (100% for coherent absorption95)) under the normal incidence with the symmetric configuration (i.e. n1 = n2 = 1.0).96,97) However, it is impossible for the continuous graphene with Amax = 50% in terahertz and infrared frequencies because the requirement of the surface conductivity including the real and imaginary parts cannot be satisfied simultaneously.

Alternatively, the structured graphene (ribbons, disks) or graphene/metasurface has more freedom degrees to tailor the surface conductivity, facilitating to realize the theoretical upper limit (Fig. 6(a)).94) The effective surface conductivity can be extracted by the S21 parameters from the numerical calculations

Fig. 6

(a) Schematic of the structured graphene sheet at the flat surface between two dielectrics. (b) Angular dependence of the peak absorbance for s- and p-polarization in arrays of graphene nano-disks. (c) Normal-incidence absorbance by graphene disk arrays supported on a dielectric-coated gold surface (inset) for various values of the dielectric film thickness. (d) Incidence-angle and photon-energy dependence of the absorbance under the conditions of (a) for d = 1.12 μm. (e) The schematic of a tunable metasurface absorber composed of an aluminum film, aluminum oxide layer, and a tunable metasurface on graphene. (f) Measured reflection spectra from the metasurface absorber for different gate voltages. (a) – (d) used with permission from94) copyright 2012 by American Physical Society. (e) and (f) used with permission from65) copyright 2014 by American Chemical Society.

(15) σseffZ0=2n1n2cosθ1cosθ2S21-(n1cosθ1+n2cosθ2)(TE polarization)
(16) σseffZ0=2n1n2S21cosθ1cosθ2-n1cosθ1-n2cosθ2(TM polarization).

where S21 parameter defined between two ports equals to tn2cosθ2/n1cosθ1. However, the theoretical study shows that the structured graphene only reach an absorption maximum in the limited incident angle due to the critical coupling (Fig. 6(b)).98) Besides, the light transferred from the vacuum to medium n1 undergoes a reflection loss in the actual situation that it is not a real perfect absorption.

To obtain the perfect absorption, the structured graphene or the graphene/metasurface on the top of a metallic ground plate separated by a thick dielectric spacer is adopted (Fig. 6(c) and (e)).65,99101) Based on the modified Fabry-Perot interference, we can get the total reflection coefficient as

(17) r=r1-t1t2exp(-2iϕ)1+r1exp(-2iϕ)

where r1 and t1(t2) are the reflection and transmission coefficients from medium n1(n2) to n2(n1) for TE or TM polarization, ϕ = k0n2d/cosθ2 is the phase accumulation upon the propagation through dielectric. Due to the thick metallic plate in the back of dielectrics, there is no transmittance (T = 0). Therefore, the absorbance is a function of the effective surface conductivity σseff and the phase accumulation ϕ. For example, the absorbance under the normal incidence reaches the maximum value when the conditions of Re(σseff)=n1/Z0 and Im(σseff)=-n2cot(ϕ)/Z0 are satisfied simultaneously which can be easily realized by adjusting the thickness of dielectrics (Fig. 6(c)).101) Importantly, the perfect absorber consisted of composite structure realizes an omnidirectional absorption (Fig. 6(d)) compared with the angle sensitive absorption from prism coupling (Fig. 6(b)). By varying the surface conductivity (or doping), the frequency and strength of the perfect absorber can be modulated in real time (Fig. 6(f)).65)

4.3. Electromagnetic applications

Except for the tunable ability of amplitude, the structured graphene is also capable to tune the phase and reshape the wavefront. For example, the graphene tailored into biaxial anisotropic nanoparticles can demonstrate different phase shifts along two principal axis due to the LSPs resonances. As a result, it can be utilized for polarization converter and optical activity.3840) Besides, the dynamic tunability of the effective surface conductivity by Fermi level further modulates the phase shift and amplitude for applications. For example, the asymmetric graphene nano-crosses (Fig. 7(a)) at Fermi level of 0.75 eV,40) demonstrate two distinct resonances with near 100° phase abrupt under the normally incidence with polarizations along x- and y-directions. At the wavelength of 7.92 μm, the transmission amplitude is equal and the phase difference approaches 90° between the x and y components, indicating a quarter wave plate under the linearly polarized wave with a polarization angle of 45°. Besides, the transmission phase difference is always 0 at 7.85 μm and 7.99 μm, indicating an optical activity. By varying the Fermi level, the structured graphene nano-crosses can realize a transition between a quarter wave plate and optical activity.

Fig. 7

(a) Graphene nano-crosses exhibit a giant optical activity (top & bottom) or a quarter wave plate (middle) at different wavelengths. (b, c) Anomalous refraction and flat lens composed of identical graphene bi-nanoribbon unit cells with different Fermi levels. (d) Anomalous reflection using graphene micro-ribbons/dielectric/metal with identical micro-ribbons but different Fermi levels and (e) An Airy beam generated by a graphene metasurface consisting of micro-ribbons/dielectric/metal with identical Fermi levels but different widths of micro-ribbons. (a) used with permission from40) copyright 2013 by The Optical Society, (b) and (c) used with permission from41) copyright 2014 by John Wiley and Sons, (d) used with permission from42) copyright 2015 by AIP Publishing, (e) used with permission from43) copyright 2015 by Nature Publishing.

In addition, the structured graphene scatters with the desired phase (varying geometric parameters or Fermi levels) can be properly arranged to reshape the wavefront and realize many functionalities in ultrathin devices, such as flat lens, anomalous reflection and refraction. The study shows that the transmission phase almost realizes a full coverage by controlling the Fermi level of the graphene ribbons for a normally incident light polarized perpendicular to the nano-ribbons.41) When two parallel graphene nano-ribbons are chosen as the unit cell, the transmission amplitude keeps uniform despite of the different Fermi levels. Hence, we can reshape the wavefront, i.e. anomalous refraction and lens (Fig. 7(b) and (c)) by arranging the distribution of nanoribbons with required Fermi levels. Besides, the structured graphene can be designed for the reflection modes42,43) (Fig. 7(d) and (e)) or the circular polarized wave based on the Pancharatnam-Berry (P-B) phase44) which is more suitable for broadband operation and multi-functions.

5. Graphene/Metasurfaces

5.1. Graphene/metallic metasurfaces

As mentioned above, the modulation depth of the graphene sheets, including frequency, amplitude and phase, usually exhibits a very limited tuning range due to the fact of a single atomic layer. One way to overcome this problem is using the structured graphene with the LSPs resonances in the THz and infrared frequencies, as mentioned above. Another effective solution is using the hybrid graphene/metallic metasurfaces because the metasurfaces offer an enhanced light-mater interaction from the terahertz to the visible frequencies. The main drawback of the metasurfaces is that they are passive, i.e. the optical properties are fixed after fabrication. Therefore, the hybrid graphene/metasurface follows from their complementarity and demonstrates a more powerful application.

In the terahertz frequencies, the surface conductivity of graphene is mainly derived from the intra-band transitions. The imaginary part of the permittivity dominates the real part, that is, the intrinsic losses caused by the intra-band transitions play an important role in the optical properties of graphene. By varying Fermi level, the surface conductivity of graphene demonstrates a largely tunable range, facilitating to realize a deep amplitude modulation when it is combined with the metallic metasurface of strongly localized in-plane field.47,48,66,101105) Fig. 8(a) shows the metal gratings with a high duty ratio, where the localized electric fields in the slits are greatly enhanced (Fig. 8(b)) and the frequency is independent under the normal incidence with polarization perpendicular to the gratings.103) The constructive interference between the coherent radiations from the oscillating local fields produces a high transmission. Therefore, the optical response of graphene/metallic gratings possess a broadband and deep modulation capacities (Fig. 8(c)). The change of transmission is mainly derived from the graphene absorption

Fig. 8

(a) Hybrid graphene/metasurface device by transferring a single layer graphene on top of a gold slit device. The incoming THz wave is polarized perpendicular to the slit orientation. (b) Electrical field enhancement factor in the slit for different surface conductivity of graphene. (c) Simulated (solid lines) and experimental (empty cycles) results of normalized THz power transmission modulation for bare graphene and the hybrid graphene-metasurface structure as a function of graphene conductivity at 1 THz. (d) SEM image of the fabricated Fano-resonant metasurface composed of asymmetric π-shape plasmonic resonators. (e) Simulated reflectivity spectra of hybrid graphene/metasurface for different Fermi levels. (f) Experimentally measured reflectivity spectra for the two samples with different geometrical dimensions before (solid curves) and after (dashed curves) graphene transfer. (g) A typical high resolution SEM image of the single gold nano-rod covered by graphene. (h) Scattering intensity (2D plot) is plotted as a function of the photon energy and gate voltage. (i) Four line cuts of the 2D plot for Rayleigh scattering spectra of an individual graphene-nano-rod hybrid structure. (a) – (i) used with permission from103,109,74) copyright 2015, 2013, 2012 by American Chemical Society.

(18) ΔT=ΔA~SEinc2η2ΔRe(σs)

where S is the graphene area, Einc is the electric field competent of incidence and η = E/Einc is the enhancement factor. From the Eq. (18), a large modulation depth can be obtained by greatly enhanced field despite of a little quality of surface conductivity. Beyond that, metallic resonators can exhibit strongly enhanced local fields for terahertz modulation by controlling the Fermi level of graphene, such as the splitring resonators (SRRs) with electric resonance,105) the sandwich structure with magnetic resonance66) and the asymmetric structures with Fano resonance.47,104)

In the mid-infrared frequencies, the graphene behave as a plasmonic material with a low loss, that is, the real part of the graphene permittivity dominates over the imaginary part. The optical response of the graphene/plasmonic metasurface is mainly reflected in the change of the resonant frequency by varying the Fermi level of graphene.45,65,67,106111) Fig. 8(d) shows the asymmetric π-shape plasmonic resonators which can exhibit high-sole asymmetric Fano resonances with a strong field enhancement when the incident polarization is parallel to the vertical resonators.109) The Fano resonance is derived from the interference between a narrow resonance (dark mode, quadrupole resonance) and a broad continuum of states (bright mode, dipole resonance), leading to an enhanced transmission and a reduced reflection.112,113) As shown in Fig. 8(e) and (f), the introduction of graphene to the metasurface significantly changes the frequency of Fano resonance, which can be evaluated by the perturbation theory and expressed as75,109)

(19) Δω=SΔσsimagEt2dsW0,

where W0=2ɛ0VɛE2dv is the stored energy of the given mode. In the mid-infrared frequencies, the surface conductivity has a largely tunable range, endowing the broadband modulation of resonant frequency and resulting in a deep modulation of amplitude or phase at a certain frequency which can be used for optical switch, modulator and photonic applications.

In the near-infrared frequencies, the inter-band transitions become the decisive factor and give rise to considerable losses which can be utilized for the amplitude modulation. As shown in Fig. 1(a), the surface conductivity of graphene has a maximum of e2/4ħ and produces a considerable losses of 2.3% at near-infrared frequencies. The limited varying range of the surface conductivity for isolated graphene sheet constrains the modulation depth. Therefore, the metasurfaces which offer a remarkable local field enhancement to the critical coupling with the tunable graphene are proposed. Fig. 8(g) shows a hybrid graphene/gold nano-rod device which can be efficiently controlled by the Fermi level due to the strongly confined the electric field in the vicinity of surface at the LSPs resonance.74) The Rayleigh scattering spectra, demonstrate the efficient electrical modulation of plasmon resonance (6%) at near infrared frequency utilizing the strong and tunable inter-band transitions in graphene (Fig. 8(h) and (i)). The change of amplitude can be obtained by75,114)

(20) ΔT-ΔA=-ΔσsrealSEt2dsc0ɛ0Einc2A,

where c0 is the speed of light, ɛ0 is the vacuum permittivity, A is the surface area of the unit cell, |Et| = η|Einc| is the tangent electric field in the graphene with an enhancement factor η compared with the incident electric Einc. To further improve the modulation depth in the near infrared frequency, the graphene should interact with high quality factor (Q-factor) resonances, such as the Fano resonances emerged in the asymmetric structures or plasmonic clusters. 46,68) Note that high Q-factor resonance usually accompany with a narrow bandwidth which demands the nearinfrared detector with a high resolution in reverse.

5.2. Graphene/dielectric metasurfaces

In the above content, we mainly focus on the hybrid graphene/metallic metasurfaces. In the near-infrared and optical frequencies, however, the non-radiative losses in metal lead to a low Q-factor resonance and a limited modulation depth.114116) Alternatively, the dielectric metasurfaces with high refractive index dielectric particles embedded in a low index matrix can support different order electric and magnetic resonances with an extremely low intrinsic losses based on the Mie theory.6163) For example, the asymmetric structures, such as asymmetric Si-based II-shape and π-shape resonators,117,118) demonstrate a high Q-factor Fano resonances and a strongly localized field enhancement in the gaps, facilitating to the depth modulation by varying the Fermi level (Fig. 9(a)). Similarly, the asymmetric periodity also has the same functions as the asymmetric structure.75) Fig. 9(b) shows an 85% modulation depth for the left-handed passband with 2π phase coverage when the Fermi level is varied from 0 to 0.8 eV. The modulation depth may be further improved by decreasing the asymmetric degree,119) but it also brings a challenge to a high machining precision and a higher Fermi level to realize high transmission, which makes the conventional gating configuration not suitable any more. Up to date, the experimental demonstration of dynamic tunability via graphene/dielectric metasurface has not yet been emerged.

Fig. 9

(a, b) Transmission spectra from the hybrid graphene/dielectric metasurface consisted of II-shape dielectric bars with structure asymmetry or dielectric blocks with periodity asymmetry for different Fermi levels. (c) SEM image of the hybrid graphene/dielectric waveguide modulator. (d) Static electro-optical response of the device for different gate voltages at 1.53 μm. (e) SEM image of the hybrid graphene/photonic crystal device. (f) Measured reflectivity from the hybrid device for different gate voltages. (a) used with permission from117) copyright 2015 by The Optical Society, (c) and (d) used with permission from16) copyright 2011 by Nature Publishing, (e) and (f) used with permission from123) copyright 2013 by American Chemical Society.

Another solution to realize a broadband or deep modulation is using dielectric waveguide or photonic crystal16,120124) which supports the guided mode or optical cavity resonance with the electric field confined in vicinity of surface. Fig. 9(c) shows the graphene/dielectric waveguide device which exhibits a large modulation depth of 0.1 dB μm1 (Fig. 9(d)) and a broadband modulation range from 1.35 μm to 1.6 μm.16) By introducing defects in the photonic crystal,123) as shown in Fig. 9(e), the localized field are greatly enhanced near the region of defects. The hybrid graphene/2D photonic crystal demonstrate a 6 dB modulation for the cavity resonance with 1.5 nm line width in reflectivity, as shown in Fig. 9(f).

In addition to the above mentioned applications of optical modulator or switch, the graphene integrated with nanoparticles with the greatly enhanced local fields are also promising applications for the optical bistability,125,126) super-scattering effect127129) and enhanced photo-detection.130,131)

6. Conclusions and Outlook

Graphene combines many superior properties from mechanical to electronic, has attracted intensive research interests from different fields. Owning to the high carrier mobility and tunable carrier concentration, graphene opens up an exciting opportunity to high performance applications, such as plasmonics, optoelectronic devices based on graphene plasmons. Different from noble metal plasmons in range of visible frequencies, graphene plasmons operate from terahertz to mid-infrared frequencies and demonstrate an extremely high confinement, facilitating to light-matter interaction for quantum optics, chemical identifications and biosensors. Importantly, the tunable surface conductivity makes graphene a wider application prospect. Graphene integrated with silicon-based optical circuit can dynamically manipulate the amplitude, phase and polarization of light, promising applications for optical communication and optical computing. Besides, graphene can couple with metasurfaces of enhanced local electric field to improve the optical modulation depth.

The behavior of graphene plasmons is tightly associated with the graphene itself. Large-scale grown high quality graphene is desirable in many applications and has been realized by chemical vapor deposition (CVD) on copper foils and films. However, the following transfer process is usually very complicated and need be optimized to minimize the damage to graphene. Another important challenge is the fabrication of structured graphene with fine control of shape and quality which will significantly influence the graphene plasmons.

Note that except for the graphene, many 2D materials, such as black phosphorus, molybdenum disulfide and topological insulators, also demonstrate the similar or even preferable electrical and optical properties which can be used to explore the tunable collective optical oscillations.

Acknowledgments

This work was supported by grants from the National Youth Top-Notch Talent Support Program and the Fundamental Research Funds for the Central Universities.

References

1. Novoselov KS, Geim AK, Morozov SV, Jiang D, Zhang Y, Dubonos SV, Grigorieva IV, Firsov AA. Electric Field Effect in Atomically Thin Carbon Films. Science 306(5696):666–69. 2004;
2. Neto AC, Guinea F, Peres NM, Novoselov KS, Geim AK. The Electronic Properties of Graphene. Rev Mod Phys 81(1):109. 2009;
3. Falkovsky L. Optical Properties of Graphene. J Phys: Conf Ser 129(1):012004. 2008;
4. Lemme MC, Echtermeyer TJ, Baus M, Kurz H. A Graphene Field-Effect Device. IEEE Electron Device Lett 28(4):282–84. 2007;
5. Bonaccorso F, Sun Z, Hasan T, Ferrari A. Graphene Photonics and Optoelectronics. Nat Photonics 4(9):611–22. 2010;
6. Lui CH, Mak KF, Shan J, Heinz TF. Ultrafast Photoluminescence from Graphene. Phys Rev Lett 105(12):127404. 2010;
7. Dawlaty JM, Shivaraman S, Strait J, George P, Chandrashekhar M, Rana F, Spencer MG, Veksler D, Chen Y. Measurement of the Optical Absorption Spectra of Epitaxial Graphene from Terahertz to Visible. Appl Phys Lett 93(13):131905. 2008;
8. Falkovsky L, Pershoguba S. Optical Far-Infrared Properties of a Graphene Monolayer and Multilayer. Phys Rev B 76(15):153410. 2007;
9. Nair RR, Blake P, Grigorenko AN, Novoselov KS, Booth TJ, Stauber T, Peres NM, Geim AK. Fine Structure Constant Defines Visual Transparency of Graphene. Science 320(5881):1308. 2008;
10. Jablan M, Buljan H, Soljačić M. Plasmonics in Graphene at Infrared Frequencies. Phys Rev B 80(24):245435. 2009;
11. Grigorenko A, Polini M, Novoselov K. Graphene Plasmonics. Nat Photonics 6(11):749–58. 2012;
12. Koppens FH, Chang DE, Garcia de Abajo FJ. Graphene Plasmonics: A Platform for Strong Light-Matter Interactions. Nano Lett 11(8):3370–77. 2011;
13. Chen J, Badioli M, Alonso-Gonzalez P, Thongrattanasiri S, Huth F, Osmond J, Spasenović M, Centeno A, Pesquera A, Godignon P. Optical Nano-Imaging of Gate-Tunable Graphene Plasmons. Nature 487(7405):77–81. 2012;
14. Kasry A, Kuroda MA, Martyna GJ, Tulevski GS, Bol AA. Chemical Doping of Large-Area Stacked Graphene Films for Use as Transparent, Conducting Electrodes. ACS Nano 4(7):3839–44. 2010;
15. Wang F, Zhang Y, Tian C, Girit C, Zettl A, Crommie M, Shen YR. Gate-Variable Optical Transitions in Graphene. Science 320(5873):206–9. 2008;
16. Liu M, Yin X, Ulin-Avila E, Geng B, Zentgraf T, Ju L, Wang F, Zhang X. A Graphene-based Broadband Optical Modulator. Nature 474(7349):64–7. 2011;
17. Nikitin AY, Guinea F, Garcia-Vidal F, Martin-Moreno L. Fields Radiated by a Nanoemitter in a Graphene Sheet. Phys Rev B 84(19):195446. 2011;
18. Bao Y, Zu S, Zhang Y, Fang Z. Active Control of Graphene-based Unidirectional Surface Plasmon Launcher. ACS Photonics 2(8):1135–40. 2015;
19. Fei Z, Rodin AS, Andreev GO, Bao W, McLeod AS, Wagner M, Zhang LM, Zhao Z, Thiemens M, Dominguez G, Fogler MM, Neto AHC, Lau CN, Keilmann F, Basov DN. Gate-Tuning of Graphene Plasmons Revealed by Infrared Nano-Imaging. Nature 487(7405):82–5. 2012;
20. Gerber JA, Berweger S, O’Callahan BT, Raschke MB. Phase-Resolved Surface Plasmon Interferometry of Graphene. Phys Rev Lett 113(5):055502. 2014;
21. Ni GX, Wang L, Goldflam MD, Wagner M, Fei Z, McLeod AS, Liu MK, Keilmann F, Ozyilmaz B, Castro Neto AH, Hone J, Fogler MM, Basov DN. Ultrafast Optical Switching of Infrared Plasmon Polaritons in High-Mobility Graphene. Nat Photonics 10(4):244–47. 2016;
22. Zhou W, Lee J, Nanda J, Pantelides ST, Pennycook SJ, Idrobo J-C. Atomically Localized Plasmon Enhancement in Monolayer Graphene. Nat Nanotechnol 7(3):161–65. 2012;
23. Fei Z, Rodin AS, Gannett W, Dai S, Regan W, Wagner M, Liu MK, McLeod AS, Dominguez G, Thiemens M, Castro Neto Antonio H, Keilmann F, Zettl A, Hillenbrand R, Fogler MM, Basov DN. Electronic and Plasmonic Phenomena at Graphene Grain Boundaries. Nat Nanotechnol 8(11):821–25. 2013;
24. Gao W, Shi G, Jin Z, Shu J, Zhang Q, Vajtai R, Ajayan PM, Kono J, Xu Q. Excitation and Active Control of Propagating Surface Plasmon Polaritons in Graphene. Nano Lett 13(8):3698–702. 2013;
25. Gao W, Shu J, Qiu C, Xu Q. Excitation of Plasmonic Waves in Graphene by Guided-Mode Resonances. ACS Nano 6(9):7806–13. 2012;
26. Lu H, Zhao J, Gu M. Nanowires-Assisted Excitation and Propagation of Mid-Infrared Surface Plasmon Polaritons in Graphene. J Appl Phys 120(16):163106. 2016;
27. Zhao Y, Chen G, Tao Z, Zhang C, Zhu Y. High Q-Factor Plasmonic Resonators in Continuous Graphene Excited by Insulator-Covered Silicon Gratings. RSC Adv 4(50):26535–42. 2014;
28. Bludov YV, Vasilevskiy M, Peres N. Mechanism for Graphene-based Optoelectronic Switches by Tuning Surface Plasmon-Polaritons in Monolayer Graphene. Europhys Lett 92(6):68001. 2011;
29. Nikitin AY, Alonso-Gonzalez P, Hillenbrand R. Efficient Coupling of Light to Graphene Plasmons by Compressing Surface Polaritons with Tapered Bulk Materials. Nano Lett 14(5):2896–901. 2014;
30. Brar VW, Jang MS, Sherrott M, Lopez JJ, Atwater HA. Highly Confined Tunable Mid-Infrared Plasmonics in Graphene Nanoresonators. Nano Lett 13(6):2541–47. 2013;
31. Hu H, Yang X, Zhai F, Hu D, Liu R, Liu K, Sun Z, Dai Q. Far-Field Nanoscale Infrared Spectroscopy of Vibrational Fingerprints of Molecules with Graphene Plasmons. Nat Commun 7:12334. 2016;
32. Ju L, Geng B, Horng J, Girit C, Martin M, Hao Z, Bechtel HA, Liang X, Zettl A, Shen YR. Graphene Plasmonics for Tunable Terahertz Metamaterials. Nat Nanotechnol 6(10):630–34. 2011;
33. Yan H, Low T, Zhu W, Wu Y, Freitag M, Li X, Guinea F, Avouris P, Xia F. Damping Pathways of Mid-Infrared Plasmons in Graphene Nanostructures. Nat Photonics 7(5):394–99. 2013;
34. Rodrigo D, Limaj O, Janner D, Etezadi D, de Abajo FJG, Pruneri V, Altug H. Mid-Infrared Plasmonic Biosensing with Graphene. Science 349(6244):165–68. 2015;
35. Garcia de Abajo FJ. Graphene Plasmonics: Challenges and Opportunities. ACS Photonics 1(3):135–52. 2014;
36. Fang Z, Thongrattanasiri S, Schlather A, Liu Z, Ma L, Wang Y, Ajayan PM, Nordlander P, Halas NJ, Garcia de Abajo FJ. Gated Tunability and Hybridization of Localized Plasmons in Nanostructured Graphene. ACS Nano 7(3):2388–95. 2013;
37. Yan H, Li X, Chandra B, Tulevski G, Wu Y, Freitag M, Zhu W, Avouris P, Xia F. Tunable Infrared Plasmonic Devices Using Graphene/Insulator Stacks. Nat Nanotechnol 7(5):330–34. 2012;
38. Guo T, Argyropoulos C. Broadband Polarizers based on Graphene Metasurfaces. Opt Lett 41(23):5592–95. 2016;
39. Fallahi A, Perruisseau-Carrier J. Design of Tunable Biperiodic Graphene Metasurfaces. Phys Rev B 86(19):195408. 2012;
40. Cheng H, Chen S, Yu P, Li J, Deng L, Tian J. Mid-Infrared Tunable Optical Polarization Converter Composed of Asymmetric Graphene Nanocrosses. Opt Lett 38(9):1567–69. 2013;
41. Lu F, Liu B, Shen S. Infrared Wavefront Control based on Graphene Metasurfaces. Adv Opt Mater 2(8):794–99. 2014;
42. Yatooshi T, Ishikawa A, Tsuruta K. Terahertz Wavefront Control by Tunable Metasurface Made of Graphene Ribbons. Appl Phys Lett 107(5):053105. 2015;
43. Li Z, Yao K, Xia F, Shen S, Tian J, Liu Y. Graphene Plasmonic Metasurfaces to Steer Infrared Light. Sci Rep 5:12423. 2015;
44. Hu X, Wen L, Song S, Chen Q. Tunable Graphene Metasurfaces by Discontinuous Pancharatnam-Berry Phase Shift. Nanotechnology 26(50):505203. 2015;
45. Dabidian N, Kholmanov I, Khanikaev AB, Tatar K, Trendafilov S, Mousavi SH, Magnuson C, Ruoff RS, Shvets G. Electrical Switching of Infrared Light Using Graphene Integration with Plasmonic Fano Resonant Metasurfaces. ACS Photonics 2(2):216–27. 2015;
46. Emani NK, Chung T-F, Kildishev AV, Shalaev VM, Chen YP, Boltasseva A. Electrical Modulation of Fano Resonance in Plasmonic Nanostructures Using Graphene. Nano Lett 14(1):78–82. 2013;
47. Lee SH, Choi M, Kim T-T, Lee S, Liu M, Yin X, Choi HK, Lee SS, Choi C-G, Choi S-Y. Switching Terahertz Waves with Gate-Controlled Active Graphene Metamaterials. Nat Mater 11(11):936–41. 2012;
48. Jadidi MM, Sushkov AB, Myers-Ward RL, Boyd AK, Daniels KM, Gaskill DK, Fuhrer MS, Drew HD, Murphy TE. Tunable Terahertz Hybrid Metal-Graphene Plasmons. Nano Lett 15(10):7099–104. 2015;
49. Ni X, Kildishev AV, Shalaev VM. Metasurface Holograms for Visible Light. Nat Commun 4:2807. 2013;
50. Zheng G, Muhlenbernd H, Kenney M, Li G, Zentgraf T, Zhang S. Metasurface Holograms Reaching 80% Efficiency. Nat Nanotechnol 10(4):308–12. 2015;
51. Wen D, Yue F, Li G, Zheng G, Chan K, Chen S, Chen M, Li KF, Wong PWH, Cheah KW, Yue Bun Pun E, Zhang S, Chen X. Helicity Multiplexed Broadband Metasurface Holograms. Nat Commun 6:8241. 2015;
52. Ni X, Ishii S, Kildishev AV, Shalaev VM. Ultra-Thin, Planar, Babinet-Inverted Plasmonic Metalenses. Light Sci Appl 2(4):e72. 2013;
53. Aieta F, Genevet P, Kats MA, Yu N, Blanchard R. Aberration-free Ultrathin Flat Lenses and Axicons at Telecom Wavelengths based on Plasmonic Metasurfaces. Nano Lett 12:4932–36. 2012;
54. Kang M, Feng T, Wang H-T, Li J. Wave Front Engineering from an Array of Thin Aperture Antennas. Opt Express 20(14):15882–90. 2012;
55. Liu C, Bai Y, Jing L, Yang Y, Chen H, Zhou J, Zhao Q, Qiao L. Equivalent Energy Level Hybridization Approach for High-Performance Metamaterials Design. Acta Mater 135(15):144–49. 2017;
56. Yu N, Genevet P, Kats MA, Aieta F, Tetienne J-P, Capasso F, Gaburro Z. Light Propagation with Phase Discontinuities: Generalized Laws of Reflection and Refraction. Science 334(6054):333–37. 2011;
57. Liu C, Bai Y, Zhao Q, Yang Y, Chen H, Zhou J, Qiao L. Fully Controllable Pancharatnam-Berry Metasurface Array with High Conversion Efficiency and Broad Bandwidth. Sci Rep 6:34819. 2016;
58. Ni X, Emani NK, Kildishev AV, Boltasseva A, Shalaev VM. Broadband Light Bending with Plasmonic Nanoantennas. Science 335(6067):427. 2012;
59. Smith DR, Padilla WJ, Vier DC, Nemat-Nasser SC, Schultz S. Composite Medium with Simultaneously Negative Permeability and Permittivity. Phys Rev Lett 84(18):4184–87. 2000;
60. Schurig D, Mock JJ, Justice BJ, Cummer SA, Pendry JB, Starr AF, Smith DR. Metamaterial Electromagnetic Cloak at Microwave Frequencies. Science 314(5801):977–80. 2006;
61. 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 14(3):1394–99. 2014;
62. Chong KE, Staude I, James A, Dominguez J, Liu S, Campione S, Subramania GS, Luk TS, Decker M, Neshev DN. Polarization-Independent Silicon Metadevices for Efficient Optical Wavefront Control. Nano Lett 15(8):5369–74. 2015;
63. Shalaev MI, Sun J, Tsukernik A, Pandey A, Nikolskiy K, Litchinitser NM. High-Efficiency All-Dielectric Metasurfaces for Ultracompact Beam Manipulation in Transmission Mode. Nano Lett 15(9):6261–66. 2015;
64. Gao W, Shu J, Reichel K, Nickel DV, He X, Shi G, Vajtai R, Ajayan PM, Kono J, Mittleman DM. High-Contrast Terahertz Wave Modulation by Gated Graphene Enhanced by Extraordinary Transmission through Ring Apertures. Nano Lett 14(3):1242–48. 2014;
65. Yao Y, Shankar R, Kats MA, Song Y, Kong J, Loncar M, Capasso F. Electrically Tunable Metasurface Perfect Absorbers for Ultrathin Mid-Infrared Optical Modulators. Nano Lett 14(11):6526–32. 2014;
66. Miao Z, Wu Q, Li X, He Q, Ding K, An Z, Zhang Y, Zhou L. Widely Tunable Terahertz Phase Modulation with Gate-Controlled Graphene Metasurfaces. Phys Rev X 5(4):041027. 2015;
67. Dabidian N, Dutta-Gupta S, Kholmanov I, Lai K, Lu F, Lee J, Jin M, Trendafilov S, Khanikaev A, Fallahazad B. Experimental Demonstration of Phase Modulation and Motion Sensing Using Graphene-Integrated Metasurfaces. Nano Lett 16(6):3607–15. 2016;
68. Fang Z, Liu Z, Wang Y, Ajayan PM, Nordlander P, Halas NJ. Graphene-Antenna Sandwich Photo-detector. Nano Lett 12(7):3808–13. 2012;
69. Fan Y, Shen NH, Zhang F, Wei Z, Li H, Zhao Q, Fu Q, Zhang P, Koschny T, Soukoulis CM. Electrically Tunable Goos-Hanchen Effect with Graphene in the Terahertz Regime. Adv Opt Mater 4(11):1824–28. 2016;
70. Kim KS, Zhao Y, Jang H, Lee SY, Kim JM, Kim KS, Ahn J-H, Kim P, Choi J-Y, Hong BH. Large-Scale Pattern Growth of Graphene Films for Stretchable Transparent Electrodes. Nature 457(7230):706. 2009;
71. Li X, Magnuson CW, Venugopal A, Tromp RM, Hannon JB, Vogel EM, Colombo L, Ruoff RS. Large-Area Graphene Single Crystals Grown by Low-Pressure Chemical Vapor Deposition of Methane on Copper. J Am Chem Soc 133(9):2816–19. 2011;
72. Castro EV, Ochoa H, Katsnelson M, Gorbachev R, Elias D, Novoselov K, Geim A, Guinea F. Limits on Charge Carrier Mobility in Suspended Graphene due to Flexural Phonons. Phys Rev Lett 105(26):266601. 2010;
73. Chen C-F, Park C-H, Boudouris BW, Horng J, Geng B, Girit C, Zettl A, Crommie MF, Segalman RA, Louie SG, Wang F. Controlling Inelastic Light Scattering Quantum Pathways in Graphene. Nature 471(7340):617–20. 2011;
74. Kim J, Son H, Cho DJ, Geng B, Regan W, Shi S, Kim K, Zettl A, Shen Y-R, Wang F. Electrical Control of Optical Plasmon Resonance with Graphene. Nano Lett 12(11):5598–602. 2012;
75. Liu C, Bai Y, Zhou J, Zhao Q, Qiao L. Large-Scale Modulation of Left-Ganded Passband in Gybrid Graphene/Dielectric Metasurface. Ann Phys 529(8):1700125. 2017;
76. Fang Z, Wang Y, Liu Z, Schlather A, Ajayan PM, Koppens FH, Nordlander P, Halas NJ. Plasmon-Induced Doping of Graphene. ACS Nano 6(11):10222–28. 2012;
77. Yan H, Li Z, Li X, Zhu W, Avouris P, Xia F. Infrared Spectroscopy of Tunable Dirac Terahertz Magneto-Plasmons in Graphene. Nano Lett 12(7):3766–71. 2012;
78. Verma VP, Das S, Lahiri I, Choi W. Large-Area Graphene on Polymer Film for Flexible and Transparent Anode in Field Emission Device. Appl Phys Lett 96(20):203108. 2010;
79. Schedin F, Lidorikis E, Lombardo A, Kravets VG, Geim AK, Grigorenko AN, Novoselov KS, Ferrari AC. Surface-Enhanced Raman Spectroscopy of Graphene. ACS Nano 4(10):5617–26. 2010;
80. Mikhailov S, Ziegler K. New Electromagnetic Mode in Graphene. Phys Rev Lett 99(1):016803. 2007;
81. Liu Y, Willis RF, Emtsev K, Seyller T. Plasmon Dispersion and Damping in Electrically Isolated Two-Dimensional Charge Sheets. Phys Rev B 78(20):201403. 2008;
82. Fei Z, Andreev GO, Bao W, Zhang LM, McLeod AS, Wang C, Stewart MK, Zhao Z, Dominguez G, Thiemens M. Infrared Nanoscopy of Dirac Plasmons at the Graphene-SiO2 Interface. Nano Lett 11(11):4701–5. 2011;
83. Barnes WL, Dereux A, Ebbesen TW. Surface Plasmon Subwavelength Optics. Nature 424(6950):824. 2003;
84. Vakil A, Engheta N. Transformation Optics Using Graphene. Science 332(6035):1291–94. 2011;
85. Otto A. Excitation of Nonradiative Surface Plasma Waves in Silver by the Method of Frustrated Total Reflection. Z Phys 216(4):398–410. 1968;
86. Kretschmann E, Raether H. Radiative Decay of Non Radiative Surface Plasmons Excited by Light. Z Naturforsch A 23(12):2135–36. 1968;
87. Alonso-Gonzalez P, Nikitin AY, Golmar F, Centeno A, Pesquera A, Velez S, Chen J, Navickaite G, Koppens F, Zurutuza A. Controlling Graphene Plasmons with Resonant Metal Antennas and Spatial Conductivity Patterns. Science 344(6190):1369–73. 2014;
88. Wang B, Zhang X, Yuan X, Teng J. Optical Coupling of Surface Plasmons between Graphene Sheets. Appl Phys Lett 100(13):131111. 2012;
89. Wang B, Zhang X, Garcia-Vidal FJ, Yuan X, Teng J. Strong Coupling of Surface Plasmon Polaritons in Monolayer Graphene Sheet Arrays. Phys Rev Lett 109(7):073901. 2012;
90. Nikitin AY, Guinea F, Garcia-Vidal FJ, Martin-Moreno L. Edge and Waveguide Terahertz Surface Plasmon Modes in Graphene Microribbons. Phys Rev B 84(16):161407. 2011;
91. Christensen J, Manjavacas A, Thongrattanasiri S, Koppens FH, Garcia de Abajo FJ. Graphene Plasmon Waveguiding and Hybridization in Individual and Paired Nanoribbons. ACS Nano 6(1):431–40. 2011;
92. Nikitin AY, Guinea F, Garcia-Vidal FJ, Martin-Moreno L. Surface Plasmon Enhanced Absorption and Suppressed Transmission in Periodic Arrays of Graphene ribbons. Phys Rev B 85(8):081405. 2012;
93. Sensale-Rodriguez B, Yan R, Kelly MM, Fang T, Tahy K, Hwang WS, Jena D, Liu L, Xing HG. Broadband Graphene Terahertz Modulators Enabled by Intraband Transitions. Nat Commun 3:780. 2012;
94. Thongrattanasiri S, Koppens FH, De Abajo FJG. Complete Optical Absorption in Periodically Patterned Graphene. Phys Rev Lett 108(4):047401. 2012;
95. Fan Y, Zhang F, Zhao Q, Wei Z, Li H. Tunable Terahertz Coherent Perfect Absorption in a Monolayer Graphene. Opt Lett 39(21):6269–72. 2014;
96. Fan Y, Shen N-H, Koschny T, Soukoulis CM. Tunable Terahertz Meta-Surface with Graphene Cut-Wires. ACS Photonics 2(1):151–56. 2015;
97. Fan Y, Wei Z, Zhang Z, Li H. Enhancing Infrared Extinction and Absorption in a Monolayer Graphene Sheet by Harvesting the Electric Dipolar Mode of Split Ring Resonators. Opt Lett 38(24):5410–13. 2013;
98. Teperik TV, De Abajo FG, Borisov A, Abdelsalam M, Bartlett P, Sugawara Y, Baumberg J. Omnidirectional Absorption in Nanostructured Metal Surfaces. Nature Photon 2(5):299–301. 2008;
99. Jang MS, Brar VW, Sherrott MC, Lopez JJ, Kim L, Kim S, Choi M, Atwater HA. Tunable Large Resonant Absorption in a Midinfrared Graphene Salisbury Screen. Phys Rev B 90(16):165409. 2014;
100. Alaee R, Farhat M, Rockstuhl C, Lederer F. A Perfect Absorber Made of a Graphene Micro-Ribbon Metamaterial. Opt Express 20(27):28017–24. 2012;
101. Andryieuski A, Lavrinenko AV. Graphene Meta-materials based Tunable Terahertz Absorber: Effective Surface Conductivity Approach. Opt Express 21(7):9144–55. 2013;
102. Chen P-Y, Huang H, Akinwande D, Alu A. Graphene-based Plasmonic Platform for Reconfigurable Terahertz Nanodevices. ACS Photonics 1(8):647–54. 2014;
103. Shi SF, Zeng B, Han HL, Hong X, Tsai HZ, Jung HS, Zettl A, Crommie MF, Wang F. Optimizing Broadband Terahertz Modulation with Hybrid Graphene/Metasurface Structures. Nano Lett 15(1):372–77. 2015;
104. Li Q, Cong L, Singh R, Xu N, Cao W, Zhang X, Tian Z, Du L, Han J, Zhang W. Monolayer Graphene Sensing Enabled by the Strong Fano-Resonant Metasurface. Nanoscale 8(39):17278–84. 2016;
105. Li Q, Tian Z, Zhang X, Xu N, Singh R, Gu J, Lv P, Luo L-B, Zhang S, Han J, Zhang W. Dual Control of Active Graphene-Silicon Hybrid Metamaterial Devices. Carbon 90:146–53. 2015;
106. Emani NK, Chung T-F, Ni X, Kildishev AV, Chen YP, Boltasseva A. Electrically Tunable Damping of Plasmonic Resonances with Graphene. Nano Lett 12(10):5202–6. 2012;
107. Yao Y, Kats MA, Genevet P, Yu N, Song Y, Kong J, Capasso F. Broad Electrical Tuning of Graphene-Loaded Plasmonic Antennas. Nano Lett 13(3):1257–64. 2013;
108. Yao Y, Shankar R, Rauter P, Song Y, Kong J, Loncar M, Capasso F. High-Responsivity Mid-Infrared Graphene Detectors with Antenna-Enhanced Photocarrier Generation and Collection. Nano Lett 14(7):3749–54. 2014;
109. Mousavi SH, Kholmanov I, Alici KB, Purtseladze D, Arju N, Tatar K, Fozdar DY, Suk JW, Hao Y, Khanikaev AB. Inductive Tuning of Fano-Resonant Metasurfaces Using Plasmonic Response of Graphene in the Mid-Infrared. Nano Lett 13(3):1111–17. 2013;
110. Smirnova DA, Miroshnichenko AE, Kivshar YS, Khanikaev AB. Tunable Nonlinear Graphene Metasurfaces. Phys Rev B 92(16):161406. 2015;
111. Yao Y, Kats MA, Shankar R, Song Y, Kong J, Loncar M, Capasso F. Wide Wavelength Tuning of Optical Antennas on Graphene with Nanosecond Response Time. Nano Lett 14(1):214–19. 2013;
112. Luk’yanchuk B, Zheludev NI, Maier SA, Halas NJ, Nordlander P, Giessen H, Chong CT. The Fano Resonance in Plasmonic Nanostructures and Metamaterials. Nature Mater 9(9):707–15. 2010;
113. Miroshnichenko AE, Flach S, Kivshar YS. Fano Resonances in Nanoscale Structures. Rev Mod Phys 82(3):2257. 2010;
114. Zhao B, Zhao J, Zhang Z. Enhancement of Near-Infrared Absorption in Graphene with Metal Gratings. Appl Phys Lett 105(3):031905. 2014;
115. Boltasseva A, Atwater HA. Low-Loss Plasmonic Metamaterials. Science 331(6015):290–91. 2011;
116. West PR, Ishii S, Naik GV, Emani NK, Shalaev VM, Boltasseva A. Searching for Better Plasmonic Materials. Laser Photonics Rev 4(6):795–808. 2010;
117. Argyropoulos C. Enhanced Transmission Modulation Based on Dielectric Metasurfaces Loaded with Graphene. Opt Express 23(18):23787–97. 2015;
118. Wu C, Arju N, Kelp G, Fan JA, Dominguez J, Gonzales E, Tutuc E, Brener I, Shvets G. Spectrally Selective Chiral Silicon Metasurfaces based on Infrared Fano Resonances. Nat Commun 27(5):3892. 2014;
119. Zhang J, Liu W, Zhu Z, Yuan X, Qin S. Strong Field Enhancement and Light-Matter Interactions with All-Dielectric Metamaterials based on Split Bar Resonators. Opt Express 22(25):30889–98. 2014;
120. Wu S, Buckley S, Jones AM, Ross JS, Ghimire NJ, Yan J, Mandrus DG, Yao W, Hatami F, Vučković J. Control of Two-Dimensional Excitonic Light Emission via Photonic Crystal. 2D Mater 1(1):011001. 2014;
121. Li R, Imran M, Lin X, Wang H, Xu Z, Chen H. Hybrid Airy Plasmons with Dynamically Steerable Trajectories. Nanoscale 9(4):1449–56. 2017;
122. Shiue R-J, Gan X, Gao Y, Li L, Yao X, Szep A, Walker D Jr, Hone J, Englund D. Enhanced Photodetection in Graphene-Integrated Photonic Crystal Cavity. Appl Phys Lett 103(24):241109. 2013;
123. Majumdar A, Kim J, Vuckovic J, Wang F. Electrical Control of Silicon Photonic Crystal Cavity by Graphene. Nano Lett 13(2):515–18. 2013;
124. Fan Y, Wei Z, Li H, Chen H, Soukoulis CM. Photonic Band Gap of a Graphene-Embedded Quarter-Wave Stack. Phys Rev B 88(24):241403. 2013;
125. Li R, Wang H, Zheng B, Dehdashti S, Li E, Chen H. Bistable Scattering in Graphene-Coated Dielectric Nanowires. Nanoscale 9:8449–57. 2017;
126. Li R, Lin X, Lin S, Zhang X, Li E, Chen H. Graphene Induced Mode Bifurcation at Low Input Power. Carbon 98:463–67. 2016;
127. Li R, Lin X, Lin S, Liu X, Chen H. Tunable Deep-Subwavelength Superscattering Using Graphene Monolayers. Opt Lett 40(8):1651–54. 2015;
128. Li R, Lin X, Lin S, Liu X, Chen H. Atomically Thin Spherical Shell-Shaped Superscatterers based on a Bohr Model. Nanotechnology 26(50):505201. 2015;
129. Li R, Zheng B, Lin X, Hao R, Lin S, Yin W, Li E, Chen H. Design of Ultracompact Graphene-Based Superscatterers. IEEE J Sel Top Quantum Electron 23(1):130–37. 2017;
130. Liu Y, Cheng R, Liao L, Zhou H, Bai J, Liu G, Liu L, Huang Y, Duan X. Plasmon Resonance Enhanced Multicolour Photodetection by Graphene. Nat Commun 2:579. 2011;
131. Echtermeyer TJ, Britnell L, Jasnos PK, Lombardo A, Gorbachev RV, Grigorenko AN, Geim AK, Ferrari AC, Novoselov KS. Strong Plasmonic Enhancement of Photovoltage in Graphene. Nat Commun 2:458. 2011;

Article information Continued

Fig. 1

(a) Real parts of intra-band (solid) and inter-band (dash) surface conductivities of graphene with EF = 0.2 eV (black) and EF = 0.4 eV (red). The carrier mobility is assumed to be μ = 1000 cm2/Vs. The inset shows the band structure of graphene. (b) The optical response of moderately doped graphene in different spectral ranges. In the terahertz, near-infrared or midinfrared frequencies, graphene exhibits a lossy (|Im(σs)/Re(σs) < 1|) or plasmonic behavior (Im(σs)/Re(σs) > 1), respectively.

Fig. 2

(a) Schematic of silicon-diffractive grating-assisted graphene plasmon excitation by guided-wave resonance (GWR). (b) Scanning electron microscopy (SEM) image of transferred graphene on a silicon grating, where the dark regions are the cracks. (c) Electrical tuning of the fabricated GWR graphene plasmonic device where arrows indicate the position of the resonant dips. (d) Schematic of the proposed coupling device, consisting of an Otto prism (couple the incident light to SPs), a tapered slab (compress the wavelength of SPs) and a thin slab (couple the SPs to SPPs in graphene). (e, f) Snapshots of the real part of the vertical component of the electric field Re(Ez) for SiC and n-GaAs waveguides. (a) – (f) used with permission from24,29) copyright by 2013, 2014 by American Chemical Society.

Fig. 3

(a) Illustration of metallic tip-induced SPP excitation and subsequent interference due to the emission of scattered and reflected SPP waves. (b, c) AFM topography (left), near-field amplitude (middle), and near-field phase (right) of a graphene wedge and a region of single-layer graphene with a high concentration of grain boundaries and folds. (d) Schematic of an antenna-based graphene SPPs launcher. (e, f) Experimental near-field images Re(Es,p) of a convex and a concave antenna extremity, respectively (λ0 = 11.06 μm). (g) Experimental near-field image Re(Es,p) of a bilayer graphene prism on a monolayer graphene (λ0 = 10.20 μm). (a) – (c) used with permission from20) copyright 2014 by American Physical Society, (d) – (g) used with permission from87) copyright 2014 by AAAS.

Fig. 4

(a) Electric field (Ex, perpendicular to graphene surface) distribution of SPPs with bias on the output arm by a gate voltage: (top) no bias as splitter, upper arm (middle) or lower arm (bottom) biased for digital switch. (b) Schematic of a monolayer graphene sheet array (MGSA). (c) Intensity distributions of SPPs in the MGSA for different initial phase differences (ϕ0) of SPPs between neighboring graphene. (a) used with permission from88) copyright 2012 by AIP Publishing, (b) and (c) used with permission from89) copyright 2012 by American Physical Society.

Fig. 5

(a) AFM and (b) SEM images of a graphene micro-ribbon array and a nano-ribbon array. (c) Control of terahertz resonance of plasmon excitations through electrical gating. Terahertz radiation was polarized perpendicular to the graphene ribbons. The inset shows extinction spectra for terahertz radiation polarized parallel to the ribbons. (d) Plasmon resonance frequency ωp as a function of |EF| or |n|1/2 for micro-ribbon arrays of different widths. (e) Normalized plasmon resonance frequency ωpw1/2 as a function of |EF| or |n|1/2 for micro-ribbon arrays of different widths. (f) SEM image (false color) of a stacked graphene/insulator micro-disk array arranged in a triangular lattice. (g) Extinction spectra for stacked plasmonic devices with one, two and five graphene layers under normal incidence. (h) Normalized plasmonic resonance frequency versus normalized total carrier density for single-layer graphene and graphene/insulator graphene stack. (i) Extinction spectra in single-layer graphene plasmonic devices with identical d (4.4 μm) but different a (4.8 and 9 μm, respectively). Inset: schematic of the graphene disk-disk interaction. (a) and (c) – (e) used with permission from32) copyright 2011 by Nature Publishing, (b) used with permission from34) copyright 2015 by AAAS, (f) – (i) used with permission from37) copyright 2012 by Nature Publishing.

Fig. 6

(a) Schematic of the structured graphene sheet at the flat surface between two dielectrics. (b) Angular dependence of the peak absorbance for s- and p-polarization in arrays of graphene nano-disks. (c) Normal-incidence absorbance by graphene disk arrays supported on a dielectric-coated gold surface (inset) for various values of the dielectric film thickness. (d) Incidence-angle and photon-energy dependence of the absorbance under the conditions of (a) for d = 1.12 μm. (e) The schematic of a tunable metasurface absorber composed of an aluminum film, aluminum oxide layer, and a tunable metasurface on graphene. (f) Measured reflection spectra from the metasurface absorber for different gate voltages. (a) – (d) used with permission from94) copyright 2012 by American Physical Society. (e) and (f) used with permission from65) copyright 2014 by American Chemical Society.

Fig. 7

(a) Graphene nano-crosses exhibit a giant optical activity (top & bottom) or a quarter wave plate (middle) at different wavelengths. (b, c) Anomalous refraction and flat lens composed of identical graphene bi-nanoribbon unit cells with different Fermi levels. (d) Anomalous reflection using graphene micro-ribbons/dielectric/metal with identical micro-ribbons but different Fermi levels and (e) An Airy beam generated by a graphene metasurface consisting of micro-ribbons/dielectric/metal with identical Fermi levels but different widths of micro-ribbons. (a) used with permission from40) copyright 2013 by The Optical Society, (b) and (c) used with permission from41) copyright 2014 by John Wiley and Sons, (d) used with permission from42) copyright 2015 by AIP Publishing, (e) used with permission from43) copyright 2015 by Nature Publishing.

Fig. 8

(a) Hybrid graphene/metasurface device by transferring a single layer graphene on top of a gold slit device. The incoming THz wave is polarized perpendicular to the slit orientation. (b) Electrical field enhancement factor in the slit for different surface conductivity of graphene. (c) Simulated (solid lines) and experimental (empty cycles) results of normalized THz power transmission modulation for bare graphene and the hybrid graphene-metasurface structure as a function of graphene conductivity at 1 THz. (d) SEM image of the fabricated Fano-resonant metasurface composed of asymmetric π-shape plasmonic resonators. (e) Simulated reflectivity spectra of hybrid graphene/metasurface for different Fermi levels. (f) Experimentally measured reflectivity spectra for the two samples with different geometrical dimensions before (solid curves) and after (dashed curves) graphene transfer. (g) A typical high resolution SEM image of the single gold nano-rod covered by graphene. (h) Scattering intensity (2D plot) is plotted as a function of the photon energy and gate voltage. (i) Four line cuts of the 2D plot for Rayleigh scattering spectra of an individual graphene-nano-rod hybrid structure. (a) – (i) used with permission from103,109,74) copyright 2015, 2013, 2012 by American Chemical Society.

Fig. 9

(a, b) Transmission spectra from the hybrid graphene/dielectric metasurface consisted of II-shape dielectric bars with structure asymmetry or dielectric blocks with periodity asymmetry for different Fermi levels. (c) SEM image of the hybrid graphene/dielectric waveguide modulator. (d) Static electro-optical response of the device for different gate voltages at 1.53 μm. (e) SEM image of the hybrid graphene/photonic crystal device. (f) Measured reflectivity from the hybrid device for different gate voltages. (a) used with permission from117) copyright 2015 by The Optical Society, (c) and (d) used with permission from16) copyright 2011 by Nature Publishing, (e) and (f) used with permission from123) copyright 2013 by American Chemical Society.