Abstract
We propose using a topological plasmonic crystal structure composed of an array of nearly parallel nanowires with unequal spacing for manipulating light. In the paraxial approximation, the Helmholtz equation that describes the propagation of light along the nanowires maps onto the Schrödinger equation of the Su-Schrieffer-Heeger (SSH) model. Using a full three-dimensional finite difference time domain solution of the Maxwell equations, we verify the existence of topological defect modes, with sub-wavelength localization, bound to domain walls of the plasmonic crystal. We show that by manipulating domain walls we can construct spatial mode filters that couple bulk modes to topological defect modes, and topological beam-splitters that couple two topological defect modes. Finally, we show that the structures are tolerant to fabrication errors with an inverse length-scale smaller than the topological band gap.
© 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
The interface between two electronic materials with topologically distinct band structures necessarily supports a topologically protected mode [1–3]. The robustness of this mode arises from global symmetries like parity, space inversion, and time reversal [4,5]. In electronic systems, the topology of the electronic band structures can be controlled relatively easily by either using materials properties (e.g. picking a suitable semiconductor [6]) or by materials engineering (e.g, by making hybrid structures like quantum wells out of suitable materials [7]). The combination of robustness of topological edge modes and the relative ease of fabrication has resulted in an explosion of interest in these systems with applications ranging from decoherence-free quantum state manipulation [8] to electronic device design [9–11].
Following the seminal work of [12], it was realized that topological band structures can also exist in photonic systems. These ideas have been explored in a number of theoretical proposals [13–19]. They have also been realized experimentally in the following photonic crystal systems: gyromagnetic photonic crystal at microwave frequencies [20], coupled whispering gallery mode optical resonators [21], optical waveguides [22,23], and plasmonic arrays [24–26].
The analysis of Refs. [22], which provides a mapping between electronic systems and paraxial light propagation in photonic and plasmonic systems, inspired a number of works. Here, we are concerned with the one-dimensional topological model of Su-Schrieffer-Heeger (SSH) [27,28], which was originally used to describe electron motion in polyacetylene chains. The SSH model hosts topologically protected modes on domain boundaries, the protection being provided by the sub-lattice (or chiral) symmetry which belongs to the AIII-class of [4]. The mapping inspired by Rechtsman et al’s work was adopted to the SSH model and explored experimentally in plasmonic nanowire arrays [24, 29] as well as photonic crystals [23]. The propagation of plasmonic topological modes in graphene nanowires was also explored theoretically in [30]. These works serve to establish the existence of optical topological defect modes.
In this work, we explore how to use SSH model plasmonic topological modes to manipulate light. We provide an in-depth analysis of the topological defect modes, similar to the ones experimentally observed in [23, 24], and propose applications of these modes in plasmonic systems. We show that the proposed setup allows us not only to guide light but also to robustly manipulate it by shifting the topological defect modes, as a function of axial position, inside the structure. Using full 3D finite-difference-time-domain (FDTD) solutions of the Maxwell equations, we demonstrate a beam splitter and a spatial mode filter (that couples light from a pair of bulk modes to a pair of topological defect modes). Further, we analyze tolerance to manufacturing defects. While topological robustness to disorder that is invariant in the axial direction (i.e. time-independent disorder in electronic systems) has been extensively studied before, here we analyze tolerance to disorder that varies in the axial direction (i.e. time-dependent disorder in electronic systems). We investigate two types of defects: (a) defects in nanowire positions, that violate sublattice symmetry exponentially weakly and (b) defects in nano-wire diameter, that violate sublattice symmetry strongly. Using the 3D FDTD calculations we show that our structures are very tolerant to type (a) defects and reasonably tolerant to type (b) defects with transverse length-scales comparable to inter-nanowire spacing, as long as the axial length-scales are larger than the inverse topological band gap.
2. Guiding light using topological defect modes
Consider a system that is almost translationally invariant along the axial, i.e. z-direction. Our goal is to describe paraxial modes, time-harmonic electromagnetic waves that propagate at small angles to the z-axis. For paraxial TM-modes it is natural to focus on the transverse components of the electric field [31,32] and to separate out the fast oscillating part, Ex = ψx(x, y, z)eiβ0z and Ey = ψy(x, y, z)eiβ0z, where β0 = ω/c. The propagation of paraxial TM-modes is governed by the paraxial Schrödinger equation (see appendix for details)
where, , V(x, y, z) describes the position of the metal nanowires, and we have neglected as per the paraxial approximation. Eq. (1) has the form of the two-dimensional time-dependent Schrödinger equation with c ∂z → ∂t, ψ{x,y} – the wave function, and – the Hamiltonian operator. Thus, stationary states of H become the TM modes of paraxial light that propagate along the z-direction and the eigenenergies ℰ of H become the z-wavenumbers βz = β0 − ∊/c. In [22], Rechtsman et al showed that this analogy can be used to gain intuition about the topological structure of electromagnetic waves by mapping solutions of the Schrödinger equation with non-trivial topology onto the Helmholtz equation.Let us now consider the SSH model, which was originally used for studying electrons in polyacetylene. The backbone of polyacetylene is a chain of carbon atoms with staggered single and double bonds, schematically depicted in Fig. 1(a). (Here, we assume that the pattern of single and double bonds does not have intrinsic dynamics, but is instead a prescribed function of time.) The discretized version of the Hamiltonian describing the hopping of spinless electrons along the backbone of polyacetylene is
where the operator ci () annihilates (creates) an electron on the i-th carbon atom and the hopping matrix element ti,i+1 = t1(t2) if the bond between sites i and i + 1 is a single (double) bond. The unit cell of this model consists of a two atom dimer. Let us assume the hopping strength within a dimer is t1 and the hopping strength between the atoms belonging to two nearby dimers is t2. Applying the Fourier transform to the Hamiltonian, i.e. and , where k is the wave vector, R is the unit cell size, and n is the unit cell index, we obtain This Hamiltonian can be written as a 2×2 matrix Hk = h⃗(k)· σ⃗, where σ⃗ = {σx, σy} is a vector of Pauli matrices and h⃗(k) = {t1 + t2 cos(kR), t2 sin(kR)}. Note that ||h⃗(k)|| is non-zero throughout the Brillouin zone, indicating that the two bands of the SSH model never touch. As k traverses the Brillouin zone in 1D, going from −π/R to π/R, h⃗(k) goes around a circle centered at t1 with radius t2. The number of times this trajectory winds around the origin is a topological invariant, being either 0 if t1 < t2 or 1 if t1 > t2, corresponding to the two distinct topological phases of the SSH chain. The winding number is directly related to the Zak phase , which is the 1D equivalence of the Berry phase in higher dimensions [33,34]. The existence of a topologically protected mode at the boundary of two topologically distinct domains is prescribed by the Atiyah-Singer index theorem [35].A domain wall in the single-bond double-bond pattern of the chain, as depicted in Fig. 1(a), is a topological defect: locally changing the bond strength around the kink cannot eliminate it as the single- double-bond pattern is a non-local property. Domain walls can, however, be eliminated in pairs by moving them towards each other, as the single-double-bond pattern far away from where the domain walls are being merged is not effected by the merger. Moreover, each domain wall must host a mid-gap state, a topological defect mode, that is localized in the vicinity of the domain wall.
We now investigate the optical equivalent of the SSH topological defect modes in plasmonic crystals. Consider a plasmonic crystal that consists of an array of parallel nanowires with staggered spacing as depicted in Fig. 1(b). The Helmholtz-Schrödinger analogy tells us that for each eigenmode of the SSH Hamiltonian Eq. (2), there is an equivalent electromagnetic mode in the plasmonic crystal. We note that the fermionic commutation relations do not play a role here as we are considering the non-interacting case; hence, we can replace the operator that creates an electron in a carbon atom atomic orbital by the operator that create a surface plasmon on the i-th nanowire. (We present the details of the connection between the continuous and discrete Hamiltonians in the appendix.) Thus the band gap in the electronic system maps onto a βz gap in the plasmonic system. Moreover, electronic states that are localized on domain walls (and appear inside the band gap) map directly onto guided plasmonic states that propagate along the domain walls in the z-direction [and appear in the βz gap, see Fig. 1(c)].
Now consider injecting a spatially truncated plane wave into a plasmonic crystal with staggered spacing but no domain walls [Fig. 1(d)]. This is equivalent to injecting an electron into polyacetylene using a local probe like an STM tip. Because the electron is being injected locally, it overlaps many k-modes, and hence the electron wave-packet will spatially spread out as time advances. Similarly, the plasmonic wave packet will expand in the transverse direction as it advances along the z-direction. To illustrate this expansion we introduce the normalized Poynting vector Pz in the x–z plane that cuts through the middle of the nanowire array
where S⃗ is the non-normalized Poynting vector and ẑ is the unit vector along the z-direction. Figure 1(e) shows the spreading out of a plasmonic wave-packet obtained using 3D FDTD simulation of the plasmonic crystal. Completing the analogy, the group velocity of the electron in polyacetylene corresponds to the opening angle of the light cone in the plasmonic crystal.The introduction of a topological defect into the plasmonic crystal [see Fig. 1(b)] gives rise to a localized mode. We expect that light injected in the vicinity of the topological defect couples both to the localized defect mode as well as to the bulk modes. We plot the results of injecting a truncated plane wave into a plasmonic crystal with a topological defect [Fig. 1(f)] in Fig. 1(g). In accord with our expectations, we observe that light coupled into the bulk modes forms a diffracting fan, while light coupled into the topological mode propagates without spreading transversely.
We comment that light guided by a topological defect mode is spatially concentrated. As a figure of merit, we consider the quantity λ2Pz, where λ is the free space wavelength, which measures how much the light is squeezed spatially as compared with diffraction limited optics (λ2Pz ≈ 1 at diffraction limit). Plasmonic confinement of light in the topologically guided mode of the structure depicted in Fig. 1(g) results in λ2Pz ≈ 16.6.
3. Manipulation of light using topological defect modes
In the electronic system, there are two well-established operations for manipulating topological defect modes: (i) shifting the position of a topological domain wall causes the associated topological defect mode to be carried along with the domain wall, and (ii) pairs of domain walls can be nucleated and pulled apart, causing two of the bulk modes to be turned into topological defect modes. In this section, we explicitly apply these operations to achieve topological manipulation of light in the plasmonic crystal of metal nanowires.
In the Helmholtz-Schrödinger correspondence, ∂t is mapped to c∂z, so that time dependent manipulation in the electronic system is mapped onto axial dependent manipulation in the plasmonic crystal. For instance, the time-dependent shifting of domain walls in the electronic picture is mapped to the shifting of the domain walls as a function of the axial position z in the plasmonic crystal, which we achieve through the variation of the nanowire spacing as a function of z.
3.1. Topological spatial mode filter
Nucleating two domain walls in the middle of the SSH chain and adiabatically moving them apart results in a spectral flow in which a pair of delocalized electronic states, one from the upper bulk band and one from the lower, are adiabatically transformed into the two mid-gap states spatially localized on the domain walls. We take advantage of this spectral flow to perform spatial mode filtering of light using a plasmonic crystal.
Specifically, we design a plasmonic crystal with a pair of domain walls that are created inside the array and shifted apart as a function of z, see Fig. 2(a). Two of the bulk modes at the input side (z = 0) of the array are mapped into the two topological defect modes on the output side (z = 200 μm) by the spectral flow, see Fig. 2(b). If one of these two bulk modes is injected into the array, after propagation, the maximum of the light intensity on the output side will be strongly localized around the kinks. To verify this behavior, we perform full 3D FDTD simulations of the structure depicted in Fig. 2(a), and plot the results in Fig. 2(c). We observe that the majority of the light flux is indeed guided into the topological defect modes. On the other hand, if we inject any other bulk mode on the input side it will be rejected by the mode filter and the output light flux in the vicinity of the domain walls will be small (see appendix). We note that the topological defect modes are not always in the center of the band-gap [see Fig. 2(e)] as we are locally slightly breaking the sub-lattice symmetry. Indeed, topological protection endows our device with robustness against exactly these types of perturbations.
3.2. Beam splitter with topological defect modes
Consider the geometry of a nanowire array depicted in Fig. 2(d): two domain walls that are well separated spatially on the input and output sides of the structure are brought close together in the middle of the structure. The propagation of the two topological defect modes can be described by the Landau-Zener Hamiltonian, in which the coupling is controlled by the spatial separation of the domain walls. The spectral flow of the array along the z direction is shown in Fig. 2(e). When the domain walls are well separated, the topological defect modes are non-interacting, and hence both are in the middle of the gap. As the domain walls are moved closer, the topological defect modes begin to interact and the degeneracy is broken. By controlling the interaction strength and the length of the interaction region it is possible to construct a 50–50 beam splitter. We perform 3D FDTD simulation of the beam splitter in which we inject light into one of the defect modes on the input side and observe an equal superposition of light in the two defect modes on the output side [see Fig. 2(f)].
4. Tolerance to perturbations in wire placement and diameter
Consider a structure in which a topological defect mode is guided by a kink that is being shifted as a function of z [see Fig. 3]. To test the tolerance of the light manipulation to perturbations (e.g. manufacturing defects), we measure the amount of light that leaks from the topological defect mode into the bulk modes due to perturbations in the plasmonic nanowire structure. This is a particularly stringent test of mode transport fidelity, as perturbations near the domain wall can affect the structure of the defect mode directly as well as how the defect mode hops from nanowire-to-nanowire as the domain wall is shifted. As it is not possible to test the tolerance to all kinds of perturbations, we focus on two specific types that are likely to occur in experimental situations: meander of a single nanowire and distortions in the diameter of a single nanowire. We note that topological robustness only protect defect modes when they are moved infinitely slowly. In a realistic plasmonic system we would rather move the defect modes as fast as possible to avoid absorption losses which has the potential to break the adiabatic approximation. It is the tolerance of light manipulation rather than topological mode robustness protected by the chiral symmetry that is examined in this section.
In Figs. 3(a) and 3(b), we compare the results of 3D FDTD simulation for two test structures in which we displace one of the nanowires by the addition of a meander M(z) = δx e−(z−z0)2/w2. In Fig. 3(a) the meander is smooth with w = 25 μm, while in Fig. 3(b) the meander is abrupt with w = 6.25 μm. The positions of unperturbed wires are indicated with black lines, the original position of the perturbed wire is indicated with black dashed line and its new position with the red line. The color scale shows Pz for light that is injected into the topological defect mode on the left side of the structure. We observe that if the meander is sufficiently smooth the topological mode remains guided along the domain wall [Fig. 3(a)]. However, if the meander is too abrupt there is significant leakage of light into the bulk modes [Fig. 3(b)].
What role does topological protection play here? The spectral flow calculation (see Appendix D) shows that the topological defect mode βz shifts as M(z) becomes large. This feature is due to the existence of next-nearest-neighbor (NNN) coupling between the nanowires which breaks the sublattice symmetry. As the NNN coupling is exponentially small, the topological defect mode is well protected from nanowire meander disorder (see Appendix D for details).
Despite the topological protection from perturbations that respects sub-lattice symmetry, these perturbations can still break adiabaticity (i.e. an abrupt meander) to give heavy leakage. This is a consequence of the fact that while topological protection ensures that the defect mode is pinned to the middle of the topological band-gap, the wave function of the topological mode still depends on the perturbation. Thus if the perturbation is turned on sufficiently abruptly, on a length scale ≈ Δβz, adiabaticity will break down. To quantify the amount of leakage caused by a perturbation, we introduce the quantity η,
where E⃗α and H⃗α are the electric and magnetic fields that correspond to the guided mode (η = 1 perfect guidance, η = 0 complete leakage). In Fig. 3(c) we plot η as a function of the smoothness parameter w. We observe a dramatic loss of adiabaticity for w ≲ 10μm.Perturbations of the nanowire diameter, as opposed to position, are potentially more problematic as these explicitly break the sublattice symmetry [by affecting the plasmon self-energy]. We investigate a series of structures similar to the one depicted in Fig. 3(a), but instead of shifting the red wire, we modify its diameter D(z) = D0 − δd e−(z−z0)2/w2, where D0 = 200 nm is the unperturbed diameter. We plot η as a function of δd in Fig. 3(d). We observe robustness to small perturbations followed by a sharp drop in η for large perturbations. The loss of transport fidelity is again associated with the loss of adiabaticity. However, due to the sublattice symmetry breaking nature of the perturbation, the loss of adiabaticity occurs prematurely as the defect mode is pushed close to the edge of the topological band-gap. That is the defect mode is indeed protected from small perturbations in the nanowire diameter (and other non-avoidable symmetry-breaking perturbations) by the topological band gap.
In summary, the device is quite tolerant to perturbations in both the nanowire position and diameter. Indeed, the device is very tolerant to nanowire position errors, which preserves the chiral symmetry if the NNN interaction is ignored. In order to see strong leakage from the topological defect mode into the bulk modes we need to displace a nanowire until it is almost touching its neighbor. The device is more sensitive to nanowire diameter which breaks the chiral symmetry, however it is still tolerant to diameter errors of ∼ 30%.
5. Outlook and summary
One of the key issue of plasmonic devices is the absorption of light due to finite optical conductivity in metals and hence a small but finite imaginary part of the dielectric constant. A useful figure of merit is FOM = Δβz Ld/(2π) the product of the decay length and the topological band gap. The FOM counts how many adiabatic operations we can do on a topologically guided mode before it decays by 1/e. The proposed topological mode filter and beam splitter both require ∼ 50 operations. For silver nanostructures with λ ∼ 900 nm (in vacuum), we find that the FOM ∼ 5 − 20 depending on the details of the structure (see appendix). From the experimental perspective, these are appealing length scales due to the availability of lasers and the ease of nano-fabrication. Moreover, for demonstration purposes the attenuation is quite reasonable. Applications would require further optimization of the structure to limit attenuation.
In summary, we have provided an in-depth analysis of the plasmonic analogue of the topological protected defect modes in the SSH model theoretically. We have explored using topological defect modes to manipulate light. Specifically, we numerically demonstrated the functionality of two devices: a mode filter and a beam splitter. Moreover, we showed that light manipulation is tolerant to fabrication errors in nanowire diameter and very tolerant to fabrication errors in nanowire displacement.
Appendix
Appendix A. The paraxial Schrödinger equations
The propagation of electromagnetic waves is governed by the Helmholtz equation
where ∊(x, y, z) is the position-dependent relative permittivity, E⃗(x, y, z), B⃗(x, y, z) are the electric and magnetic field components, and the right hand side encodes the boundary conditions at the metal-air interface. For the case of dielectric waveguides, the dielectric constant tends to vary gently, and hence the terms on the right hand side of Eq. (6) can be neglected. Therefore it is natural to obtain the paraxial Schödinger equation from the Helmholtz equation for Ez. For metallic structures variations of ∊ cannot be neglected. However, for structures that are translationally invariant in z-direction, ∇∊ only has transverse components. Hence, it is natural to focus on the transverse components Ex and Ey, as the Helmholtz equations for those components have closed form. For TM-modes propagating at small angles to the z-axis it is natural to implement the paraxial approximation, which leads us to the paraxial Schrödinger equation, Eq. (1). The potential energy operator, that appears in Eq. (1), is given by where the second term accounts for the boundary conditions at the metal-air interface.Appendix B. Connecting the continuous and discrete Helmholtz equations
In this appendix, we connect the continuous description of electromagnetic waves Eq. (1) with the discrete description of the SSH model given by the Hamiltonian
This Hamiltonian is an extension of Eq. (2), where we have added the on-site energy vi which will be used to describe the plasmon self-energy. We note that this connection is only precise in the limit of weak coupling (i.e. when the distance between nanowires is sufficiently large compared to the wavelength of light). In the strong coupling limit the topological properties of the Helmholtz equation remain intact but the tight-binding model can no longer be used to accurately describe light propagation. Consequently we use full 3D FDTD solutions of the Helmholtz equation throughout the main text.Our strategy to make the connection is to (1) describe the plasmon “self-energy” by modeling a single nanowire, and (2) describe the plasmon hopping by modeling two nanowires. We begin by considering a single nanowire of the type that makes up the plasmonic crystal. The single nanowire has a well defined plasmon mode, i.e. a radially symmetric solution of Eq. (1) of the form ψ(x, y, z) = ψ1(x, y)eiβ1z. We can capture the plasmon self-energy β1, by setting vi = c(β0 − β1) in Eq. (7).
Next, we consider the plasmon modes of a system of two parallel nanowires separated by distance s. The plasmon spectrum is now composed of two modes ψ± with eigenvalues β±. These are approximately the symmetric and the antisymmetric combinations of the single nanowire modes. Comparing this spectrum with the spectrum of the two site discrete model, we identify t = (β+ − β−)/2. In Fig. 4(a) we plot βz of the symmetric and antisymmetric modes obtained using the Helmholtz equations as a function of s. We observe that the splitting of the symmetric and antisymmetric modes with respect to the single-nanowire β1 (black dash line) is even for s > 0.5 μm and hence we can extract the tight binding parameter t. For s < 0.5 μm the splitting becomes uneven signaling the breakdown of the tight binding model. We plot the extracted tight-binding parameter t = (β+ − β−)/2, which is applicable for s > 0.5 μm, in Fig. 4(b).
Appendix C. Mode filtering
In Fig. 5 we demonstrate an example of mode rejection by the mode filter. We inject a mode that is not adiabatically connected to the two topological defect modes in the structure. The light flux spreads out over the whole structure except for the area in the vicinity of the two topological defects.
Appendix D. Tolerance to fabrication error
In this set of figures we provide additional data for the same set of structures that were used to construct Fig. 3 of the main text. The data demonstrate the tolerance of the topological defect mode manipulation against perturbation in nanowire position and diameter. The top panel of Fig 6(a) shows the light flux computed using 3D FDTD in a structure with a single topological defect shifting across the nanowire array. Light is injected into the topological defect mode on the left side of the array. The light is guided in the middle of the nanowire array and remains confined to the topological defect until it exits the array on the right side. The bottom panel of Fig 6(a) shows the spectral flow of the βz spectrum as a function of position along the wire (similar to Figs. 2(b) and 2(e) of the main text). The spectrum plot shows that the defect mode is well separated from the bulk modes throughout the structure.
Next, we test tolerance to perturbations by displacing the red nanowire by a Gaussian with maximum displacement of δx and a width w. In Figs. 6(b)–6(f) we plot the light flux, wire displacement, and the spectral flow of the βz spectrum as a function of z for δx = 150 nm and w ranging from 25 μm to 6.25 μm. From the light flux plots (top panels), we observe that the light intensity is well guided when w ≳ 12.5 μm. For w ≲ 10.00 μm, light flux heavily leaks into the bulk modes. The reason for this, is the breaking of adiabaticity as w ≲ 1/Δβz.
Figure 7 is similar to Fig. 6, except we perturb the diameter D of the red nanowire. The unperturbed nanowire has a diameter of 200 nm, and we shrink it by 25.0nm to 100.0 nm in Figs. 7(a)–7(f). The perturbation profiles have a fixed width of 25 μm. Topological mode guidance starts to break down when the diameter shrinks by 87.5 nm as the βz of the guided mode is approaching the bulk spectrum.
We note that the topological defect modes in Fig. 6 are shifted slightly from the middle of the band gap when we introduce disorder on the nanowire position. This shift is caused by the exponentially weak next-nearest-neighbor (NNN) hopping of the plasmons, which breaks the sub-lattice symmetry. Indeed, the effect of the NNN interaction on the energy of the topological defect mode (in polyacetylene) has been studied before [36–39].
Appendix E. Decay of the topological defect mode
In this section, we compute the decay length of the topological defect modes using complex susceptibility for silver nanowire structures. We use Lumerical Mode Solution numerical eigenmode solver to obtain complex βz and extract the decay length for two types of topological defects (see Fig. 8 inset). Specifically, we simulate a nanowire array consisting of 7 nanowires with periodic boundary conditions. The nanowires have a diameter of either 400, 800, or 1200 nm and staggered surface-to-surface spacing of 50 nm and 350 nm. The key figure of merit is FOM = Δβz Ld/(2π). Here, Δβz is the topological band gap, which sets the length scale for adiabatic manipulation of topological defects; Ld = 1/Im[βz] is the decay length for the topologically guided mode.
In Fig. 8 we plot the FOM as a function of free space wavelength (λ0) for two flavors of topological defects: domain wall with major spacing (solid lines) and domain wall with minor spacing (dashed lines). We observe that FOM for topological defects with major spacing is larger than the FOM for defects with minor spacing. The reason for this is the large penetration of the electric field into the nanowires for the kink with minor spacing [see plots of the mode structure, Fig. 9]. We also observe that the FOM has a broad maximum around 900 nm, which defines an optimal wavelength for silver structures. The maximum is a result of the competition between the topological band gap (which decreases with increasing wavelength) and the decay length (which increases with increasing wavelength). Finally, we observe that the FOM for domain wall with major spacing is larger for larger diameter nanowires, which is a consequence of the suppression of Ld as the wire diameter becomes smaller than wavelength.
Appendix F. Phase matching to improve coupling to the topological defect mode
Finally, we make a remark regarding [24], which claims to experimentally detect topologically guided modes in the type of structures that we propose. We found that the analysis carried out in [24] is flawed. Specifically, (i) there is breakdown of the coupled-mode theory, which predicts velocities of modes that are faster than the speed of light (ii) there is zero overlap of the injected mode with the topological mode in minor-spacing defect array due to symmetry reasons.
The the self-energy in the coupled mode theory results in the real part of mode index β = 1.047k0, where k0 is the free space wavenumber. According to their coupled mode theory calculation, the lower band edge reaches βz ≈ β − 0.225 k0, which corresponds to βz ≈ 0.822k0. This value of the mode index indicates that the mode is propagating faster than the speed of light, which signals the breakdown of the coupled mode theory.
The topological protected defect modes in major-spacing case is symmetric while minor-spacing case tends to have an anti-symmetric nature [see Fig. 9]. In real experiment, if the plasmonic mode is injected directly into the beginning of the nanowire array without any modification of the plasmonic field, the mode is likely to be symmetric and fails to couple to the topological defect mode in minor-spacing case. Here we proposed a way of modifying the length of the nanowires at the beginning of the nanowire array to induce phase difference between the nanowires [see Fig. 10], in order to improve topological mode coupling efficiency. If we assume a truncated plane wave injected into the plasmonic nanowire arrays shown in Fig. 10 from above and the plasmon is excited at the edges of each nanowire, in order the achieve the π-phase shift between the defect nanowire and the nearby nanowires [see Fig. 9], the nanowire length difference should be Δl = π/β, where β is the wave-number of the plasmonic mode. The different staggered pattern in minor-spacing defect array [Fig. 10(a)] and major-spacing defect array [Fig. 10(b)] is because of the distinct phase profile of the topological defect modes in these two nanowire arrays.
Funding
Charles E. Kaufman Foundation (KA2016-85221 and KA2014-73919).
Acknowledgments
C. L. was partially supported by a Pittsburgh Quantum Institute graduate student fellowship. We acknowledge the support from the Charles E. Kaufman Foundation.
References and links
1. C. L. Kane and E. J. Mele, “Quantum spin hall effect in graphene,” Phys. Rev. Lett. 95, 226801 (2005). [CrossRef] [PubMed]
2. B. A. Bernevig and S.-C. Zhang, “Quantum spin hall effect,” Phys. Rev. Lett. 96, 106802 (2006). [CrossRef] [PubMed]
3. L. Fu and C. L. Kane, “Topological insulators with inversion symmetry,” Phys. Rev. B 76, 045302 (2007). [CrossRef]
4. A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Ludwig, “Classification of topological insulators and superconductors in three spatial dimensions,” Phys. Rev. B 78, 195125 (2008). [CrossRef]
5. A. Kitaev, “Periodic table for topological insulators and superconductors,” AIP Conference Proc. 1134, 22–30 (2009). [CrossRef]
6. M. Z. Hasan and C. L. Kane, “Colloquium : Topological insulators,” Rev. Mod. Phys. 82, 3045–3067 (2010). [CrossRef]
7. M. König, S. Wiedmann, C. Brüne, A. Roth, H. Buhmann, L. W. Molenkamp, X.-L. Qi, and S.-C. Zhang, “Quantum spin hall insulator state in hgte quantum wells,” Science 318, 766–770 (2007). [CrossRef] [PubMed]
8. C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. Das Sarma, “Non-abelian anyons and topological quantum computation,” Rev. Mod. Phys. 80, 1083–1159 (2008). [CrossRef]
9. F. Xiu, L. He, Y. Wang, L. Cheng, L.-T. Chang, M. Lang, G. Huang, X. Kou, Y. Zhou, X. Jiang, Z. Chen, J. Zou, A. Shailos, and K. L. Wang, “Manipulating surface states in topological insulator nanoribbons,” Nat. Nano. 6, 216–221 (2011). [CrossRef]
10. K. Nakayama, K. Eto, Y. Tanaka, T. Sato, S. Souma, T. Takahashi, K. Segawa, and Y. Ando, “Manipulation of topological states and the bulk band gap using natural heterostructures of a topological insulator,” Phys. Rev. Lett. 109, 236804 (2012). [CrossRef]
11. A. R. Mellnik, J. S. Lee, A. Richardella, J. L. Grab, P. J. Mintun, M. H. Fischer, A. Vaezi, A. Manchon, E.-A. Kim, N. Samarth, and D. C. Ralph, “Spin-transfer torque generated by a topological insulator,” Nature 511, 449–451 (2014). Letter. [CrossRef] [PubMed]
12. F. D. M. Haldane and S. Raghu, “Possible realization of directional optical waveguides in photonic crystals with broken time-reversal symmetry,” Phys. Rev. Lett. 100, 013904 (2008). [CrossRef] [PubMed]
13. A. B. Khanikaev, S. Hossein Mousavi, W.-K. Tse, M. Kargarian, A. H. MacDonald, and G. Shvets, “Photonic topological insulators,” Nat. Mater. 12, 233–239 (2013). [CrossRef]
14. L. Lu, L. Fu, J. D. Joannopoulos, and M. Soljacic, “Weyl points and line nodes in gyroid photonic crystals,” Nat. Photonics 7, 294–299 (2013). [CrossRef]
15. S. A. Skirlo, L. Lu, and M. Soljačić, “Multimode one-way waveguides of large chern numbers,” Phys. Rev. Lett. 113, 113904 (2014). [CrossRef] [PubMed]
16. A. Poddubny, A. Miroshnichenko, A. Slobozhanyuk, and Y. Kivshar, “Topological majorana states in zigzag chains of plasmonic nanoparticles,” ACS Photonics 1, 101–105 (2014). [CrossRef]
17. L. Lu, C. Fang, L. Fu, S. G. Johnson, J. D. Joannopoulos, and M. Soljačić, “Symmetry-protected topological photonic crystal in three dimensions,” Nat. Phys. 12, 337–340 (2016). [CrossRef]
18. S. Yu, X. Piao, J. Hong, and N. Park, “Metadisorder for designer light in random systems,” Science Advances2 (2016). [CrossRef] [PubMed]
19. S. Longhi, “Non-hermitian bidirectional robust transport,” Phys. Rev. B 95, 014201 (2017). [CrossRef]
20. Z. Wang, Y. Chong, J. D. Joannopoulos, and M. Soljacić, “Observation of unidirectional backscattering-immune topological electromagnetic states,” Nature 461, 772–775 (2009). [CrossRef] [PubMed]
21. M. Hafezi, S. Mittal, J. Fan, A. Migdall, and J. M. Taylor, “Imaging topological edge states in silicon photonics,” Nat. Photonics 7, 1001–1005 (2013). [CrossRef]
22. M. C. Rechtsman, J. M. Zeuner, Y. Plotnik, Y. Lumer, D. Podolsky, F. Dreisow, S. Nolte, M. Segev, and A. Szameit, “Photonic floquet topological insulators,” Nature 496, 196–200 (2013). [CrossRef] [PubMed]
23. A. Blanco-Redondo, I. Andonegui, M. J. Collins, G. Harari, Y. Lumer, M. C. Rechtsman, B. J. Eggleton, and M. Segev, “Topological optical waveguiding in silicon and the transition between topological and trivial defect states,” Phys. Rev. Lett. 116, 163901 (2016). [CrossRef] [PubMed]
24. Q. Cheng, Y. Pan, Q. Wang, T. Li, and S. Zhu, “Topologically protected interface mode in plasmonic waveguide arrays,” Laser Photonics Rev. 9, 392–398 (2015). [CrossRef]
25. I. S. Sinev, I. S. Mukhin, A. P. Slobozhanyuk, A. N. Poddubny, A. E. Miroshnichenko, A. K. Samusev, and Y. S. Kivshar, “Mapping plasmonic topological states at the nanoscale,” Nanoscale 7, 11904–11908 (2015). [CrossRef] [PubMed]
26. W. Zhang, X. Chen, and F. Ye, “Plasmonic topological insulators for topological nanophotonics,” Opt. Lett. 42, 4063–4066 (2017). [CrossRef] [PubMed]
27. A. J. Heeger, S. Kivelson, J. R. Schrieffer, and W. P. Su, “Solitons in conducting polymers,” Rev. Mod. Phys. 60, 781–850 (1988). [CrossRef]
28. W. P. Su, J. R. Schrieffer, and A. J. Heeger, “Solitons in polyacetylene,” Phys. Rev. Lett. 42, 1698–1701 (1979). [CrossRef]
29. F. Bleckmann, Z. Cherpakova, S. Linden, and A. Alberti, “Spectral imaging of topological edge states in plasmonic waveguide arrays,” Phys. Rev. B 96, 045417 (2017). [CrossRef]
30. S. Ke, B. Wang, H. Long, K. Wang, and P. Lu, “Topological edge modes in non-hermitian plasmonic waveguide arrays,” Opt. Express 25, 11132–11143 (2017). [CrossRef] [PubMed]
31. W. Huang, C. Xu, S. T. Chu, and S. K. Chaudhuri, “The finite-difference vector beam propagation method: analysis and assessment,” J. Lightwave Technol. 10, 295–305 (1992). [CrossRef]
32. W. P. Huang and C. L. Xu, “Simulation of three-dimensional optical waveguides by a full-vector beam propagation method,” IEEE J. Quantum Electron. 29, 2639–2649 (1993). [CrossRef]
33. M. V. Berry, “Quantal phase factors accompanying adiabatic changes,” Proc. R. Soc. Lond. A 392, 45–57 (1984). [CrossRef]
34. J. Zak, “Berry’s phase for energy bands in solids,” Phys. Rev. Lett. 62, 2747–2750 (1989). [CrossRef] [PubMed]
35. M. F. Atiyah and I. M. Singer, “The index of elliptic operators: I,” Annals Math. 87, 484–530 (1968). [CrossRef]
36. X. Yi, “Note on the electron-hole symmetry breaking of the energy band structure of polyacetylene,” Commun. Theoretical Phys. 14, 209 (1990). [CrossRef]
37. Z. J. Li, Z. An, K. L. Yao, and Z. G. Li, “Effects of the second and third neighbor hopping interactions on the electronic states of soliton in polyacetylene,” Zeitschrift für Physik B Condensed Matter 92, 195–197 (1993). [CrossRef]
38. L. Li, Z. Xu, and S. Chen, “Topological phases of generalized su-schrieffer-heeger models,” Phys. Rev. B 89, 085111 (2014). [CrossRef]
39. B. Narayan, “End modes in arrays of modulated su–schrieffer–heeger chains,” Pramana 87, 25 (2016). [CrossRef]