Translator Disclaimer
13 August 2020 Nanopolaritonic second-order topological insulator based on graphene plasmons
Author Affiliations +

Ultrastrong confinement, long lifetime, and gate-tunability of graphene plasmon polaritons (GPPs) motivate wide-ranging efforts to develop GPP-based active nanophotonic platforms. Incorporation of topological phenomena into such platforms will ensure their robustness as well as enrich their capabilities as promising test beds of strong light–matter interactions. A recently reported approach suggests an experimentally viable platform for topological graphene plasmonics by introducing nanopatterned gates—metagates. We propose a metagate-tuned GPP platform emulating a second-order topological crystalline insulator. The metagate imprints its crystalline symmetry onto graphene by modulating its chemical potential via field-effect gating. Depending on the gate geometry and applied voltage, the resulting two-dimensional crystal supports either one-dimensional chiral edge states or zero-dimensional midgap corner states. The proposed approach to achieving the hierarchy of nontrivial topological invariants at all dimensions lower than the dimension of the host material paves the way to utilizing higher-order topological effects for on-chip and ultracompact nanophotonic waveguides and cavities.



Nanoscale plasmonic excitations in graphene1 have proven their versatile potential for nanophotonic device applications2,3 and fundamental studies of strong light–matter interaction4,5 owing to their field-effect tunability and ultrasubwavelength confinement.6,7 In high-quality graphene encapsulated between hexagonal boron nitride (hBN) layers, these plasmon polaritons can propagate over long distances with minimal attenuation.810 In the past two years, researchers have sought new ways of combining the established concepts of topological photonics1117 with this attractive nanophotonic platform in order to further enhance its robustness. The earliest proposals for topological graphene plasmonics were based on nanopatterned graphene with broken time-reversal symmetry (TRS) via external DC magnetic fields.18,19 However, DC magnets are not compatible with overall device miniaturization and even present a hindrance to the plasmon detection tools such as scanning near-field microscopes.1,310 Moreover, graphene patterning, even if done over hBN-encapsulated structure for maximal protection,20 would considerably degrade the carrier mobility and introduce severe damping of the graphene plasmon polaritons (GPPs). Therefore, it is imperative to find an approach that preserves both the graphene quality and the TRS.

Although TRS-preserving photonic topological insulator (PTI) phases resembling spin-Hall electronic insulators have been realized by emulating spin-1/2 Kramers pairs using right/left circulating21 or transverse electric/magnetic22 modes, this approach is not applicable here because GPPs are essentially scalar-like excitations due to their quasielectrostatic nature. The remaining two options for TRS-invariant topological phases are then valley-Hall insulators2325 and topological crystalline insulators.2628 Earlier proposals29,30 for realizing these topological phases with GPPs suggested specific spatial landscapes of graphene optical conductivity, and our recent work31 suggested a practical path to realizing a GPP-based valley-Hall insulator using metagate-tuned graphene. A metagate, a designer nanopatterned planar electrode, imprints its structural symmetry onto the static carrier density landscape in graphene through the field-effect gating. Modulation of the carrier density translates into modulation of the optical conductivity, which creates a photonic crystal structure for GPPs. A similar method, which imparts a pattern to the dielectric spacer between graphene and backgate, has been experimentally shown to create a photonic crystal structure for GPPs while maintaining their long propagation distance.32 These approaches are not deleterious to the graphene quality because all patterning is done on the dielectric (or metallic) substrate while leaving the graphene itself intact.

In this paper, we present a metagate-graphene geometry that emulates topological crystalline insulators with GPPs. The employed photonic structure corresponds to a honeycomb lattice with a Kekulé distortion,33 which is known to produce two topologically distinct phases depending on the distortion directions.26,27 Orbital-locked chiral states arise at the boundary between two distinct domains, and these confined edge modes are shown to propagate freely around sharp corners with no reflection. Remarkably, a slightly modified structure can be used to emulate topological midgap corner states34,35 of second-order topological insulators (SOTIs), d-dimensional crystals supporting topological states of a dimension less than d1.3643 We show that our metagate geometry successfully realizes these phenomena by providing the first-principle calculations for the GPPs and qualitatively comparing the results to a generic tight binding model for Kekulé-textured graphene-like structures. While second-order topological effects have been demonstrated in macroscopic (photonic and phononic)3942 and waveguide-based34,43 systems, where the lattice spacing is much greater than the photon/phonon wavelengths, our work presents the first implementation of subwavelength polaritonic SOTI phases hosted by a two-dimensional (2-D) material. We also estimate how plasmonic band structures are affected by quantum nonlocal effects of Dirac electrons.4



Chemical potential (Fermi level) landscapes from field-effect gating using a metagate were obtained by electrostatic simulations in COMSOL. In order to incorporate the quantum capacitance effect of graphene, we considered the electrochemical equilibrium condition, which was done by setting a graphene sheet as a surface charge density element.31 The static surface charge density on graphene was given as a function of the local potential [n=(eV/vF)2/π]; therefore, an iterative solver was used. Photonic bands of the GPPs were calculated using both the in-house (Fourier-based quasielectrostatic) code31 and the finite-elements (COMSOL) simulations. The two codes agree well with each other, and we used them according to their computational efficiency: the in-house code for obtaining the results plotted in Figs. 1 and 4(c) and Fig. 6, COMSOL for the large-domain simulation shown in Figs. 4(d) and 5.

