Completeness and soundness of the off-diagonal Bethe ansatz solution of anisotropic spin-1/2 chains with arbitrary boundary fields: the open XYZ chain

The off-diagonal Bethe ansatz (ODBA) exactly solved the anisotropic spin-1/2 chains with arbitrary boundary fields — open XXZ and XYZ chains with generic non-diagonal boundary magnetic fields, which break the U(1) spin-rotation symmetry required by the conventional Bethe ansatz — by deriving inhomogeneous T-Q relations whose pole-cancellation conditions give Bethe ansatz equations (BAEs). The source work verified completeness, i.e. that the BAEs reproduce all $2^N$ transfer-matrix eigenvalues, only numerically, for the XXZ chain at $N=3,4$. In the subsequent literature the XXZ/XXX (six-vertex) case was proven complete in 2014-2015 by separation-of-variables methods and by a Bezout-counting argument, and the model extensions and thermodynamic-limit program posed by the source work have been substantially carried out. The surviving open core is the eight-vertex case: a rigorous completeness and soundness proof (or refutation) of the ODBA Bethe equations for the open XYZ chain with generic non-diagonal boundary fields, for both parities of the chain length and including the homogeneous limit, together with a criterion separating physical from spurious solutions. Literature of 2025-2026 explicitly states that a complete description of the general open XYZ transfer-matrix spectrum in terms of Q-function solutions of a homogeneous T-Q equation is still unknown, and no completeness or soundness proof of the ODBA Bethe equations for the XYZ chain exists.

Background

Integrable quantum spin chains are lattice models whose dynamics is generated by a one-parameter family of mutually commuting transfer matrices. The transfer matrix is built from an $R$-matrix — a solution of the Yang-Baxter equation that guarantees the factorized consistency of scattering — and, for open (finite) chains, from boundary matrices ($K$-matrices) satisfying Sklyanin's reflection equation, which encodes integrability-compatible boundary conditions (Sklyanin 1988). The spin-1/2 XXZ chain (the anisotropic Heisenberg chain with anisotropy $\Delta=\cos\eta$, where $\eta$ is the crossing parameter) and the fully anisotropic XYZ chain with such boundary fields have Hamiltonians that commute with the transfer matrix, so the transfer-matrix spectrum encodes the full set of conserved quantities. Generic integrable boundary fields are non-diagonal in the $S^z$ basis and depend on six boundary parameters: in the trigonometric (XXZ) case these are the parameters $\alpha_\pm$, $\beta_\pm$, $\theta_\pm$ of the source work, whereas in the elliptic (XYZ) case the source's $K$-matrices (its equations (5.1)-(5.3)) are parametrized by $\alpha_l^{(-)}$ and $\alpha_l^{(+)}$, $l=1,2,3$ (the symbol $\theta_j$ is reserved there for the inhomogeneity parameters).

The central difficulty is that the conventional algebraic Bethe ansatz (ABA) — the standard exact-solution method, which builds eigenstates from a reference state (a 'local vacuum' such as all spins up, annihilated by the off-diagonal entries of the monodromy matrix) — fails for generic non-diagonal boundary fields: they break the U(1) symmetry of conservation of total $S^z$, and no reference state exists. Between Sklyanin's 1988 construction and 2013, exact solutions were known only when the boundary parameters satisfy a constraint (a Nepomechie-type condition relating the six boundary parameters), under which a modified Bethe ansatz works and its completeness was verified numerically (Nepomechie and Ravanini 2003).

The off-diagonal Bethe ansatz (ODBA), introduced by Cao, Yang, Shi and Wang, circumvents the reference state: it uses only functional relations satisfied by every transfer-matrix eigenvalue $\Lambda(u)$ — operator identities at the inhomogeneity points $\theta_j$ (auxiliary site-dependent offsets of the spectral parameter, a standard technical device), of the type $\Lambda(\theta_j)\Lambda(\theta_j-\eta)=\Delta_q(\theta_j)$ where $\Delta_q$ is the quantum determinant (a known scalar function), together with the asymptotics (XXZ) or quasi-periodicity (XYZ) of $\Lambda(u)$. These imply a generalized (inhomogeneous) T-Q relation: a functional equation expressing $\Lambda(u)$ through an auxiliary polynomial $Q(u)$ — Baxter's Q-function, whose zeros are the Bethe roots — plus an extra off-diagonal term proportional to $1/Q(u)$ with no analogue when U(1) symmetry is present. Requiring the right-hand side to be pole-free yields the Bethe ansatz equations (BAEs) for the roots; for open chains the BAEs differ structurally for even and odd chain length $N$. With this method the source work exactly solved the open XXZ and XYZ chains with arbitrary boundary fields (Cao, Yang, Shi and Wang 2013, arXiv:1307.2023), following the isotropic XXX case (Cao, Yang, Shi and Wang 2013, arXiv:1306.1742); the XYZ solution was further elaborated by Cao, Cui, Yang, Shi and Wang (2014, arXiv:1307.0280).

