The open spin-1/2 XXZ chain in the massive antiferromagnetic regime (Δ>1\Delta>1) with longitudinal boundary fields supports boundary modes that appear in the Bethe-ansatz solution as isolated complex roots, called boundary roots, exponentially localized at an edge. Grijalva, De Nardis and Terras (SciPost Phys. 7, 023 (2019)) showed that these roots drive a ground-state double degeneracy for equal fields ∣h±∣<Δ−1|h_\pm|<\Delta-1 and explicitly posed the questions restated here: is the operatorial strong zero mode the creation/annihilation operator of the boundary-root excitation, and can it be constructed within the Bethe-ansatz/integrability framework; do representative finite-temperature Bethe/TBA states contain boundary roots, with two thermodynamically distinct thermal states distinguished by opposite boundary-root deviations; how should boundary-root data be incorporated into generalized Gibbs ensembles describing post-quench steady states near the edge, given that bulk-extensive conserved charges do not encode them; and whether supersymmetric structures connect to boundary roots. Recent work (2024-2026) constructed exact strong zero modes from commuting transfer matrices and partially identified their action on Bethe eigenstates, but the identification is restricted to one vanishing boundary field, is proved only in the ground-state sector, and is unproven at general energy densities; the finite-temperature, non-equilibrium and supersymmetric questions are untouched. The problem, namely the operatorial and statistical-mechanical fate of the boundary root beyond the ground state, remains open.

Background

The spin-1/2 XXZ chain is the paradigmatic interacting one-dimensional quantum magnet. In this problem it is taken on an open (finite) chain of LL sites with longitudinal magnetic fields applied only at the two ends,

H=∑j=1L−1[σjxσj+1x+σjyσj+1y+Δ(σjzσj+1z−1)]+h−σ1z+h+σLz,H = \sum_{j=1}^{L-1}\left[\sigma^x_j\sigma^x_{j+1}+\sigma^y_j\sigma^y_{j+1}+\Delta\left(\sigma^z_j\sigma^z_{j+1}-1\right)\right] + h_-\sigma^z_1 + h_+\sigma^z_L,

where σjx,y,z\sigma^{x,y,z}_j are Pauli matrices acting on site jj, Δ>1\Delta>1 (conventionally Δ=cosh⁡η\Delta=\cosh\eta) is the anisotropy, and h±h_\pm are the boundary fields. For Δ>1\Delta>1 the chain is in the massive antiferromagnetic regime: the bulk spectrum is gapped and the ground state has antiferromagnetic Néel order. Open integrable chains remain exactly solvable: eigenstates in the sector with NN down spins are parametrized by rapidities obeying Sklyanin's boundary Bethe equations, the boundary-integrable form of the algebraic Bethe ansatz. Besides real roots and bound-state strings (patterns of nearly parallel rapidities), in this regime the boundary phase factors possess zeros that can attract isolated complex roots, called boundary roots, which converge to a zero of the boundary factor exponentially fast in the chain length LL and describe an excitation exponentially localized at one edge, i.e. a boundary bound state of magnons. The source paper (Grijalva, De Nardis and Terras, 2019) showed that for even LL and ∣h±∣<hcr(1)=Δ−1|h_\pm|<h_{\rm cr}^{(1)}=\Delta-1 the boundary root is present in the Bethe solutions of both the ground state and the first excited state; for equal fields h−=h+=hh_-=h_+=h the spectrum is gapped and the ground state is doubly degenerate up to corrections exponentially small in LL, the two states being distinguished by opposite deviations of the boundary root, and the boundary root itself being delocalized between the two edges. As consequences, the magnetization at one edge depends on the value of the field at the opposite edge even in the half-infinite-chain limit, with a discontinuity at h−=h+h_-=h_+, and the zero-temperature spin autocorrelation at the edge saturates at long times to a finite value fixed by the boundary-root contribution. A strong zero mode is an operator localized near one edge that commutes with the Hamiltonian up to corrections exponentially small in LL and anticommutes with a discrete symmetry; it maps every eigenstate in one symmetry sector to a nearly degenerate partner in the conjugate sector, so that the near-degeneracy persists at all energy densities rather than only in the ground state, and it endows edge spins with anomalously long coherence times (its breakdown at a critical coupling defines an eigenstate phase transition). For this model class a strong zero mode was constructed by Fendley (2016) for the massive XYZ chain, which contains the free-boundary XXZ chain as a special case; it was later recast as a series in 1/Δ1/\Delta and, recently, in exact transfer-matrix form. The Bethe-ansatz picture and the operatorial picture describe the same edge physics from two sides, yet the boundary root is a property of individual eigenstates while the strong zero mode is an operator on the whole Hilbert space, and the source paper left their exact relation open. Three statistical-mechanical questions complete the problem. First, finite temperature: the thermodynamic Bethe ansatz (TBA), introduced by Takahashi, rewrites the finite-temperature thermodynamics of a Bethe-solvable model as integral equations for rapidity densities, with representative states characterizing a macrostate; it should say whether thermal representative states contain boundary roots, and whether the strong-zero-mode degeneracy at all energies translates into two thermodynamically distinct representative thermal states distinguished by opposite boundary-root deviations, in analogy with the ground-state pair. The source also predicts a non-decaying connected edge autocorrelation at any temperature, whose verification within TBA requires precisely this input. Second, non-equilibrium: after a quantum quench (a sudden change of a Hamiltonian parameter), integrable systems relax to a generalized Gibbs ensemble (GGE), the statistical ensemble built from all extensive conserved charges; since these charges do not encode boundary-root data, and Bertini and Fagotti showed that steady states emerging from a defect require information beyond the bulk charges, the way boundary-root information enters the post-quench steady state near the edge is open. Third, structure: the open XXZ chain admits lattice supersymmetry in other parameter regimes (Hagendorf and Liénardy at Δ=−1/2\Delta=-1/2; Weston and Yang via the algebraic Bethe ansatz), and whether supersymmetric structures connect to boundary roots in the gapped regime was suggested but not resolved.

