Pith. sign in

REVIEW 3 major objections 4 minor 3 cited by

This paper derives exact edge operators that commute with the Hamiltonian of any integrable spin-S chain, guaranteeing infinitely long edge coherence times even where locality is weak.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-03 17:51 UTC pith:VNCUMAWT

load-bearing objection Solid arbitrary-spin ESZM construction with a real weak-locality story; the main gap is an unproved spectral assertion on the auxiliary transfer matrix. the 3 major comments →

arxiv 2512.07742 v1 pith:VNCUMAWT submitted 2025-12-08 cond-mat.stat-mech math-phmath.MPquant-ph

Strong zero modes in integrable spin-S chains

classification cond-mat.stat-mech math-phmath.MPquant-ph
keywords strong zero modesintegrable spin chainsXXZ chainmatrix product operatorstransfer matricesedge coherenceBethe ansatzboundary strings
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper sets out to show that exact strong zero modes (ESZMs) — operators that commute exactly with the Hamiltonian and are localized near the edge — exist not just for the spin-1/2 XXZ chain but for every integrable spin-S XXZ chain with open boundaries and a boundary field. Because the integrable spin-S chain has 2S+1 degenerate ground states, and an odd number of them for integer S, the new ESZMs cannot have the same sharp locality or the property of squaring to the identity that the spin-1/2 mode has. Instead the authors prove these properties hold in a weaker, Hilbert-Schmidt sense, and argue that this is still enough to force edge autocorrelation functions to plateau at nonzero values — infinite edge coherence times. The construction uses the family of commuting transfer matrices of the model and gives explicit matrix-product forms for the modes. A sympathetic reader should care because this gives exact, parameter-free statements about long-lived edge physics in a class of strongly interacting many-body systems, and identifies why integer spin is structurally different.

Core claim

On the paper's own terms: for any spin S and real anisotropy η, the derivative of the auxiliary-spin-1/2 transfer matrix at the special value u* = iπ/2 yields a Hermitian operator Ψ that commutes exactly with the integrable open spin-S Hamiltonian. In the Hilbert-Schmidt norm, Ψ is normalized to 1 and squares to the identity up to corrections exponentially small in system size, and each local piece Ψ_j of its MPO decomposition decays exponentially in j; however, for S ≥ 1 the spectral norm of Ψ_j does not decay, so locality is strictly weaker than the spin-1/2 case. The authors explain this weakening as forced by the odd number of degenerate ground states for integer S, and show numerically

What carries the argument

The central object is the double-row transfer matrix T^{(S,1/2)}(u) built from quantum-group generators S^z, S^± and diagonal reflection matrices; it commutes with the Hamiltonian for all u. The special value u* = iπ/2 makes the L-operator satisfy L_j(u*) σ^z_a L_j(u*) = σ^z_a, which collapses the transfer matrix to zero and turns its derivative into a sum of local pieces Ψ_j. Writing Ψ as a matrix-product operator with a four-dimensional ancilla, the Hilbert-Schmidt norms are evaluated by a transfer matrix U = (1/(2S+1))Tr(A⊗A); the claim that U has a unique eigenvalue 1 with all others smaller in modulus is what makes the mode normalizable and almost idempotent.

Load-bearing premise

The entire normalizability and near-inversion of Ψ rests on an unproved spectral assertion: that the auxiliary transfer matrix U (and its fourth-power analogue V) has a unique eigenvalue of modulus 1 and a gap below it; the paper says this was 'verified by direct calculation' but gives no derivation.

What would settle it

Compute the spectrum of the four-dimensional transfer matrix U (or V) for a specific S and real η and find an eigenvalue of modulus greater than or equal to 1 other than the Perron eigenvalue; then the claimed exponential decay of ∥Ψ_j∥ and the limit ∥Ψ²−1∥²→0 fail. Alternatively, for the spin-1 chain with boundary conditions of Eq. (59), compute the edge autocorrelator plateau at larger L and show it decays to a value below |c_1|²+|c_2|², contradicting the claim that additional ESZMs exist.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • For every S, the open integrable spin-S XXZ chain with appropriate boundary fields has an exactly conserved edge operator, so local observables overlapping it show nonzero infinite-temperature autocorrelation plateaus.
  • The weaker locality for S≥1 is not a technical accident but forced by the odd number of degenerate ground states for integer S, so any attempt to construct a sharp SZM in those models must fail.
  • The SZM (non-exact) version is recovered when both boundaries are set to the integrable spin-flip-invariant value, with commutator exponentially small in Hilbert-Schmidt norm but not spectral norm.
  • In the spin-1/2 case, the ESZM's action on energy eigenstates is exactly characterized by the Bethe-ansatz boundary strings: eigenvalue -1 precisely on states containing the right-boundary string.
  • Numerical edge-autocorrelator plateaus for the spin-1 chain require additional ESZMs beyond the single constructed Ψ, indicating a family of conserved edge operators.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the spectral-gap assertion holds, the construction likely extends by fusion to higher auxiliary spins, and the paper's conjectured S=3/2 ESZM from T^{(3/2,3/2)}(iπ/2) would give a conventional strongly localized mode; a direct check of its full MPO would settle that.
  • The 'weak locality but normalizable' pattern suggests a general principle: when the ground-state manifold has odd degeneracy, edge operators can still be exact conserved charges but must act nontrivially on an exponentially small (in Hilbert-Schmidt measure) fraction of the Hilbert space, so disorder-free prethermal-like plateaus survive.
  • The boundary-string characterization for S=1/2 suggests a dictionary between ESZM eigenvalues and root structure that, if extended to higher S, would let one compute exact edge-coherence plateaus directly from Bethe ansatz data without diagonalizing.
  • Perturbing away from integrability while keeping the U(1) symmetry, the plateau in edge autocorrelators appears robust at accessible sizes; a quantitative prediction is that the plateau height for large L is set by the sum of squares of overlaps with all ESZMs, which could be tested against exact numerics.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper constructs, for integrable open spin-S XXZ chains with a Z2-invariant left boundary and generic right boundary, an operator Ψ obtained from the derivative of the transfer matrix T^{(S,1/2)}(u) at u=iπ/2. The authors show that Ψ commutes exactly with the Hamiltonian, has a Hilbert-Schmidt norm equal to 1 up to exponentially small corrections, nearly squares to the identity in the same norm, and is weakly localized near the left edge. They argue that these properties imply finite late-time plateaus in edge autocorrelation functions, present numerical evidence for S=1/2 and S=1, and discuss connections to boundary bound states in the Bethe ansatz for S=1/2.

