Multiple reference states and complete spectrum of the $Z_n$ Belavin model with open boundaries

The multiple reference state structure of the $\Z_n$ Belavin model with non-diagonal boundary terms is discovered. It is found that there exist $n$ reference states, each of them yields a set of eigenvalues and Bethe Ansatz equations of the transfer matrix. These $n$ sets of eigenvalues together constitute the complete spectrum of the model. In the quasi-classic limit, they give the complete spectrum of the corresponding Gaudin model.

💡 Research Summary

The paper addresses a long‑standing problem in the theory of integrable lattice models: how to obtain the complete set of eigenvalues for the Zₙ Belavin model when the system possesses non‑diagonal open‑boundary conditions. The Belavin model is an elliptic‑type quantum spin chain whose bulk interactions are encoded in an R‑matrix satisfying the Yang‑Baxter equation. Open boundaries are introduced through K‑matrices that solve the reflection equation; when these K‑matrices are non‑diagonal the usual algebraic Bethe ansatz (ABA) fails because the standard ferromagnetic vacuum is no longer an eigenstate of the double‑row transfer matrix.

The authors resolve this difficulty by discovering a “multiple reference state” structure. Exploiting the underlying Zₙ symmetry, they construct n distinct reference vectors |Ω^{(α)}⟩ (α = 1,…,n), each associated with a different weight sector of the quantum space. For each reference state they perform a full ABA analysis: the double‑row monodromy matrix is decomposed into creation operators B_i^{(α)}(u), annihilation operators C_i^{(α)}(u) and diagonal operators A^{(α)}(u), D^{(α)}(u). The commutation relations close thanks to the quasi‑periodicity of elliptic theta functions and the crossing symmetry of the R‑matrix. Consequently, for every α a set of Bethe roots {λ_j^{(α)}} emerges, satisfying a self‑consistent Bethe Ansatz equation (BAE) of the form

\