  Recent work shows that highly excited many-body localized eigenstates can
exhibit broken symmetries and topological order, including in dimensions where
such order would be forbidden in equilibrium. In this paper we extend this
analysis to discrete symmetry protected order via the explicit examples of the
Haldane phase of one dimensional spin chains and the topological Ising
paramagnet in two dimensions. We comment on the challenge of extending these
results to cases where the protecting symmetry is continuous.
