Topological Spaser

Jhih-Sheng Wu,∗ Vadym Apalkov,† and Mark I. Stockman‡ 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)∇ϕνk(r)] = sνk∇ ϕν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: ∇ϕνk∇ϕν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 †Electronic address: [email protected] electric field components in the xy-plane of the lattice ‡Electronic 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)∇ϕm (r) = sm ∇ ϕ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 † † X (α ) (α) HSP = ~ωνK aˆνKaˆνK +a ˆνK0 aˆνK0 , (6) Hαm,α0m0 (k)Cνkm0 = sνkCνkm . (4) α0m0 wherea ˆ† 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  (α) 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 †  ∇ϕ r − R ∇ϕ r − R 0 d r, α 6= α ,  m j m j E(r, t) = − Aνk∇ϕνk(r)(ˆaνk +a ˆνk), (7) (5) k=K,K0 3

FIG. 4: Dynamics of SP population per unit supercell. (a) Temporal dynamics of the SP population. Initially the SPs in both the valleys are present, which is case (iii) – see the text. The pumping rate is g = 1013 s−1. (b) The stationary number of SPs in both valleys as a function of pumping rate.

classically by replacing the creation and annihilation op- erators of the SP’s by the respective complex c-number ∗ ∗ amplitudes aνK, aνK, and aνK0 , aνK0 . The correspond- ing SP population numbers per the composite unit cell 2 2 are NK = |aνK| and NK0 = |aνK0 | . The coupled system of equations [3], which describes both the SPs and the gain medium, takes the following form X (p)∗ ˜ (p)∗ a˙ νK = −γspaνK + iNg ρ¯ ΩνK , (9) p X (p)∗ ˜ (p)∗ a˙ νK0 = −γspaνK0 + iNg ρ¯ ΩνK0 , (10) p   (p) (p) X ˜ (p) (p) n˙ = −4Im ρ¯ Ωνk aνk  + k=K,K0 0 FIG. 3: (a)-(d) The valence band SPs in the K- and K - g(1 − n(p)) − γ (1 + n(p)), (11) valleys. For both the valleys, only the SPs at the sublattice A 2 (p) (p) (p) X (p)∗ (p)∗ are excited. For the K-valley, these are the SPs with m = 1, ρ¯˙ = −Γ12ρ¯ + in Ω˜ a , (12) 0 νk νk while for the K -valley, the SP with m = −1 are generated. k=K,K0 Their phase shifts from site to site are 2π/3 and −2π/3, re- spectively. The SP local fields rotate in the directions shown where Ng is the number of electrons of the gain medium by the red arrows.(e) Schematic illustration of the system of inside each shell, γsp is the SP relaxation rate, γ2 is 0 gain medium and metal array. The SPs at the K- and K - the non-radiative decay rate of the level |2i of the chro- points compete for the same gain. mophore, Γ12 is polarization relaxation rate for the |2i → |1i transition of the chromophores, g is the ex- citation (pumping) rate per a chromophore, and Ω˜ (p) = where A = p4π s /  s0 , s0 = νk νk ~ νk d νk νk 1 (p) (p) Aνk∇ϕνk(r )d is the Rabi frequency. In compu- Re[ds(ω)/dω|ω=ωνk ]. ~ The gain medium is assumed to be achiral with a linear tations, we set: d = 2, d = 10 debye, Ng = 514, γsp = 13 −1 14 −1 12 −1 transition dipole moment d; this dipole equally couples 4.1×10 s ,Γ12 = 2.1×10 s , and γ2 = 4×10 s . to both the K and K0 SP modes. The gain medium is de- We used dielectric data [35] for silver as the nanoshells’ scribed quantum mechanically using density matrix ρ(p) metal. ˜ (p) ˜ (p) for a pth nanoshell. Within the rotating wave approx- The Rabi frequencies, ΩνK and ΩνK0 , are periodic with imation (RWA), the non-diagonal elements of the den- respect to the position vector, r(p), with the periods that (p) (p) 0 sity matrix can be written as (ρ )12 =ρ ¯ exp(iωνKt), are determined by the crystal momenta K and K , re- while the diagonal elements determine the population in- spectively. Any solution of Eqs. (9)-(12), which contains 0 (p) (p) (p) both the K and K components, will be periodic on a version, n = ρ22 − ρ11 . The interaction between the gain medium and the SP system is determined by Hamil- superlattice with an eighteen-nanoshell unit cell (“super- tonian cell”) shown in Fig. 1. Then the summation over p in Eqs. (9) and (10) is extended over the corresponding eigh- X (p) (p) Hint = − E(r , t)d . (8) teen points. p As initial conditions for Eqs. (9)-(12), we set (p) n (t)|t=0 = −1, which implies that initially all the chro- Following Ref. 3, we treat the plasmonic field quasi- mophores of the gain medium are in the ground state. 4