Problem Statement

The questions are posed for the open spin-1/2 XXZ chain with Hamiltonian as above, in the massive antiferromagnetic regime Δ>1\Delta>1, for even chain length LL and boundary fields in the boundary-root regime ∣h±∣<Δ−1|h_\pm|<\Delta-1. (1) Operatorial identity of the boundary root. Determine whether the strong zero mode is the creation/annihilation operator of the boundary excitation corresponding to the boundary root at an edge of the chain, and construct this operator within the Bethe-ansatz/integrability framework. Existing exact constructions impose one vanishing boundary field; the question remains open for generic boundary fields, in particular the equal-field case h−=h+h_-=h_+ (including h=0h=0, where Fendley's original strong zero mode lives), in which the boundary root is delocalized between the two edges and the ground state is quasi-degenerate. (2) Finite-temperature fate. Determine whether representative finite-temperature Bethe/TBA states of the open chain contain boundary roots, and whether there exist two thermodynamically distinct representative thermal states distinguished by two opposite deviations of the boundary root, analogous to the ground-state pair at h−=h+h_-=h_+; decide accordingly whether the strong-zero-mode-induced degeneracy at all energies translates into a two-fold structure of the thermal description, and give the corresponding prediction for the long-time limit of the connected edge autocorrelation at finite temperature. (3) Non-equilibrium bookkeeping. Formulate how boundary-root information is incorporated into the ensemble description of post-quench steady states near the edge, for example as a generalized Gibbs ensemble supplemented by boundary data beyond the bulk-extensive conserved charges, and characterize the resulting edge observables in the steady state. (4) Related structural question. Decide whether the open XXZ chain with boundary roots admits supersymmetric properties connecting the boundary modes, in analogy with the lattice supersymmetry found in other regimes of the open XXZ chain. Items (1)-(3) form one scientific objective, and settling that objective resolves the task: to promote the boundary root from a zero-temperature, ground-state diagnostic of the Bethe equations to a fully identified operatorial and statistical degree of freedom of the gapped open chain. Item (4) is a structural probe of the same boundary modes, explicitly suggested by the source work; it may be resolved in addition, but a resolution of the task does not require it. Each item may be resolved by explicit construction, by proof, or by a controlled demonstration that the putative structure is absent (for example, a proof that the thermal description admits only one representative state).

References

  1. S. Grijalva, J. De Nardis, V. Terras, Open XXZ chain and boundary modes at zero temperature, SciPost Phys. 7, 023 (2019), DOI 10.21468/SciPostPhys.7.2.023, arXiv:1901.10932, https://arxiv.org/abs/1901.10932
  2. P. Fendley, Strong zero modes and eigenstate phase transitions in the XYZ/interacting Majorana chain, J. Phys. A: Math. Theor. 49, 30LT01 (2016), DOI 10.1088/1751-8113/49/30/30LT01, arXiv:1512.03441, https://arxiv.org/abs/1512.03441
  3. E. K. Sklyanin, Boundary conditions for integrable quantum systems, J. Phys. A: Math. Gen. 21, 2375 (1988), DOI 10.1088/0305-4470/21/10/015, https://doi.org/10.1088/0305-4470/21/10/015
  4. M. Takahashi, One-Dimensional Heisenberg Model at Finite Temperature, Prog. Theor. Phys. 46, 401 (1971), DOI 10.1143/PTP.46.401, https://doi.org/10.1143/PTP.46.401
  5. B. Bertini and M. Fagotti, Determination of the Nonequilibrium Steady State Emerging from a Defect, Phys. Rev. Lett. 117, 130402 (2016), DOI 10.1103/PhysRevLett.117.130402, https://doi.org/10.1103/PhysRevLett.117.130402
  6. E. Vernier, H.-C. Yeh, L. Piroli, A. Mitra, Strong Zero Modes in Integrable Quantum Circuits, Phys. Rev. Lett. 133, 050606 (2024), DOI 10.1103/PhysRevLett.133.050606, arXiv:2401.12305, https://arxiv.org/abs/2401.12305
  7. P. Fendley, S. Gehrmann, E. Vernier, F. Verstraete, XYZ integrability the easy way, J. Stat. Phys. 193, 79 (2026), DOI 10.1007/s10955-026-03645-8, arXiv:2511.04674, https://arxiv.org/abs/2511.04674
  8. F. H. L. Essler, P. Fendley, E. Vernier, Strong zero modes in integrable spin-S chains, SciPost Phys. 20, 164 (2026), DOI 10.21468/SciPostPhys.20.6.164, arXiv:2512.07742, https://arxiv.org/abs/2512.07742
  9. S. Gehrmann and F. H. L. Essler, Exact strong zero modes in quantum circuits and spin chains with non-diagonal boundary conditions, SciPost Phys. 20, 142 (2026), DOI 10.21468/SciPostPhys.20.5.142, arXiv:2511.05490, https://arxiv.org/abs/2511.05490
  10. P. R. Pasnoori, Y. Tang, J. Lee, J. H. Pixley, N. Andrei, P. Azaria, Spin fractionalization and zero modes in the spin-1/2 XXZ chain with boundary fields, Phys. Rev. B 112, 075121 (2025), DOI 10.1103/thlq-h58t, arXiv:2312.05970, https://arxiv.org/abs/2312.05970
  11. P. R. Pasnoori, Y. Tang, J. Lee, J. H. Pixley, P. Azaria, N. Andrei, Complete boundary phase diagram of the spin-1/2 XXZ chain with boundary fields in the antiferromagnetic gapped regime, Phys. Rev. B 113, 054424 (2026), DOI 10.1103/5j68-1pqh, arXiv:2507.00386, https://arxiv.org/abs/2507.00386
  12. P. R. Pasnoori, J. Lee, J. H. Pixley, N. Andrei, P. Azaria, Boundary quantum phase transitions in the spin-1/2 Heisenberg chain with boundary magnetic fields, Phys. Rev. B 107, 224412 (2023), DOI 10.1103/PhysRevB.107.224412, arXiv:2212.14832, https://arxiv.org/abs/2212.14832
  13. P. Kattel, P. R. Pasnoori, J. H. Pixley, N. Andrei, Edge modes and boundary impurities in the anisotropic Heisenberg spin chain, Phys. Rev. B 111, 174430 (2025), DOI 10.1103/PhysRevB.111.174430, arXiv:2412.00194, https://arxiv.org/abs/2412.00194
  14. K. K. Kozlowski and V. Terras, Multi-point correlation functions in the boundary XXZ chain at finite temperature, Ann. Henri Poincaré 25, 1007 (2024), DOI 10.1007/s00023-023-01310-4, arXiv:2212.08824, https://arxiv.org/abs/2212.08824
  15. C. Abetian, N. Kitanine, V. Terras, Boundary overlap in the open XXZ spin chain, SciPost Phys. 18, 026 (2025), DOI 10.21468/SciPostPhys.18.1.026, arXiv:2409.15194, https://arxiv.org/abs/2409.15194
  16. C. Hagendorf and J. Liénardy, Open spin chains with dynamic lattice supersymmetry, J. Phys. A: Math. Theor. 50, 185202 (2017), DOI 10.1088/1751-8121/aa67ff, arXiv:1612.02951, https://arxiv.org/abs/1612.02951
  17. R. Weston and J. Yang, Lattice supersymmetry in the open XXZ model: an algebraic Bethe ansatz analysis, J. Stat. Mech. (2017) 123104, DOI 10.1088/1742-5468/aa9f42, arXiv:1709.00442, https://arxiv.org/abs/1709.00442
  18. T. He and X. Wu, Quantum wetting transition in the gapped antiferromagnetic XXZ chain, Phys. Rev. B 112, 064412 (2025), DOI 10.1103/bdh8-vhr9, https://doi.org/10.1103/bdh8-vhr9