Significance. If the construction is correct, it provides exact strong zero modes for arbitrary spin S in an integrable setting, extending the earlier MPO/transfer-matrix construction for S=1/2. The exact commutation with the Hamiltonian is a strong structural result, and the operator-norm formulation of locality is a useful generalization. The paper also gives a clear physical interpretation of the degenerate ground states as a first-order transition line and presents reproducible numerics for edge autocorrelators. The weakening of locality for S≥1 and the explanation in terms of an odd number of ground states for integer S is an interesting conceptual contribution. However, the central claims rely on a spectral-gap assertion that is stated but not proved, and for S=1 the constructed operator does not by itself explain the numerically observed plateau heights.

major comments (3)
  1. [§4.2.1, Eqs. (91)–(96)] The normalizability and near-inversion properties of Ψ rest on the assertion, made after Eq. (91), that the ancillary transfer matrix U (and similarly V in Eq. (95)) has a unique eigenvalue equal to 1 and all other eigenvalues strictly smaller in modulus for all S and η∈R. The text says this was 'verified by direct calculation', but no derivation or explicit computation is given. This is load-bearing: Eqs. (93) and (87), and hence the plateau argument leading to Eq. (102), all follow from the exponential decay of subleading eigenvalues. If a subleading eigenvalue had modulus 1, the exponentially small corrections would become O(1) and the ESZM would lose its normalizability and inversion properties. Please provide a proof, or at least a complete and reproducible algebraic verification, including a bound on the second-largest eigenvalue.
  2. [§4.2.2, Fig. 11 and Eq. (103)] For S=1 the overlap of s^z_1 with the constructed Ψ is |c1|²=0.2637 at η=1, whereas the numerical plateau in Fig. 11 is substantially larger. The paper correctly concludes that 'there must be additional ESZM operators', but none are constructed. This is not a cosmetic gap: the abstract and conclusion claim that the constructed ESZMs imply the observed infinite coherence times, but for S=1 the single explicitly constructed operator cannot account for the plateau. The argument needs either an explicit additional conserved edge operator, or a more careful statement that the plateau demonstrates the existence of additional ESZMs whose construction is left open.
  3. [§4.4, Eq. (116)] The claimed S=3/2 ESZM Ψ_{3/2}=T^{(3/2,3/2)}(iπ/2) is a conjecture: only the first two terms of its large-η expansion are matched with the perturbative construction. Since the paper uses this to argue that a 'conventional' ESZM exists for S=3/2, this part is not established. This does not affect the main T^{(S,1/2)} construction, but the text should clearly separate conjecture from proof, especially in the Conclusions where it is used as evidence.
minor comments (4)
  1. [Throughout] There are several typographical issues: 'spin-Schains' in the title line, 'cf.section' in §3.1, 'withre' in §4.2, and missing spaces around parentheses. These should be corrected.
  2. [§5, Eqs. (126)–(127)] The conjecture that the ESZM eigenvalue is +1 without a right-boundary string and −1 with one is supported by numerical results for L≤12, N≤6 in Appendix D. This is reasonable evidence, but the statement should explicitly note the limited system sizes, as the text already does implicitly in Table 1.
  3. [§4.2.2, Eq. (102)] The lower bound in Eq. (102) assumes that the only conserved edge-localized operators are the ESZMs. For S=1 the paper itself shows there must be more, so the inequality is not tight. This should be stated directly where the equation appears, not only later.
  4. [References] The construction relies on the companion papers [39,45,74], two of which are preprints. If the journal allows, it would help readers if the key MPO identities from [45] were summarized in an appendix, since the present paper's normalizability proof depends on them.

Circularity Check

0 steps flagged

No significant circularity: the ESZM construction is explicit and algebraic; the main caveat is an unproved spectral-gap assertion, which is a completeness risk rather than a circular step.

full rationale