NK0 , decays exponentially for t & 200 fs. This occurs due to a strong competition of the K- and K0-modes for the common gain. At the same time, the decay to zero, NK0 → 0, implies that there is no cross-mode talk, which is due to the topological protection: these two SP modes carry conserved topological charges QT = ±1. Due to the T -symmetry, the valleys are symmetric: starting with the K0-valley with a dominant SP population, the final population will only be in that valley. Thus the topologi- cal spaser is a topologically-protected symmetric bistable device. The kinetics of the topological spaser, i.e., the depen- dence of the stationary SP population numbers, NK and NK0 for t → ∞, as functions of the pumping rate g, is illustrated in Fig. 4(b) for case (i) when only the K- valley has a non-zero initial population: NK = 0.01 and NK0 =0. There is a pronounced threshold after which the SP population in the K-valley grows linearly with the pumping rate. The population in the K0-valley remains zero due to the topological protection. Similarly, in case (ii), the SP population will be amplified and remain only in the K0-valley. The dynamics of the real-space distributions of the local fields in the topological spaser is shown in Fig. 6 FIG. 5: Dynamics of K-valley SP eigenmode. (a)-(d) The (a)-(e) for the case of the K-valley generation (see also distributions of the local field modulus are plotted for the unit Appendix for the K0 mode). The fields are normalized cell in the real space. The phases of the spaser oscillation are to one SP in the unit supercell. This distribution con- indicated for each panel. The color-coded scale of the local tains three dipolar fields localized on the A-sublattice fields is given at the left-hand side for one SP quantum per the nanoshells, which rotate clockwise in the optical cycle. unit supercell. The local fields rotate in time in the direction This field distribution is generally chiral and compliant of the black arrow, i.e., clockwise. (e) Edge fields for K-valley 0 mode. The instantaneous dipole vectors are indicated by the with the E representation of the D3h group. At a given green arrows; the direction of the field rotation is depicted by moment of time, this field distribution can also be de- the curved black arrows; crystal momentum k of the mode scribed as a chiral lattice wave propagating in the clock- is shown by the straight black arrow. Each such a rotating wise direction along the unit cell boundary. The 2d E0 dipole generates a current that is normal to this dipole. representation also describes a mode with the opposite (counter-clockwise) chirality, which can be obtained by 0 an application of either T or σv symmetry operations. Similarly, we assume that initially there is no polariza- The present topologically-charged E0 spasing mode is tion,ρ ¯(t)|t=0 = 0. With respect to the initial SP am- indeed dark: obviously, it does not possess a net dipole plitudes, aνK and aνK0 , we consider three cases: (i) The moment. However, there is a non-zero chiral current initial SP amplitude is small and located only at the K- propagating unidirectionally clockwise (for QT = 1) or point where we set aνK = 0.1 and aνK0 = 0; (ii) The counterclockwise (for QT = −1) within the A-sublattice same as the previous case but with the SP amplitude lo- unit cell. Such currents for the neighboring cells cancel 0 cated at the K -point, aνK = 0 and aνK0 = 0.1; (iii) The out each other. However, there will be an uncompen- SP’s are initially present at both the K- and K0-points sated chiral current at the edge of the 2d plasmonic lat- but with different amplitudes, aνK = 1 and aνK0 = 0.2. tice, which is illustrated in Fig. 6 (e). The corresponding We start with case (iii) where initially both the K- local fields resemble a field of a plasmon polariton prop- and K0-valleys are weakly populated, with a higher pop- agating clockwise along the edge despite the absence of ulation at the K-point, which emulates a random initial the edge states, similar to the edge plasmon polaritons condition created by quantum fluctuations. The pump- in Ref. [32]. ing begins at t = 0 with a rate of g = 1013 s−1. The There is the second E0 mode of the same frequency, dynamics of the SP populations in both the valleys as a which is related to the above-discussed right-rotating one 0 function of time t after the beginning of the pumping is by the application of either T or σv symmetry operations. displayed in Fig. 4(a). As we see, after the initial period That mode is generated in the K0-valley; it is rotating of decay and relaxation oscillations, the population of the counter-clockwise within the unit cell and propagating dominant K-valley starts growing exponentially and then counter-clockwise along the edge. The choice of which levels off becoming stationary with the SP population one of these modes will be generated is random and is a 0 NK ≈ 150. At the same time, the minor SP population, spontaneous violation of the T and σv symmetries. 5

