Supersymmetric Many-Body Systems from Partial Symmetries: Integrability, Localization and Scrambling

Partial symmetries are described by generalized group structures known as symmetric inverse semigroups. We use the algebras arising from these structures to realize supersymmetry in (0+1) dimensions and to build many-body quantum systems on a chain. This construction consists in associating appropriate supercharges to chain sites, in analogy to what is done in spin chains. For simple enough choices of supercharges, we show that the resulting states have a finite non-zero Witten index, which is invariant under perturbations, therefore defining supersymmetric phases of matter protected by the index. The Hamiltonians we obtain are integrable and display a spectrum containing both product and entangled states. By introducing disorder and studying the out-of-time-ordered correlators (OTOC), we find that these systems are in the many-body localized phase and do not thermalize. Finally, we reformulate a theorem relating the growth of the second Renyi entropy to the OTOC on a thermal state in terms of partial symmetries.


Introduction
In quantum mechanics, a symmetry is implemented by requiring a system to be invariant under a certain set of transformations. These transformations must form a group, whose action on the associated Hilbert space is realized by unitary or anti-unitary operators, as stated by Wigner's theorem [1]. If the symmetry is obeyed only by part of the system, one can still implement it by the use of sets of transformations, but these have to be taken to form an inverse semigroup rather than a group. This generalization leads to the notion of partial symmetries.
An inverse semigroup (S, * ) is a pair formed by a set S and an associative binary operation * , such that every element has a unique inverse. So, for a given x ∈ S, there exists a unique y ∈ S such that y * x * y = y and x * y * x = x. (1.1) There is no single identity on an inverse semigroup, but only idempotents or projectors which can be thought of as partial identities. The elements of S are called partial symmetries.
Since they do not act on the entire Hilbert space, but only on parts of it, there is no unitary representation of these operators. Such structures do arise in quantum mechanics but have often been discarded without getting much consideration. It turns out, however, that inverse semigroups are relevant for physics as they provide the precise description for invariances that underlie certain physical systems [2]. To convince the readers, we illustrate this point in two instances of physical interest.
The first instance concerns a complete classification of tilings of R n . It is known that, for every such tiling, there exists an inverse semigroup associated with it [3]. Whereas periodic crystal structures -familiar to n-dimensional crystallography -are well described by group theory, aperiodic structures like quasicrystals [4][5][6] -see Fig. 1 -are associated to aperiodic tilings described by inverse semigroups [3,[8][9][10]. Hamiltonian [11], works on gap-labelling theorems for such systems and, more generally, for substitution sequence Hamiltonians [12,13].
The second instance concerns operators in Hilbert space. In a finite-dimensional vector space V, a complex matrix A is always decomposable to polar factorization, A = U R or A = LU , where L, R are non-negative Hermitian matrices and U is a unitary matrix. In an In this paper, we shall work with a particular type of inverse semigroup, known as the symmetric inverse semigroup (SIS), which can be thought of as the analog of the permutation group in this context. Just like any group can be embedded into a symmetric group, according to the Wagner-Preston representation theorem [14], any inverse semigroup can be embedded into a SIS such that it is isomorphic to some sub SIS [2]. This means that we can think of the SISs as building blocks of an arbitrary inverse semigroup. Thus, by working with SISs, there is no loss of generality. This result is the counterpart of Cayley's representation theorem in group theory stating that any group is isomorphic to the group of permutations or some subgroup of it.
When exploring a physical system, it is always a good idea to endow it with extra symmetries that might constrain the dynamics and bring the evaluation of observables under computation control. A symmetry that is proven notoriously capable of doing so is supersymmetry. This is a relation between bosons and fermions, which was first introduced in the context of relativistic quantum field theory as an extension of the Poincaré group, see for example [15][16][17][18][19] and [20] for a review. In this paper, we shall focus on (0 + 1)-dimensional supersymmetric systems, viz. supersymmetric quantum mechanics [21,22]. This system has proven to be a very useful arena, where problems with wide range of potentials could be solved exactly [23].
A supersymmetric system is equipped with a multiplet of fermionic charges, called supercharges, Q, that generate the supersymmetry transformations. By construction, the vacuum expectation value of top auxiliary components of supermultiplets is the order parameter of spontaneous supersymmetry breaking. After eliminating these auxiliary components, the system's Hamiltonian is given by H = {Q † , Q}, implying that the vacuum expectation value of the Hamiltonian is positive semidefinite. If the supersymmetry is unbroken, Q|0 is zero and 0|H|0 = 0. If supersymmetry is spontaneously broken, Q|0 is non-zero and 0|H|0 > 0.
Thus, the ground state energy provides an easily calculable order parameter of spontaneous supersymmetry breaking.
Another related order parameter is the Witten index [24], ∆ = N fermion , the number of bosonic zero-energy ground states minus the number of fermionic zero-enegy ground states. If this quantity is non-zero, then there ought to be unequal numbers of bosonic and fermionic zero-energy ground states. If the quantity is zero, then there are either equal pairs of bosonic and fermionic zero-energy ground states or the ground states have positive energy.
In the former case, the system preserves supersymmetry. In the latter case, the system breaks supersymmetry. So, provided the system's energy spectrum is discrete, the Witten index is a quantized quantity that cannot be continuously changed by supersymmetry preserving deformations.
Topological phases of matter are described by topological invariants, examples of which include the quantum double models of Kitaev described by the genus of a surface [25], or the symmetry protected states of matter described by the cohomology of corresponding symmetry groups [26]. These invariants are often computable as the partition functions of systems or as the trace of corresponding transfer matrices [27]. Along these lines, we propose to use the Witten index, which can be thought of as a twisted partition function of systems or as a supertrace 2 of the corresponding transfer matrix formed out of supersymmetric Hamiltonian, to also be a useful indicator for interesting phases of many-body systems, protected by the global supersymmetry. 3 Following these lines of reasonings, we are naturally led to construct supersymmetric systems with partial symmetries, or supersymmetric SISs. We will do just that in the following, focusing in particular on supersymmetric many-body systems on a lattice, where the internal degrees of freedom are built from SIS algebra elements. In this way, we implement a supersymmetric algebra on the Hilbert space spanned by these elements, which can be thought of as realizing supersymmetric algebras out of partial symmetries. Many-body supersymmetry preserving zero-energy ground states can be constructed as eigenstates of many-body Hamiltonians by considering supersymmetry in the non-relativistic setting, as in statistical physics, see [33,34] and citations therein. Supersymmetric many-body systems have been considered in the past [35], where their lattice formulation was done by realizing the global supercharges on the total Hilbert space of the system [36][37][38][39][40][41]. Exact solutions for these models were obtained in [42]. Entanglement entropy in such systems have also been considered [43].
An important aspect of the models we construct with partial symmetries is that they are generically quantum integrable: all of these models have as many local integrals of motion as the number of sites, as we shall see. This is the case when we construct our supersymmetric system with a unique grading of the Hilbert space spanned by the elements of the SIS algebra. 4 Given this, we ask the question about scrambling properties in such supersymmetric systems, which forms the subject of the remaining part of the main text.
Recently, the topic of scrambling has proven to be of great interest in many fields, from the physics of black holes to quantum information. Concretely, scrambling can be used to identify the onset for quantum chaos [44] and the propagation of entanglement [45]. From a gravitational/holographic perspective, black holes are conjectured to be the fastest scramblers in Nature [46] and, in the same spirit, it has been conjectured that chaos cannot grow faster than in Einstein gravity [47,48]. These works provide several hints connecting fast scrambling and thermalized phases, which are embedded into a large program concerned with the full 2 We use the twisted sum or the supertrace as the Hilbert space is graded. 3 The Witten index also has a deep connection with the topology of the bundle spaces upon which the Hamiltonian acts. This is the content of the Atiyah-Pataudi-Singer index theorem [28][29][30][31]. Generalizations of this index have also been considered [32]. 4 However, we can work in a scenario where the system is constructed out of two different gradings of the Hilbert space and in such a case we obtain non-integrable systems. This is discussed in App. B.
understanding of thermalization in complex quantum systems.
On the opposite extreme, and the focus of this work, are systems where time evolution results in a localized phase and no thermalization occurs at all. Roughly speaking, thermalization is the process by which a system evolves to a point in which it can be well described by few thermodynamical quantities. This seems to conflict with the quantum mechanical intuition that evolution is unitary and the system is not expected to lose information about the initial state. The resolution to this puzzle is the eigenstate thermalization hypothesis (ETH) [49], which says that a quantum system thermalizes if it does so irrespective of the initial state it was prepared in, thereby also giving a proper definition of quantum thermalization. This brings up the possibility of the existence of systems that have states where such a process of equilibration does not occur and, consequently, the ETH is not satisfied. A classic example of such a system includes integrable models whose local integrals of motion prevent the transport of conserved quantities hindering the reach of equilibrium. The presence of disorder is another source of localization and loss of thermalization. In situations in which the phenomenon is described by one-particle effects, this is called Anderson localization (AL) [50], while in cases in which it is described by the many-particle effects, it is known as many-body localization (MBL). A recent review on thermalization and localization is provided by [51].
The supersymmetric models we present here fall into the MBL class. The presence of "local" integrals of motion in these systems prevent the transport of the conserved quantities and thus we expect to find MBL states in these systems. Since any such supersymmetric system constructed out of a unique grading of the SIS algebra is expected to have this property, we hereby developed a method of constructing supersymmetric MBL states using the SIS algebras.
The system we study is a long-ranged interacting one and we emphasize that this is done just for simplicity. In App. A, we construct short-ranged interacting supersymmetric systems built out of a unique grading of the SIS algebra, which continue to have local integrals of motion and are thus expected to possess MBL states.
We shall deduce the MBL property of these systems by studying the out-of-time-order correlators (OTOC) [44,45,47]. OTOCs have been recently used for analyzing localized systems, in particular in distinguishing MBL and AL phases in [52][53][54][55][56]. This leads to connections between the OTOC's, scrambling, and localization. We briefly summarize the procedure used for this purpose.
To probe scrambling using OTOCs, one has to consider two local, arbitrary operators at different times, say V (0) ≡ V and W (t) = e iHt W (0)e −iHt , and the squared commutator between these probe operators in some state of interest, which can be taken to be a thermal state with temperature β −1 2) or, alternatively, the out-of-time order correlator that is contained in Eq. (1.2). The behavior of these quantities as a function of time has crucial information about the system. For example, thermalized, MBL and AL phases can be diagnosed, respectively, by an exponential decay, a power-law decay and a constant behavior of the OTOC. We will refer to the systems presenting less than exponential OTOC as slow scramblers. It should be noted that OTOC have been considered long ago in the context of semi-classical methods in superconductivity [57].
With this setup we organize the paper as follows. In Sec. 2, we review the basics of partial symmetries and SIS and construct single-particle supersymmetric systems out of SISs.
In Sec. 3, we construct supersymmetric models in (0 + 1) dimensions based on SIS with supercharges leading both to free and to interacting Hamiltonians. In Sec. 4, we apply this setup to investigate scrambling in interacting many-body systems. We will present a toy model with quenched disorder that exhibits a supersymmetric MBL phase and argue that systems generated using supersymmetric realizations of SIS should be slow scramblers. An outlook for future work and a discussion of the scope of these supersymmetric SIS models is given in Sec. 5.
We also include several appendices, where we discuss other possible supersymmetric systems based on SISs. In App. A, we show other supercharges that exhibit a spectrum which is more complicated than the ones considered in the main body of the paper. App. B sets the ground for possible ways of constructing non-integrable supersymmetric many-body systems with partial symmetries. We construct supersymmetric systems that have entangled eigenstates in App. C. Finally, in App. D we sketch how to obtain para-supersymmetric quantum mechanical systems from partial symmetries.
2 Partial symmetries and supersymmetric quantum mechanics In this section, we begin with the simplest quantum system with partial symmetries: a oneparticle or one-site quantum mechanical system defined on a finite-dimensional Hilbert space.
Let S n = {1, 2 · · · , n}. Consider the set of all partial bijections on S n together with the usual composition rule, which is binary and associative. This pair forms a SIS, denoted by S n = (S n , * ). We can also form a class of SIS by choosing subsets of order p ≤ n, S p , and considering the set of partial bijections on this subset. We will refer to the resulting SIS as S n p . Taking the elements of S n p as basis, we will construct a Hilbert space, whose inner product is defined by the natural pairing of the basis. In what follows, we will work with such SISs and Hilbert spaces.

