Kaleidoscope Yang-Baxter Equation for Gaudin's Kaleidoscope models

TL;DR

The work extends the coordinate Bethe ansatz to Gaudin's kaleidoscope models by introducing the Kaleidoscope Yang-Baxter Equation (KYBE) as a sector-aware integrability condition tied to $D_N$ symmetry. It analyzes two-dimensional LQZC-type systems and Gaudin’s $D_6$ kaleidoscope, showing Bethe ansatz solvability is highly sensitive to symmetry subspaces, with BA equations being overconstrained for $D_{2N}$ ($N\ge4$) while certain sectors remain solvable. A detailed algebraic framework is developed, revealing that KYBE imposes a rich structure linked to the quantum torus algebra and yielding new quantum identities. Numerical validation via finite element methods confirms symmetry-resolved integrability in specific sectors and highlights missing BA states (brick modes), illustrating the nuanced interplay between boundary conditions, symmetry, and solvability in kaleidoscopic quantum systems.

Abstract

Recently, researchers have proposed the Asymmetric Bethe ansatz method - a theoretical tool that extends the scope of Bethe ansatz-solvable models by "breaking" partial mirror symmetry via the introduction of a fully reflecting boundary. Within this framework, the integrability conditions which were originally put forward by Gaudin have been further generalized. In this work, building on Gaudin's generalized kaleidoscope model, we present a detailed investigation of the relationship between DN symmetry and its integrability. We demonstrate that the mathematical essence of integrability in this class of models is characterized by a newly proposed Kaleidoscope Yang-Baxter Equation. Furthermore, we show that the solvability of a model via the coordinate Bethe ansatz depends not only on the consistency relations satisfied by scattering matrices, but also on the model's boundary conditions and the symmetry of the subspace where solutions are sought. Through finite element method based numerical studies, we further confirm that Bethe ansatz integrability arises in a specific symmetry sector. Finally, by analyzing the algebraic structure of the Kaleidoscope Yang-Baxter Equation, we derive a series of novel quantum algebraic identities within the framework of quantum torus algebra.

Figures (15)

• Figure 1: A schematic of the $D_6$ group. Here, ${}_{ r_m}\vec{k}$ and ${}_{s_n}\vec{k}$ denote the points obtained by applying rotation ${ r_m}$ or reflection $s_n$ operations, respectively, to the initial object (for example, the momentum vector $\vec{k}$, which corresponds to the point $r_0 := \mathrm{id}$), $\theta$ denotes the initial angle.
• Figure 3: Schematics of the consistency conditions. The left panel illustrates the equivalence of different paths in the free propagation region, which follows directly from the fact that the matrices $\bm{\beta}_{\vec{x}}$ are diagonal. The middle panel depicts the equivalence between two different paths starting from point A: one path crosses the potential barrier first and then propagates along the barrier, while the other path first propagates along the barrier and then crosses the barrier. This equivalence is also manifested as the commutativity of the matrices $\bm{\beta}_{\vec{x}_k}$ and $\bm{S}_k$. The right panel represents the consistency of the boundary conditions, i.e., a point A that satisfies the boundary condition will continue to satisfy the boundary condition after propagation along the hard wall direction. This property can likewise be established through the commutativity of the matrices $\bm{\beta}_{\vec{x}_k}$ and $\bm{\Gamma}_k$.
• Figure 4: (a) Schematic phase diagram of the Liu-Qi-Zhang-Chen model jackson2024. The original 1D two-body problem is equivalent to a 2D single-particle problem, i.e. a quantum particle confined in a box of size $\frac{1}{2}L \times \frac{\sqrt{3}}{2}L$. (b) Consider the section $A$ of the vector bundle over this region. For computational convenience, we rotated the model clockwise by $90^\circ$.
• Figure 5: Example of numerical solution of the equation (\ref{['eqBAESep']}): for a fixed set of $(n_1, n_2, n_3) = (-2, -1, 1)$, we only need to search for $w=f(x_i)$ such that $x_1 + x_2 + x_3 = 0$ holds, thereby obtaining the unique solution of the BA equations associated with this set of integers.
• Figure 6: The numerical solutions of the LQZC model are obtained using the FEM method. We set the parameters $L = 1$ and $c = 10.0$. The three contour plots in the figure show the ground-state solutions of the three different modes.

Theorems & Definitions (1)

• Proposition 1