In conclusion, we propose a topological spaser con- Appendix A: Dipole Modes of metallic spherical stituted by a honeycomb lattice of spherical metal shell nanoshells containing a gain medium. The two sub- lattices, A and B, are built from two different types The surface plasmon eigenmodes of the ith nanoshell of nanoshells. Such a spaser generates one of two chi- are described by the equation ral plasmonic modes characterized by topological charge QT = ±1 whose local fields rotate clockwise or counter- clockwise, respectively. Which of these two modes is gen- h i ∇ θ(i)(r)∇ϕ(i) = s(i)∇2ϕ(i). (A1) erated is determined by a spontaneous violation of sym- m m m metry as defined by the initial conditions. The chirality of the spasing mode is stable and topologically protected. Function θ(i)(r) of the nanoshell is equal to 1 inside the The macroscopic fields of the spasing mode are localized metal of the ith nanoshell and 0 elsewhere. The m = ±1 along the edge of the lattice despite the absence of the dipole eigenmode is given by edge states. These edge fields propagate clockwise (for  r QT = 1) or counter-clockwise (QT = −1). Due to a very c1 , 0 < r < R1, R1 0  2 high lattice momentum k = K or k = K , they are dark imφ r R1 ϕ = sin θe c2 , +c3 2 ,R1 < r < R2, (A2) R1 r in contrast to the original lasing spaser of Ref. [7]. Never- 2  R1 theless, they can be outcoupled by using a corresponding c4 r2 ,R2 < r. coupler, e.g., grating. Application-wise, such a topolog- (i) One finds the eigenvalue sm for these eigenmodes as ical spaser is a symmetric bistable that can be used for ultrafast optical storage and processing of information. 1 ± 1 p1 + 8η3 s(i) = 3 , (A3) m 2 Acknowledgments where

Major funding was provided by Grant No. DE-FG02- R η = 1 . (A4) 11ER46789 from the Materials Sciences and Engineering R2 Division of the Office of the Basic Energy Sciences, Office of Science, U.S. Department of Energy. Numerical simu- We impose a normalization condition lations have been performed using support by Grant No. Z DE-FG02-01ER15213 from the Chemical Sciences, Bio- (i) 2 |∇ϕm | dv = 1 . (A5) sciences and Geosciences Division, Office of Basic Energy All Space Sciences, Office of Science, US Department of Energy. The work of V.A. was supported by NSF EFRI NewLAW There are two solutions corresponding to the ± signs in Grant EFMA-17 41691. Support for J.W. came from a Eq. (A3). We choose one corresponding to the − sign, MURI Grant No. N00014-17-1-2588 from the Office of which describes a lower frequency mode. Substituting it, Naval Research (ONR). we find:

1/2 1/2 " (i) (i) 1 # " 3 p # 3 (1 − sm )(sm − ) −1 + 4η + 1 + 8η3 c = 3 = , (A6) 1 (i) p 16πR1 1 16πR 1 + 8η3 (sm − 2 ) 1 (i) 2 ! sm − c = 3 c , (A7) 2 (i) 1 sm − 1 1  −1  c = c , (A8) 3 (i) 1 3 sm − 1 ! 1 −1 c4 = c1. (A9) 3 (i) 1 sm − 3

The electric field of the SP mode inside the gain core, where er, eθ, eφ and ex, ey are the corresponding spheri- i.e., for r ≤ R1, is cal and Cartesian unit vectors, and

imφ s c1e (i) E = −A (e sin θ + e cos θ + ie ) 4π~sm R r θ φ A = . (A11) 1  s(i)0 c1 h m = −A (ex + imey), (A10) R1 6

Appendix B: SPP of Plasmonic Array The eigenvectors are orthogonal and normalized as,

