## Many-body-localization induced protection of symmetry-protected-topological order in an XXZ spin model

Yoshihito Kuno

Department of Physics, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan (Dated: November 1, 2019)

There is a counter-intuitive expectation proposed by Huse et al [Phys. Rev. B 88 , 014206 (2013)], and Chandran et al [Phys. Rev. B 89,144201 (2014)]: Localization protects quantum order of groundstate even in high excited eigenstates. In this work, we numerically investigate the localization protection for symmetry-protected-topological (SPT) order by considering a modified XXZ spin model, related to an interacting Su-Schrieffer-Heeger (SSH) model. We systematically study how Many-body-localization (MBL) protects SPT for different disorder types and different forms of interactions. A certain disorder leads the clear degeneracy of the low-lying entanglement spectrum of each excited many-body eigenstates. This fact indicates that existence of edge modes, which is a hallmark of the SPT order, is protected by the MBL in excited many-body eigenstates. In addition, we also report how the MBL protected edge modes in excited many-body eigenstates fade away for a diagonal type disorder.

Introduction.Many-body-localization (MBL) is one of the most active topics in condensed matter physics. Theoretically, beyond the classical studies of the Anderson localization [1, 2], rich physical phenomena have been expected so far [3-9]. Experimentally, recent cold atom experiments have claimed to capture the ergodicity breaking dynamics of the MBL [10-13]. On the other hand, topological state of matter is also one of the important topics in condensed matter physics. Various artificial topological states of matter have been realized in recent experimental systems [14, 15] and theoretical predictions have been simulated. In this work, we study the interplay of MBL and topological state. Recently, Huse et al. [16] and Chandran et al. [17] proposed a unique conjecture about the Anderson localization and MBL: A groundstate quantum order is protected even in excited many-body eigenstates, i.e., the MBL protects the groundstate quantum order. Shortly after proposing the conjecture, Kjall et al. [18] and Pekker et al. [19] have tested the MBL induced protection of ferromagnetic Z 2 -symmetry breaking order by considering a onedimensional (1D) quantum random Ising model. Bahri et al. [20] have also considered a 1D disordered cluster spin model as a test model and confirmed the MBL induced protection of the symmetry-protected-topological phase (SPT)[21, 22]. This remarkable result indicates that the SPT order under the MBL could in principle exist out of equilibrium and hence, cooling to detect the SPT order in real experiments is unnecessary[17, 23]. Furthermore, recently, for the 1D disordered cluster spin model, the detailed phase diagram including co-existence of the MBL and SPT order has been investigated [24]. In addition, the relationship between the MBL and SPT phase for a time-reversal symmetric model has been extensively discussed from the perspective of the matrix product state [25, 26]. However, the following issues require further investigation: (I) Does the conjecture stated in [16, 17] hold universally for various condensed matter models? (II) To what extent in the whole spectrum does the mechanism of the MBL protection of quantum order work? (III) How strong does the quantum order of groundstate remain in many-body excited states? In this work, we consider a modified XXZ spin model relating to the SuSchrieffer-Heeger (SSH) model [27], which is one of the simple models exhibiting 1D non-trivial topological insulating phase corresponding to the SPT order [28], and numerically investigate the conjecture in [16, 17].

In addition, whether the SPT phase in an SSH type model or the AKLT model is protected by the MBL remains an open question [17]. Recently, a study [29] indicated that an XXZ model with a particle-hole symmetric disorder and an SSH-type coupling exhibits no protection of the SPT phase by the MBL. Here, for various conditions of disorder and interactions, we numerically evaluate the protection of the SPT phase and also how effective the arguments in [29] are.

Modified XXZ spin model and interacting SSH model.We consider a 1D XXZ spin model with dimerized coupling:

$$\frac { \alpha \mu \, \epsilon } { \text{er spin} } \quad H _ { X X Z } \ = \ \sum _ { i } \left [ J _ { i } ( S _ { i } ^ { + } S _ { i + 1 } ^ { - } + S _ { i } ^ { - } S _ { i + 1 } ^ { + } ) + U _ { i } S _ { i } ^ { z } S _ { i + 1 } ^ { z } + h _ { i } S _ { i } ^ { z } \right ],$$

where S +( -) i and S z i are spin-1/2 raising (lowering) operator and spin-1/2 z -component operator, i denotes lattice site taking i = 1 , · · · , L , where L is a system size. The model parameters are shown in Fig 1 (a), for i ∈ odd , J i → J o i = J 1 + δJ o i , for i ∈ even , J i → J e i = J 2 + δJ e i , where δJ o e ( ) i is a uniform disorder, [ -δJ 1(2) ; δJ 1(2) ], and h i is also a uniform disorder, [ -δh δh ; ]. U i is regarded as an interaction in a fermion picture and in this work acts only on the i ∈ even links with a strength U . This interaction pattern U i may be implemented by considering a combination of a double well and a spin dependent optical lattice [30]. Using the Jordan-Wigner transformation, the XXZ model of H XXZ can be mapped into

