arxiv: v2 [cond-mat.str-el] 6 Sep 2016

Transcription

1 Abelin Floquet symmetry-protected topologicl phses in one dimension Rhul Roy nd Fenner Hrper Deprtment of Physics nd Astronomy, University of Cliforni, Los Angeles, Cliforni 90095, USA Dted: September 11, 018) rxiv: v [cond-mt.str-el] 6 Sep 016 Time-dependent systems hve recently been shown to support novel types of topologicl order tht cnnot be relized in sttic systems. In this pper, we consider rnge of time-dependent, intercting systems in one dimension tht re protected by n Abelin symmetry group. We clssify the distinct topologicl phses tht cn exist in this setting nd find tht they my be described by bulk invrint ssocited with the unitry evolution of the closed system. In the open system, nontrivil phses correspond to the ppernce of edge modes, which hve signtures in the mnybody qusienergy spectrum nd which relte to the bulk invrint through form of bulk-edge correspondence. We introduce simple models which relize nontrivil dynmicl phses in number of cses, nd outline loop construction tht cn be used to generte such phses more generlly. I. INTRODUCTION In recent yers, topologicl order hs become concept of fundmentl importnce in condensed mtter physics. The most unifying work in this field hs been the production of detiled clssifiction schemes tht systemticlly ctegorize the vriety of topologicl orders tht my exist in physicl systems under certin conditions. The originl clssifiction schemes considered gpped, nonintercting fermionic systems, leding to the so-clled tenfold wy of topologicl insultors nd superconductors [1, ]. More recently, clssifiction schemes hve been developed for intercting systems tht re protected by globl symmetry, G, but which hve no intrinsic longrnge) topologicl order, known s symmetry protected topologicl phses SPTs) [3 0]. Using group cohomology nd other techniques, schemes hve been developed tht clssify SPT phses for systems with mny different symmetry groups nd in vrious dimensions. Topologicl or SPT order ws originlly thought only to exist in strict sense t zero temperture, since mny topologicl phses brek down in the presence of therml excittions [1, ]. Recently, however, this belief hs been overturned by the theoreticl prediction of mnybody loclistion MBL), phenomenon wherein strong disorder induces non-therml behvior in mny-body system, which my in turn exhibit ordered phses tht correspond to high tempertures [3 3]. In systems undergoing MBL, topologicl order cn persist bove zero temperture [9, 33]. Moreover, in some MBL SPT systems with sufficiently strong disorder, the entire spectrum cn disply signtures of topologicl order, extending the notion of clssifiction beyond the study of the ground stte [34 37]. SPT order hs lso been studied in the context of Floquet systems: systems described by time-dependent Hmiltonin Ht) tht vries periodiclly in time with period T. In Floquet system, to be described more fully below, the relevnt sttes re the eigensttes of the unitry evolution opertor t the end of the driving cycle, UT ) [38]. Using periodic driving, it is possible to force system into topologicl phse see Ref. [39] for review), nd in prticulr, to generte stroboscopic SPT phses [40]. Remrkbly, Floquet systems hve lso been shown to exhibit novel types of topologicl order tht cnnot exist in sttic systems, including the existence of edge modes when the bulk Hmiltonin is trivil [41 50]. More recently, topologicl order in Floquet systems hs been shown to persist in the presence of disorder [51, 5]. The study of intercting phses, however, hs proved to be more difficult endevour. Driven, intercting systems genericlly led to heting, nd the resulting eigensttes might nively be expected to fll into trivil, infinite-temperture, ergodic phses [53]. This my be prevented by introducing strong disorder to the system, which leds to loclized phses tht void the problems of infinite heting [54 60]. Moreover, s demonstrted in Ref. [60] in model of one-dimensionl 1D) spin chins with Ising symmetry, MBL llows shrp phses to be defined by considering order in the Floquet eigensttes. Two of the phses identified in this work hve no sttic nlog, nd re in fct connected to Floquet topologicl insultors. In this pper, we ttempt to clssify intercting, symmetry-protected Floquet systems, focusing minly on Abelin symmetry groups nd systems of bosons in one dimension, nd building on intuition gined from studying the nonintercting problem [50] nd the disordered Mjorn Floquet chins in Ref. [44, 60]. Our pproch considers the decomposition of generl unitry evolution into constnt Hmiltonin component, which my generte topologicl order nlogous to sttic SPT phse, nd loop component, which cn introduce dynmicl SPT order tht is only possible in driven systems. We present expressions for topologicl invrints nd illustrte the bulk-edge correspondence for these phses, furnishing our results with exmples from clss D nd Z N Z N SPTs. Severl recent ppers consider the clssifiction of Floquet SPT phses using different methods [61 63], nd we find tht our results re consistent with these. The structure of this pper is s follows. We begin in Sec. II by discussing some generl fetures of timedependent systems nd introducing the concept of unitry loop evolution, which we will require in the min text. In Sec. III we consider fermionic systems in clss