Completeness is the foundational question for any Bethe-ansatz solution: the Hilbert space of $N$ spin-1/2 sites has dimension $2^N$, the scheme is complete if BAE solutions reproduce all $2^N$ transfer-matrix eigenvalues, and sound if every BAE solution is physical (yields a true eigenvalue). For periodic six- and eight-vertex models completeness was proven by Baxter (2002). For open chains with non-diagonal boundary fields the analogous statement long remained unproven; the source work verified it numerically for the XXZ chain at $N=3,4$ only, and left open (i) a rigorous completeness proof for arbitrary $N$, (ii) the classification of admissible T-Q parametrizations, (iii) extensions to higher spin and other algebras, and (iv) the thermodynamic limit.

Substantial progress followed for the six-vertex (XXZ/XXX) case: a Bezout-counting proof that the inhomogeneous T-Q functional relations completely characterize the spectrum for generic non-diagonal boundary matrices, both parities of $N$, the homogeneous chain and the spin-$s$ generalization (Cao, Yang, Shi and Wang 2015); and a separation-of-variables (SoV) program — Sklyanin's method representing the transfer matrix in a basis where it acts diagonally — proving completeness of an inhomogeneous Baxter T-Q equation for generic inhomogeneous open XXZ chains with the most general boundaries (Faldella, Kitanine and Niccoli 2014; Kitanine, Maillet and Niccoli 2014; Maillet and Niccoli 2019). In sharp contrast, the eight-vertex (XYZ) case with generic boundary fields remains without any completeness proof of the ODBA Bethe equations: Niccoli and Terras (2025) solve the open XYZ chain with arbitrary boundary fields by the new SoV approach and, at generic inhomogeneities, characterize the full transfer-matrix spectrum by an inhomogeneous system of $N$ quadratic equations, but a description of the spectrum in terms of Q-function solutions of a T-Q equation is obtained only under a boundary-parameter constraint, and they state explicitly that a complete description of the general open XYZ transfer-matrix spectrum in terms of Q-function solutions of a homogeneous T-Q equation is still unknown; and the recently developed rational Q-system / Wronskian Bethe-equation techniques — a reformulation of the BAEs as polynomial recursion relations whose solutions are automatically physical — exist only for six-vertex-type models, with both their completeness proof and their extension to XYZ explicitly listed as open (Jiang, Liu, Miao and Tan 2026).

Problem Statement

The problem, in the form posed by the source work and surviving the later literature, is the completeness and soundness of the off-diagonal Bethe ansatz solution of the anisotropic spin-1/2 chains with arbitrary boundary fields. The completeness of the open XXZ/XXX (six-vertex) case is settled in substance (see previous progress), although the soundness of its BAEs — the exclusion of 'unphysical' solutions — and the treatment of degenerate spectra remain open even there; the unresolved target of this statement is the eight-vertex case, in two linked parts.

(1) Completeness for the open XYZ chain with generic boundary fields. Prove — or refute by an explicit finite counterexample — that the inhomogeneous T-Q relations (equations (5.33) and (5.39) of Cao, Yang, Shi and Wang, arXiv:1307.2023) together with their Bethe ansatz equations ((5.36)-(5.37) for even $N$, (5.42)-(5.43) for odd $N$), taken with the prescribed number of Bethe roots, account for the full set of $2^N$ transfer-matrix eigenvalues, counted with multiplicity: every eigenvalue is obtained from at least one solution (completeness). Injectivity of the solution-to-eigenvalue map is not part of the criterion: distinct solutions or parametrizations may correspond to the same eigenvalue, as already happens in the six-vertex case, where infinitely many equivalent T-Q parametrizations of a given eigenvalue are known. The statement must cover both parities of the chain length $N$, generic values of the six elliptic boundary parameters $\alpha_l^{(\pm)}$ ($l=1,2,3$) and of the anisotropy (the elliptic modulus or half-period ratio), the homogeneous limit, and must either include or explicitly delimit the degenerate hyperplanes of boundary-parameter space (in particular the constraint loci $\sum_{\sigma=\pm}\sum_{i=1}^{3}\epsilon_i^{\sigma}\alpha_i^{(\sigma)}=(N-2M-1)\eta$ (with signs $\epsilon_i^{\sigma}\in\{\pm1\}$ satisfying $\prod_{\sigma,i}\epsilon_i^{\sigma}=1$, and $M\in\{0,\dots,N\}$ the integer degree parameter of the factorized Q-function), where Q-function solutions of a homogeneous T-Q equation and other solution methods exist).