Diagrammatics for SISs
We shall first introduce a convenient diagrammatic way of representing the elements that renders the binary operations transparent and defines an algebra over S n p .
To illustrate this in a transparent manner, let us start with the simplest example, S 2 1 , whose diagrammatics are shown in Fig. 2. The partial symmetry elements of S 2 1 are denoted by x i,j , with i, j ∈ {1, 2}, and obey the following composition rule (2.1) The indices i and j can be thought of, respectively, as the domain and range of the partial symmetry operation. The product between these elements is null when the range of the first element is different from the domain of the second element it is being composed with. Note that this product is non-commutative.
The next step up in complexity is illustrated by S 3 1 , which is made up of nine elements, x i,j with i, j ∈ {1, 2, 3}. These partial symmetries of S 3 1 are depicted in Fig. 3. Further moving up, we can construct a SIS with arbitrary n and p, S n p . For the sake of definiteness, we will focus on S 3 1 , but the formulation is straightforwardly generalizable to any (p, n).
We next associate Hilbert spaces to every SISs. The Hilbert spaces we consider are spanned by the partial symmetries of the chosen SIS, meaning that we are working with the algebra of Composition on S 2 1 Figure 2: Diagrammatic representation of S 2 1 . The composition rules are obtained by tracing arrows connecting two elements. If one cannot trace a continuous arrow, the product is 0. this SIS. Therefore, the SIS acquires a vector space structure and an inner product is naturally defined by For instance, an arbitrary element of the Hilbert space spanned by the elements of S 2 1 is given This is equivalent to working in the regular representation of the chosen SIS.