FIG. 1: (a) Modified XXZ spin model with open boundary condition. The red-dotted circles on either ends represent the edge sites of the chain. In SPT phase ( J 2 &gt; J 1 ), the edge mode appears as free spin-1/2 degree of freedom, corresponding to zero energy excitation (in thermodynamic limit). (b) Lowest to 20-th energy spectrum for L = 18 system. (c) 81th to 100-th energy spectrum with the same parameters as in (b).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00297v2-figure-d882c369-8c45-4fa1-9992-a806bc57d89c.jpg)

an interacting SSH model [30, 31] with disorder

$$\begin{array} { c c c } \lambda & \cdots & \cdots & \cdots \\ H _ { S S H } & = & \sum _ { i } \left [ J _ { i } ( f _ { i } ^ { \dagger } f _ { i + 1 } + \text{h.c.} ) \\ & + & U _ { i } ( n _ { i } - 1 / 2 ) ( n _ { i + 1 } - 1 / 2 ) + h _ { i } n _ { i } \right ], \\ \cdots & \cdots & \cdots \end{array}$$

where f ( † ) i is an annihilation (creation) spinless fermion and n i = f † i f i is a number operator. Based on this model, let us note the topological properties. In the clean limit and U = 0 non-interacting case, the model belongs to the BDI class in the topological classification [33, 34]. At half-filling, the SSH model exhibits a non-trivial topological phase. The non-trivial topological phase is known to be the SPT phase protected by the chiral symmetry. The SPT phase remains even for a finite U i [30, 31, 35-37]. The model has time-reversal and particle-hole symmetry [29], then the interacting SSH model is invariant under the chiral (sub-lattice) transformation, f † i → -( 1) i f i , f i → -( 1) i f † i [38]. In a disordered case, since δJ o i and δJ e i disorders do not break the chiral symmetry, the groundstate SPT phase is expected to be fairly robust to them. In contrast, if the disorder h i is diagonal so that the chiral symmetry is broken, the groundstate SPT becomes fragile. Since the XXZ model of H XXZ inherits the topological properties of the SSH model of H SSH [39], we numerically treat the XXZ model of H XXZ . In what follows, for all numerical simulations, we fix J 1 = 0 1, . J 2 = 1, and δh = 0 01[41] except for Fig. 4. .

Phenomenology: Localization protects SPT order.The studies by Huse et al. [16] and Chandran et al. [17] state that a quantum order of groundstate in an MBL system with weak interactions is protected in low energy excited states. The low-energy excitations above the groundstate are expected to be localized due to a disorder. Since very weak interacting (almost non-interacting) excitation particles are in Anderson localization, the excitations are suppressed to extend. As a concrete example, for 1D quantum random Ising model the groundstate is ferromagnetic. The excitation is a kink (domain wall) excitation. The excitation kinks are localized owing to disorder, and therefore, the ferromagnetic domains are stable. This leads to a spin-glass order even in high excited eigenstates [18, 19].

The above phenomenon may apply to the model of H XXZ . For L = even clean system with J 1 /lessmuch J 2 and open boundary, the groundstate of the system exhibits a valence-bond-solid order in the bulk, where the singlet state is formed on J e i -links. The groundstate is in the SPT phase because on the left and right edge sites a free spin-1/2 degree of the freedom exists (See Fig. 1 (a)) corresponding to zero-energy edge mode [40]. The edge modes induce fourfold degeneracy of the groundstate. Here, let us consider 2 J 2 = U case and a single excitation on a single J e i -link. The singlet groundstate on the link excites to three triplet states. The excitation may be regarded as a triplon as in the AKLT model [17]. In low energy excited eigenstates, such triplon excitations are dominant. If we regard the triplon excitation as a quasiparticle, the quasiparticles may be localized under a finite disorder. It is expected that the localization of the quasiparticles does not affect (interact with) the edge modes, and hence, the edge modes are protected by localization of the quasiparticles. This means that the SPT order of the groundstate survives even in excited eigenstates. Here, under the change of U , disorder type, and the target sector of the excited eigenstates, it is interesting how the SPT order sustains or not.

