Entanglement of exact excited states of Affleck-Kennedy-Lieb-Tasaki models: Exact results, many-body scars, and violation of the strong eigenstate thermalization hypothesis

We obtain multiple exact results on the entanglement of the exact excited states of nonintegrable models we introduced in Phys. Rev. B 98, 235155 (2018). We first discuss a general formalism to analytically compute the entanglement spectra of exact excited states using matrix product states and matrix product operators and illustrate the method by reproducing a general result on single-mode excitations. We then apply this technique to analytically obtain the entanglement spectra of the infinite tower of states of the spin-S AKLT models in the zero and finite energy density limits. We show that in the zero energy density limit, the entanglement spectra of the tower of states are multiple shifted copies of the ground-state entanglement spectrum in the thermodynamic limit. We show that such a resemblance is destroyed at any nonzero energy density. Furthermore, the entanglement entropy S of the states of the tower that are in the bulk of the spectrum is subthermal S proportional to log L, as opposed to a volume law S proportional to L, thus indicating a violation of the strong eigenstate thermalization hypothesis (ETH). These states are examples of what are now called many-body scars. Finally, we analytically study the finite-size effects and symmetry-protected degeneracies in the entanglement spectra of the excited states, extending the existing theory.