Supersymmetric systems from SISs
We now turn to the construction of a single-site, one-particle supersymmetric system. The first step is to realize the supersymmetry algebra from the SISs. Start with S 2 1 and define q = x 1,2 and q † = x 2,1 , As q and q † are nilpotent, they can be thought of as supercharges. Out of these supercharges, we construct the Hamiltonian H in the usual manner In the regular representation, this is just the identity operator and hence S 2 1 leads to a trivial spectrum. The Hilbert space has a Z 2 -graded structure, as shown in Fig. 4. These two halves may be dubbed as the "bosonic" and the "fermionic" halves of the space corresponding to the fermion number operator, F = x 2,2 . In this case this turns out to be the projector P as well. Consequently, the M and P operators are the "bosonic" and "fermionic" parts of the Hamiltonian. Note that this system admits no zero-energy ground state, for if |z were such a state, it must satisfy that q |z = q † |z = 0 and there is no such state in H. The Witten index for this system is zero.
We can construct a system with non-empty zero-energy ground states, the first nontrivial example being using the partial symmetries of S 3 1 with dim[H(S 3 1 )] = 9. We choose the one-site supercharges as and we straightforwardly confirm that they are nilpotent. The Hamiltonian is now given by where Figure 4: The grading of the Hilbert space spanned by the partial symmetries of S 2 1 . H b and H f are the "bosonic" and "fermionic" parts of this grading. The grading operator is given by consists of two subspaces, the "bosonic" and "fermionic" subspaces, H b and H f respectively.
First, there is the three-dimensional subspace H 0 comprised of the normalized zero-energy ground states We introduce the fermion number operator F as Clearly, this has eigenvalue 1 upon acting on the states in H 0 , so that the ground states are all fermionic. This makes H 0 ⊂ H f . Second, there is the three-dimensional subspace consisting of the excited fermionic states We denote this space by H ef ⊂ H f , where the label stands for "excited fermionic". Finally, there is the three-dimensional subspace H b , consisting of bosonic states (the fermion number operator F acting on these states gives zero). These states are given by The supercharges q, q † pair up bosonic states in H b and excited fermionic states in H ef , as dictated by the product rules of the underlying SIS. They are excited states with positive energy eigenvalues. Note that so the M, P act as projection operators on bosonic and fermionic subspaces, respectively.
Correspondingly, the energy eigenvalue of the excited states is 1. The fermionic zero-energy ground states are unpaired. Summing up, we have The Witten index ∆ is computed as thus counting correctly the three fermionic zero-energy states in H 0 .
The As a remark we only consider operators that are even under the Z 2 grading, they then form superselection sectors of the theory which are essentially the "bosons" and the "fermions". We cannot have physical states that are a superposition of the two sectors, as they do not have a well defined fermion number. However, it is possible to construct many-body supersymmetric systems with the above grading where entangled states exist, a subject discussed in App. C.

Supersymmetric deformations and Witten index
We are interested in classifying deformations of the Hamiltonian H while preserving supersymmetry. We do not expect these deformations to mix different zero-energy ground states nor to lift them to excited states. As such, we require the Witten index to be invariant under supersymmetry preserving deformations added to the Hamiltonian.
We will classify the supersymmetry preserving deformations into two parts, those that are obtained by deforming the supercharges q to q d resulting in deformed supersymmetric Hamiltonians, H d , and those that are added to the original supersymmetric Hamiltonian H as just perturbations that can possibly lift the ground state degeneracy. We show that the Witten index is left invariant by both these types of deformations.

Deformed supercharges
Consider deformed supercharges of the form with a, b ∈ C. The resulting deformed Hamiltonian is given by with both M d , P d being orthogonal projectors as in the undeformed Hamiltonian, which is recovered with a = b = 1. Now we can compute the zero modes for this deformed system and we find them to be precisely

Supersymmetry preserving perturbations
The simplest class of supersymmetry preserving deformations corresponds to local perturbations that can be added to the Hamiltonian while at the same time preserving the supercharges, q and q † , and commuting with the Klein operator, (1 − 2F ). The only operator that preserves the supercharges in Eq. (2.6) and the Klein operator is provided by the M + P , which is the Hamiltonian itself. Clearly this does neither mix the three zero modes in Eqs. (2.9)-(2.11) nor does it create an energy gap among these states.
We could ask for nontrivial deformations which do not preserve the supercharges but still keep the Witten index unchanged. It turns out that the system is stable to any perturbation built as a linear combination of x 1,1 , x 1,2 and x 1,3 and their hermitian conjugates. This can be seen by their action on the zero-energy ground states in Eqs. (2.9)-(2.11). Accordingly, the system keeps the Witten index intact. Such deformation extends to more nontrivial partial symmetries.
For S n 1 , with arbitrary n, we find that the Witten index ∆ is given by and remains invariant under the class of deformations specified above.
This number for the Witten index can be seen by considering the following supercharge for an arbitrary S n 28) and its conjugate, q † . The zero modes are now given by This index is again left invariant by deformed supercharges of the form The argument goes just as in the case of n = 3. Apart from these deformed supercharges, the system is also invariant under the local perturbations which are a linear combination of x 1,1 , x 1,2 , · · · , x 1,n and their hermitian conjugates, as these operators do not lift the degeneracy of the ground states just as in the n = 3 case.