The central construction defines Ψ as a derivative of the commuting transfer matrix T^{(S,1/2)}(u) at u*=iπ/2, so exact commutation [H,Ψ]=0 follows from the standard transfer-matrix commutation (63), which the paper explicitly uses via the L-operator (66), K-matrices (68), and the key identity (72). The locality and normalizability results, Eqs. (93) and (87), do rely on the assertion in §4.2.1 that U and V have a unique eigenvalue of modulus 1 and all others strictly smaller; the paper states this was 'verified by direct calculation' but gives no derivation. That is an unproved computational assumption and a genuine load-bearing gap, but it is not circular: the spectrum of U is not an input that is later renamed as the prediction. The overlap coefficients c1 in Eqs. (102)-(103) are computed from the constructed operator, not fitted to the numerical plateaus; indeed for S=1 the authors note that the observed plateaus exceed |c1|^2 and therefore require additional ESZM operators, which is an acknowledged incompleteness rather than a circular explanation. The self-citations to Refs. [39,45,74] point to earlier treatments of the S=1/2 case and to the plateau/overlap argument, but the present paper rederives the needed MPO form, recursion relations and norm computations explicitly, so those citations are not load-bearing. The unproved spectral assertion and the S=1 completeness issue should be counted as correctness/completeness risks, not as circularity.

Axiom & Free-Parameter Ledger

0 free parameters · 7 axioms · 2 invented entities

The central construction is self-contained down to standard integrability, but several load-bearing statements are asserted or conjectural: the spectral gap of U, the existence of additional S=1 ESZMs, the S=3/2 ESZM, and the Bethe-state eigenvalue conjecture. No parameters are fitted to data.