Fig. 1

(a) Schematic of the platform: graphene over a metagate. Left: lateral view of the entire structure and right: unit cell geometry of the metagate. (b) Metagate-induced Fermi energy EF(x,y) profile inside graphene (degenerate case: r1=r2=120  nm). Solid (dashed) line: unit cell for nondegenerate (degenerate) cases. (c) BZ and its high-symmetry points. Solid (dashed) hexagons: BZ boundaries for the nondegenerate (degenerate) cases. (d) Unfolded and (e) artificially folded BZs. Parameters: a=600  nm, h1=15  nm, h2=h3=10  nm, and V0=1  V throughout this paper.





Metagate-Induced Plasmonic Crystal

The dispersion relation q(ω) of the GPPs propagating in homogeneous graphene with a uniform chemical potential EF surrounded by an embedding medium with a dielectric constant ε is given by

Eq. (1)

where the semiclassical Drude conductivity,44 σDrude=i(e2/π2)|EF|/(ω+iγ), of graphene is assumed. The finite temperature effect would be negligible as long as kBT|EF|. Here q and ω are the GPP’s 2-D wave vector and frequency, respectively, γ is the phenomenological Drude loss corresponding to free-carrier scattering, e is the electron charge, and is Planck’s constant. Graphene’s chemical potential (corresponding to the Fermi level of its free carriers) can be varied within ±1  eV by tuning the carrier density n through field-effect gating: EF=sgn(n)vFπ|n|.45 Throughout this paper, we use the n-dependent renormalized expression for the Fermi velocity vF.46

Equation (1) implies that a periodic modulation in EF generates a periodic texture in the effective refractive index for the propagating GPPs on graphene, thereby creating a photonic crystal. A wide range of periodic modulation of EF can be achieved using a metagate;31 see Fig. 1(a) for the metagate geometry employed in the rest of this work. The metagate is a thin slab of metal (perfect electric conductor or PEC) perforated by a periodic array of unit cells containing a circular hole of radius r1 at its center and two identical holes of radius r2 at the periphery. When graphene is capacitively doped by this structured gate, the area above the empty holes attracts a smaller carrier density than the area above the unperforated regions. A periodic landscape of EF produced by metagate-induced doping is shown in Fig. 1(b) (a specific case of r1=r2 is depicted). With an eye on future experimental realizations, our electrostatic simulation assumes that graphene is hBN-encapsulated and accounts for graphene’s quantum capacitance self-consistently under electrochemical equilibrium in graphene.31,46

For the r1=r2 case shown in Fig. 1(b), the metagate is patterned as a honeycomb lattice with a smaller unit cell consisting of a single hole, as shown by the dashed line in Fig. 1(b). The photonic band structure (PBS) of the terahertz (THz)-frequency GPPs is shown in Fig. 1(d); it is calculated along the irreducible boundary of the Brillouin zone (BZ) defined from the smaller unit cell [dashed line in Fig. 1(c)]. As expected, the honeycomb photonic lattice produces a twofold degenerate conical dispersion at the two equivalent (k and k) valleys. The presence of a metallic gate in close proximity of the graphene layer is captured by the linear dispersion relation of the so-called acoustic GPPs9 at small momenta |q|a1 (i.e., near the γ-point of the BZ). Also, it is previously reported that the screening through this periodically structured gate, along with EF landscape, contributes to stronger modulation of PBS.31 We also note the subwavelength nature of the metagate (and, correspondingly, of the imprinted photonic crystal): the unit cell periodicity a=600  nm is much smaller than the free-space wavelength of the THz radiation. For example, at the k(k) point, we obtain an ultrahigh confinement ratio λvacuum/λGPP140.

The folded band structure in Fig. 1(e) is identical to the one in Fig. 1(d), except that it is drawn along the smaller BZ defined by the original unit cell for generic (r1r2) cases. The double degeneracy of the Dirac cones at Γ originates from the trivial band folding from the “true” (minimal) first BZ to an artificially shrunken BZ [i.e., from the dashed hexagon to the solid hexagon in Fig. 1(c)]. The four-fold degenerate modes at Γ shown here are the fundamental building blocks26,27 for the topological crystalline phases discussed below.


Kekulé-Textured Honeycomb Lattice

To highlight the underlying physics, we provide a well-known tight-binding (TB) model for a Kekulé-textured graphene-like system.26 Figure 2(a) illustrates the setting of the TB model and its mapping onto the metagate-tuned graphene structure shown in Fig. 1. The hexagonal unit cell of our metagate can be divided into six equivalent segments [see the shaded triangle in Fig. 2(a)] supported by C6 point group symmetry. If r1>r2 (r1<r2), the distance between the center-of-mass positions of the two adjacent metagate segments is shorter (greater) when they are positioned across the unit cell boundary versus within the same unit cell. Qualitatively, if we consider each segment to be an effective TB site, then the effective hopping strength within a unit cell tin is weaker (stronger) than that across adjacent unit cells tout in the case of r1>r2 (r1<r2).

Fig. 2

