  Many-body localized systems in which interactions and disorder come together
defy the expectations of quantum statistical mechanics: In contrast to ergodic
systems, they do not thermalize when undergoing nonequilibrium dynamics. What
is less clear, however, is how topological features interplay with many-body
localized phases as well as the nature of the transition between a topological
and a trivial state within the latter. In this work, we numerically address
these questions, using a combination of extensive tensor network calculations,
specifically DMRG-X, as well as exact diagonalization, leading to a
comprehensive characterization of Hamiltonian spectra and eigenstate
entanglement properties.