axioms (7)
  • standard math U_q(sl2) quantum-group algebra (Eq. 67) underlies the L-operator and transfer-matrix commutation.
    Standard integrability framework; used to construct commuting family (63).
  • domain assumption Transfer matrices T^{(S,S')}(u) commute for all u,v (Eq. 63).
    Core of the construction; taken from Yang-Baxter integrability, not re-proven.
  • ad hoc to paper The ancillary transfer matrix U (and V) has a unique eigenvalue 1 and all others strictly smaller in modulus for η∈R.
    Stated as 'verified by direct calculation' after Eq. (91); no derivation. Load-bearing for normalizability and Ψ²≈1.
  • ad hoc to paper Additional ESZM operators exist for the integrable spin-1 chain with boundary conditions (59).
    Inferred in §4.2.2 from plateau values exceeding the overlap of the constructed Ψ; operators not constructed.
  • ad hoc to paper For S=3/2, Ψ_{3/2}=T^{(3/2,3/2)}(iπ/2) is an ESZM with good locality.
    Conjecture in §4.4; only first two strong-coupling terms match; full evaluation beyond scope.
  • domain assumption The integrable spin-S XXZ chain has 2S+1 degenerate ground states along η>0, forming a first-order transition line.
    Based on corner transfer matrices and scattering [53,62], plus explicit large-η states; not fully derived for all η.
  • ad hoc to paper For S=1/2, the ESZM eigenvalue is +1 on Bethe states without a right-boundary string and −1 with one (Eqs. 126-127).
    Stated as conjecture, supported by numerics up to L=12 and analytic large-L ground-state-sector calculation.
invented entities (2)
  • Additional ESZM operators for the S=1 integrable chain no independent evidence
    purpose: Account for the infinite-temperature autocorrelation plateaus, which exceed the overlap of the explicitly constructed Ψ.
    Postulated in §4.2.2 on the basis of numerical plateaus; no explicit construction or independent falsifiable handle is given.
  • Conventional ESZM for S=3/2 from T^{(3/2,3/2)}(iπ/2) no independent evidence
    purpose: Produce a zero mode with strong spatial locality, unlike the weak-locality Ψ from T^{(S,1/2)}.
    Conjectured in §4.4; only the first two terms in the large-η expansion are checked.

pith-pipeline@v1.3.0-alltime-deepseek · 38583 in / 18014 out tokens · 152384 ms · 2026-08-03T17:51:59.323237+00:00 · methodology

0 comments
read the original abstract

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.

Figures

Figures reproduced from arXiv: 2512.07742 by Eric Vernier, Fabian H.L. Essler, Paul Fendley.

Figure 1
Figure 1. Figure 1: Triple-well Landau-Ginzburg potential. One involved but powerful method for deriving this degeneracy in integrable models such as ours is to use the corner transfer matrix [51, 52]. When a spin chain is the limit of an integrable two-dimensional classical model, this approach allows the computation of one-point functions in the presence of various boundary conditions. From these one can infer the ground-st… view at source ↗
Figure 2
Figure 2. Figure 2: Ground state phase diagram of the integrable S = 1 chain perturbed by a single-ion anisotropy term (36). 2.3 Spin 3 2 and beyond We expect that the integrable XXZ chain describes a first-order transition line for all η > 0 for any S > 1 2 . Corner transfer matrices [53] and exact scattering calculations indicate that the ground-state degeneracy becomes 2S + 1-fold [62]. Some intuition into this result come… view at source ↗
Figure 3
Figure 3. Figure 3: Landau-Ginzburg potential describing the four ground states for S = 3/2. Low-lying excitations consist of domain walls between between neighbouring minima. when not symmetry related, i.e. the gap for the wall between ground states 2 and 3 is the same as that for the wall between 1 and 2 (and 3 and 4). Note as well that the energies of other nearest-neighbor pairs from (41) also can be read off from the pic… view at source ↗
Figure 4
Figure 4. Figure 4: Ground state phase diagram of the model (36) for S = 3/2. We thus arrive at the ground state phase diagram shown in [PITH_FULL_IMAGE:figures/full_fig_p014_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Projected infinite-temperature autocorrelation function f0(t) for S = 1 2 , ∆ = 4.5, hL = 0 and system sizes L = 10, 12, 14. Here C n m denotes a binomial coefficient. Even when breaking this pairing by taking hL ̸= 0, the projected autocorrelator f0(t) dominates the infinite-temperature autocorrelator. We start with free boundary case hL = 0, which possesses a SZM [23]. Results for the projected autocorre… view at source ↗
Figure 6
Figure 6. Figure 6: Left panel: Projected infinite-temperature autocorrelation function f0(t) for S = 1 2 , ∆ = 3.5, system size L = 10 and several values of hL. At late times f0(t) relaxes to a finite constant due to the presence of an ESZM. Right panel: same for hL = 0.1 and three different system sizes. The dashed red line is the prediction from the ESZM projection (102) for L ≫ 1. left panel shows that as hL becomes small… view at source ↗
Figure 7
Figure 7. Figure 7: Projected infinite-temperature autocorrelation function f0(t) for S = 1 2 , ∆ = 4.5, hL = 0, J x 2 = 0.25 and system sizes L = 10, 12, 14. We see the emergence of plateau at intermediate times, the duration of which grows with L. feature is the plateau at intermediate times, the duration of which grows with L. At very late times the projected autocorrelator decays to a very small value. Both features are i… view at source ↗
Figure 8
Figure 8. Figure 8: Projected infinite-temperature autocorrelation function f0(t) for S = 1 2 , ∆ = 4.5, hL = 0.1, J x 2 = 0.25 and system sizes L = 10, 12, 14. with L. Moreover, at large times, the projected autocorrelator now does not decay to 17 [PITH_FULL_IMAGE:figures/full_fig_p017_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Infinite-temperature edge autocorrelation function C s z 1 0 (t) for the non-integrable, open S = 1 XXZ chain (57) with ∆ = 3.7622 and system sizes L = 6, 8, 10. The solid lines are the infinite time averages (58). The inset shows early times on a double logarithmic scale. 3.4 Edge coherence in the integrable spin-1 chain We now turn to the integrable spin-1 chain (2), (17). We first consider open boundary… view at source ↗
Figure 10
Figure 10. Figure 10: Infinite-temperature edge autocorrelation function C s z 1 (t) for the S = 1 chain with Hamiltonian density (17) with η = 1, vanishing boundary fields b (1) 1 = 0 = b (1) L and system sizes L = 6, 8. The solid lines are the infinite-time averages (58). 3.4.1 Integrable, spin-flip invariant boundary conditions I We now fix the boundary terms in the integrable spin-1 Hamiltonian to be b (1) 1 = [PITH_FULL_… view at source ↗
Figure 11
Figure 11. Figure 11: Infinite-temperature autocorrelation functions (50) for O1 = s z 1 (upper two curves) and O2 = (1 − (s z 1 ) 2 )s z 2 (lower two curves) in the integrable open spin-1 chain with η = 1, boundary potentials (59), and system sizes 6 ≤ L ≤ 9. These results suggest the presence of an ESZM operator localized around the left boundary. autocorrelation functions in [PITH_FULL_IMAGE:figures/full_fig_p019_11.png] view at source ↗
Figure 5
Figure 5. Figure 5 [PITH_FULL_IMAGE:figures/full_fig_p020_5.png] view at source ↗
Figure 12
Figure 12. Figure 12: Infinite-temperature autocorrelation functions (50) for O1 = s z 1 in the integrable open spin-1 chain with η = 1, boundary potentials (60), and system sizes 6 ≤ L ≤ 9. These results suggest the presence of a SZM operator localized around the left boundary. 3.4.2 Integrable, spin-flip invariant boundary conditions II Spin-flip invariant boundary conditions are also obtained by flipping the two boundaries … view at source ↗
Figure 13
Figure 13. Figure 13: Infinite-temperature autocorrelation functions (50) for O1 = s z 1 (upper two curves) and O2 = (1−(s z 1 ) 2 )s z 2 (lower two curves) in the integrable open spin-1 chain with η = 1, boundary potentials (61), and system sizes L = 6 and L = 8. In the parametrization used in section 4 these correspond to ξ + = 0 = ξ −. The infinite￾temperature autocorrelation function C s z 1 0 (t) for these boundary condit… view at source ↗
Figure 15
Figure 15. Figure 15: Infinite-temperature edge autocorrelation function C O1,2 (t) for O1 = s z 1 (upper two curves) and O2 = (1 − (s z 1 ) 2 )s z 2 in the integrable spin-1 chain with η = 1, boundary potentials (59) perturbed by (36) with D = 0.25 and system sizes 6 ≤ L ≤ 9. 4 Exact strong zero modes from transfer matrices We explained in section (2.2) that the structure of SZMs in the integrable spin-1 chain cannot be as el… view at source ↗
Figure 16
Figure 16. Figure 16: The transfer matrix T (S,S′ ) (u) acting on L sites. The physical spin S are represented in red, and the auxiliary spin S ′ in blue. the L-operators, and the K± S′(u) are reflection matrices [68, 69] that act on the auxiliary space and encode the boundary conditions. The transfer matrices thus can be written as the trace over the 2S ′ + 1-dimensional auxiliary space, namely T (S,S′ ) (u) = tra  K+ a (u)T… view at source ↗
Figure 17
Figure 17. Figure 17: Hibert-Schmidt norm of Ψj from (85), for S = 1 2 (left panel) and S = 1 (right panel). left boundary. While the latter action for j large is non-trivial on an exponentially large number of states, their number is still negligible compared to the full Hilbert space. Hence this behaviour is still consistent with the vanishing norms derived above. This lack of exact localisation for S > 1 2 can be quantified… view at source ↗
Figure 18
Figure 18. Figure 18: Spectral norm [PITH_FULL_IMAGE:figures/full_fig_p030_18.png] view at source ↗
Figure 10
Figure 10. Figure 10: 4.4 ESZM operators from transfer matrices T (S,S′> 1 2 ) (u) The ESZM described above is constructed from the transfer matrix T (S,S′ ) (u), with auxil￾iary spin S ′ = 1 2 and physical spin S. A natural question is whether one can obtain other ESZMs from other values of S ′ . Our work indicates that only for some values of (S, S′ ) it is possible. For the simplest case (S, S′ ) = ( 1 2 ,1) the answer appe… view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 3 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Exact strong zero modes are generic in integrable spin systems with large anisotropy

    quant-ph 2026-05 unverdicted novelty 7.0

    Exact strong zero modes arise generically in integrable spin systems with large anisotropy from quasi-periodicity of the R-matrix and tracelessness of the K-matrix.

  2. Exact strong zero modes are generic in integrable spin systems with large anisotropy

    quant-ph 2026-05 unverdicted novelty 7.0

    Exact strong zero modes arise generically in integrable anisotropic spin models from quasi-periodicity of R-matrices and tracelessness of K-matrices, unifying known cases and predicting new ones.

  3. Quasi-local Edge Mode in XXX Spin Chain/Circuit with Interaction Boundary Defect

    cond-mat.stat-mech 2026-03 conditional novelty 7.0

    Above a critical boundary interaction, the XXX spin chain/circuit hosts an explicit matrix-product quasi-local edge mode that yields a nonzero boundary Drude weight and a transition to ergodic boundary dynamics below ...

Reference graph

Works this paper leans on

86 extracted references · 12 canonical work pages · cited by 2 Pith papers

  1. [1]

    A. B. Zamolodchikov and V. Fateev,Model factorized s-matrix and an integrable spin-1 heisenberg chain, Sov. J. Nucl. Phys.(Engl. Transl.);(United States)32(2) (1980)

  2. [2]

    P. P. Kulish, N. Y. Reshetikhin and E. K. Sklyanin,Yang-baxter equation and repre- sentation theory: I, Letters in Mathematical Physics5, 393 (1981)

  3. [3]

    Kulish and N

    P. Kulish and N. Y. Reshetikhin,Generalized heisenberg ferromagnet and the gross- neveu model, Sov. Phys. JETP53(1), 108 (1981)

  4. [4]

    Takhtajan,The picture of low-lying excitations in the isotropic heisenberg chain of arbitrary spins, Physics Letters A87(9), 479 (1982)

    L. Takhtajan,The picture of low-lying excitations in the isotropic heisenberg chain of arbitrary spins, Physics Letters A87(9), 479 (1982)

  5. [5]

    Babujian,Exact solution of the one-dimensional isotropic heisenberg chain with arbitrary spins s, Physics Letters A90(9), 479 (1982)

    H. Babujian,Exact solution of the one-dimensional isotropic heisenberg chain with arbitrary spins s, Physics Letters A90(9), 479 (1982)

  6. [6]

    Babujian,Exact solution of the isotropic heisenberg chain with arbitrary spins: thermodynamics of the model, Nuclear Physics B215(3), 317 (1983)

    H. Babujian,Exact solution of the isotropic heisenberg chain with arbitrary spins: thermodynamics of the model, Nuclear Physics B215(3), 317 (1983)

  7. [7]

    Sogo,Ground state and low-lying excitations in the heisenberg xxz chain of arbi- trary spin s, Physics Letters A104(1), 51 (1984)

    K. Sogo,Ground state and low-lying excitations in the heisenberg xxz chain of arbi- trary spin s, Physics Letters A104(1), 51 (1984)

  8. [8]

    A. N. Kirillov and N. Y. Reshetikhin,Exact solution of the heisenberg xxz model of spin s, Journal of Soviet Mathematics35, 2627 (1986)

  9. [9]

    A. N. Kirillov and N. Y. Reshetikhin,Exact solution of the integrable xxz heisenberg model with arbitrary spin. i. the ground state and the excitation spectrum, Journal of Physics A: Mathematical and General20(6), 1565 (1987)

  10. [10]

    Kirillov and N

    A. Kirillov and N. Y. Reshetikhin,Exact solution of the integrable xxz heisenberg model with arbitrary spin. ii. thermodynamics of the system, Journal of Physics A: Mathematical and General20(6), 1587 (1987)

  11. [11]

    Frahm, N.-C

    H. Frahm, N.-C. Yu and M. Fowler,The integrable xxz heisenberg model with arbitrary spin: construction of the hamiltonian, the ground-state configuration and conformal properties, Nuclear Physics B336(3), 396 (1990)

  12. [12]

    Suzuki,Excited state nonlinear integral equations for an integrable anisotropic spin-1 chain, Journal of Physics A: Mathematical and General37(50), 11957 (2004)

    J. Suzuki,Excited state nonlinear integral equations for an integrable anisotropic spin-1 chain, Journal of Physics A: Mathematical and General37(50), 11957 (2004)

  13. [13]

    Deguchi and C

    T. Deguchi and C. Matsui,Correlation functions of the integrable higher-spin xxx and xxz spin chains through the fusion method, Nuclear Physics B831(3), 359 (2010)

  14. [14]

    Mezincescu, R

    L. Mezincescu, R. I. Nepomechie and V. Rittenberg,Bethe ansatz solution of the fateev-zamolodchikov quantum spin chain with boundary terms, Physics Letters A 147(1), 70 (1990)

  15. [15]

    De Vega and A

    H. De Vega and A. Gonzalez-Ruiz,Exact solution of the suq (n)-invariant quantum spin chains, Nuclear Physics B417(3), 553 (1994)

  16. [16]

    Affleck and F

    I. Affleck and F. Haldane,Critical theory of quantum spin chains, Physical Review B36(10), 5291 (1987). 44 SciPost Physics Submission

  17. [17]

    L¨ auchli, G

    A. L¨ auchli, G. Schmid and S. Trebst,Spin nematics correlations in bilinear-biquadratic spin chains, Physical Review B74(14) (2006), doi:10.1103/physrevb.74.144426

  18. [18]

    A. Y. Kitaev,Unpaired majorana fermions in quantum wires, Physics-uspekhi 44(10S), 131 (2001)

  19. [19]

    Gu and X.-G

    Z.-C. Gu and X.-G. Wen,Tensor-entanglement-filtering renormalization approach and symmetry-protected topological order, Physical Review B80(15) (2009), doi:10.1103/physrevb.80.155131

  20. [20]

    Pollmann, E

    F. Pollmann, E. Berg, A. M. Turner and M. Oshikawa,Symmetry protection of topological phases in one-dimensional quantum spin systems, Physical Review B 85(7) (2012), doi:10.1103/physrevb.85.075125

  21. [21]

    Fendley,Parafermionic edge zero modes in zn-invariant spin chains, Journal of Statistical Mechanics: Theory and Experiment2012(11), P11020 (2012)

    P. Fendley,Parafermionic edge zero modes in zn-invariant spin chains, Journal of Statistical Mechanics: Theory and Experiment2012(11), P11020 (2012)

  22. [22]

    Alicea and P

    J. Alicea and P. Fendley,Topological phases with parafermions: theory and blueprints, Annual Review of Condensed Matter Physics7(1), 119 (2016)

  23. [23]

    Fendley,Strong zero modes and eigenstate phase transitions in the xyz/interacting majorana chain, Journal of Physics A: Mathematical and Theoretical49(30), 30LT01 (2016)

    P. Fendley,Strong zero modes and eigenstate phase transitions in the xyz/interacting majorana chain, Journal of Physics A: Mathematical and Theoretical49(30), 30LT01 (2016)

  24. [24]

    Moran, D

    N. Moran, D. Pellegrino, J. Slingerland and G. Kells,Parafermionic clock models and quantum resonance, Physical Review B95(23), 235127 (2017)

  25. [25]

    L. M. Vasiloiu, F. Carollo, M. Marcuzzi and J. P. Garrahan,Strong zero modes in a class of generalized ising spin ladders with plaquette interactions, Physical Review B 100(2), 024309 (2019)

  26. [26]

    J. Kemp, N. Y. Yao and C. R. Laumann,Symmetry-enhanced boundary qubits at infinite temperature, Physical Review Letters125(20), 200506 (2020)

  27. [27]

    J. Kemp, N. Y. Yao, C. R. Laumann and P. Fendley,Long coherence times for edge spins, Journal of Statistical Mechanics: Theory and Experiment2017(6), 063105 (2017)

  28. [28]

    D. V. Else, P. Fendley, J. Kemp and C. Nayak,Prethermal strong zero modes and topological qubits, Physical Review X7(4), 041062 (2017)

  29. [29]

    Jinet al.,Topological prethermal strong zero modes on superconducting processors, Nature645(8081), 626 (2025), doi:10.1038/s41586-025-09476-z,2501.04688

    F. Jinet al.,Topological prethermal strong zero modes on superconducting processors, Nature645(8081), 626 (2025), doi:10.1038/s41586-025-09476-z,2501.04688

  30. [30]

    Thakurathi, A

    M. Thakurathi, A. A. Patel, D. Sen and A. Dutta,Floquet generation of ma- jorana end modes and topological invariants, Phys. Rev. B88, 155133 (2013), doi:10.1103/PhysRevB.88.155133

  31. [31]

    Chandran, V

    A. Chandran, V. Khemani, C. R. Laumann and S. L. Sondhi,Many-body local- ization and symmetry-protected topological order, Phys. Rev. B89, 144201 (2014), doi:10.1103/PhysRevB.89.144201

  32. [32]

    Bahri, R

    Y. Bahri, R. Ronen and E. Altman,Localization and topology protected quantum coherence at the edge of hot matter, Nature Communications6, 7341 (2015). 45 SciPost Physics Submission

  33. [33]

    Iadecola, L

    T. Iadecola, L. H. Santos and C. Chamon,Stroboscopic symmetry-protected topological phases, Physical Review B92(12), 125107 (2015)

  34. [34]

    Khemani, A

    V. Khemani, A. Lazarides, R. Moessner and S. L. Sondhi,Phase struc- ture of driven quantum systems, Phys. Rev. Lett.116, 250401 (2016), doi:10.1103/PhysRevLett.116.250401

  35. [35]

    Sreejith, A

    G. Sreejith, A. Lazarides and R. Moessner,Parafermion chain with 2π/k floquet edge modes, Physical Review B94(4), 045127 (2016)

  36. [36]

    Potirniche, A

    I.-D. Potirniche, A. C. Potter, M. Schleier-Smith, A. Vishwanath and N. Y. Yao, Floquet symmetry-protected topological phases in cold-atom systems, Phys. Rev. Lett. 119, 123601 (2017), doi:10.1103/PhysRevLett.119.123601

  37. [37]

    Kumar, P

    A. Kumar, P. T. Dumitrescu and A. C. Potter,String order parameters for one- dimensional floquet symmetry protected topological phases, Phys. Rev. B97, 224302 (2018), doi:10.1103/PhysRevB.97.224302

  38. [38]

    Mukherjee, R

    B. Mukherjee, R. Melendrez, M. Szyniszewski, H. J. Changlani and A. Pal,Emergent strong zero mode through local floquet engineering, Physical Review B109(6), 064303 (2024)

  39. [39]

    Vernier, H.-C

    E. Vernier, H.-C. Yeh, L. Piroli and A. Mitra,Strong zero modes in integrable quantum circuits, Phys. Rev. Lett.133, 050606 (2024), doi:10.1103/PhysRevLett.133.050606

  40. [40]

    D. J. Yates, F. H. L. Essler and A. Mitra,Almost strong (0,π) edge modes in clean interacting one-dimensional floquet systems, Physical Review B99(20), 205419 (2019)

  41. [41]

    D. J. Yates, A. G. Abanov and A. Mitra,Lifetime of almost strong edge-mode opera- tors in one-dimensional, interacting, symmetry protected topological phases, Physical Review Letters124(20), 206803 (2020)

  42. [42]

    Klobas, P

    K. Klobas, P. Fendley and J. P. Garrahan,Stochastic strong zero modes and their dynamical manifestations, Physical Review E107(4), L042104 (2023)

  43. [43]

    C. T. Olund, N. Y. Yao and J. Kemp,Boundary strong zero modes, arXiv preprint arXiv:2305.16382 (2023)

  44. [44]

    R. J. Baxter,Exactly solved models in statistical mechanics, doi:https://physics.anu.edu.au/research/ftp/mpg/baxter book.php (1982)

  45. [45]

    Fendley, S

    P. Fendley, S. Gehrmann, E. Vernier and F. Verstraete,XYZ integrability the easy way(2025),2511.04674

  46. [46]

    Hagendorf,Spin chains with dynamical lattice supersymmetry, Journal of Statis- tical Physics150(4), 609 (2013)

    C. Hagendorf,Spin chains with dynamical lattice supersymmetry, Journal of Statis- tical Physics150(4), 609 (2013)

  47. [47]

    Yang and P

    X. Yang and P. Fendley,Non-local spacetime supersymmetry on the lattice, Journal of Physics A: Mathematical and General37(38), 8937 (2004)

  48. [48]

    Fendley, K

    P. Fendley, K. Sengupta and S. Sachdev,Competing density-wave orders in a one-dimensional hard-boson model, Phys. Rev. B69(7), 075106 (2004), doi:10.1103/PhysRevB.69.075106

  49. [49]

    A. B. Zamolodchikov,Thermodynamic Bethe ansatz for RSOS scattering theories, Nucl. Phys. B358, 497 (1991), doi:10.1016/0550-3213(91)90422-T. 46 SciPost Physics Submission

  50. [50]

    O’Brien and P

    E. O’Brien and P. Fendley,Lattice supersymmetry and order-disorder coexis- tence in the tricritical Ising model, Phys. Rev. Lett.120(20), 206403 (2018), doi:10.1103/PhysRevLett.120.206403,1712.06662

  51. [51]

    G. E. Andrews, R. J. Baxter and P. J. Forrester,Eight vertex SOS model and generalized Rogers-Ramanujan type identities, J. Statist. Phys.35, 193 (1984), doi:10.1007/BF01014383

  52. [52]

    D. A. Huse,Exact exponents for infinitely many new multicritical points, Phys. Rev. B30, 3908 (1984), doi:10.1103/PhysRevB.30.3908

  53. [53]

    E. Date, M. Jimbo, A. Kuniba, T. Miwa and M. Okado,One-dimensional Configu- ration Sums in Vertex Models and Affine Lie Algebra Characters, Lett. Math. Phys. 17, 69 (1989), doi:10.1007/BF00420017

  54. [54]

    Di Vecchia and S

    P. Di Vecchia and S. Ferrara,Classical solutions in two-dimensional supersymmetric field theories, Nuclear Physics B130(1), 93 (1977)

  55. [55]

    Ahn,Complete s-matrices of supersymmetric sine-gordon theory and perturbed superconformal minimal model, Nuclear Physics B354(1), 57 (1991)

    C. Ahn,Complete s-matrices of supersymmetric sine-gordon theory and perturbed superconformal minimal model, Nuclear Physics B354(1), 57 (1991)

  56. [56]

    C. Ahn, R. I. Nepomechie and J. Suzuki,Finite size effects in the spin-1 xxz and su- persymmetric sine-gordon models with dirichlet boundary conditions, Nuclear Physics B767(3), 250 (2007)

  57. [57]

    Tsvelik,Field-theory treatment of the heisenberg spin-1 chain, Physical Review B 42(16), 10499 (1990)

    A. Tsvelik,Field-theory treatment of the heisenberg spin-1 chain, Physical Review B 42(16), 10499 (1990)

  58. [58]

    F. H. L. Essler and I. Affleck,Haldane-gap chains in a magnetic field, Journal of Statistical Mechanics: Theory and Experiment2004(12), P12006 (2004)

  59. [59]

    A. O. Gogolin, A. A. Nersesyan and A. M. Tsvelik,Bosonization and strongly corre- lated systems, Cambridge university press (2004)

  60. [60]

    Reshetikhin and F

    N. Reshetikhin and F. Smirnov,Hidden Quantum Group Symmetry and Integrable Perturbations of Conformal Field Theories, Commun. Math. Phys.131, 157 (1990), doi:10.1007/BF02097683

  61. [61]

    P. H. Ginsparg,APPLIED CONFORMAL FIELD THEORY, InLes Houches Summer School in Theoretical Physics: Fields, Strings, Critical Phenomena(1988), hep-th/9108028

  62. [62]

    Bernard and A

    D. Bernard and A. Leclair,The Fractional supersymmetric Sine-Gordon models, Phys. Lett. B247, 309 (1990), doi:10.1016/0370-2693(90)90901-H

  63. [63]

    Dupont and J

    M. Dupont and J. E. Moore,Universal spin dynamics in infinite-temperature one- dimensional quantum magnets, Physical Review B101(12), 121106 (2020)

  64. [64]

    Y. kui Zhou,Fusion hierarchy and finite-size corrections of uq[sl(2)]-invariant vertex models with open boundaries, Nuclear Physics B453(3), 619 (1995), doi:https://doi.org/10.1016/0550-3213(95)00293-2

  65. [65]

    Mezincescu and R

    L. Mezincescu and R. I. Nepomechie,Fusion procedure for open chains, Journal of Physics A: Mathematical and General25(9), 2533 (1992), doi:10.1088/0305- 4470/25/9/024. 47 SciPost Physics Submission

  66. [66]

    Yu-kui Zhou,Row transfer matrix functional relations for Baxter’s eight-vertex and six-vertex models with open boundaries via more general reflection matrices, Nuclear Physics B458(3), 504 (1996), doi:https://doi.org/10.1016/0550-3213(95)00553-6

  67. [67]

    E. K. Sklyanin,Boundary conditions for integrable quantum systems, Journal of Physics A: Mathematical and General21(10), 2375 (1988)

  68. [68]

    I. V. Cherednik,Factorizing particles on a half-line and root systems, Teoreticheskaya i Matematicheskaya Fizika61(1), 35 (1984)

  69. [69]

    De Vega and A

    H. De Vega and A. Gonz´ alez-Ruiz,Boundary k-matrices for the xyz, xxz and xxx spin chains, Journal of Physics A: Mathematical and General27(18), 6129 (1994)

  70. [70]

    Piroli, B

    L. Piroli, B. Pozsgay and E. Vernier,From the quantum transfer matrix to the quench action: the loschmidt echo in xxz heisenberg spin chains, Journal of Statisti- cal Mechanics: Theory and Experiment2017(2), 023106 (2017), doi:10.1088/1742- 5468/aa5d1e

  71. [71]

    Ilievski, J

    E. Ilievski, J. De Nardis, B. Wouters, J.-S. Caux, F. H. Essler and T. Prosen,Complete generalized gibbs ensembles in an interacting theory, Physical review letters115(15), 157201 (2015)

  72. [72]

    Ilievski, M

    E. Ilievski, M. Medenjak, T. Prosen and L. Zadnik,Quasilocal charges in integrable lattice systems, Journal of Statistical Mechanics: Theory and Experiment2016(6), 064008 (2016)

  73. [73]

    Saleur and J.-B

    H. Saleur and J.-B. Zuber,Integrable lattice models and quantum groups, Spring School on String Theory and Quantum Gravity Trieste, Italy23, 0001 (1990)

  74. [74]

    Gehrmann and F

    S. Gehrmann and F. H. L. Essler,Exact strong zero modes in quantum circuits and spin chains with non-diagonal boundary conditions, arXiv preprint arXiv:2511.05490 (2025)

  75. [75]

    Kapustin and S

    A. Kapustin and S. Skorik,Surface excitations and surface energy of the antiferro- magnetic xxz chain by the bethe ansatz approach, Journal of Physics A: Mathematical and General29(8), 1629 (1996)

  76. [76]

    Asakawa and M

    H. Asakawa and M. Suzuki,Finite-size corrections in the xxz model and the hubbard model with boundary fields, Journal of Physics A: Mathematical and General29(2), 225 (1996)

  77. [77]

    F. H. Eßler and H. Frahm,X-ray edge singularity in integrable lattice models of correlated electrons, Physical Review B56(11), 6631 (1997)

  78. [78]

    Bed¨ urftig and H

    G. Bed¨ urftig and H. Frahm,Spectrum of boundary states in the open hubbard chain, Journal of Physics A: Mathematical and General30(12), 4139 (1997)

  79. [79]

    Deguchi, R

    T. Deguchi, R. Yue and K. Kusakabe,A gapless charge mode induced by the boundary states in the half-filled hubbard open chain, Journal of Physics A: Mathematical and General31(36), 7315 (1998)

  80. [80]

    Grijalva, J

    S. Grijalva, J. D. Nardis and V. Terras,Open XXZ chain and boundary modes at zero temperature, SciPost Phys.7, 023 (2019), doi:10.21468/SciPostPhys.7.2.023

Showing first 80 references.