(a) Schematic of the TB model on a Kekulé lattice. Shaded triangular segment of a metagate maps onto a lattice site of the TB model. (b) PBS calculated for tin=1.1t/tout=0.9t (left) and tin=0.9t/tout=1.1t (right). Color: ratio between the projection amplitudes onto the dipole-like orbitals |p1,2 and onto the quadrupole-like orbitals |d1,2. (c) PBS calculated for a structure interfacing the two configurations from (b) along the zigzag-oriented domain wall. Bandgap-crossing purple lines: dispersion curves for the edge states propagating along the domain wall.


Any eigenstates on six hexagonal vertices can be described by linear combinations of six orbital states: a monopole-like (|s=[1;1;1;1;1;1]), two dipole-like (|p1=[1;1;0;1;1;0], |p2=[1;1;2;1;1;2]), two quadrupole-like (|d1=[1;1;0;1;1;0], |d2=[1;1;2;1;1;2]), and an octopole-like (|f=[1;1;1;1;1;1]).26 The monopole and octopole modes are at high energies, playing negligible roles for low-energy physics (high/low compared to tin and tout). When two hopping strengths are the same, the dipole and quadrupole modes are degenerate at zero Bloch momentum just like in Fig. 1(e). This is because, if tin=tout, |p1,2 and |d1,2 become equivalent with respect to the redefinition of the unit hexagon (around the center hole or around a periphery hole). This fourfold degeneracy is lifted when tintout, as shown in Fig. 2(b). Specifically, the dipole (quadrupole) modes can be shown to reside below (above) the bandgap for tout<tin. Band inversion occurs between these four orbitals around the Γ-point for the tout>tin case, thus creating an insulating phase that is topologically distinct from the “topologically trivial” insulating phase realized in the case of tout<tin.

In electronic systems, such topologically nontrivial insulators are sometimes referred to as quantum orbital-Hall insulators26 because their underlying cause—band inversion of orbital-locked states—is reminiscent of that in quantum spin-Hall insulators. Indeed, the low-energy effective four-band Hamiltonian around the Γ point is characterized by a Z2 topological invariant16,26,27 (Z2=0 for tout<tin and Z2=1 for tout>tin) as in quantum spin-Hall insulators (see the Supplementary Material). Thus following the bulk-edge correspondence,1113 linear edge states should arise at the interface between the topological and trivial domains. Indeed, we observe two expected domain wall states in the bandgap, as shown in Fig. 2(c). The forward (backward) propagating mode is linked to the orbital states such that the relative phase of the lattice sites within a unit cell is increasing in the counter-clockwise (clockwise) direction.


Topological Crystalline Phase in Graphene Plasmons

Next, we confirm that our metagate-tuned GPPs indeed exhibit the two topologically distinct phases introduced above. This is done by comparing the band structures and examining the plasmonic mode profiles above and below the bandgap. Figure 3(a) provides the GPP PBS calculated for two different cases: r1=90  nm and r2=150  nm (topologically trivial case), and the topologically nontrivial case of r1=160  nm and r2=110  nm. Here r1 and r2 are chosen such that any adjacent holes are clearly separated by at least a 40-nm gap. Plasmonic bandgaps of similar sizes are opening for both configurations, and the overall shape of the bands agrees with the band structures of the TB model shown in Fig. 2(b). The dynamic charge density in graphene associated with the propagating GPPs is plotted in Fig. 3(b) for a single unit cell.

Fig. 3

(a) Plasmonic band structures calculated for r1=90  nm/r2=150  nm (left: topologically trivial domain) and r1=160  nm/r2=110  nm (right: topologically nontrivial domain). Line color: dipole-to-quadrupole ratios (see text). (b) Color: AC surface charge density in graphene associated with the GPPs at the Γ-point (below the bandgaps). Parameters: same as in Fig. 1.


Next, we color-coded [Fig. 3(a)] each eigenstate according to the ratio between the dipole and quadrupole strengths. By comparing the left and right plots in Fig. 3(a), we note that the band inversion predicted by the TB model is indeed observed for the GPPs propagating through the modulated two-dimensional Fermi energy landscape EF(x,y) in graphene. As depicted in Fig. 3(b), the plasmonic modes near the Γ point display dipole (quadrupole) charge distributions below the bandgap when r1 is smaller (greater) than r2. Although the visualized dipole-quadrupole ratio does not show a perfect agreement between the plasmonic band structure and the TB model, this is due to the fact that the actual plasmonic modes are not strongly bounded as in the idealized TB model. To account for weak localization more precisely, the effective TB description would need to include the second (or higher) nearest neighbor couplings. These minor discrepancies, however, can be neglected since the simplest TB model can capture all the topological features of our system.


Topologically Protected Plasmonic Edge States

To observe the topologically protected edge states propagating along the domain wall between the trivial and nontrivial topological domain, it is necessary to carry out additional engineering of the interface between the two domains. This is needed because the microscopic field profiles of the dipole and quadrupole modes are not identical inside the two domains owing to the change of the unit cell structure. Therefore, when the two distinct metagate domains are interfaced, these abrupt changes across the domain wall can introduce a small bandgap between two edge states.26,27 To remove this gapping of the edge states, we engineer the domain wall by setting the radius r0 of the peripheral metagate holes (i.e., the holes at the domain wall) to be the average of the two r2 values from each domain: r¯=(r2e+r2s)/2. Here the subscripts s and e stand for the shrunken (topologically trivial) and expanded (topologically nontrivial) configuration, respectively. This additional interface engineering mitigates the local C6 symmetry breaking at the interface.

