Matrix Product Symmetries and Breakdown of Thermalization from Hard Rod Deformations

We construct families of exotic spin-1/2 chains using a procedure called "hard rod deformation." We treat both integrable and nonintegrable examples. The models possess a large noncommutative symmetry algebra, which is generated by matrix product operators with a fixed small bond dimension. The symmetries lead to Hilbert space fragmentation and to the breakdown of thermalization. As an effect, the models support persistent oscillations in nonequilibrium situations. Similar symmetries have been reported earlier in integrable models, but here we show that they also occur in nonintegrable cases.