2 D, which ct s bridge between free nd intercting Floquet systems. We then study bosonic systems with Z N Z N symmetry in Sec. IV nd more generl Abelin groups in Sec. V, where we lso discuss bulk invrints nd bulk-edge connections. We end with discussion in Sec. VI. II. PRELIMINARY DISCUSSION The evolution due to locl, time-dependent Hmiltonin Ht) cn be described by the unitry opertor t ) Ut) = T exp i Ht )dt, 1) where T is the time ordering opertor nd t usully runs from 0 to T. We ssume tht the Hmiltonin Ht) t every instnt is invrint under group of symmetry trnsformtions, G. We write the instntneous eigensttes of Ut) s ψt) nd note tht these correspond to instntneous qusienergies ɛt) through Ut) ψt) = e iɛt) ψt), where π < ɛt) π. It is conventionl to define the effective Floquet Hmiltonin t the end of the evolution s H F, obtined through UT ) = exp ih F T ), when this gives meningful i.e., locl nd symmetry-preserving) Hmiltonin. For the purposes of this pper, we restrict our discussion to unitry evolutions which, for closed system, cn be expressed s exp ih F T ) for some Floquet Hmiltonin H F. As we will see below, this implies tht the unitry for the open system cn be written s product of two fctors, where one fctor is of the form exp ihf o T ) nd the other is n effective edge unitry opertor. We wnt to understnd nd clssify the distinct types of robust effective edge opertor tht cn rise from such unitry evolutions. The requirement tht the effective edge unitry hve observble physicl consequences plces some constrints on the form of H F, which we do not study in detil here. One expects tht integrbility or MBL is necessry, lthough likely not sufficient. One wy to generte edge modes protected from bulk trnsitions by qusienergy gp or n effective mobility gp in n MBL phse, is by evolving constnt SPT Hmiltonin in time. Phses produced in this wy hve protected edge modes, resulting in qusienergy spectrum with degenerte multiplets t every level. The edge modes, which re connected to this vector spce, then correspond to projective representtion of the group, nd the different phses correspond to equivlence clsses of projective representtions [10]. However, free fermion studies point to set of more interesting possibilities, which we will find correspond to nontrivil unitry loops. We define loop to be unitry evolution of the form in Eq. 1) tht, for closed system, stisfies UT ) = e iφ I for some phse φ. In much of wht follows, we will ssume φ = 0 without loss of generlity. In this wy, the qusienergies t the end of unitry loop 0 coincide, nd the effective Floquet Hmiltonin H F = 0. Nevertheless, fter the loop evolution of n open system, there my be edge modes present t different qusienergy vlues, which leds to richer set of possible topologicl phses. The utility of studying loops becomes evident when we consider n end point unitry for closed system tht cn be expressed s UT ) = exp ih F T ), where H F is sttic nd locl Hmiltonin. If this is the cse, then we cn continuously deform the unitry evolution into composition of unitries of the form U = C L, where L is loop nd C is the evolution due to the sttic Hmiltonin H F. We define the composition of two unitries U 1 U in the usul wy, s the complete evolution due to U 1 followed by the complete evolution due to U, with both components ppropritely rescled so tht the evolution runs from t = 0 to t = T see Ref. 50 for n explicit construction). If the closed system unitry Ut) results from the evolution due to the time-dependent Hmiltonin Ht), nd hs the end-point form UT ) = exp ih F T ), the continuous deformtion into the decomposed form U = C L my be obtined through the homotopy Hmiltonin +H F s 0 t T s/ H t, s) = 4H F s H t 3sT/4) 1 3s/4 ) T s/ t 3sT/4 3sT/4 t T. When s = 0, this Hmiltonin genertes the originl Ut), but when s = 1, the first hlf of the evolution is constnt scled) evolution corresponding to H F, nd the second hlf of the evolution corresponds to loop Lt). 1 The end point of the evolution remins fixed throughout. In nonintercting system, this unitry decomposition leds to rigorous clssifiction of unitry evolutions [50]. In n open system, the unitry evolution will differ from tht of the closed system t the edge, giving rise to n effective edge unitry. Notbly, for nontrivil stte, it my not be possible to write the endpoint unitry for the open system s U o T ) = exp ihf o T ) for locl, symmetry-preserving Hmiltonin HF o : some terms could connect sites t the two edges of the chin. In this cse, following the decomposition of Eq. ), the open system unitry tkes the form U o T ) = U eff exp ihf o T ), where U eff is the effective edge unitry tht derives from the loop component of the evolution. Edge unitries of loops hve some specil properties, s we shll see below, nd these re likely to persist when the loops re followed by unitries with locl integrls of motion. When the constnt Hmiltonin is MBL, the SPT clssifiction of 1D ground sttes [1, 16 18] cn be ) 1 We ssume in this rticle tht the constnt evolution precedes the loop evolution in the decomposition of generl unitry. Similr rguments my be obtined using the opposite ordering.

3 3 extended to ll eigensttes [33, 34] with detectble consequences [34, 35]. The effective edge unitry of loop followed by constnt MBL Hmiltonin of some SPT phse will not in generl hve the direct product structure of bulk times edge. However, it is likely tht the edge modes of the loop will continue to hve some observble signtures, such s those considered in Ref. 35. The clssifiction of possible phses of such systems, s obtined from studying the distinct edge unitries, then hs nturl product structure, nl) H G, U1)) where the second fctor origintes from the usul sttic clssifiction of SPT phses [7, 10] nd the first fctor counts the number of distinct loops. Since the constnt evolution component is well understood, we will focus our study on loops. As noticed previously, generic Floquet systems driven with intercting Hmiltonins suffer from heting, nd this is certinly true for loops, which re fetureless nd hence, lredy t infinite temperture. However, one cn construct n MBL Floquet drive which voids this problem by piring loop with constnt MBL Hmiltonin. More generlly, loops re of interest nd my be defined even in systems without MBL. In this work, we restrict ourselves to one dimension where progress cn be mde most esily. We study nontrivil loops in vriety of contexts nd chrcterize them by bulk topologicl invrints. III. CLASS D SYSTEMS To motivte the ides described in this pper, we begin by studying the exmple of clss D free fermionic systems [43, 44]. These hve simple description in terms of nonintercting Mjorn fermions, nd llow us to develop the intuition required for the mny-body picture tht we will use in lter sections. A. Single-prticle Picture As motivted in Sec. II, we will focus on forming nontrivil unitry loops, which will in turn llow us to generte systems with dynmicl topologicl order tht re unique to time-dependent systems. In free-fermion system, loops of this kind cn lwys be generted by evolving with nontrivil, flt-bnd SPT Hmiltonin for the first hlf of cycle, nd trivil, flt-bnd Hmiltonin for the reminder of the cycle. We will illustrte this construction with specific exmple below. A generl clss D Hmiltonin my be written in terms of Mjorn fermions [65] s H = i i A ij j, 3) 4 ij We note, however, tht robust edges my lso exist for systems without MBL [64]. where i = i re Mjorn fermions nd A ij is n ntisymmetric mtrix. In the sttic system, this Hmiltonin is clssified by studying the symmetric momentum points, k = 0 nd k = π, bout which the Hmiltonin hs the full symmetry of clss D. By considering the pproprite homotopy groups, it cn be shown tht sttic systems in clss D hve twofold clssifiction described by n element of Z [65]. In the Floquet cse, we form nd define) the Mjorn unitry evolution opertor through Ot) = T exp t 0 A ijt )dt ). Since A ij is ntisymmetric, nd O0) is the identity, we see tht Ot) belongs to the specil orthogonl group SOn), if n is the number of bnds in the BdG Hmiltonin. We cn clssify the Floquet loop by gin studying the behvior t the symmetric momentum points, k = 0 nd k = π. Since we re considering loops, Ot) must end t the identity mtrix, nd so the clssifiction is given by the fundmentl group π 1 SOn)), which for n > 1 is Z. We see tht there is Z invrint ssocited with both k = 0 nd k = π, yielding four possible phses in totl. The set of loops is therefore clssified by e 1, e ) Z Z for n > 1, with ech e i tking vlues from {0, 1}. The topologicl invrint tht determines the existence or lck) of nontrivil boundry modes for n open system is given by the sum modulo ) e 1 + e Z. The extr Z element, which does not correspond to the edge modes, my be regrded s wek topologicl invrint of the drive. 3 The nontrivil Floquet loop is chrcterized by single Mjorn π mode in the open system, which, in the bsence of dditionl sttic SPT order, is lso ccompnied by Mjorn zero mode. For unitry evolution tht is not loop, there my lso be edge mode contributions from subsequent constnt Hmiltonin evolution, s described in Sec. II. We now illustrte these ides by describing specific model in clss D, with further detils given in Appendix A. We consider one-dimensionl fermionic chin of length K tht hs two-stte Hilbert spce t ech site, with the corresponding nnihiltion opertors on site j being j nd b j. 4 Following Kitev [65], we cn define two sets of Mjorn opertors through j 1 = j + j, b j 1 = b j + b j, j = j j i b j = b j b j, 4) i 3 Note tht for n = 1, the clssifiction of loops is Z Z, but the Z invrint tht describes the edge is still given by e 1 + e ) mod. 4 We choose two-stte Hilbert spce so tht the resulting phse hs nontrivil dynmicl SPT order but trivil sttic SPT order. This requires Mjorn mode t ɛ = 0 nd Mjorn mode t ɛ = π.