Next, we calculate the PBS of the GPPs propagating along the zigzag-shaped domain wall oriented along the x direction. The computational domain is comprised of 18 unit cells on each side of the domain wall. The results are shown in Fig. 4(a), where the black thin lines correspond to the projected bulk modes (i.e., having the same value of kx), whereas the red (blue) solid lines correspond to the forward (backward) propagating chiral edge states. The edge states have been found to be exponentially localized at the domain wall, as shown in Figs. 4(c) and 4(d), while the bulk modes are delocalized (not shown). The effectiveness of the interface engineering described above is apparent from the fact that the gap between the forward- and backward-propagating edge modes is extremely small, as can be seen in Fig. 4(a). The existence of an extremely small but finite gapping of the edge modes due to symmetry breaking at the domain wall has been noted earlier.27

Fig. 4

(a) Projection of the PBS of the GPPs supported by a graphene nanoribbon gated by a metagate comprised of the two domains shown in Fig. 3. Each domain contains 18 unit cells, zigzag-type domain wall runs parallel to the x-direction, the radius of peripheral holes at the domain wall: r¯=130  nm. Black thin lines: x-projected bulk bands; red (blue) lines and circles: forward (backward)-propagating chiral edge states. (b) Metagate-induced Fermi energy profile on graphene near the domain boundary (white dotted line). (c) AC plasma charge density profiles for the forward- and backward-propagating edge states. 2-D color code: the magnitude and phase of the charge density. The gray arrows: guiding to the eye for the directions of the phase increase. (d) Reflection-free propagation of topological edge modes along the domain wall with sharp corners. The forward (backward) mode on the left (right) is excited at f=4.05  THz with a surface current density (green arrows) J+=x^+iy^ (J=x^iy^). Color plots are absolute squares of the out-of-plane electric field. Parameters: same as in Figs. 1 and 3.


As predicted by the TB model, the forward-propagating mode (red circle) has the AC surface charge density distribution with the phase increasing around the center holes in the counterclockwise direction, see Fig. 4(c). Due to the mirror symmetry in the nanoribbon, the backward-propagating mode (blue circle) is a mirror image of the forward-propagating mode. Therefore, these two modes form a chiral pair. We next investigated the topological robustness of these chiral edge states by modeling their propagation along a complex pathway with several types of sharp corners along the domain wall. For this simulation, each mode was selectively excited by circularly polarized external surface currents [green lines in Fig. 4(d): J±=x^±iy^]. The modeling results shown in Fig. 4(d) clearly indicate that the edge states are indeed immune to backscattering at all—60 deg, 90 deg, or 120 deg— lattice-preserving sharp turns.


Midgap Topological Corner States

In addition to the Z2 invariant introduced above for characterizing the quantum orbital-Hall effect, the C6 crystalline insulator carries another binary invariant associated with topological corner charges:35 qc0 (mod 1) for tout<tin and 1/2 for tout>tin (see the Supplementary Material). The nontrivial half corner charge, given that the chiral symmetry is respected both inside the domain and at its termination, gives rise to a corner-localized mode in the middle of the bandgap.34 This midgap corner state is a key-signature of second-order topological effects, where a band topology of a d-dimensional crystal is manifested in the (d2)-dimensional observable. For 2-D crystalline PTIs, such observables are corner states that can be viewed as topologically robust and highly localized photonic cavities. Below we demonstrate that such nanocavities for GPPs can be developed using the metagate-graphene platform that enables us to embed a topologically nontrivial PTI inside a topologically trivial one.

One challenge in designing corner states in 2-D plasmonic systems is that, unlike the case of a waveguide-based platform34,43 that is naturally described within the framework of a TB model, our plasmonic crystal does not possess strongly localized/bounded lattice sites. Therefore, any abrupt termination of a domain causes large shifts in effective on-site potentials, i.e., breaks the local chiral symmetry around a corner of the structure. To be specific, when graphene is set to be undoped or removed outside the domain boundary of our plasmonic crystal, the resulting structure is mapped to a TB model with severe redshifts in on-site potentials of boundary lattice sites (see the Supplementary Material). To avoid such unwanted perturbation along the boundary, we choose to imbed the topological domain within a trivial C6 crystalline insulator. In other words, while we are utilizing the same configuration as the one used for demonstrating bandgap-spanning chiral edge states in Fig. 4, we are now pursuing a different objective: to produce a sufficient bandgap between the two edge states in order to spectrally isolate the midgap state inside it. Therefore, we again used interface engineering to maximize the gapping of the edge states. As shown in Fig. 5(a), the size of the edge bandgap can be tuned from almost negligible [as was done in Figs. 4(a) and 4(b), where peripheral hole radii at the domain wall were chosen as r¯=(r2s+r2e)/2] to large. The latter objective is accomplished by staggering the radii of the peripheral holes between r2s and r2e, as shown on the right side of Fig. 5(a). Such interface engineering amounts to edge roughening and enables us to increase the edge bandgap in such a way that the corner states appear near the middle of the bandgap.