In our H XXZ model, we observe a certain typical signal of the protection of the SPT order by localization. Using the exact diagonalization [44-47], we numerically calculated the low energy spectrum for L = 18 system with open boundary, U = 2, δJ 2 = 3 [42], where edge modes are expected to appears on both sides of the edge. In the clean system, the groundstate is fourfold degenerate. See Fig. 1 (b)-(c), we plot the lowest to 28-th and the 81-th to 100-th spectrum. We captured fairly fourfold degeneracy of the spectrum. As observed in the inset figure in Fig. 1 (b), the deviation is negligibly small despite the existence of the finite size effect. Accordingly, the edge mode degeneracy seems to be protected in low energy excited states.

Entanglement spectrum and the degeneracy.We obtained the signal of the protection of the SPT order from the spectrum degeneracy. To verify this phenomenon in detail, we calculate the entanglement spectrum (ES). Consider dividing the system into two equal halves (A and B partial system). The ES { Λ n /lscript } is defined by Λ n /lscript = -2 log | λ n /lscript | , each λ n /lscript is a singular value of the Schmidt decomposition for the n -th eigenstate, then the

FIG. 2: Low-lying entanglement spectrum for low energy excited eigenstates. (a) Lowest and next-lowest entanglement spectrum in L = 14 clean system with U = 0, where the partial system cut has one free 1/2 spin on the left edge. (b) Lowest and nextlowest entanglement spectrum in L = 14 disordered system with U = 2 and δJ 2 = 3. (c) Lowest to fourth entanglement spectrum in L = 16 clean system with U = 0, where the partial system cut has two free 1/2 spins on both left and right edges in the SPT order. (d) Lowest and fourth entanglement spectrum in L = 16 disordered system with U = 2 and δJ 2 = 3.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00297v2-figure-2b14e4f4-4962-476e-b295-d58004fc590b.jpg)

reduced density-matrix of n -th eigenstate is given by ρ n A = ∑ /lscript ( λ n /lscript ) 2 | /lscript 〉 A 〈 /lscript | A . In what follows, we assume that { Λ n /lscript } is ascending-ordered for /lscript . The degeneracy of the low-lying ES characterizes the SPT order. It corresponds to the number of the edge modes [32]. For an interacting SSH model in the clean limit, the ES of the groundstate has been extensively studied to characterize the SPT groundstate by the ES [31, 43]. In the clean limit, we confirmed that the SPT groundstate in the SSH model of H SSH with J 2 &gt; J 1 and open boundary has two or fourfold degeneracy of the lowest ES depending on the system size L .

Let us now consider the disordered system. We focus on effects of δJ e i disorder in the low excited energy sector, and therefore, switch off δJ o i and h i disorders. A periodic boundary condition is considered. In the calculation of the ES, we virtually cut the system into two equal halves with the system size L/ 2. Here, for L/ 2 = odd , it is necessary to cut one J o i -link and one J e i -link to obtain two equal halves. If the system is in the SPT phase, the partial system has a single free spin-1/2 edge mode on either right or left edge sites. This leads to twofold degeneracy of the lowest ES [31]. Meanwhile, for L/ 2 = even , we cut two J e i -links to obtain two equal halves. In the SPT phase, this leads to two free spin-1/2 edge modes on both right and left edge sites corresponding to fourfold degeneracy of the lowest ES [31]. For L = 14 system with U = 0 and δJ 2 = 0, we first calculated the ES by using many-body eigenstates from groundstate to 100-th excited state. The lowest and second lowest ES, Λ 1 and Λ 2 are plotted in Fig. 2 (a). While Λ 1 and Λ 2 are degenerate for the groundstate, Λ 1 and Λ 2 for many other excited states are not degenerate. This indicates that the SPT order vanishes in low energy excited states except for the groundstate. In contrast, see an interacting disorder case δJ 2 = 3 as shown in Fig. 2 (b). This is a single-shot result for a particular δJ e i disorder distribution. Interestingly, we observed the proliferation of twofold degeneracy of the ES even in low energy excited eigenstates. As shown in inset Fig. 2 (b), the deviation of the degeneracy is negligibly small even in the finite system size. The result indicates the survival of the edge modes even in low energy excited states due to δJ e i disorder. Furthermore, we also calculated the ES of L = 16 system for a noninteracting clean case in Fig. 2 (c). The lowest to fourth lowest ES, Λ , Λ 1 2 Λ 3 and Λ 4 are plotted. As observed in Fig. 2 (a), for the groundstate, the fourfold degeneracy from Λ 1 to Λ 4 appears, while for many other excited states the fourfold degeneracy is not observed. However, as shown in Fig. 2 (d), we calculated the ES for an interacting disorder case δJ 2 = 3 and found the proliferation of the fourfold degeneracy of the ES even in low energy excited states. This indicates for L/ 2 = even system the survival of the edge modes even in low energy excited eigenstates. Thus, the ES calculation indicates that the δJ e i disorder protects the SPT order in the low energy excited eigenstates.