4 4 which stisfy =. Let H 1 = j j 1 ) = i j 1 j) j j H = 1 ) j j+1 j+1 j + j j+1 + j+1 j = i j j j j+1 ), 5) with H 1b, H b defined identiclly in terms of the b opertors nd b Mjorns. H 1 nd H respectively correspond to the trivil nd nontrivil phses of the 1D clss D superconductor. We now evolve the system with H 1 = H 1 + H 1b for 0 t π nd with H = H + H b for π t π. We notice tht for n open system, H nd H b ech hve Mjorn edge modes. The evolution by the trivil Hmiltonin H 1 pushes the Mjorn mode of subsystem to qusienergy ɛ = π, while the Mjorn mode of subsystem b hs been pushed to qusienergy ɛ = π 0. Thus, we see tht t t = π, Mjorn mode t both zero nd π emerge, which persist until the end of the evolution. On the other hnd, evolving the closed system until t = π leds to unitry tht is the identity up to n overll phse fctor). In this wy, the two-prt drive described by Hmiltonins H 1 nd H produces nontrivil unitry loop with dynmicl SPT order. This order is mnifested in the properties of the unitry evolution opertor. At the end of the loop, the unitry describing the closed system differs from the unitry describing the open system through terms tht ct t the ends of the chin. In the Mjorn lnguge, it is cler tht the term i K 1 from H nd the term i b K b 1 from H b cnnot ct in the open system, since they would connect the two edges. If we define new fermion opertors through d = 1 K + i 1 ), d = 1 K i 1 ) e = 1 b K + i1 b ), e = 1 b K i1 b ), 6) we see tht the open nd closed unitry opertors U op nd U cl, respectively) re relted through U op π) = e [ π K 1 +b K b 1)] Ucl π) 7) = e [ iπ d d dd )+iπe e ee )] Ucl π), where U cl π) is the identity up to some overll phse fctor. The detils of this clcultion re presented in Appendix A. Thus, we see tht the effective unitry t the edge is U eff π) = e [ iπ d d dd )+iπe e ee )], 8) where the second term in the exponent merely provides n overll fctor of 1. We will use the concept of n effective edge unitry throughout this rticle to discuss topologicl order nd define invrints. We will be prticulrly interested in its behvior under the symmetry trnsformtions of the relevnt group G, which for clss D is fermion prity. We denote the prity opertor by ˆP, which cn be expressed s product of the onsite prity opertors of ll sites, ˆP = ˆP j j, with ˆP j = ) ) ij 1 j i b j 1 j b. Now, if we define n effective prity opertor ˆP L t the left edge of the open chin, which is product over finite set of onsite prities in the vicinity of this edge, then for the effective edge unitry in Eq. 8), we find { ˆP L, U eff } = 0. 9) For trivil effective edge unitry, we would insted find [ ˆP L, U eff ] = 0. In the discussion bove, we hve considered specific model Hmiltonin, whose unitry evolution exhibits cler signtures of the underlying dynmicl SPT order. While this model ws chosen for its simplicity, the resulting signtures will pply more generlly. First, Mjorn π mode will exist t ny point fter dynmicl phse trnsition hs occurred until nother equivlent phse trnsition occurs). In the model considered here, the existence of the π mode could be shown nlyticlly t the specil point t = π. At other times, nd for more generl models, the presence of Mjorn π mode my be inferred from the properties of the mnybody spectrum, which we discuss below in Sec. III B 3. We my lso consider perturbing the system with locl short-rnge) Hmiltonin, H L, which cts ner the left edge t the end of the loop nd which hence commutes with suitbly defined ˆP L ). After time t, the effective edge Hmiltonin will hve evolved to U eff = e ihlt U eff. However, it is cler tht the new effective unitry still stisfies { ˆP L, U eff } = 0. In this wy, the effective edge unitry is chrcterized by its commuttion or nticommuttion properties with the prity opertor, ˆP L, nd is robust under locl perturbtions tht preserve symmetry. We cn lso consider threding π-flux through the closed system, nd the effects of this will provide nother signture of the SPT order which is robust to perturbtions. By explicit clcultion, detiled in Appendix A, it cn be shown tht the unitry of the closed system t the end of the cycle with flux differs from the unitry without flux by n overll fctor of e iπ. We write this reltionship s U 0 U π = e iπ I. B. Mny-body Picture We re now in position to formulte the problem in mny-body setting, returning to the generl set of systems in clss D. Although we will t first not introduce ny dditionl interctions, we note tht much of the discussion tht follows is pplicble to both free nd