Supersymmetric systems on a chain
So far, we constructed a system whose supercharge is defined on a single site, so it could be thought of as a one-particle supersymmetric quantum mechanics. Our next step is to extend the construction to a many-body supersymmetric system on a chain, described by a globally defined supercharge. For simplicity, we will choose the homogenous chain such that all sites where each site supports one and the same Hilbert space spanned by the partial symmetries of S 3 1 . In total, dim(H) = 9 N . The Hamiltonian is determined once the supercharges are specified.
We will first need to choose the grading of H, which again can be chosen from many possibilities. We continue adopting the convention that, locally, the sector I is "bosonic" and the sectors II+III are "fermionic". We will now present several examples of many-body systems that can be obtained by selecting different supercharges for this choice of grading.

Non-interacting supersymmetric chain
As a warm-up, we start by considering a simple system in which different lattice sites do not interact with one another. This corresponds to taking the supercharge as where, by definition, θ's are anticommuting and nilpotent variables, {θ i , θ j } = 0, that are built out of q's and q † 's. We can concretely realize these variables using the q's in Eq. (2.6) as follows where the fermion number operator F j was defined in Eq. (2.12). The variable Eq. (3.2) can be thought of as the well-known non-local Jordan-Wigner transformation of the local q i variables.
The purpose of this procedure is to ensure that the θ i 's on adjacent sites anticommute. We will use these variables repeatedly in the rest of this paper.
The Hamiltonian defined by the supercharge Eq. (3.1) is The total fermion number operator F in this case is given by F = N j=1 F j . Therefore, the Z 2 -grading operator is given by It is easy to see that this operator commutes with the Hamiltonian H, thus forming superselection sectors as in the one-particle case. It also anticommutes with the supercharges Q and Q † only when N is odd. Henceforth we assume that N is odd.
Clearly, the N -site chain Hamiltonian Eq.
For simplicity, we assume the chain to be homogeneous and set a i = 1 for all sites i = 1, · · · , N . The form of H i in terms of commuting orthogonal projectors, M i and P i (recall that M 2 i = M i and P 2 i = P i ), results in integer eigenvalues between 0 and N . Then, we can write the energy spectrum of the chain as These spectra exhaust all the possible states of H as N j=0 Degeneracy which, as we have seen, is the dimension of the total Hilbert space H.  At finite temperature β −1 , the partition function is given by

Long-range interacting supersymmetric chain
We now demonstrate how a model of an interacting N -site system can be constructed. The model is associated with long-ranged supercharges and Hamiltonian. We can, however, also construct interacting models but with local supercharges and Hamiltonians. We relegate them to App. A. We emphasize that the MBL property studied for the supercharges in this section is also shared by the supercharges in App. A as those systems continue to possess the local integrals of motion that are possessed by the long-range interacting supercharges in this section. We study the long-ranged interacting supercharges merely for the simplicity of computations.
Consider the following choice of supercharge which is just a product of the local supercharges at each site. This is clearly a nilpotent operator and hence generates a supersymmetry algebra. The resulting Hamiltonian is given by Though interacting, the resulting Hamiltonian is integrable; there are N local integrals of motion given by We can organize the Hilbert space in terms of the cohomologies of nilpotent Q and Q † .
The subspace H 0 of zero-energy states is spanned by solutions of Q|Z = Q † |Z = 0. We see that they are labelled by the product states of the following types. The first type of ground states have at least one local zero-energy state on an individual site: where the ellipses denote any of single-site boson or fermion excited states. There are 9 N − 6 N many such states. The second type of ground states is built from the mixture of single-site boson and fermion excited states with at least one local fermion excited state: 14) The number of such states is precisely 3 N + 3 N for H b and H f , all with eigenvalue 1. Although they are entangled eigenstates of the operator Q + Q † , note that they are not entangled as eigenstates of the Hamiltonian. This can be understood as arising due to the fact that Eq. (3.14) is a superposition of a bosonic and a fermionic state (except in the even N case, when this state is an eigenstate of the fermion number operator). The Hamiltonian of this long-range interacting system then has only product states as eigenstates. 5 The total number of eigenstates is the number of ground states plus the number of excited states, which is equal to 9 N , the total dimension of the Hilbert space, dim H = 9 N . Note that the spectrum of this system is independent of N and is given by the two eigenvalues 0 and 1.
The partition function can be easily computed for this system and is found to be

Supersymmetric deformations and Witten index
The Z 2 -grading Klein operator W is The supercharges Q and Q † anticommute with W .
The Witten index ∆ is defined as the trace of the Klein operator. We can count this index from the ground states we identified above and find precisely −3 N for arbitrary N .
This can be easily seen by considering the form of the states enumerated in Eq. (3.12) and Eq. (3.13). In each of these product states, there are an equal number of bosonic and fermionic states. The only state which is unpaired is the product state |z 1 , z 2 , · · · z N −1 , z N made of oneparticle ground state at every site. As each of these local zero modes are fermionic (recall from Eqs. (2.9)-(2.11)), all these states are fermionic. One can easily confirm that the excited states are paired between bosonic and fermionic states, with multiplicity one. 5 It is however possible to construct supercharges resulting in supersymmetric Hamiltonians that do preserve the Klein operator, have entangled eigenstates and a Witten index different from −3 N . An example will be discussed in App. C.
As in the one-particle case, we are interested in classifying supersymmetry preserving deformations in the many-body setting. Such deformations are defined by continuous perturbations of the Hamiltonian that commute with the supercharge Q given by Eq. (3.10), its adjoint Q † , and the Z 2 grading Klein operator W in Eq. (3.16). We split these up into those that can be added as perturbations to the supersymmetric Hamiltonian and those that are obtained by deforming the supercharge Q as in the one particle case.

Local and quasi-local supersymmetry preserving perturbations
On a chain, we can deform the system in a variety of manners. First, we can deform the system on each site. Such deformations are given by where C(i) is a site-dependent function. Obviously, the system is invariant under these singlesite deformations as they do not change the Witten index ∆.
Next, we can also deform the system over two sites. These deformations take the form Continuing in a similar manner, we can also deform the system over multiple sites; these deformations are supported on several sites and take the form, where k i=1 α i = 0 and the site-dependent coefficient function C k (i 1 , · · · , i k ) is taken to be suitably quasi-local. Such operators are again diagonal in the eigenbasis of this system and hence the Witten index is left invariant. These operators account for all the allowed deformations to this system.