Further, we investigate how the low-lying ES degeneracy depends on δJ 2 and the regime of excited energy density. To quantify the degree of the degeneracy of the low-lying ES, we define the degree of the two or fourfold degeneracy of the ground ES, for L/ 2 = odd system ∆Λ 1 -2 ( n ) ≡ Λ 2 -Λ , 1 for L/ 2 = even system ∆Λ 1 -4 ( n ) ≡ 3Λ 1 -Λ 2 -Λ 3 -Λ , and then define the 4 following measure:

$$O _ { 2 E S } = \frac { 1 } { N. }$$

$$O _ { 4 E S } = \frac { 1 } { N _ { s } } \sum _ { n \in \mathbb { N } _ { E } } \theta ( \epsilon - \Delta \Lambda _ { 1 - 4 } ( n ) ),$$

$$O _ { 2 E S } & = \frac { 1 } { N _ { s } } \sum _ { n \in \mathbb { N } _ { E } } \theta ( \epsilon - \Delta \Lambda _ { 1 - 2 } ( n ) ), \\ O _ { 4 E S } & = \frac { 1 } { \pi \nu } \sum _ { \cdot } \theta ( \epsilon - \Delta \Lambda _ { 1 - 4 } ( n ) ),$$

where N s is a number of eigenstate in the target regime of the spectrum denoted by N E , θ X ( ) is the Heaviside function, and /epsilon1 represents a practical numerical criteria of the degeneracy of the ES [48]. We calculate the disorderaveraged values of O 2 ES and O 4 ES denoted by 〈 O 2 ES 〉 and 〈 O 4 ES 〉 . In the disordered system, 〈 O 2 ES 〉 and 〈 O 4 ES 〉 can be regarded as an order parameter to charac-

FIG. 3: The behavior of the ES degenerate order parameters with varying disorder strength δJ 2 and an interaction form. The behavior of 〈 O 2 ES 〉 in L = 14 system for the low excited energy sector (a) and for the middle excited energy sector (b) The behavior of 〈 O 4 ES 〉 in L = 12 system for the low excited energy sector (c) and for the middle excited energy sector (d) As a whole, the interaction U enhances the saturation values of 〈 O 2 ES 〉 and 〈 O 4 ES 〉 for large δJ 2 . (e) The effects of an odd-link interaction U 1 for L = 14 system with U = 2, δJ 2 = 3 and δJ 1 = 0 for the low excited energy sector. (f) Modulated-disordered interaction case. U i = 2 J i . L = 14 system with δJ 2 = 3 and δJ 1 = 0 for the low excited energy sector. For all cases, 10 2 disorder samples were used. Error bars are small and are included in the point labels.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00297v2-figure-22944ece-832c-4f11-8d19-55f68d49839d.jpg)

terize the SPT order in the excitation energy parameter space. Figure 3 shows various behaviors of 〈 O 2 ES 〉 and 〈 O 4 ES 〉 as varying δJ 2 with δJ 1 = 0.

See Fig. 3 (a), 〈 O 2 ES 〉 is of L = 14 system and we set N E to 10% excited eigenstates from the groundstate in all eigenstates. We found that for various U , 〈 O 2 ES 〉 increases with the increase of δJ 2 , and that interestingly 〈 O 2 ES 〉 saturates up to ∼ 0 7 for large . U . This implies that almost 70% excited eigenstates in the set of N E exhibit a single edge mode on the edge site. The low energy excited eigenstates are fairly ordered in the SPT order with the increase of δJ 2 . Furthermore, for large δJ 2 , 〈 O 2 ES 〉 of large U case tends to be larger than that of small U . The result indicates that the interaction U tends to promote the SPT order for excited eigenstates. We further calculated 〈 O 2 ES 〉 with the same condition but for the middle 10% excited eigenstates in all eigenstates. Remarkably, 〈 O 2 ES 〉 also increases as increasing δJ 2 , the saturate value of 〈 O 2 ES 〉 is somewhat small compared to that in the low energy excited sector. This indicates that the protection of the SPT order is weak in the middle excited energy sector. However, interestingly for the high excited energy sector, the protection of the SPT order tends to be strong again, although the interaction dependence varies. See the supplemental material [49]. Let us turn on the disorder dependent behavior of 〈 O 4 ES 〉 as shown in Fig. 3 (c)-(d). The results are of L = 12 system size. The results in Fig. 3 (c)-(d) have the same tendency as those in Fig. 3 (a)-(b), and therefore, we expect L/ 2 = even system to clearly exhibit the protection of the SPT order in the low of the excited states. In addition, we checked the system size dependence for 〈 O 2 ES (4 ES ) 〉 , see [49].

