## Topological Phases of Many-Body Localized Systems: Beyond Eigenstate Order

David M. Long 1, ∗ and Dominic V. Else 2

1

Condensed Matter Theory Center and Joint Quantum Institute, Department of Physics, The University of Maryland, College Park, Maryland 20742, USA 2 Perimeter Institute for Theoretical Physics, Waterloo, Ontario N2L 2Y5, Canada (Dated: August 23, 2024)

Many-body localization (MBL) lends remarkable robustness to nonequilibrium phases of matter. Such phases can show topological and symmetry breaking order in their ground and excited states, but they may also belong to an anomalous localized topological phase (ALT phase). All eigenstates in an ALT phase are trivial, in that they can be deformed to product states, but the entire Hamiltonian cannot be deformed to a trivial localized model without going through a delocalization transition. Using a correspondence between MBL phases with short-ranged entanglement and locality preserving unitaries-called quantum cellular automata (QCA)-we reduce the classification of ALT phases to that of QCA. This method extends to periodically (Floquet) and quasiperiodically driven ALT phases, and captures anomalous Floquet phases within the same framework as static phases. The QCA framework further generalizes to include symmetry-enriched ALT phases (SALT phases), and provides a large class of soluble models suitable for realization in quantum simulators. In systematizing the study of ALT phases, we both greatly extend the classification of interacting nonequilibrium systems and clarify a confusion in the literature which implicitly equates nontrivial Hamiltonians with nontrivial ground states.

Many-body localization (MBL) [1-11] can protect robust, topological phases even in systems which are far from equilibrium [12, 13]. For instance, localized systems can support spontaneous symmetry breaking [1215], topological order [16-21], and nontrivial phases with external driving [22-32], including time crystals [3340]. In fact, any Hamiltonian built out of commuting terms, of the kind frequently employed to study topological order [41-48], becomes robustly localized by choosing strongly disordered coefficients for the commuting terms [20].

tems, ALT phases have only been discussed in the singleparticle context [50]. Nonetheless, we propose that the static interacting models of Refs. [47, 48, 55-57] in fact represent ALT phases upon disordering the Hamiltonian.

Two MBL systems that can be deformed into one another without going through a delocalization transition belong to the same localized phase [23-25, 29, 31, 49, 50]. Two Hamiltonians certainly belong to distinct localized phases when they support distinct eigenstate order [1215]-their eigenstates cannot be deformed into one another without their correlation length diverging. However, in this Letter, we argue eigenstate order is not the right concept to obtain a general understanding of localized topological phases. For one thing, not all ground state topological phases are compatible with MBL. More importantly, there exist anomalous localized topological phases 1 (ALT phases) [22-24, 31, 49, 50]. ALT phases have no eigenstate order, but nonetheless cannot be deformed to a trivial model-in which the commuting terms have disjoint support [51]-without going through a delocalization transition.

The existence of ALT phases in driven systems is recognized [22, 23, 26, 27, 31, 49, 52-54], but in static sys-

∗ dmlong@stanford.edu

1 The label 'anomalous' descends from Refs. [22, 24]. It is not related to other usages of the word 'anomalous' in the context of topological phases of matter.

This Letter systematizes the study of ALT phases in strongly interacting systems. Specifically, we focus on localized topological phases satisfying a condition we refer to as short-ranged entanglement . We will argue that instead of eigenstate order, the right way to characterize all topological localized phases with short-ranged entanglement is via quantum cellular automata (QCA) [56, 5860]-the locality preserving unitaries which map physical degrees of freedom to local integrals of motion (LIOMs) [Fig. 1(a)].

z

FIG. 1. (a) Many-body localization is characterized by commuting LIOMs, τ z j [6]. The LIOMs of an MBL model with short-ranged entanglement can be prepared from trivial, single-site LIOMs σ z j by a QCA, τ z j = V σ z j V † . ALT phases cannot be deformed back to a trivial model without delocalizing, even though they have no eigenstate order. (b) Multitone driven ALT phases correspond to loops or tori of QCA. The Ω -spectrum property QCA d -1 ≃ ΩQCA d (5) implies that n -tone-driven ALT phases in d dimensions correspond to ( d -1) -dimensional phases with n -1 drives. The translation QCA in d = 1 [58] corresponds with the anomalous Floquet phase in d = 2 [23], which in turn corresponds with a two-tone-driven phase in d = 3 with circulating surface states.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00825v2-figure-d23d9e4a-1b2f-47a0-bb5d-50e30a7f9274.jpg)

As QCA in spin systems have been classified in one [58] and two [61] dimensions, this correspondence allows for a classification of ALT phases in low dimensional systems. The same formalism extends to periodically (Floquet [62]) and quasiperiodically driven ALT phases [31], which are classified by loops or tori of QCA, respectively. Thus, we can treat static and n -tone-driven ALT phases in d dimensions on the same footing. We also considerably develop the study of the topology of QCA, which allows us to classify all ALT phases with d -n ≤ 2 [Fig. 1(b)].

Symmetry-enriched ALT phases (SALT phases) can also be classified by QCA, and we compute the classification for any symmetry group G and d -n ≤ 1 . Our methods naturally provide explicit soluble models of driven (S)ALT phases which, being finite depth circuits, are well suited to realization in quantum simulators [63-65].

ALT phases .-We aim to classify many-body localized phases. A Hamiltonian H is many-body localized when it possesses a complete set of commuting LIOMs τ z j such that [6]