Consider the problem of two media with permittivities hm|ni = δmn. (B7) d and m. The eigenmodes of SP’s in the quasistatic limit can be obtained from Consider that there are N metal particles (same metal).

2 ∇ [θ(r)∇ϕ] = s∇2ϕ, (B1) ∇ [θ(r)∇ϕn] = sn∇ ϕn, (B8) N X (i) with s = d/(d − m) and θ(r) is 1 inside metal and zero θr) = θ (r). (B9) elsewhere. The geometric properties of the eigenmodes i=1 are characterized by a generalized eigenvalue problem, (i) 2 We can solve Eq. (B8) in terms of the local basis ϕm , ∇ [θ(r)∇ϕn] = sn∇ ϕn, (B2) defined as where sn and ϕn are the eigenvalues and eigenfunctions, h (i) (i)i (i) 2 (i) respectively. The conditions on the boundary S are ∇ θ (r)∇ϕm = sm ∇ ϕm , (B10)

(i) ϕn · n∇ϕn = 0, (B3) where ϕm are the eigenmodes of the ith metal particle. We expand the total field ϕ in the local basis where n is the normal vector to S. Multiplying Eq. (B2) n by ϕn and integrating by parts, we obtain X (i) (i) ϕn = Cnmϕm , (B11) i,m sn = hn|θ(r)|ni, (B4)

|ni = ∇ϕn, (B5) (i) Z where Cnm are the coefficients to be determined. Sub- 3 ∗ stituting Eq. (B11) in Eq. (B8) and using Eq. (B10) we hm|ni = d r (∇ϕm)(∇ϕn) , (B6) obtained

X X (i) h (i)i X (i) (i) 2 (i) X (i) 2 (i) Cnm∇ θj(r)∇ϕm + sm Cnm∇ ϕm = sn Cnm∇ ϕm , (B12) j i6=j,m i,m i,m

where the first term describes the inter-particle interac- or equivalently tion. To deal with the first term, we expand the field (i) (j) ∇ϕm over the other local eigenstates ∇ϕm0 (i 6= j). The (i) field ∇ϕm inside the jth particle can be expanded over (j) (j) the local basis ∇ϕm0 , since the local basis ∇ϕm0 forms a complete set of field inside the jth metal. X X (i) (j) 2 (j) Cnmβijmm0 sm0 ∇ ϕm0 . (B15) (i) X (j) j,m0 i6=j,m ∇ϕm = βijmm0 ∇ϕm0 . (B13) m0 The interaction can be written as

X X (i) h (j)i Cnmβijmm0 ∇ θj(r)∇ϕm0 (B14) j,m0 i6=j,m The eigenmode equation becomes

X X (i) (j) 2 (j) X (i) (i) 2 (i) X (i) 2 (i) Cnmβijmm0 sm0 ∇ ϕm0 + sm Cnm∇ ϕm = sn Cnm∇ ϕm . (B16) j,m0 i6=j,m i,m i,m

(j) Multiplying φm on the both sides, integrating by parts tain and utilizing the othogonality of the eigenmodes, we ob- X (i) (k) (k) (k) Cnm0 βikm0msm = (sn − sm )Cnm (B17) i6=k,m0 7

More compactly, of βijmm0 , i.e.,

X (k) (i) (k) βikm0msm Cnm0 = snCnm. (B18) βijmm0 = βjim0m (B19) i,m0

Equation (B18) defines an eigenproblem, which yields the is protected by the Green’s reciprocity. The interaction (k) egienvalue sn and the coefficients Cnm. The Hermiticity coefficeints βijmm0 are given by an integral

Z Z 1 (j)  (j)∗  (i) 3 1 (i)  (j)∗  (i) 3 β 0 = θ (r) ∇ϕ 0 ∇ϕ d r = θ (r) ∇ϕ 0 ∇ϕ d r (B20) ijmm (j) m m (i) m m sm0 sm Z  (j)∗  (i) 3 = ∇ϕm0 ∇ϕm d r. (B21) All Space