〈 O 2 ES 〉 is decreasing. The MBL protection is decreasing. The uniform interaction case tends to be close to the claim in [29]. Second, U i = 2 J i case was also calculated as shown in Fig. 3. The parameters are same with Fig. 3 (e). From the result, as a whole 〈 O 2 ES 〉 is smaller than those of the interacting cases in Fig. 3 (a) and (c). Our numerics shows that the MBL protection is subtle.

Other interaction conditions were also calculated. Effect of the odd-link interaction U i = odd = U 1 is shown in Fig. 3 (e). Here, other parameters are U = 2 and δJ 2 = 3. As approaching a uniform interaction U i = U for all i ,

Effects of diagonal disorder and MBL signal.We investigated effects of δh disorder. This disorder breaks the chiral symmetry in the system of H SSH . Hence, we expect that this disorder strongly affects the edge mode in the SPT order and that the SPT order tends to vanish with a certain strength of δh . Figure. 4 (a) shows a typical δh dependence of 〈 O 2 ES 〉 for L = 14 system with δJ 2 = 3, δJ 1 = 0. For the low and middle sector of the excited energy, 〈 O 2 ES 〉 decreases with the increase of δh in various U . There is no clear spectrum phase transition, but a crossover towards the conventional MBL without the SPT order is observed. The protection of the SPT is fairly robust, which is almost suppressed at δh ∼ 1, comparable to J 2 energy scale. We also calculated effects of δJ 1 disorder and obtained a global phase diagram [49]. In addition, we investigated the stability of MBL. To characterize the MBL for each excited eigenstates, we employ the average level spacing ratio 〈 r 〉 [8, 50, 51]. Particularly, to obtain the clear behavior of 〈 r 〉 , we focused on the ∑ L i =1 S z i = 0 sector in the Hilbert space and on the middle excited eigenstates [50]. In the calculation of 〈 r 〉 , consider the spectrum { E i } (in ascending order). Each level spacing in { E i } is given as r k = [min( δ ( k ) , δ ( k +1) )] / [max( δ ( k ) , δ ( k +1) )], where δ ( k ) = E k +1 -E k , and then the value of 〈 r 〉 is obtained by averaging over disorder samples. From the value of 〈 r 〉 , we can estimate localization in the target sector of the excited eigenstates: For the MBL, 〈 r 〉 ∼ 0 386, . corresponding to the Poisson random matrix ensemble. The results of 〈 r 〉 are shown in Fig. 4 (b)-(c), where we

FIG. 4: (a) δh dependence of 〈 O 2 ES 〉 for L = 14 system with J 1 = 0 1 and . J 2 = 1. Error bars are small and are included in the point labels. (b) δJ 2 -dependence of the average gap ratio 〈 r 〉 for the middle excited energy sector. (c) δh -dependence of the average gap ratio 〈 r 〉 for the middle excited energy sector. For all results, we averaged over 10 2 disorder samples.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00297v2-figure-2474e2ab-32f4-43dd-90ce-cc1255550547.jpg)

employed L = 16 system. Figure. 4 (b) is δJ 2 -dependence of 〈 r 〉 for a finite moderate δh case ( δh = 0 01, 0 1 and . . 0 2) . and U = 2. For small δJ 2 , 〈 r 〉 takes ∼ 0 53 . the signal of an extended state, corresponding to Gaussian orthogonal ensemble [50, 51], and with the increase of δJ 2 , the values of 〈 r 〉 approach ∼ 0 386. . δJ 2 -disorder induces MBL, i.e., MBL is exhibited when the eigenstates turn into the SPT order as increasing δJ 2 from zero [52]. Figure. 4 (c) displays δh -dependence of 〈 r 〉 with δJ 2 = 3. Even in varying δh where the SPT order is suppressed in excited eigenstates, the MBL signal 〈 r 〉 ∼ 0 386 re-. mains. This indicates that during the crossover behavior in Fig. 4 (a), the MBL sustains.

Conclusion.For a simple modified XXZ model, we numerically investigated the MBL protection of the SPT phase. We found that for a certain interaction pattern, the MBL clearly protects the SPT order in many-body excited eigenstates. However, a uniform interaction situation and the diagonal disorder weaken the tendency of the MBL protection. Although the limitation in this work is the small system sizes accessible to exact diagonalization, our results indicate that the SPT order in the XXZ model possibly appears even in high temperature or out of equilibrium.

