Strong zero modes in integrable spin-S chains

We derive exact strong zero mode (ESZM) operators for integrable spin-S chains with open boundary conditions and a boundary field. Their locality properties are generally weaker than in the previously known cases, but they still imply infinite coherence times in the vicinity of the edges. We explain how such integrable chains possess multiple ground states describing a first-order quantum phase transition, and that the odd number of such states for integer S makes the weaker locality properties necessary. We make contact with more traditional approaches by showing how the ESZM for S=1/2 acts on energy eigenstates given by solutions of the Bethe equations.