---
title: Bethe Eigenlines in Integrable Systems
url: https://www.emergentmind.com/topics/bethe-eigenlines
type: topic
---

# Bethe Eigenlines in Integrable Systems

Searching arXiv for recent and directly relevant papers on Bethe eigenlines and closely related usages.
Bethe eigenlines are one-dimensional common eigenspaces arising in Bethe-ansatz-integrable spectral problems. In the strictest sense represented in the Gaudin-model literature, they are the joint eigenspaces of commutative Gaudin algebras acting on tensor products or multiplicity spaces, and they vary in families over real moduli spaces of marked curves [2507.12829]. In the algebraic Bethe ansatz for spin chains, the same term is naturally attached to the line spanned by an on-shell Bethe vector, since the vector is defined only up to normalization and depends on Bethe-root data rather than on a canonical scale [1304.7978]. In numerical and functional treatments of Bethe–Salpeter equations, the closest analogue is a structured eigenline or paired eigendirection of a parameter-dependent matrix problem, often tracked through eigenvalue curves or preserved by structure-aware eigensolvers [1009.1535], [1501.03830]. Across these settings, the common theme is that Bethe data do not merely produce eigenvalues: they determine distinguished one-dimensional spectral objects whose geometry, regularization, monodromy, and numerical computation are central to the theory.

## 1. Bethe eigenlines in algebraic Bethe ansatz for spin chains

