Closure relations and complete separation-of-variables spectrum characterization for twisted $gl_{\mathcal{M}|\mathcal{N}}$ and Hubbard models

Maillet, Niccoli and Vignoli (SciPost Phys. 9, 060 (2020)) constructed quantum separation-of-variables (SoV) bases for the fundamental inhomogeneous $gl_{\mathcal{M}|\mathcal{N}}$ supersymmetric spin chains and for the inhomogeneous Hubbard model with quasi-periodic twisted boundary conditions given by twist matrices of simple spectrum, proving non-degeneracy and diagonalizability of the transfer matrices in these bases. However, the complete solution of the spectral problem under general twisted boundary conditions was left open: the closure relations (fusion relations plus inner-boundary and null out-boundary conditions) were conjectured, but proven only for the special class of non-invertible simple-spectrum twists in the $gl_{1|\mathcal{N}}$ cases (established in detail for $gl_{1|2}$ ), and the fusion/closure relations for the Hubbard transfer matrices were not derived at all, with both deferred to a future publication that has not appeared. The open problem is to derive the complete fusion/closure relations for the Hubbard transfer matrices and to prove that the resulting finite functional system completely characterizes the transfer-matrix spectrum for both model classes under general simple-spectrum twists, yielding the complete SoV spectral solution (eigenvalues and factorized eigenvectors). The 2025 review by Levkovich-Maslyuk reports that for supersymmetric chains only partial SoV results are known; the eigenvalue-completeness content for fundamental $gl(m|n)$ chains was, however, independently resolved by Chernyak, Leurent and Volin (Comm. Math. Phys. 391 (2022) 969) via Wronskian Bethe equations, by a route that does not prove the SoV closure statements.

Background

A large family of quantum integrable lattice models is built from a solution $R(\lambda)$ of the Yang-Baxter equation, an operator on the tensor product of two auxiliary spaces that encodes the integrable structure. Monodromy matrices are formed as products of $R$ -matrices acting on a chain of $N$ quantum sites, and the transfer matrix — the trace (in the auxiliary space, with grading signs for supersymmetric models) of the monodromy, taken with a twist matrix $K$ inserted before tracing — generates commuting conserved charges. Quasi-periodic (twisted) boundary conditions correspond to a non-trivial $K$ ; a twist matrix has simple spectrum if its eigenvalues are pairwise distinct. For supersymmetric chains based on the Lie superalgebra $gl_{\mathcal{M}|\mathcal{N}}$ the auxiliary space is $\mathbb{C}^{\mathcal{M}+\mathcal{N}}$ with a $\mathbb{Z}_2$ grading, and signs from permuting fermionic indices modify the trace. The Hubbard model — the paradigmatic model of correlated electrons on a lattice — is integrable through Shastry's $R$ -matrix, whose two-particle scattering entangles a pair of spectral parameters, a structure with no direct analogue in the rational or trigonometric $gl_n$ chains.

The quantum separation of variables (SoV) method replaces the $N$ -body spectral problem of the transfer matrix by $N$ one-variable problems: one constructs a basis of the Hilbert space in which a family of commuting operators (the conserved charges contained in the transfer matrix) take prescribed separate values, so that transfer-matrix eigenvectors have factorized wavefunctions and the spectral problem reduces to a finite system of functional equations. Maillet and Niccoli (SciPost Phys. 6, 071 (2019)) realized this program completely for the non-graded rational $Y(gl(n))$ chains, where $Y(gl(n))$ is the Yangian, the algebra of commuting charges generated by the rational monodromy matrices: the SoV basis is generated by multiple actions of transfer matrices on a generic co-vector, and the spectrum is characterized by the fusion relations — functional relations among transfer matrices $T_{\mathcal{V}}(\lambda)$ fused over symmetric and antisymmetric auxiliary representations $\mathcal{V}$ , indexed by Young diagrams — supplemented by boundary conditions on the Young-diagram domain (the inner-boundary condition and the vanishing of transfer matrices outside the domain, the null out-boundary conditions). Maillet, Niccoli and Vignoli (SciPost Phys. 9, 060 (2020)) extended the SoV basis construction to the fundamental inhomogeneous $gl_{\mathcal{M}|\mathcal{N}}$ chains — fundamental meaning that each quantum space carries the defining $(\mathcal{M}+\mathcal{N})$ -dimensional representation, inhomogeneous that the $R$ -matrix at site $a$ is evaluated at a site-dependent spectral shift $\xi_a$ — and to the inhomogeneous Hubbard model with simple-spectrum twists, proving that the transfer matrices then have non-degenerate, simple spectra and are diagonalizable, with factorized SoV wavefunctions. For the spectrum itself they proved the complete characterization only for the $gl_{1|2}$ chain with a special class of non-invertible simple-spectrum twists (zero determinant), where the closure of the fusion relations reduces the spectral problem to a single polynomial functional equation (a quantum spectral curve) and the completeness of the Bethe ansatz — the statement that the eigenstates obtained from the standard ansatz of plane waves parametrized by quasimomenta (the Bethe roots) span the whole Hilbert space — follows as a corollary, and the same special-class characterization extends to the $gl_{1|\mathcal{N}}$ non-invertible twists; for general twists the characterization was stated as Conjecture 3.1, verified on two- and three-site chains, and the complete solution for both model classes was deferred to a future publication. The Hubbard fusion/closure relations were explicitly not derived, the fusion structure being, in the authors' words, rather involved.