Fig. 5

Designer interfaces between topologically distinct PTIs enable nanoscale waveguides and cavities. (a) Left: a uniform domain wall configuration used in Fig. 4 for supporting bandgap-spanning chiral edge modes. Right: staggered domain wall for supporting cavity-like corner modes. Small icons above the arrow denote the 1-D band structures for each domain wall type. (b) Left: gapping of the edge states supported by the staggered domain wall (r1e=165  nm, r2e=120  nm, r1s=100  nm, and r2s=155  nm) enables corner states. Right: DOS for the topological (expanded) domain (10×10 unit cells) imbedded into a trivial (shrunken) domain. Inset: the field profile of a midgap corner state. (c) A TB model qualitatively mimics the imbedded structure with a staggered domain wall. Small circles: sublattice sites, color-coded links: hopping amplitudes. (d) Same as (b), but for the TB model shown in (c).


Figure 5(b) shows the emergence of the expected topological corner state. An eigensimulation for a structure with 120-deg-angled corners reveals an eigenstate per corner that is well inside the edge bandgap and localized around the corner. Also, the mode profile of this corner state [Fig. 5(b) inset] features dominant amplitudes at two effective lattice sites that are the second corner-most, which exactly matches the reported result in an earlier work.34 The density of states (DOS) was defined as Im[n1ωnω], where ωn is the complex eigenfrequency (energy).

The proposed approach of embedding an SOTI into a topologically trivial structure and appropriately engineering the interface between them can be applied to any subwavelength platform. For example, we show in Figs. 5(c) and 5(d) that a generic TB model predicts the same results. Figure 5(c) describes the way we incorporate the edge roughening effect into the model. The relative strengths of hopping amplitudes are determined by considering the widths of metagate junctions, at which two adjacent effective lattice sites [Fig. 2(a)] meet together. Comparing to the setting drawn on the right side of Fig. 5(a), the hopping strengths are assigned to be proportionate to the junction widths. Then the following results in Fig. 5(d) confirm that the midgap corner state is guaranteed in this scheme. The reason for this is that our plasmonic crystal is not strictly chiral symmetric, as it is observed in Fig. 3(a) that the lower two bands are not exact mirror images of the upper two bands with respect to the midgap frequency. The exact emulation of the chiral symmetric Hamiltonian is not possible due to the weak localization of effective lattice sites and the presence of additional bands at higher frequencies.

To ensure that the DOS peak corresponding to a corner state is clearly standing out from the DOS due to the edge states, intrinsic graphene losses must be taken into account. At cryogenic temperatures (T=20 to 100 K), the intrinsic plasmonic loss in graphene is primarily due to the pseudomagnetic field induced by electron–phonon scattering:10 γ/2π=[1.5×(EFineV)×(TinK)]  GHz (see the Supplementary Material). Thus, for the DOS calculation in Fig. 5(b) we used γ/2π=20  GHz from EF=0.15  eV and T=90  K. Then the peak from the corner state contribution is clearly separated from other peaks, indicating that the observation of the corner state would be experimentally feasible at cryogenic settings.


Nonlocal Optical Response of Dirac Electrons

So far, we have used the semiclassical Drude model [Eq. (2)] to describe graphene optical conductivity. Several experimental works3,4,8,9 have, however, revealed that the quantum nonlocal response of Dirac electrons needs to be considered for the precise prediction of graphene plasmon dispersion. The nonlocal optical conductivity in homogeneously doped graphene is described accurately enough at the random phase approximation (RPA) level24,31

Eq. (2)

where f(z) is a dimensionless correction factor that describes the nonlocal effect,

Eq. (3)


Equation (2) is valid only when vF|q|ω<1 and ω|EF|. Note that f(z) is an increasing function in z bounded below by 1. This momentum-dependent rise of conductivity stems from velocity matching between the Dirac electrons, traveling at vF, and plasmons, traveling at vp=ω/|q|; as z=(vF/vp)21. The electrons interact longer with each plasma wavefront and, therefore, exhibit an enhanced response.4

For nonlocal RPA conductivity under a spatially varying Fermi level, as in our metagate-tuned graphene, we adopt a semiphenomenological extension from the homogeneous case31,47

Eq. (4)

where f[n](z) is the Taylor expansion of f(z) in Eq. (3) up to the n’th order term, and E˜F(q) is the Fourier transform of EF(r). Equation (4) reduces back to Eq. (2) under a constant Fermi level: E˜F(q)=EFδq,0σpheno(ω,q,q)=δq,qσRPA(ω,q). Also, the local Drude model is equivalent to the case n=0. The effectiveness of Eq. (4) with n=1 has been verified in a recent experimental work.3