5 5 intercting systems. This is in much the sme wy tht the clssifiction of sttic clss D free-fermion systems crries over to the intercting cse [65]. 1. Bulk Invrint To find bulk invrint in the mny-body picture, we consider the unitry evolution corresponding to the generl Hmiltonin given in Eq. 3). A nontrivil Mjorn Floquet loop stisfying OT ) = I corresponds to fermionic unitry loop stisfying UT ) = I. For nontrivil evolution, we therefore see tht U 0 T )U πt ) = e iπ I, which grees with our findings for the specific model considered in the previous section. For trivil loop, the equivlent product yields U 0 T )U πt ) = I, nd so we cn define the bulk invrint, k {0, 1} Z of n intercting clss D system through U 0 T )U πt ) = e ikπ I. We expect this invrint to chrcterize loop even in the presence of prity-preserving interctions, nd we will rgue tht this is indeed the cse below.. Edge Picture We found previously tht nontrivil unitry loop corresponds to the ppernce of zero nd π Mjorn mode in the single-prticle spectrum of n open system. In the mny-body spectrum, this is mnifested s twofold degenercy t every qusienergy, ccompnied by π-degenercy. By π-degenercy, we men tht every stte in the spectrum hs counterprt seprted by qusienergy π, which differs only by boundry Mjorn mode. This implies tht there re two sttes of differing prity t ech energy, which we expect to persist in the presence of prity-preserving interctions. Without loss of generlity, we will restrict our discussion to nontrivil loops tht cn be written s two-prt drive, comprising evolution by trivil Hmiltonin followed by evolution by n SPT Hmiltonin more generl unitry loops cn be considered s compositions of such two-prt drives). In these cses, the unitry for the closed system is the identity up to phse), while the nontrivil behvior t the edge of the open system cn be cptured by n effective unitry of the form given in Eq. 8), U eff = e [ iπ d d dd )+iπe e ee )]. 10) In the generl intercting cse, the opertors d nd e cn be mny-body fermion opertors tht depend nontrivilly on the opertors of the originl Hmiltonin, but which trnsform under prity in the sme wy s in the single-prticle cse. We regrd the π- nd 0-Mjorn modes s rising out of dynmicl phse trnsition tht does not chnge the sttic SPT order of the system. 3. Phse Trnsitions We now consider the phse trnsitions tht must ccompny chnges in the dynmicl topologicl order of the system. Formlly, we note tht single unitry evolution Ut) for t T cn lso be regrded s fmily of unitries U t, s) with U t, s) = Ust) for st T nd where s is continuously vried. If, s in the model system bove, the closed system unitry U T, s) cn be expressed s the exponent of nturl Floquet Hmiltonin for vlues of s excluding some point s c, but the open system unitry does not hve such nturl Floquet Hmiltonin expression for s > s c, then we my regrd this s phse trnsition of the unitry tht occurs t s c. The trnsition point U T, s c ) is then equivlent to UT/s c ). We will loosely refer to this s the existence of criticl point in the originl Floquet evolution Ut) t the point t = T/s c t c. In other words, we ssume tht for t > t c the open system unitry tkes the form U o t > t c ) = U eff exp ihf o t), with the dynmicl order chrcterized by the effective edge unitry U eff. We ssume tht the system undergoes such trnsition t t = t c, nd consider the instntneous mny-body qusienergy spectrum of Ut) cross t c. A similr pproch ws used to study Floquet trnsitions in the nonintercting cse in Ref. [48]. We ssume tht the sttic SPT order of the eigensttes of the unitry evolution opertors before nd fter the trnsition is trivil, so tht there re no Mjorn modes before t c nd one Mjorn mode ech t zero nd π fter t c. For simplicity, we will lso ssume tht there re no other phse trnsitions during the evolution. These ssumptions suggest tht there is no difference between the closed nd open systems for t < t c, since in this region the evolution is topologiclly trivil. At the trnsition point t c, the emergence of Mjorns is ssocited with bulk 0, π)-trnsition, where ech bulk eigenstte ψ becomes prt of multiplet of four sttes tht my be written { ψ, d ψ, e ψ, d e ψ } fter the trnsition. If two of these sttes hve mnybody qusienergy ɛt c ), the other two hve qusienergy ɛt c ) + π modulo π). The two sttes in ech pir hve opposite prity. This degenercy pttern in the mnybody spectrum will persist until nother phse trnsition occurs. We now consider wht hppens in the open system if we reconnect the edges t some point just beyond t c with Hmiltonin of the form H = h d d d dd ) + h e e e ee ). 11) The effective edge unitry t time t > t c is then given by U efft) = e [ iπ d d dd ) it t c)h dd d dd )] e [iπe e ee ) it t c)h ee e ee )], 1)