Independently, Chernyak, Leurent and Volin (Comm. Math. Phys. 391 (2022) 969) proved completeness of the Wronskian Bethe equations — functional equations among Baxter Q-functions (polynomial eigenvalue functions of mutually commuting Q-operators) written as Wronskian determinants — for rational $gl(m|n)$ spin chains in the defining representation, for generic diagonalizable twists with distinct eigenvalues and for periodic chains, showing the solution count equals the Hilbert-space dimension. This resolves the eigenvalue-completeness question for fundamental supersymmetric chains by a Q-function route, but it does not prove the SoV closure conjecture, does not construct the SoV eigenvectors, and does not address the Hubbard model. For ordinary $gl(n)$ chains the SoV program advanced further (SoV in any compact representation by Ryan and Volin; determinant formulae for correlators by Gromov, Levkovich-Maslyuk and Ryan), while a 2025 review by Levkovich-Maslyuk records that for supersymmetric chains only partial SoV results exist and that closure of the algebra of SoV creation operators is more nontrivial there. The Q-operator construction and determinant formulae for scalar products (overlaps of SoV basis elements and eigenstates) and form factors (matrix elements of local operators between transfer-matrix eigenstates) in the graded and Hubbard SoV frameworks were outlined as expected consequences and remain undone.

Problem Statement

Derive the complete fusion/closure relations for the transfer matrices of the inhomogeneous Hubbard model defined with quasi-periodic twisted boundary conditions given by a general simple-spectrum twist matrix, and prove, for the fundamental inhomogeneous $gl_{\mathcal{M}|\mathcal{N}}$ chains and the inhomogeneous Hubbard model alike, that the resulting finite system of functional and discrete equations — the fusion relations among the fused transfer-matrix eigenvalues together with the inner-boundary condition and the null out-boundary conditions, with the trivial solution excluded in the invertible-twist case — is in bijection with the transfer-matrix spectrum for arbitrary chain length $N$ (Conjecture 3.1 of Maillet, Niccoli and Vignoli and its $gl_{\mathcal{M}|\mathcal{N}}$ generalization). A complete answer supplies, for both model classes and general simple-spectrum twists, the explicit closure system, the proof that every transfer-matrix eigenvalue solves it and that it admits no spurious solutions, and thereby the complete SoV solution of the spectral problem: eigenvalues and the factorized eigenvector wavefunctions in the SoV basis.

References

J.-M. Maillet, G. Niccoli, and L. Vignoli, Separation of variables bases for integrable $gl_{\mathcal{M}|\mathcal{N}}$ and Hubbard models, SciPost Phys. 9, 060 (2020). DOI: 10.21468/SciPostPhys.9.4.060, arXiv:1907.08124. https://arxiv.org/abs/1907.08124
D. Chernyak, S. Leurent, and D. Volin, Completeness of Wronskian Bethe equations for rational $gl(m|n)$ spin chains, Communications in Mathematical Physics 391 (2022) 969-1045. DOI: 10.1007/s00220-021-04275-9, arXiv:2004.02865. https://arxiv.org/abs/2004.02865
F. Levkovich-Maslyuk, Separation of variables for higher rank integrable models: review of recent progress, Journal of Physics A: Mathematical and Theoretical 58 (2025) 243001. DOI: 10.1088/1751-8121/ade103, arXiv:2503.15398. https://arxiv.org/abs/2503.15398
J.-M. Maillet, G. Niccoli, and L. Vignoli, On scalar products in higher rank quantum separation of variables, SciPost Phys. 9, 086 (2020). DOI: 10.21468/SciPostPhys.9.6.086, arXiv:2003.04281. https://arxiv.org/abs/2003.04281
J.-M. Maillet and G. Niccoli, Complete spectrum of quantum integrable lattice models associated to $Y(gl(n))$ by separation of variables, SciPost Phys. 6, 071 (2019). DOI: 10.21468/SciPostPhys.6.6.071, arXiv:1810.11885. https://arxiv.org/abs/1810.11885
P. Ryan and D. Volin, Separation of variables for rational $gl(n)$ spin chains in any compact representation, via fusion, embedding morphism and Bäcklund flow, Communications in Mathematical Physics 383 (2021) 311-343. DOI: 10.1007/s00220-021-03990-7, arXiv:2002.12341. https://arxiv.org/abs/2002.12341
N. Gromov, F. Levkovich-Maslyuk, and P. Ryan, Determinant form of correlators in high rank integrable spin chains via separation of variables, Journal of High Energy Physics 2021 (2021) 169. DOI: 10.1007/JHEP05(2021)169, arXiv:2011.08229. https://arxiv.org/abs/2011.08229
K. Lu and E. Mukhin, Jacobi-Trudi identity and Drinfeld functor for super Yangian, International Mathematics Research Notices 2021 (2021) 16751-16810. DOI: 10.1093/imrn/rnab023, arXiv:2007.15573. https://arxiv.org/abs/2007.15573
K. Lu, On Bethe eigenvectors and higher transfer matrices for supersymmetric spin chains, Journal of High Energy Physics 2023 (2023) 120. DOI: 10.1007/JHEP04(2023)120, arXiv:2209.14416. https://arxiv.org/abs/2209.14416
N. Gromov and F. Levkovich-Maslyuk, New compact construction of eigenstates for supersymmetric spin chains, Journal of High Energy Physics 2018 (2018) 085. DOI: 10.1007/JHEP09(2018)085, arXiv:1805.03927. https://arxiv.org/abs/1805.03927
C. Huang, K. Lu, and E. Mukhin, Solutions of $gl_{m|n}$ XXX Bethe ansatz equation and rational difference operators, Journal of Physics A: Mathematical and Theoretical 52 (2019) 375204. DOI: 10.1088/1751-8121/ab1960, arXiv:1811.11225. https://arxiv.org/abs/1811.11225
L. Vignoli, Quantum separation of variables in higher rank and supersymmetric integrable models, PhD thesis (2021). https://theses.hal.science/tel-03550854 (INSPIRE: https://inspirehep.net/literature/2028004)