(k) 0 Because matrix β 0 sm in Eq. (B18) is not Hermitian, E σh 3C3 3S3 3U2 3σv ikm m 0 the eigenvectors are not orthogonal, and the normaliza- E 2 2 -1 -1 0 0 P (i) 2 tion condition is not |Cnm| = 1. In actuality, the 0 i,m TABLE I: Characters of the E representation of the D3h normalization condition is given by point symmetry group, adapted from Ref. 31. The elements of D are the identity E, the horizontal-plane reflection σ , X (k)∗ (i) 3h h 0 Cnm βikm mCnm0 = 1. (B22) the three rotations C3, the three rotary-reflections S3, the 0 i,k,m ,m three rotations about the horizontal axes U2, and the three 0 vertical-plane reflections σv. We can convert Eq. (B18) into a Hermitian eigenproblem q ˜(i) (i) (i) by a substitution, Cnm = sm Cnm. With this substi- given lattice site and in space from a site to a site. A basis tution, we rewrite Eq. (B18) as for such a field can be constituted, e.g., by two spherical harmonics, Y corresponding to left- and right-handed X ˜(i) ˜(k) 1±1 Hikm0mCnm0 = snCnm , (B23) rotations. i,m0 An example of the local fields of the two eigenmodes, q (i) (j) which are mutually time-reversed and should be de- H 0 = s s 0 β 0 . (B24) ikm m m m ikm m scribed by the same doublet representation, is provided by our solution of the spaser equations – see the main Because Hikm0m is Hermitian, the eigenvectors are or- thogonal and the normalization is simply text of the article. Here we show such two modes: the K-mode rotating right [Fig. 6] and the T -reversed K0- X |C˜(i) |2 = 1. (B25) mode rotating left [Fig. 7]. nm To unambiguously identify the corresponding represen- i,m tation, we note that the underlying material system is a sublattice A whose unit cell has the third-oder vertical (z) Appendix C: Characters of Topological SP axis (3C3 symmetry elements), three second-order axes Eigenmodes as Irreducible Representations of the (in the xy plane) (3U2 symmetries), also 3S3 elements 0 D3h Point Symmetry Group (mirror-rotational symmetries), and 3σv elements (reflec- tions in the vertical planes). Spherical harmonics Y1±1 The spasing eigenmodes are surface plasmons (SPs) are spherical representations of a polar vector in the xy that should be classified corresponding to irreducible rep- plane: resentations of the point symmetry group of the lattice 1r 3 unit cell, which is D3h [31]. This group has six such rep- Y1±1 = ∓ (x ± iy) . (C1) resentations: four of them are singlet or one-dimensional 2 2π 0 0 (1d), which are A1 (scalar or fully-symmetric), A2 (a The character of the trivial E operation is always the 00 pseudovector in the z-direction), A1 (pseudoscalar), and representation dimensionality, i.e., χ(E) = 2. To find 00 A2 (a polar vector in the z-direction), and two are dou- χ(σh) for our eigenvectors we note that the σh opera- blet or 2d-representations: E0 (a polar vector in the xy- tion is reflection in the xy mirror plane (i.e., a trans- plane) and E00 (a pseudovector in the xy-plane). We seek formation x → x, y → y, and z → −z), which leaves for a solution where the local electric field is in the plane both of the basis functions of Eq. (C1) unchanged. Cor- of the lattice (the xy-plane) and rotates in time for a respondingly, χ(σh) = 2. Operation C3 is a rotation 8

by a 2π/3 angle whose character, obviously, is χ(C3) = 2 cos (2π/3) = −1. An operation S3 is a product of σh, which leaves the basis unchanged, and C3; thus its char- acter is χ(S3) = χ(C3) = −1. A U2 rotation, say about the y axis of symmetry, causes y → y, x → −x; this results in Y11 ↔ Y1−1. Thus the basis functions are exchanged, and the corresponding transformation ma- trix is purely off-diagonal; correspondingly, χ(U2) = 0. 0 An operation σv acts for our purposes exactly like U2; 0 thus χ(σv) = 0. These characters that we have obtained 0 are exactly those of the E representation of the D3h group summarized in Table I. This unambiguously iden- tifies that our spasing eigenmodes transform accordingly 0 FIG. 6: Dynamics of K-valley SP eigenmode. (a)-(d) The dis- to the E doublet representation of the D3h point sym- tributions of the local field modulus are plotted for the unit metry group. Note that that the two components of this cell in the real space. The phases of the spaser oscillation are doublet describe field distribution rotating in time and indicated for each panel. The color-coded scale of the local space in the opposite directions. They transform one fields is given at the left=hand side for one SP quantum per into another by either of the following transformations: the unit supercell. The local fields rotate in time in the di- 0 T , σ , or U2. rection of the black arrow, i.e., clockwise. (e) Edge fields for v K-valley mode. The instantaneous dipole vectors are indi- cated by the green arrows; the direction of the field rotation is depicted by the curved black arrows; crystal momentum k of the mode is shown by the straight black arrow. Each such a rotating dipole generates a current that is normal to this dipole.