6 6 which is the effective unitry from Eq. 10) multiplied by the evolution due to the edge-closing Hmiltonin, Eq. 11). We note tht if we close the system with ntiperiodic boundry conditions, the sign of the coefficients h d nd h e chnges reltive to closing the system with periodic boundry conditions. This llows us to mke the following sttements bout the spectrum immeditely fter t c : 1. In the system with π flux which is equivlent to imposing ntiperiodic boundry conditions), the spectrum fter t c my be obtined from the system with zero flux by shifting the entire spectrum by π.. If the zero flux spectrum is shifted in this wy, then sttes re mpped onto sttes with the sme prity. 3. If h e 0, there is prity shift in the ground stte or ny other prticulr stte) of the multiplet in the system with flux. These spectrl fetures will rise in generic clss D system with dynmicl SPT order. In Fig. 1, we show the splitting of the multiplet schemticlly for both periodic nd nti-periodic boundry conditions, ssuming the complete evolution forms unitry loop. We note tht if there re no other phse trnsitions, then U π T )U 0 T ) = I, demonstrting the correspondence between the bulk invrint nd the behvior of the edge modes. Finlly, we note tht 0, π)-trnsition of this form cn be viewed s fourfold energy level crossing in the folded qusienergy spectrum, which plots mny-body qusienergy vlues modulo π in the rnge 0 ɛt) π. Notbly, ordinry phse trnsitions would pper s pirwise crossings in both the folded nd the unfolded spectr. IV. SYSTEMS WITH Z N Z N SYMMETRY Let us now consider Floquet phses of bosonic systems tht re symmetric under the group Z N Z N, following similr strtegy to the discussion of clss D. We will first construct pproprite flttened in this context, integervlued) SPT nd trivil Hmiltonins, which hve useful properties under time evolution. Constructing loops out of these, we will find nontrivil topologicl phses where there re qusienergy degenercies t multiples of π/n, relted to edge modes in the open system. In the sttic cse, the group cohomology clssifiction predicts N different SPT phses, which we lbel by some integer c {0, 1,,..., N 1}. As shown in Appendix B, we cn construct integer-vlued SPT Hmiltonins H,b c, defined on spin chin with K sites, which generte these sttic SPT phses nd which hve the following properties under time evolution: For closed system, U π N ) πi = e N Hc,b) = [V To )] [V T e )] b, 13) ϵt) π π 0 - π -π ϵt) π π 0 - π -π e ψ ψ e d ψ d ψ t 0 e ψ ψ e d ψ d ψ t 0 ) b) FIG. 1. Schemtic splitting of Mjorn multiplet t 0, π)-phse trnsition s the system is closed with ) periodic nd b) ntiperiodic boundry conditions. Blue lines correspond to sttes with even prity, while red lines correspond to sttes with odd prity. Spectrum b) cn be obtined from spectrum ) through shift by π. See min text for detils. where T o, T e re genertors of the two Z N groups corresponding to odd or even sites), nd V T o/e ) re the globl unitry opertors corresponding to T o/e defined explicitly in Appendix B). We cn lso form trivil Hmiltonins, H,b triv which hve the property for both open nd closed systems) tht U π N ) = e πi N Htriv,b ) = [V To )] [V T e )] b. 14) By construction, we cn form loop by evolving the system in time with the Hmiltonin { H c Ht) =,b for 0 t π/n, 15) for π/n t 4π/N, H triv ā, b where ā = N is the inverse element of. Choosing c = 1 we drop the superscript c for the rest of this section), we find tht the result of this evolution for n open system is described by the effective edge unitry [ U eff = exp i π ] N H 1 + bh K ), 16) where H 1 nd H K re opertors cting on the two ends of the open system, whose explicit form is given in Ap- T T

7 7 pendix B. Then, every stte in the open system with qusi-energy ɛ is ccompnied by sttes tht differ by edge modes nd which hve qusienergies ɛ+πk 1 /N + k b/n), where k 1/ {1,,... N 1}. This is signture of the topologicl phse tht will pper in the mny-body spectrum. Anlogous to the clss D discussion, we my lso consider the effect of twisting the boundry conditions of the closed system. In our model, the effect of T o -twisted boundry conditions cn be relized by trnsforming the opertors H 1, H K by H 1 X 1 H 1 X 1, H K X 1 H K X 1, where X 1 is the unitry opertor corresponding to the symmetry T o t site 1. We now consider evolving the open system until t = 4π/N, before reconnecting the ends of the chin with n opertor H. The effective unitry t the edge fter this point is then given by U efft) = e [ it 4π/N)H ] e [ iπ N H1+bH K)]. 17) On the other hnd, if we reconnect the system with T o -twisted boundry conditions, the Hmiltonin H chnges to H X 1 H X 1 nd the resulting effective unitry chnges to U eff t) U T o t) where U To t) = X 1 e ] [ it 4π/N)H X iπ 1e[ N H1+bH K)]. 18) Since detx 1 ) is non-vnishing, nd X 1 H KX 1 = H K I + NP 0, 19) where P 0 is the zero-energy subspce of H K, we see tht U To t) hs the sme spectrum s Ũ To t) = e [ it 4π/N)H ] e [ iπ N H1+bH K I))], 0) which is shifted by phse πb/n compred to the untwisted unitry opertor. Similrly, under T e -twisted boundry conditions, the spectrum of the effective unitry is shifted by π/n, nd generl effective unitry my be chrcterized by these phses. A relted clcultion tells us tht evolving the closed system through the entire loop with T o -twisted boundry conditions gives unitry e iπb/n I, while evolving with T e -twisted boundry conditions gives unitry e iπ/n I. As we will rgue in the next section, these phses define bulk topologicl invrint of the evolution. Since there re N possible vlues for ech phse, there re totl of N distinct unitry loops, corresponding to N distinct topologicl phses. V. GENERAL ABELIAN GROUPS For generl finite Abelin group with p fctors, G = Z n1 Z n... Z np, we follow the sme strtegy s before to construct models for loops with dynmicl topologicl order. We first form n expnded group, G, which contins G s subgroup, for which nontrivil SPT cn be constructed whose ground stte multiplet contins sttes tht cn be lbeled by ll the irreducible representtions of the group G. SPTs in one dimension re clssified by the second cohomology group H G, U1)). Using the result [66] H G, U1)) = 1 j<k p Z n j,n k ), 1) where n j, n k ) indictes the highest common fctor of n j nd n k, we see tht the group G = Z n1 Z n... Z np Zñ, where ñ is the lowest common multiple LCM) of the set {n 1,... n p }, hs SPT sttes with ground stte multiplets tht contin ll the irreducible representtions of G. When coupled with pproprite trivil Hmiltonins to form loops, in much the sme wy s in the Z N Z N cse, these led to systems whose effective edge unitries for evolution on n open chin cn be chrcterized by the phses {α j }, where ech α j lbels n irreducible representtion of T j, the genertor corresponding to the group fctor Z nj. These models re chrcterized by set of robust edge modes in the open system, which cn be described by n effective edge unitry U eff. For loop, since H F = 0, this implies tht U = U eff. Following our previous discussion, we introduce twisted boundry conditions s follows: first imgine cutting closed chin t some site, nd lbeling the two edges L nd R. There exists set of terms in the Hmiltonin for the closed chin tht connect the two edges L nd R, whose sum we lbel by H edge. For ech genertor T j, there is then n opertor X L,j = i L V it j ), which is product over the onsite unitry opertors corresponding to T j over finite set of sites on the left edge, such tht the effect of the twisted boundry conditions is equivlent to trnsforming H edge s H edge X L,j H edge X L,j ). A. Effective Edge Unitries nd Spectrl Shift For the models we hve just introduced, one cn show, generlizing the rguments of the Z N Z N cse, tht the effective edge unitries stisfy the property X L,j ) U eff X L,j = e iαj U eff. ) We will rgue below tht this property is true more generlly: Any effective edge unitry of loop must stisfy n eqution of the form of Eq. ) for some X L,j which hs support only t the left edge. Eq. ) implies tht if ɛ is qusienergy of U eff then so too is ɛ + nα j, for ny integer n. Furthermore, we see tht α j must in generl be multiple of π/n j, where n j is the size of the Abelin group ssocited with T j, for ny U eff tht stisfies Eq. ). We cn lbel n effective unitry by the set of phses {α j } it is multiplied by under trnsformtion with the {X L,j }, writing U eff {α j }). We cn thus obtin distinct U eff for every irreducible representtion of the group G. The opertor U eff must lso commute with V T j ) = i V it j ) for ny j, where the product is now over ll sites.