In the isotropic periodic spin-\(\frac12\) Heisenberg chain, the standard ABA begins from the Hamiltonian
\[
H=\frac14 \sum_{n=1}^{N}\left(\vec{\sigma}_n\cdot \vec{\sigma}_{n+1}-1\right), \qquad \vec{\sigma}_{N+1}=\vec{\sigma}_1,
\]
with monodromy matrix
\[
T_a(\lambda)=L_{Na}(\lambda)\cdots L_{1a}(\lambda) = \begin{pmatrix} A(\lambda) & B(\lambda)\\ C(\lambda) & D(\lambda) \end{pmatrix},
\]
transfer matrix
\[
t(\lambda)=A(\lambda)+D(\lambda),
\]
and pseudovacuum
\[
|0\rangle=\binom10^{\otimes N}.
\]
For distinct Bethe roots \(\lambda_1,\dots,\lambda_M\), the standard Bethe vector is
\[
|\lambda_1,\dots,\lambda_M\rangle = B(\lambda_1)\cdots B(\lambda_M)|0\rangle,
\]
with roots constrained by the Bethe equations
\[
\left(\frac{\lambda_k+\frac{i}{2}}{\lambda_k-\frac{i}{2}}\right)^N = \prod_{\substack{j=1\\j\neq k}^{M} \frac{\lambda_k-\lambda_j+i}{\lambda_k-\lambda_j-i}, \qquad k=1,\dots,M.
\]
The transfer-matrix action has the standard wanted term plus unwanted terms, and the Bethe equations are exactly the conditions that remove all unwanted contributions [1304.7978].

For generic nonsingular roots, this gives the usual correspondence: a solution of the Bethe equations determines a Bethe vector and therefore a common eigenline of the commuting transfer matrices and of the Hamiltonian [1304.7978]. The paper on singular solutions makes explicit that the eigenline, not the unnormalized vector, is the invariant object. This becomes decisive when the vector formula degenerates.

A closely related situation appears for the open XXX chain with general boundaries. There the transfer matrix contains both \(\mathscr B(u)\) and \(\mathscr C(u)\), so \(U(1)\) symmetry is broken and the standard fixed-magnetization ABA no longer closes. The generalized construction replaces the usual Bethe vector by
\[
|\Phi^N(\lambda_1,\lambda_2,\dots,\lambda_N)\rangle=
\overline{\mathscr{B}(\lambda_1)}\, \overline{\mathscr{B}(\lambda_2)}\cdots \overline{\mathscr{B}(\lambda_N)}|\Omega\rangle,
\]
where the modified creation operator is
\[
\overline{\mathscr{B}(\lambda)}=
\mathscr{B}(\lambda)+\frac{\rho}{\xi^-}\left(\frac{2\lambda}{2\lambda+1} \mathscr{A}(\lambda)-\mathscr{D}(\lambda)\right)-\left(\frac{\rho}{\xi^-}\right)^2\mathscr{C}(\lambda),
\]
with
\[
\rho=1-\sqrt{1+\xi^+\xi^-}.
\]
When the corresponding Bethe equations hold, the off-shell transfer-matrix action reduces to a wanted term, so the vector becomes an eigenvector. Since it is meaningful only up to an overall factor, the associated object is again an eigenline [1309.6165].

This suggests a broad ABA principle: Bethe equations parametrize eigenlines, while explicit Bethe vectors are representatives of those lines whose normalization is secondary. A plausible implication is that “Bethe eigenline” is often the conceptually correct formulation whenever Bethe vectors are only canonically defined projectively.

## 2. Singular solutions and the failure of naive root-to-eigenline correspondence

The clearest breakdown of the naive Bethe-solution/eigenvector dictionary occurs for singular solutions of the XXX chain containing
\[
\lambda_1=\frac{i}{2}, \qquad \lambda_2=-\frac{i}{2}.
\]
For \(N>4\), this pair is an exact formal solution of the Bethe equations, but it is singular in two ways. First, the energy formula
\[
E=-\frac12\sum_{k=1}^{M}\frac{1}{\lambda_k^2+\frac14}
\]
diverges because \(\lambda^2+\frac14=0\) at \(\lambda=\pm \frac{i}{2}\). Second, the Bethe vector becomes ill-defined: in the conventions of the paper, \(B(\tfrac{i}{2})\) is finite while \(B(-\tfrac{i}{2})\) is singular, producing components of indeterminate \(0/0\) type [1304.7978].

The central conceptual point is that a formal solution of the Bethe equations containing \(\pm \frac{i}{2}\) does not automatically define a genuine transfer-matrix eigenvector, or equivalently a genuine Bethe eigenline [1304.7978]. Root data alone are insufficient at the singular locus.

A naive regularization,
\[
\lambda_1^{\text{naive}}=\frac{i}{2}+\epsilon, \qquad \lambda_2^{\text{naive}}=-\frac{i}{2}+\epsilon,
\]
recovers the correct energy limit \(E=-1\) for the two-root state, but not the correct eigenvector. The reason is that the unwanted coefficients in the off-shell action have the wrong asymptotics:
\[
F_1\sim \epsilon^N, \qquad F_2\sim 1,
\]
whereas one needs
\[
F_1(\lambda,\{\lambda\})\sim \epsilon^{N+1}, \qquad F_2(\lambda,\{\lambda\})\sim \epsilon.
\]
The cure is the higher-order regularization
\[
\lambda_1=\frac{i}{2}+\epsilon+c\,\epsilon^N, \qquad \lambda_2=-\frac{i}{2}+\epsilon,
\]
with \(c\) fixed by generalized consistency conditions [1304.7978].

For even \(N\), both asymptotic conditions are satisfied by
\[
c=2i(-1)^{N/2}.
\]
For odd \(N\), they cannot be satisfied simultaneously, so the singular two-string is unphysical [1304.7978]. In the general singular configuration
\[
\left\{\frac{i}{2},-\frac{i}{2},\lambda_3,\dots,\lambda_M\right\},
\]
compatibility reduces to the practical criterion
\[
\left[ -\prod_{k=3}^{M} \frac{\lambda_k+\frac{i}{2}}{\lambda_k-\frac{i}{2}} \right]^N=1.
\]
A singular solution is physical iff this condition holds, equivalently iff one can choose \(c\) so that the regularized Bethe vector becomes a genuine eigenvector [1304.7978].

This criterion sharply refines completeness questions. For \(N=6\), \(M=3\), there are \(5\) singular solutions of the formal Bethe equations, but only \(1\) is physical. For \(N=8\), \(M=4\), there are \(21\) singular solutions, but only \(3\) are physical [1304.7978]. The practical conclusion is that the space of admissible Bethe eigenlines is strictly smaller than the space of formal Bethe-root solutions once singular roots are allowed.

## 3. Geometric and categorical Bethe eigenlines in Gaudin models

In the Gaudin-model setting, Bethe eigenlines are treated directly as geometric objects. For dominant weights \(\lambda_1,\dots,\lambda_n\), write
\[
V(\underline\lambda):=V(\lambda_1)\otimes\cdots\otimes V(\lambda_n).
\]
For regular \(\chi\in\mathfrak h\), the inhomogeneous Gaudin algebra
\[
\mathcal A_\chi(C)\subset U(\mathfrak g)^{\otimes n}
\]
acts on \(V(\underline\lambda)\), preserving weight spaces, and the paper denotes by
\[
\mathcal E_\chi(C,\underline\lambda)
\]
the set of eigenlines of \(\mathcal A_\chi(C)\) on \(V(\underline\lambda)\) [2507.12829]. In the homogeneous case over \(\overline M_{n+1}\), one similarly considers eigenlines of \(\mathcal A(C)\) on multiplicity spaces
\[
\operatorname{Hom}_{\mathfrak g}(V(\mu),V(\underline\lambda)).
\]

These eigenlines vary over real cactus flower moduli spaces. The families
\[
\mathcal E_n=\bigsqcup_{\underline\lambda}\mathcal E_\chi(\underline\lambda),\qquad
\mathcal X_{n+1}=\bigsqcup_{\underline\lambda,\mu}\mathcal X(\underline\lambda,\mu)
\]
form a \(\Lambda_+\)-coloured operadic covering of \(\overline F_n(\mathbb R)\) [2507.12829]. The key structural reason is factorization of Gaudin algebras on boundary strata, for example
\[
\mathcal A_\chi(\alpha_k(C_1,C_2)) = \mathcal A_\chi(C_1)\otimes \mathcal A_\chi(C_2),
\]
together with analogous formulas for mixed and homogeneous boundary gluings [2507.12829].

The resulting conceptual statement is stronger than a mere parametrization of eigenvectors. Isomorphism classes of operadic coverings of \(\overline F_n(\mathbb R)\) are naturally one-to-one with equivalence classes of concrete coboundary monoidal categories satisfying certain semisimplicity and finiteness conditions, and the covering defined by Bethe eigenlines reconstructs the concrete coboundary category of Kashiwara \(\mathfrak g\)-crystals [2507.12829]. The paper proves an equivalence
\[
\Psi:\mathcal C \xrightarrow{\sim} \mathfrak g\text{-}\mathtt{Crys},
\]
sending simple objects \(L(\lambda)\) to crystals \(B(\lambda)\) [2507.12829].

This yields a monodromy interpretation. Over \(\overline F_n(\mathbb R)\), Bethe eigenlines carry an action of the virtual cactus group \(vC_n\), and there are bijections
\[
\mathcal E_\chi(\infty,\underline\lambda)\cong \mathcal B(\lambda_1)\times\cdots\times \mathcal B(\lambda_n)
\]
compatible with the \(vC_n\)-actions [2507.12829]. Over split and compact real loci of the trigonometric deformation, the monodromy is described by the mirabolic cactus group \(MC_n\) and the extended affine cactus group \(\widetilde{AC}_n\), acting respectively by elementary transpositions or cyclic rotations together with the usual cactus-group action on tensor-product crystals [2507.12829].

In this geometric regime, Bethe eigenlines are literally points in topological covering spaces. Their importance lies not only in encoding spectra of Gaudin Hamiltonians, but in carrying the operadic and monodromic data from which crystal tensor combinatorics can be recovered. This suggests that “Bethe eigenline” is, in this context, a genuinely geometric notion rather than a shorthand for an eigenspace.

## 4. Bethe eigenlines, Jordan degeneration, and generalized eigenvectors

A different use of the eigenline perspective appears in non-diagonalizable limits of integrable spin chains. In the eclectic spin chain, obtained as a strong-twist limit of a finitely twisted \(su(3)\) XXX chain,
\[
\hat H_{\xi_1,\xi_2,\xi_3}=\lim_{\epsilon\to 0}\,\epsilon\, H\!\left(\frac{\xi_1}{\epsilon},\frac{\xi_2}{\epsilon},\frac{\xi_3}{\epsilon}\right),
\]
the limiting Hamiltonian becomes defective and develops Jordan blocks [2112.13883].

For finite twist, one has ordinary Bethe eigenvectors from nested coordinate Bethe ansatz. In the strong-twist limit, several of these eigenvectors collapse onto the same limiting “locked state”
\[
|(k)\rangle=\sum_{l=1}^L e^{2\pi i mkl/L}\,U^l\,|1,\ldots,1,2,\ldots,2,3,\ldots,3\rangle,
\]
so distinct finite-\(\epsilon\) Bethe eigenlines coalesce [2112.13883]. The missing spectral information is not lost; it is encoded in suitable singular linear combinations of the colliding Bethe vectors, which produce generalized eigenvectors and full Jordan chains.

The general mechanism is formulated for a family \(M(\epsilon)\) of diagonalizable matrices approaching a defective limit \(M(0)\). If eigenvectors \(v_i(\epsilon)\) coalesce, then limits of singular combinations such as
\[
\lim_{\epsilon\to0}\frac{\alpha_1v_1+\alpha_2v_2}{\lambda_1-\lambda_2}
\]
recover rank-2 generalized eigenvectors, and higher-rank constructions follow recursively [2112.13883]. The paper proves that, when the limiting Jordan blocks are distinguishable, these limits reconstruct the full generalized eigenspace.

This is relevant to Bethe eigenlines because it shows that a family of ordinary Bethe eigenlines can degenerate into non-semisimple spectral data. The limit is not simply a set of surviving eigenlines; it carries extra extension data detectable only through the singular asymptotics of the Bethe vectors. A plausible implication is that the projective viewpoint is especially natural near exceptional points: the collision of eigenlines is the primary event, and Jordan structure emerges from how these lines approach one another.

## 5. Bethe–Salpeter eigenlines as structured eigendirections

In Bethe–Salpeter theory, the phrase “Bethe eigenline” is not standard, but several papers treat closely related objects: one-dimensional eigendirections, eigenvalue branches, and paired invariant subspaces of the Bethe–Salpeter Hamiltonian.

After discretization, the homogeneous Bethe–Salpeter equation takes the matrix form
\[
\lambda(P^2)\,\vec F_{[h]}(P^2)=\mathbf K(P^2)\cdot \vec F_{[h]}(P^2),
\]
and physical masses satisfy
\[
\lambda_n(P^2)=1 \quad \text{at} \quad P^2=-M_n^2.
\]
Each \(\lambda_n(P^2)\) is an eigenvalue trajectory, and its crossing with \(1\) marks a bound state [1009.1535]. In this sense, a “Bethe eigenline” is the \(P^2\)-dependent eigendirection \(\vec F_n(P^2)\) associated with an eigenvalue branch.

The inhomogeneous equation,
\[
(\mathbf 1-\mathbf K)\vec F=\vec F_0,\qquad \vec F=(\mathbf 1-\mathbf K)^{-1}\vec F_0,
\]
detects the same spectrum through poles of the resolvent [1009.1535]. Poles appear precisely when some eigenline satisfies \(\lambda_i(P^2)=1\). Thus eigenline crossings and poles are two views of the same spectral event [1009.1535].

For the finite-dimensional Bethe–Salpeter Hamiltonian
\[
H=\begin{bmatrix} A & B \\ -\overline{B} & -\overline{A} \end{bmatrix},\qquad A=A^H,\quad B=B^T,
\]
the spectrum has the symmetry
\[
\lambda,\ -\lambda,\ \overline{\lambda},\ -\overline{\lambda},
\]
and under the definiteness condition
\[
\begin{bmatrix} A & B \\ \overline{B} & \overline{A} \end{bmatrix}\succ 0
\]
all eigenvalues are real and occur in \(\pm\lambda\) pairs [1501.03830]. The corresponding right and left eigenvector matrices can be written
\[
X=\begin{bmatrix} X_1 & \overline{X_2}\\ X_2 & \overline{X_1} \end{bmatrix},\qquad
Y=\begin{bmatrix} X_1^H & -X_2^H\\ -X_2^T & X_1^T \end{bmatrix},
\]
with
\[
X_1^H X_1-X_2^H X_2=I_n.
\]
Thus the positive eigendirections are not isolated objects; each determines companion negative and dual directions through the BSE symmetry [1501.03830].

Later work recasts this geometry in explicit indefinite-inner-product language. For a definite Bethe–Salpeter Hamiltonian
\[
H=C_n\Omega,\qquad
C_n=\begin{bmatrix}I_n&0\\0&-I_n\end{bmatrix},
\]
one has a structured spectral decomposition
\[
H= \begin{bmatrix} X & \overline{Y}\\ Y & \overline{X} \end{bmatrix}
\begin{bmatrix} \Lambda & 0\\ 0 & -\Lambda \end{bmatrix}
\begin{bmatrix} X & -\overline{Y}\\ -Y & \overline{X} \end{bmatrix}^{\!\herm},
\]
with
\[
\Phi(X,-Y)^\herm \Phi(X,Y)=I_{2n},
\qquad
\Phi(U,V)=\begin{bmatrix} U & \overline{V}\\ V & \overline{U} \end{bmatrix}.
\]
If
\[
z_i=\begin{bmatrix}x_i\\y_i\end{bmatrix}
\]
is an eigenvector for \(+\lambda_i\), then
\[
\widetilde z_i=\begin{bmatrix}\overline{y_i}\\ \overline{x_i}\end{bmatrix}
\]
is the corresponding eigenvector for \(-\lambda_i\) [2603.07910]. The relevant eigenspaces lie in
\[
\mathcal C(2n,2k)\cap \mathcal H(2n,2k),
\]
meaning they admit bases that are simultaneously of Bethe–Salpeter block form and \(C_n\)-orthonormal [2603.07910].

The algorithmic literature is built around preserving this structure. Structure-preserving dense and parallel algorithms reduce the BSE to a real Hamiltonian problem and then to a skew-symmetric eigenproblem [1501.03830], [1912.04062]. Structure-preserving \(\Gamma\)QR and \(\Gamma\)-Lanczos methods keep the \(\Pi^{-}\)-Hermitian form intact, so computed eigenvalues and eigenvectors retain the symmetry pairing of the original BSE matrix [1801.00884]. A recent structure-preserving LOBPCG solver computes a few smallest positive eigenpairs while maintaining the \(C_n\)-orthogonal paired geometry [2603.07910].

The precise terminology in these papers is “eigenvectors,” “paired eigendirections,” “invariant subspaces,” and “eigenvalue trajectories,” not “Bethe eigenlines.” Still, the notion is effectively present: physically relevant spectral data are organized into one-dimensional eigendirections or their paired counterparts, and preserving that organization is the central numerical concern.

## 6. Reconstruction, orthogonality, and computational access to Bethe eigenlines

Several papers address how Bethe-type spectral lines are reconstructed when standard ABA machinery is unavailable or when explicit orthogonality/completeness must be established by other means.

For integrable models solved by off-diagonal Bethe Ansatz, the inhomogeneous \(T\)-\(Q\) relation gives eigenvalues but not automatically eigenvectors. The retrieval method constructs an orthogonal basis of the Hilbert space from commuting monodromy entries and determines an eigenstate by its overlaps with this basis. In the XXZ spin torus model, the candidate Bethe state is
\[
|\lambda_1,\ldots,\lambda_N\rangle =
\prod_{j=1}^{N}\left\{\frac{D(\lambda_j)}{d(\lambda_j)}\right\} |\Omega;\{\theta_j\}\rangle,
\]
with a reconstructed reference state
\[
|\Omega;\{\theta_j\}\rangle =\sum_{l=0}^{N}\frac{(B^-)^l}{[l]_q!}|0\rangle.
\]
For the open XXX chain, the reconstructed Bethe states are
\[
{}_B\langle \lambda_1,\ldots,\lambda_N| =\langle 0|\prod_{j=1}^{N}\bar C(\lambda_j), \qquad
|\lambda_1,\ldots,\lambda_N\rangle_B=\prod_{j=1}^{N}\bar B(\lambda_j)|0\rangle.
\]
In both cases the state is determined uniquely up to an overall scalar factor by its overlaps with a complete orthogonal basis, so what is reconstructed is naturally an eigenline [1407.5294].

Orthogonality of Bethe eigenfunctions is the main result in a different continuum setting: the Laplacian on the hyperoctahedral Weyl alcove with Robin boundary conditions. The Bethe Ansatz eigenfunctions
\[
v(\xi_\lambda,x),\qquad \lambda\in\Lambda
\]
form an orthogonal basis of \(L^2(A,dx)\) for \(g,g_+,g_->0\) [1602.02152]. These functions arise as the continuum limit of an orthogonal basis of algebraic Bethe Ansatz eigenfunctions for a finite \(q\)-boson system. In the finite model, the \(n\)-particle Bethe states
\[
\Psi(v_1,\dots,v_n):= \hat{\mathcal B}_m(v_1;a_-)\cdots \hat{\mathcal B}_m(v_n;a_-)\,|0\rangle
\]
diagonalize the commuting transfer operator and Hamiltonian, and distinct spectral points give orthogonal states [1602.02152]. This is precisely the setting in which each admissible Bethe solution determines a distinguished spectral line.

Recent computational work accesses Bethe-root-parameterized eigenstates directly. In the XXZ chain, Bethe states are used as the variational ansatz in VQE, with Bethe roots treated as variational parameters. The coordinate Bethe state
\[
|B^L_M (\vec{k})\rangle= \sum_{w\in P(L,M)}f(w, \vec{k})\,|w\rangle
\]
is optimized by minimizing either energy or variance, thereby estimating Bethe roots corresponding to actual eigenstates [2404.18244]. The paper explicitly notes that this approach is not limited to real Bethe roots and that, when optimization converges to an eigenstate, the resulting roots are automatically physical [2404.18244]. It does not formulate a geometric eigenline theory, but it gives a practical parameter-space route to projective Bethe states.

A different computational direction appears in the 2026 study of three periodic spin-\(\tfrac12\) chains solved by coordinate and nested Bethe Ansatz. For the scalar model \(Y1\), the eigenstate is parametrized by momenta \(\{p_j\}\); for the nested model \(Y2\), by \(\{k_j\}\) and auxiliary roots \(\{\lambda_\alpha\}\); for the model \(Y3\), by first-level roots \(\{z_k\}\) together with an eigenline of a non-\(U(1)\)-symmetric free-fermionic auxiliary transfer matrix [2603.29932]. The last case is especially notable because the auxiliary problem is free-fermionic and lacks \(U(1)\)-invariance, so the spectral line at the nesting level is not encoded by a conventional auxiliary-magnon number sector [2603.29932].

These examples show that the construction of Bethe eigenlines is not confined to standard ABA. They can be reconstructed from functional relations, proved orthogonal through finite-volume limits, or searched variationally through root-parameterized trial states.

## 7. Conceptual synthesis and common misconceptions

A recurring misconception is that a solution of Bethe equations automatically determines a physical eigenstate. The singular-root analysis of the XXX chain is the clearest counterexample: formal Bethe solutions may fail to define finite, nonzero transfer-matrix eigenvectors, and additional generalized consistency conditions are required [1304.7978]. Likewise, in non-diagonalizable limits, ordinary Bethe eigenvectors may collide, and their limits alone do not exhaust the spectral data; generalized eigenvectors must be extracted from singular combinations [2112.13883].

Another misconception is that Bethe eigenlines are always isolated one-dimensional objects with no extra structure. In Gaudin models they form coverings with operadic gluing and nontrivial monodromy, recovering coboundary tensor-category data [2507.12829]. In Bethe–Salpeter problems they typically appear in paired \(\pm\)-symmetric families and are naturally embedded in structured invariant subspaces rather than treated as independent lines [1501.03830], [2603.07910], [1801.00884].

A third misconception is that the term applies uniformly across all “Bethe” contexts. The data suggest a more careful picture. In the strict Gaudin and ABA literature, “Bethe eigenline” refers quite directly to the line spanned by an on-shell Bethe vector or to the common eigenline of a commutative Bethe algebra [1304.7978], [2507.12829]. In Bethe–Salpeter literature, the same wording is best understood as an interpretive bridge rather than established terminology: the actual objects are structured eigendirections, eigenvalue branches, and paired invariant subspaces [1009.1535], [1501.03830].

Taken together, the literature supports a robust definition. A Bethe eigenline is the one-dimensional spectral object singled out by Bethe data—roots, auxiliary roots, or more general commuting-algebra eigenconditions—once normalization is quotiented out. Depending on context, that line may require regularization, may sit inside a monodromic family over moduli space, may degenerate into Jordan structure, or may come paired with a symmetry-related companion. The unifying point is that Bethe methods do not merely compute eigenvalues: they organize spectral theory around these line-level objects.

Source: https://www.emergentmind.com/topics/bethe-eigenlines