$$H = \sum _ { j } h _ { j } \tau _ { j } ^ {$$

Here, j indexes LIOMs, ' · · · ' indicates higher-body interactions, h jk...l are real numbers which decay with the maximum distance between any of τ z j , . . . , τ z l , and we consider local degrees of freedom on the integer lattice in d dimensions, Z d .

Equation (1) identifies MBL Hamiltonians with commuting models, which are independently used as exactly solvable models of topological order [41-45] and symmetry-protected topological (SPT) phases [16, 46, 48, 66-70]. In the MBL context, the LIOMs may have exponential tails.

MBL arising from strong random disorder is believed to be asymptotically unstable for d &gt; 1 [71], and doubt has arisen regarding the LIOM description even for d = 1 [11, 72-77]. The case of correlated disorder remains open [78-82]. Regardless, the Hamiltonian (1) continues to function as an effective model for astronomical prethermal timescales [10, 11, 72, 76, 77, 83]. Indeed, the timescale for topological response in a localized model tends to be parametrically separate from the thermalization time [23, 25, 29-32, 50]. This justifies treating the topological classification of LIOM Hamiltonians (1) and the eventual thermalization process as separate problems. This Letter addresses the former.

We say two MBL Hamiltonians H 0 and H 1 are in the same localized topological phase if there is a continuous path of MBL Hamiltonians H t ( t ∈ [0 , 1] ) which connects them. Ignoring long-range resonances in the spectrum of H t , it is possible to construct a continuous path of locality preserving unitaries V t (that is, V AV t † t is localized near the support of A for all operators A ) such that V t conjugates the LIOMs of H 0 into LIOMs of H t [84-89]. Indeed, this unitary simply follows the adiabatic evolution of the many-body eigenstates. Including long-range resonances prevents an adiabatic construction of V t [90], but it is expected that sufficiently dilute resonances can be incorporated such that there still exists a V 1 that conjugates LIOMs of H 0 into those of H 1 and is continuously connected to the identity [9]. Thus, if H 0 and H 1 are in the same localized phase, their LIOMs are related by a continuous path of locality-preserving unitary operators. Locality-preserving unitaries are also called quantum cellular automata (QCA), and we will use this terminology in the rest of the paper.

An immediate consequence is that if H 0 and H 1 have distinct eigenstate order, they cannot be in the same localized phase. However, the converse is not true: even if they have the same eigenstate order, it need not be possible to deform H 0 to H 1 . This phenomenon is exemplified by ALT phases, which have no eigenstate order, but which cannot be deformed to a trivial MBL model with LIOMs supported on a single lattice site. We denote these trivial LIOMs σ z j , even when the local Hilbert space dimension is not 2.

An example of a nontrivial ALT phase is the threedimensional 3-fermion Walker-Wang model [47, 48, 56]. In the absence of symmetry, this model is believed to have no eigenstate order [48, 91], but nonetheless it supports topologically ordered surface states protected by localization. Phrased in the language of MBL, [56, Theorem II.4] showed that the 3-fermion model, and, in fact, any MBL Hamiltonian H with short-ranged entanglement -there are local unitaries τ x j which create and destroy excitations of τ z j without affecting other τ z k (Appendix A1)can be prepared from a trivial model by a QCA. That is, there is a QCA V such that H = V H trivial V † and H trivial has single-site LIOMs σ z j (Appendix A2). We note that the property we call short-range entanglement for MBL Hamiltonians implies all eigenstates are short-range entangled in the usual sense.

Hence, short-range entangled MBL models can be parameterized by QCA as V H trivial V † . Paths of nondegenerate MBL Hamiltonians induce paths of unitaries as described above, so Hamiltonians in the same phase can be prepared by QCA which are connected by a path. Thus, the path-connected components of QCA can be used to classify phases of MBL Hamiltonians. An additional subtlety is that this parameterization is not unique-any QCA V and V P which differ by a permutation of the sites P prepare identical LIOMs, and so represent Hamiltonians in the same phase (Appendix A 3).

We further allow stabilization by the stacking (tensor product) of additional ancillaeV 0 and V 0 ⊗ ✶ are said to be stably equivalent, and are considered to prepare Hamiltonians in the same phase [51, 58, 92, 93]. The space of stable equivalence classes of d -dimensional QCA is denoted QCA d , and its connected components π 0 (QCA ) d form an abelian group under tensor products of QCA [58, 94]. Short-range entangled localized phases are classified by the quotient of this group by permutation QCA.

The QCA formalism also lets us straightforwardly

identify which phases are anomalous. For any simultaneous eigenstate | 0 ⟩ of all σ z j , V | 0 ⟩ is an eigenstate of H . This defines a map e 0 ( V ) = V | 0 ⟩ from QCA d to the space of invertible states [95] (short-range entangled states) in d -dimensions, inv d . There is an induced homomorphism

$$e _ { 0 } ^ { * } \colon \pi _ { 0 } ( \text{$$

where addition of states is again defined by stacking. The crucial observation is that the homomorphism e ∗ 0 is, in general, neither surjective nor injective. The lack of surjectivity means that not all invertible states can arise as eigenstates of commuting models (for instance, chiral states [96, 97]). The lack of injectivity (even accounting for equivalence by permutations) precisely means that nontrivial ALT phases exist. Indeed, ALT phases may be identified as the nontrivial phases in the kernel of e ∗ 0 .

Most examples of ALT phases in the literature are driven systems [23-25, 29, 31, 49, 50]. We specialize to the case of multi-tone time dependence [28, 98, 99] (which includes Floquet driving as a special case), such that the Hamiltonian can be parameterized H t ( ) = H θ ( 1 ( ) t , . . . , θ n ( )) t , where θ k ( ) t = ω t k is a drive phase defined modulo 2 π . It is convenient to assemble the drive phases into a vector ⃗ θ t = ( θ 1 ( ) t , . . . , θ n ( )) t . Then we say that the system is MBL with short-ranged entanglement if it has a locality preserving generalized Floquet decomposition ( ℏ = 1 ) [28, 99]

$$U ( t ) = V ( \vec { \theta } _ { t } ) \exp \$$

Here, U t ( ) is the time evolution operator. Thus we have a correspondence between multi-tone driven localized phases and parameterized families of QCA, V ( ⃗ θ ) . That is, homotopy classes of continuous maps from the n -dimensional torus, T n := [ R / (2 π Z )] n , to QCA d , which we denote π 0 [ T n , QCA ] d . Ignoring weak invariants [100, 101] corresponding to low dimensional cross sections of T n , we may instead consider the n th homotopy group π n (QCA ) d .

Once more, there is an induced homomorphism e ∗ 0 : π n (QCA ) d → π n (inv d ) , and ALT phases are identified as the nontrivial elements of ker e ∗ 0 . As in the case of static systems, the parameterization by QCA is not unique, so that the classification of ALT phases is properly a quotient of ker e ∗ 0 . Further details regarding driven MBL phases are provided throughout Appendix A.

Classification .-The classification of static ALT phases in d dimensions is determined by the classification of QCA. That is, by π 0 (QCA ) d . There has been great recent progress in classifying QCA in spin systems-a complete classification now exists in one [58] and two [61] dimensions, and several nontrivial phases are known in three dimensions [55-57]. The classification of ALT phases inherits this progress.

For Floquet systems, there is another correspondence known between localized phases and QCA. Reference [23] showed that a periodically driven localized phase can be

FIG. 2. The swindle map S : QCA d -1 → ΩQCA d maps a QCA v to a loop of QCA S v, t ( ) in one higher dimension. S v, t ( ) is a finite depth circuit built from v , v † (squares), and a path from the identity to a swap gate, SWAP( ) t (line with crosses).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00825v2-figure-44630102-139e-4cb3-bc03-a5a1a6b150ac.jpg)

classified by its action on the edge when placed on open boundary conditions. This action is described by a QCA of one dimension lower, π 0 (QCA d -1 ) . On the other hand, both our construction and an equivalent formulation of unitary loops [49] show that Floquet localized phases should be classified by π 1 (QCA ) d . We deduce that these groups must be isomorphic,

$$\begin{array} {$$

The relationship Eq. (4) is reminiscent of a widely held conjecture that invertible states form an Ω -spectrum [95, 96, 102-104]. That is, there is a homotopy equivalence [105] inv d -1 ≃ Ωinv d , where Ω X is the space of loops that begin and end at a specified base point in the space X . In particular π n -1 (inv d -1 ) = ∼ π n (inv d ) for all n . We make the conjecture that QCA also forms an Ω -spectrum,

$$\underset { \begin{array} { c } \omega \\ \text{local} \$$

The base point for the loops in ΩQCA d is the identity.

Equation (5) is interpreted through a picture of pumping [95, 103, 104]. Taking a loop γ ∈ ΩQCA d and restricting it to open boundary conditions, the action in the bulk must still be trivial at the end of the loop, due to the finite range of the QCA. However, on the ( d -1) -dimensional edge, there may be a nontrivial action by some v ∈ QCA d -1 [23]. We say that γ pumps v across the system, and Eq. (5) hypothesizes that the homotopy structure of the loop space ΩQCA d is completely characterized by the pumped QCA v .

More formally, we argue that the swindle map S : QCA d -1 → ΩQCA d is a homotopy equivalence. Given some v ∈ QCA d -1 , the swindle S v, t ( ) is defined through a finite depth quantum circuit [59, 60, 95, 106] (Fig. 2). Let swap jk ( ) t be a path from the identity to the swap gate in the unitary group on sites j and k . We stack these gates together in a codimension 1 plane,

$$\text{SWAP} _ { x, x + 1 } ( t ) = \bigotimes _ { \vec { y } \in \mathbb { Z } ^ { d - 1 } } \text{swap} _ { ( x, \vec { y } ), ( x + 1, \vec { y } ) } ( t ). \quad ( 6 ) \\ \text{...} \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

Then the path of QCA on a two-layer system

$$\stackrel { \iota - \mu } { \text{ be } } \quad u _ { x, x + 1 } ( v, t ) \coloneqq v _ { x } S W A P _ { x, x + 1 } ( t ) v _ { x } ^ { \dagger } S W A P _ { x, x + 1 } ( t ) ^ { \dagger } \ ( 7 )$$

TABLE I. The homotopy groups π n (QCA ) =: d h q are a function of q = d -n , as is π n (inv d ) =: h q inv . For spin systems h 3 is conjectured to be the Witt group of modular tensor categories [55-57, 108, 109]; h 2 is trivial [61]; h 1 is isomorphic to log Q × [58]; and h q for q ≤ 0 is trivial for even q , Q Z / for q = -1 , and Q for odd q ≤ -3 (Appendix C). In all groups in the table, the group operation is addition.

| q       | 3   | 2   | 1       |   0 | - 1   | - 2   | - 3   |   - 4 | - 5   |
|---------|-----|-----|---------|-----|-------|-------|-------|-------|-------|
| h q     | ?   | 0   | log Q × |   0 | Q / Z | 0     | Q     |     0 | Q     |
| h q inv | 0   | Z   | 0       |   0 | 0     | Z     | 0     |     0 | 0     |

TABLE II. The classification of strong invariants in n -tone driven d -dimensional localized phases of spins with short-range entanglement. The classifying group differs from π n (QCA ) d for ( d, n ) ∈ { (0 1) (1 0) , , , } due to the nonuniqueness of the parameterization by QCA. All nontrivial elements of this table correspond to ALT phases.

|   d n | 0   | 1       | 2       | 3     | 4     | 5   |
|-------|-----|---------|---------|-------|-------|-----|
|     0 | 0   | 0       | 0       | Q     | 0     | Q   |
|     1 | 0   | 0       | Q / Z   | 0     | Q     | 0   |
|     2 | 0   | log Q × | 0       | Q / Z | 0     | Q   |
|     3 | ?   | 0       | log Q × | 0     | Q / Z | 0   |

interpolates between the identity and v x v † x +1 , where v x acts as v on the ( d -1) -dimensional surface { x } × Z d -1 . Finally, the swindle map is

$$S ( v, t ) = \left ( \prod _ { x \in 2 \mathbb { Z } + 1 } u _ { x, x + 1 } ( v, t ) \right ) \left ( \prod _ { x \in 2 \mathbb { Z } } u _ { x, x + 1 } ( v, t ) \right ), \quad \begin{array} { c } \text{incide} \\ \text{$a lgeb$} \\ \text{To} \\ \text{specific} \end{array} \end{array}$$

which is continuous in t and v and equal to the identity for ( t, v ) ∈ { 0 1 , } × QCA d -1 (Fig. 2).

Assuming that any loop of QCA can be restricted to open boundary conditions such that it remains a QCA, we show that S is a homotopy equivalence in Appendix B. Formal details are deferred to Ref. [107].

Iteratively applying Eq. (5), we reduce the classification of driven localized phases in any dimension to either the classification of static QCA, or a driven classification of zero-dimensional QCA,

$$\pi _ { n } ( Q C A _ { d } ) \cong \begin{cases} \pi _ { 0 } ( Q C A _ { d - n } ) & \text{for $d\geq n$}, \quad ( 9 ) & \text{other} \\ \pi _ { n - d } ( Q C A _ { 0 } ) & \text{for $d<n$}. \quad ( 9 ) & \text{tary} \\ & \text{of pa} \end{cases}$$

The space QCA 0 is a stabilization of the projective unitary group by stacking ancillae. We find its homotopy groups in Appendix C. Using the known classification of one- and two-dimensional QCA, we have all the homotopy groups π n (QCA ) d for spin systems with d -n ≤ 2 . These are listed in Tab. I, and the corresponding classification of localized phases (excluding permutation QCA) is listed in Tab. II. Comparing to the classification of invertible states shows that all the nontrivial phases in the tables are ALT phases.

SALT phases .-Two systems belong to the same symmetric localized phase if they can be deformed into one another without breaking some symmetry or delocalizing. If the short-ranged entangled MBL Hamiltonian H commutes with an on-site unitary abelian symmetry, then it can still have stable and complete LIOMs [46, 54, 67, 110], and the correspondence with QCA continues to hold. Thus, symmetry-enriched ALT phases (SALT phases, symmetric localized phases without eigenstate order) can be classified by symmetric QCA. (We do not consider weak symmetries, where the symmetry only exists on average [111, 112].) In fact, at the level of QCA, there is no need to restrict to abelian groups.

The conjectured classification of symmetry-enriched invertible states can be stated in terms of just the space inv d and the abstract symmetry group G [66, 95, 96, 102104, 113]. However, the classification of symmetric QCA depends not just on G , but on the specific representation ρ of the symmetry action [93, 114]. We compute the classification of symmetric QCA, and thus SALT phases, with d -n ≤ 1 and any on-site unitary representation of any group in Appendix D. The classification of symmetric QCA closely mirrors Tab. I, except that it includes additional subgroups in d -n ≥ 0 [given by H d -n +1 ( G U ; (1)) for finite G ] related to nontrivial SPT states or pumps thereof; moreover, rational numbers p/q ∈ Q in Tab. I are replaced with formal ratios of representations, ρ/σ . Addition and multiplication for ρ/σ is defined through the same rules as for ratios of integers, but replacing ' + ' with ' ⊕ ' and ' × ' with ' ⊗ '. In an appropriate coincidence, this mathematical structure is known as an algebraic localization .

To define the algebraic localization precisely, we must specify the set of ancillae representations,

$$\varphi \quad \mathcal { R } = \{ \rho \in \text{Rep} ( G ) \, \colon \, \text{some ancillae carry } \rho \}, \quad ( 1 0 )$$

which is closed under tensor product-if ρ 1 and ρ 2 can be added then ρ 1 ⊗ ρ 2 corresponds to adding them both. The algebraic localization with denominators in R is the quotient

$$\text{her} _ { \text{ion} } & & \mathcal { R } ^ { - 1 } R _ { G } = ( R _ { G } \times \mathcal { R } ) / \sim \quad \text{where} \ ( \rho _ { 1 }, \sigma _ { 1 } ) \sim ( \rho _ { 2 }, \sigma _ { 2 } ) \\ & \quad \Leftrightarrow \exists \rho \in \mathcal { R } \ \text{such that } \rho \otimes \rho _ { 1 } \otimes \sigma _ { 2 } \simeq \rho \otimes \rho _ { 2 } \otimes \sigma _ { 1 }, \quad ( 1 1 )$$

where R G is the representation ring of G , and ' ≃ ' is unitary equivalence of representations. Equivalence classes of pairs ( ρ, σ ) are denoted ρ/σ .

The representation theory of any group which is of relevance to physics is understood in extraordinary detail, so calculating R -1 R G , and thus the classification of lowdimensional SALT phases, becomes straightforward. For instance, the classification of ( d, n ) = (2 1) , U (1) -SALT phases in Ref. [26] by rational polynomials follows from observing that R U (1) ∼ = Z [ x, x -1 ] is isomorphic to the ring of integer Laurent polynomials. The fact that including the regular representation, ρ reg , for a finite group in R leads to a collapse in the classification [93] follows

FIG. 3. Sketch of the pulse schedule for the models defined by the swindle map, with resulting unitary evolution ˜ S . Sites are coupled along even bonds (red, t ∈ [0 , 1 3 ) ) and then along odd bonds (blue, t ∈ [ 1 3 , 2 3 ) ). Finally, an on-site disorder potential may be applied to promote localization (gray, t ∈ [ 2 3 , 1] ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00825v2-figure-4d4edb59-9f58-4f05-8474-c4c609e32990.jpg)

from the fact that ρ reg is a zero divisor in R G . This means Eq. (11) can be satisfied even if ρ 1 ⊗ σ 2 ≃ ρ 2 ⊗ σ 1 .

̸

Models .-While the swindle map [Eq. (8), Fig. 2] was introduced as an abstract tool, it also defines a class of soluble models of n -tone-driven (S)ALT phases given any ( n -1) -parameter family of QCA in one dimension lower. As these models are finite depth circuits, they are also well suited for realization in quantum simulator experiments [63-65] (Fig. 3).

As an example, Fig. 1(b) illustrates a threedimensional phase with two driving tones which pumps anomalous Floquet phases to its boundary, and thus has a quantized surface magnetization [27, 52]. This is obtained by applying the swindle map to any model of the anomalous Floquet phase [22, 23].

As three-dimensional systems are typically still challenging for quantum simulators, in Appendix E we construct a soluble model of the quasiperiodic energy pump (QP pump) [31, 32, 53], a one-dimensional two-tonedriven SALT phase with U (1) symmetry. The QP pump supports edge states which pump energy between the two driving fields. The classification of SALT phases shows that the QP pump is stable to strong interactions, provides a nearest-neighbor qubit model for this phase, and indicates that forgetting the U (1) symmetry puts the model in the phase of Refs. [29, 30].

Discussion .-In the context of exactly soluble com-

- [1] P. W. Anderson, Absence of diffusion in certain random lattices, Physical review 386 , 1492 (1958).
- [2] I. V. Gornyi, A. D. Mirlin, and D. G. Polyakov, Interacting electrons in disordered wires: Anderson localization and lowt transport, Phys. Rev. Lett. 95 , 206603 (2005).
- [3] D. M. Basko, I. L. Aleiner, and B. L. Altshuler, Metal-insulator transition in a weakly interacting manyelectron system with localized single-particle states, Annals of Physics 321 , 1126 (2006), arXiv:0506617 [condmat].
- [4] V. Oganesyan and D. A. Huse, Localization of inter-

muting models and driven phases of localized systems, authors have repeatedly expressed confusion regarding the position of their models within existing classifications [23, 26, 29, 91]. Most of these are classifications of states [28, 66, 115], while localized models (including any commuting projector model) impose nontrivial structure on the entire spectrum. This structure need not be captured by any single state, making such classifications of limited applicability. This point is implicit in the literature [23, 31, 50, 91], but the lack of a systematic and general study seems to have allowed confusion to persist.

Indeed, both the study of ground state phases and localized phases are rich mathematical problems with important implications for disparate physical systems. While the two problems are related (by e ∗ 0 , with shortranged entanglement), neither can be regarded as a special case of the other.

The general relationship between (S)ALT phases and QCA we have demonstrated, as well as the Ω -spectrum property which relates static (S)ALT phases to driven phases, extends to fermionic systems essentially without modification (Appendix A-B) [56]. The Ω -spectrum structure of the classification can also be exploited to classify SALT phases with lattice symmetries [103, 104, 116]. The more challenging direction of generalization is to include topological order (long-ranged entanglement), where the correspondence with QCA breaks down [20]. This case is of particular interest given recent proposals of using driven topologically ordered localized models to perform quantum computation [117-120]. The notion of measurement QCA introduced by Ref. [121] may provide a route towards classifying localized phases with topological order, which we will explore in future work.

The authors are particularly grateful to R. Kobayashi for collaboration during the early stages of this project. The authors also thank M. Barkeshli, A. Chandran, L. Fidkowski, S. Pace, R. Thorngren, Y.-T. Tu, B. Ware, and C. Zhang for useful discussions. This work is supported by the Laboratory for Physical Sciences (DML). Research at Perimeter Institute is supported in part by the Government of Canada through the Department of Innovation, Science and Economic Development, and by the Province of Ontario through the Ministry of Colleges and Universities (DVE).

acting fermions at high temperature, Phys. Rev. B 75 , 155111 (2007).

- [5] A. Pal and D. A. Huse, Many-body localization phase transition, Phys. Rev. B 82 , 174411 (2010).
- [6] D. A. Huse, R. Nandkishore, and V. Oganesyan, Phenomenology of fully many-body-localized systems, Phys. Rev. B 90 , 174202 (2014).
- [7] M. Schreiber, S. S. Hodgman, P. Bordia, H. P. Lüschen, M. H. Fischer, R. Vosk, E. Altman, U. Schneider, and I. Bloch, Observation of many-body localization of interacting fermions in a quasirandom optical lattice, Science

## 349 , 842 (2015), arXiv:1501.05661.

- [8] J. Smith, A. Lee, P. Richerme, B. Neyenhuis, P. W. Hess, P. Hauke, M. Heyl, D. A. Huse, and C. Monroe, Many-body localization in a quantum simulator with programmable random disorder, Nature Physics 12 , 907 (2016).
- [9] J. Z. Imbrie, On Many-Body Localization for Quantum Spin Chains, Journal of Statistical Physics 163 , 998 (2016), arXiv:1403.7837.
- [10] W. D. Roeck, F. Huveneers, B. Meeus, and O. A. Prośniak, Rigorous and simple results on very slow thermalization, or quasi-localization, of the disordered quantum chain (2023), arXiv:2308.08462 [math-ph].
- [11] P. Sierant, M. Lewenstein, A. Scardicchio, L. Vidmar, and J. Zakrzewski, Many-body localization in the age of classical computing (2024), arXiv:2403.07111 [condmat.dis-nn].
- [12] D. A. Huse, R. Nandkishore, V. Oganesyan, A. Pal, and S. L. Sondhi, Localization-protected quantum order, Physical Review B - Condensed Matter and Materials Physics 88 , 1 (2013), arXiv:1304.1158.
- [13] S. A. Parameswaran and R. Vasseur, Many-body localization, symmetry and topology, Reports on Progress in Physics 81 , 082501 (2018).
- [14] D. Pekker, G. Refael, E. Altman, E. Demler, and V. Oganesyan, Hilbert-glass transition: New universality of temperature-tuned many-body dynamical quantum criticality, Phys. Rev. X 4 , 011052 (2014).
- [15] J. A. Kjäll, J. H. Bardarson, and F. Pollmann, Manybody localization in a disordered quantum ising chain, Phys. Rev. Lett. 113 , 107204 (2014).
- [16] X.-G. Wen, Colloquium: Zoo of quantum-topological phases of matter, Rev. Mod. Phys. 89 , 041004 (2017).
- [17] B. Bauer and C. Nayak, Area laws in a many-body localized state and its implications for topological order, Journal of Statistical Mechanics: Theory and Experiment 2013 , P09005 (2013).
- [18] K. Slagle, Z. Bi, Y.-Z. You, and C. Xu, Many-body localization of symmetry protected topological states (2015), arXiv:1505.05147 [cond-mat.str-el].
- [19] A. Chan and T. B. Wahl, Classification of symmetryprotected topological many-body localized phases in one dimension, Journal of Physics: Condensed Matter 32 , 305601 (2020).
- [20] T. B. Wahl and B. Béri, Local integrals of motion for topologically ordered many-body localized systems, Phys. Rev. Res. 2 , 033099 (2020).
- [21] F. Venn, T. B. Wahl, and B. Béri, Many-bodylocalization protection of eigenstate topological order in two dimensions (2022), arXiv:2212.09775 [cond-mat.disnn].
- [22] M. S. Rudner, N. H. Lindner, E. Berg, and M. Levin, Anomalous edge states and the bulk-edge correspondence for periodically driven two-dimensional systems, Phys. Rev. X 3 , 031005 (2013).
- [23] H. C. Po, L. Fidkowski, T. Morimoto, A. C. Potter, and A. Vishwanath, Chiral floquet phases of many-body localized bosons, Phys. Rev. X 6 , 041070 (2016).
- [24] P. Titum, E. Berg, M. S. Rudner, G. Refael, and N. H. Lindner, Anomalous floquet-anderson insulator as a nonadiabatic quantized charge pump, Phys. Rev. X 6 , 021013 (2016).
- [25] F. Nathan, D. Abanin, E. Berg, N. H. Lindner, and M. S. Rudner, Anomalous floquet insulators, Phys. Rev.
- B 99 , 195133 (2019).
- [26] C. Zhang and M. Levin, Classification of interacting floquet phases with u (1) symmetry in two dimensions, Phys. Rev. B 103 , 064302 (2021).
- [27] F. Nathan, D. A. Abanin, N. H. Lindner, E. Berg, and M. S. Rudner, Hierarchy of many-body invariants and quantized magnetization in anomalous floquet insulators, arXiv:1907.12228 (2019).
- [28] D. V. Else, W. W. Ho, and P. T. Dumitrescu, Longlived interacting phases of matter protected by multiple time-translation symmetries in quasiperiodically driven systems, Phys. Rev. X 10 , 021032 (2020).
- [29] A. J. Friedman, B. Ware, R. Vasseur, and A. C. Potter, Topological edge modes without symmetry in quasiperiodically driven spin chains, Phys. Rev. B 105 , 115117 (2022).
- [30] P. T. Dumitrescu, J. G. Bohnet, J. P. Gaebler, A. Hankin, D. Hayes, A. Kumar, B. Neyenhuis, R. Vasseur, and A. C. Potter, Dynamical topological phase realized in a trapped-ion quantum simulator, Nature 607 , 463 (2022).
- [31] D. M. Long, P. J. D. Crowley, and A. Chandran, Nonadiabatic topological energy pumps with quasiperiodic driving, Phys. Rev. Lett. 126 , 106805 (2021).
- [32] F. Nathan, R. Ge, S. Gazit, M. Rudner, and M. Kolodrubetz, Quasiperiodic floquet-thouless energy pump, Phys. Rev. Lett. 127 , 166804 (2021).
- [33] K. Sacha, Modeling spontaneous breaking of timetranslation symmetry, Phys. Rev. A 91 , 033617 (2015).
- [34] V. Khemani, A. Lazarides, R. Moessner, and S. L. Sondhi, Phase structure of driven quantum systems, Phys. Rev. Lett. 116 , 250401 (2016).
- [35] D. V. Else, B. Bauer, and C. Nayak, Floquet time crystals, Phys. Rev. Lett. 117 , 090402 (2016).
- [36] C. W. von Keyserlingk, V. Khemani, and S. L. Sondhi, Absolute stability and spatiotemporal long-range order in floquet systems, Phys. Rev. B 94 , 085112 (2016).
- [37] S. Choi, J. Choi, R. Landig, G. Kucsko, H. Zhou, J. Isoya, F. Jelezko, S. Onoda, H. Sumiya, V. Khemani, C. von Keyserlingk, N. Y. Yao, E. Demler, and M. D. Lukin, Observation of discrete time-crystalline order in a disordered dipolar many-body system, Nature 543 , 221 (2017).
- [38] V. Khemani, R. Moessner, and S. L. Sondhi, A brief history of time crystals (2019), arXiv:1910.10745 [condmat.str-el].
- [39] D. V. Else, C. Monroe, C. Nayak, and N. Y. Yao, Discrete time crystals, Annual Review of Condensed Matter Physics 11 , 467 (2020).
- [40] X. Mi et al. , Time-crystalline eigenstate order on a quantum processor, Nature 601 , 531 (2022).
- [41] A. Y. Kitaev, Fault-tolerant quantum computation by anyons, Annals of Physics 303 , 2 (2003).
- [42] M. A. Levin and X.-G. Wen, String-net condensation: A physical mechanism for topological phases, Phys. Rev. B 71 , 045110 (2005).
- [43] C. Chamon, Quantum glassiness in strongly correlated clean systems: An example of topological overprotection, Phys. Rev. Lett. 94 , 040402 (2005).
- [44] S. Bravyi, B. Leemhuis, and B. M. Terhal, Topological order in an exactly solvable 3d spin model, Annals of Physics 326 , 839 (2011).
- [45] J. Haah, Local stabilizer codes in three dimensions without string logical operators, Phys. Rev. A 83 , 042330

(2011).

- [46] M. Levin and Z.-C. Gu, Braiding statistics approach to symmetry-protected topological phases, Phys. Rev. B 86 , 115109 (2012).
- [47] K. Walker and Z. Wang, (3+1)-TQFTs and topological insulators, Frontiers of Physics 7 , 150 (2012).
- [48] F. J. Burnell, X. Chen, L. Fidkowski, and A. Vishwanath, Exactly soluble model of a three-dimensional symmetry-protected topological phase of bosons with surface topological order, Physical Review B Condensed Matter and Materials Physics 90 , 1 (2014), arXiv:1302.7072.
- [49] R. Roy and F. Harper, Floquet topological phases with symmetry in all dimensions, Phys. Rev. B 95 , 195128 (2017).
- [50] B. Lapierre, T. Neupert, and L. Trifunovic, Topologically localized insulators, Phys. Rev. Lett. 129 , 256401 (2022).
- [51] M. B. Hastings, Classifying quantum phases with the kirby torus trick, Phys. Rev. B 88 , 165114 (2013).
- [52] F. Nathan, M. S. Rudner, N. H. Lindner, E. Berg, and G. Refael, Quantized magnetization density in periodically driven systems, Phys. Rev. Lett. 119 , 186801 (2017).
- [53] D. M. Long, P. J. D. Crowley, and A. Chandran, Coupled layer construction for synthetic Hall effects in driven systems, Phys. Rev. B 106 , 144203 (2022).
- [54] C. Zhang, Topological invariants for SPT entanglers (2022), arXiv:2210.02485 [cond-mat.str-el].
- [55] W. Shirley, Y.-A. Chen, A. Dua, T. D. Ellison, N. Tantivasadakarn, and D. J. Williamson, Three-dimensional quantum cellular automata from chiral semion surface topological order and beyond, PRX Quantum 3 , 030326 (2022).
- [56] J. Haah, L. Fidkowski, and M. B. Hastings, Nontrivial Quantum Cellular Automata in Higher Dimensions, Communications in Mathematical Physics 398 , 469 (2023).
- [57] L. Fidkowski and M. B. Hastings, Pumping chirality in three dimensions (2023), arXiv:2309.15903 [condmat.str-el].
- [58] D. Gross, V. Nesme, H. Vogts, and R. F. Werner, Index Theory of One Dimensional Quantum Walks and Cellular Automata, Communications in Mathematical Physics 310 , 419 (2012).
- [59] T. Farrelly, A review of Quantum Cellular Automata, Quantum 4 , 368 (2020).
- [60] P. Arrighi, An overview of quantum cellular automata, Natural Computing 18 , 885 (2019).
- [61] M. Freedman and M. B. Hastings, Classification of Quantum Cellular Automata, Communications in Mathematical Physics 376 , 1171 (2020).
- [62] G. Floquet, Sur les équations différentielles linéaires à coefficients périodiques, Ann. Sci. École Norm. Sup. 12 , 47 (1883).
- [63] A. Eckardt, Colloquium: Atomic quantum gases in periodically driven optical lattices, Rev. Mod. Phys. 89 , 011004 (2017).
- [64] C. Gross and I. Bloch, Quantum simulations with ultracold atoms in optical lattices, Science 357 , 995 (2017).
- [65] M. Kjaergaard, M. E. Schwartz, J. Braumüller, P. Krantz, J. I.-J. Wang, S. Gustavsson, and W. D. Oliver, Superconducting qubits: Current state of play, Ann. Rev. Condens. Matter Phys. 11 , 369 (2020).
- [66] X. Chen, Z.-C. Gu, and X.-G. Wen, Classification of gapped symmetric phases in one-dimensional spin systems, Phys. Rev. B 83 , 035107 (2011).
- [67] X. Chen, Z.-X. Liu, and X.-G. Wen, Two-dimensional symmetry-protected topological orders and their protected gapless edge excitations, Phys. Rev. B 84 , 235141 (2011).
- [68] X. Chen, Z.-C. Gu, Z.-X. Liu, and X.-G. Wen, Symmetry-protected topological orders in interacting bosonic systems, Science 338 , 1604 (2012).
- [69] X. Chen, Z.-C. Gu, Z.-X. Liu, and X.-G. Wen, Symmetry protected topological orders and the group cohomology of their symmetry group, Phys. Rev. B 87 , 155114 (2013).
- [70] X. Chen, Y.-M. Lu, and A. Vishwanath, Symmetryprotected topological phases from decorated domain walls, Nature Communications 5 , 3507 (2014).
- [71] W. De Roeck and F. Huveneers, Stability and instability towards delocalization in many-body localization systems, Phys. Rev. B 95 , 155129 (2017).
- [72] J. Šuntajs, J. Bonča, T. c. v. Prosen, and L. Vidmar, Quantum chaos challenges many-body localization, Phys. Rev. E 102 , 062144 (2020).
- [73] M. Schulz, S. R. Taylor, A. Scardicchio, and M. Žnidarič, Phenomenology of anomalous transport in disordered one-dimensional systems, Journal of Statistical Mechanics: Theory and Experiment 2020 , 10.1088/1742-5468/ab6de0 (2020), arXiv:1909.09507.
- [74] D. Sels and A. Polkovnikov, Dynamical obstruction to localization in a disordered spin chain, Phys. Rev. E 104 , 054105 (2021).
- [75] D. Sels and A. Polkovnikov, Thermalization of dilute impurities in one-dimensional spin chains, Phys. Rev. X 13 , 011041 (2023).
- [76] A. Morningstar, L. Colmenarez, V. Khemani, D. J. Luitz, and D. A. Huse, Avalanches and many-body resonances in many-body localized systems, Phys. Rev. B 105 , 174205 (2022).
- [77] D. Sels, Bath-induced delocalization in interacting disordered spin chains, Phys. Rev. B 106 , L020202 (2022).
- [78] V. Khemani, D. N. Sheng, and D. A. Huse, Two universality classes for the many-body localization transition, Phys. Rev. Lett. 119 , 075702 (2017).
- [79] U. Agrawal, S. Gopalakrishnan, and R. Vasseur, Universality and quantum criticality in quasiperiodic spin chains, Nature Communications 11 , 2225 (2020).
- [80] U. Agrawal, R. Vasseur, and S. Gopalakrishnan, Quasiperiodic many-body localization transition in dimension d ≥ 1 , Phys. Rev. B 106 , 094206 (2022).
- [81] Y.-T. Tu, D. Vu, and S. Das Sarma, Avalanche stability transition in interacting quasiperiodic systems, Phys. Rev. B 107 , 014203 (2023).
- [82] Y.-T. Tu, D. M. Long, and S. D. Sarma, Interacting quasiperiodic spin chains in the prethermal regime (2024), arXiv:2405.01622 [cond-mat.dis-nn].
- [83] D. M. Long, P. J. D. Crowley, V. Khemani, and A. Chandran, Phenomenology of the prethermal manybody localized regime, Phys. Rev. Lett. 131 , 106301 (2023).
- [84] M. V. Berry and J. M. Robbins, Chaotic classical and half-classical adiabatic reactions: geometric magnetism and deterministic friction, Proc. R. Soc. Lond. A 442 , 659-672 (1993).
- [85] C. Jarzynski, Geometric phases and anholonomy for a

- class of chaotic classical systems, Phys. Rev. Lett. 74 , 1732 (1995).
- [86] M. B. Hastings, Locality in quantum systems (2010), arXiv:1008.5137 [math-ph].
- [87] P. W. Claeys, M. Pandey, D. Sels, and A. Polkovnikov, Floquet-engineering counterdiabatic protocols in quantum many-body systems, Phys. Rev. Lett. 123 , 090602 (2019).
- [88] M. Pandey, P. W. Claeys, D. K. Campbell, A. Polkovnikov, and D. Sels, Adiabatic eigenstate deformations as a sensitive probe for quantum chaos, Phys. Rev. X 10 , 041017 (2020).
- [89] S. Gopalakrishnan, M. Müller, V. Khemani, M. Knap, E. Demler, and D. A. Huse, Low-frequency conductivity in many-body localized systems, Phys. Rev. B 92 , 104202 (2015).
- [90] V. Khemani, R. Nandkishore, and S. L. Sondhi, Nonlocal adiabatic response of a localized system to local manipulations, Nature Physics 11 , 560 (2015).
- [91] L. Fidkowski, J. Haah, and M. B. Hastings, Exactly solvable model for a 4 + 1D beyond-cohomology symmetry-protected topological phase, Phys. Rev. B 101 , 155124 (2020).
- [92] L. Piroli and J. I. Cirac, Quantum cellular automata, tensor networks, and area laws, Phys. Rev. Lett. 125 , 190402 (2020).
- [93] Z. Gong, C. Sünderhauf, N. Schuch, and J. I. Cirac, Classification of matrix-product unitaries with symmetries, Phys. Rev. Lett. 124 , 100402 (2020).
- [94] M. Freedman, J. Haah, and M. B. Hastings, The Group Structure of Quantum Cellular Automata, Communications in Mathematical Physics 389 , 1277 (2022).
- [95] A. Kitaev, On the classification of short-range entangled states (2013), talk at the Simons Center.
- [96] A. Kitaev, Anyons in an exactly solved model and beyond, Annals of Physics 321 , 2 (2006), january Special Issue.
- [97] A. Kapustin and L. Fidkowski, Local Commuting Projector Hamiltonians and the Quantum Hall Effect, Communications in Mathematical Physics 373 , 763 (2020).
- [98] T.-S. Ho, S.-I. Chu, and J. V. Tietz, Semiclassical manymode floquet theory, Chem. Phys. Lett. 96 , 464 (1983).
- [99] D. M. Long, P. J. D. Crowley, and A. Chandran, Manybody localization with quasiperiodic driving, Phys. Rev. B 105 , 144204 (2022).
- [100] A. Kitaev, Periodic table for topological insulators and superconductors, AIP Conf. Proc. 1134 , 22 (2009).
- [101] R. Roy and F. Harper, Periodic table for floquet topological insulators, Phys. Rev. B 96 , 155118 (2017).
- [102] A. Kitaev, Homotopy-theoretic approach to SPT phases in action: Z 16 classification of three-dimensional superconductors (2015), talk at the Symmetry and Topology in Quantum Matter Workshop, Institute for Pure and Applied Mathematics, University of California, Los Angeles, California.
- [103] C. Z. Xiong, Minimalist approach to the classification of symmetry protected topological phases, Journal of Physics A: Mathematical and Theoretical 51 , 445001 (2018).
- [104] D. Gaiotto and T. Johnson-Freyd, Symmetry protected topological phases and generalized cohomology, Journal of High Energy Physics 2019 , 7 (2019).
- [105] A. Hatcher, Algebraic Topology , Algebraic Topology (Cambridge University Press, 2002).
- [106] P. Arrighi, V. Nesme, and R. Werner, Unitarity plus causality implies localizability, Journal of Computer and System Sciences 77 , 372 (2011), adaptivity in Heterogeneous Environments.
- [107] D. M. Long and D. V. Else, Homotopy theory of quantum cellular automata, in preparation.
- [108] J. Haah, Clifford quantum cellular automata: Trivial group in 2D and Witt group in 3D, Journal of Mathematical Physics 62 , 092202 (2021).
- [109] J. Haah, Invertible Subalgebras, Communications in Mathematical Physics 403 , 661 (2023).
- [110] I. V. Protopopov, W. W. Ho, and D. A. Abanin, Effect of su(2) symmetry on many-body localization and thermalization, Phys. Rev. B 96 , 041122 (2017).
- [111] C. de Groot, A. Turzillo, and N. Schuch, Symmetry Protected Topological Order in Open Quantum Systems, Quantum 6 , 856 (2022).
- [112] R. Ma and C. Wang, Average symmetry-protected topological phases, Phys. Rev. X 13 , 031016 (2023).
- [113] D. V. Else and C. Nayak, Classifying symmetryprotected topological phases through the anomalous action of the symmetry on the edge, Phys. Rev. B 90 , 235137 (2014).
- [114] C. Zhang, Note on quantum cellular automata and strong equivalence (2023), arXiv:2306.03171 [quant-ph].
- [115] D. V. Else and C. Nayak, Classification of topological phases in periodically driven interacting systems, Phys. Rev. B 93 , 201103 (2016).
- [116] R. Thorngren and D. V. Else, Gauging spatial symmetries and the classification of topological crystalline phases, Phys. Rev. X 8 , 011040 (2018).
- [117] H. C. Po, L. Fidkowski, A. Vishwanath, and A. C. Potter, Radical chiral floquet phases in a periodically driven kitaev model and beyond, Phys. Rev. B 96 , 245116 (2017).
- [118] M. B. Hastings and J. Haah, Dynamically Generated Logical Qubits, Quantum 5 , 564 (2021).
- [119] D. Aasen, Z. Wang, and M. B. Hastings, Adiabatic paths of hamiltonians, symmetries of topological order, and automorphism codes, Phys. Rev. B 106 , 085122 (2022).
- [120] A. Paetznick, C. Knapp, N. Delfosse, B. Bauer, J. Haah, M. B. Hastings, and M. P. da Silva, Performance of planar floquet codes with majorana-based qubits, PRX Quantum 4 , 010310 (2023).
- [121] D. Aasen, J. Haah, Z. Li, and R. S. K. Mong, Measurement quantum cellular automata and anomalies in floquet codes (2023), arXiv:2304.01277 [quant-ph].
- [122] P. Hall, On representatives of subsets, Journal of the London Mathematical Society s1-10 , 26 (1935).
- [123] M. Hall, Combinatorial Theory , Wiley Classics Library (Wiley, 1998).
- [124] H. Sambe, Steady states and quasienergies of a quantum-mechanical system in an oscillating field, Phys. Rev. A 7 , 2203 (1973).
- [125] P. M. Blekher, H. R. Jauslin, and J. L. Lebowitz, Floquet spectrum for two-level systems in quasiperiodic time-dependent fields, J. Stat. Phys. 68 , 271 (1992).
- [126] A. Verdeny, J. Puig, and F. Mintert, Quasi-periodically driven quantum systems, Z. Naturforsch., A: J. Phys. Sci. 71 , 897 (2016), arXiv:1603.03923.
- [127] I. Martin, G. Refael, and B. Halperin, Topological frequency conversion in strongly driven quantum systems, Phys. Rev. X 7 , 041008 (2017).
- [128] P. J. D. Crowley, I. Martin, and A. Chandran, Topolog-

ical classification of quasiperiodically driven quantum systems, Phys. Rev. B 99 , 064306 (2019).

- [129] A. Y. Kitaev, Unpaired majorana fermions in quantum wires, Physics-Uspekhi 44 , 131 (2001).
- [130] S. Yao, Z. Yan, and Z. Wang, Topological invariants of floquet systems: General formulation, special properties, and floquet topological defects, Phys. Rev. B 96 , 195303 (2017).
- [131] J. P. C. Greenlees and J. P. May, Chapter 8 - equivariant stable homotopy theory, in Handbook of Algebraic Topology , edited by I. M. James (North-Holland, Amsterdam, 1995) pp. 277-323.
- [132] J. I. Cirac, D. Perez-Garcia, N. Schuch, and F. Verstraete, Matrix product unitaries: structure, symmetries, and topological invariants, Journal of Statistical Mechanics: Theory and Experiment 2017 , 083105 (2017).
- [133] D. M. Long, P. J. D. Crowley, A. J. Kollár, and A. Chandran, Boosting the quantum state of a cavity with floquet driving, Phys. Rev. Lett. 128 , 183602 (2022).
- [134] M. H. Kolodrubetz, F. Nathan, S. Gazit, T. Morimoto, and J. E. Moore, Topological floquet-thouless energy pump, Phys. Rev. Lett. 120 , 150601 (2018).

## Appendix A: From ALT phases to QCA

Our primary tool for studying localized topological phases is a mapping from quantum cellular automata (QCA) to many-body localized (MBL) Hamiltonians with short-ranged entanglement. Explicitly, we define a locality preserving unitary which maps physical on-site degrees of freedom to the conserved local integrals of motion (LIOMs).

We formalize our definitions of 'MBL Hamiltonian' and 'short-ranged entanglement' thereof in Appendix A1. The main point of this section and Appendix A2 is that these physical notions reproduce the formal requirements of [56, Theorem II.4]. Thus, we realize abstract restrictions on commuting projector models as physical properties of localized systems. [56, Theorem II.4] then states that the local integrals of motion (LIOMs) of the MBL Hamiltonian can be prepared by a QCA. We make an equivalent proof in Appendix A 2. The QCA in question is not unique, and we characterize the redundancy in Appendix A 3. We also discuss localized phases of driven systems, which we also characterize (non-uniquely) by parameterized families of QCA.

While our discussion is phrased for spin systems, the results we review all hold for fermions as well [56].

## 1. Setup

While all classification results in this work are for spin systems on the integer lattice Z d , the correspondence between localized phases and QCA holds for either spins or fermions and any locally finite set of sites. That is, any subset of space with finite diameter contains finitely many degrees of freedom. We will call the of points in space supporting some degree of freedom a lattice, though it does not need to have any further regularity.

Formally, we consider a metric space X and a lattice K ⊆ X (which may include the same point multiple times) with associated finite dimensional local Hilbert spaces H k . This specializes our discussion to spin systems, though a parallel discussion holds for systems with fermions [56]. The total Hilbert space is the tensor product of the local Hilbert spaces,

$$\mathcal { H } = \bigotimes _ { k \in K } \mathcal { H } _ { k }. \quad \quad ( \text{A} 1 ) \quad \underset { \text{su} } { \text{su} } \underset { \text{t} } { \text{in} } \text{$\mathcal{H}$}$$

(We will not need to be concerned with precise definitions of the infinite tensor product.)

A Hamiltonian on X is a Hermitian operator H on H . We will consider fully localized Hamiltonians with strictly local LIOMs. The generic case is for the LIOMs to be quasilocal (with exponential tails), but for the purpose of mathematical proofs it is convenient to require that LIOMs have finite support on the lattice.

Definition A.1. The Hamiltonian H is ranger fully many-body localized if there is a lattice of points J ⊆ X

(points may occur more than once) and Hermitian operators { τ z j } j ∈ J , called LIOMs, of the form

$$\tau _ { j } ^ { z } \in \mathbf 1 \otimes \bigotimes _ { k \in B _ { r } ( j ) \cap K } \text{End} ( \mathcal { H } _ { k } ), \quad ( \mathbf A 2 ) \\ \vdots B _ { r } ( j ) \text{ is the ball of radius } r \text{ around } j \Big / \text{ such that}$$

[where B r ( j ) is the ball of radius r around j ] such that [ τ z j , H ] = [ τ z j , τ z j ′ ] = 0 . Further, the set of LIOMs must be maximal, such that a local operator A commutes with all LIOMs only if A ∈ ⟨ τ z j ⟩ j ∈ J .

This definition should be compared to [56, Definition II.1].

We further say that H is short-range entangled if all its excitations can be created and destroyed locally. This excludes models with nontrivial anyons, where there are local excitations which cannot be removed except by annihilating them with other, potentially distant, excitations. It also excludes symmetry breaking, where domain walls also cannot be destroyed locally. We formalize this notion with the following convenient definition.

Definition A.2. The Hamiltonian H is ranger fully many-body localized with short-ranged entanglement if it is ranger fully many-body localized and additionally there are unitary operators

$$\tau _ { j } ^ { x } \in 1 \otimes \bigotimes _ { k \in B _ { r } ( j ) \cap K } \text{End} ( \mathcal { H } _ { k } ), \quad \text{(A3)} \\ \text{that } [ \tau _ { j } ^ { x }, \tau _ { j ^ { \prime } } ^ { z } ] \ = \ [ \tau _ { j } ^ { x }, \tau _ { j ^ { \prime } } ^ { x } ] \ = \ 0 \ f o r \ j \ \neq \ j ^ { \prime } \ a n d \ \tau _ { j } ^ { x }$$

such that [ τ x j , τ z j ′ ] = [ τ x j , τ x j ′ ] = 0 for j = j ′ and τ x j cyclically permutes the eigenspaces of τ z j . That is, denoting the ordered set of distinct eigenvalues of τ z j as ( λ , . . . λ j 0 j N j -1 ) and the associated eigenspaces eig( λ , τ j l z j ) , we have

̸

$$\tau _ { j } ^ { x } \text{eig} ( \lambda _ { l } ^ { j }, \tau _ { j } ^ { z } ) = \text{eig} ( \lambda _ { l + 1 \bmod N _ { j } } ^ { j }, \tau _ { j } ^ { z } ). \quad ( A 4 )$$

This definition should be compared to [56, Definition II.3]. One can also compare this notion of 'short-ranged entanglement' for Hamiltonians to the usual notion of short-ranged entanglement for states. A state | ψ ⟩ may be said to be short-range entangled if, for any other state | ψ ′ ⟩ which differs from | ψ ⟩ only in some ball B r ( x ) [ ⟨ ψ A ψ | | ⟩ = ⟨ ψ A ψ ′ | | ′ ⟩ for any local A not supported in B r ( x ) ] there is a local unitary τ x ψψ ′ supported in B r ( x ) such that τ x ψψ ′ | ψ ′ ⟩ = | ψ ⟩ . Any eigenstate of an H satisfying Definition A.2 is short-range entangled in this sense (the required τ x ψψ ′ can be constructed from the τ x j operators), but the converse-that τ x j operators exist for any MBL Hamiltonian H whose eigenstates are all shortrange entangled-is not obvious.

Note that a commuting model is characterized less by the Hamiltonian H , and more by the algebra of the operators { τ z j , τ x j } j ∈ J . Indeed, focusing on LIOM operator algebras will be much more useful than the Hamiltonian itself in this context. In this way, we also avoid needing to be specific regarding the resolution of degeneracies in the Hamiltonian. A classification more directly

based on the topology of the space of MBL Hamiltonians would need a definition of MBL which excludes highly degenerate Hamiltonians, such as the zero Hamiltonian, from being MBL. Otherwise all MBL Hamiltonians H 0 and H 1 would be connected by a continuous path with H t = (1 -2 ) t H 0 for t ∈ [0 , 1 2 ) and H t = (2 t -1) H 1 for t ∈ [ 1 2 , 1] . Focusing on LIOMs captures the continuous change of local properties of the model while allowing for accidental degeneracies, so we will not address how to formulate a definition of MBL which avoids this issue.

For convenience, we define some abbreviated notation for operator algebras. For S ⊆ X define

$$\mathcal { O } _ { S } = 1 \otimes \bigotimes _ { k \in S \cap K } \text{End} ( \mathcal { H } _ { k } ) \cong \bigotimes _ { k \in S \cap K } \text{Mat} _ { M _ { k } } ( \mathbb { C } ), \ \ ( A 5 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \text{where } M _ { k } = \dim ( \mathcal { H } _ { k } ) \text{ and } \text{Mat} _ { M _ { k } } ( \mathbb { C } ) \text{ is the algebra of } \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where M k = dim( H k ) and Mat M k ( C ) is the algebra of M k × M k matrices over the complex numbers. Similarly define the algebra generated by the τ operators within some region as

$$\mathcal { T } _ { S } = \langle \tau _ { j } ^ { z }, \tau _ { j } ^ { x } \rangle _ { j \in S \cap J }. \quad \quad ( \text{A6} ) \quad \text{$r^{\dots}$}$$

Definition A.3. A unitary operator V on H is a ranger quantum cellular automaton (QCA) if

$$V \mathcal { O } _ { S } V ^ { \dagger } \subseteq \mathcal { O } _ { B _ { r } ( S ) } \quad \quad ( \mathbf A 7 ) \quad P r _ { \substack { \text{se} } }$$

for all S ⊆ X .

Here, B r ( S ) is the ball of radius r around the set S . Properly, a QCA should be a projective unitary, and this distinction will be important in later appendices, but not here.

Reference [51] defines a commuting model to be 'without intrinsic topological order' when the LIOMs τ z j can be prepared from single site operators

σ z k ∈ ✶ ⊗ End( H k ) (A8) by a QCA V . Any such Hamiltonian also clearly has short-ranged entanglement by our Definition A.2, as we can construct τ x k = V σ x k V † from a matrix σ x k ∈ End( H k ) which acts as a cyclic permutation on the eigenvectors of σ z k . The result of Appendix A 2 is that the QCA V always exists given Definition A.2 as well, and as such the two notions short-ranged entanglement and lack of intrinsic topological order are equivalent.

## 2. Existence

In this appendix, we prove the following theorem.

Theorem A.1. If H is ranger many-body localized with short-ranged entanglement, then there is a range2 r QCA V such that each V † τ z j V is a sum of products of commuting single-site operators σ z k .

In fact, after a redefinition of LIOMs to another equivalent set, we will have that V † τ z j V is a single site operator.

The proof is essentially the same as [56, Theorem II.4]. The main insight here is that the abstract requirements of that theorem correspond to natural physical requirements for an MBL Hamiltonian.

The strategy will be to show that each of the algebras

$$\mathcal { T } _ { j } = \langle \tau _ { j } ^ { z }, \tau _ { j } ^ { x } \rangle \quad \quad ( A 9 ) \\ \cdot \cdot \cdot \cdot \cdot \quad \text{--} \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot$$

is isomorphic to a matrix algebra. This implies that there is an isomorphism between T j and O k whenever they have the same dimension. [All isomorphisms will be ∗ -isomorphisms-that is an isomorphism α such that α A ( † ) = α A ( ) † for all operators A . All such automorphisms of this kind between matrix algebras can be expressed as a conjugation by a unitary.] To find a QCA V which implements this isomorphism simultaneously for all sites, we need to match the algebras O k to T j ( k ) locally and bijectively. This can be phrased as a graph problem, and Hall's marriage theorem [122, 123] implies that this problem always has a solution.

Lemma A.2. There is an isomorphism

$$\mathcal { T } _ { j } \cong \text{Mat} _ { N _ { j } } ( \mathbb { C } ) \text{ \quad \ \ } ( A 1 0 )$$

for each j ∈ J .

Proof. Wedderburn's classification of finite dimensional semisimple algebras implies that T j is isomorphic to a full matrix algebra if and only if it is central. That is, A ∈ T j commutes with all of T j only if A ∈ C .

✶ Suppose A ∈ Z ( T j ) is in the center of T j . Then it must commute with all the LIOMs-it commutes with τ z j by assumption and it commutes with the other τ z l because it is contained in T j . Then, by Definition A.1, A is a sum of products of LIOMs. However, as A ∈ T j , it must be a polynomial in τ z j . In particular,

$$\hat { A } = \bigoplus _ { l = 1 } ^ { N _ { j } - 1 } a _ { l } \mathbb { 1 } _ { \text{eig} ( \lambda _ { l } ^ { j }, \tau _ { j } ^ { z } ) } \quad ( A 1 1 ) \\ \hat { \iota } _ { l } \in \mathbb { C }. \text{ Evaluating the commutator with $\tau_{j}^{x}$,}$$

for some a l ∈ C . Evaluating the commutator with τ x j , we have

τ x † j Aτ x j = N j -1 ⊕ l =1 a l +1 mod N j ✶ eig( λ ,τ j l z j ) = A, (A12) so all the a l = a are equal and A ∈ C ✶ . It is convenient to require that all the N j and M k = dim( H k ) are prime. This can always be achieved by redefining the lattices J and K . Indeed, focusing on N j , we can make a prime factorization

$$N _ { j } = p _ { 0 } p _ { 1 } \cdots p _ { m _ { j } - 1 } \quad \quad ( A 1 3 ) \\.. \. \. \. \. \.$$

(not all necessarily distinct) and express

$$\mathcal { T } _ { j } \cong \text{Mat} _ { N _ { j } } ( \mathbb { C } ) \cong \bigotimes _ { n = 0 } ^ { m _ { j } - 1 } \text{Mat} _ { p _ { n } } ( \mathbb { C } ). \quad \quad ( A 1 4 )$$

Then we have that T j = ∏ m -1 n =0 T j n with each tensor factor being of prime dimension. Further, we can make this factorization such that the eigenbasis of τ z j becomes a product basis. Indeed, we have that any l ∈ { 0 , . . . , N j - } 1 can be uniquely expressed as a sum of integers

$$l = \sum _ { n = 0 } ^ { m _ { j } - 1 } \mu _ { n } ^ { l } \prod _ { q = 0 } ^ { n - 1 } p _ { q } \, \quad \, \quad \, ( A 1 5 )$$

with µ l n ∈ { 0 , . . . , p n -1 } . (Empty sums are taken to be 0 and empty products are 1.) Then we can define Hermitian projectors π j µn (with µ ∈ { 0 , . . . , p n - } 1 and n ∈ { 0 , . . . , m j - } 1 ) onto the direct sum of eigenspaces eig( λ , τ j l z j ) such that µ l n = µ . Note that all π j µn with the same n have the same rank, specifically N /p j n . As T j is a full matrix algebra, we can find τ x j n such that

$$\tau _ { j _ { n } } ^ { x } \pi _ { \mu n } ^ { j } \tau _ { j _ { n } } ^ { x \dagger } = \pi _ { ( \mu + 1 \bmod p _ { n } ) n } ^ { j } \quad \quad ( A 1 6 ) \quad \text{tha} \\ \text{nen}$$

and such that τ x j n commutes with all the other projectors. Indeed, τ x j n is defined by cyclically permuting the value of µ l n in the decomposition of l for eig( λ , τ j l z j ) . By essentially the same calculation as Lemma A.2, ⟨ τ x j n , π j µn µ ⟩ p n -1 =0 is isomorphic to a full matrix algebra, and so occurs as a tensor factor in T j . Then we make the decomposition of T j as

$$\mathcal { T } _ { j } \cong \bigotimes _ { n = 0 } ^ { m _ { j } - 1 } \langle \tau _ { j _ { n } } ^ { x }, \pi _ { \mu n } ^ { j } \rangle _ { \mu = 0 } ^ { p _ { n } - 1 }. \quad \quad ( A 1 7 ) \quad \underset { \text{dim} } { \text{du} }$$

We pick new LIOMs τ z j n = ∑ p n -1 µ =0 λ j n µ π j µn with λ j n µ ∈ R all distinct. These commute with τ z j , and so are each actually functions of τ z j by Definition A.1. On the other hand, they generate all diagonal operators in T j in the τ z j eigenbasis, so τ z j can also be recovered from these new LIOMs. Thus, the new LIOMs are equivalent to the old LIOM.

Replace the lattice J with a new lattice J ′ which consists of the old lattice sites j each repeated m j times, and associated LIOMs τ z j n and algebras T j n . This lattice has the properties we need. We henceforth drop primes, and just use J . Similarly, we can ensure that each of the local Hilbert spaces H k is of prime dimension.

As T j and O k are both full matrix algebras, there is an isomorphism between them whenever they have the same dimension. To build these individual isomorphisms into a QCA, we need to make a matching between T j and O k algebras which are supported near one another. Define the bipartite graph B = (( J, K ) , E ) with edges

$$( j, k ) \in E \iff N _ { j } = M _ { k } \text{ and dist} ( j, k ) \leq r. \quad ( A 1 8 )$$

↦

Then we are looking for a perfect matching in this graph. That is, a bijection k → j ( k ) such that ( j ( k , k ) ) ∈ E for all k ∈ K . Hall's marriage theorem [122] (specifically, the infinite version due to Hall [123]-no relation) gives necessary and sufficient conditions for a perfect matching to exist. Define the neighborhood of S ⊆ B as

$$\overset { \prime } { y } \quad \partial S = \{ b \in B \ \colon ( s, b ) \in E \ \text{for some} \ s \in S \}. \quad ( A 1 9 )$$

Then the marriage theorem states that if: (1) every vertex has finite degree, (2) for every finite S ⊆ J we have

$$| S | \leq | \partial S |, \text{ \quad \ \ } ( A 2 0 )$$

and similarly (3) for each finite S ⊆ K we have | S | ≤ | ∂S | , then the required perfect matching bijection exists. The criterion (A20) is called the marriage condition .

Lemma A.3. The graph B has a perfect matching.

Proof. We verify the marriage condition, Eq. (A20).

Note that every vertex in B has finite degree by the assumption that K and J are locally finite. Further observe that the graph B splits into distinct connected components labelled by primes p such that N j = M k = p . Thus, we can consider each component separately. Henceforth, we assume that we are working in a given component.

Let S ⊆ J be a finite subset of LIOM sites. Then we have that T S is a subalgebra in O B r ( S ) . As both T S and O B r ( S ) are full matrix algebras, O B r ( S ) must be a tensor product of T S and its commutant in O B r ( S ) . Comparing dimensions then shows that ∏ j ∈ S N j divides ∏ k ∈ B r ( S ) M k . Then, all of the primes N j = p must occur among the primes M k , repeated according to their multiplicity. The sites in B r ( S ) with local Hilbert space dimension M k = p are precisely the elements of ∂S in this component. Thus, we conclude | S | ≤ | ∂S | .

✶ We can now prove Theorem A.1.

In the other direction, let S ⊆ K be a finite subset of lattice sites. Then we also have that O S is a subalgebra of T B r ( S ) , whereupon the proof proceeds as in the previous paragraph. Indeed, suppose that A ∈ O S ⊆ O B 2 r ( S ) . As T B r ( S ) ⊆ O B 2 r ( S ) occurs as a tensor factor in O B 2 r ( S ) , A is a sum of products of operators in T B r ( S ) and its commutant. We need to show that if A is in the commutant, then it is a scalar. Indeed, as their supports do not overlap, A must also be in the commutant of T K B \ r ( S ) . Thus, A commutes with all of the LIOMs, and by Definition A.1 we must have A ∈ ⟨ τ z j ⟩ j ∈ J . The only such operator which also commutes with all the τ x j operators is a scalar, A = a .

Proof of Theorem A.1. Let j : K → J be a perfect matching as in Lemma A.3. By Lemma A.2 there is an isomorphism of algebras

$$\alpha _ { k } \colon \mathcal { O } _ { k } \to \mathcal { T } _ { j ( k ) }. \quad \quad ( A 2 1 )$$

Then

$$\alpha = \bigotimes _ { k \in K } \alpha _ { k } \quad \quad ( A 2 2 )$$

is an algebra isomorphism between ⊗ k ∈ K O k and ⊗ j ∈ J T j . It is bijective because the matching k j ( ) is perfect, and it is a homomorphism because all the algebras T j commute with one another.

All such complex matrix algebra homomorphisms are of the form α A ( ) = V AV † for some unitary V , which is the QCA we need. It maps all LIOMs to single-site operators, thus meeting the requirements of the theorem.

Further, α is 2 r -locality preserving: we have dist( j ( k , k ) ) ≤ r by the definition of the graph B , and then α ( O k ) = T j ( k ) ⊆ O B 2 r ( k ) by the fact that T j ( k ) is supported in O B r ( j ( k )) and the triangle inequality.

Define the single-site operators σ z k = V † τ z j ( k ) V . Then we have that

$$V ^ { \dagger } H V = \sum _ { k } h _ { k } \sigma _ { k } ^ { z } + \sum _ { k l } h _ { k l } \sigma _ { k } ^ { z } \sigma _ { l } ^ { z } + \cdots \quad ( A 2 3 ) \quad \text{erat} \quad \text{In} \quad \text{(A)}$$

is a sum of products of the commuting single-site operators σ z k (by Definition A.1). We note that this is just a diagonalization of H in the simultaneous eigenbasis of all the σ z j . The nontrivial conclusion of this section is that this diagonalization can be made locally.

This characterization of localized system with shortranged entanglement by QCA motivates a similar definition for driven systems. We consider an n -tone-driven system to be localized with short-ranged entanglement when it has a generalized Floquet decomposition [28, 99]

$$U ( t ) = V ( \vec { \theta } _ { t } ) e ^ { - i t H _ { F } } V ( \vec { \theta } _ { 0 } ) ^ { \dagger }, \quad \ \ ( A 2 4 ) \quad \text{tati}$$

where the static Floquet Hamiltonian

$$H _ { F } = \sum _ { k } h _ { k } \sigma _ { k } ^ { z } + \sum _ { k l } h _ { k l } \sigma _ { k } ^ { z } \sigma _ { l } ^ { z } + \cdots \quad \text{ (A2.5)} \quad \text{ the } 1 \\ \text{ As s} \text{ s}$$

is a sum of products of the fixed σ z k operators and the micromotion operator V ( ⃗ θ ) is a continuous parameterized family of QCA. By a similar proof to the above, we see that this definition is equivalent to the existence of a maximal set of LIOMs τ z j ( ⃗ θ ) with explicit continuous ⃗ θ dependence and a family of conjugate τ x j ( ⃗ θ ) operators also with continuous ⃗ θ dependence. Alternatively, the definition is also equivalent to the quasienergy operator [98, 124-127] having a complete set of LIOMs [99] and conjugate operators which create or destroy excitations locally.

The requirement of the family τ x j ( ⃗ θ ) being continuous is nontrivial. For instance, even for 2-tone-driven qubits we may construct continuous families of LIOMs such that each family of LIOM eigenstates has a nonzero Chern number [127, 128]. However, there is an obstruction to continuously defining a family of conjugate operators, as the eigenstates cannot be given a smooth gauge. We do not consider such systems to be localized as they have continuous spectrum [128], and so do not enjoy the same stability to generic perturbations we expect of localized systems [99]. Indeed, in a frequency lattice picture, such systems are delocalized [128].

We emphasize that Floquet localized phases with short-ranged entanglement correspond to loops of QCAnot loops in the space of QCA modulo generalized permutations, even if QCA which differ by generalized permutations correspond to the same LIOMs. Indeed, a loop of QCA modulo permutations could correspond to a path from the identity QCA to a permutation. This path of QCA does not define a legitimate micromotion operator V ( ⃗ θ ) .

## 3. Non-uniqueness

A QCA which maps H to a sum of product of σ z k operators always exists, but this QCA is not unique.

Indeed, the QCA V is a unitary operator which diagonalizes H in the simultaneous eigenbasis of { σ z k } k ∈ K . Even if H is nondegenerate, such unitaries are only unique up to reordering of eigenvalues and multiplying eigenstates by arbitrary phases. Specifically, if P is a unitary generalized permutation matrix (a unitary matrix with exactly one unit modulus complex number in each row and column) in the eigenbasis of { σ z k } k ∈ K , then ( V P ) † H V P ( ) is still diagonal. If H is nondegenerate, the full redundancy of the correspondence between localized models with short-ranged entanglement and QCA is found by restricting P to be a QCA.

Translations form an important class of QCA permutation operators [58, 61]. In one dimension, all QCA are given by translation operators, up to a finite depth circuit. Thus, while one-dimensional QCA are classified by log Q × [58], all these nontrivial QCA land within the redundancy of the description of localized phases. As such, all one-dimensional localized phases of spins without symmetry are trivial. The same is not true of fermions, where the Kitaev chain [129] represents a nontrivial localized phase which can be prepared by a QCA [117]. It has invertible eigenstate order, and so is not an ALT phase.

Mathematically, the classification of localized phases with short-ranged entanglement is given by a quotient of π 0 (QCA ) d . Denoting the set of generalized permutation QCA by GP d , there is a canonical inclusion map ι : GP d → QCA d , and an induced homomorphism ι ∗ : π 0 (GP ) d → π 0 (QCA ) d . We write the image of this homomorphism as I 0 = ι ∗ ( π 0 (GP )) d . Then the quotient which classifies localized topological phases is

$$\pi _ { 0 } ( \text{QCA} _ { d } ) / I _ { 0 }. \text{ \quad \ \ } ( A 2 6 )$$

This group is isomorphic to π 0 (QCA d / GP ) d . Indeed, there is a fibration sequence

$$\text{GP} _ { d } \stackrel { \iota } { \rightarrow } \text{QCAT} _ { d } \stackrel { q } { \rightarrow } \text{QCAT} _ { d } / \text{GP} _ { d } \quad \text{ (A27)}$$

(where q is the quotient map which sends a QCA V to its coset V GP d ) with associated long exact sequence of

$$\dots \to \pi _ { 0 } ( \text{GP} _ { d } ) \stackrel { \iota ^ { * } } { \longrightarrow } \pi _ { 0 } ( \text{QCA} _ { d } ) \stackrel { q ^ { * } } { \longrightarrow } \pi _ { 0 } ( \text{QCA} _ { d } / \text{GP} _ { d } ) \to 0, \quad \begin{array} { c } \text{gener} \\ ( A 2 8 ) & \text{spond} \end{array} \end{array}$$

and π 0 (QCA d / GP ) = d ∼ π 0 (QCA ) d /I 0 follows from the first isomorphism theorem applied to q ∗ .

In the driven case, the redundancy is captured by a family of generalized permutation matrices P ⃗ ( θ ) with elements

$$P _ { \alpha \beta } ( \vec { \theta } ) = \delta _ { \pi ( \alpha ) \beta } e ^ { - i \vec { m } _ { \beta } \cdot \vec { \theta } - i \phi _ { \beta } }, \quad \ \ ( A 2 9 ) \quad \$$

where π is a permutation of the simultaneous eigenstates of the { σ z k } k ∈ K operators and ⃗ m α ∈ Z n are vectors of integers. [Such P ⃗ ( θ ) are generalized permutation matrices in a frequency lattice picture [31, 99].] Any P of this form defines an identical unitary evolution (A24) with

$$V ^ { \prime } ( \vec { \theta } ) = V ( \vec { \theta } ) P ( \vec { \theta } ), & & ( A 3 0 ) \quad \overset { \text{a} \, \varepsilon } { \text{in} \, t }$$

$$H _ { F } ^ { \prime } = H _ { F } - i [ \partial _ { t } P ( \vec { \theta } _ { t } ) ] P ( \vec { \theta } _ { 0 } ) ^ { \dagger }. \quad ( A 3 1 ) \quad \text{$o gy$} \\. \quad \text{$\text{$spa$} }$$

Note that i ∂ P [ t ( ⃗ θ t )] P ⃗ ( θ t ) † is a static diagonal operator with elements

$$\begin{array} {$$

where ⃗ ω = ( ω , . . . , ω 1 n ) is the vector of drive frequencies. If P ⃗ ( θ ) is also a QCA for all ⃗ θ , then V ′ ( ⃗ θ ) also defines an acceptable generalized Floquet decomposition.

The only example for which this redundancy affects the classification of localized phases for d ≤ 3 is the classification of Floquet ( n = 1 ) phases in zero dimensions [31]. Looking ahead to Appendix C 2, we have that π 1 (QCA ) = 0 Q Z / , and that a complete index is given by the winding number ˜ [ w 1 V ] . In zero dimensions, all unitaries are QCA, so the diagonal unitary with elements

$$P _ { m } ( \theta ) _ { \alpha \beta } = e ^ { i$$

defines an acceptable V ′ ( θ ) = V ( θ P ) m ( θ ) . However, we can compute that

$$\tilde { w } _ { 1 } [ V P _ { m } ] = \tilde$$

where N is the Hilbert space dimension. Thus, multiplication by P m ( θ ) can take ˜ [ w 1 V ′ ] to any allowed value, and we conclude that all loops of zero dimensional QCA (loops of projective unitaries) define trivial localized phases.

More systematically, we can again classify driven localized phases by a quotient of the group π 0 [ T n , QCA ] d . The generalized permutations P ⃗ ( θ ) are particular maps from T n to QCA d , and we denote the subgroup of π 0 [ T n , QCA ] d induced by the inclusion of the generalized permutations in the full space of maps from T n to QCA d by I n . Then the classification of n -tone driven localized phases in d -dimensions with short-ranged entanglement is given by

$$\pi _ { 0 } \left [ \mathbb { T } ^ { n }, \text$$

This group is not isomorphic to π 0 [ T n , QCA d / GP ] d in general. Taking, for instance, n = 1 , there are loops in π 0 [ T 1 , QCA d / GP d ] ∼ = π 1 (QCA d / GP ) d which correspond to paths from the identity to a locally generated permutation. This is not a legitimate micromotion operator V ( θ ) , as it is not periodic as a (projective) unitary.

## Appendix B: The Ω -spectrum property

In topology, an Ω -spectrum refers to a sequence of topological spaces ( Y d ) d ∈ Z together with (weak) homotopy equivalences

$$Y _ { d - 1 } \simeq \Omega Y _ { d } \text{$$

where Ω Y d is the space of loops in Y d which begin and end at a specified base point. Ω -spectra play a special role in topology, as they serve to define generalized cohomology theories -functors from the category of topological spaces to abelian groups obeying certain axioms [95, 102104]. The functor in question is

$$h _ { Y } ^ { d } ( X ) = \pi _ { 0 } \left [ X, Y _ { d } \right ]. \text{ \quad \ \ } \text{(B2)}$$

The Ω -spectrum property comes with substantial structure, and highly constrains the properties of the spaces Y d .

We conjecture that the spaces of quantum cellular automata on the integer lattice Z d stabilized by adding ancillae, QCA d , form an Ω -spectrum,

$$Q C A _ { d - 1 } \simeq \Omega Q C A _ { d }. \text{ \quad \ \ } ( B 3 )$$

For d &lt; 0 , we define QCA d = Ω -d QCA 0 , making the negative degrees of this statement trivial. In positive degree, we propose that the swindle map S [Eq. (8)] provides the homotopy equivalence. Thus, we propose there is a generalized cohomology theory

$$h ^ { d } ( X ) = \pi _ { 0 } \left [ X, Q C A _ { d } \right ]. \text{ \quad \ \ } \text{(B4)}$$

In this appendix, we provide arguments in support of this conjecture. In Appendix B1, we recall that if loops of QCA can be put on open boundary conditions such that they act by QCA on their edges-known as blending with the identity -then loops of QCA can all be parameterized as finite depth circuits. We also argue that this condition holds for loops of one-dimensional QCA. Then, in Appendix B2, we construct an explicit deformation retract from the space of circuits to the image of the swindle map. As the swindle map is injective, this shows

$$Q C A _ { d - 1 } \cong \lim ( S ) \simeq \Omega Q C A _ { d } \quad \quad ( B 5 )$$

under the assumptions of Appendix B 1.

Making these arguments rigorous requires being more specific about the topology on QCA d , even in d = 1 . This can be done, but we defer such technicalities to Ref. [107].

While it may be useful to keep spin systems in mind, the arguments we present in this appendix are agnostic to local operators being fermionic or spin operators. Provided that the additional ancillae used to stabilize a system may support fermionic operators, everything in the following sections also applies to systems with fermions.

## 1. Blending and circuits

A blending between α and β ∈ QCA d is a third QCA γ such that

$$\gamma ( A ) & = \alpha ( A ) \quad \text{for supp} ( A ) \subseteq ( - \infty, x _ { l } ] \times \mathbb { Z } ^ { d - 1 }, \quad ( B 6 ) \quad \text{QCA}, \\ \psi ( A \vee A \vee A \vee \frac { \epsilon } { \epsilon } \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \ dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots$$

$$\gamma ( A ) = \beta ( A ) \text{ \ for supp} ( A ) \subseteq [ x _ { r }, \infty ) \times \mathbb { Z } ^ { d - 1 }, \text{ \ \ } ( B 7 ) \text{ \ \ is not\ \ } \text{the\ } \text{$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$ \quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\ \text{$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$$$

where supp( A ) is the smallest subset of Z d on which the operator A acts non-trivially, and x l,r are positions along the first axis in Z d . That is, γ agrees with α to the left of some cut, and agrees with β to the right of the cut. We can similarly define a blending γ θ ( ⃗ ) between parameterized families of QCA α ⃗ ( θ ) and β ⃗ ( θ ) . Two (families of) QCA are said to blend into each other if a blending exists between them.

Two QCA blend into one another if and only if they are related by a finite depth circuit [94]. That is, α blends with β at any position of the cut x if and only if αβ -1 is equal to a finite depth unitary circuit. (Actually, it is sufficient for a blending to exist at any position [109].)

Indeed, it is sufficient to consider a blending between γ = αβ -1 and the identity. Then γ should be equal to a finite depth circuit. Let the blending at the cut across { x }× Z d -1 be ˜ γ x , and coarse grain such that ˜ γ agrees with γ for x ′ &lt; x and acts as the identity for x ′ &gt; x . Then we have that ˜ γ γ x ˜ -1 x +1 and ˜ γ -1 x ˜ γ x +1 are QCA supported on { x, x +1 } × Z d -1 .

We can split γ into a depth two circuit by iteratively splitting off gates of the form ˜ γ γ x ˜ -1 x +1 or ˜ γ -1 x ˜ γ x +1 . (The gates are themselves QCA in d -1 dimensions.) First, observe that we can split γ into a product of two semiinfinite QCA

$$\gamma = \tilde { \gamma } _ { 0 } ( \tilde { \gamma } _ { 0 } ^ { - 1 } \gamma ). \quad \quad \quad ( B 8 )$$

Then we can express the ˜ γ 0 factor as

$$\tilde { \gamma } _ { 0 } = ( \tilde { \gamma } _ { 0 } \tilde { \gamma } _ { - 1 } ^ { - 1 } ) \tilde { \gamma } _ { - 1 }. \quad \quad ( B 9 ) \quad T t$$

The first factor is the gate we need. We can proceed similarly with the second factor, decomposing ˜ γ -1 = ˜ γ -2 (˜ γ -1 -2 ˜ γ -1 ) and so on, alternating whether we perform the multiplication by ˜ γ -1 x ˜ γ x on the right or ˜ γ γ x ˜ -1 x on the left. Mirroring the same procedure for positive x , we obtain a depth-two circuit decomposition

$$\gamma = \left ( \prod _ { x \in 2 \mathbb { Z } + 1 } w _ { x, x + 1 } \right ) \left ( \prod _ { x \in 2 \mathbb { Z } } w _ { x, x + 1 } \right ). \quad ( B 1 0 ) \quad u _ { 0, 1 }$$

The converse direction for the equivalence of blending and circuit equivalence is straightforward. If γ is a finite depth circuit, we can simply delete all gates to one side of any cut to obtain a blending with the identity.

It is intuitive that a blending should exist for any path of QCA beginning at the identity. Such a path of unitary operators γ t ( ) can be generated by a local Hamiltonian H t ( ) if it is differentiable. Then tracing over the sites to the right of any cut gives a new Hamiltonian Tr ( x, ∞ × ) Z d -1 H t ( ) which, assuming it continues to generate a QCA, generates a blending between γ t ( ) and the identity. This is certainly true in one dimension, where any unitary evolution at the zero-dimensional edge is a QCA, but in higher dimensions the new Hamiltonian evolution may develop exponential tails along the edge. This is not much of a problem physically, but such tails are an obstacle to performing proofs.

Nonetheless, we expect that the path-space P QCA d can be continuously parameterized in terms of a strictly local path of circuits in dimensions higher than one too. The loop space ΩQCA d can then be parameterized in terms of circuits which begin and end at the identity. Note that this does not require that the individual gates must return to the identity, only that they cancel at the end of the path.

In Ref. [107] we study this structure in greater detail, and characterize the circuit liftings implicit here as a fibration on the space of QCA.

## 2. Deformation to the swindle

We now argue that Eq. (B5) holds. First, we show that S : QCA d -1 → ΩQCA d is injective, such that im( S ) = QCA ∼ d -1 is homeomorphic to QCA d -1 . That is, S v, t ( ) = S w,t ( ) for all t only if v = w . Indeed, multiplying S v, t ( ) on the left by the inverse of the top layer in S w,t ( ) and on the right by the inverse of the bottom layer (Fig. 2), we have

$$\text{st}, \\ \text{ni} - \left ( \prod _ { x \in \mathbb { Z } + 1 } u _ { x, x + 1 } ( w, t ) ^ { \dagger } u _ { x, x + 1 } ( v, t ) \right ) \\ \times \left ( \prod _ { x \in \mathbb { Z } } u _ { x, x + 1 } ( v, t ) u _ { x, x + 1 } ( w, t ) ^ { \dagger } \right ) = 1. \ \ ( B 1 1 ) \\ \text{ This equation can only be satisfied if, for all } t \in [ 0, 1 ], \\.$$

This equation can only be satisfied if, for all t ∈ [0 , 1] ,

$$\underset { - } { \ e d } \quad u _ { 2 x, 2 x + 1 } ( v, t ) u _ { 2 x, 2 x + 1 } ( w, t ) ^ { \dagger } = q _ { 2 x } ( t ) q _ { 2 x + 1 } ( t ) \quad ( B 1 2 )$$

$$\rho _ { \text{rm} } ^ { = } \sum _ { u _ { 2 x - 1, 2 x } ( w, t ) ^ { \dagger } u _ { 2 x - 1, 2 x } ( v, t ) = q _ { 2 x - 1 } ( t ) ^ { \dagger } q _ { 2 x } ( t ) ^ { \dagger } \ ( B 1 3 ) }$$

for some sequence of q x ( ) t supported on { x }× Z d -1 . However, the gates u x,x +1 are translationally invariant, so considering x = 0 gives

$$u _ { 0, 1 } ( v, t ) u _ { 0, 1 } ( w, t ) ^ { \dagger } = q _ { 0 } ( t ) q _ { 1 } ( t ) = u _ { 0, 1 } ( v, t ) ^ { \dagger } u _ { 0, 1 } ( w, t ), \\ ( B 1 4 )$$

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00825v2-figure-ac6a01a3-5c8b-4a16-be37-acf34f6bb0d7.jpg)

-→

-→

FIG. 4. The loop of finite depth circuits γ t ( ) ⊗ S v ( γ 0 , s ) † is deformed to the constant loop at the identity as a function of s ∈ [0 , 1] . (a) At s = 0 , the circuit consists of two layers. Between s = 0 and s = 1 2 / the gates u † x,x +1 ( v γ 0 , t ) in the top layer (red) are swapped to the bottom layer while having their order of multiplication reversed. (b) The gate u † 0 1 , ( v γ 0 , t ) w 0 1 , ( ) t [ u † -1 0 , ( v γ 0 , t ) w -1 0 , ( ) t ] acts as the identity on site 0 for t ∈ { 0 1 , } . It can thus be swapped (over s ∈ [ 1 2 , 3 4 ] ) so that it acts on site 1 [ -1 ] and its ancilla while ensuring that the whole circuit is still a loop. (c) Between s = 3 4 and s = 1 the circuit (red) is 'rolled up', such that it acts as the identity (gray dots) in any finite region.

where the second equality comes from translating Eq. (B13) by one unit towards positive x and inverting the equation. We see that u 0 1 , ( v, t ) 2 = u 0 1 , ( w,t ) 2 . Making a spectral decomposition of both sides, we conclude that u 0 1 , ( v, t ) = u 0 1 , ( w,t ) up to multiplying eigenvalues by ± 1 . However, these paths both start at t = 0 with the same eigenvalues λ = +1 , which remain nonzero throughout all t . Thus, all the discrete choices of ± 1 must be +1 , and u 0 1 , ( v, t ) = u 0 1 , ( w,t ) . Then comparing these at t = 1 gives v v 0 † 1 = w w 0 † 1 , and hence v = w up to a phase. We regard QCA as projective unitaries, so v = w .

We now assume that all loops of QCA γ t ( ) can be represented by paths of circuits (Appendix B 1)

$$\gamma ( t ) = \left ( \prod _ { x \in 2 \mathbb { Z } + 1 } w _ { x, x + 1 } ( t ) \right ) \left ( \prod _ { x \in 2 \mathbb { Z } } w _ { x, x + 1 } ( t ) \right ) \quad ( B 1 5 ) \quad \stackrel { \text{rows} } { s \in [ 1 } } \\ \tilde { \Phi } _ { x } \tilde { \Phi } _ { x } \tilde { \Phi } _ { x } \tilde { \Phi } _ { x } \tilde { \Phi } _ { x } \tilde { \Phi } _ { x } \tilde { \Phi } _ { x } \tilde { \Phi } _ { x } \tilde { \Phi } _ { x } \tilde { \Phi } _ { x } \tilde { \Phi } _ { x } \tilde { \Phi } _ { x } \tilde { \Phi } _ { x } \tilde { \Phi } _ { x } \tilde { \Phi } _ { x } \tilde { \Phi } _ { x } \tilde { \Phi } _ { x } \tilde { \Phi                                                                                                                                                                                                         } }$$

such that w x,x +1 (0) = ✶ and the circuit is the identity at t = 1 . This requires that w x,x +1 (1) = v γ x ⊗ v γ † x +1 for some sequence v γ x ∈ QCA d -1 which are not necessarily all equal.

Then, up to stable equivalence, we have

γ t ( ) = γ t ( ) ⊗ ✶ = γ t ( ) ⊗ S v ( γ 0 , t ) † S v ( γ 0 , t ) . (B16) We now make a sequence of deformations which cancel the γ t ( ) and S v , t ( 0 ) † parts of Eq. (B16), leaving just ✶ ⊗ S v , t ( 0 ) . That is, we define a continuous function E : ΩQCA d × [0 , 1] → ΩQCA d (B17)

$$( \gamma ( t ), s ) \mapsto E ( \gamma, t, s ) \quad ( B 1 8 ) \quad h \,.$$

↦

such that E γ,t, ( 0) = γ t ( ) ⊗ ✶ and E γ,t, ( 1) = ✶ ⊗ S v ( γ 0 , t ) . The deformations all take the form of swapping the positions of gates in the circuit representation of Eq. (B16) (Fig. 4). The S v ( γ 0 , t ) factor is a spectator, and so we will focus on γ t ( ) ⊗ S v , t ( 0 ) † .

First, in s ∈ [0 , 1 2 ] , we swap the gates u † x ,x ′ ′ +1 ( v γ 0 , t ) (we denote ancillae sites with a prime, x ′ ) from the second tensor factor-the second layer -to the first, such that the new gates in the first layer become

$$w _ { x, x + 1 } ^ { \prime } ( t ) = u _ { x, x + 1 } ^ { \dagger } ( v _ { 0 } ^ { \gamma }, t ) w _ { x, x + 1 } ( t ). \quad \text{(B19)}$$

This also requires reversing the order of multiplication in

$$S ^ { \dagger } = \left ( \prod _ { x ^ { \prime } \in 2 \mathbb { Z } } u _ { x ^ { \prime }, x ^ { \prime } + 1 } ^ { \dagger } \right ) \left ( \prod _ { x ^ { \prime } \in 2 \mathbb { Z } + 1 } u _ { x ^ { \prime }, x ^ { \prime } + 1 } ^ { \dagger } \right ), \quad ( B 2 0 ) \\ \sim \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \delta \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \dotdet \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det\det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det$$

so that it matches the circuit structure of γ t ( ) . This can be done by adding the gates SWAP x,x ′ (2 s ) to the positions marked in Fig. 4(a).

Observe that w ′ 0 1 , (1) = ✶ ⊗ v γ 0 v γ † 1 and w ′ -1 0 , (1) = v γ † 0 v γ -1 ⊗ ✶ (B21) both act as the identity on site x = 0 at t = 1 . Thus, in s ∈ [ 1 2 , 3 4 ] , we can introduce SWAP(4( s -1 2 )) gates to the positions shown in Fig. 4(b) which 'roll up' the circuitthe gates previously acting on site 0 now act on-site at x = ± 1 between ancillae. Explicitly, for s ∈ [ 1 2 , 3 4 ) , the gate in E γ,t, s ( ) acting across the (00 ′ , 11 ) ′ edge is

$$\underset { \text{tbing} } { \text{tbing} } \quad \text{SWAP} _ { 0, 1 ^ { \prime } } ( 4 ( s - \frac { 1 } { 2 } ) ) w _ { 0, 1 } ^ { \prime } ( t ) \text{SWAP} _ { 0, 1 ^ { \prime } } ( 4 ( s - \frac { 1 } { 2 } ) ) ^ { \dagger }. \ ( B 2 2 ) \\ \text{taken} \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \text{, } \quad \end{$$

This rolling up procedure can be iterated. At s = 3 4 , we see that the gate w ′ 1 2 , ( ) t w ′ 1 ′ , 1 ( ) t is the identity on sites 1 and 1 ′ when t = 1 . Thus, by adding ancillae to site 2 , this gate can also be rolled up within s ∈ [ 3 4 , 7 8 ] [Fig. 4(c)]. We continue rolling up sites x and -x within s ∈ [1 -2 - -x 1 , 1 -2 - -x 2 ] , such that, as s → 1 , we have E γ,t, s ( ) → ⊗ ✶ S v ( γ 0 , t ) on any finite patch of the lattice. Technically, this makes E γ,t, s ( ) continuous in the strong operator topology.

As w ′ 0 1 , (1) acts as the identity on site 0, the introduced swap gates cancel at t = 1 , and the entire circuit remains a loop for all s . By s = 3 4 , the gate acts as the identity on site 0 for all t ∈ [0 , 1] . We repeat an analogous deformation on the other gate with support on site 0 . The part of the circuit which is supported on x = 0 is thus rolled up to act on the ancillae at sites x = ± 1 .

The deformation is completed by adding a parallel layer of swap gates which takes ✶ ⊗ S v ( γ 0 , t ) to S v ( γ 0 , t ) ⊗ ✶ ∈ im( S ) . If the circuit decomposition Eq. (B15) gives w x,x +1 = u x,x +1 ( v ) [as it does when using the Tr ( x, ∞ × ) Z d H t ( ) construction of Appendix B 1] then E is a deformation retract, and Eq. (B5) holds.

## Appendix C: Homotopy groups of QCA

The Ω -spectrum property reduces the calculation of homotopy groups π n (QCA ) d of QCA to either finding

the connected components π 0 (QCA d -n ) or the homotopy groups of a zero-dimensional QCA, π n -d (QCA ) 0 . In this appendix, we summarize progress in the literature in computing π 0 (QCA d -n ) for spin systems (Appendix C1), and compute all the groups π n -d (QCA ) 0 (Appendix C2). We briefly comment on the classifcation of QCA on lattices other than Z d in Appendix C3, but leave any technical arguments to Ref. [107].

Unlike the previous two appendices, the results of this appendix are specific for spin systems. The homotopy groups of the space of fermionic QCA will differ from the groups presented here.

## 1. Connected components

In general, the connected components π 0 (QCA ) d are not known. However, these groups have been calculated for d ∈ { 1 2 , } (and d = 0 , which is trivial), and the groups are known to be nontrivial in some higher dimensions.

For one-dimensional QCA, the complete classification [that is, π 0 (QCA ) 1 ] was obtained by Ref. [58], which also explicitly defined an invariant which diagnoses which connected component a given QCA belongs to. All onedimensional QCA are path connected to a shift -a QCA which acts as a translation on some subset of the degrees of freedom. For example, in a chain of qubits, each QCA is connected to a QCA which translates all qubits by r sites either right or left. For qubit chains, QCA are classified by the group Z .

However, we stabilize QCA by allowing ancillae of any Hilbert space dimension to be added to any site by tensor product. Thus, we must assume that the chain also has qutrit (three-state) degrees of freedom, and other higher level qudits. While a 4-level qudit can be decomposed into two qubits, and so does not give new invariants compared to the qubit system, qutrits and other prime qudits all produce invariants. A shift of qutrits cannot be connected to a shift of qubits.

The resulting classifying group is succinctly presented as

$$\pi _ { 0 } ( Q C A _ { 1 } ) \cong \log \mathbb { Q } ^ { \times }, \text{ \quad \ \ } ( C 1 )$$

where Q × is the set of positive rationals with multiplication as the group operation, and the group operation in log Q × is addition. Connecting to our discussion, this can be written as

$$\log \mathbb { Q } ^ { \times } \cong \bigoplus _ { p \text{ prime} } \mathbb { Z } \log p. \quad \quad \quad ( C 2 )$$

Each prime Hilbert space dimension is labelled by log p , and the integer coefficients give the shift index for that qudit species.

The classification of two-dimensional QCA on spheres was shown to be trivial in Ref. [61],

$$\pi _ { 0 } ( Q C A _ { 2 } ) \cong 0. \quad \ \ \ ( C 3 ) \quad m i$$

There is some subtlety in this statement. For QCA on a general two-dimensional manifold-that is, the lattice is obtained from a triangulation of a manifold-the classification may be nontrivial. However, all invariants relate to one-dimensional cycles in the manifold, and so belong to a lower dimensional part of the classification. For appropriate choices of topology on QCA 2 [107] (where the manifold is R 2 ), this allows all two-dimensional QCA on the plane to be connected by continuous paths.

The classification of three dimensional QCA is still open. It is conjectured that π 0 (QCA ) 3 is isomorphic to the Witt group of modular tensor categories [55, 56, 108, 109]-that is, it has the same group structure as stacking two-dimensional topological orders, modulo an equivalence relation which identifies phases which may share a gapped edge. Several nontrivial components of π 0 (QCA ) 3 have been identified, which are believed to form a subgroup [55],

$$\pi _ { 0 } ( \text{QCA} _ { 3 } ) \geq \mathbb { Z } _ { 8 } \oplus \mathbb { Z } _ { 2 }.$$

The elements of this subgroup can, indeed, be identified with two-dimensional topologically ordered phases. In fact, the ALT phases corresponding to these elements exhibit that topological order (specifically, the associated anyonic excitations) on their two dimensional surfaces [48, 55, 56, 109].

## 2. Zero dimensions

The space of zero dimensional QCA with finite Hilbert space dimension, say N , is naturally identified with the projective unitary group,

$$\text{QCA} _ { 0 } ( N ) \cong P U ( N ). \text{ \quad \ \ } ( C 5 )$$

The complication in finding the homotopy groups of QCA 0 is the stabilization by adding ancillae. That is, the N →∞ limit.

## a. Unstable invariants

For finite but large N , the homotopy groups of PU N ( ) are all known,

$$\begin{array} { c c c } \text{this} \\ \pi _ { n } ( P U ( N ) ) \cong \begin{cases} \mathbb { Z } _ { N } & \text{for $n=1$,} \\ 0 & \text{for $n<2N$ and even,} \end{cases} \quad ( \mathbb { C } 6 ) \\ \mathbb { Z } & \text{for $n<2N$ and odd.} \end{array}$$

There are even integral formulae for winding number invariants indexing the elements of these groups. Given any continuous map V : S n → PU N ( ) (where S n is the n -dimensional sphere) it is possible to find a lift ˜ : V S n \ pt → U N ( ) from the punctured sphere ( pt is the topological space with a single point) to the usual unitary group. That is, if q : U N ( ) → PU N ( ) is the quotient map from the unitary group to the projective unitaries,

we have qV ˜ = V | S n \ pt . This is possible because the punctured sphere is contractible, and so it is possible to give a smooth gauge to ˜ V away from the puncture. Defining A = (d ˜ ) ˜ V V † , where d is the exterior derivative on S n \ pt , the winding number invariants are given by [31, 130]

$$W _ { n } [ \tilde { V } ] = C _ { n } \int _ { S ^ { n } \vee \text{pt} } \text{Tr} \left [ A ^ { \wedge n } \right ], \quad \quad ( C 7 )$$

where C n is a constant independent of N and A ∧ n = A ∧ · · · ∧ A is the n -fold wedge product of A with itself. For even n we have W n = 0 , while for odd n ≥ 3 the winding number is an integer which only depends on the homotopy class of V , so we can denote them W V n [ ] . For n = 1 , only

$$w _ { 1 } [ V ] = W _ { 1 } [ \tilde { V } ] \bmod N \in \mathbb { Z } _ { N } \quad \ \ ( C 8 ) \quad \text{Th} \quad \text{trait}$$

acts as an invariant for V . Indeed, redefining ˜ V → e iθ ˜ V does not affect V , but takes W 1 → W 1 + N .

↦

## b. Stable invariants

The complication is that these winding numbers are invariant under the direct sum with an ancilla, not the tensor product. Indeed, denoting ✶ M as the identity in PU M ( ) , we have from Eq. (C7)

$$W _ { n } [ V \otimes \mathbb { 1 } _ { M } ] = M W _ { n } [ V ]. \quad \quad ( C 9 ) \quad \underset { \text{acc} } { \text{fac} }.$$

✶ More generally, for Q ∈ PU M ( ) ,

$$W _ { n } [ V \otimes Q ] = M W _ { n } [ V ] + N W _ { n } [ Q ]. \quad \ \ ( C 1 0 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

To find invariants which are stable under the addition of ancillae, we simply need to rescale the winding numbers as

$$\tilde { W } _ { n } [ V ] = C _ { n } \int _ { S ^ { n } \vee p t } \frac { \text{Tr} \left [ A ^ { \wedge n } \right ] } { \text{Tr} \left [ 1 \right ] }. \quad \ \ ( C 1 1 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

✶ These are no longer integer invariants, but rather rational invariants [modulo 1 for n = 1 ]. Then Eq. (C10) shows that ˜ W V n [ ⊗ Q ] = ˜ W V n [ ] + ˜ W Q n [ ] .

## c. Classification

At a more formal level, for all N,M ∈ N we have injective homomorphisms

$$\iota _ { N \to N M } \colon P U ( N ) & \to P U ( N M ) \quad \text{ (C12)} \\ V & \mapsto V \otimes \mathbb { 1 } _ { M } \quad \text{ the} \quad \text{ has}$$

↦

✶ which define a directed system of topological groups. QCA 0 is the direct limit

$$Q C A _ { 0 } = \varinjlim P U ( N ) \text{ \quad \ \ } ( C 1 3 )$$

with respect to this system. Strictly, we require that the multiplicity n p of every prime p in the factorization N = 2 n 2 3 n 3 . . . p n p . . . approaches infinity, n p →∞ . The homotopy groups of QCA 0 can then be computed from the induced system

$$\iota _ { N \to N M } ^ { * } \colon \pi _ { n } ( P U ( N ) ) & \to \pi _ { n } ( P U ( N M ) ) \quad \text{ (C14)} \\ W _ { n } [ V ] & \mapsto M W _ { n } [ V ]. \\ \tau \quad \cdot \quad. \quad. \quad. \quad \tau \quad. \quad.$$

↦

By rescaling the winding numbers, we find that

$$\pi _ { n } ( \text{QCA} _ { 0 } ) \cong \varinjlim _ { \ell \cdot \cdot \cdot } \pi _ { n } ( P U ( N ) )$$

$$\mathfrak { Q } \text{CA} _ { 0 } ) & \cong \lim \pi _ { n } ( P U ( N ) ) \quad \quad ( C 1 5 ) \\ & \cong \begin{cases} \mathbb { Q } / \mathbb { Z } & \text{for $n=1$}, \\ 0 & \text{for $n$ even}, \quad \quad ( C 1 6 ) \\ \mathbb { Q } & \text{for $n$ odd}. \end{cases}$$

The group operation is given by the usual addition of rational numbers, and corresponds to the tensor product of QCA.

## d. Limited ancillae

We obtained a classification by rationals because we allowed ancillae of any dimension to be appended to the system. Thus, all Hilbert space dimensions can occur in the denominator of ˜ W n . If, instead, we considered only a limited collection of Hilbert space dimensions for the ancillae, we would obtain a different classification. In fact, this classification will be a useful warm-up for the case of QCA with symmetry.

Say our available ancillae dimensions are given by

$$0 ) \quad \mathcal { R } = \{ \rho \in \mathbb { N } \ \colon \, \text{acilla of dimension } \rho \text{ exists} \}, \quad ( C 1 7 )$$

and the system is built out of degrees of freedom with these dimensions. Then R forms a submonoid (similar to a group, except without inverses) of N under multiplication, as if ✶ ρ and ✶ σ can be ancillae, then so can ✶ ρσ = ✶ ρ ⊗ ✶ σ . R also contains 1 , which corresponds to adding no ancillae.

A similar directed system to Eqs. (C12, C13) defines the limit QCA R 0 describing the stabilization of QCA by the ancillae in R . The only difference is that now the system proceeds through increasing ρ ∈ R , rather than N ∈ N . Thus, not all denominators may occur in the rescaled winding numbers. Defining the algebraic localization (not to be confused with localization in physical space)

$$\mathcal { R } ^ { - 1 } \mathbb { Z } \coloneqq \left \{ \frac { p } { q } \in \mathbb { Q } \ \colon \, p \in \mathbb { Z }, q \in \mathcal { R } \right \}, \quad \ \ ( C 1 8 )$$

then R -1 Z is an additive group (indeed, a ring) and we have

$$\pi _ { n } ( Q C A _ { 0 } ^ { \mathcal { R } } ) \cong \begin{cases} \begin{array} { c c } ( \mathcal { R } ^ { - 1 } \mathbb { Z } ) / \mathbb { Z } & \text{for $n=1$,} \\ 0 & \text{for $n$ even,} \end{array} \end{cases} ( C 1 9 )$$

## 3. QCA on manifolds

The Ω -spectrum property discussed in Appendix B assumes that QCA are defined on the integer lattice Z d . By coarse graining, any locally finite set of points in R d can be deformed to Z d , so the QCA is more properly based on the manifold R d . This distinction is important, as the classification of QCA based on the d -dimensional manifold M , QCA M , depends on the homology of M [61] and not just its dimension.

It is natural to conjecture that QCA on an orientable manifold M can all be patched together from QCA defined on the open charts in an atlas for M . That is, we construct a map

$$Q C A _ { M } \to [ M, Q C A _ { d } ] \quad \quad ( C 2 0 ) \quad \underset { \longrightarrow } { \text{res} }$$

(where [ M, QCA ] d is the space of continuous functions from M to QCA d ) by identifying the action of V ∈ QCA M near each point x ∈ M as a QCA on R d via a chart in M containing x . We suspect that a homotopy equivalence QCA M ≃ [ M, QCA ] d can be constructed in this way.

This leads to the hypothesis that, for orientable M [116],

$$\pi _ { 0 } \left [ X, Q C A _ { M } \right ] & \cong \pi _ { 0 } \left [ X \times M, Q C A _ { d } \right ] \cong h ^ { d } ( X \times M ). \quad \text{ if the} \\ \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots$$

We defer the majority of work to Ref. [107], but we remark here that Eq. (C21) reproduces the classification of QCA on two-dimensional orientable manifolds in Ref. [61]. For non-orientable manifolds, one must consider twisted generalized cohomology-that is, QCA M would be homotopy equivalent to spaces of sections in a QCA d -bundle over M [103, 104].

## Appendix D: Homotopy groups of symmetric QCA

The QCA formalism for classifying localized phases can be adapted to incorporate symmetry constraints, and thus provide a classification of symmetry-enriched ALT phases (SALT phases), and more generally of symmetryenriched MBL phases with short-ranged entanglement. In this appendix, we show that the Ω -spectrum property of QCA respects any on-site unitary group action, from which it follows that symmetric QCA also form an Ω -spectrum (Appendix D1). While the Ω -spectrum property of invertible states is linked to a classification of symmetry-enriched phases by BG -parameterized families of states ( BG is the classifying space for the symmetry group), this construction fails for QCA (Appendix D 2). Nonetheless, we identify complete invariants for all homotopy groups of symmetric QCA for d ≤ 1 , and give a presentation of the group structure (Appendices D 3, D 4).

An important observation is that the classification of symmetric QCA depends on the representation of the symmetry action [26, 93, 114]. This can be formalized within the stabilization procedure of QCA. Each added ancilla carries some representation of the symmetry group G . We do not assume that all representations of G are necessarily available. Define a countable set of representations which occur in the ancillae,

$$\mathcal { R } \coloneqq \{ \rho \in \text{Rep} ( G ) \, \colon \text{some ancillae carry } \rho \}, \quad ( D 1 )$$

where Rep( G ) is the tensor category of finite dimensional unitary representations of G . The set R is a multiplicative submonoid of Rep( G ) under tensor product, as if ρ and σ are representations carried by some ancillae, then ρ ⊗ σ is the representation of the ancilla formed by stacking both ancillae. R also contains the trivial one-dimensional representation 1 , which corresponds to appending no ancillae. This definition of R generalizes Eq. (C17), for which we can regard each ancilla as a representation of the trivial group. (One can compare R to the closely related notion of a G -universe [131], which is defined similarly, except that a G -universe is usually required to be closed under direct sum, rather than tensor product.)

We say that two QCA V 1 and V 2 which carry a group action by G are R -equivalent if ancillae from R can be appended to each of them such that they become equal, and carry the same representation. In the literature, two QCA have been called strongly equivalent if they are connected by a symmetric path modulo R -equivalence [93, 114].

## 1. Symmetric swap

In Appendix B, we argued that the swindle map S : QCA d -1 → ΩQCA d is a homotopy equivalence. Both the swindle map and the deformations which demonstrated this equivalence were defined in terms of a path swap jk ( ) t from the identity map to a swap gate between two sites with the same local Hilbert space dimension. In the context of QCA with an on-site G -action, we can repeat all the same arguments, provided we can find a path swap ρ jk ( ) t which commutes with the representation ρ .

Consider two sites with Hilbert space H ⊗H j k and a representation ρ ⊗ ρ of the group G . It is a well known fact that the representation ρ ⊗ ρ decomposes into a symmetric and alternating subspace. Indeed, we have

$$\text{swap} _ { j k } ( 1 ) \rho _ { g } \otimes \rho _ { g } = \rho _ { g } \otimes \rho _ { g } \text{swap} _ { j k } ( 1 ) \quad ( D 2 )$$

for all g ∈ G . Thus, all ρ g ⊗ ρ g become block diagonal in an eigenbasis of the swap gate, swap jk (1) . Decomposing

$$\mathcal { Y } _ { 1 } \quad \mathcal { H } _ { j } \otimes \mathcal { H } _ { k } \cong \text{Sym} ( \mathcal { H } _ { j } \otimes \mathcal { H } _ { k } ) \oplus \text{Alt} ( \mathcal { H } _ { j } \otimes \mathcal { H } _ { k } ), \quad ( \text{D3} ) \\ \mathcal { Y } _ { 1 } \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots$$

where Sym( H ⊗H j k ) is the (+1) -eigenspace of swap jk (1) and Sym( H ⊗H j k ) its ( -1) -eigenspace, we have

$$\rho \otimes \rho \cong \rho _ { \text{Sym} } \oplus \rho _ { \text{Alt} }. \quad \text{(D4)}$$

Then the path of operators

$$\text{swap} _ { j k } ^ { \rho } ( t ) = \mathbb { 1 } _ { \text{Sym} } \oplus ( e ^ { i \pi t } \mathbb { 1 } _ { \text{Alt} } ) \quad \text{(D5)}$$

is continuous, connects ✶ to the swap gate, and commutes with ρ g ⊗ ρ g for all g ∈ G and t ∈ [0 , 1] .

Then, given some on-site group action ρ on QCA d -1 ∋ v , we can construct a loop of QCA on a doubled Hilbert space

$$u ^ { \rho } _ { x, x + 1 } ( v, t ) \coloneqq & v _ { x } \text{SWAP} ^ { \rho } _ { x, x + 1 } ( t ) v ^ { \dagger } _ { x } \text{SWAP} ^ { \rho } _ { x, x + 1 } ( t ) ^ { \dagger }, \\ \dots & \quad \text{cuda: no } \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \ \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \cdots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \ dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \ Dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \Dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \ \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots$$

where SWAP ρ x,x +1 ( ) t is a product of swap ρ jk ( ) t paths, as in Eq. (6). Multiplying these together gives an instance of the swindle map

$$S ^ { \rho } ( v, t ) = & \left ( \prod _ { x \in 2 \mathbb { Z } + 1 } u _ { x, x + 1 } ^ { \rho } ( v, t ) \right ) \left ( \prod _ { x \in 2 \mathbb { Z } } u _ { x, x + 1 } ^ { \rho } ( v, t ) \right ), \quad \text{the fu} \\ \text{which reconnect the room action in the case that}$$

which respects the group action in the sense that

$$\left ( \prod _ { x } \rho _ { x } \right ) S ^ { \rho } ( v, t ) \left ( \prod _ { x } \rho _ { x } ^ { \dagger } \right ) = S ^ { \rho } ( \rho v \rho ^ { \dagger }, t ). \quad ( D 8 ) \quad \underset { \text{deep} } { \text{det} }$$

That is, the group action on the d -dimensional loop agrees with the action on the ( d -1) -dimensional QCA. Further, as the deformation from ΩQCA d to the image of S v, t ( ) in Appendix B proceeded exclusively through interleaving finite depth circuits with swap gates, we can make all the same deformations with SWAP ρ x,x +1 ( ) t . Thus, S ρ gives a G -homotopy equivalence-a homotopy equivalence which respects the G action.

This makes QCA a G -spectrum , or more specifically a (naive) Ω -spectrum with G -action. Formally, this defines a (Bredon) equivariant cohomology theory. On general grounds, it is known that the fixed point subspaces of the spectrum-QCA which commute with the group action-themselves form an Ω -spectrum [131]. (An equivariant map must map fixed points to fixed points.) We write the fixed point spaces as QCA R d , where we use R as the decoration rather than ρ to denote the stabilization by appending ancillae. Also define QCA R d for negative d by QCA R d = Ω -d QCA R 0 . Then we have

$$\pi _ { n } ( Q C A _ { d } ^ { \mathcal { R } } ) \cong \pi _ { 0 } ( Q C A _ { d - n } ^ { \mathcal { R } } ). \quad \quad ( D 9 ) \quad \text{gro}$$

More generally, there is a generalized cohomology theory

$$h _ { \mathcal { R } } ^ { d } ( X ) \coloneqq \pi _ { 0 } \left [ X, Q C A _ { d } ^ { \mathcal { R } } \right ] \quad \quad ( D 1 0 ) \quad \underset { \text{eva} } { \text{$\ }$\ } \\ \text{$\ }$\ } \text{$\ } k n c$$

which classifies parameterized families of symmetric QCA.

## 2. Symmetric QCA are not classified by BG -parameterized QCA

Invertible states are also conjectured to form an Ω -spectrum [95, 96, 102-104]. Further, it is believed that the classification of G -symmetry-enriched invertible states can be deduced just from the cohomology theory h d inv ( X ) := π 0 [ X, inv d ] and the abstract group G . In particular, denoting inv G d for the space of fixed points under the G -action, it is conjectured that

$$\pi _ { 0 } ( \text{inv} _ { d } ^ { G } ) \cong h _ { \text{inv} } ^ { d } ( B G ), \text{ \quad \ \ } ( D 1 1 )$$

where BG is the classifying space for the group G . That is, connected components of symmetric QCA are in bijection with families of not-necessarily-symmetric QCA parameterized by BG . This conjecture can has a natural interpretation in which the map in π 0 [ BG, inv d ] encodes the fusion rules of domain walls [95, 96, 102-104, 113].

The corresponding hypothesis for QCA is, in general, false:

$$\pi _ { 0 } ( \text{QCA} _ { d } ^ { \mathcal { R } } ) \not \cong h ^ { d } ( B G ). \text{ \quad \ \ } ( D 1 2 )$$

̸

Indeed, the hypothetical classification by h d ( BG ) is independent of the symmetry representation ρ , and as the strong classification of QCA (with restricted representations of ancillae) is known to depend on the representation, it cannot be the case that symmetric QCA are classified by h d ( BG ) . In fact, Ref. [26] computed the classification of U (1) -symmetric QCA in one dimension with an arbitrary symmetry representation. In this case, while the classification depends on the symmetric QCA V and the representation ρ , the invariant associated to the pair ( V, ρ ) does not depend on which stabilizing ancillae can be included, R . Thus, some of the subtleties associated to stacking ancillae are absent in this case. Regardless, we already see the ρ -independent h 1 ( BG ) cannot be a complete classification.

By explicit computation, we show

$$\pi _ { 0 } ( Q C A _ { 1 } ^ { \mathcal { R } } ) \geq h ^ { 1 } ( B U ( 1 ) ) \quad \quad ( D 1 3 )$$

̸

for any collection R of representations of U (1) . Indeed, the left hand side was found in Ref. [26]. We will repeat this classification with greater generality in Appendix D 4, but for now all we need is that the classifying group is free abelian, so no element has finite order.

On the right hand side of Eq. (D13), we can compute h 1 ( BU (1)) by solving the Atiyah-Hirzebruch spectral sequence (AHSS) [103, 104]. The AHSS is a method for evaluating generalized cohomology theories h d ( X ) given knowledge of h p (pt) for all p ≤ d +1 and the cohomology groups of X H , p ( X M ; ) . In general, it gives that h d ( X ) is a quotient of a subgroup of

$$\bigoplus _ { p = 0 } ^ { \infty } H ^ { p } ( X ; h ^ { d - p } ( \text{pt} ) ). \quad \quad ( D 1 4 )$$

Identifying the right subgroup and quotient is difficult in general, but for X = BU (1) ≃ CP ∞ and d = 1 it happens to be easy-the spectral sequence stabilizes on the E 2 page, so the subgroup is the entire direct sum and

the quotient is trivial. We have

$$h ^ { 1 } ( B U ( 1 ) ) \cong \bigoplus _ { p = 0 } ^ { \infty } H ^ { p } ( \mathbb { C P } ^ { \infty } ; h ^ { 1 - p } ( \mathbf p t ) ) \quad \quad ( D 1 5 )$$

$$\cong H ^ { 0 } ( \mathbb { C P } ^ { \infty } ; \log \mathbb { Q } _ { + } ) \oplus 0 \oplus H ^ { 2 } ( \mathbb { C P } ^ { \infty } ; \mathbb { Q } / \mathbb { Z } ) \\ \sim \partial \cap \partial \cap H ^ { 4 } / \mathbb { Q } ^ { \infty }. \, \partial \cup \cap \partial \dots \quad / \mathbb { Q } 1 \epsilon \uplus \quad \text{an}.$$

$$\cong \log \mathbb { Q } ^ { \times } \oplus \mathbb { Q } / \mathbb { Z } \oplus \mathbb { Q } \oplus \cdots. \quad ( D 1 7 )$$

$$\oplus 0 \oplus H ^ { 4 } ( \mathbb { C P } ^ { \infty } ; \mathbb { Q } ) \oplus 0 \cdots \ ( D 1 6 ) \\ \infty \cap \mathbb { m } ^ { \times } \ \Delta \cap \varPi / \varPi \ \Delta \cap \Delta \ \Delta \ \dots \quad / \mathfrak { m } 1 7 \ \mathbb { m } 1 7 \ \mathbb { m }$$

We used that BU (1) = CP ∞ and H n ( CP ∞ ; M ) = M for even n and is trivial otherwise. In particular, we have h 1 ( BU (1)) ≥ Q Z / , which is not a free group-all its elements have finite order. Thus, π 0 (QCA R 1 ) cannot contain h 1 ( BU (1)) .

In the case of invertible states, the classification by h d inv ( BG ) is justified through a decorated domain wall construction [103, 104]. We speculate that this fails for QCA because, while QCA form a spectrum, that spectrum is non-connective -it has nontrivial homotopy groups h d (pt) in arbitrarily negative degree d &lt; 0 . Thus, the putative domain wall structure encoded in h d ( BG ) includes algebraic fusion rules for domain walls which exist in negative dimension, which are not physical. This may also be why single-particle symmetric phases are not classified by h d K ( BG ) even though the phases without symmetry are classified by the K -theory spectrum h d K (which has nontrivial homotopy groups in all even degrees) [100].

It seems plausible that all symmetric QCA can be classified through a decorated domain wall construction. Certainly, all examples in the literature have this structure [26, 91, 93]. However, when correctly implemented, this approach must treat the zero-dimensional defects which terminate such constructions [113] with greater care than h d ( BG ) .

## 3. Zero dimensions with symmetry

As in the case without symmetry, the homotopy groups π n (QCA R d ) with d ≤ n are determined by π n -d (QCA R 0 ) . Thus, we make significant progress on characterizing symmetric QCA in higher dimensions through an indepth analysis of zero-dimensional QCA with symmetry. Our strategy is the same as the case without symmetry: we identify unstable invariants (at fixed Hilbert space dimension) and subsequently find a stable presentation of these invariants which does not change upon adding ancillae. Finally, we characterize the abstract classifying group.

## a. Unstable invariants

We consider a Hilbert space H with finite dimension N carrying a continuous unitary representation ρ : G → U N ( ) of an arbitrary topological group G . Then ρ is a direct sum of irreducible representations (irreps) r each with multiplicity n r ,

and we have

$$\rho = \bigoplus _ { r \, \text{irrep of } G } r \otimes 1 _ { n _ { r } }, \text{ \quad \ \ (D 1 8 ) } \\ \text{ve}$$

$$N = \sum _ { r \text{ irrep of } G } n _ { r } d _ { r }, \text{ \quad \ \ } ( D 1 9 )$$

where d r is the dimension of the irrep r . In particular, all but finitely many n r are zero.

We first describe the classification of symmetric unitaries, before explaining how this classification is modified for projective unitaries. Any family of unitaries ˜ V : S n → U N ( ) which commutes with ρ must decompose as

$$\tilde { V } = \bigoplus _ { r \, \text{irrep of } G } \mathbb { 1 } _ { d _ { r } } \otimes \tilde { V } _ { r }. \quad \text{(D2)} \\ \text{rmatio of this family which continues to com-} \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot$$

Any deformation of this family which continues to commute with ρ must maintain this structure. Thus, the families ˜ ( V 1 ⃗ θ ) and ˜ ( V 2 ⃗ θ ) are homotopic through a symmetric homotopy if and only if ˜ V 1 r ( ⃗ θ ) and ˜ V 2 r ( ⃗ θ ) are homotopic for all r . We have a complete set of invariants given by the winding numbers of each ˜ V r ,

$$W _ { n } ^ { \rho } [ \tilde { V } ] = ( W _ { n } [ \tilde { V } _ { r } ] ) _ { r \, \text{irrep of } G }. \quad \ \ ( D 2 1 )$$

With some foresight, we choose to identify W V ρ n [ ˜ ] with an element of the module

$$\det _ { \text{cts} } \quad \text{Rep} ( G ) [ \pi _ { n } ( U ( N ) ) ] = \left \{ \sum _ { \rho \in \text{Rep} ( G ) } \rho \cdot W _ { n } ^ { \rho } \\ \quad \colon W _ { n } ^ { \rho } \in \pi _ { n } ( U ( N ) ) \right \}, \quad ( D 2 2 )$$

which consists of finite formal sums of winding numbers with coefficients in Rep( G ) , regarded as a rig (a ri n g without n egatives). That is, addition in Rep( G ) is direct sum, multiplication is tensor product (which distributes over direct sums) and elements are regarded up to unitary equivalence of representations, so that addition and multiplication are commutative. We do not append the formal negatives used to define the representation ring R G . However, as π n ( U N ( )) ∈ { 0 , Z } for sufficiently large N , the module Rep( G π )[ n ( U N ( ))] is isomorphic to either 0 or R G .

Then we define the complete invariant for symmetric unitaries

$$W _ { n } ^ { \rho } [ \tilde { V } ] = \sum _ { r \text{ irrep of } G } r \cdot W _ { n } [ \tilde { V } _ { r } ]. \quad \text{(D23)}$$

QCAare properly projective unitaries, rather than unitaries. W V ρ n [ ˜ ] continues to function as a complete invariant for n ≥ 2 , as in the case without symmetry. We

may thus denote them W V ρ n [ ] , where V : S n → PU N ( ) is a family of symmetric projective unitaries which lie in the connected component of the identity. However, there are symmetric projective unitaries V which cannot be expressed as symmetric proper unitaries (requiring a larger group of n = 0 invariants), and there are loops V ( θ ) of symmetric proper unitaries which become equivalent when regarded as symmetric projective unitaries (requiring a smaller group of n = 1 invariants).

We first recall that connected components of symmetric projective unitaries [ n = 0 ] are labelled by onedimensional representations of G . Given a symmetric projective unitary V , we can arbitrarily choose a proper unitary ˜ V in the preimage of V under the quotient by global phases. The fact that V is symmetric only implies that ˜ V is symmetric up to a global phase,

$$\rho ( g ) ^ { \dagger } \tilde { V } \rho ( g ) = \chi [ V ] ( g ) \tilde { V }. \text{ \quad \ \ } ( D 2 4 )$$

Note that χ V [ ] is the same when ˜ V is multiplied by a global phase, justifying the notation χ V [ ] rather than χ V [ ˜ ] . Using that ρ is a representation, we see that χ V [ ] : G → U (1) is a one-dimensional representation of G . From Eq. (D18), we see that ˜ V must decompose as

$$\tilde { V } = \sum _ { r \text{ irrep of } G } \sigma _ { r \rightarrow \chi [ V ] r } \otimes \tilde { V } _ { r }, \quad \text{ (D25)} \quad \text{for} \\ \sum _ { r }$$

where σ r → χ V r [ ] is a matrix with a single nonzero element equal to 1, such that ˜ V maps the irrep r block of ρ to the irrep χ V [ ] r block. Any path of projective unitaries V ( ) t lifts to a path of proper unitaries ˜ ( V t ) (as the unit interval is contractible), and the space of permutations σ r → χ V r [ ] is discrete, so the permutation label χ V [ ] is an invariant for path components of projective unitaries. Any unitaries with the same χ V [ ] in the decomposition Eq. (D25) can be connected by a path, simply by connecting each ˜ V r by a path. Thus, χ V [ ] is a complete invariant.

However, not all χ V [ ] ∈ Hom( G → U (1)) =: X G can occur in a finite Hilbert space. For the ˜ V in Eq. (D25) to be unitary, we need that ∑ r σ r → χ V r [ ] is a finitedimensional permutation matrix, and so in particular χ V [ ] has finite order: there is an n such that χ V [ ]( g ) n = 1 for all g . Thus, χ V [ ] must belong to the torsion subgroup of X G (the subgroup of elements of finite order, with the group operation in X G defined pointwise), which we denote X T G . Conversely, given any χ V [ ] ∈ X T G with order n , we can take

$$\rho = \bigoplus _ { p = 0 } ^ { n - 1 } \chi [ V ] ^ { p } \quad \quad ( D 2 6 ) \quad \begin{pmatrix} S t \\ \ t h \\ \ t i c \end{pmatrix}$$

and then the unitary ∑ n -1 p =0 σ χ V [ ] p → χ V [ ] p +1 achieves the invariant χ V [ ] .

As an aside related to the discussion in Appendix D 2, we note that

$$X _ { G } ^ { T } \cong H ^ { 1 } ( G ; \mathbb { Q } / \mathbb { Z } ) \cong H ^ { 1 } ( G ; \pi _ { 1 } ( \mathbb { Q } \mathbb { A } _ { 0 } ) ) \quad ( \mathbb { D } 2 7 )$$

for the vast majority of groups G relevant to physics. For instance, a sufficient condition is that the component group G/G 0 (where G 0 is the identity component of the topological group G ) is finitely generated.

For n = 1 , we can multiply a loop of proper unitaries ˜ ( V θ ) by e imθ and change all the winding numbers, while still forming the same loop of projective unitaries after the quotient by global phases. Indeed, we find

$$W _ { 1 } ^ { \rho } [ e ^ { i m \theta } \tilde { V } ] = \sum _ { r \, \text{irrep of } G } r \cdot W _ { 1 } [ e ^ { i m \theta } \tilde { V } _ { r } ]$$

$$e ^ { i m \theta } \tilde { V } ] & = \sum _ { r \, \text{irrep of } G } r \cdot W _ { 1 } [ e ^ { i m \theta } \tilde { V } _ { r } ] \quad \text{(D28)} \\ & = \sum _ { r \, \text{irrep of } G } r \cdot ( m n _ { r } + W _ { 1 } [ \tilde { V } _ { r } ] ) \quad \text{(D29)}$$

$$& r \text{ irrep of } G \\ & = \left ( \sum _ { r \text{ irrep of } G } r \otimes \mathbb { 1 } _ { n _ { r } } \right ) \cdot m + W _ { 1 } ^ { \rho } [ \tilde { V } ] \\ & - \alpha \,. \, m \pm W ^ { \rho } [ \tilde { V } ] & / \mathbb { D } { 3 } 1 \, \mathbb { U }$$

$$= \rho \cdot m + W _ { 1 } ^ { \rho } [ \tilde { V } ],$$

where we used r · n r = r · (1 + · · · +1) = r + · · · + r = r ⊗ ✶ n r · 1 (D32) for n r ≥ 0 and the decomposition of ρ into irreps, ρ = ∑ r r ⊗ n r . Thus, we must consider the coset

$$\int _ { r } \in \mathcal { C } \, \in \, \text{and} \, \text{and} \, \text{composition} \, \in \, \text{$\rho$} \, \in \, \text{$\rho$,} \, \rho = \\ \sum _ { r } r \otimes \mathbb { 1 } _ { n _ { r } }. \, \text{Thus, we must consider the cost} \\ \text{ele-} _ { \text{$\mathbf k$ of} } \, \\ \cdot \, \quad \, \, \,$$

which is now a well-defined invariant for loops of projective unitaries V ( θ ) .

## b. Stable invariants

With the unstable invariants identified, we must now consider the stabilization process of appending ancillae. We have a directed system among symmetric projective unitary groups PU N ( ) ρ [the subset of PU N ( ) which commutes with the representation ρ ] given by

$$\iota _ { \rho \to \rho \sigma } \colon P U ( N ) ^ { \rho } & \to P U ( N M ) ^ { \rho \sigma }, \quad ( D 3 4 ) \\ V & \mapsto V \otimes 1 \\. & \quad - \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots$$

↦

✶ where we take ρ, σ ∈ R , N = dim( ) ρ , and M = dim( σ ) . The stabilized space QCA R 0 is the direct limit

$$Q C A _ { 0 } ^ { \mathcal { R } } = \varinjlim P U ( N ) ^ { \rho }. \quad \text{(D3.5)} \\ \intertext { a n o u s i o n t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h t h$$

Starting at n = 0 , we find the invariant χ V [ ] = e iω ρ [ V ] of π 0 ( PU N ( ) ρ ) is actually already stable to appending ancillae. With ω ρ [ V ] ∈ R / (2 π Z ) = ∼ U (1) defined modulo 2 π , we can see from Eq. (D24) that

Strictly, we require that the representation ρ grows through the directed system such that each representation in R eventually occurs as a tensor factor in ρ .

$$\iota _ { \rho \rightarrow \rho \sigma } ^ { * } ( \omega ^ { \rho } [ V ] ) = \omega ^ { \rho \sigma } [ V \otimes 1 ] = \omega ^ { \rho } [ V ] \quad \text{(D36)}$$

and more generally

$$\omega ^ { \rho \sigma } [ V \otimes Q ] = \underset { \dots } { \omega ^ { \rho } } [ V ] + \omega ^ { \sigma } [ Q ]. \quad \underset { \dots } { \omega ^ { \dagger } } \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots$$

To compute the action of the induced maps on the remaining unstable invariants, ι ∗ ρ → ρσ ( W V ρ n [ ]) = W ρσ n [ V ⊗ ✶ ] , we suppose that we have Clebsch-Gordan multiplicities C k rr ′ ∈ Z such that, for r, r , k ′ irreducible,

$$r r ^ { \prime } = \sum _ { k } C _ { k } ^ { r r ^ { \prime } } k \in \text{Rep} ( G ), \quad \text{(D38)} \quad \text{(Woo)}$$

where C rr ′ k k = k + · · · + k = C rr ′ k ⊗ k . Taking

$$\text{where $C_{k}^{r}^{r}k=k+\cdots+k=\mathbb{1}_{C_{k}^{r}^{r}^{\prime}\otimes k.$ Taking } \quad \text{$two\quad}$ this} \\ \rho = \sum _ { r } n _ { r } ^ { \rho } r \quad \text{and} \quad \sigma = \sum _ { r ^ { \prime } } n _ { r ^ { \prime } } ^ { \sigma } r ^ { \prime }, \quad \text{(D3)} \\. \quad \dots \quad. \quad. \. \.$$

we have, with summation implied,

$$\rho \sigma = n _ { r } ^ { \rho } n _ { r ^ { \prime } } ^ { \sigma } C _ { k } ^ { r r ^ { \prime } } k. \quad \quad ( D 4 0 ) \quad \ \.$$

Then V ⊗ ✶ decomposes into the irreducible subspaces of ρσ (possibly after multiplication by a permutation of the irrep subspaces) as

$$V \otimes 1 \cong \bigoplus _ { k } \left [ \mathbb { 1 } _ { d _ { k } } \otimes \left ( \bigoplus _ { r, r ^ { \prime } } V _ { r } \otimes \mathbb { 1 } _ { n ^ { \sigma } _ { r ^ { \prime } } } C ^ { r r ^ { \prime } } _ { k } \right ) \right ], \quad ( D 4 1 ) \quad & \text{The $\cosh$ an} \\ \tilde { w } _ { 1 } ^ { \rho } [ \mathbb { I } \\ \text{and the invariant becomes} \\ \cdot$$

and the invariant becomes

$$W _ { n } ^ { \rho \sigma } [ V \otimes \mathbb { 1 } ] = \sum _ { k, r, r ^ { \prime } } k \cdot W _ { 1 } [ V _ { r } \otimes \mathbb { 1 } _ { n _ { r ^ { \prime } } ^ { \sigma } C _ { k } ^ { r r ^ { \prime } } } ] \quad ( \Gamma$$

$$\text{invariant becomes} ^ { \text{ (} \sigma ) } _ { \iota } [ V \otimes 1 ] & = \sum _ { k, r, r ^ { \prime } } k \cdot W _ { 1 } [ V _ { r } \otimes 1 _ { n _ { r ^ { \prime } } ^ { \sigma } C _ { k } ^ { r r ^ { \prime } } } ] \quad ( D { \text{4} } { \text{as} } \quad \text{iso} \\ & = \sum _ { k, r, r ^ { \prime } } ( n _ { r ^ { \prime } } ^ { \sigma } C _ { k } ^ { r r ^ { \prime } } k ) \cdot W _ { 1 } [ V _ { r } ] \quad ( D { \text{4} } { \text{as} } \quad \text{in} \quad \text{$\quad$W_{n}$} \\ & = \sum _ { r } ( \sigma r ) \cdot W _ { 1 } [ V _ { r } ] \quad ( D { \text{4} } { \text{and} } \quad \text{are} \\ & = \sigma \cdot W _ { n } ^ { \rho } [ V ]. \quad ( D { \text{4} } { \text{and} } \quad \text{as} \quad \text{$\quad$V$} \\ \text{ogous calculation produces the generalization} \quad \text{$\quad$am}.$$

An analogous calculation produces the generalization

$$\sup _ { \dots } \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdot$$

Then the induced directed system is

$$\iota _ { \rho \to \rho \sigma } ^ { * } \colon \text{Rep} ( G ) [ \pi _ { n } ( U ( N ) ) ] & \to \text{Rep} ( G ) [ \pi _ { n } ( U ( N M ) ) ] \quad \text{also} \\ W _ { n } ^ { \rho } [ V ] & \mapsto \sigma \cdot W _ { n } ^ { \rho } [ V ]. \quad ( D 4 7 ) \\ \intertext { W o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o w o } = \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } +\partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 0 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \odot _ { 1 } + \partial _ { 1 } + \partial _ { 1 } + \odot _ { 1 } + \partial _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odot _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odot _ { 1 } + \odots _ { 1 } + \odot _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odot _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odts _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odts _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odts _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odts _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odelta _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odelta _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odts _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odets _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odots _ { 1 } + \odets _ { 1 } + \odots _ { 1 } + \od$$

↦

We recover Eq. (C14) by taking G = 0 to be the trivial group.

Thus, we can identify a stable invariant by formally dividing W V ρ n [ ] by the representation ρ ,

$$\tilde { W } _ { n } ^ { \rho } [ V ] = \frac { W _ { n } ^ { \rho } [ V ] } { \rho } \in \mathcal { R } ^ { - 1 } \text{Rep} ( G ) [ \pi _ { n } ( U ( N ) ) ], \quad ( D 4 8 ) \quad \text{This} \ \,.$$

which then obeys ˜ W ρσ n [ V ⊗ ✶ ] = ˜ W V ρ n [ ] . Here, R -1 Rep( G π )[ n ( U N ( ))] is the algebraic localization of the module Rep( G π )[ n ( U N ( ))] , defined in the next section.

For n = 1 , the invariant is a coset

$$\tilde { w } _ { 1 } ^ { \rho } [ V ] = \tilde { W } _ { 1 } ^ { \rho } [ V ] + 1 \cdot \mathbb { Z } \in \mathcal { R } ^ { - 1 } \text{Rep} ( G ) [ \pi _ { 1 } ( U ( N ) ) ] / ( 1 \cdot \mathbb { Z } ). \\ ( D 4 9 )$$

## c. Classification

If M is a module over R and S ⊆ R is a multiplicatively closed subset of R , then the algebraic localization of M is the set S -1 M := ( M × S / ) ∼ , where the equivalence relation is

$$( m _ { 1 }, s _ { 1 } ) \sim ( m _ { 2 }, s _ { 2 } ) \iff \exists \sigma \in S \text{ such that} \\ \sigma \cdot ( s _ { 1 } \cdot m _ { 2 } - s _ { 2 } \cdot m _ { 1 } ) = 0 \in M. \quad ( D 5 0 ) \\ \text{W-out with to a problem alone at} \, ( m _ { 1 } \cdot \partial \, m / a \, \partial \, m / a \, \partial \, \text{and} \, \text{coin}$$

We write the equivalence class of ( m,s ) as m/s and equip this set with the addition operation and R -action

$$\frac { m _ { 1 } } { s _ { 1 } } + \frac { m _ { 2 } } { s _ { 2 } } = \frac { s _ { 2 } \cdot m _ { 1 } + s _ { 1 } \cdot m _ { 2 } } { s _ { 1 } s _ { 2 } }, \quad \text{(D51)}$$

$$\nu _ { 1 } \quad \nu _ { 2 } \quad \rho \cdot \frac { m } { s } = \frac { \rho \cdot m } { s }.$$

The formal ratios ˜ W V ρ n [ ] are properly defined as elements of the localization R -1 Rep( G π )[ n ( U N ( ))] , which from Eqs. (D46, D51) gives

$$\tilde { W } _ { n } ^ { \rho \sigma } [ V \otimes Q ] = \tilde { W } _ { n } ^ { \rho } [ V ] + \tilde { W } _ { n } ^ { \sigma } [ Q ]. \quad ( D 5 3 ) \\.. \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot$$

The one-dimensional stable symmetric winding number is an additive coset in

$$\tilde { w } _ { 1 } ^ { \rho } [ V ] \in \mathcal { R } ^ { - 1 } \text{Rep} ( G ) [ \pi _ { 1 } ( U ( N ) ) ] / ( 1 \cdot \pi _ { 1 } ( U ( N ) ) ). \ \text{(D54)} \\ \bullet \ \cdot \ \cdot \ \, \, \,$$

Again, as π n ( U N ( )) ∈ { 0 , Z } , these groups are always isomorphic to either the trivial group or R -1 R G .

Note that Eq. (D50) implies that ˜ W ρ 1 n [ V ] = ˜ W ρ 2 n [ Q ] if and only if there are representations σ ⊗ ρ 1 and σ ⊗ ρ 2 in R such that the unstable winding numbers

$$W _ { n } [ V \otimes 1 _ { \sigma \rho _ { 2 } } ] = \sigma \rho _ { 2 } \cdot W _ { n } [ V ] = \sigma \rho _ { 1 } \cdot W _ { n } [ Q ] = W _ { n } [ Q \otimes 1 _ { \sigma \rho _ { 1 } } ] \\ ( D 5 5 ) \\ \text$$

are equal. Thus, V and Q are strongly equivalent, exactly as we want.

We have identified a complete set of stable invariants among the elements of the localization R -1 R G , but not all elements in R -1 R G necessarily correspond to achievable invariants with a given R . The subset which can occur is that where the numerator involves irreps which also occur in the denominator. We define the additive subgroup of R -1 R G

̸

$$\begin{array} { c }. \\ \rho = s + r \text{ for some } s \in \text{Rep} ( G ) \end{array} \right \}. \ \ ( D 5$$

$$\text{ival} \quad & \mathcal { K } ^ { \mathcal { R } } = \left \{ \frac { \sum _ { r } \, \text{irrep} \, n _ { r } r } { \rho } \in \mathcal { R } ^ { - 1 } R _ { G } \, \colon \, n _ { r } \neq 0 \, \Longrightarrow \\ \text{ally$$

This is precisely the subgroup of winding number invariants which can occur.

To conclude, the homotopy groups of symmetric zerodimensional QCA are

$$\pi _ { n } ( Q C A _ { 0 } ^ { \mathcal { R } } ) \cong \begin{cases} \begin{array} { c c } \leq X _ { G } ^ { T } & n = 0, \\ \mathcal { K } ^ { \mathcal { R } } / \mathbb { Z } & n = 1, \\ 0 & n \geq 2 \text{ even}, \\ \mathcal {$$

where by Z we mean the subgroup generated by 1 ∈ K R . We have also not precisely characterized the restriction on π 0 (QCA R 0 ) ≤ X T G imposed by the limited ancillae in R , as our focus is on the less familiar higher homotopy groups.

↦

As an example, we consider G = U (1) . Then all irreducible representations are one-dimensional, and given by U (1) ∋ e iθ → e imθ (in multiplicative notation). Identifying this irrep with the polynomial x m , the representation ring is isomorphic to the ring of integer Laurent polynomials in x ,

$$R _ { U ( 1 ) } \cong \mathbb { Z } [ x, x ^ { - 1 } ]. \text{ \quad \ \ } ( D 5 8 )$$

We have that X U (1) ∼ = Z , so X T U (1) ∼ = 0 is trivial. Let us suppose that the monoid of ancillae R contains all finite dimensional representations of U (1) , so R is the full representation rig and thus isomorphic to Laurent polynomials with non-negative coefficients,

$$\mathcal { R } \cong \mathbb { N } [ x, x ^ { - 1 } ]. \quad \quad ( D 5 9 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

Then K R is the additive group of rational polynomials with positive denominators where all powers which occur in the numerator also occur in the denominator,

$$\mathcal { K } ^ { \mathcal { R } } & \cong \left \{ \frac { p ( x ) } { q ( x ) } \in \text{Frac} \, \mathbb { Z } [ x ] \, \colon \, q ( x ) \in \mathbb { N } [ x ], \quad \begin{array} { c c } \text{for all} \\ \text{on th} \\ \text{Th} \\ \text{Eq. (} \\ \text{are } p ^ { \prime } \\ \text{$\mathbb{ }$} \, \text{$\mathbb{ }$} \end{array} \right \}$$

We used that one can divide the numerator and denominator of a ratio of Laurent polynomials by the most negative power of x in either to obtain a ratio of usual polynomials. Then the homotopy groups of QCA R 0 with this set of ancillae is given by Eq. (D57). We note that the rational polynomial structure which appeared in Ref. [26] is again occurring here, as it is a consequence of the representation theory of U (1) .

## 4. One dimension with symmetry

The study of the connected components π 0 (QCA R 1 ) has been considered in Refs. [26, 93, 114], where complete classifications were found for U (1) and finite cyclic groups of prime order. In this appendix, we use the representation theory technology of the previous section to generalize these classifications to arbitrary topological groups G and countable collections of ancillae representations R .

## a. Unstable invariants

Two one-dimensional QCA V 1 and V 2 are path equivalent if and only if V V 1 † 2 is a finite depth circuit [58, 59]. Leveraging this fact simplifies the definition of invariants and proofs of classifications. Indeed, all one-dimensional

QCA can be expressed as circuits of isometries [132]. That is, possibly after coarse graining by blocking sites together, any QCA V on a one-dimensional chain with sites x and local Hilbert space dimensions N x is equal to a product

$$V = \left ( \prod _ { x \in 2 \mathbb { Z } + 1 } u _ { x, x + 1 } ^ { \dagger } \right ) \left ( \prod _ { x \in 2 \mathbb { Z } } u _ { x, x + 1 } \right ) \quad ( \text{D61} )$$

where

$$u _ { x, x + 1 } \in U ( N _ { x } N _ { x + 1 } \to N _ { x } ^ { \prime } N _ { x + 1 } ^ { \prime } ). \quad \text{(D62)}$$

That is, V can be written as a finite depth circuit where the internal legs of the circuit do not necessarily have the same dimension as the external legs-the gates of the circuit are isometries, rather than unitaries. However, for the entire circuit to be unitary, we must have N N x x +1 = N N ′ x ′ x +1 .

Then V commutes with the on-site symmetry ∏ x ρ x if and only if

$$\mathfrak { r } \quad u _ { x, x + 1 } \rho _ { x } ( g ) \rho _ { x + 1 } ( g ) = v _ { x } ( g ) v _ { x + 1 } ( g ) u _ { x, x + 1 } \quad ( \mathrm D 6 3 )$$

for all g ∈ G and x ∈ Z , where v x ( g ) ∈ U N ( ′ x ) is a unitary on the internal leg [132].

The product of two consecutive v x v x +1 is, by Eq. (D63), unitarily equivalent to the representation ρ ρ x x +1 . However, this only implies that v x and v x +1 are projective representations of G . The classifications of Ref. [26, 93, 114] rely on the following observation: V is equal to a finite depth circuit with symmetric unitary gates (which implies they are connected by a symmetric path) if and only if there are functions ϕ x : G → U (1) such that e iϕ x v x is unitarily equivalent to ρ x for all x . That is, the projective representation v x can be lifted to an ordinary unitary representation, and that representation is precisely ρ x . A more concise statement of the condition is that v x and ρ x are isomorphic when both regarded as projective representations.

Indeed, suppose that

$$e ^ { i \phi _ { x } ( g ) } w _ { x } v _ { x } ( g ) w _ { x } ^ { \dagger } = \rho _ { x } ( g ) \text{ \ for all } g \in G. \text{ \ \ (D64)}$$

From the fact that v x v x +1 is unitarily equivalent to ρ ρ x x +1 , we conclude that ϕ x = -ϕ x +1 . Further, N x = N ′ x as w x is unitary, so we can write V as a circuit of unitary gates

$$V = \left ( \prod _ { x \in 2 \mathbb { Z } + 1 } u _ { x, x + 1 } ^ { \prime \dagger } \right ) \left ( \prod _ { x \in 2 \mathbb { Z } } u _ { x, x + 1 } ^ { \prime } \right ) \quad \text{(D65)}$$

with

$$\begin{array} { l l } { { 1 } { \dots } } & { { u ^ { \prime } _ { x, x + 1 } = w _ { x } w _ { x + 1 } u _ { x, x + 1 } } } \\ { { \dots } } & { { = ( e ^ { i \phi _ { x } } w _ { x } ) ( e ^ { i \phi _ { x + 1 } } w _ { x + 1 } ) u _ { x, x + 1 } } } & { ( D 6 6 ) } \end{array}$$

such that

$$u _ { x, x + 1 } ^ { \prime } \rho _ { x } \rho _ { x + 1 } & = w _ { x } w _ { x + 1 } v _ { x } v _ { x + 1 } u _ { x, x + 1 } & ( D 6 7 ) \quad \text{ dimer} \\ & = ( e ^ { i \phi _ { x } } w _ { x } v _ { x } w _ { x } ^ { \dagger } ) ( e ^ { i \phi _ { x + 1 } } w _ { x + 1 } v _ { x + 1 } w _ { x + 1 } ^ { \dagger } ) \\ & \quad \times w _ { x } w _ { x + 1 } u _ { x, x + 1 } & ( D 6 8 ) \\ & = \rho _ { x } \rho _ { x + 1 } u _ { x, x + 1 } ^ { \prime } & ( D 6 9 )$$

is symmetric.

Conversely, if V can be written as a circuit with symmetric unitary gates, as in Eq. (D65), then we must have that

$$\left ( \prod _ { x \in 2 \mathbb { Z } + 1 } u _ { x, x + 1 } ^ { \prime \dagger } \right ) \left ( \prod _ { x \in 2 \mathbb { Z } } u _ { x, x + 1 } ^ { \prime } u _ { x, x + 1 } ^ { \dagger } \right ) = 1. \quad \begin{array} { c c } \text{me } v \vdots \\ \text{(D70)} \end{array} \quad \begin{array} { c } \text{me } v \vdots \\ \text{(de)} \\ \text{We see by, for instance, the fact that the operator entan-} \end{array}$$

We see by, for instance, the fact that the operator entanglement across the ( x, x +1) bond in this product must vanish, that

$$u _ { x, x + 1 } ^ { \prime } u _ { x, x + 1 } ^ { \dagger } = w _ { x } w _ { x + 1 } \quad \quad ( D 7 1 ) \quad \ m u$$

must be a tensor product of isometries w x ∈ U N ( x → N ′ x ) . Further, as the whole circuit is the identity, we have w w † x x = ✶ and thus N ′ x ≥ N x for all x . However, N N x x +1 = N N ′ x ′ x +1 , and these conditions can only be simultaneously satisfied if N x = N ′ x for all x , and thus w x is a unitary gate.

Then we have

$$v _ { x } v _ { x + 1 } u _ { x, x + 1 } = u _ { x, x + 1 } \rho _ { x } \rho _ { x + 1 } & & ( D 7 2 )$$

$$\begin{array} { c } \cdot \\ = w _ { x } w _ { x + 1 } u ^ { \prime } _ { x, x + 1 } \rho _ { x } \rho _ { x + 1 } \end{array} \quad \text{(D7 3)}$$

$$= w _ { x } w _ { x + 1 } \rho _ { x } \rho _ { x + 1 } u ^ { \prime } _ { x, x + 1 } \quad \text{(D7 4)} \quad \text{So}$$

$$= ( w _ { x } \rho _ { x } w _ { x } ^ { \dagger } ) ( w _ { x + 1 } \rho _ { x + 1 } w _ { x + 1 } ^ { \dagger } ) u _ { x, x + 1 }, \quad \text{ an} ( D 7 5 ) \quad \text{ an}$$

so that v x = e -iϕ x w ρ w x x † x , as we needed.

Thus, an invariant for symmetric QCA in one dimension must somehow measure the difference between v x and ρ x . We split this into two parts. First, we need v x to lift to a usual unitary representation. It is well known that the obstruction to lifting a projective unitary representation to a unitary representation is captured by group cohomology with coefficients in U (1) , so there is an

$$\omega _ { x } \in H ^ { 2 } ( G ; U ( 1 ) ) \text{ \quad \ \ } ( \text{D7} ) \text{ \ \ } s a i$$

such that v x lifts to a unitary representation if and only if ω x = 0 . (We use an additive notation for ω x .) Further, e -iω x v x always lifts to a unitary representation, say ˜ v x . The unitary representation ˜ v x is unique up to multiplication by a one dimensional representation χ x : G → U (1) . Thus, we can identify the isomorphism class of ˜ v x , with a coset in the representation rig,

$$\tilde { v } _ { x } X _ { G } \in \text{Rep} ( G ) / X _ { G }, \text{ \quad \ \ ( D 7 7 ) \quad v _ { 0 } / }$$

where X G is the multiplicative subgroup of units (one dimensional representations) in Rep( G ) . Then V is a symmetric circuit if and only if

$$( \omega _ { x }, \tilde { v } _ { x } X _ { G } ) = ( 0, \rho _ { x } X _ { G } ) \text{ \quad \ \ } ( D 7 8 )$$

for all x .

## b. Stable invariants

While the ω x invariant is stable to stacking ancillae, the ˜ v x X G invariant is not, and both apparently depend on position x . Stabilizing the ˜ v x X G invariant through the algebraic localization construction fixes both of these issues.

We modify the ˜ v x X G invariant to

$$\frac { v _ { x } } { \rho _ { x } } X _ { G } \in \mathcal { R } ^ { - 1 } \text{Rep} ( G ) / X _ { G }, \text{ \quad \ \ } ( \text{D79} )$$

which we now regard as a multiplicative group via the multiplication in Rep( G ) .

The quotient by X G may refer to either the quotient by X / G 1 ⊆ R -1 Rep( G ) or performing the X G quotient first, then the localization. The results are isomorphic.

When appending additional ancillae to site x with representation σ x , V ⊗ ✶ acts trivially on the ancilla, so v x just becomes v x ⊗ σ x . Then v x σ /ρ σ x x x = v x /ρ x is unchanged.

Further, because v x v x +1 = ρ ρ x x +1 ∈ R , we have

$$\omega _ { x } + \omega _ { x + 1 } = 0 \ \text{ and } \ \frac { \tilde { v } _ { x } } { \rho _ { x } } \frac { \tilde { v } _ { x + 1 } } { \rho _ { x + 1 } } = 1. \quad \ \text{(D80)}$$

So the invariant for x +1 is just the inverse of that for x , and has no new information. Thus, we can just focus on x = 0 , and henceforth drop the x index.

## c. Classification

We have a complete set of invariants in H 2 ( G U ; (1)) × R -1 Rep( G /X ) G . Again, not all ( ω, vX ˜ G /ρ ) can actually be achieved with a fixed R . We focus primarily on ˜ vX /ρ G .

The only numerators ˜ v which may occur in the invariant ˜ vX /ρ G are those which occur as tensor factors in ρ for some ρ ∈ R . That is, the numerators belong to the saturation

$$\text{nly} \quad \hat { \mathcal { R } } = \{ v \in \text{Rep} ( G ) \ \colon \, \exists \rho \in \text{Rep} ( G ) \\ \text{er}, \quad \text{such that } v \otimes \rho \in \mathcal { R } \}. \quad ( D 8 1 ) \\ \hat { v } _ { x }. \quad \hat { \nu } \quad \hat { \nu }$$

Note that ˆ R is closed under multiplication by R , so the localization R -1 ˆ R is well defined and a multiplicative group under the natural multiplication operation of fractions.

All the elements of ˆ R may occur as invariants. Suppose v /ρ 0 0 ∈ R -1 ˆ R . Then, as v 0 belongs to ˆ R , there are

representations v 1 ∈ Rep( G ) and ρ 1 ∈ R and a unitary u ′ 0 1 , such that

$$v _ { 0 } \otimes v _ { 1 } u _ { 0, 1 } ^ { \prime } = u _ { 0, 1 } ^ { \prime } \rho _ { 1 }. \text{ \quad \ \ \ } ( D 8 2 )$$

We tensor in ρ 0 to get

$$v _ { 0 } \otimes ( v _ { 1 } \otimes \rho _ { 0 } ) u _ { 0, 1 } = u _ { 0, 1 } \rho _ { 0 } \otimes \rho _ { 1 }. \quad ( D 8 3 )$$

Then we define a QCA V which commutes with the representation ∏ x ∈ 2 Z ρ 0 x ⊗ ρ 1( x +1) by

$$V = \left ( \prod _ { x \in 2 \mathbb { Z } + 1 } u _ { x, x + 1 } ^ { \dagger } \right ) \left ( \prod _ { x \in 2 \mathbb { Z } } u _ { x, x + 1 } \right ), \quad ( D 8 4 ) \quad \text{there} \quad \text{there} \quad \text{Unde} \quad \text{ }$$

where u x,x +1 is a translation of u 0 1 , from Eq. (D83). The invariant for the QCA V is, by construction, ( v X 0 G ) /ρ 0 , or its inverse ( v ρ X 1 0 G ) /ρ 1 .

With regard to the invariant described by Eq. (D76), let us briefly address whether all elements of H 2 ( G,U (1)) can actually be achieved by a QCA [recall that this was not necessarily the case for the H 1 ( G,U (1)) invariant discussed in d = 0 , for infinite groups G ]. This invariant characterizes the obstruction in lifting the projective representation v x to a unitary representation. For finite groups, all such invariants can be achieved. Taking the local Hilbert space to be the span of | g ⟩ with g ∈ G , and given any 2-coycle α : G × G → U (1) , we define

$$v _ { 0 } ( g ) = \rho _ { \text{reg} } ( g ) D ( g ) \text{ \quad \ \ } ( D 8 5 )$$

where ρ reg is the regular representation and D g ( ) is the diagonal matrix with elements ⟨ h D g | ( ) | h ⟩ = α g, h ( ) . Then v 0 is known to be a projective representation with obstruction [ α ] ∈ H 2 ( G U ; (1)) . By a similar construction to the previous paragraph, we can construct a QCA (indeed, a finite depth circuit) from the unitary gate u 0 1 , which conjugate ρ reg ⊗ ρ reg to v 0 ⊗ v 1 , where v 1 has obstruction -[ α ] . This gate exists by the vanishing of the obstruction for v 0 ⊗ v 1 . We leave open the precise restrictions on ω ∈ H 2 ( G U ; (1)) imposed by limited ancillae representations or by G being infinite.

The full classifying group is then a subgroup

$$\pi _ { 0 } ( \text{QCA} _ { 1 } ^ { \mathcal { R } } ) \leq H ^ { 2 } ( G ; U ( 1 ) ) \oplus \log \mathcal { R } ^ { - 1 } \hat { \mathcal { R } } / X _ { G }, \ \ ( D 8 6 ) \quad \text{represent}$$

with only the first summand possibly allowing for this to be a proper subgroup. That is, the projection to the second summand is surjective. We express the second summand with a logarithm just to make the group operation addition. The induced homomorphism from evaluation on a product state, e ∗ 0 , just corresponds to projecting to the first summand, and thus reproduces the group cohomology classification of invertible states in one dimension [66, 93].

The case of G = U (1) again serves as a simple example [26, 27]. We take R ∼ = N [ x, x -1 ] to consist of all representations of U (1) . Then ˆ = R ∼ N [ x, x -1 ] as well, and the algebraic localization is isomorphic to the multiplicative group of rational polynomials with positive coefficients,

$$\mathcal { R } ^ { - 1 } \hat { \mathcal { R } } \cong \text{Frac} \, \mathbb { N } [ x ] ^ { \times } \,. \quad \quad ( D 8 7 )$$

The group of units X U (1) consists of monomials, x m . Expressing the quotient in the additive notation, this is

$$& ) \quad \log \mathcal { R } ^ { - 1 } \hat { \mathcal { R } } / X _ { U ( 1 ) } \cong \log \text{Frac} \, \mathbb { N } [ x ] ^ { \times } / \mathbb { Z } \log x. \quad ( D 8 8 ) \\ &$$

As H 2 ( U (1); U (1)) ∼ = 0 is trivial, this is the complete classification of U (1) -symmetric QCA in one dimension.

It has also been observed that including the regular representation for a finite group G in R , ρ reg ∈ R , reduces the classification to just H 2 ( G U ; (1)) ⊕ log Q × [93]. Indeed, in the current presentation it is clear that this is because R G has zero divisors when G is a finite group. If ρ reg is the regular representation, then in the representation rig we have ρ reg σ = dim( σ ρ ) reg , and so

$$\rho _ { \text{reg} } [ \sigma \dim ( \sigma ^ { \prime } ) - \sigma ^ { \prime } \dim ( \sigma ) ] = 0, \quad \ ( D 8 9 )$$

for any σ, σ ′ ∈ Rep( G ) . Thus, σ/σ ′ = dim( σ / ) dim( σ ′ ) ∈ R -1 ˆ R , and any two representations of the same dimension define the same element of the algebraic localization, which is then isomorphic to ratios of positive integers, R -1 ˆ R /X G ∼ = Q × . (More correctly, this should itself be a localization of the natural numbers, dim( R ) -1 dim( ˆ ) R .)

## 5. Functorality

Any classification of symmetric objects should, ideally, be compatible with forgetting the symmetry . That is, it should be possible to map the classification of symmetric objects into the classification of objects without symmetry in a consistent way. This is generalized by the notion of the classification being functorial [103].

Suppose that we have a homomorphism of topological groups

$$\phi \colon H \to G, \text{ \quad \ \ } ( \mathbb { D } 9 0 )$$

and a representation ρ : G → U N ( ) of G . Then ρϕ is a representation of H ,

$$\rho \phi \colon H \stackrel { \phi } { \rightarrow } G \stackrel { \rho } { \rightarrow } U ( N ). \text{ \quad \ \ } ( D 9 1 )$$

↦

Further, we have that ( ρ 1 ⊕ ρ 2 ) ϕ ≃ ρ ϕ 1 ⊕ ρ ϕ 2 and ( ρ 1 ⊗ ρ 2 ) ϕ ≃ ρ ϕ 1 ⊗ ρ ϕ 2 , so that ρ → ρϕ defines an induced rig homomorphism

$$\phi _ { \text{Rep} } ^ { * } \colon \text{Rep} ( G ) \to \text{Rep} ( H ). \quad \quad ( \text{D92} )$$

Note that the order of G and H are reversed. The map is contravariant.

What's more, any symmetric V ∈ QCA R d (or more generally a QCA subject to some G action, not necessarily symmetric) can also be regarded as being symmetric with

respect to the H action ρϕ . Thus, we have a continuous inclusion map

$$\phi _ { \text{QCA} } ^ { * } \colon \text{QCA} _ { d } ^ { \mathcal { R } } \to \text{QCA} _ { d } ^ { \mathcal { R } ^ { \prime } } \text{ \quad \ \ (D 9 3 ) }.$$

where ϕ ∗ Rep ( R ⊆ R ) ′ . In turn, this continuous map induces homomorphisms

$$\phi _ { n } ^ { * } \colon \pi _ { n } ( \text{QCA} _ { d } ^ { \mathcal { R } } ) \to \pi _ { n } ( \text{QCA} _ { d } ^ { \mathcal { R } ^ { \prime } } ). \quad ( \text{D94} ) \quad.$$

In Appendices D3 and D4, we classified homotopy groups of symmetric QCA with d -n ≤ 1 . More concretely, we constructed isomorphisms between π n (QCA R d ) and certain groups defined in terms of Rep( G ) . One should want these isomorphisms to be functorial-the induced map ϕ ∗ n should be reproduced by applying ϕ ∗ Rep (or more correctly the induced map on the localization or quotient) to the invariant in R -1 Rep( G ) .

It is not difficult to see that this is the case. Indeed, as the map ϕ ∗ QCA does not change the QCA V , only the representation, carrying through the calculation of the invariant just involves replacing all the representations ρ which occur with ρϕ , which is precisely the action of ϕ ∗ Rep .

An important example is completely forgetting the symmetry, corresponding to the inclusion of the identity, ϕ : 1 → G . Then we have that all representations reduce to the identity operator ρϕ = ✶ dim( ρ ) . In the rig Rep( G ) , this corresponds to the dimension map,

$$\dim \colon \rho \mapsto \dim ( \rho ), \quad \quad ( D 9 5 ) \quad S ( \begin{matrix} \\ \\ \text{pli} \end{matrix}$$

↦

where the right hand side is to be interpreted as the sum of dim( ρ ) ∈ N copies of the trivial one-dimensional representation 1 . For G = U (1) , identifying Rep( U (1)) with N [ x, x -1 ] , the dimension map just corresponds to evaluating the rational polynomial invariant f ( x ) at x = 1 [26],

$$\dim ( f ( x ) ) = f ( 1 ). \text{ \quad \ \ } ( D 9 6 ) \text{ \ \ an}$$

## 6. Lattice symmetries

It is believed that the generalized cohomology classification of invertible states also captures lattice symmetries and antiunitary symmetries [103, 104]. We also believe this to be true for QCA, and in this section we outline the structure that this should take. We leave detailed explorations to Ref. [107] and future work.

In the picture of Appendix C3, a QCA on some d -dimensional orientable manifold M can be identified by a map from M to QCA d

$$V \in Q C A _ { M } \longleftrightarrow v \, \colon M \to Q C A _ { d }. \quad \text{ (D9 7)} \quad \text{ ever}.$$

More precisely, we suspect QCA M and [ M, QCA ] d to be homotopy equivalent. Lattice symmetries can be treated in this picture as actions of the symmetry group G on the base manifold M [116]. Then, identifying the QCA V with the map v , the QCA is symmetric if

$$v ( x ) = \rho _ { g } v ( g \cdot x ) \rho _ { g } ^ { \dagger } \text{ \quad \ \ } \text{(D98)}$$

for all g ∈ G and x ∈ M . Then QCA with lattice symmetries would be fixed points of the mapping space

$$[ M, Q C A _ { d } ] ^ { \mathcal { R } } \, \quad \, \quad \, ( D 9 9 )$$

under the simultaneous action of G on QCA d and M . The action on QCA d will generically not be on-site. For instance, a reflection should also act to reflect the QCA v x ( ) .

Abstractly, we are lead to equivariant Bredon cohomology, which classifies such fixed points of mapping spaces when the target space is an Ω -spectrum [131]. The cohomology may, in general, need to be twisted when either M is not orientable or, as would be more common for physics, the symmetry action is orientation reversing (involves reflections) or antiunitary.

## Appendix E: Models

The swindle map S v, t ( ) is used in the main text and Appendix B as a theoretical tool. In fact, it also provides a useful construction for explicit models of d -dimensional driven (S)ALT phases, given a QCA (or family of QCA) v in one dimension lower. If v is nearest-neighbor, so is S v, t ( ) (in the ℓ ∞ distance-it may involve diagonal couplings in the lattice). Even better, a slightly modified model only involves coupling two qubits at a time in one dimension. Thus, these models are well suited to experimental realization in a variety of quantum simulation platforms [63-65].

In Appendix E 1, we define the modified swindle map which only couples disjoint pairs of qubits at any time, and show how disorder may be introduced into the swindle map without spoiling solubility [24]. Then we explore the particular example of the quasiperiodic energy pump (QP pump) [31, 32, 53] in Appendix E 2. This is a one-dimensional U (1) -SALT phase with two drive tones which, thus far, has not been experimentally realized. The model we obtain through the swindle construction is particularly simple, paving the way for experiments on this phase.

## 1. Modified swindle

If v ∈ QCA d -1 has a range r = 1 ( vAv † is localized within one site of A ), then S v, t ( ) also has range 1 . However, say in one dimension, it typically couples a site j to both j -1 and j +1 simultaneously. This makes experimental implementation of the swindle models in some quantum simulation platforms more complicated. Further, the swindle model is translationally invariant, and

so does not enjoy the robustness of disordered systems to perturbations.

We address both of these issues by defining a modified swindle map ˜ ( S v, t ) . We split the loop defining the swindle unitary into three segments, t ∈ [0 , 1 3 ) , t ∈ [ 1 3 , 2 3 ) , and t ∈ [ 2 3 , 1] . Within the first segment, sites coupled along even numbered bonds along the first dimension of the lattice are coupled, while odd numbered bonds are coupled in the second segment. The third segment introduces on-site disorder to the model.

We pick arbitrary smooth functions η 1 2 3 , , ( ) t which respectively map t ∈ [0 , 1 3 ] , [ 1 3 , 2 3 ] , and [ 2 3 , 1] to η ∈ [0 , 1] . Sensible choices for these functions are monotonic increasing. We also choose disordered ( d -1) -dimensional evolution operators u dis j ( η ) , which are required to commute with any desired symmetry of the model. Then, on open boundary conditions with L ∈ 2 Z sites along the first dimension, we have

$$\tilde { S } ( v, t ) = \prod _ { j = 0 } ^ { L / 2 - 1 } u _ { 2 j, 2 j + 1 } ( v, \eta _ { 1 } ( t ) ) \quad \text{(E1)}$$

[where u j,j +1 is defined in Eq. (7)] for t ∈ [0 , 1 3 ) ,

$$\tilde { S } ( v, t ) = v _ { 0 } \left ( \prod _ { j = 1 } ^ { L / 2 - 1 } u _ { 2 j - 1, 2 j } ( v, \eta _ { 2 } ( t ) ) v _ { 2 j - 1 } ^ { \dagger } v _ { 2 j } \right ) v _ { L } ^ { \dagger } \text{ (E2)} \quad \text{FIG. } \frac { \text{fig. } \frac { \text{f} } { \text{where} } }$$

for t ∈ [ 1 3 , 2 3 ) , and

$$\tilde { S } ( v, t ) = u _ { j } ^ { \text{dis} } ( \eta _ { 3 } ( t ) ) v _ { 0 } v _ { L } ^ { \dagger } \text{ \quad \ \ } ( \text{E3} ) \text{ \ \ } \text{das}$$

for t ∈ [ 2 3 , 1] (Fig. 3).

The Hamiltonian associated to the unitary ˜ S can be obtained by differentiation,

Away from the boundaries, ˜ S acts as the identity at t = 2 3 . The LIOMs are then determined by u dis j . At the edges, ˜ S does not act as the identity, and instead acts as v or v † .

$$H ( t ) = [ i \partial _ { t } \tilde { S } ( v, t ) ] \tilde { S } ( v, t ) ^ { \dagger }. \quad \quad ( \mathbb { E } 4 ) \quad \text{syn}$$

## 2. QP pump

We carry through the modified swindle construction explicitly for the QP pump. This phase is of particular interest, as its edge states can be used to prepare highly excited non-classical states in quantum cavities [133]. Thus, an experimental realization of this model is especially desirable.

Our model is based on a one-dimensional chain of qubits with open boundary conditions and length L ∈ 2 Z . The construction for odd numbers of qubits is essentially the same. The resulting model will be qualitatively similar to that in Refs. [32, 134].

FIG. 5. (a) The power pumped into the θ drive in the state ρ 0 -, ( L -1)+ (E20). The average value is quantized to ωW , where W is the winding number of the unitary used to construct the QP pump model and ω is the frequency of the θ drive. The power is only nonzero when t mod 1 ∈ [0 , 1 3 ) . (b) The quantized power appears as an average slope (red dashed) in the energy absorbed by the θ drive (blue). In fact, the modified swindle model has E ± ( ) t = ∓ ωWt at integer times. Parameters: ω = (1+ √ 5) / 2 , W = 1 , η 1 = sin (3 2 πt/ 2) .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00825v2-figure-cac841d6-f49a-4554-a3c6-01f85fdff50b.jpg)

The model has a U (1) symmetry generated by

$$Q = \sum _ { j } Q _ { j } = \sum _ { j } \frac { 1 } { 2 } ( 1 + \sigma _ { j } ^ { z } ), \quad \text{(E5)} \\ \text{choose a path to the swap gate which preserves this}$$

so we choose a path to the swap gate which preserves this symmetry. Indeed, we may take

$$\text{swap} _ { j, j + 1 } ( t ) = \exp \left [ - i \frac { 2 \pi } { 8 } t \left ( \vec { \sigma } _ { j } \cdot \vec { \sigma } _ { j + 1 } \right ) \right ] \quad \text{(E6)}$$

up to a global phase, where ⃗ σ j · ⃗ σ j +1 = σ σ x j x j +1 + σ σ y j y j +1 + σ σ z j z j +1 .

The lower dimensional QCA to which the swindle map will be applied will actually be a loop of QCA in zero dimensions. That is, a loop of unitary matrices,

$$v _ { j } ( \theta ) = \exp \left [ - i \frac { W } { 2 } \theta \sigma _ { j } ^ { z } \right ], \text{ \quad \ \ } \text{(E7)}$$

where W ∈ Z is the winding number and θ ∈ [0 , 2 π ) parameterizes the loop. Note that all W ∈ Z give distinct winding numbers with U (1) symmetry. Without

this symmetry, all odd (even) winding numbers would be connected by homotopies (Appendix C 2).

̸

As the model is quasiperiodically driven, we will take θ t ( ) = ωt mod 2 π to be a linear function of t with ω ∈ Q irrational.

We take the disorder in the model to be on-site and diagonal in the σ z basis,

$$u _ { j } ^ { \text{dis} } ( t ) = \exp \left [ - i t \delta _ { j } \sigma _ { j } ^ { z } \right ], \quad \quad ( \text{E8} ) \quad \text{per} \\ \text{ter}$$

where δ j is, say, uniformly random in [ -∆ ∆] , .

assists in stabilizing localization [28, 99].

The Hamiltonian is then defined piecewise as

$$\begin{array} { c } \text{II take} \\ \omega \not \in \mathbb { Q } \\ \text{e and} \end{array} \begin{cases} \begin{cases} \omega h ^ { 1 } ( \eta _ { 1 } ( t ), \theta ) + \dot { \eta } _ { 1 } ( t ) h ^ { 2 } ( \eta _ { 1 } ( t ), \theta ), & t \in [ 0, \frac { 1 } { 3 } ), \\ \omega h ^ { 3 } ( \eta _ { 2 } ( t ), \theta ) + \dot { \eta } _ { 2 } ( t ) h ^ { 4 } ( \eta _ { 2 } ( t ), \theta ), & t \in [ \frac { 1 } { 3 }, \frac { 2 } { 3 } ), \\ \omega h ^ { 5 } + \dot { \eta } _ { 3 } ( t ) h ^ { 6 }, & t \in [ \frac { 2 } { 3 }, 1 ], \\ \end{cases} \\ \text{where} \, \dot { n } _ { n } = \partial n _ { n } \text{ and } H ( t, \theta ) \text{ has neriod } T = 1 \text{ in } t \text{ and} \end{cases}$$

where ˙ η k = ∂ η t k and H t, θ ( ) has period T = 1 in t and period 2 π in θ . Defining σ ± j = σ x j ± iσ y j , the individual terms for the first pulse are

The final elements of the modified swindle model we need are the functions η 1 2 3 , , ( ) t . These are essentially the pulse shapes which implement each unitary as a function of t , and their specific choice should be tailored to the experimental platform in mind. We leave them unspecified, beyond requiring that they are smooth, monotonic, and their derivatives vanish at their end points. This ensures that there are no discontinuities in the pulses, which also

$$h ^ { 1 } ( \eta _ { 1 }, \theta ) = \sum _ { j = 0 } ^ { L / 2 - 1 } h _ { 2 j, 2 j + 1 } ^ { 1 } ( \eta _ { 1 }, \theta ) \quad \text{(E10)}$$

$$h ^ { 2 } ( \eta _ { 1 }, \theta ) = \sum _ { j = 0 } ^ { L / 2 - 1 } h _ { 2 j, 2 j + 1 } ^ { 2 } ( \eta _ { 1 }, \theta ), \quad \text{(E11)}$$

where

$$h _ { j, j + 1 } ^ { 1 } ( \eta _ { 1 }, \theta ) = \frac { W } { 4 } \left [ ( 1 - \cos ( \pi \eta _ { 1 } ) ) ( \sigma _ { j } ^ { z } - \sigma _ { j + 1 } ^ { z } ) - \frac { i } { 2 } \sin ( \pi \eta _ { 1 } ) e ^ { - i W \theta } \sigma _ { j } ^ { + } \sigma _ { j + 1 } ^ { - } + \text{H.c.} \right ], \text{ and}$$

$$h _ { j, j + 1 } ^ { 2 } ( \eta _ { 1 }, \theta ) = \frac { 2 \pi } { 8 } \stackrel { \bullet } { \left [ \sin ( W \theta ) \sin ( \pi \eta _ { 1 } ) ( \sigma _ { 2 j } ^ { z } - \sigma _ { 2 j + 1 } ^ { z } ) + \frac { e ^ { - i W \theta } } { 2 } ( 1 - \cos ( W \theta ) - i \sin ( W \theta ) \cos ( \pi \eta _ { 1 } ) ) \sigma _ { j } ^ { + } \sigma _ { j + 1 } ^ { - } + \text{H.c.} \right ] }.$$

$$\cdot _ { 1 } ) & - \frac { \imath } { 2 } \sin ( \pi \eta _ { 1 } ) e ^ { - i W \theta } \sigma _ { j } ^ { + } \sigma _ { j + 1 } ^ { - } + \text{H.c.} \Big | \,, \text{ and } \\ \sigma _ { 2 i + 1 } ^ { z } & + \frac { e ^ { - i W \theta } } { \partial } ( 1 - \cos ( W \theta ) - i \sin ( W \theta ) \cos ( \pi \eta _ { 1 } ) ) \sigma _ { i } ^ { + } \sigma _ { i + 1 } ^ { - } + \text{H.c.} \Big | \,.$$

For the second pulse, we have

$$h ^ { 3 } ( \eta _ { 2 }, \theta ) = \frac { W } { 2 } ( \sigma _ { 0 } ^ { z } - \sigma _ { L - 1 } ^ { z } ) + \sum _ { j = 1 } ^ { L / 2 - 1 } h _ { 2 j - 1, 2 j } ^ { 3 } ( \eta _ { 2 }, \theta ), \quad \text{and} \quad h ^ { 4 } ( \eta _ { 2 }, \theta ) = \sum _ { j = 1 } ^ { L / 2 - 1 } h _ { 2 j - 1, 2 j } ^ { 2 } ( \eta _ { 2 }, \theta ),$$

where

$$h _ { j, j + 1 } ^ { 3 } ( \eta _ { 2 }, \theta ) & = \frac { W } { 4 } \left [ [ 1 + \cos ( W \theta ) + \cos ( 2 \pi \eta _ { 2 } ) ( 1 - \cos ( W \theta ) ) ] ( \sigma _ { j } ^ { z } - \sigma _ { j + 1 } ^ { z } ) \\ & \quad - \frac { i } { 2 } e ^ { - i W \theta } [ \sin ( \pi \eta _ { 2 } ) ( 1 - 2 i \sin ( W \theta ) ) + \sin ( 2 \pi \eta _ { 2 } ) ( 1 - \cos ( W \theta ) ) ] \sigma _ { j } ^ { + } \sigma _ { j + 1 } ^ { - } + H. c. \right ].$$

Finally, the disorder pulse has at the cost of introducing a time-independent term

$$h ^ { 5 } = \frac { W } { 2 } ( \sigma _ { 0 } ^ { z } - \sigma _ { L - 1 } ^ { z } ), \text{ \ and \ } h ^ { 6 } = \sum _ { j = 0 } ^ { L - 1 } \delta _ { j } \sigma _ { j } ^ { z }.$$

Additionally, rotating frame transformations generated by combinations of σ z j can be used to eliminate some of the complex time dependence of the coupling terms. For instance, making the transformation

$$\tilde { S } ( t ) \mapsto e ^ { i \omega t \frac { W } { 2 } \sum _ { j } j \sigma _ { j } ^ { z } } \tilde { S } ( t ) \text{ \quad \ \ } ( \text{E17} ) \text{ \ \ \ } \text{one}$$

↦

eliminates the e -iWθ factors in all the expressions above,

$$H ^ { \prime } = - \omega \frac { W } { 2 } \sum _ { j = 0 } ^ { L - 1 } j \sigma _ { j } ^ { z }$$

to the Hamiltonian. That is, the time-dependent vector potential in the model can be traded for a linear scalar potential. On the other hand, if implementing time-dependent detunings σ z j -σ z j +1 is more challenging in some architecture, moving to a frame co-rotating with those terms of the Hamiltonian moves this time dependence to the σ + j σ -j +1 exchange term.

The main interesting feature of the QP pump is its edge states, which transfer energy from one drive to the

other. To compute the pumped power in the edge state, we evaluate the expectation value of

$$P ( t, \theta ) = - \omega \partial _ { \theta } H ( t, \theta ), \quad \quad ( \text{E19} ) \quad \text{$\rho$.} _ { \text{by}$}$$

which is the work done by the qubit chain on the drive with phase θ . We evaluate this operator in a mixture of quasienergy states, with the initial density matrix being

$$\rho _ { 0 \pm, ( L - 1 ) \pm } = \frac { 1 } { 2 ^ { L } } ( 1 \pm \sigma _ { 0 } ^ { z } ) ( 1 \pm \sigma _ { L - 1 } ^ { z } ). \quad ( E 2 0 ) & \quad \langle. \\ \text{This is maximally mixed in the bulk LION states, but} \\ \text{has a definite state at both edges. The expectation value} & \quad \text{Firt}$$

This is maximally mixed in the bulk LIOM states, but has a definite state at both edges. The expectation value of the power operator is thus

$$\langle P ( t, \theta ) \rangle _ { \pm \pm } = \text{Tr} \left [ P \tilde { S } \rho _ { 0 \pm, ( L - 1 ) \pm } \tilde { S } ^ { \dagger } \right ] \quad \text{(E21)} \quad \text{(E)}$$

The model is constructed so that τ j ( t = 0) = σ z j , and the Heisenberg picture operator σ z H j ( ) t is a function of ⃗ θ . We can thus find the edge LIOMs within a period by 'reverse Heisenberg evolution' [28] of the σ z 0 operator. We have, for t ∈ [0 , 1 3 )

$$\inf _ { \mathbb { G } } \tau _ { 0 } ^ { z } ( t, \theta ) = u _ { 0, 1 } ( v ( \theta ), \eta _ { 1 } ( t ) ) \sigma _ { 0 } ^ { z } u _ { 0, 1 } ^ { \dagger } ( v ( \theta ), \eta _ { 1 } ( t ) ) \quad ( \mathbb { E } 2 2 )$$

and τ z 0 = σ z 0 otherwise. Equation (E21) then gives that

$$\langle P ( t, \theta ) \rangle _ { \pm \pm } = \frac { 1 } { 2 ^ { L } } \text{Tr} \left [ P ( 1 \pm \tau _ { 0 } ^ { z } ) ( 1 \pm \tau _ { L - 1 } ^ { z } ) \right ]. \quad ( \text{E23} ) \\ \intertext { a rther using that the power is a traceless local operator, } P = \sum _ { \lambda } P _ { i, i + 1 }, \text{ and that } \tau _ { 0 } ^ { z } _ { L - 1 } \text{ have disjoint support, we }$$

Further using that the power is a traceless local operator, P = ∑ j P j,j +1 , and that τ z 0 ,L -1 have disjoint support, we can expand Eq. (E23) as where we suppressed the dependence on t and θ in the right hand side.

$$\int _ { \text{the} } ^ { 2 1 } \left \langle P ( t, \theta ) \right \rangle _ { \pm \pm } = \pm \frac { 1 } { 2 ^ { L } } \text{Tr} \left [ P \tau _ { 0 } ^ { z } \right ] \pm \frac { 1 } { 2 ^ { L } } \text{Tr} \left [ P \tau _ { L - 1 } ^ { z } \right ].$$

Calculating the LIOMs from the definition of u 0 1 , gives

$$u _ { 0, 1 } ( v ( \theta ), \eta _ { 1 } ) \sigma _ { 0 } ^ { z } u _ { 0, 1 } ^ { \dagger } ( v ( \theta ), \eta _ { 1 } ) & = \frac { 1 } { 2 } \left [ \sigma _ { 0 } ^ { z } + \sigma _ { 1 } ^ { z } + \frac { 1 } { 2 } [ 1 + \cos ( W \theta ) + \cos ( 2 \pi _ { 1 } ) ( 1 - \cos ( W \theta ) ) ] ( \sigma _ { 0 } ^ { z } - \sigma _ { 1 } ^ { z } ) \\ & \quad + \frac { e ^ { - i W \theta } } { 4 } [ 2 \sin ( \pi \eta _ { 1 } ) \sin ( W \theta _ { 1 } ) + i \sin ( 2 \pi \eta _ { 1 } ) ( 1 - \cos ( W \theta ) ) ] \sigma _ { 0 } ^ { + } \sigma _ { 1 } ^ { - } + H. c. \right ].$$

Similarly, τ z L -1 ( t, θ ) is given by σ z L -1 except during the first pulse, which is determined by

$$u _ { L - 2, L - 1 } ( v ( \theta ), \eta _ { 1 } ) \sigma _ { L - 1 } ^ { z } u _ { L - 2, L - 1 } ^ { \dagger } ( v ( \theta ), \eta _ { 1 } ) & = \frac { 1 } { 2 } \left [ \sigma _ { L - 2 } ^ { z } + \sigma _ { L - 1 } ^ { z } - \frac { 1 } { 2 } [ 1 + \cos ( W \theta ) + \cos ( 2 \pi \eta _ { 1 } ) ( 1 - \cos ( W \theta ) ) ] ( \sigma _ { L - 2 } ^ { z } - \sigma _ { L - 1 } ^ { z } ) \\ & - \frac { e ^ { - i W \theta } } { 4 } [ 2 \sin ( \pi \eta _ { 1 } ) \sin ( W \theta _ { 1 } ) + i \sin ( 2 \pi \eta _ { 1 } ) ( 1 - \cos ( W \theta ) ) ] \sigma _ { 0 } ^ { + } \sigma _ { 1 } ^ { - } + H. c. ]. \quad ( E 2 6 )$$

The power pumped by the left (smaller j ) edge state is

$$\frac { 1 } { 2 L } \, \text{Tr} \left [ P \tau _ { 0 } ^ { \tau } \right ] & = - \frac { \omega } { 2 L } \, \text{Tr} \left [ ( \omega \partial _ { \theta } h _ { 0, 1 } ^ { 1 } + \dot { \eta } _ { 0 } \partial _ { \theta } h _ { 0, 1 } ^ { 2 } ) \tau _ { 0 } ^ { \tau } \right ] = \frac { \omega W } { 4 } \left [ \frac { \omega W } { 2 } \sin ( W \theta ) ( 1 - \cos ( 2 \pi \eta _ { 1 } ) ) \\ & \quad + \pi \dot { \eta } _ { 1 } ( \sin ( 2 \pi \eta _ { 1 } ) ( 1 - \cos ( W \theta ) ) - \sin ( \pi \eta _ { 1 } ) ) \right ],$$

$$( 1 - \cos ( W \theta ) ) - \text{si}$$

for t ∈ [0 , 1 3 ) , and Tr[ Pτ z 0 ] = 0 otherwise. Further, Tr[ Pτ z 0 ] = -Tr [ Pτ z L -1 ] , so the pumped power is nonzero only when the two edge states are opposite. We will subsequently let ± refer only to the occupation of the LIOM at site j = 0 , and assume the other edge LIOM is in the opposite state. Then the total pumped power is

$$\langle P ( t, \theta ) \rangle _ { \pm } = \pm \frac { \omega W } { 2 } \left [ \frac { \omega W } { 2 } \sin ( W \theta ) \left ( 1 - \cos ( 2 \pi \eta _ { 1 } ( t ) ) \right ) + \pi \dot { \eta } _ { 1 } ( t ) \left ( \sin ( 2 \pi \eta _ { 1 } ( t ) ) \right ) ( 1 - \cos ( W \theta ) ) - \sin ( \pi \eta _ { 1 } ( t ) ) \right ) \right ] \quad ( E 2 8 )$$

for t ∈ [0 , 1 3 ) . This function is plotted as a function of t in Fig. 5.

grating ⟨ P t, θ ( ) ⟩ ± over t and θ . This gives

We can compute the average pumped power by inte-

$$\nearrow \quad \bar { P } _ { \pm } = \frac { 1 } { 2 \pi } \int _ { 0 } ^ { 2 \pi } \mathrm d \theta \int _ { 0 } ^ { 1 } \mathrm d t \left \langle P ( t, \theta ) \right \rangle _ { \pm } = \mp \omega W. \quad ( \mathrm E 2 9 )$$

That is, W photons of frequency ω are pumped into (or out of, if W &lt; 0 ) the θ drive per period of t (which we have set to 1) when ⟨ τ z 0 ⟩ = -1 and ⟨ τ z L -1 ⟩ = +1 .

## a. Invariant

Finally, we briefly explain where this instance of the QP pump sits in the classification of Appendix D3. Specifically, we find its associated invariant.

The QP pump is U (1) -symmetric and has d = 1 , n = 2 . It thus has the same classifying group as loops of unitaries in d = 0 with U (1) symmetry. Indeed, this is the object from which our model was constructed. From Appendix D3, to get the classification of v θ ( ) we must find its one-dimensional winding number in each symmetry sector. We have

$$v ( \theta ) = e ^ { - i W \theta / 2 } \begin{pmatrix} 1 & 0 \\ 0 & e ^ { i W \theta } \end{pmatrix}, \ \ Q _ { 0 } = \begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix}. \ \ ( E 3 0 ) \quad \ m u l a \omega \\ \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

Thus, up to a global phase, v θ ( ) has winding number W in the charge 0 sector, and winding number 0 in the charge 1 sector.

The invariant associated to v θ ( ) is a coset [Eq. (D49)] of formal ratios of representations with the numerator weighted by the winding numbers. We identify the charge 0 representation with 1 and the charge 1 representation with x ∈ N [ x, x -1 ] ∼ = Rep( U (1)) . Then the invariant is

$$\tilde { w } _ { 1 } ^ { Q } [ v ] = \frac { 1 \cdot W + x \cdot 0 } { 1 + x } + \mathbb { Z } = \frac { W } { 1 + x } + \mathbb { Z }. \quad ( \mathbb { E } 3 1 ) \quad \text{emen}$$

We see that all distinct values of W indeed define inequivalent values of the invariant.

Suppose now that we explicitly broke the U (1) symmetry. Then the charge sector 1 ↔ x should be mapped to having charge 0 . The new invariant can be found by evaluating the rational function ˜ w Q 1 [ v ] at x = 1 (Appendix D 5). We find

$$\tilde { w } _ { 1 } ^ { Q } [ v ] \mapsto \frac { W } { 2 } + \mathbb { Z } \in \mathbb { Q } / \mathbb { Z }.$$

↦

Now, all even W are equivalent to the trivial phase. Interestingly, the odd W phases remain non-trivial. They no longer support a quantized edge current, as we can deform the model with W = +1 to W = -1 , but there is apparently still some nontrivial dynamics at the edge.

̸

Indeed, we believe the odd W phase without symmetry to be the two-tone-driven phase of Refs. [29, 30]. While the solvable limit of the model in Ref. [29] is highly singular, which makes computing invariants from integral formulae inconvenient, the edge phenomenology described in Ref. [29] is consistent with what we expect from having ˜ w 1 = 0 mod 1 . Namely, the edge dynamics was argued to be anomalous in that it could not be achieved by a zerodimensional localized [in the sense that a continuous micromotion operator exists as in Eq. (A24)] system which did not occur at the edge of an extended one-dimensional system. As the model in Ref. [29] is based around a version of the AKLT model, which can be trivialized when stacked with itself, it is thus reasonable to suppose that ˜ w 1 = 1 2 mod 1 / for that model, as this is the only element of Q Z / of order 2. Thus, the abstract classification reveals a relationship between the interacting QP pump and Refs. [29, 30] which was not clear previously. It also extends the arguments for stability in Ref. [29] to the nonperturbative regime, where a picture of emergent symmetry (as employed in that reference) becomes inapplicable.
David M. Long Dominic V. Else
Topological Phases of Many-Body Localized Systems: Beyond Eigenstate Order
Many-body localization (MBL) lends remarkable robustness to nonequilibrium
phases of matter. Such phases can show topological and symmetry breaking order
in their ground and excited states, but they may also belong to an anomalous
localized topological phase (ALT phase). All eigenstates in an ALT phase are
trivial, in that they can be deformed to product states, but the entire
Hamiltonian cannot be deformed to a trivial localized model without going
through a delocalization transition. Using a correspondence between MBL phases
with short-ranged entanglement and locality preserving unitaries - called
quantum cellular automata (QCA) - we reduce the classification of ALT phases to
that of QCA. This method extends to periodically (Floquet) and
quasiperiodically driven ALT phases, and captures anomalous Floquet phases
within the same framework as static phases. The QCA framework further
generalizes to include symmetry-enriched ALT phases (SALT phases), and provides
a large class of soluble models suitable for realization in quantum simulators.
In systematizing the study of ALT phases, we both greatly extend the
classification of interacting nonequilibrium systems and clarify a confusion in
the literature which implicitly equates nontrivial Hamiltonians with nontrivial
ground states.