FIG. 7: Dynamics of the K0-valley SP eigenmode. Same as Fig. 6 but the dynamics is T -reversed: fields rotate left and propagate left.

[1] D. J. Bergman and M. I. Stockman, Surface plasmon am- [5] D. Y. Fedyanin, Toward an electrically pumped spaser, plification by stimulated emission of radiation: Quantum Opt. Lett. 37, 404 (2012). generation of coherent surface plasmons in nanosystems, [6] D. G. Baranov, A. P. Vinogradov, A. A. Lisyansky, Y. M. Phys. Rev. Lett. 90, 027402 (2003). Strelniker, and D. J. Bergman, Magneto-optical spaser, [2] M. I. Stockman and D. J. Bergman, Surface plasmon am- Opt. Lett. 38, 2002 (2013). plification by stimulated emission of radiation (spaser), [7] N. I. Zheludev, S. L. Prosvirnin, N. Papasimakis, and US Patent 7,569,188 (2009). V. A. Fedotov, Lasing spaser, Nat. Phot. 2, 351 (2008). [3] M. I. Stockman, The spaser as a nanoscale quantum [8] M. A. Noginov, G. Zhu, A. M. Belgrave, R. Bakker, V. M. generator and ultrafast amplifier, Journal of Optics 12, Shalaev, E. E. Narimanov, S. Stout, E. Herz, T. Sutee- 024004 (2010). wong, and U. Wiesner, Demonstration of a spaser-based [4] K. Li, X. Li, M. I. Stockman, and D. J. Bergman, Sur- nanolaser, Nature 460, 1110 (2009). face plasmon amplification by stimulated emission in [9] R. F. Oulton, V. J. Sorger, T. Zentgraf, R.-M. Ma, nanolenses, Phys. Rev. B 71, 115409 (2005). C. Gladden, L. Dai, G. Bartal, and X. Zhang, Plas- 9