Figure 6 shows how much the topological bandgap in our metagate-tuned plasmonic crystal gets shifted by the quantum nonlocal effect. We consider two cases: (a) V0=1  V, as assumed so far, and (b) V0=0.05  V, to demonstrate a case where nonlocal responses become significant. With V0=10.05  V, GPP frequencies are roughly halved; ωV01/4 [Eq. (1)]. As stated earlier, Eqs. (2) and (4) are valid only when ρvF|q|ω<1. Here the dimensionless quantity ρ is a measure of the importance of the nonlocal quantum effects.31 Bloch eigenstates in a periodic system contain many wavevector components, but here we take only the primary wavevector for simplicity in estimating ρ.18 The primary wavevector for dipole and quadrupole states around the plasmonic bandgap shown in Fig. 3(a) is Γ with the first band folding; that is, its magnitude is |q|=4π3a. For case (a), ω/2π is around 4 THz, and thereby we get ρ0.5. As a result, the resonance shift shown in Fig. 6(a) is insignificant compared to the bandgap size, and the nonlocal effect converges rapidly at the first two orders. For case (b), we get ρ1.0 from ω/2π2  THz. With ρ approaching near unity, the nonlocal correction already exceeds the bandgap size at the first order and keeps increasing as the higher-order terms are added [see Fig. 6(b)]. Thus Eq. (4) is applicable only under high gate voltages and fails to be an effective description as ρ1.

Fig. 6

Shifts of plasmonic resonance frequencies due to nonlocal responses; eigenfrequencies are evaluated at Γ above and below topological bandgap, using σpheno[n] in Eq. (4) for n=0 (local Drude limit) and n=1,2 (nonlocal corrections up to the first or the second order). (a) With V0=1  V, the bandgap is located around 4 THz and its size is about 0.3 THz. (b) With V0=0.05  V, the bandgap is located around 2 THz and its size is about 0.15 THz.




We emphasize that our metagate-graphene structure was designed for experimentally feasible low-loss propagation of topological GPP excitations. Thus our structure is designed to operate at a few THz frequency range, where the dielectric losses from hBN or SiO2 phonon absorptions are vanishing. Then the 1/e decay length of the topological chiral GPPs in Fig. 4 can be estimated to be ωγ1|q|30a if we assume γ/2π=20  GHz. Considering that the lateral confinement of the edge states shown in Fig. 4(c) is 3 to 4 unit cells, it should be possible to observe the propagation along the sharp turns shown in Fig. 4(d). Also, we expect that the second-order topological corner state can be efficiently excited, featuring a high quality factor (Q400). Again, our metagate scheme plays an essential role in enabling this low-loss operation through avoiding direct patterning of graphene.

Our metagate-tuned graphene system serves as an interesting case where topological photonics becomes a useful tool for unraveling condensed matter properties. Figure 6 suggests that we can electrically tune the proposed system from a condition in which Eq. (4) is effective, to another in which a modification is required. In this sense, our platform can be used to study quantum nonlocal effects in periodically doped graphene. For example, it is known that when the nonlocal response strength parameter ρ becomes greater than 1, graphene plasmons experience severe Landau damping.4,18 The Bloch states in our plasmonic crystal are described by multiple ρ values from multiple wavevector components. It is an intriguing problem to ask whether the reflection-free propagation of 1-D topological Bloch states would be disabled by this selective Landau damping acting on a particular set of wavevectors. Thus exotic optical properties of Dirac electrons under strongly modulated Fermi surfaces can be revealed by examining how our plasmonic crystal deviates from the predictions based on local descriptions or trivial nonlocal corrections.

In conclusion, we have designed a nanophotonic platform for achieving two types of topologically robust optical elements based on graphene plasmons: a topological crystalline insulator for making reflections-free waveguides and a second-order PTI for making robust nanocavities. The platform consists of a single-layer graphene that is electrically biased using a designer metagate: an array of holes forming a honeycomb lattice with a Kekulé distortion. We demonstrate that the change from the propagating GPPs guided in two dimensions to cavity-like plasmons localized in all three dimensions can be accomplished by engineering the interface between two plasmonic PTIs: trivial and topological. We envision that the metagate geometry introduced in this work can also be used for realizing fractionalized Kekulé-vortex states33 with graphene plasmons as well as many other topological phenomena such as strain-induced pseudomagnetic fields.48 Also, by adding the second gate (a top metagate or a simple backgate beneath the existing metagate), the electric control over topological plasmons can be greatly enhanced to enable switching between distinct low-dimensional topological phases.49 This would constitute a major advance over the single-gate geometry analyzed here, e.g., switching between two domain wall configurations in Fig. 5(a). Our metagate approach, combined with engineered domain walls between topologically different PTIs, paves the way to novel methods of manipulating nanolight in graphene and for developing topologically robust devices.


This work was supported by the U.S. Army Research Office (ARO) under Grant No. W911NF-16-1-0319 and the National Science Foundation (NSF) under Grant Nos. DMR-1741788 and DMR-1719875. M. J. also acknowledges the support from the Kwanjeong Fellowship from the Kwanjeong Educational Foundation.



S. S. Sunku et al., “Photonic crystals for nano-light in moiré graphene superlattices,” Science, 362 (6419), 1153 –1156 (2018). SCIEAS 0036-8075 Google Scholar


X. Lin et al., “All-angle negative refraction of highly squeezed plasmon and phonon polaritons in graphene–boron nitride heterostructures,” Proc. Natl. Acad. Sci. U. S. A., 114 (26), 6717 –6721 (2017). Google Scholar


