Piroli, Pozsgay and Vernier (2017) defined an initial state of an integrable lattice model to be integrable if it is annihilated by all local conserved charges odd under space reflection, and proved that this class contains every state built from boundary K-matrices (solutions of the boundary Yang-Baxter reflection equations) via the rotated-channel construction, including matrix product states of arbitrary bond dimension generated by the action of fused transfer matrices. The problem is to complete this classification: prove, or disprove by an explicit finite-bond-dimension counterexample, that the boundary construction exhausts all integrable matrix product states (equivalently, that the annihilation conditions imply the boundary Yang-Baxter-type equations, currently known only under an unproven irreducibility condition); and classify all finite-dimensional operator-valued solutions of these equations for trigonometric (XXZ) and higher-spin chains. Later literature resolved the nested SU(3) boundary interpretation and the factorized-overlap formulas for rational chains, while the exhaustiveness, trigonometric and higher-spin classification questions remain open.
Background
In a global quantum quench, a many-body system is prepared in a simple, far-from-equilibrium initial state and then evolves unitarily under a Hamiltonian whose dynamics is exactly tractable. The natural laboratory is one-dimensional integrable lattice models such as the XXZ spin chain, whose algebraic Bethe ansatz structure provides a commuting family of transfer matrices (traces, over all lattice sites, of products of the two-particle scattering matrix, the -matrix) generating local conserved charges ; the eigenstates (Bethe states) are parametrized by sets of rapidities (quasi-momenta) . The same framework covers the XXX, XXZ and XYZ chains, called rational, trigonometric and elliptic according to whether the -matrix depends on the rapidity rationally, trigonometrically or through elliptic (doubly periodic) functions. Because of these extensive conserved charges the system does not thermalize, and the exact late-time steady state can be computed by the quench action method, whose essential input is the overlap between the initial state and the Bethe eigenstates. Before 2017, exact quench-action computations were restricted to a handful of initial states for which such overlaps could be guessed and verified.
Piroli, Pozsgay and Vernier introduced a definition unifying these cases: an initial state is integrable if it is annihilated by all local conserved charges that are odd under space reflection (parity reversal of the chain). They proved that this class contains every state related to an integrable boundary condition in a rotated channel, the lattice analogue of the Ghoshal-Zamolodchikov boundary states of integrable quantum field theory, in which the Euclidean time direction is exchanged with space. On the lattice, integrable boundaries are described by -matrices solving the boundary Yang-Baxter (reflection) equation; when the -matrix satisfies a crossing relation, a -matrix produces a two-site block whose repetition along the chain is an integrable initial state. Acting on such boundary states with fused transfer matrices, whose auxiliary space is a fused (higher-spin) representation, generates integrable matrix product states (MPS), i.e. states of the form with finite-dimensional auxiliary matrices (the auxiliary dimension is the bond dimension), of arbitrary bond dimension. For the spin-1/2 XXZ chain every two-site state arises this way and is integrable; for spin only a subset does, because the fused -matrices carry only three free parameters, too few to parametrize the -dimensional space of two-site blocks. In nested higher-rank models such as the -invariant Lai-Sutherland chain, solved by a nested Bethe ansatz with two quasiparticle species, the -matrix lacks the ordinary crossing relation; the boundary construction then passes to twisted (soliton-non-preserving) reflection equations, whose operator-valued solutions are representations of twisted Yangians, which are deformations of the Yangian symmetry algebra attached to symmetric pairs of Lie algebras such as . A direct sequel to the 2017 work showed that for matrix product states with a single repeating block the annihilation conditions are equivalent to a linear intertwining relation for the two-site block, called the square root relation because iterating it twice yields the quadratic reflection equation; the square root relation implies the (twisted) reflection equations whenever the dressed matrix product state, obtained by acting with the fused transfer matrices, is irreducible, meaning that its auxiliary-space matrices admit no common proper invariant subspace. This irreducibility condition has so far only been verified numerically, and the equivalence without it is not established.
In every case where exact overlaps are known they are factorized: they vanish unless the Bethe rapidities form opposite pairs (the pair structure, with exceptions at special non-generic points such as the XX chain at vanishing anisotropy), and are then given by products of one-particle factors and a ratio of Gaudin determinants, the determinants built from the derivatives of the Bethe equations that also give the norms of Bethe states. The 2017 paper left open whether its boundary construction produces all integrable states or only a proper subclass, and how much of this overlap structure follows from the definition itself.
Problem Statement
The problem is to complete the classification of integrable initial states of Bethe-ansatz integrable lattice models, where 'integrable' means annihilated by all local conserved charges odd under space reflection. The common objective is a complete characterization of this class of states; the requested results are its two linked layers, exactly as split by the source paper into a necessity part and a classification part.
-
Exhaustiveness and necessity. Prove, or disprove by an explicit finite-bond-dimension counterexample, that every translationally invariant matrix product state with finite bond dimension annihilated by all odd charges is generated from boundary two-site blocks by the action of (fused) transfer matrices; equivalently, that the annihilation conditions imply the boundary Yang-Baxter-type equations, which in the rational case means the linear square-root intertwiner relation and the twisted reflection equation. The currently known implication assumes an irreducibility condition on the dressed matrix product states that has only been verified numerically; an analytic proof of the implication without the irreducibility assumption, or a delimitation of the class where exhaustiveness fails that exhibits an explicit finite-bond-dimension counterexample, is required.
-
Classification of solutions. Enumerate all finite-dimensional operator-valued solutions of the boundary Yang-Baxter-type equations that produce integrable matrix product states. The rational case is organized by the twisted-Yangian symmetric pairs such as , and (with and the orthogonal and symplectic Lie algebras, and the pairs referring to the symmetric pairs introduced above); a corresponding classification for trigonometric (XXZ) R-matrices, and a complete classification of the integrable two-site blocks for spin chains (where the three fused-K parameters cover only a proper subset of the block space), does not exist.
A complete answer consists of the necessity theorem (or explicit counterexample) together with the classification of all finite-dimensional operator-valued solutions of the boundary Yang-Baxter-type equations; the two items are successive layers of one and the same question, namely what the annihilation definition of an integrable quench actually comprises. The factorized-overlap and pair-structure material described in the Background is a source-discussed consequence of this definition and of later work on it, not a further requested deliverable.
References
- L. Piroli, B. Pozsgay, E. Vernier, What is an integrable quench?, Nucl. Phys. B 925 (2017) 362-402, DOI 10.1016/j.nuclphysb.2017.10.012, arXiv:1709.04796, https://arxiv.org/abs/1709.04796
- B. Pozsgay, L. Piroli, E. Vernier, Integrable Matrix Product States from boundary integrability, SciPost Phys. 6 (2019) 062, DOI 10.21468/SciPostPhys.6.5.062, arXiv:1812.11094, https://arxiv.org/abs/1812.11094
- L. Piroli, E. Vernier, P. Calabrese, B. Pozsgay, Integrable quenches in nested spin chains I: the exact steady states, J. Stat. Mech. (2019) 063103, DOI 10.1088/1742-5468/ab1c51, arXiv:1811.00432, https://arxiv.org/abs/1811.00432
- L. Piroli, E. Vernier, P. Calabrese, B. Pozsgay, Integrable quenches in nested spin chains II: fusion of boundary transfer matrices, J. Stat. Mech. (2019) 063104, DOI 10.1088/1742-5468/ab1c52, arXiv:1812.05330, https://arxiv.org/abs/1812.05330
- T. Gombor, On the classification of rational K-matrices, J. Phys. A 53 (2020) 135203, DOI 10.1088/1751-8121/ab7602, arXiv:1904.03044, https://arxiv.org/abs/1904.03044
- B. Pozsgay, Overlaps with arbitrary two-site states in the XXZ spin chain, J. Stat. Mech. (2018) 053103, DOI 10.1088/1742-5468/aabbe1, arXiv:1801.03838, https://arxiv.org/abs/1801.03838
- Y. Jiang, B. Pozsgay, On exact overlaps in integrable spin chains, JHEP 06 (2020) 022, DOI 10.1007/JHEP06(2020)022, arXiv:2002.12065, https://arxiv.org/abs/2002.12065
- T. Gombor, B. Pozsgay, On factorized overlaps: Algebraic Bethe Ansatz, twists, and Separation of Variables, Nucl. Phys. B 967 (2021) 115390, DOI 10.1016/j.nuclphysb.2021.115390, arXiv:2101.10354, https://arxiv.org/abs/2101.10354
- T. Gombor, On exact overlaps for gl(N) symmetric spin chains, Nucl. Phys. B 983 (2022) 115909, DOI 10.1016/j.nuclphysb.2022.115909, arXiv:2110.07960, https://arxiv.org/abs/2110.07960
- T. Gombor, Exact overlaps for all integrable two-site boundary states of gl(N) symmetric spin chains, JHEP 05 (2024) 194, DOI 10.1007/JHEP05(2024)194, arXiv:2311.04870, https://arxiv.org/abs/2311.04870
- T. Gombor, Exact Overlaps for All Integrable Matrix Product States of Rational Spin Chains, Phys. Rev. Lett. 135 (2025) 150402, DOI 10.1103/4vy2-8cnk, arXiv:2410.23282, https://arxiv.org/abs/2410.23282
- T. Gombor, Derivations for the MPS overlap formulas of rational spin chains, JHEP 10 (2025) 035, DOI 10.1007/JHEP10(2025)035, arXiv:2505.20234, https://arxiv.org/abs/2505.20234
- T. Gombor, C. Kristjansen, V. Moustakis, X. Qian, On exact overlaps of integrable matrix product states: inhomogeneities, twists and dressing formulas, JHEP 02 (2025) 100, DOI 10.1007/JHEP02(2025)100, arXiv:2410.23117, https://arxiv.org/abs/2410.23117
- T. Gombor, C. Kristjansen, Overlaps for matrix product states of arbitrary bond dimension in ABJM theory, Phys. Lett. B 834 (2022) 137428, DOI 10.1016/j.physletb.2022.137428, arXiv:2207.06866, https://arxiv.org/abs/2207.06866
- C. Rylands, B. Bertini, P. Calabrese, Integrable quenches in the Hubbard model, J. Stat. Mech. (2022) 103103, DOI 10.1088/1742-5468/ac98be, arXiv:2206.07985, https://arxiv.org/abs/2206.07985
- X. Qian, X. Zhang, Generalized Integrable Boundary States in XXZ and XYZ Spin Chains, arXiv:2601.16605 (2026), https://arxiv.org/abs/2601.16605
- Y. Liu, J. Wu, Chiral integrable boundary states in the SU(4) alternating spin chain, Sci. China Phys. Mech. Astron. 69 (2026) 231011, DOI 10.1007/s11433-025-2834-3, arXiv:2507.03489, https://arxiv.org/abs/2507.03489
- N. Bai, H. Yang, M.-Z. Shao, Solving for the integrable boundary states of the ABJM spin chain from KT-relations, arXiv:2607.04203 (2026), https://arxiv.org/abs/2607.04203
- Y. Zhang, Y.-H. Wu, M. Cheng, H.-H. Tu, SO(n) Affleck-Kennedy-Lieb-Tasaki states as conformal boundary states of integrable SU(n) spin chains, SciPost Phys. 20 (2026) 067, DOI 10.21468/SciPostPhys.20.3.067, arXiv:2508.13114, https://arxiv.org/abs/2508.13114
- B. Bertini, E. Tartaglia, P. Calabrese, Entanglement and diagonal entropies after a quench with no pair structure, J. Stat. Mech. (2018) 063104, DOI 10.1088/1742-5468/aac73f, arXiv:1802.10589, https://arxiv.org/abs/1802.10589
- M. de Leeuw, C. Kristjansen, S. Mori, AdS/dCFT one-point functions of the SU(3) sector, Phys. Lett. B 763 (2016) 197-202, DOI 10.1016/j.physletb.2016.10.044, https://doi.org/10.1016/j.physletb.2016.10.044