Acknowledgments.Y. K. acknowledges the support of the Grant-in-Aid for JSPS Fellows (No.17J00486). Y. K. is very grateful to Jakub Zakrzewski for his critical reading of my manuscript and for providing numerous helpful suggestions that improved my draft, and thanks an anonymous referee for valuable suggestions.

- [1] P. W. Anderson, Phys. Rev. 109 1492 (1958).
- [2] A. Lagendijk, B. Van Tiggelen, and D. S.Wiersma, Phys. Today 62 , 24 (2009).
- [3] R. Nandkishore and D. A. Huse, Annu. Rev. Condens. Matter Phys. 6 , 15 (2015).
- [4] D. M. Basko, I. L. Aleiner, and B. L. Altshuler, Ann. Phys. 321 , 1126 (2006).
- [5] D. A. Abanin and Z. Papic, Annalen derPhysik 529 , 1700169 (2017).
- [6] F. Alet and N. Laflorencie, Comptes Rendus Physique 19 , 498 (2018).
- [7] D. A. Abanin, E. Altman, I. Bloch, and M. Serbyn, Rev. Mod. Phys. 91 , 021001 (2019).
- [8] V. Oganesyan and D. A. Huse, Phys. Rev. B 75 , 155111 (2007).
- [9] S. Iyer, V. Oganesyan, G. Refael, and D. A. Huse, Phys. Rev. B 87 , 134202 (2013).
- [10] M. Schreiber, S. S. Hodgman, P. Bordia, H. P. Luschen, M. H. Fischer, R. Vosk, E. Altman, U. Schneider, and I. Bloch, Science 349 , 842 (2015).
- [11] J.-Y. Choi, S. Hild, J. Zeiher, P. Schaus, A. RubioAbadal, T. Yefsah, V. Khemani, D. A. Huse, I. Bloch, and C. Gross, Science 352 , 1547 (2016).
- [12] A. Lukin, M. Rispoli, R. Schittko, M. E. Tai, A. M. Kaufman, S. Choi, V. Khemani, J. Le`nard, and M. Greiner, o Science 364 , 256 (2019).
- [13] M. Rispoli, A. Lukin, R. Schittko, S. Kim, M. E. Tai, J. Le` onard and M. Greiner, arXiv: 1812.06959 (2018).
- [14] N. R. Cooper, J. Dalibard, and I. B. Spielman, Rev. Mod. Phys. 91 , 15005 (2019).
- [15] T. Ozawa, H. M. Price, A. Amo, N. Goldman, M. Hafezi, L. Lu, M. Rechtsman, D. Schuster, J. Simon, O. Zilberberg, and I. Carusotto, Rev. Mod. Phys. 91 , 015006 (2019).
- [16] D. A. Huse, R. Nandkishore, V. Oganesyan, A. Pal, and S. L. Sondhi, Phys. Rev. B 88 , 014206 (2013).
- [17] A. Chandran, V. Khemani, C. R. Laumann, and S. L. Sondhi, Phys. Rev. B 89 ,144201 (2014).
- [18] J. A. Kjall, J. H. Bardarson, and F. Pollmann, Phys. Rev. Lett. 113 ,107204 (2014).
- [19] D. Pekker, G. Refael, E. Altman, E. Demler, and V. Oganesyan, Phys. Rev. X 4 , 011052 (2014).
- [20] Y. Bahri, R. Vosk, E. Altman, and A. Vishwanath, Nat. Commun. 6 , 7341 (2015).
- [21] F. Pollmann, A. M. Turner, E. Berg, and M. Oshikawa, Phys. Rev. B 81 , 064439 (2010).
- [22] F. Pollmann, E. Berg, A. M. Turner, and M. Oshikawa, Phys. Rev. B 85 , 075125 (2012).
- [23] S. A. Parameswaran and R. Vasseur, Reports Prog. Phys. 81 , 082501 (2018).
- [24] K. S. C. Decker, D. M. Kennes, J. Eisert, and C. Karrasch, arXiv: 1902.02259 (2019).
- [25] T. B. Wahl, Phys. Rev. B 98 , 054204 (2018).
- [26] A. Chan, T. B. Wahl, arXiv:1808.05656 (2019).
- [27] W. Su, J. Schrieffer and A. Heeger, Phys. Rev. Lett. 42 , 1698 (1979).
- [28] J. K. Asb´th, L. Oroszl´ny, and A. P´lyi, A Short Course o a a on Topological Insulators: Band Structure and Edge States in One and Two Dimensions (Springer International Publishing, New York, 2016), Lecture Notes in Physics, Vol. 919 .
- [29] R. Vasseur, A. J. Friedman, S. A. Parameswaran, and A. C. Potter, Phys. Rev. B 93 , 134207 (2016).
- [30] Y. Kuno, Phys. Rev. B 99 064105 (2019).
- [31] J. Sirker, M. Maiti, N. P. Konstantinidis, and N. Sedlmayr, J. Stat. Mech. P10032 (2014).
- [32] L. Fidkowski, Phys. Rev. Lett. 104 , 130502 (2010).

