Arxiv:1812.01183V2 [Cond-Mat.Str-El] 10 Dec 2018 Tically to a Photonic Mass [26–28] in the Low-Energy Physics Talline Structure

Non-local topological electromagnetic phases of matter

Todd Van Mechelen and Zubin Jacob∗ Birck Nanotechnology Center and Purdue Quantum Center, Department of Electrical and Computer Engineering, Purdue University, West Lafayette, Indiana 47907, USA

In 2+1D, nonlocal topological electromagnetic phases are defined as atomic-scale media which host photonic monopoles in the bulk band structure and respect bosonic symmetries (e.g. time-reversal T 2 = +1). Additionally, they support topologically protected spin-1 edge states, which are fundamentally different than 1 1 spin- ⁄2 and pseudo-spin- ⁄2 edge states arising in fermionic and pseudo-fermionic systems. The striking feature of the edge state is that all electric and magnetic field components vanish at the boundary - in stark contrast to analogs of Jackiw-Rebbi domain wall states. This surprising open boundary solution of Maxwell’s equations, dubbed the quantum gyroelectric effect [Phys. Rev. A 98, 023842 (2018)], is the supersymmetric partner of the topological Dirac edge state where the spinor wave function completely vanishes at the boundary. The defining feature of such phases is the presence of temporal and spatial dispersion in conductivity (the linear response function). In this paper, we generalize these topological electromagnetic phases beyond the continuum approximation to the exact lattice field theory of a periodic atomic crystal. To accomplish this, we put forth the concept of microscopic photonic band structure of solids - analogous to the traditional theory of electronic band structure. Our definition of topological invariants utilizes optical Bloch modes and can be applied to naturally occurring crystalline materials. For the photon propagating within a periodic atomic crystal, our theory shows that besides the Chern invariant C ∈ Z, there are also symmetry-protected topological (SPT) invariants ν ∈ ZN which are related to the cyclic point group CN of the crystal ν = C mod N. Due to the rotational symmetries of light R(2π) = +1, these SPT phases are manifestly bosonic and behave very differently from their fermionic counterparts R(2π) = −1 encountered in conventional condensed matter systems. Remarkably, the nontrivial bosonic phases ν 6= 0 are determined entirely from rotational (spin-1) eigenvalues of the photon at high-symmetry points in the Brillouin zone. Our work accelerates progress towards the discovery of bosonic phases of matter where the electromagnetic field within an atomic crystal exhibits topological properties.

I. INTRODUCTION long wavelength approximation to the exact lattice field theory of optical Bloch waves. In this regime, we must consider From a material science standpoint, all known topological not only the first spatial component σxy(ω, k) = σxy(ω, k, 0) phases of matter to date have been characterized by electronic but all spatial harmonics of the crystal g 6= 0, to infinite order, phenomena [1, 2]. This is true for both time-reversal broken Hall X phases - often called Chern insulators [3–6] and time-reversal Jx (ω, k) = σxy(ω, k, g)Ey(ω, k + g). (1) unbroken phases - known as topological insulators [7–9]. The g signature of time-reversal broken phases is the quantum Hall 2 g · R ∈ 2π R conductivity σxy = ne /h, which is quantized in terms of Z are the reciprocal lattice vectors and is the the electronic Chern invariant n ∈ Z [10–12]. e being the primitive vector of the crystal. In this case, Ei is the micro- elementary charge of the electron and h the Planck constant. scopic electric field. The electromagnetic field must be de- Only recently has the idea of bosonic Hall conductivity and scribed to the same scale as the electronic wave functions, topological bosonic phases been put forth [13–21]. i.e. for photon momenta on the order of the lattice constant However, it should be emphasized that the traditional Hall ka = π, with a ≈ 5 A.˚ Since topological invariants are conductivity [22, 23] only has topological significance, with fundamentally global properties, these astronomically deep respect to the electron, in the static ω = 0 and local k = 0 lim- subwavelength fields actually play a role in the topological 2 its of the electromagnetic field σxy(0, 0) = ne /h. At high physics. frequency ω 6= 0 and short wavelength k 6= 0, the Hall con- The idea of lattice topologies in electromagnetism was first ductivity σxy(ω, k) acquires new physical meaning. We have proposed by Haldane [29, 30] in the context of photonic crys- shown that the electromagnetic field itself becomes topolog- tals [31–38]. These are artificial materials composed of two or ical [24, 25] and nonlocal Hall conductivity functions iden- more different constituents which form a macroscopic crys-

arXiv:1812.01183v2 [cond-mat.str-el] 10 Dec 2018 tically to a photonic mass [26–28] in the low-energy physics talline structure. A few important examples are gyrotropic ω ≈ 0. These topological electromagnetic phases of matter photonic crystals [31–33], Floquet topological insulators [39– depend on the global behavior of σxy(ω, k), over all frequen- 41] and bianisotropic metamaterials [42–56]. Instead, we fo- cies and wave vectors. cus on the microscopic domain and utilize the periodicity of As of yet, only the continuum topological theory of the the atomic lattice itself. Thus, the topological invariants in our aforementioned quantum gyroelectric effect (QGEE) has been theory are connected to the microscopic atomic lattice and not solved [24, 25]. Our goal is to extend this concept beyond the artiﬁcially engineered macroscopic structures. We stress that in the microscopic case, the electromagnetic theory is manifestly bosonic [57–60] (e.g. time-reversal T 2 = +1) and characterizes topological phases of matter fundamentally dis- ∗ [email protected] tinct from known fermionic and pseudo-fermionic phases. 2