(2) Soundness and physicality criterion. Prove — or refute by an explicit finite counterexample — that every solution of the inhomogeneous BAEs, taken with the prescribed number of Bethe roots, yields a transfer-matrix eigenvalue (soundness); if spurious solutions exist, give a criterion distinguishing the physical solutions of these equations from them. This component is unresolved even for the six-vertex case, whose completeness proof explicitly allows that not all T-Q relations built from BAE solutions need correspond to right eigenvalues because the BAEs may admit 'unphysical' solutions (Cao, Yang, Shi and Wang 2015, footnote to Proposition 2); in the elliptic case no such criterion exists at all. The relations among the different inhomogeneous T-Q parametrizations used in the source work and its successors (different root numbers and zero-root patterns) remain relevant background for formulating such a criterion, but a full classification theorem for admissible parametrizations is not part of the target: in the six-vertex case all parametrizations of a given eigenvalue are known to be equivalent.

The common scientific objective of the two deliverables is to establish the off-diagonal Bethe ansatz as a complete and sound exact solution of the anisotropic spin-1/2 chains with arbitrary boundary fields — the last unsolved flagship case being the eight-vertex chain with U(1)-breaking boundaries — rather than a formally derived but unverified spectral characterization.

References

1. J. Cao, W.-L. Yang, K. Shi, Y. Wang, 'Off-diagonal Bethe ansatz solutions of the anisotropic spin-1/2 chains with arbitrary boundary fields', Nucl. Phys. B 877 (2013) 152-175. DOI: 10.1016/j.nuclphysb.2013.10.001; arXiv:1307.2023. URL: https://arxiv.org/abs/1307.2023
2. J. Cao, W. Yang, K. Shi, Y. Wang, 'Off-diagonal Bethe ansatz solution of the XXX spin-chain with arbitrary boundary conditions', Nucl. Phys. B 875 (2013) 152-165. DOI: 10.1016/j.nuclphysb.2013.06.022; arXiv:1306.1742. URL: https://arxiv.org/abs/1306.1742
3. J. Cao, S. Cui, W.-L. Yang, K. Shi, Y. Wang, 'Spin-1/2 XYZ model revisit: general solutions via off-diagonal Bethe ansatz', Nucl. Phys. B 886 (2014) 185-201. DOI: 10.1016/j.nuclphysb.2014.06.026; arXiv:1307.0280. URL: https://arxiv.org/abs/1307.0280
4. J. Cao, W.-L. Yang, K. Shi, Y. Wang, 'On the complete-spectrum characterization of quantum integrable spin chains via inhomogeneous T-Q relation', J. Phys. A 48 (2015) 444001. DOI: 10.1088/1751-8113/48/44/444001; arXiv:1409.5303. URL: https://arxiv.org/abs/1409.5303
5. S. Faldella, N. Kitanine, G. Niccoli, 'The complete spectrum and scalar products for the open spin-1/2 XXZ quantum chains with non-diagonal boundary terms', J. Stat. Mech. (2014) P01011. DOI: 10.1088/1742-5468/2014/01/P01011; arXiv:1307.3960. URL: https://arxiv.org/abs/1307.3960
6. N. Kitanine, J.-M. Maillet, G. Niccoli, 'Open spin chains with generic integrable boundaries: Baxter equation and Bethe ansatz completeness from separation of variables', J. Stat. Mech. (2014) P05015. DOI: 10.1088/1742-5468/2014/05/P05015; arXiv:1401.4901. URL: https://arxiv.org/abs/1401.4901
7. J.-M. Maillet, G. Niccoli, 'On separation of variables for reflection algebras', J. Stat. Mech. (2019) 094020. DOI: 10.1088/1742-5468/ab357a; arXiv:1904.00852. URL: https://arxiv.org/abs/1904.00852
8. R. I. Nepomechie, 'An inhomogeneous T-Q equation for the open XXX chain with general boundary terms: completeness and arbitrary spin', J. Phys. A 46 (2013) 442002. DOI: 10.1088/1751-8113/46/44/442002; arXiv:1307.5049. URL: https://arxiv.org/abs/1307.5049
9. Y. Jiang, S. Cui, J. Cao, W.-L. Yang, Y. Wang, 'Completeness and Bethe root distribution of the spin-1/2 Heisenberg chain with arbitrary boundary fields', unpublished preprint (2013). arXiv:1309.6456. URL: https://arxiv.org/abs/1309.6456
10. R. I. Nepomechie, F. Ravanini, 'Completeness of the Bethe Ansatz solution of the open XXZ chain with nondiagonal boundary terms', J. Phys. A 36 (2003) 11391-11401. DOI: 10.1088/0305-4470/36/45/003; arXiv:hep-th/0307095. URL: https://arxiv.org/abs/hep-th/0307095
11. R. J. Baxter, 'Completeness of the Bethe ansatz for the six and eight vertex models', J. Stat. Phys. 108 (2002) 1-48. DOI: 10.1023/a:1015437118218. URL: https://doi.org/10.1023/a:1015437118218
12. E. K. Sklyanin, 'Boundary conditions for integrable quantum systems', J. Phys. A 21 (1988) 2375-2389. DOI: 10.1088/0305-4470/21/10/015. URL: https://doi.org/10.1088/0305-4470/21/10/015
13. G. Niccoli, V. Terras, 'The open XYZ spin 1/2 chain: Separation of Variables and scalar products for boundary fields related by a constraint', SciPost Phys. 19 (2025) 090. DOI: 10.21468/SciPostPhys.19.4.090; arXiv:2402.04112. URL: https://arxiv.org/abs/2402.04112
14. G. Niccoli, V. Terras, 'Correlation functions for open XXZ spin 1/2 quantum chains with unparallel boundary magnetic fields', J. Phys. A 55 (2022) 405203. DOI: 10.1088/1751-8121/ac7ca1; arXiv:2202.12870. URL: https://arxiv.org/abs/2202.12870
15. Y. Jiang, Y.-C. Liu, Y. Miao, Z.-X. Tan, 'Rational Q-systems for integrable spin chains without U(1) symmetry', SciPost Phys. 20 (2026) 175. DOI: 10.21468/SciPostPhys.20.6.175; arXiv:2512.01551. URL: https://arxiv.org/abs/2512.01551
16. Y.-Y. Li, J. Cao, W.-L. Yang, K. Shi, Y. Wang, 'Thermodynamic limit and surface energy of the XXZ spin chain with arbitrary boundary fields', Nucl. Phys. B 884 (2014) 17-27. DOI: 10.1016/j.nuclphysb.2014.04.010; arXiv:1401.3045. URL: https://arxiv.org/abs/1401.3045
17. R. I. Nepomechie, C. Wang, 'Boundary energy of the open XXX chain with a non-diagonal boundary term', J. Phys. A 47 (2014) 032001. DOI: 10.1088/1751-8113/47/3/032001; arXiv:1310.6305. URL: https://arxiv.org/abs/1310.6305
18. B. Pozsgay, O. Rákos, 'Exact boundary free energy of the open XXZ chain with arbitrary boundary conditions', J. Stat. Mech. (2018) 113102. DOI: 10.1088/1742-5468/aae5a5; arXiv:1804.09992. URL: https://arxiv.org/abs/1804.09992
19. P. Sun, Z.-R. Xin, Y. Qiao, K. Hao, L. Cao, J. Cao, T. Yang, W.-L. Yang, 'Surface energy and elementary excitations of the XXZ spin chain with arbitrary boundary fields', J. Phys. A 52 (2019) 265201. DOI: 10.1088/1751-8121/ab2259; arXiv:1901.01514. URL: https://arxiv.org/abs/1901.01514
20. Z. Xin, J. Cao, W.-L. Yang, Y. Wang, 'Surface energy and elementary excitations of the XYZ spin chain with integrable open boundary fields', J. High Energ. Phys. 2024 (2024) 090. DOI: 10.1007/JHEP10(2024)090; arXiv:2404.19240. URL: https://arxiv.org/abs/2404.19240
21. Z. Xin, J. Cao, W.-L. Yang, Y. Wang, 'Exact physical quantities of the XYZ spin chain in the thermodynamic limit', J. High Energ. Phys. 2025 (2025) 086. DOI: 10.1007/JHEP02(2025)086; arXiv:2411.12200. URL: https://arxiv.org/abs/2411.12200
22. X. Zhang, Y.-Y. Li, J. Cao, W.-L. Yang, K. Shi, Y. Wang, 'Bethe states of the XXZ spin-1/2 chain with arbitrary boundary fields', Nucl. Phys. B 893 (2015) 70-88. DOI: 10.1016/j.nuclphysb.2015.01.022; arXiv:1412.6905. URL: https://arxiv.org/abs/1412.6905
23. Y. Wang, W.-L. Yang, J. Cao, K. Shi, 'On the inhomogeneous T-Q relation for quantum integrable models', preprint (2015). arXiv:1506.02512. URL: https://arxiv.org/abs/1506.02512