8 8 Let us now study the implictions of the effective unitry for the system. Suppose the effective edge unitry t some time is U eff. We then consider connecting the system t t this point nd evolving it for n infinitesiml time t. The resulting evolution is cptured by the trnsformtion U eff e ih t U eff, where H is Hmiltonin tht hs support on both L nd R. The effect of crrying out the sme process with T j -periodic boundry conditions is U eff X L,j e ih t X L,j U eff. From Eq. ), this hs the sme spectrum s e ih t U eff but is shifted by α j. Since the unitry of the full open system cn be written s U o T ) = U eff exp ihf o T ), the spectrum of the full system under T j -periodic boundry conditions is shifted reltive to the untwisted cse by α j. Now consider the effect of perturbing the system with some locl short-rnged) term with support only on L), described by Hmiltonin H L. At time t fter the perturbtion is pplied, the effective unitry U eff hs evolved to U eff = exp ih Lt)U eff. Since H L is symmetric under V T j ) = i V it j ), then V T j )e ihlt [V T j )] = e ihlt. However, since H L only hs support on L, we see tht It follows tht V T j )e ih Lt [V T j )] = X L,j e ih Lt X L,j = e ih Lt. 3) X L,j U effx L,j = e iαj U eff. 4) Thus we see tht the defining property of the effective unitry given in Eq. )) is stble ginst locl perturbtions. B. Phse Trnsitions Suppose we hve unitry evolution such tht for some rnge of times t, the system cn be described by effective edge unitries U eff {α j }, t) nd instntneous bulk Floquet Hmiltonins H F t)) such tht the effective edge unitry is nontrivil t some time t 1 nd trivil t lter time t. As we now rgue, bulk phse trnsition must then hve occurred t some point t c between t 1 nd t. 5 We first consider evolving the open system until time t 1, before connecting the edges of the system nd evolving it further for n infinitesiml time δ, with both T j - twisted nd regulr boundry conditions. It follows tht the spectrum of the trnsformed edge unitry t t 1 + δ in the twisted cse is shifted by α j reltive to the untwisted cse. On the other hnd, if we evolve the open 5 We re considering the idelized sitution of phse trnsition occurring t single point. More generlly, extended sttes my persist for some intervl of time, but our discussion cn be esily extended to this more generl cse. system until t, nd then gin evolve it for n infinitesiml time δ, closing with twisted nd regulr boundry conditions, there is no reltive spectrl shift between the two boundry conditions. This chnge in the boundry behvior implies bulk phse trnsition t some point t c between these two times, where by continuity, the bulk spectrum t t c is symmetric under shifts by α j. At the phse trnsition point, the edge modes t the two ends of the open system re connected through bulk modes. This trnsition is of the type nontrivil to trivil, but it is cler tht similr process occurs for trnsitions of the type trivil to nontrivil. At phse trnsition where nontrivil U eff chrcterized by {α j } ppers or disppers, there must be bulk modes tht re in one-to-one correspondence with, nd which hve the sme spectrl gps s, the edge modes represented by U{α j }). From our previous discussion following Eq. ), eigensttes of the bulk spectrum t the phse trnsition point must be ccompnied by other modes t qusienergies which re nα j prt. In generl we my consider phse trnsitions of the type {α j } {β j } by decomposing them into sequences of phse trnsitions of the type considered bove. C. Bulk Topologicl Invrint We now consider evolving the closed system through loop so tht the unitry evolution opertor t the end of the evolution is I up to phse fctor e iφ, which we tke to be 1 without loss of generlity below. If we now evolve the sme system with T j twisted boundry conditions, then we clim tht the unitry opertor t the end of the evolution hs to be of the form e iθj I. To see this, consider the following rgument: the T j -twisted periodic boundry conditions cn be pplied by trnsforming the Hmiltonin term H edge,i cross prticulr link i of the closed chin s ), X Li),j H edge,i XLi),j where Li) indictes the left side of n edge cut t site i. The Hmiltonins obtined by performing these trnsformtions t different sites i of the chin re relted by unitry trnsformtions. In the cses we consider, the unitry loop cn be regrded s finite depth quntum circuit or n LU trnsformtion) [10]. The effect of this trnsformtion in sufficiently long chin t prticulr site i must led to n effective unitry which is of the form U i, where U i cts nontrivilly only in the vicinity of i [67]. The unitry obtined by crrying out this trnsformtion t distnt site j, U j, must be relted to U i through unitry trnsformtion. Thus, U i must be fetureless nd cn only contribute n overll phse i.e. it must be of the form e iθj I. Further, introducing n j such twists t distnt points long the chin must hve the sme effect s introducing T j ) nj twist, which brings us bck to the untwisted cse. Thus, e inθj = 1, nd the discrete llowed vlues