mon lasers at deep subwavelength scale, Nature 461, 629 L. Dai, and R. M. Ma, Lasing enhanced surface plasmon (2009). resonance sensing, Nanophotonics 6, 472 (2017). [10] R.-M. Ma, R. F. Oulton, V. J. Sorger, G. Bartal, [22] E. I. Galanzha, R. Weingold, D. A. Nedosekin, M. Sa- and X. Zhang, Room-temperature sub-diffraction-limited rimollaoglu, J. Nolan, W. Harrington, A. S. Kuchyanov, plasmon laser by total internal reflection, Nat. Mater. 10, R. G. Parkhomenko, F. Watanabe, Z. Nima, A. S. Biris, 110 (2010). A. I. Plekhanov, M. I. Stockman, and V. P. Zharov, [11] R. A. Flynn, C. S. Kim, I. Vurgaftman, M. Kim, J. R. Spaser as a biological probe, Nat. Commun. 8, 15528 Meyer, A. J. M¨akinen,K. Bussmann, L. Cheng, F. S. (2017). Choa, and J. P. Long, A room-temperature semiconduc- [23] W. Zhou, M. Dridi, J. Y. Suh, C. H. Kim, D. T. Co, tor spaser operating near 1.5 micron, Opt. Express 19, M. R. Wasielewski, G. C. Schatz, and T. W. Odom, Las- 8954 (2011). ing action in strongly coupled plasmonic nanocavity ar- [12] M. J. H. Marell, B. Smalbrugge, E. J. Geluk, P. J. van rays, Nature Nano 8, 506 (2013). Veldhoven, B. Barcones, B. Koopmans, R. N¨otzel,M. K. [24] M. I. Stockman, Lasing spaser in two-dimensional plas- Smit, and M. T. Hill, Plasmonic distributed feedback monic crystals, NPG Asia Mater. 5, e71 (2013). lasers at telecommunications wavelengths, Opt. Express [25] M. V. Berry, Quantal phase factors accompanying adia- 19, 15109 (2011). batic changes, Proc. Royal Soc. London Ser. A 392, 45 [13] F. v. Beijnum, P. J. v. Veldhoven, E. J. Geluk, M. J. A. (1984). d. Dood, G. W. t. Hooft, and M. P. v. Exter, Surface [26] D. Xiao, M.-C. Chang, and Q. Niu, Berry phase effects plasmon lasing observed in metal hole arrays, Phys. Rev. on electronic properties, Reviews of Modern Physics 82, Lett. 110, 206802 (2013). 1959 (2010). [14] Y.-J. Lu, C.-Y. Wang, J. Kim, H.-Y. Chen, M.-Y. Lu, Y.- [27] R. Guo, M. Necada, T. K. Hakala, A. I. Vakevainen, and C. Chen, W.-H. Chang, L.-J. Chen, M. I. Stockman, C.- P. Torma, Lasing at k points of a honeycomb plasmonic K. Shih, and S. Gwo, All-color plasmonic nanolasers with lattice, Phys. Rev. Lett. 122, 013901 (2019). ultralow thresholds: Autotuning mechanism for single- [28] Y. M. You, X. X. Zhang, T. C. Berkelbach, M. S. Hy- mode lasing, Nano Lett. 14, 43814388 (2014). bertsen, D. R. Reichman, and T. F. Heinz, Observation of [15] Q. Zhang, G. Li, X. Liu, F. Qian, Y. Li, T. C. Sum, C. M. biexcitons in monolayer wse2, Nat. Phys. 11, 477 (2015). Lieber, and Q. Xiong, A room temperature low-threshold [29] K. S. Novoselov, A. Mishchenko, A. Carvalho, and ultraviolet plasmonic nanolaser, Nat. Commun. 5, 4953 A. H. C. Neto, 2d materials and van der Waals het- (2014). erostructures, Science 353, 461 (2016). [16] B. T. Chou, Y. H. Chou, Y. M. Wu, Y. C. Chung, W. J. [30] D. N. Basov, M. M. Fogler, and F. J. G. de Abajo, Po- Hsueh, S. W. Lin, T. C. Lu, T. R. Lin, and S. D. Lin, laritons in van der waals materials, Science 354, aag1992 Single-crystalline aluminum film for ultraviolet plasmonic (2016). nanolasers, Sci. Rep. 6, 19887 (2016). [31] L. D. Landau and E. M. Lifshitz, Quantum Mechanics: [17] C.-J. Lee, H. Yeh, F. Cheng, P.-H. Su, T.-H. Her, Y.- Non-Relativistic Theory (Pergamon Press, Oxford and C. Chen, C.-Y. Wang, S. Gwo, S. R. Bank, C.-K. Shih, New York, 1965). and W.-H. Chang, Low-threshold plasmonic lasers on [32] J. C. W. Song and M. S. Rudner, Chiral plasmons with- a single-crystalline epitaxial silver platform at telecom out magnetic field, Proceedings of the National Academy wavelength, ACS Photonics 4, 1431 (2017). of Sciences 113, 4658 (2016). [18] S. Sun, C. Zhang, K. Wang, S. Wang, S. Xiao, and [33] D. J. Bergman and D. Stroud, Properties of macro- Q. Song, Lead halide perovskite nanoribbon based uni- scopically inhomogeneous media, in Solid State Physics, form nanolaser array on plasmonic grating, ACS Photon- Vol. 46, edited by H. Ehrenreich and D. Turnbull (Aca- ics 4, 649656 (2017). demic Press, Boston, 1992) pp. 148–270. [19] R.-M. Ma, S. Ota, Y. Li, S. Yang, and X. Zhang, Explo- [34] M. I. Stockman, S. V. Faleev, and D. J. Bergman, Lo- sives detection in a lasing plasmon nanocavity, Nature calization versus delocalization of surface plasmons in Nanotechnology 9, 600 (2014). nanosystems: Can one state have both characteristics?, [20] S. Wang, B. Li, X. Y. Wang, H. Z. Chen, Y. L. Wang, Phys. Rev. Lett. 87, 167401 (2001). X. W. Zhang, L. Dai, and R. M. Ma, High-yield plas- [35] P. B. Johnson and R. W. Christy, Optical constants of monic nanolasers with superior stability for sensing in noble metals, Phys. Rev. B 6, 4370 (1972). aqueous solution, ACS Photonics 4, 1355 (2017). [21] X. Y. Wang, Y. L. Wang, S. Wang, B. Li, X. W. Zhang,