Deformed supercharges
We can introduce a deformation of the supercharge as follows where each of the local deformed supercharges are given by Eq.

The spreading of quantum information
So far, we focused on the spectrum and Witten index of the supersymmetric system on a chain. Here, we dwell on the time evolution of many-body entanglement. This is captured by correlations functions of various time-orderings. More specifically, we will compute outof-time-order correlators (OTOC) and study whether the system scrambles and equilibrates, and, if so, how it does it. We do this for a prototype model, consisting of an interacting disordered system built from Eq. (3.10) that exhibits a many-body localized phase which is supersymmetric. The solvability of this model, as we have seen for the spectral analysis in the previous section, is a remarkable feature brought by the supersymmetric nature of the SIS we utilized in the construction. This will allow us to proceed with analytic computations for the OTOC.

Slow scrambling
First of all, we introduce a quenched disorder in the system by dressing the supercharge in Eq. (3.10) as where J i are real-valued time-independent random variables that can be thought of as analogous to a static random on-site potential. One could generalize this choice by restricting the product to subsets of sites, rather than including all sites, as we shall present in Appendix A.
For now, however, we work with this simpler choice. The results should not depend on this choice.
The many-body interacting Hamiltonian built out of this supercharge is To probe scrambling behavior of this model, we should compute the correlator in Eq. (1.2) and study its time dependence, as discussed in the Introduction. We choose local supercharges as local operators, W (t) = q i (t) and V = q j (0), with i = j. Moreover, we will set β = 0 for simplicity, since we are primarily interested in highly excited states, We prefer to compute C(t) rather than the OTOC (usually considered in this context), but of course the results are independent of this choice.
It is possible to show that the time evolved operator q i (t) is given by and, consequently, All the information about disorder is contained in the argument of the cosine, as a result of the "on-site disorder". In such circumstances, it is reasonable to absorb the effect of randomness into a single variable, defined with its probability measure as where J is a constant. Performing the disorder average leads to (4.7) The behavior of Eq. (4.7) and Eq. (4.8) are shown in Fig. 6. Notice that no N -dependence appears in the time-dependence, apart from the trivial one in the normalization factor of the commutator.
Instead of the quenched quantity we have computed, C(t) J , one could first average over realizations of the couplings and then take the expectation value. While the two procedures generally lead to different results, one can readily perform the computation on the reverse order and see that in this simple case they provide the same answer, that is, the averages commute.
Under what conditions the two procedures are (in)equivalent is an interesting question, that we leave for future studies.
In both choices of the ensemble, the early-time behavior is given by

OTOC-EE theorem from partial symmetries
Recently, a connection between the decay of the OTOC taken in a thermal equilibrium state and the growth of a certain entanglement entropy was proposed for a system quenched by an arbitrary operator [53]. To formulate the exact statement, assume a system described by (4.10) We can state a modified result in terms of our formalism that will supply the story we are developing with further insights. This is accomplished by requiring the quench operator O to act only on the subspace spanned by the partial symmetries of S 3 1 . In other words, we demand {x i,j } to form a complete set for the quench operators we may consider. This is certainly a restriction, since an arbitrary operator acting on the full Hilbert space cannot in general be expressed in terms of partial symmetries. However, this restriction will also provide hints for expecting slow scrambling in any supersymmetric model constructed out of SIS algebras.
As an example, we will verify the partial symmetry version of the OTOC-EE theorem for Eq. (4.2). To this end, we assume again, for simplicity, that the system is at infinite temperature such that the right-hand side of Eq. (4.10) is reduced to As a first step, suppose the initial state to be a maximally mixed state, where ρ(0) ∝ 1. We then quench the system at the first site with To trace out B, we first need to compute the eigenvectors of the local Hamiltonian h j ∝ (P j + M j ). Recall from previous sections that these are given by (4.14) With this in mind, it is straightforward to show that (4. 15) In order to compute the OTOCs in Eq. (4.11), we use the hypothesis that the partial symmetries form a complete set on B, i.e. any quench operator in our model can be expressed in Thus, it is a straightforward, albeit tedious, exercise to show that Adding all the above contributions, where the constant c can be set equal to one by choosing a convenient normalization for the partial symmetries. This completes the check of the theorem.