- [33] A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Ludwig, Phys. Rev. B 78 , 195125 (2008).
- [34] A. Kitaev, in Advances in Theoretical Physics: Landau Memorial Conference, edited by V. Lebedev and M. Feigelfman, AIP Conf. Proc. No. 1134 (AIP, Melville, NY, 2009), p. 22.
- [35] M. Yahyavi, L. Saleem, and B. Hetenyi, J. Phys. Condens. Matter 30 , 445602 (2018).
- [36] B. Het´ enyi and B. D´ra, Phys. Rev. B o 99 , 085126 (2019)
- [37] M. Weber, F. P. Toldin, and M. Hohenadler, arXiv:1905.05218 (2019).
- [38] S. R. Manmana, A. M. Essin, R. M. Noack, and V. Gurarie, Phys. Rev. B 86 , 205119 (2012).
- [39] M. P. Estarellas, I. DfAmico, and T. P. Spiller, Sci. Rep. 7 , 42904 (2017).
- [40] In the thermodynamic limit, the edge mode exhibits exact zero-energy excitation. For a finite size case, the energy-excitation is exponentially small with the system size L .
- [41] In our numerics, if we set δh = 0, degeneracy in the spectrum is proliferated. This may originate from the symmetries of the topological model and the non-conservation of the total spin-z component in the system. Such a degeneracy can lead difficulty in evaluating the degeneracy of the entanglement spectrum. Once we inject small but finite δh , the degree of the degeneracy is expected to be lessened.
- [42] In all numerical simulation except for Fig. 4 (b)-(c), the Hilbert space includes all spin z-component sectors. The total spin z-component S z total = i S z i is not conserved.
- ∑ [43] B. T. Ye, L. Z. Mu, and H. Fan, Phys. Rev. B 94 , 165167 (2016).
- [44] M. Noack, S. R. Manmana, AIP Conf. Proc. 789 , 93-163 (2005).
- [45] M. Zhang and R. X. Dong, Eur. J. Phys. 31 , 591 (2010).
- [46] P. Prelovsek and J. Bonca, Strongly Correlated Systems: Numerical Methods (Springer, Berlin, 2013), Springer Series in Solid-State Sciences, Vol. 176.
- [47] D. Raventos, T. Gras, M. Lewenstein, and B. Julia-Diaz, J. Phys. B 50 , 113001 (2017).
- [48] As shown in the inset figures of Fig. 2 (b) and (d), since the deviation of the degeneracy coming from finite size effects is negligibly small, we set /epsilon1 = 0 1 for the calculation . of 〈 O 2 ES 〉 and /epsilon1 = 0 2 for the calculation of . 〈 O 4 ES 〉 .
- [49] Supplemental material:
- [50] D. J. Luitz, N. Laorencie, and F. Alet, Phys. Rev. B 91 , 081103 (2015).
- [51] J. Janarek, D. Delande, J. Zakrzewski, Phys. Rev. B 97 , 155133 (2018).
- [52] For δh = 0 case, 〈 r 〉 indicated less than 0 386. This comes . from the fact that the level spacing ratio is difficult to capture the signal of the MBL when the high degeneracy of the entire spectrum of the model of H XXZ without h i disorder remains. However, we expect since very small δh = 0 01 case clearly . indicates the MBL signal and therefore even for the limit δh → 0, the system with only δJ 2 disorder is in the MBL.

## Supplemental Material

## DEGENERACY OF THE ES IN HIGH EXCITED EIGENSTATES

In the main text, we observed the MBL induced protection of the edge mode in the low and middle excited energy sector. Here, let us observe the protection behavior for in the high excited energy sector. Figure S1 shows various behaviors of 〈 O 2 ES 〉 and 〈 O 4 ES 〉 as varying δJ 2 with J 1 = 0 1, . J 2 = 1, δh = 0 01, and . δJ 1 = 0. We set N E to 10% excited eigenstates below the highest eigenstate in all eigenstates. See Fig. S1 (a), 〈 O 2 ES 〉 is of L = 14 system with J 1 = 0 1, . J 2 = 1. As same with the low and middle excited energy sector, we found that for various U , 〈 O 2 ES 〉 increases with the increase of δJ 2 . However, compared with Fig. 3 (a)-(b) in the main text, the interaction U does not enhance 〈 O 2 ES 〉 for large δJ 2 . This tendency is also true for the behavior of 〈 O 4 ES 〉 with L = 12 system as shown in Fig. S1 (b). As a whole, the protection of the SPT in the high excited energy sector is stronger than that of the middle excited energy sectors. Therefore, the MBL protection tends to be somewhat strong at band edges compared to the center of the band.