With that in mind, this paper is dedicated to solving two conductivity) and show unidirectional, completely transverse longstanding problems, which is of interest to both photonics electro-magnetic (TEM) edge states [24, 25]. However, time- and condensed matter physics. The first, is developing the rig- reversal symmetric topological phenomena can arise in higher orous theory of optical Bloch modes in natural crystal solids. dimensional systems in the context of nonlocal magnetoelec- This problem gained significant interest in the 60’s and 70’s tricity [84]. These time-reversal symmetric phases possess in the context of spatial dispersion (nonlocality) as it lead to counter-propagating TEM edge states and are interpreted as qualitatively new phenomena - such as natural optical activity two copies of a bosonic Chern insulator. Features of topo- (gyrotropy) [61–64]. The current paper builds on our recent logical phenomena, such as spin-momentum locking [85–89], discovery of the quantum gyroelectric effect [24, 25] where have also been reported in conventional surface state prob- we have shown that nonlocality is also essential for topologi- lems - surface plasmon-polaritons (SPPs), Dyakonov waves, cal phenomena and is a necessary ingredient in any long wave- etc. However, these traditional surface properties are not con- length theory. However, since topological field theories are nected to any topologically protected edge states or nontrivial global constructs, a complete picture can only be achieved in phases. the microscopic domain of Bloch waves. Most of the founda- Note: Due to the frequency of integral formulas, all dif- tions have been summarized by Agronovich and Ginzburg in ferential elements are assumed√ to be correctly normalized, their seminal monograph on crystal optics [65]. Nevertheless, such as√ Fourier transforms dω/ 2π → dω and Fourier series topological properties have never been tackled to date and a eig·r/ V → eig·r, where V is the unit cell area of a crystal. few fundamental quantities, such as the Bloch energy density, have not been defined. This leads to the second problem - deriving the electro- II. LATTICE ELECTROMAGNETISM magnetic topological invariants of these systems given only the atomic lattice. We solve this problem and also provide a A. 2+1D electrodynamics systematic bosonic classification of all 2+1D topological photonic matter. Utilizing the optical Bloch modes, we show that In this paper we focus on two-dimensional materials and a Chern invariant C ∈ can be found for any two-dimensional Z the topological electromagnetic phases associated with them. crystal and characterizes distinct topological phases. We then The preliminaries for 2+1D electromagnetism can be found in go one step further and classify these topological phases with Appendix A of Ref. [24]. Conveniently, the restriction to 2D respect to the symmetry group of the crystal - the cyclic point limits the degrees of freedom of both the electromagnetic field groups C . These are known as symmetry-protected topo- N and the induced response of the material, such that strictly logical (SPT) phases [66–77] and the spin of the photon is transverse-magnetic (TM) waves propagate. The correspond- critical to their definition. The rotational symmetries of light ing wave equation reads, R(2π) = +1 impart an intrinsically bosonic nature to these phases, which are fundamentally different than their fermionic E  D  counterparts R(2π) = −1 encountered in conventional con- x x H f = i∂ g, f = E , g = D . (2) densed matter systems. We illustrate this fact by directly com- 0 t  y   y H B paring SPT bosonic and fermionic phases side-by-side. Our z z rigorous formalism of microscopic photonic band structure f is the TM polarization state of the electromagnetic field and provides an immediate parallel with the traditional theory of the material response is captured by the displacement field g. electronic band structure in crystal solids. H0(p) = p · S are the vacuum Maxwell equations in real This article is organized as follows. In Sec. II we develop space and describe the dynamics of the free photon, the general formalism of 2+1D lattice electromagnetism. First we derive the generalized linear response function - account-   0 0 −py ing for spatiotemporal dispersion to infinite order in the crys- ˆ ˆ H0(p) = pxSx + pySy =  0 0 px  . (3) tal’s spatial harmonics g. Thereafter, we find the equivalent −py px 0 Hamiltonian that governs all light-matter Bloch excitations of the material. In Sec. III we study the discrete rotational p = −i∇ is the two-dimensional momentum operator. Sˆx symmetries (point groups) of the crystal and the implications ˆ on spin-1 quantization [78–83] of the photon. The following and Sy are spin-1 operators that satisfy the angular momentum ˆ ˆ ˆ Sec. IV discusses the electromagnetic Chern number and its algebra [Si, Sj] = iijkSk, relationship to symmetry-protected topological (SPT) bosonic   phases. The bosonic classification of each phase is related di- 0 −i 0 ˆ rectly to integer quantization of the photon [Tab. I] and this Sz = i 0 0 . (4) is compared alongside their fermionic counterparts [Tab. II]. 0 0 0 Sec. V presents our conclusions. ˆ The focus of this paper is 2+1D topological electromag- Here, (Sz)ij = −iijz is the generator of rotations in the x- netic (bosonic) phases of matter C 6= 0 which requires break- y plane and is represented by the antisymmetric matrix. In ing time-reversal symmetry. These bosonic Chern insulators two dimensions, Sˆz governs all rotational symmetries of the are ultimately related to nonlocal gyrotropic response (Hall electromagnetic field. 3

TABLE I. Summary of 2+1D topological electromagnetic (bosonic) phases. Symmetry-protected topological (SPT) bosonic phases exist in all cyclic point groups CN=2,3,4,6. The continuous group C∞ describes the long wavelength theory k ≈ 0. The topological phases are characterized by their Chern invariant C ∈ Z and SPT invariant ν ∈ ZN . These numbers are not independent - but intimately related by the symmetries of the crystal: ν = C mod N. ν is protected by N-fold rotational symmetry and determines the Chern number up to a factor of N N. The bosonic classification of ν represents the direct product of rotational eigenvalues (ηN ) = +1 (roots of unity) of the electromagnetic field at high-symmetry points (HSPs) in the Brillouin zone. For the spin-1 photon, this classification is more intuitively understood in terms of 2π modulo integers mN ∈ ZN , which determine the N possible eigenvalues of ηN = exp i N mN . N Point group, CN Symmetry, ZN Bosonic classification, (ηN ) = +1 Boson SPT invariant, ν = C mod N C1 --- C2 Z2 exp (i2πC/2) = η2(Γ)η2(X)η2(Y)η2(M) ν = m2(Γ) + m2(X) + m2(Y) + m2(M) mod 2 0 0 C3 Z3 exp (i2πC/3) = η3(Γ)η3(K)η3(K ) ν = m3(Γ) + m3(K) + m3(K ) mod 3 C4 Z4 exp (i2πC/4) = η4(Γ)η4(M)η2(Y) ν = m4(Γ) + m4(M) + 2m2(Y) mod 4 C6 Z6 exp (i2πC/6) = η6(Γ)η3(K)η2(M) ν = m6(Γ) + 2m3(K) + 3m2(M) mod 6 ∗ C∞ Z exp (iθC) = ηθ(0)ηθ (∞), ηθ = exp(iθm) ν = C = m(0) − m(∞)

TABLE II. Summary of 2+1D SPT fermionic phases for comparison. The fermionic classification of ν represents the direct product of N 1 rotational eigenvalues (ζN ) = −1 (roots of negative unity) of the spinor field at HSPs in the Brillouin zone. For the spin- ⁄2 electron, this 1 classification is more intuitively understood in terms of modulo half-integers mN ∈ ZN + /2, which determine the N possible eigenvalues 2π of ζN = exp i N mN . N Point group, CN Symmetry, ZN Fermionic classification, (ζN ) = −1 Fermion SPT invariant, ν = C mod N C1 --- C2 Z2 exp (i2πC/2) = ζ2(Γ)ζ2(X)ζ2(Y)ζ2(M) ν = m2(Γ) + m2(X) + m2(Y) + m2(M) mod 2 0 0 3 C3 Z3 exp (i2πC/3) = −ζ3(Γ)ζ3(K)ζ3(K ) ν = m3(Γ) + m3(K) + m3(K ) + /2 mod 3 C4 Z4 exp (i2πC/4) = −ζ4(Γ)ζ4(M)ζ2(Y) ν = m4(Γ) + m4(M) + 2m2(Y) + 2 mod 4 C6 Z6 exp (i2πC/6) = −ζ6(Γ)ζ3(K)ζ2(M) ν = m6(Γ) + 2m3(K) + 3m2(M) + 3 mod 6 ∗ C∞ Z exp (iθC) = ζθ(0)ζθ (∞), ζθ = exp(iθm) ν = C = m(0) − m(∞)

B. 2+1D linear response theory Equation (7) implies energy conservation in Hermitian systems ω0 = ω. However, a crystal is not translationally invari- 0 The effective electromagnetic properties of a material are ant in space - momentum is not conserved k 6= k. Instead, the very accurately described by a linear response theory - assum- crystal is periodic and possesses discrete translational symme- ing nonlinear interactions are negligible. This is true for low try [65, 90], 8 intensity light |f| / 10 V/m that is sufficiently weak com- M(ω, r, r0) = M(ω, r + R, r0 + R), (8) pared to the atomic fields governing the binding of the crystal itself. Our goal is to characterize the entire topological field where R is the primitive lattice vector of the crystal. This theory in this regime. With this in mind, the most general lin- admits a Fourier decomposition in the spatial harmonics of ear response of a 2D material is nonlocal in both space and the crystal g, time coordinates, 0 X 0 −ir0·g M(ω, r, r ) = Mg(ω, r − r )e , (9) t g g(t, r) = d2r0 dt0M(t, t0, r, r0)f(t0, r0). (5) ˆ ˆ−∞ with g · R ∈ 2π Z arbitrary integer combinations of the reciprocal lattice vectors. M is the response function and compactly represents the con- Due to nonlocality, it is necessary to convert to the recipro- stitutive relations in space-time, cal space,   εxx εxy χx 0 2 2 0 0 −ik·r ik0·r0 0 0 M(ω, k, k ) = d rd r M(ω, r, r )e e . (10) M(t, t , r, r ) = εyx εyy χy . (6) ¨ ζ ζ µ x y M(ω, k, k0) determines the linear transformation properties 0 Note that M is a 3 × 3 dimensional matrix and we include all of an input wave with momentum k to an output wave with possible material responses as a generalization, for instance momentum k. In a periodic crystal, the momentum is conserved up to a reciprocal vector k0 = k + g and represents a magnetism µ and magnetoelectricity χi, ζi. If the properties of the crystal are not changing temporally discrete spectrum, (no external modulation), the response function is translation- X M(ω, k, k0) = M (ω, k)δ2(k + g − k0). (11) ally invariant in time, g g 0 0 0 0 M(t, t , r, r ) = M(t − t , r, r ) δ2(k + g − k0) is the momentum conserving delta function. 0 (7) Each Fourier element of the response function M (ω, k) de- = dωM(ω, r, r0)e−iω(t−t ). g ˆ termines the polarization dependent scattering amplitude from 4 k + g → k. These are essentially the photonic structure fac- with M(0, k, k0) = M¯ (0, k, k0) at zero frequency. To en- tors of the two-dimensional crystal. sure the electromagnetic field is real-valued, i.e. represents a In this case, k is the crystal momentum and is only uniquely neutral particle, we always require the reality condition, defined within the Brillouin zone (BZ). Hence, the electro- 0 ∗ 0 magnetic eigenstates of the medium are Bloch waves, M(ω, k, k ) = M (−ω, −k, −k ). (18)

2 0 0 Furthermore, the response is transparent at high frequency H (k)f = ω d k M(ω, k, k )f 0 0 k ˆ k ω → ∞, as the material cannot respond to sufficiently fast X (12) temporal oscillations, = ω Mg(ω, k)fk+g. g 0 1 2 lim M(ω, k, k ) = 3δk−k0 . (19) ω→∞ H (k) = k · S are the vacuum Maxwell equations in mo- 0 1 2 2 0 mentum space. The Bloch photonic wave function f(k, r) = 3 is the 3 × 3 identity matrix and δk−k0 = δ (k − k ) is the momentum conserving delta function. Lastly, the response hr|fki corresponds to the net propagation of all k+g scattered waves in the medium, must be causal and satisfy the Kramers-Kronig relations. Combining all the above criteria, we find that M can al- X ig·r ways be decomposed as a discrete summation of oscillators f(k, r) = fk+ge , (13) g [29, 52, 91], † C C 00 0 where f(k, r + R) = f(k, r) is periodic in the crystal lattice. 0 1 2 X 2 00 αk00k αk k M(ω, k, k ) = 3δk−k0 − d k . ˆ ω 00 (ω − ω 00 ) Note that Eq. (12) and (13) reduce to the continuum theory α αk αk [24, 25] when considering only the 0th order harmonic g = 0. (20) Any Hermitian (lossless) response function can be expressed in this form. Equation (20) is easily extended to 3D materi- C. Generalized response function als but our focus is on 2D topological field theories. In this case, α labels an arbitrary bosonic excitation in the material, Nevertheless, Eq. (12) poses a few serious problems; it does such as an exciton or phonon, which couples linearly to the not represent a proper first-order in time Hamiltonian since all electromagnetic fields via the 3 × 3 tensor, harmonics of the response function Mg(ω, k) depend on the eigenvalue ω. Moreover, it is not evident that the Bloch waves 0 2 2 0 0 −ik·r ik0·r0 Cα(k, k ) = d rd r Cα(r, r )e e . (21) in Eq. (13) are normalizable, as the system contains complex ¨ spatial and temporal dispersion. Due to these issues, it is ad- ω is the resonant energy of the oscillator and corresponds vantageous to return to the more general form of M(ω, k, k0) αk to a first-order pole of the response function. Notice that M without assuming discrete translational symmetry. This will itself contains an integral over k00. Microscopically, this con- allow us to derive very robust properties of the response func- stitutes the overlap with the electronic momentum to infinites- tion that can also be applied to amorphous materials or qua- imally small scale k → ∞. sicrystals. Substituting Eq. (20) into Eq. (15), we can exchange the First, we demand Hermiticity, order of integration U(ω) = d2kU(ω, k) and define, M(ω, k, k0) = M†(ω, k0, k), ´ (14) 2 2 X 2 0 Cαkk0 fk0 U(ω, k) = |fk| + d k > 0, (22) such that the response is lossless. To account for normalizable ˆ (ω − ω ) α αk electromagnetic waves, the energy density must be positive definite for all ω, which is positive definite for all ω and k. Equation (22) is the generalized inner product for the electromagnetic field and 2 2 0 † ¯ 0 represents the energy density at an arbitrary frequency and U(ω) = d kd k fkM(ω, k, k )fk0 > 0, (15) ¨ wave vector. We will now show that Eq. (20) is derived from a first-order in time Hamiltonian. where M¯ describes the inner product space in a dispersive medium, ∂ D. Generalized Hamiltonian M¯ (ω, k, k0) = [ωM(ω, k, k0)] . (16) ∂ω To find the corresponding Hamiltonian, we expand the re- Notice that U(ω) = U ∗(ω) is only real-valued when M is sponse function M in terms of three-component matter oscil- Hermitian. For realistic materials, the energy density is also lators ψα. These represent internal polarization and magneti- stable at static equilibrium ω = 0, zation modes of the material,

2 2 0 † 0 2 0 U(0) = d kd k f M(0, k, k )f 0 > 0, (17) ωψ = ω ψ + d k C 0 f 0 . (23) ¨ k k αk αk αk ˆ αkk k 5

Substituting Eq. (23) and (20) into Eq. (12) we obtain, The ﬁrst two terms on the right hand side of Eq. (24) represent the vacuum equations and self-energy of the electromagnetic 2 00 2 0 X d k d k † ﬁeld. The third term is the linear coupling to the oscillators. ωfk = H0(k)fk+ Cαk00kCαk00k0 fk0 ¨ ω 00 Combining Eq. (23) and (24) into a single algebraic matrix, α αk (24) we write the generalized Hamiltonian H(k, k0) as, X 2 0 † + d k C ψ 0 . ˆ αk0k αk α

 2 P d2k00 † † †  H0(k)δ 0 + C 00 Cαk00k0 C 0 C 0 ... k−k α ωαk00 αk k 1k k 2k k 2 ´ 0 0  C1kk ω1kδk−k0 0 ... H(k, k ) =  2  , (25)  C 0 0 ω δ 0 ...  2kk 2k k−k  ......

which is manifestly Hermitian H(k, k0) = H†(k0, k). generally periodic in k, since they correspond to energy gaps We now define uk as the generalized state vector of the in the electronic band structure. However, we do not need to electromagnetic problem; accounting for the photon fk and assume this to define the optical Bloch excitations. A periodic all possible internal excitations ψαk, coupling is sufficient. The quasiparticle eigenstates of this Hamiltonian describe   fk the complete spectrum of Bloch waves, ψ1k 2 0   d k Hkk0 uk0 = ωuk, uk = ψ2k , (26) X ˆ Hg(k)unk+g = ωnkunk, ωn(k + g) = ωn(k),  .  . g (30) which is a first-order wave equation. Notice that contraction and the eigenenergies ωnk are periodic Bloch bands. n la- of uk naturally reproduces the energy density [Eq. (22)] upon bels a particular energy band of the material with its associ- summation over all degrees of freedom, ated Bloch eigenstate |unki. The total wave function |unki contains the photon |fnki and all internal degrees of freedom † 2 X 2 ukuk = |fk| + |ψαk| = U(ω, k) describing the linear response |ψnαki. This is expressed com- α pactly in the Fourier basis un(k, r) = hr|unki, 2 (27) 2 X 2 0 Cαkk0 fk0 = |fk| + d k .  fnk+g  ˆ (ω − ωαk) α ψn1k+g X ig·r   un(k, r) = unk+ge , unk+g = ψ  , The complete set of eigenvectors and eigenvalues is repre-  n2k+g g . sented by uk. We must define all relevant electromagnetic . quantities in terms of this generalized state vector. (31) where un(k, r + R) = un(k, r) is periodic in the crystal lattice. In this basis, |unki is normalized to the energy density E. Crystal Hamiltonian as,

X † We are now ready to enforce crystal periodicity. Instead of 1 = hunk|unki = unk+gunk+g expanding M directly, we utilize the periodicity of the cou- g 0 0 pling tensors Cα(r, r ) = Cα(r + R, r + R), which is a dis- ! X † X † crete spectrum in g, = fnk+gfnk+g + ψnαk+gψnαk+g (32) g α 0 X 2 0 Cα(k, k ) = Cαg(k)δ (k + g − k ). (28) X † ¯ 0 0 g = fnk+gMg −g(ωnk, k + g)fnk+g . gg0 Cαg(k) tells us the scattering amplitude of a photon fk+g with momentum k + g into an internal mode of the material ψαk The bra-ket notation h|i implies integration over the 2D unit at momentum k, and vice versa. The crystal Hamiltonian ac- cell and we have utilized the linear response theory to express counts for all such scattering events, ψα in terms of the driving field f, P 0 X 2 0 0 Cαg0 (k + g)fnk+g0+g H(k, k ) = Hg(k)δ (k + g − k ), (29) g ψnαk+g = . (33) g ωnk − ωαk+g † ¯ with Hermiticity Hg(k) = H−g(k + g) satisfied by defini- Mg(ω, k) = ∂ω [ωMg(ω, k)] is the contribution to the en- tion. Note, the resonant energies ωα(k + g) = ωα(k) are ergy density arising from each spatial harmonic of the crystal. 6

Finally, the eigenenergies ωnk are the n nontrivial roots of the characteristic wave equation,

X H0(k)fnk = ωnk Mg(ωnk, k)fnk+g, (34) g which generates all possible photonic bands of the crystal. Note, the response function Mg(ω, k) is now expressed in terms of Cαg(k) and describes the net summation of all scattering and back-scattering events in the material,

† 0 0 X Cα−g0 (k + g )Cαg−g0 (k + g ) Mg(ω, k) = 13δg − . ωαk+g0 (ω − ωαk+g0 ) αg0 (35) This proves that the wave equation is derived from a ﬁrst-order Hamiltonian, has real eigenvalues ω = ωnk for all momenta, and is normalizable in terms of |unki.

III. DISCRETE ROTATIONAL SYMMETRY

A. Point groups in 2D

Point groups are the discrete analogs of continuous rota- FIG. 1. Brillouin zone of each cyclic point group CN . (a), (b), tions and reflections. They represent the number of ways the (c), (d), and (e) correspond to N = 2, 3, 4, 6, and ∞ respectively. atomic lattice can be transformed into itself [92, 93]. Due to Due to rotational symmetry, the total Brillouin zone is equivalent the crystallographic restriction theorem (CRT), there are ten to N copies of the irreducible Brillouin zone (IBZ), which is rep- such point groups in 2D. The first five are the cyclic groups resented by the blue quadrant. For continuous symmetry N = ∞, CN , this is simply a line. The yellow circles label high-symmetry points Rki = ki where the crystal Hamiltonian is invariant under a certain rotation Rˆ. At these specific momenta, a Bloch photonic wave func- C1,C2,C3,C4,C6. (36) tion Rˆ N |f(ki)i = ηN (ki)|f(ki)i is an eigenstate of a an N-fold 2π rotation ηN (ki) = i N mN (ki) such that the photon possesses For instance, C3 implies threefold cyclic symmetry while C1 quantized integer eigenvalues mN (ki) ∈ ZN . Since mN are dis- is no symmetry. The last five are the dihedral groups DN , crete quantum numbers, their values cannot vary continuously if the crystal symmetry is preserved - they can only be changed at a topological phase transition. D1,D2,D3,D4,D6. (37)

The dihedral group DN contains CN plus reflections. How- ever, it can be proven that the Chern number for all DN point B. Spin-1 discrete symmetries groups vanish [69]. Therefore, we concern ourselves with C only the cyclic groups N . The Brillouin zone of each point If the two-dimensional crystal belongs to a cyclic point group is displayed in Fig. 1. group CN , the Hamiltonian possesses discrete rotational sym- The defining characteristic of each cyclic group is the metry about the z-axis, fermionic or bosonic representation. When we rotate the fields by 2π, we take the particle into itself and acquire a phase, −1 R HRg(Rk)R = Hg(k), ωn(Rk) = ωn(k), (39) R(2π) = (−1)F . (38) where R is any rotation in CN . It is important to note that F is twice the total spin of particle, or equivalently, the R is diagonal in u, meaning the photon and each oscillator is fermion number. Fermions with half-integer spin are anti- rotated individually, f → Rf and ψα → Rψα. This implies symmetric under rotations R(2π) = −1, while bosons with there is no mixing of fields. The symmetries of the Hamilto- integer spin are symmetric R(2π) = +1. Depending on the nian are endowed by the coupling tensors, which dictates the symmetries of the lattice, the topology fundamentally changes degrees of freedom of the material response, for fermions and bosons. We will understand the implications −1 this has for spin-1 photons. R CαRg(Rk)R = Cαg(k), ωα(Rk) = ωα(k). (40) 7

After summation over all Cαg(k), we can prove that the re- Importantly, there are speciﬁc points in the Brillouin zone sponse function transforms identically under such a rotation, where k is invariant under a discrete rotation,

−1 R MRg(ω, Rk)R = Mg(ω, k). (41) Rki = ki. (47)

Therefore, the photon inherits all symmetries of the crystal. This is because the crystal momentum only differs by a lattice In this case, the R matrix represents a discrete rotation and translation at these points Rki = ki +g, which leaves a Bloch can be expressed as the exponential of the spin-1 generator wave function unchanged, ˆ (Sz)ij = −iijz, iRki·r i(ki+g)·r e un(Rki, r) = e un(ki + g, r) (48)  2π 2π  iki·r cos N sin N 0 = e un(ki, r). 2π ˆ 2π 2π RN = exp i Sz = − sin N cos N 0 , (42) N 0 0 1 These are called high-symmetry points (HSPs); they occur at the center and certain vertices of the Brillouin zone. The crys- where 2π/N is an N-fold rotation. We stress that every cyclic tal Hamiltonian is rotationally invariant at these momenta - ˆ group for the photon is a vector representation, which is i.e. it commutes with R. Therefore, the wave functions are ˆ bosonic, simultaneous eigenstates of R at HSPs,

N R|ˆ u (k )i = η (k )|u (k )i, R(2π) = R(NθN ) = (RN ) = +13. (43) n i n i n i (49)

The electromagnetic ﬁeld returns in phase under cyclic revo- which immediately implies, lution. R|ˆ fn(ki)i = ηn(ki)|fn(ki)i. (50)

Here, ηn(ki) is the eigenvalue of Rˆ at ki for the nth band. C. High-symmetry points

The Bloch eigenstates |unki are essentially a collection of D. Spin-1 eigenvalues periodic vector fields. To rotate the fields, we must perform an operation on both the coordinates r and the polarization Depending on the point group and the precise HSP, ˆ states f and ψα. In real space, the operation of a rotation R is ηn(ki) = ηN,n(ki) can represent any Nth root of unity cor- preformed as, responding to the rotation operator RˆN , hr|R|ˆ u i = Ru (k, R−1r) = η (k)u (Rk, r), (44) 2π nk n n n η (k ) = exp i m (k ) , (η )N = +1. N,n i N N,n i N,n where R is a discrete rotation defined in Eq. (42). This implies (51) the Fourier coefficients obey, mN,n(kp) ∈ ZN is a modulo integer - it labels the N possible spin-1 eigenvalues at ki. In C4 for example, the Γ and Runk+R−1g = ηn(k)unRk+g. (45) M points are invariant under Rˆ4 rotations, while the X and ˆ It follows from symmetry that the operation of Rˆ takes a wave Y points are invariant under R2 rotations (inversion). This means there are 4 possible spin-1 charges located at m4,n(Γ) function at k to Rk with the same energy ωn(k) = ωn(Rk) - 2 & m4,n(M) ∈ 4 respectively and 2 possible charges located but with a possibly different phase |ηn(k)| = 1. Utilizing the Z at m2,n(X) = m2,n(Y) ∈ 2. A visualization of these topo- linear response theory, we notice that the phase factor ηn(k) Z is governed entirely by the photon, logical charges is presented in Fig. 2 and this is contrasted with their fermionic counterparts in Fig. 3. In Sec. IV we will −1 P 0 connect these rotational eigenvalues directly to the topological g0 RCαg (k + R g)fnk+g0+R−1g Rψnαk+R−1g = invariants. ωnk − ωαk+R−1g P 0 CαRg0 (Rk + g)Rfnk+g0+R−1g = g IV. TOPOLOGICAL ELECTROMAGNETIC (BOSONIC) ωnk − ωαk+R−1g P PHASES OF MATTER 0 Cαg0 (Rk + g)ηn(k)fnRk+g0+g = g ωnk − ωαRk+g A. Electromagnetic Chern number = ηn(k)ψnαRk+g. (46) The Berry connection for a band n is found by varying the total Bloch wave function |unki with respect to the momen- This is an incredibly convenient simplification and implies the tum, precise coordinates of the matter oscillations ψα are superflu- ous when discussing symmetries. The electromagnetic field f X † An(k) = −ihunk|∂kunki = −i unk+g∂kunk+g. (52) tells us everything. g 8

4 FIG. 2. The collection of spin-1 (bosonic) charges for the C4 point group. (a) Fourfold rotations (R4) = +1; there are four unique 2π 4 eigenvalues η4 = exp i 4 m4 corresponding to the roots of unity (η4) = +1. These represent the modulo 4 integers m4 ∈ Z4. Note that 2 m4 = 3 = −1 can also be interpreted as a left-handed eigenstate. (b) Bosonic inversion (R2) = +1; there are two unique eigenvalues 2π 2 η2 = exp i 2 m2 corresponding to the roots of unity (η2) = +1. These represent the modulo 2 integers m2 ∈ Z2.

1 4 FIG. 3. The collection of spin- ⁄2 (fermionic) charges for the C4 point group. (a) Fourfold rotations (R4) = −1; there are four unique 2π 4 eigenvalues ζ4 = exp i 4 m4 corresponding to the roots of negative unity (ζ4) = −1. These represent the modulo 4 half-integers 1 7 1 3 1 5 1 m4 ∈ Z4 + /2. Note that m4 = /2 = − /2 can be interpreted as a spin-down fermion while m4 = /2 = /2 + 1 and m4 = /2 = − /2 + 3 2 2π constitute a fermion plus a boson. (b) Fermionic inversion (R2) = −1; there are two unique eigenvalues ζ2 = exp i 2 m2 corresponding 2 1 3 1 to the roots of negative unity (ζ2) = −1. These represent the modulo 2 half-integers m2 ∈ Z2 + /2. Note that m2 = /2 = − /2 can also be interpreted as a spin-down fermion under modulo 2.

This can be simplified slightly to obtain, As can be seen from Eq. (53), the Berry connection is only defined within the Brillouin zone Ank+g = Ank + ∂kχnk, X † ¯ An(k) = − i fnk+gMg0−g(ωnk, k + g)∂kfnk+g0 up to a possible U(1) gauge. Hence, the gauge invariant Berry 0 gg curvature is periodic Fnk+g = Fnk, X † + fnk+gAg0−g(ωnk, k + g)fnk+g0 . Fn(k) = ˆz · [∂k × An(k)]. (55) gg0 (53) The Chern number is found by integrating the Berry curvature The first term gives the Berry connection of the photon, while over the two-dimensional Brillouin zone, the second term Ag(ω, k) arises solely from the matter oscil- 1 2 lations, Cn = Fn(k)d k, Cn ∈ Z, (56) 2π ˆBZ † 0 0 X Cα−g0 (k + g )∂kCαg−g0 (k + g ) A which determines the winding number of the collective light- Ag(ω, k) = −i 2 . (ω − ωαk+g0 ) 2 1 1 αg0 matter excitations over the torus T = S ×S . Equation (56) (54) is one of the central results of this paper. An electromagnetic Due to nonlocality, Eq. (54) does not generally vanish. This Chern invariant can be found for any 2D crystal and charac- additional contribution to the Berry phase corresponds to vor- terizes distinct topological phases of matter Cn 6= 0. tices in the response function itself - independent of the Berry gauge of the photon. This means the Chern number can be B. Symmetry-protected topological bosonic phases nonzero Cn 6= 0 even if the winding of electromagnetic field is trivial. However, we will show in the proceeding sections that all symmetry constraints on the Chern number can be es- Nevertheless, even if we knew the specifics of the mate- tablished entirely in terms of the photon. rial, evaluating the Chern number by brute force would be a 9

FIG. 4. Examples of SPT bosonic phases in a crystal with C4 symmetry. These phases are characterized by their SPT invariant ν = m4(Γ) + m4(M) + 2m2(Y) mod 4 which determines the electromagnetic Chern number up to a multiple of 4. Here, m4 ∈ Z4 and m2 ∈ Z2 are modulo integers. (a), (b), (c) and (d) correspond to SPT bosonic phases of ν = 3, 2, 1 and 0 respectively. For bosons, we simply add up all the integer charges within the irreducible Brillouin zone. For instance, the ν = 2 phase has eigenvalues of m4(Γ) = 1 at the center and m4(M) = 3 = −1 at the vertices, with inversion eigenvalues of m2(Y) = m2(X) = 1 at the edge centers: ν = 1 + 3 + 2 × 1 = 2.

FIG. 5. Examples of SPT fermionic phases in a crystal with C4 symmetry. These phases are characterized by their SPT invariant ν = 1 m4(Γ) + m4(M) + 2m2(Y) + 2 mod 4 which determines the electronic Chern number up to a multiple of 4. In this case, m4 ∈ Z4 + /2 1 and m2 ∈ Z2 + /2 are modulo half-integers. (a), (b), (c) and (d) correspond to SPT fermionic phases of ν = 3, 2, 1 and 0 respectively. The problem is more complicated for fermions because the charges are fractional and we must also account for the antisymmetric phases of a 1 spinor wave function. As an example, the ν = 2 phase has eigenvalues of m4(Γ) = m4(M) = /2 at the center and vertices, with inversion 3 1 1 1 3 eigenvalues of m2(Y) = m2(X) = /2 = − /2 at the edge centers: ν = /2 + /2 + 2 × /2 + 2 = 2. herculean task. Instead, we invoke constraints of the point contributes a Berry phase of 2π, and by symmetry, there are groups, which constitute a type of symmetry-protected topo- N such vortices within the total Brillouin zone Cn → Cn +N. logical (SPT) phase [66–77]. SPT phases are protected by the However, this has no effect on νn → νn. Only the vortices N-fold rotational symmetry of CN and this gives rise to an lying at HSPs contribute to νn because these come in fractions additional topological invariant νn ∈ ZN . Remarkably, νn of 2π. is classified entirely from ηn(ki) eigenvalues at HSPs and re- In the following sections we will discuss the bosonic classi- quires no complicated integration to compute. This invariant fication of νn for each cyclic point group and the SPT phases is related to the Chern number up to a multiple of N, associated with them. We do not present the full derivations here since the rigorous proofs have been carried out by others ν = C mod N, C ∈ N + ν . n n n Z n (57) (see Ref. [69]) - we simply state the salient results. For completeness, in Appendix A we also discuss the SPT fermionic The interpretation of νn is quite simple - it tells us the geometric phase around the irreducible Brillouin zone (IBZ) of phases associated with each point group. We do this to em- the crystal, phasize that fermionic and bosonic systems represent distinct topological field theories, with fundamentally different inter- pretations. These differences are highlighted with a few ex- 2π 2 exp i Cn = exp i Fn(k)d k amples [Fig. 4 and 5]. N ÎBZ (58) = exp i An(k) · dk , ˛∂IBZ where ∂IBZ is the path around IBZ. This follows from ro- C. Twofold (inversion) symmetry: C2 tational symmetry of the Berry curvature Fn(k) = Fn(Rk). For instance, the path in C4 is ∂(IBZ)4 = ΓXMYΓ. Applying For the C2 point group, or simply inversion symmetry, the the logarithm, νn is equivalent to, SPT phase is related to the Chern number by νn = Cn mod 2 which is a invariant. There is only one nontrivial SPT N Z2 νn = An(k) · dk mod N. (59) phase and it can be found modulo 2 from, 2π ˛∂RBZ As we will see more explicitly, ν is tied entirely to η . The 2π n n exp i C = η (Γ)η (X)η (Y)η (M). (60) reason is subtle - any vortex within the interior of the IBZ 2 n 2,n 2,n 2,n 2,n 10

Applying the logarithm, this classiﬁcation can be expressed This is equivalent to the summation of spin-1 eigenvalues at equivalently in terms of m2,n ∈ Z2 inversion eigenvalues, the HSPs, ν = m (Γ) + m (X) + m (Y) + m (M) mod 2. n 2,n 2,n 2,n 2,n νn = m6,n(Γ) + 2m3,n(K) + 3m2,n(M) mod 6, (67) (61) If the summation of m eigenvalues is odd, the SPT phase 2,n where m6,n(Γ) ∈ Z6 is a modulo 4 integer, m6,n(K) ∈ Z3 is nontrivial ν = 1 and corresponds to an odd-valued Chern n is a modulo 3 integer and m2,n(M) ∈ Z2 is a modulo 2 inte- number. Likewise, νn = 0 is an even-valued Chern number. ger. This completes the classiﬁcation of all 2+1D topological electromagnetic (bosonic) phases of matter which is summarized in Tbl. I. These are compared alongside their fermionic D. Threefold symmetry: C3 counterparts in Tbl. II.

C3 is unique because it is the only point group with an odd rotational symmetry - i.e. it lacks inversion symmetry. This G. Continuous symmetry: C∞ means the parity of Chern number (odd or even) is not re- stricted by the symmetries of the crystal. For C3, the SPT To finish, we briefly discuss the continuum limit g = 0 phase is νn = Cn mod 3 which is a Z3 invariant. There are two nontrivial SPT phases and they can be found modulo 3 and the topological phases that can be described by a long from, wavelength theory k ≈ 0. The physics is significantly more tractable here and exactly solvable models are possible 2π 0 [24, 25]. In this limit, the rotational symmetry of the crystal exp i Cn = η3,n (Γ) η3,n (K) η3,n K . (62) 3 is approximately continuous C∞. The SPT invariant νn and Chern number Cn are thus equivalent, This classification is expressed equivalently in terms of quantized modulo 3 integers m ∈ at HSPs, 3,n Z3 νn = Cn = mn(0) − mn(∞). (68) 0 νn = m3,n(Γ) + m3,n(K) + m3,n(K ) mod 3. (63) Note that νn ∈ Z and mn ∈ Z are not modulo integers in Note though, odd and even phases are not distinct ν = −2 = this limit and do not have the same interpretation as the lat- 1 = 4 under modulo 3. tice theory. This is because we have gained the full rotational symmetry in the continuum approximation. Clearly though, the eigenvalues must change at HSPs mn(0) 6= mn(∞) for a E. Fourfold symmetry: C4 nontrivial phase to exist Cn 6= 0. In the continuum regulariza- tion, ki = 0 represents the Γ point and ki = ∞ is interpreted as mapping the vertices of the Brillouin zone into one another. For the C4 point group, the SPT phase is related to the Chern number by νn = Cn mod 4 which is a Z4 invariant. There are three nontrivial SPT phases and they can be found modulo 4 from, V. CONCLUSIONS 2π exp i Cn = η4,n (Γ) η4,n (M) η2,n (Y) . (64) In summary, we have developed the complete 2+1D lat- 4 tice field theory describing all symmetry-protected topolog- The classification is expressed equivalently in terms of spin-1 ical bosonic phases of the photon. To accomplish this, we an- eigenvalues, alyzed the electromagnetic Bloch waves in microscopic crystals and derived the Chern invariant of these light-matter ex- νn = m4,n(Γ) + m4,n(M) + 2m2,n(Y) mod 4, (65) citations. Thereafter, the rotational symmetries of the crystal were examined extensively and the implications these have where m4,n(Γ) & m4,n(M) ∈ Z4 are modulo 4 integers and on photonic spin. We have studied all two dimensional point m2,n(Y) ∈ 2 is a modulo 2 integer. Examples of all SPT Z groups CN with nonvanishing Chern number C 6= 0 and phases of the C4 point group are displayed in Fig. 4 and these linked the topological invariants directly to spin-1 quantized are compared with their fermionic counterparts in Fig. 5. eigenvalues of the electromagnetic field - establishing the bosonic classification for each topological phase.

F. Sixfold symmetry: C6

ACKNOWLEDGEMENTS For the C6 point group, the SPT phase is νn = Cn mod 6 which is a Z6 invariant. There are ﬁve nontrivial SPT phases and they can be found modulo 6 from, This research was supported by the Defense Advanced Re- search Projects Agency (DARPA) Nascent Light-Matter Inter- 2π actions (NLM) Program and the National Science Foundation exp i C = η (Γ) η (K) η (M) . (66) 6 n 6,n 3,n 2,n (NSF) [Grant No. EFMA-1641101]. 11

2π APPENDIX exp i C = −ζ (Γ) ζ (K) ζ K0 , (A3b) 3 3 3 3 Appendix A: Symmetry-protected topological fermionic phases

For completeness, we examine the SPT fermionic phases associated with each point group CN and highlight their es- 2π exp i C = −ζ (Γ) ζ (M) ζ (Y) , (A3c) sential differences from bosons. The most important distinc- 4 4 4 2 tion is how they transform under rotations; half-integer particles are antisymmetric R(2π) = −1. In terms of discrete rotations RˆN about the z-axis, the eigenstates of a Bloch spinor particle satisfy, 2π exp i C = −ζ6 (Γ) ζ3 (K) ζ2 (M) . (A3d) RˆN |Ψ(ki)i = ζN (ki)|Ψ(ki)i, (A1) 6 where the eigenvalues at HSPs are related by, Although the classification appears similar, the SPT fermionic 2π phases constitute very different physics than their bosonic ζ (k ) = exp i m (k ) , (ζ )N = −1. (A2) N i N N i N counterparts, which is alluded to by the antisymmetric phase factors R(2π) = −1. We illustrate this with an example in 1 mN (ki) ∈ ZN + /2 is a modulo half-integer and labels the C4. Applying the logarithm - the classification for the SPT 1 N possible spin- ⁄2 eigenvalues. Notice that ζN represents fermionic phase ν = C mod 4 can be expressed as, the Nth roots of negative unity which is characteristic of a fermionic field. ν = m4(Γ) + m4(M) + 2m2(Y) + 2 mod 4, (A4) The single-particle fermionic classification for C2, C3, C4 and C6 respectively is [69, 76], 1 2π where m4(Γ) & m4(M) ∈ Z4 + /2 are modulo 4 half-integers exp i C = ζ (Γ)ζ (X)ζ (Y)ζ (M), (A3a) and m (Y) ∈ + 1/2 is a modulo 2 half-integer. 2 2 2 2 2 2 Z2

[1] M. Z. Hasan and C. L. Kane, “Colloquium: Topological insula- based on quantized hall resistance,” Phys. Rev. Lett. 45, 494– tors,” Rev. Mod. Phys. 82, 3045–3067 (2010). 497 (1980). [2] Shinsei Ryu, Andreas P Schnyder, Akira Furusaki, and An- [11] R. B. Laughlin, “Quantized hall conductivity in two dimen- dreas W W Ludwig, “Topological insulators and superconduc- sions,” Phys. Rev. B 23, 5632–5633 (1981). tors: tenfold way and dimensional hierarchy,” New Journal of [12] D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. den Physics 12, 065010 (2010). Nijs, “Quantized hall conductance in a two-dimensional peri- [3] F. D. M. Haldane, “Model for a quantum hall effect with- odic potential,” Phys. Rev. Lett. 49, 405–408 (1982). out landau levels: Condensed-matter realization of the ”parity [13] Xie Chen, Zheng-Xin Liu, and Xiao-Gang Wen, “Two- anomaly”,” Phys. Rev. Lett. 61, 2015–2018 (1988). dimensional symmetry-protected topological orders and their [4] D Jaksch and P Zoller, “Creation of effective magnetic fields in protected gapless edge excitations,” Phys. Rev. B 84, 235141 optical lattices: the hofstadter butterfly for cold neutral atoms,” (2011). New Journal of Physics 5, 56 (2003). [14] Xie Chen, Zheng-Cheng Gu, Zheng-Xin Liu, and Xiao-Gang [5] Jens Koch, Andrew A. Houck, Karyn Le Hur, and S. M. Girvin, Wen, “Symmetry-protected topological orders in interacting “Time-reversal-symmetry breaking in circuit-qed-based photon bosonic systems,” Science 338, 1604–1606 (2012). lattices,” Phys. Rev. A 82, 043811 (2010). [15] Yuan-Ming Lu and Ashvin Vishwanath, “Theory and classifica- [6] Gregor Jotzu, Michael Messer, Remi´ Desbuquois, Martin Le- tion of interacting integer topological phases in two dimensions: brat, Thomas Uehlinger, Daniel Greif, and Tilman Esslinger, A chern-simons approach,” Phys. Rev. B 86, 125119 (2012). “Experimental realization of the topological haldane model [16] Ashvin Vishwanath and T. Senthil, “Physics of three- with ultracold fermions,” Nature 515, 237 EP – (2014). dimensional bosonic topological insulators: Surface- [7] C. L. Kane and E. J. Mele, “Quantum spin hall effect in deconfined criticality and quantized magnetoelectric effect,” graphene,” Phys. Rev. Lett. 95, 226801 (2005). Phys. Rev. X 3, 011016 (2013). [8] B. Andrei Bernevig, Taylor L. Hughes, and Shou-Cheng [17] Max A. Metlitski, C. L. Kane, and Matthew P. A. Fisher, Zhang, “Quantum spin hall effect and topological phase tran- “Bosonic topological insulator in three dimensions and the sta- sition in hgte quantum wells,” Science 314, 1757–1761 (2006). tistical witten effect,” Phys. Rev. B 88, 035131 (2013). [9] Joseph Maciejko, Taylor L. Hughes, and Shou-Cheng Zhang, [18] T. Senthil and Michael Levin, “Integer quantum hall effect for “The quantum spin hall effect,” Annual Review of Condensed bosons,” Phys. Rev. Lett. 110, 046801 (2013). Matter Physics 2, 31–53 (2011). [19] Zeren Lin and Zhirong Liu, “Spin-1 dirac-weyl fermions pro- [10] K. v. Klitzing, G. Dorda, and M. Pepper, “New method tected by bipartite symmetry,” The Journal of Chemical Physics for high-accuracy determination of the fine-structure constant 143, 214109 (2015). 12

[20] Tian Lan, Liang Kong, and Xiao-Gang Wen, “The- lidi (RRL) Rapid Research Letters 7, 101–108 (2013). ory of (2+1)-dimensional fermionic topological orders and [41] Mikael C. Rechtsman, Julia M. Zeuner, Yonatan Plotnik, fermionic/bosonic topological orders with symmetries,” Phys. Yaakov Lumer, Daniel Podolsky, Felix Dreisow, Stefan Nolte, Rev. B 94, 155113 (2016). Mordechai Segev, and Alexander Szameit, “Photonic floquet [21] Tiantian Zhang, Zhida Song, A. Alexandradinata, Hongming topological insulators,” Nature 496, 196–200 (2013), letter. Weng, Chen Fang, Ling Lu, and Zhong Fang, “Double-weyl [42] Alexander B. Khanikaev and Gennady Shvets, “Two- phonons in transition-metal monosilicides,” Phys. Rev. Lett. dimensional topological photonics,” Nature Photonics 11, 763– 120, 016401 (2018). 773 (2017). [22] Y. Ikebe, T. Morimoto, R. Masutomi, T. Okamoto, H. Aoki, [43] Alexey Slobozhanyuk, S. Hossein Mousavi, Xiang Ni, Daria and R. Shimano, “Optical hall effect in the integer quantum hall Smirnova, Yuri S. Kivshar, and Alexander B. Khanikaev, regime,” Phys. Rev. Lett. 104, 256802 (2010). “Three-dimensional all-dielectric photonic topological insula- [23] Yisong Zheng and Tsuneya Ando, “Hall conductivity of a two- tor,” Nat Photon 11, 130–136 (2017), article. dimensional graphite system,” Phys. Rev. B 65, 245420 (2002). [44] Cheng He, Xiao-Chen Sun, Xiao-Ping Liu, Ming-Hui Lu, Yulin [24] Todd Van Mechelen and Zubin Jacob, “Quantum gyroelectric Chen, Liang Feng, and Yan-Feng Chen, “Photonic topological effect: Photon spin-1 quantization in continuum topological insulator with broken time-reversal symmetry,” Proceedings of bosonic phases,” Phys. Rev. A 98, 023842 (2018). the National Academy of Sciences 113, 4924–4928 (2016). [25] Todd Van Mechelen and Zubin Jacob, “Photonic [45] Stanislav B. Glybovski, Sergei A. Tretyakov, Pavel A. Belov, Dirac monopoles and skyrmions: spin-1 quantization,” Yuri S. Kivshar, and Constantin R. Simovski, “Metasurfaces: arXiv:1806.09879 (2018). From microwaves to visible,” Physics Reports 634, 1 – 72 [26] G. V. Dunne, “Aspects of chern-simons theory,” in Aspects (2016), metasurfaces: From microwaves to visible. topologiques de la physique en basse dimension. Topologi- [46] I. V. Lindell and A. H. Sihvola, “Realization of the pemc cal aspects of low dimensional systems, edited by A. Comtet, boundary,” IEEE Transactions on Antennas and Propagation 53, T. Jolicœur, S. Ouvry, and F. David (Springer Berlin Heidel- 3012–3018 (2005). berg, Berlin, Heidelberg, 1999) pp. 177–263. [47] Rasoul Alaee, Mohamed Farhat, Carsten Rockstuhl, and Falk [27] D. Boyanovsky, R. Blankenbecler, and R. Yahalom, “Physi- Lederer, “A perfect absorber made of a graphene micro-ribbon cal origin of topological mass in 2 + 1 dimensions,” Nuclear metamaterial,” Opt. Express 20, 28017–28024 (2012). Physics B 270, 483 – 505 (1986). [48] Nikitas Papasimakis, Zhiqiang Luo, Ze Xiang Shen, [28] S. A. R. Horsley, “Topology and the optical dirac equation,” Francesco De Angelis, Enzo Di Fabrizio, Andrey E. Niko- Phys. Rev. A 98, 043837 (2018). laenko, and Nikolay I. Zheludev, “Graphene in a photonic [29] S. Raghu and F. D. M. Haldane, “Analogs of quantum-hall- metamaterial,” Opt. Express 18, 8353–8359 (2010). effect edge states in photonic crystals,” Phys. Rev. A 78, 033834 [49] Junfei Li, Chen Shen, Ana D´ıaz-Rubio, Sergei A. Tretyakov, (2008). and Steven A. Cummer, “Systematic design and experimental [30] F. D. M. Haldane and S. Raghu, “Possible realization of di- demonstration of bianisotropic metasurfaces for scattering-free rectional optical waveguides in photonic crystals with broken manipulation of acoustic wavefronts,” Nature Communications time-reversal symmetry,” Phys. Rev. Lett. 100, 013904 (2008). 9, 1342 (2018). [31] Zheng Wang, Yidong Chong, J. D. Joannopoulos, and [50] Yu Guo, Meng Xiao, and Shanhui Fan, “Topologically pro- Marin Soljacic, “Observation of unidirectional backscattering- tected complete polarization conversion,” Phys. Rev. Lett. 119, immune topological electromagnetic states,” Nature 461, 772– 167401 (2017). 775 (2009). [51] Kun Ding, Guancong Ma, Meng Xiao, Z. Q. Zhang, and C. T. [32] Zheng Wang and Shanhui Fan, “Optical circulators in two- Chan, “Emergence, coalescence, and topological properties of dimensional magneto-optical photonic crystals,” Opt. Lett. 30, multiple exceptional points and their experimental realization,” 1989–1991 (2005). Phys. Rev. X 6, 021007 (2016). [33] Ling Lu, Liang Fu, John D. Joannopoulos, and Marin Soljacic, [52] Mario´ G. Silveirinha, “Chern invariants for continuous media,” “Weyl points and line nodes in gyroid photonic crystals,” Nat Phys. Rev. B 92, 125153 (2015). Photon 7, 294–299 (2013). [53] Wenlong Gao, Mark Lawrence, Biao Yang, Fu Liu, Fengzhou [34] M. Hafezi, S. Mittal, J. Fan, A. Migdall, and J. M. Taylor, Fang, Benjamin Beri,´ Jensen Li, and Shuang Zhang, “Topolog- “Imaging topological edge states in silicon photonics,” Nat Pho- ical photonic phase in chiral hyperbolic metamaterials,” Phys. ton 7, 1001–1005 (2013), article. Rev. Lett. 114, 037402 (2015). [35] Zheng-Cheng Gu, Zhenghan Wang, and Xiao-Gang Wen, [54] Dafei Jin, Ling Lu, Zhong Wang, Chen Fang, John D. “Classification of two-dimensional fermionic and bosonic topo- Joannopoulos, Marin Soljacic, Liang Fu, and Nicholas X. logical orders,” Phys. Rev. B 91, 125149 (2015). Fang, “Topological magnetoplasmon,” Nature Communications [36] Torsten Karzig, Charles-Edouard Bardyn, Netanel H. Lind- 7, 13486 EP – (2016), article. ner, and Gil Refael, “Topological polaritons,” Phys. Rev. X 5, [55] S. A. Hassani Gangaraj, M. G. Silveirinha, and G. W. Hanson, 031001 (2015). “Berry phase, berry connection, and chern number for a con- [37] Yakir Hadad, Vincenzo Vitelli, and Andrea Alu, “Solitons and tinuum bianisotropic material from a classical electromagnet- propagating domain walls in topological resonator arrays,” ACS ics perspective,” IEEE Journal on Multiscale and Multiphysics Photonics 4, 1974–1979 (2017). Computational Techniques 2, 3–17 (2017). [38] Giuseppe De Nittis and Max Lein, “Symmetry Classification of [56] Li-kun Shi and Justin C. W. Song, “Plasmon geometric phase Topological Photonic Crystals,” arXiv:1710.08104 (2017). and plasmon hall shift,” Phys. Rev. X 8, 021020 (2018). [39] Netanel H. Lindner, Gil Refael, and Victor Galitski, “Floquet [57] Iwo Bialynicki-Birula and Zofia Bialynicka-Birula, “Berry’s topological insulator in semiconductor quantum wells,” Nature phase in the relativistic theory of spinning particles,” Phys. Rev. Physics 7, 490 EP – (2011), article. D 35, 2383–2387 (1987). [40] Jerome Cayssol, Balazs Dora, Ferenc Simon, and Roderich [58] Michael Stone, “Berry phase and anomalous velocity of weyl Moessner, “Floquet topological insulators,” physica status so- fermions and maxwell photons,” International Journal of Mod- 13

ern Physics B 30, 1550249 (2016). netic space groups,” Science Advances 4 (2018). [59] Michael Stone, “Topology, spin, and light,” Science 348, 1432– [78] Barry Bradlyn, Jennifer Cano, Zhijun Wang, M. G. Vergniory, 1433 (2015). C. Felser, R. J. Cava, and B. Andrei Bernevig, “Beyond dirac [60] O. El Gawhary, T. Van Mechelen, and H. P. Urbach, “Role and weyl fermions: Unconventional quasiparticles in conven- of radial charges on the angular momentum of electromagnetic tional crystals,” Science 353 (2016). fields: Spin-3/2 light,” Phys. Rev. Lett. 121, 123202 (2018). [79] Yan-Qing Zhu, Dan-Wei Zhang, Hui Yan, Ding-Yu Xing, and [61] J Hubbard, “The dielectric theory of electronic interactions in Shi-Liang Zhu, “Emergent pseudospin-1 maxwell fermions solids,” Proceedings of the Physical Society. Section A 68, 976 with a threefold degeneracy in optical lattices,” Phys. Rev. A (1955). 96, 033634 (2017). [62] David S. Falk, “Effect of the lattice on dielectric properties of [80] Haiping Hu and Chuanwei Zhang, “Spin-1 topological an electron gas,” Phys. Rev. 118, 105–109 (1960). monopoles in the parameter space of ultracold atoms,” Phys. [63] David R. Penn, “Wave-number-dependent dielectric function of Rev. A 98, 013627 (2018). semiconductors,” Phys. Rev. 128, 2093–2097 (1962). [81] Haiping Hu, Junpeng Hou, Fan Zhang, and Chuanwei Zhang, [64] Lev Davidovich Landau, JS Bell, MJ Kearsley, LP Pitaevskii, “Topological triply degenerate points induced by spin-tensor- EM Lifshitz, and JB Sykes, Electrodynamics of Continuous momentum couplings,” Phys. Rev. Lett. 120, 240401 (2018). Media, Vol. 8 (Elsevier, New York, 2013). [82] I. C. Fulga, L. Fallani, and M. Burrello, “Geometrically pro- [65] Vladimir M Agranovich and Vitaly Ginzburg, Crystal optics tected triple-point crossings in an optical lattice,” Phys. Rev. B with spatial dispersion, and excitons, Vol. 42 (Springer Science 97, 121402 (2018). & Business Media, New York, 2013). [83] Xinsheng Tan, Dan-Wei Zhang, Qiang Liu, Guangming Xue, [66] Liang Fu, “Topological crystalline insulators,” Phys. Rev. Lett. Hai-Feng Yu, Yan-Qing Zhu, Hui Yan, Shi-Liang Zhu, and 106, 106802 (2011). Yang Yu, “Topological maxwell metal bands in a superconduct- [67] Taylor L. Hughes, Emil Prodan, and B. Andrei Bernevig, ing qutrit,” Phys. Rev. Lett. 120, 130503 (2018). “Inversion-symmetric topological insulators,” Phys. Rev. B 83, [84] Todd Van Mechelen and Zubin Jacob, “Dirac-Maxwell cor- 245132 (2011). respondence: Spin-1 bosonic topological insulator,” arXiv: [68] Robert-Jan Slager, Andrej Mesaros, Vladimir Juricic, and Jan 1708.08192 (2017). Zaanen, “The space group classification of topological band- [85] Xiaobo Yin, Ziliang Ye, Junsuk Rho, Yuan Wang, and Xiang insulators,” Nature Physics 9, 98 EP – (2012), article. Zhang, “Photonic spin hall effect at metasurfaces,” Science 339, [69] Chen Fang, Matthew J. Gilbert, and B. Andrei Bernevig, “Bulk 1405–1407 (2013). topological invariants in noninteracting point group symmetric [86] Todd Van Mechelen and Zubin Jacob, “Universal spin- insulators,” Phys. Rev. B 86, 115112 (2012). momentum locking of evanescent waves,” Optica 3, 118–126 [70] Jorrit Kruthoff, Jan de Boer, Jasper van Wezel, Charles L. Kane, (2016). and Robert-Jan Slager, “Topological classification of crystalline [87] Farid Kalhor, Thomas Thundat, and Zubin Jacob, “Universal insulators through band structure combinatorics,” Phys. Rev. X spin-momentum locked optical forces,” Applied Physics Letters 7, 041069 (2017). 108, 061102 (2016). [71] Barry Bradlyn, L. Elcoro, Jennifer Cano, M. G. Vergniory, Zhi- [88] Sarang Pendharker, Farid Kalhor, Todd Van Mechelen, Saman jun Wang, C. Felser, M. I. Aroyo, and B. Andrei Bernevig, Jahani, Neda Nazemifard, Thomas Thundat, and Zubin Jacob, “Topological quantum chemistry,” Nature 547, 298 EP – “Spin photonic forces in non-reciprocal waveguides,” Opt. Ex- (2017), article. press 26, 23898–23910 (2018). [72] Hao Song, Sheng-Jie Huang, Liang Fu, and Michael Hermele, [89] Konstantin Y. Bliokh, Daria Smirnova, and Franco Nori, “Topological phases protected by point group symmetry,” Phys. “Quantum spin hall effect of light,” Science 348, 1448–1451 Rev. X 7, 011020 (2017). (2015). [73] Tigran A. Sedrakyan, Victor M. Galitski, and Alex Kamenev, [90] Max Born and Kun Huang, Dynamical theory of crystal lattices “Topological spin ordering via chern-simons superconductiv- (Clarendon Press, Oxford, 1954). ity,” Phys. Rev. B 95, 094511 (2017). [91] S A R Horsley and T G Philbin, “Canonical quantization of [74] Hoi Chun Po, Ashvin Vishwanath, and Haruki Watanabe, electromagnetism in spatially dispersive media,” New Journal “Symmetry-based indicators of band topology in the 230 space of Physics 16, 013030 (2014). groups,” Nature Communications 8, 50 (2017). [92] David Hestenes, “Point groups and space groups in geometric [75] Ryan Thorngren and Dominic V. Else, “Gauging spatial sym- algebra,” in Applications of Geometric Algebra in Computer metries and the classification of topological crystalline phases,” Science and Engineering, edited by Leo Dorst, Chris Doran, Phys. Rev. X 8, 011040 (2018). and Joan Lasenby (Birkhauser¨ Boston, Boston, MA, 2002) pp. [76] Akishi Matsugatani, Yuri Ishiguro, Ken Shiozaki, and Haruki 3–34. Watanabe, “Universal relation among the many-body chern [93] John Frederick Nye, Physical properties of crystals: their rep- number, rotation symmetry, and filling,” Phys. Rev. Lett. 120, resentation by tensors and matrices (Oxford University Press, 096601 (2018). New York, 1985). [77] Haruki Watanabe, Hoi Chun Po, and Ashvin Vishwanath, “Structure and topology of band structures in the 1651 mag-