9 9 of θ j imply tht it is topologicl invrint tht cnnot chnge under continuous deformtions of the loop. Thus, the set {θ j } my be used s bulk topologicl invrint for the dynmicl phse of the system. The number of distinct invrints is clerly equl to the number of distinct irreducible representtions of the group G. It lso follows tht if we consider two loops with invrints {θ j } nd {β j }, the loop obtined by running these unitry evolutions in succession hs n invrint {θ j + β j }. A similr rgument cn be used to confirm tht for generl unitry loop of fermions with clss D symmetry, U 0 U π = e ikπ I for integer k, s we found in Sec. III. For those systems, π flux insertion cn be effected by the pproprite Peierls substitution for set of edge terms in the vicinity of some site i. Agin, we cn rgue s we did bove tht this leds to unitry tht differs from tht of the closed loop only in the vicinity of i nd tht this difference must be n overll phse. To rgue tht the phse must be n integer multiple of π, we note tht introducing π flux brings us bck to the cse with no flux insertion nd tht the phse ccompnying the effective unitry must hve doubled compred to the π flux cse. D. Bulk-edge Connection Consider unitry loop with bulk topologicl invrint {θ j }, which we evolve s follows. We first evolve the system with open boundry conditions long the entire loop. This leds to some effective unitry U eff t time T. We then reconnect the edges with some Hmiltonin H edge nd evolve it for some dditionl time T, such tht the net unitry t T +T is the identity. This implies tht e ih edget U eff = I. 5) We my imgine sequence of moves through which we continuously deform the originl loop of the closed system to the bove evolution, thus preserving the bulk invrint. If we now consider the relted evolution where we evolve it with open boundry conditions until time T nd with T j -twisted boundry conditions from time T to T + T, the unitry t the end of the evolution must be X L,j e ih edget X L,j U eff = e iθj I. 6) This cn be rewritten using Eq. 5 s X L,j U effx L,j = e iθj U eff, 7) which hs the sme form s Eq. ), with {θ j } identified with {α j }. This rgument lso shows us tht ny effective edge unitry hs to stisfy n eqution of the form in Eq. ), s climed before. We cn lso use this thought experiment to confirm the bulk-edge connection in the clss D cse. We imgine evolving the system with open boundry conditions long the entire loop, resulting in n effective unitry, U eff t time T. We then reconnect it with some edge Hmiltonin H edge, s in our previous discussion, nd evolve it until some further time T + T such tht the finl unitry of the entire evolution is the identity, leding to the condition of Eq. 5). As before, we rgue tht we cn continuously deform the originl evolution of the closed system through the loop to the one just described. The effect of redoing the evolution until T + T, but now with π flux, must then led to the unitry : P L e ih edget P L U eff = e ikπ e ih edget U eff. 8) where P L is product of on-site prity opertors on the left edge of the system. This cn be rerrnged to give P L U eff P L = e ikπ U eff, 9) which shows tht for non-trivil loop, U eff must stisfy Eq. 9) nd must commute with P L for trivil loop. VI. DISCUSSION In this rticle, we hve considered the clssifiction of Abelin Floquet SPT phses in one dimension. Using systems with clss D symmetry s n intuitive exmple, we went on to discuss bosonic SPT phses with Z N Z N symmetry nd then, more generlly, symmetry under ny finite Abelin group. Our pproch ws to consider unitry evolutions s composition of constnt Hmiltonin evolution with loop. The first component, being in one-to-one correspondence with sttic Hmiltonin, reproduced the clssifiction scheme of sttic SPT phses. More interesting dynmicl SPT phses could be obtined by constructing nontrivil loop evolutions. In the cses discussed in this pper, we constructed explicit models tht could be used to generte such nontrivil loops, nd outlined how these models my be extended to more generl symmetry groups. The loop construction llowed us to identify dynmicl topologicl order by studying the behvior t the ends of n open system, which we encoded in n effective unitry tht cted only t the edge. Effective edge unitries of this form re distinguished by discrete set of chrcteristic trnsformtion properties tht re protected from locl perturbtions. This leds to spectrl piring for the mny-body eigensttes provided the effective edge unitry commutes with the unitry of the bulk, s one might expect for MBL systems. For closed systems, we rgued tht unitry loops could only chnge by certin quntized phses under grouptwisted boundry conditions. We interpreted these phses s bulk invrints: robust quntities, stble under locl perturbtions, tht could be used to count nd lbel the dynmicl Floquet SPT phses. The distinct sets of U1) phses obtined under such twists correspond to distinct one-dimensionl irreducible representtions of the group. Moreover, these bulk invrints were shown to be relted to the effective edge unitry, providing bulk-edge correspondence in these systems. We rgued tht invrints of this form cn only chnge through bulk phse trnsition, nd we enumerted the