FIG. S1: The behavior of the ES degenerate order parameters in varying the disorder strength. (a) The behavior of 〈 O 2 ES 〉 in L = 14 system for the high excited energy sector. We used 10 2 disorder samples. (b) The behavior of 〈 O 4 ES 〉 in L = 12 system for the high excited energy sector. We used 10 2 disorder samples. For the calculation of 〈 O 2 ES 〉 , /epsilon1 = 0 1 and . for the calculation of 〈 O 4 ES 〉 , /epsilon1 = 0 2. Error bars are small . and are included in the point labels.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00297v2-figure-2ea405ce-fd6f-4020-adbb-c1487f2ce848.jpg)

## FINITE SYSTEM SIZE EFFECTS IN THE BEHAVIOR OF THE ENTANGLEMENT SPECTRUM DEGENERACY

We show a typical system size dependence of the degeneracy of the ES. We calculated 〈 O 2 ES 〉 and 〈 O 4 ES 〉 for various system size. The definition of 〈 O 2 ES 〉 and 〈 O 4 ES 〉 is shown in the main text. We set the parameters J 1 = 0 1, . J 2 = 1, U = 2, δJ 1 = 0, δJ 2 = 3 and δh = 0 01, where the . system exhibits the MBL protection for the SPT phase and we used the lowest 10% many-body eigenstates. The result is shown in Fig. S2. Here, L =8,10,12,14 and 16 systems were employed. When the system size is increased, the decrease in the value of 〈 O 2 ES 〉 and 〈 O 4 ES 〉 seems small. This fact indicates that the ES degeneracy characterized by a finite value of 〈 O 2 ES 〉 and 〈 O 4 ES 〉 survives for fairly large system size.

FIG. S2: System size dependence of the degeneracy of the ES. The number nearby the labels represents the number of the disorder samples. The red line is a fitting line determined by all data. The red line is 0 1502 ln(1 . /L ) + 0 9201 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00297v2-figure-85b8b53a-8548-46a3-b162-65bba2119a0b.jpg)

## δJ 2 AND δh DISORDER PHASE DIAGRAM

In the main text, we observed the δh disorder dependence of 〈 O 2 ES 〉 with δJ 2 = 3. Here we observe effects of δJ 1 disorder by calculating a global phase diagram. We employed L = 10 system with J 1 = 0 1, . J 2 = 1, U = 2, and δJ 2 = 3, set N E to 10% excited eigenstates from the groudstate in all eigenstates, and calculated 〈 O 2 ES 〉 by averaging over disorder samples. The result is shown in Fig. S3. For small δh , the δJ 1 disorder does not suppress the value of 〈 O 2 ES 〉 so much. In contrast, the δh disorder more suppresses the value of 〈 O 2 ES 〉 . Therefore, the MBL protection of the SPT order are fairly robust to effect of the δJ 1 disorder. It is natural because the δJ 1 disorder does not break the chiral symmetry of the system of H XXZ and therefore the SPT order is hardly affected.

FIG. S3: δJ 2 and δh disorder phase diagram. We used 10 2 disordered samples.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00297v2-figure-44f06fb0-a2f8-45c6-9e77-3f1f4c21abae.jpg)
Yoshihito Kuno
Many-body-localization induced protection of symmetry-protected topological order in a XXZ spin model
There is a counter-intuitive expectation proposed by Huse {\it et al} [Phys.
Rev. B {\bf 88}, 014206 (2013)], and Chandran {\it et al} [Phys. Rev. B
89,144201 (2014)]: Localization protects quantum order of groundstate even in
high excited eigenstates. In this work, we numerically investigate the
localization protection for symmetry-protected-topological (SPT) order by
considering a modified XXZ spin model, related to an interacting
Su-Schrieffer-Heeger (SSH) model. We systematically study how
Many-body-localization (MBL) protects SPT for different disorder types and
different forms of interactions. A certain disorder leads the clear degeneracy
of the low-lying entanglement spectrum of each excited many-body eigenstates.
This fact indicates that existence of edge modes, which is a hallmark of the
SPT order, is protected by the MBL in excited many-body eigenstates. In
addition, we also report how the MBL protected edge modes in excited many-body
eigenstates fade away for a diagonal type disorder.