Topological Spaser Jhih-Sheng Wu,∗ Vadym Apalkov,y and Mark I. Stockmanz Center for Nano-Optics (CeNO) and Department of Physics and Astronomy, Georgia State University, Atlanta, Georgia 30303 (Dated: October 8, 2019) We theoretically introduce a topological spaser, which consists of a hexagonal array of plasmonic metal nanoshells containing an achiral gain medium in their cores. Such a spaser can generate two mutually time-reversed chiral surface plasmon modes in the K- and K0-valleys, which carry the opposite topological charges, ±1, and are described by a two-dimensional E0 representation of the D3h point symmetry group. Due to the mode competition, this spaser exhibits a bistability: only one of these two modes generates, which is a spontaneous symmetry breaking. Such a spaser can be used for an ultrafast all-optical memory and information processing The concept of the surface plasmon amplification by of two different sublattices is provided by the transition stimulated emission of radiation (spaser, also called plas- metal dichalcogenide (TMDC) crystals [28{30]. monic nanolaser) [1{3] has recently been experiencing The spasing eigenmodes are surface plasmons (SPs) rapid development. Many different types of spasers that should be classified corresponding to irreducible rep- have been proposed [4{7] and demonstrated [8{18]. The resentations of the symmetry point group of the lattice spasers were also applied to various problems includ- unit cell, which is D3h [31]. This group has six represen- ing explosives detection [19], monitoring of the nano- tations, of which E0 is a two-dimensional (doublet) rep- environment [20, 21], cancer therapeutics and diagnostics resentation with the desired properties { see Appendix. (theranostics) [22]. It describes two degenerate modes time-reversed to each A particular type of the spasers is represented by plas- other, which carry topological charges of QT = ±1 defin- monic crystals that include gain media [7, 13, 23, 24]. ing their chirality; their local fields rotate in time and Such spasers belong to the class of lasing spasers [7], space in the opposite directions. which are nanostructured plasmonic metasurfaces con- For the proposed topological spaser, we will show 0 sisting of a periodic lattice of individual spasers. Due to that the degenerate eigenmodes at the K- and K -points interactions in the near-field, the individual spasers lock strongly compete with each other, making such a topo- in phase to generate temporaly- and spatially-coherent logical spaser bistable: either QT = 1 or QT = −1 mode fields. However, such a fundamental question as the can generate. These spasing modes are topological Berry effects of topological properties (the Berry curvature) plasmons [32]. However, in Ref. 32, the system's chirality [25, 26] of the plasmonic Bloch bands of such crystals was due to the induced valley polarization. In a sharp has not yet been investigated. contrast, in our case the system is originally achiral and Recently, a groundbreaking work has been carried out T -reversible: no integral Berry curvature or magnetic aimed at obtaining a topological lasing in a plasmonic- field are present; the chirality is self-organized due to the mode competition causing a spontaneous violation of the photonic (diffractive) lattice of honeycomb symmetry 0 [27]. The plasmonic lasing was observed at the K-points T -reversal and σv-reflection symmetries. 0 The quasistatic SPs eigenmodes [33] are electric poten- but only in a mode of the A1 symmetry, which is a sin- glet (scalar) representation. Therefore, this mode does tials, 'νk(r), characterized by lattice momentum k and not possess a chiral topological charge. band index ν. These eigenmodes can be found from the quasistatic equation [34] In this Letter, we propose a topological spaser that 2 geometrically is a deeply-subwavelength two-dimensional r [Θ(r)r'νk(r)] = sνkr 'νk(r); (1) (2d) crystal (metasurface) with a honeycomb symmetry where 1 > sνk > 0 are the eigenvalues, and Θ(r) is the built of two different triangular sublattices, A and B { see characteristic function, which is 1 inside and 0 outside Fig. 1. A metaatom of such a lattice is a plasmonic metal the metal. These eigenmodes satisfy an orthonormality nanoshell containing an achiral gain medium, similar to R ∗ 3 condition: r'νkr'ν0k0 d r = δνν0 δkk0 . arXiv:1909.11113v2 [cond-mat.mes-hall] 4 Oct 2019 the spaser geometry of Ref. 3. The A- and B-sublattices We employ the tight-binding approximation where the of such a topological spaser differ in size and shape of the on-site states are the modes of isolated metal nanoshells, constituent nanoshells (see the caption to Fig. 1), so their and the coupling between them is the nearest-neighbor individual spasers have different eigenfrequencies. Note dipole interaction. For each nanoshell, we consider two that a natural example of a honeycomb lattice consisting eigenmodes constituting a basis for the E0 doublet repre- (α) sentation of the D3h point symmetry group, 'm / Y1m, where Y1m are the spherical harmonics, m = ±1, and ∗Electronic address: [email protected] α = A or B. These eigenmodes with m = ±1 have the yElectronic address: [email protected] electric field components in the xy-plane of the lattice zElectronic address: [email protected] and are coupled by the dipole interaction. 2 FIG. 2: (a) Band structure of the honeycomb array of metal nanoshells. Due to the broken inversion symmetry, i.e., sub- lattices A and B being inequivalent, there are band gaps at FIG. 1: Honeycomb array of metal nanoparticles. It consists the K- and K0-points. The spasing SP modes at the K- and of two inequivalent sublattices: A and B. Sublattice A consists K0-points are indicated by open circles; the arrows indicate of plasmonic metal nanoshells with the inner radius of 8 nm the direction of the rotation of the local modal fields for the and the outer radius of 14 nm while sublattice B consists corresponding valleys. (b) Profile of the electric field of an of similar nanoshells with the inner radius of 8 nm and the SP at the valence band in the K-valley at a certain instance outer radius of 16 nm. The gain medium is placed inside of time. Only the SPs at sublattice A are excited. The color the nanoshells. The lattice bond length is set as 50 nm. The bar to the right codes the modal field amplitude. primitive unit cell of the honeycomb crystal structure is shown by the red parallelogram. The supercell, which describes the periodicity of the SP's at both the K- and K0-points, see Eqs. 0 (9)- (12), is marked by the yellow parallelogram. where j is an arbitrary lattice site, and j are the nearest- neighbor sites to it. The solution of Eq. (4) produces four bands whose dis- (α) The eigenmodes 'm of an isolated nanoshell of a sub- persions are shown in Fig. 2(a). Due to the broken inver- lattice α = A; B satisfy the following equation [cf. Eq. sion symmetry owing to the sublattices A and B being (1)] inequivalent, there are band gaps opened at the K- and K0-points. h (α) (α) i (α) 2 (α) r Θ (r)r'm (r) = sm r 'm (r) : (2) The structure of SP excitations at these points is the following. At either the K-point or the K0-point, the where Θ(α)(r) is 1 inside the metal of nanoshell α and 0 SP modes are excited in the real space only on a single elsewhere, and r is relative to the center of this nanoshell. sublattice, A { see an illustration in Fig. 2(b). The SPs The solution of Eq. (2) can be expressed as an expansion at the K- and K0-points are mutually time-reversed: the over the spherical harmonics, Y1m, see Appendix. magnetic quantum numbers are m = ±1 corresponding With the known eigenmodes of the isolated nanoshells, to the topological charges QT = ±1, and the phase shifts we express the quasistatic potential for a mode with between the nearest sites on the A sublattice are ±2π=3, quantum numbers ν; k as a sum over the lattice and respectively, as illustrated in Fig. 3(a)-(d). Below we m = ±1, assume that the frequency of the optical transitions in q the gain medium is equal to the frequency of SPs in the X (α) (α) (α (α) (α) valence band at the K- and K0-points { see a schematic 'νk(r) = Cνkm sm exp ikRj 'm r − Rj ; jαm in Fig. 3(e). 0 (3) The SP's at the K- and K -points have equal frequen- (α) cies, !νK = !νK0 , as protected by the T symmetry. We where Rj is the lattice vector of the nanoshell j center (α) write down the Hamiltonian of the plasmonic system in in sublattice α. Expansion coefficients Cνkm satisfy the the second quantization as following tight-binding equations 0 y y X (α ) (α) HSP = ~!νK a^νKa^νK +a ^νK0 a^νK0 ; (6) Hαm,α0m0 (k)Cνkm0 = sνkCνkm : (4) α0m0 wherea ^y anda ^ are the SP creation and annihilation Here the nearest neighbor tight-binding Hamiltonian is νk νk operators. Here the frequency !νk is found from equa- Hαm,α0m0 = tion Re[s(!νk)] = sνk, where s(!) = d=(d − m(!)) is the Bergman spectral parameter, d and m(!) are the 8 (α) 0 >sm δmm0 ; α = α ; dielectric functions of the host material and the metal, > >q q respectively. The corresponding electric field is given by < (α) (α0)X α α sm sm0 exp[ik(Rj0 − Rj )]× [1] > j0 > 0 0 >R (α)∗ (α) (α ) (α ) 3 0 X y > r' r − R r' r − R 0 d r; α 6= α ; : m j m j E(r; t) = − Aνkr'νk(r)(^aνk +a ^νk); (7) (5) k=K;K0 3 FIG.