10 10 distinctive properties nd gp structures tht rise in the qusienergy spectrum t such points. Overll, we expect tht unitry evolutions protected by symmetry G re clssified ccording to the product nl) H G, U1)), where the second fctor clssifies the sttic component, nd the first fctor clssifies the loop component. Distinct loop components re given by the different irreducible representtions of the group G. It is likely tht the loop construction would lso permit the clssifiction of dynmicl SPT phses under more generl conditions thn those discussed here. We leve discussion of more exotic symmetries, including higher dimensions, time-reversl symmetry nd non- Abelin groups, to future work. Finlly, we note tht our clssifiction method nd discussion differs from relted work [61 63] in number of significnt wys the formultion in terms of loops, the bulk invrints, the bulk-edge correspondence nd dynmicl phse trnsitions), lthough the counting of the phses is consistent with these other results. Notbly, the product structure of our clssifiction is reminiscent of the product structure of the clssifiction given in Refs. 6 nd 63 for G with unitry symmetries only. We expect these clsses of SPT phses to be in one-to-one correspondence. Although we frme the discussion in the context of Floquet systems, our pproch does not require the system to be time periodic. Rther, we consider the unitry evolution due to generl time-dependent Hmiltonin, nd study the instntneous topologicl order tht my exist t given point in time. In order to obtin the loop component from generl unitry evolution, we ssumed tht in the closed system, the end-point unitry could be written s UT ) = exp ih F T ), with H F sttic nd locl Hmiltonin. While this is likely true for mny systems exhibiting MBL, it my lso hold more generlly. In this wy, we did not need to explicitly invoke ides from MBL, lthough such considertions would undoubtbly be importnt in periodiclly driven systems. Notbly, experimentl probes of SPT order, such s the observtion of persistent edge modes [34, 35], would seem to require MBL. We leve discussion of these issues to future work. ACKNOWLEDGMENTS We thnk C. W. von Keyserlingk nd S. L. Sondhi for shring copy of Ref. [61] in dvnce of publiction nd for mny fruitful discussions. We lso thnk P. Pdmnbhn nd A. Vishwnth for shring useful insights nd the Princeton Centre for Theoreticl Science for hospitlity. R. R. nd F. H. cknowledge support from the NSF under CAREER DMR nd the Alfred P. Slon foundtion. Note During the course of this work we becme wre of relted works in Refs. [61 63], whose results re consistent with those presented here. Appendix A: Further Detils on the Clss D Model In this ppendix, we give dditionl detils on the clcultions outlined in Sec. III. We begin with the Hmiltonin of Eq. 5), nd consider the evolution described in the min next. In the first prt of the evolution, for 0 t π, the Hmiltonins for open nd closed systems re the sme, nd the unitry evolution opertor is given by U o/c 1 π) = exp π K K j 1 j) + b j 1 j) b = K j 1 j ), where we hve mde use of the identity exp [t j k ] = cos t) + sin t) j k. The second prt of the evolution differs depending on whether we re considering the closed or open system. For the closed system, we find Uπ) c = exp π K j j+1) K + b j j+1) b = K j j+1 ). while for the open system we find U o π) = K 1 K 1 ) j j+1) K j j+1 ). Then, for the complete closed loop, we obtin U c π)u c 1π) = K j j+1) K j 1 j) = 1, mking use of the identity K 1 K 1 j j = 1) K+1. For the complete open loop, on the other hnd, we obtin U o π)u o 1 π) = K 1 U eff π)u c π)u c 1π). It my be verified tht U eff π) = K 1 is equivlent to the expressions given for the effective edge unitry in Eq. 8). We now consider threding π-flux through the closed system s suggested t the end of Sec. III A. In the first

11 11 prt of the drive, this does not ffect the unitry evolution, since the Hmiltonin only consists of on-site terms. In the second prt of the drive, we implement the new boundry conditions by chnging the sign of the term tht connects site K to site 1. Specificlly, the Hmiltonin for the second prt of the drive becomes H = i K 1 j j+1) i K 1, nd similrly for H b. This now corresponds to the unitry evolution Ũπ) c = exp π K 1 j ) j+1 K 1 K 1 + b j j+1 b ) K 1 = K j j+1 ). Finlly, for the complete closed loop, we obtin Ũ c π)u c 1π) = = +1, K j j+1) K j 1 j) which hs n overll phse difference of e iπ compred to the flux-free system. Appendix B: Integer-vlued SPT Hmiltonins for Systems with Z N Z N Symmetry In this ppendix, we outline the construction of integer-vlued Z N Z N SPT Hmiltonins. We strt from the model Z N Z N Hmiltonins given in Ref. [68]. These consist of chin of Z N vribles, divided into odd nd even sites, which hve the on-site terms Ho i,c = 1 { c } Z i 1) Xi Z i+1 ) c + H.c. B1) He i,c = 1 ) c } {Z i 1 ) c X i Z i+1 + H.c., where o e) stnds for odd even) sites. In these expressions, Z nd X re Z N generlistions of the Puli sigm mtrices, given explicitly in Ref. [69]. Different vlues for c modulo N) llow one to obtin the N different sttic SPT phses, nd the genertors of the Z N Z N symmetry re given by the opertors V T o ) = i o X i V T e ) = i e X i, B) where T o nd T e re the genertors of the globl Z N symmetries on the odd nd even sites, respectively. We emphsize tht ech term in Eq. B1) commutes with ny other term with the sme vlue of c, nd tht tht these lso commute with the symmetry genertors V T o ), V T e ). Ech term H i,c o/e hs eigenvlues ) ) ω + ω ω + ω ) E = 1,,,... B3) where ω = e πi/n is root of unity. We cn therefore rewrite Ho i,c s Ho i,c = [P 0,i ωn + ω ) n ] P n,i +..., B4) nd similrly for the even sites, where P n,i is the projector to the eigenstte subspce of Ho i,c with eigenvlue ωn +ω ) n. It is now possible to write down modified on-site terms through N 1 H i,c int,o = np n,i n=1 N 1 H i,c int,e = n P n,i, n =1 B5) which ech hve integer eigenvlues E = 0, 1,..., N 1. The complete modified integer-vlued Hmiltonins re written Hcl c = i o Hop c = i o H i,c int,o + H i,c int,e i e H i,c int,o + i e H i,c int,e, B6) where the primes on the summtions in the open system indicte tht the boundry terms corresponding to i = 1 nd i = K) should be excluded. Below, we will leve the open/closed lbel implicit if the mening is cler. Since ech term in the Hmiltonin commutes, we cn lbel ech eigenstte of the system by its eigenvlues under ech H i,c int,o/e. For the closed system, we write prticulr stte s λ 1, λ,..., λ K, where ech λ i tkes vlues from {0, 1,..., N 1}. Since there re N K sttes nd N K eigenvlue combintions, every stte is fully determined by its lbels. It is simple to verify tht the stte {λ i } hs energy E {λi} = i λ i nd is trnsformed under the symmetry genertors ccording to V T o ) {λ i } = e πi/n) i o λi {λ i } V T e ) {λ i } = e πi/n) i e λi {λ i }. B7) In the open system, there re still N K sttes but only N K eigenvlue combintions, since λ 1 nd λ K re not specified. This leves n N -fold degenercy in every eigenstte of the system. However, the degenerte