A. Woessner et al., “Electrical 2π phase control of infrared light in a 350-nm footprint using graphene plasmons,” Nat. Photonics, 11 421 –424 (2017). NPAHBY 1749-4885 Google Scholar


M. B. Lundeberg et al., “Tuning quantum nonlocal effects in graphene plasmonics,” Science, 357 (6347), 187 –191 (2017). SCIEAS 0036-8075 Google Scholar


D. A. Iranzo et al., “Probing the ultimate plasmon confinement limits with a van der Waals heterostructure,” Science, 360 (6386), 291 –295 (2018). SCIEAS 0036-8075 Google Scholar


J. Chen et al., “Optical nano-imaging of gate-tunable graphene plasmons,” Nature, 487 77 –81 (2012). NATUAS 0028-0836 Google Scholar


Z. Fei et al., “Gate-tuning of graphene plasmons revealed by infrared nano-imaging,” Nature, 487 82 –85 (2012). NATUAS 0028-0836 Google Scholar


A. Woessner et al., “Highly confined low-loss plasmons in graphene–boron nitride heterostructures,” Nat. Mater., 14 421 –425 (2015). NMAACR 1476-1122 Google Scholar


P. Alonso-Gonzalez et al., “Acoustic terahertz graphene plasmons revealed by photocurrent nanoscopy,” Nat. Nanotechnol., 12 31 –35 (2017). NNAABX 1748-3387 Google Scholar


G. X. Ni et al., “Fundamental limits to graphene plasmonics,” Nature, 557 530 –533 (2018). NATUAS 0028-0836 Google Scholar


A. B. Khanikaev et al., “Photonic topological insulators,” Nat. Mater., 12 233 –239 (2013). NMAACR 1476-1122 Google Scholar


L. Lu, J. D. Joannopoulos and M. Soljačić, “Topological photonics,” Nat. Photonics, 8 821 –829 (2014). NPAHBY 1749-4885 Google Scholar


C. He et al., “Photonic topological insulator with broken time-reversal symmetry,” Proc. Natl. Acad. Sci. U. S. A., 113 (18), 4924 –4928 (2016). Google Scholar


K. Lai et al., “Experimental realization of a reflections-free compact delay line based on a photonic topological insulator,” Sci. Rep., 6 28453 (2016). SRCEC3 2045-2322 Google Scholar


A. B. Khanikaev and G. Shvets, “Two-dimensional topological photonics,” Nat. Photonics, 11 763 –773 (2017). NPAHBY 1749-4885 Google Scholar


M. A. Gorlach et al., “Far-field probing of leaky topological states in all-dielectric metasurfaces,” Nat. Commun., 9 909 (2018). NCAOBW 2041-1723 Google Scholar


F. Gao et al., “Topologically protected refraction of robust kink states in valley photonic crystals,” Nat. Phys., 14 140 –144 (2018). NPAHAX 1745-2473 Google Scholar


D. Jin et al., “Infrared topological plasmons in graphene,” Phys. Rev. Lett., 118 (24), 245301 (2017). PRLTAO 0031-9007 Google Scholar


D. Pan et al., “Topologically protected Dirac plasmons in a graphene superlattice,” Nat. Commun., 8 1243 (2017). NCAOBW 2041-1723 Google Scholar


B. S. Jessen et al., “Lithographic band structure engineering of graphene,” Nat. Nanotechnol., 14 340 –346 (2019). NNAABX 1748-3387 Google Scholar


M. Hafezi et al., “Robust optical delay lines with topological protection,” Nat. Phys., 7 907 –912 (2011). NPAHAX 1745-2473 Google Scholar


T. Ma et al., “Guiding electromagnetic waves around sharp corners: topologically protected photonic transport in metawaveguides,” Phys. Rev. Lett., 114 (12), 127401 (2015). PRLTAO 0031-9007 Google Scholar


L. Ju et al., “Topological valley transport at bilayer graphene domain walls,” Nature, 520 650 –655 (2015). Google Scholar


T. Ma and G. Shvets, “All-Si valley-hall photonic topological insulator,” New J. Phys., 18 025012 (2016). NJOPFM 1367-2630 Google Scholar


J. Lu et al., “Observation of topological valley transport of sound in sonic crystals,” Nat. Phys., 13 369 –374 (2017). NPAHAX 1745-2473 Google Scholar


L.-H. Wu and X. Hu, “Topological properties of electrons in honeycomb lattice with detuned hopping energy,” Sci. Rep., 6 24347 (2016). SRCEC3 2045-2322 Google Scholar


L.-H. Wu and X. Hu, “Scheme for achieving a topological photonic crystal by using dielectric material,” Phys. Rev. Lett., 114 (22), 223901 (2015). PRLTAO 0031-9007 Google Scholar


R. Süsstrunk and S. D. Huber, “Classification of topological phonons in linear mechanical metamaterials,” Proc. Natl. Acad. Sci. U. S. A., 113 (33), E4767 –E4775 (2016). Google Scholar


P. Qiu et al., “Topologically protected edge states in graphene plasmonic crystals,” Opt. Express, 25 (19), 22587 –22594 (2017). OPEXFF 1094-4087 Google Scholar