Discussion: a supersymmetric MBL phase
We now present a heuristic explanation for the results encountered above. First of all, we should remark that, by a localized phase, we mean that decoherence does not occur, that is, long time dynamics does not hide all the information about the initial state. This implies that out-of-equilibrium "atypical" initial states do not evolve into equilibrium "typical" states.
The possibility of keeping coherence is determined by the partial freeze on the process of scrambling, induced by the disorder in all sites of the chain.
The localization phenomenon is clearly a breakdown or violation of the ETH [58,49,[59][60][61], which asserts that all many-body eigenstates of a given system are thermal if all its initial states are able to thermalize, which is supposed to provide a quantum version of thermalization.
Since our Hamiltonian is integrable, it was expected to violate the ETH due to the existence of many local integrals of motion [51,[62][63][64][65], therefore defining a (supersymmetric) many-body localized phase. We stress that, as our system is finite-dimensional and hence does not allow mean field approximation, the localization is the MBL driven by many-particle effects, not the Anderson localization driven by single-particle effects.
Moreover, the supersymmetric many body system we consider is an example of a full-MBL (FMBL) phase, which is a term reserved for a system where all the initial states fail the ETH.
By the statement of ETH, we require that where O is a local operator, · eq denotes the equilibrium value, and the labels C and M C denote the expectation value in a canonical and a microcanonical ensemble, respectively. By choosing the local operator to be M 1 in our setup, we can easily verify that the equilibrium value is not the same as the average in a canonical or a microcanonical ensemble for any initial state, ρ = |ψ ψ|. This ratifies our claim that the SUSY many body systems we consider are an example of FMBL systems.
It is instructive to analyze the situation from other points of view. In Section 4.1, the power-law decay of the OTOC or, alternatively, the growth of C(t) being faster than t in Eqs. (4.7) and (4.8), is a feature encountered in models that present a many-body localized phase [52][53][54][55][56]. The absence of a term linear in time slows down the growth of local operators, preventing the system from achieving thermalization.
It is also worth noting that, as observed in [53], the OTOC can distinguish between MBL and AL phases. For Anderson localization, the OTOC is constant, a feature that can be easily seen in our setup: the non-interacting Hamiltonian in Eq. (3.3) commutes with any local supercharge, hence q i (t) = q i and [q i (t), q j ] is independent of time, leading to a constant OTOC.
To complete the picture, we try to give a heuristic explanation for the absence of fast scrambling in our models. From the analysis on Sec. 4.2, it was possible to relate the Rényi entropy with correlators involving partial symmetries, a fact that came from the restriction of O to be written out of the {x i,j }. Notice that a generic operator acting on the Hilbert space of the system has many more degrees of freedom, by which we mean that the space of linear operators on a local Hilbert space cannot be generated or spanned using only the elements of the partial symmetries, S 3 1 . Nevertheless, the supersymmetric Hamiltonians we can consider are ultimately written in terms of partial symmetries, which form a closed algebra. This means that the entanglement between all possible degrees of freedom one could generally have is largely reduced due to the presence of partial symmetries, which results in the possibility to retain all the information about the initial state as time evolution goes on. The net result is a many-body localized phase as we would naturally expect in this setting.
Finally, we want to address the question of supersymmetry breaking due to finite temperature and quenched disorder effects. To avoid confusion, we stress that, although we are borrowing the terminology of equilibrium physics, the concept of temperature is not well-defined since the system does not achieve thermal equilibrium. Thus, β −1 should be considered as a characteristic energy scale such that β = ∞ probes the ground state and β = 0 probes highly excited states. As discussed in the Introduction, the ground state energy can be thought of as an order parameter of spontaneous supersymmetry breaking, since after eliminating auxiliary components, it results in a vanishing vacuum expectation value for the Hamiltonian if supersymmetry is unbroken, or a positive value if supersymmetry is broken. At finite β, in analogy to equilibrium statistical mechanics, we may define a "thermal" vacuum by utilizing the thermo field double formalism in such a way that the vacuum expectation value of an operator is equal to its thermal average. For an good overview on this subject, we refer to [66].
At quenched disorder, however, J is a random variable and so physical observables should be averaged over the disorder ensemble. Since we are interested in highly excited states, we will consider the case where β is small enough, which leads to This result is valid for both measures we considered previously. It is clear that the disorder effects also contribute to supersymmetry breaking along with the finite temperature effects.
Nonetheless, this consideration shows that supersymmetry is unbroken at infinite temperature, like what happens at zero temperature.
All of these discussions point to the picture that the many-body system we constructed is a slow scrambler and can be used for studying MBL phases. In addition, we emphasize that this localized phase is also supersymmetric and protected by the corresponding Witten index that is independent of β due to the finite-dimensionality of Hilbert space: 7

Discussion and Outlook
In this paper, we have introduced how symmetric inverse semigroups can be used to realize supersymmetric algebras by constructing supercharges out of their elements. By choosing the SIS to be S 2 1 , the simplest example, we recover the usual fermions and the supersymmetric many-body systems constructed in [35,42]. 8 Other choices of SIS lead to non-trivial, novel models. The elements of the SIS algebra allowed to consider graded Hilbert spaces and, for a given SIS other than S 2 1 , there are several choices for this grading. For a given grading, the many-body system constructed is integrable and characterized by an invariant Witten index protected by supersymmetry. We can also construct non-integrable systems by using different gradings, as shown explicitly in App. B. We emphasize that we do not know of any nonsupersymmetric analogs of these systems. We can also realize para-supersymmetric systems using the SIS algebras, as done in App. D.
The fact that these many-body systems are integrable provides a favorable hunting ground for non-thermal states. In this spirit, by introducing disorder we find many-body localized states, which were diagnosed by the behavior of an out-of-time order correlator of local operators. The supersymmetry helped us to compute the OTOC analytically and establish the MBL property, as opposed to the usual MBL literature where most of the computations are numerical.
The supersymmetry algebras and hence the supersymmetric many-body systems con- 7 Recall that N was chosen to be odd. 8 This is shown in detail in App. A. structed here are made up from the elements of S 3 1 . We can build a whole class of models by using the elements of other SISs, namely S n p with n > p. By going to higher n, we increase the dimension of the Hilbert space as the dimension of the algebra made from the elements of S n p is given by n p 2 · p!. This makes the computations tedious but could nevertheless yield surprising results.
We would like to see if these systems have anything to do with quasicrystals. This can be especially seen once we have a hold on the specific inverse semigroups that describe a chosen quasicrystal. Given that our system models MBL phases by exhibiting slow scrambling, we could ask the question if this paves a way to experimentally realize MBL phases on quasicrystals. 9 This problem, however, requires a more systematic study where we construct a system that is invariant under the relevant inverse semigroup corresponding to the quasicrystal. The fact that this system describes the dynamics of a graded Hilbert space and thus a supersymmetric phase, and the fact that it is built very much like a spin system, makes us think of a possible way to realize such a system in the lab. This would account for a table-top experiment for graded Hilbert spaces and a phase of matter characterized by the Witten index and protected by supersymmetry.
Another question to be explored in this setting is whether there is a way to introduce chaotic behavior in such systems and possibly turn it into a toy model for holographyà la SYK model [70] (or its supersymmetric version [71,72]). Toward this end, we are currently looking at introducing open systems or considering subsystems inside a closed system in this setup to induce a linear growth in the OTOC's for early times and an exponential decay for later times.
The supercharge we worked with in this paper is a very simple one which allowed for simple computations. However, as shown in App. B, there are a number of non-trivial supercharges that produce local interactions and are non-integrable. The study of the OTOC's in these systems could show signatures of thermalization. Such considerations can also help us study the transition from localized to thermalized phases in this setting. We are currently working on these issues.
Integrability features frequently in the construction of these supercharges using the SIS 9 MBL phases on quasi-periodic systems have been studied recently [67,68]. This is an extension of a long list of papers about the emergence of localization on quasi-periodic potentials on the one-dimensional chain [69].
algebras. 10 We can ask the question if it is possible to find a Lax pair to describe this system as it is done for the Heisenberg spin chain systems [73]. This will possibly give us solutions to the Yang-Baxter equation and shed more light on these systems.
We could also explore higher-dimensional versions of the models presented in this paper to find more interesting features. The robustness of MBL states to decoherence make them good candidates for engineering quantum memory devices, which might mark a testing ground for supersymmetry, partial symmetries, and SIS. We hope to come back to these explorations in the future.

A Other supersymmetric chains
In supersymmetric systems, the Hamiltonian is specified uniquely by the choice of supercharges. In Sec. 3, we discussed two models of supersymmetric chains: one non-interacting and another long-range-interacting. There are actually many more possible choices of supercharges that lead to interacting systems. In this appendix, we illustrate this point by constructing three alternatives.

A.1 Supercharges of type I
The first example concerns three-site interactions in the total supercharge 10 Non-integrable supercharges can also be obtained, as discussed in App. B.
which can also be recast in terms of the underlying q i variables as where P i = q † i q i . Here, we recall that the q i 's on adjacent sites commute with each other. Using the relations q 2 i = 0, q i P i = q i and P i q i = 0 for all sites i, it is easy to check that Q 2 I = 0.
One finds that the many-body Hamiltonian is a local expression containing terms of the form Note that for a supersymmetric system, we can solve for the spectrum by looking at the Hilbert space on the i-th site. Since they are defined on a site, we will dub these modes as local zero modes, local bosons and local fermions . The many-particle states are then labeled by the local bosons, the local fermions, and the local zero modes. An example of an N -particle state is We can now compute the (global) zero modes, which are denoted by |Z and satisfy There are many such zero modes for this system, arising in three different ways. We could have product states of the form where some of the sites are filled up with local zero modes and the rest are filled up with local bosons or local fermions, or we could have product states where all the sites are filled up with either local bosons or local fermions. Finally, we could have entangled states with the sites filled up with just the local bosons and local fermions.
We now list all these states for the simple case in which the number of sites is N = 4. The total dimension of the Hilbert space in this case is 9 4 . In what follows, we list out all possible forms of the ground states.
Product states like (A.7) These states can take the following forms We also have A similar set of states exist with z 4 on site 4. Thus the total number of such states is 12 · 3 4 .
There are no entangled states, that are ground states and contain local zero modes on the sites.
Product and entangled states like (A.8) In order to find these states, we will consider the most general state as an Ansatz and allow Q I and Q † I to act on it. When the number of sites is N = 4 this state takes the form with coefficients living in C. This provides both the entangled state and the product states.
The former are given by The product states are those which are not part of the entangled sector and there are (2 4 − 6) × 3 4 of them. Thus we have accounted for all the ground states of the system and we find this number to be 5913 states.
The total Hilbert space dimension is 9 4 = 6561. Out of these, as we have just seen, 5913 states are ground states with zero energy. The remaining 648 states are excited states, which we will now write down. The key point is that they are all entangled and have two energy eigenvalues, E = 1 and E = 2. E = 1 states These are given by for a total of 4 · 3 4 = 324 states. As eigenstates of the Hamiltonian for this system these states are not entangled but they are as eigenstates of the operator Q I +Q † I . Another way of seeing is that these entangled states are not eigenstates of the Klein operator which commutes with the Hamiltonian. Thus the four eigenstates with E = 1 are given by |z 1 , b 2 , b 3 , b 4 , |z 1 , f 2 , f 3 , f 4 , |b 1 , b 2 , b 3 , z 4 and |f 1 , f 2 , f 3 , z 4 . As can be easily seen two of these states are fermionic paired with two bosonic states. These product states now have a well defined fermion number and hence are eigenstates of the Klein operator as it should be. E = 2 states These are given by for a total of 4 · 3 4 = 324 states. These eigenstates do not have a well defined fermion number and so the eigenstates of the Hamiltonian with a well defined fermion number can be written as . Two of them are fermionic states paired with two bosonic states as expected.
We have thus accounted for all the 648 excited states and obtained the full spectrum of the theory for the system on N = 4 sites. The MBL property We expect the supercharges of type I to result in systems that possess MBL states as they have local integrals of motion as we have seen above. These provide a shortranged version of the MBL systems as opposed to the unrealistic long-ranged supercharges studied in the main text. The same holds for the supercharges of types II and III presented below.

The Witten index
Relation to [35] One of the earliest works on supersymmetric spin systems can be found in [35]. These are fermion lattice models described by the fermion creation and annihilation operators a * i and a i , respectively, acting on a site i. The supercharge is given by The nilpotency of the supercharge Q follows from the usual algebra of the fermion creation and annihilation operators. A supercharge similar to this can also be written in our framework of supersymmetry from SISs and goes as follows with θ being the non-local, anticommuting variable given by Eq. (3.2). When we choose S 2 1 as our SIS, we recover the supercharge of [35], which is then a special case of our models.
The systems built out of these supercharges continue to be integrable with as many local integrals of motion as the number of sites and are hence expected to be in the MBL phase as is the case with our other integrable models.
Remark After the first version of this article appeared on the arxiv, we found out that H.
Moriya had independently come to the same conclusion [74] of the existence of MBL phases in the supersymmetric spin models of [35], using the techniques of C * -algebras.

A.2 Supercharges of type II
Another possibility is given by The anticommuting variables θ i are given by Eq. (3.2). The Hamiltonian from this supercharge is given by Here we have suppressed the disorder terms of the supercharge, which can however be reintroduced straightforwardly. Relation to [75] It should also be noted that this model is precisely the same as the 1-dimensional lattice fermion model with N = 2 supersymmetry introduced in [75], if we choose the input SIS to be S 2 1 . The mapping is immediate when we set the disorder coefficients to 1 and choose the fermion operator of [75] to be the local supercharge, q i here. For the S 2 1 case this is precisely taken as q i = (x 1,2 ) i to give a many-body supercharge where we have adopted the notation in [75] < i > to denote the sites neighboring site i. This differs from the previous supercharge in (A.12) at the boundaries. More precisely We also note that projector M i is the same as (1 i − F i ) which is precisely the hardcore fermion condition as F i in this setting is the fermion number operator. The local Hamiltonian is  16) which is precisely the same as Eq. (6) of [75]. Thus, our models based on other SISs represent a generalization of the fermion supercharges considered in [75] and subsequent works. As the models in [75] were shown to be related to lattice versions of the (1+1)-dimensional Thirring models with four-fermi interactions we can ask the question what the above models mean for higher SISs. This interesting direction also vindicates our position on the use of SISs to create supersymmetry algebras and supersymmetric many-body systems as they possibly provide lattice realizations of superconformal field theories in the continuum. We will explore these connections in papers to follow.

A.3 Supercharges of type III
Finally, let us introduce yet another supercharge that results in an interacting system with properties quite similar to the ones in the previous sections. This is given by The θ i 's are again given by (3.2) and Again, we have suppressed the disorder terms of the supercharge. This Hamiltonian is integrable, just like the previous two cases, with conserved quantities given by The resulting models can be thought of as hardcore boson models as opposed to the hardcore fermion models introduced in [75] and for the Hamiltonian in (A.13).

B Non-integrable systems
Here we see how one can obtain a non-integrable system using the SIS algebra. A step in this direction is important from the point of view of thermalizing the integrable supersymmetric systems considered in the main text. To do this, we first look at the possible gradings of the Hilbert space spanned by the elements of the SIS algebra. We will work with the S 3 1 case as before and later we will look at the higher n case of S 4 1 . Note that we should start with at least n = 3 as there is only one possible grading with n = 2.

Different gradings
The choice of the supercharge in Eq. (2.6) corresponds to the grading of the Hilbert space shown in Fig. 5. This is one of three possible choices of gradings of the Hilbert space spanned by the elements of the S 3 1 algebra. We could have instead considered a grading shown in Fig.  7, which gives the superchargẽ Figure 7: Another grading of the Hilbert space spanned by S 3 1 elements. H b and H f again denote the "bosonic" and "fermionic" parts respectively.
Note that we have changed what we call a "boson" and a "fermion" by changing the grading of this Hilbert space. The Hamiltonian that we obtain by using this supercharge is with the projectors to the new "bosonic" and "fermionic" sectors given bỹ The other possible grading corresponds to calling the sector III as the "bosonic" part and the sector I+II as the "fermionic" sector. This comes equipped with the corresponding supercharge. For our purposes it is enough to look atq defined in (B.1).  I + II  III + IV   I + III  II + IV   I + IV  II + III   I  II + III + IV   II  I + III + IV   III  I + II + IV   IV  I + II + III There is a supercharge corresponding to each of the gradings listed in the table.
Constructing the non-integrable many-body system We now show a choice of supercharge that helps us obtain the non-integrable supersymmetric system. Until now we have only shown supercharges that resulted in an integrable system with the conserved quantities being given by the local Hamiltonians, h i for sites i ∈ {1, · · · , N }.
By writing down terms that include interactions between the supercharges obtained by the different gradings of the Hilbert space we can obtain the desired supercharges.
We go back to the S 3 1 example, with an analogous procedure for the case of other n and p. Consider the supercharge that results in the non-interacting system, given by Eq. where the anticommutingθ i variables in terms ofq's in (B.1) are as follows with the "fermionic" projector ,P k =q † kq k . We can now consider the following family of supercharges with F an invertible element written as a function of the local supercharges q i that describes the grading corresponding to the one shown in Fig. 5. It is readily seen that (Q ) 2 = 0. This supercharge involves an interaction between the two types of grading and there is no longer a unique local Hamiltonian that forms an integral of motion thus rendering the system to be non-integrable.
An example of a possible choice for F is given by F = e aQ = 1 + aQ, (B.8) with Q corresponding to one of the many supercharges of the grading corresponding to the one in Fig. 5 such as the one shown in (A.1).

Another route
The above method of similarity transforming the superchargeQ with F (Q) is not the only way to obtain an interaction between the two different gradings of S 3 1 . Here we present an example of a supercharge that leads to a non-integrable system but is not obtained using the above procedure.
Consider the following supercharge where q andq are given by Eq.(2.6) and B.1 respectively. It is easy to see that this operator is nilpotent and satisfies the supersymmetric QM algebra. We have chosen just two sites for simplicity, it is easy to generalize to an arbitrary number of sites. The Hamiltonian obtained from this supercharge and its adjoint is given by H = q † 1 q 1 (x 2,2 ) 2 +q † 1q 1 (x 3,3 ) 2 + M 1 M 2 + (x 2,2 ) 1 M 2 . (B.10) Other than the trivial identity operator there is no conserved quantity on site 1. Thus this system is non-integrable. It should also be noted that there is no well defined notion of a "boson" and "fermion" on this Hilbert space as we have mixed two different gradings of the same Hilbert space.
We will use f to denote the two sets together. E = 1 states These are given by giving a total of 2 × 3 N non-zero energy states. Note that the entangled states are fermionic eigenstates of the Klein operator. Thus the bosonic and fermionic states get paired by the supercharges. E = 0 states There are several types. Those that are a mix of local "bosons" and "fermions" are |f 1 · · · , b i 1 , · · · , b i k , · · · , f N , N k where the f i 's now denote one of the 6 basis elements of the local II + III sector. These amount to 9 N − 3 N − 6 N states.

D Para-supersymmetric systems from SISs
Supersymmetric systems can be thought of as dynamics on graded Hilbert spaces. There is no reason why we should stop at just a single grading of the Hilbert space. This leads us to the notion of para-supersymmetry, which is dynamics on a Hilbert space which includes more than one grading. The simplest is given by r = 2 para-supersymmetry which is a generalization of the usual r = 1 supersymmetric systems, and where r defines the para-supersymmetric algebra through its generator q r+1 = 0.
These types of algebras for r = 2 have been studied in the past [76,77] and generalized to higher r in [78]. Such systems have a corresponding generalization of the Witten index and topological interpretation for them [79][80][81].
In the following we recall the version of the r = 2 para-supersymmetry algebra as given in [76] and show how they can be obtained from the SISs algebra.
The r = 2 algebra acts on a Hilbert space, H with the grading which naturally generalizes for arbitrary r. The para-supersymmetry is generated by q, q † and H with the relations q 3 = q †3 = 0, (D.1) [H, q] = H, q † = 0, (D.2) Their actions on the graded Hilbert space H is given by are the zero modes annihilated by both q and q † . The remaining are excited states and are related by the p-supercharges, q and q † .
We can use similar gradings for SISs with other n and p. We cannot use S 2 1 as there is not enough room for the grading. We could start with S 3 1 , but that gives a trivial Hamiltonian just as S 2 1 gave a trivial Hamiltonian in the supersymmetric case.
This realization was done on a single point and can be extended to the many-body case on N sites by a simple global supercharge Q = q 1 · · · q N . (D.7) However we can construct more non-trivial ones as done for the supersymmetric case in App.
A, but we do not do this here.
The MBL property These systems continue to be integrable as they have as many local integrals of motion as the number of sites in the system. Thus we expect them to exhibit the MBL phase just as the supersymmetric systems discussed previously.