P. Qiu et al., “Pseudospin dependent one-way transmission in graphene-based topological plasmonic crystals,” Nano. Res. Lett., 13 113 (2018). NRLAAD 1556-276X Google Scholar


M. Jung, Z. Fan and G. Shvets, “Midinfrared plasmonic valleytronics in metagate-tuned graphene,” Phys. Rev. Lett., 121 (8), 086807 (2018). PRLTAO 0031-9007 Google Scholar


L. Xiong et al., “Photonic crystal for graphene plasmons,” Nat. Commun., 10 4780 (2019). NCAOBW 2041-1723 Google Scholar


C.-Y. Hou, C. Chamon and C. Mudry, “Electron fractionalization in two-dimensional graphenelike structures,” Phys. Rev. Lett., 98 (18), 186809 (2007). PRLTAO 0031-9007 Google Scholar


J. Noh et al., “Topological protection of photonic mid-gap defect modes,” Nat. Photonics, 12 408 –415 (2018). NPAHBY 1749-4885 Google Scholar


W. A. Benalcazar, T. Li and T. L. Hughes, “Quantization of fractional corner charge in cn-symmetric higher-order topological crystalline insulators,” Phys. Rev. B, 99 (24), 245151 (2019). Google Scholar


W. A. Benalcazar, B. A. Bernevig and T. L. Hughes, “Quantized electric multipole insulators,” Science, 357 (6346), 61 –66 (2017). SCIEAS 0036-8075 Google Scholar


W. A. Benalcazar, B. A. Bernevig and T. L. Hughes, “Electric multipole moments, topological multipole moment pumping, and chiral hinge states in crystalline insulators,” Phys. Rev. B, 96 (24), 245115 (2017). Google Scholar


C. W. Peterson et al., “A quantized microwave quadrupole insulator with topologically protected corner states,” Nature, 555 346 –350 (2018). NATUAS 0028-0836 Google Scholar


M. Serra-Garcia et al., “Observation of a phononic quadrupole topological insulator,” Nature, 555 342 –345 (2018). NATUAS 0028-0836 Google Scholar


F. Schindler et al., “Higher-order topological insulators,” Sci. Adv., 4 (6), eaat0346 (2018). STAMCV 1468-6996 Google Scholar


H. Xue et al., “Acoustic higher-order topological insulator on a kagome lattice,” Nat. Mater., 18 108 –112 (2019). NMAACR 1476-1122 Google Scholar


X. Ni et al., “Observation of higher-order topological acoustic states protected by generalized chiral symmetry,” Nat. Mater., 18 113 –120 (2019). NMAACR 1476-1122 Google Scholar


S. Mittal et al., “Photonic quadrupole topological phases,” Nat. Photonics, 13 692 –696 (2019). NPAHBY 1749-4885 Google Scholar


M. Jablan, H. Buljan and M. Soljacic, “Plasmonics in graphene at infrared frequencies,” Phys. Rev. B, 80 (24), 245435 (2009). Google Scholar


D. K. Efetov and P. Kim, “Controlling electron-phonon interactions in graphene at ultrahigh carrier densities,” Phys. Rev. Lett., 105 (25), 256805 (2010). PRLTAO 0031-9007 Google Scholar


R. J. G. L. Yu et al., “Interaction phenomena in graphene seen through quantum capacitance,” Proc. Natl. Acad. Sci. U. S. A., 110 (9), 3282 –3286 (2013). Google Scholar


I. Torre et al., “Lippmann–Schwinger theory for two-dimensional plasmon scattering,” Phys. Rev. B, 96 (3), 035433 (2017). Google Scholar


C. Brendel et al., “Pseudomagnetic fields for sound at the nanoscale,” Proc. Natl. Acad. Sci. U. S. A., 114 (17), E3390 –E3395 (2017). Google Scholar


Z. Fan et al., “Electrically defined topological interface states of graphene surface plasmons based on a gate-tunable quantum Bragg grating,” Nanophotonics, 8 (8), 1417 –1431 (2019). Google Scholar


Minwoo Jung received his master’s degree in physics from Cornell University in 2020. He is currently a PhD student in the Department of Physics at Cornell University. His research topic mainly focuses on utilizing nanopolaritons in two-dimensional materials as platforms for advanced photonics.

Ran G. Gladstone received his master’s degree in applied physics from Cornell University in 2019. He is currently a PhD student at the School of Applied Physics at Cornell University. His main research focus is implementing theoretical topological insulator models in realistic photonic systems.

Gennady Shvets received his PhD in physics from MIT in 1995. He is a professor of applied and engineering physics at Cornell University. His research interests include nanophotonics, optical and microwave metamaterials and their applications (including biosensing, optoelectronic devices, and vacuum electronics), topological concepts in photonics, and ultraintense laser–matter interactions. He is the author or co-author of more than 200 papers in refereed journals. He is a fellow of the American Physical Society, Optical Society of America, and SPIE.

© The Authors. Published by SPIE and CLP under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Minwoo Jung, Ran G. Gladstone, and Gennady B. Shvets "Nanopolaritonic second-order topological insulator based on graphene plasmons," Advanced Photonics 2(4), 046003 (13 August 2020).
Received: 28 May 2020; Accepted: 22 July 2020; Published: 13 August 2020

Back to Top