The generalized non-linear Schrodinger model on the interval

The generalized (1+1)-D non-linear Schrodinger (NLS) theory with particular integrable boundary conditions is considered. More precisely, two distinct types of boundary conditions, known as soliton preserving (SP) and soliton non-preserving (SNP), are implemented into the classical $gl_N$ NLS model. Based on this choice of boundaries the relevant conserved quantities are computed and the corresponding equations of motion are derived. A suitable quantum lattice version of the boundary generalized NLS model is also investigated. The first non-trivial local integral of motion is explicitly computed, and the spectrum and Bethe Ansatz equations are derived for the soliton non-preserving boundary conditions.

💡 Research Summary

The paper investigates the generalized (1+1)-dimensional nonlinear Schrödinger (NLS) model with a $gl_N$ internal symmetry when the system is confined to a finite interval. The authors introduce two distinct integrable boundary conditions—soliton preserving (SP) and soliton non‑preserving (SNP)—by means of reflection matrices $K^{\pm}(\lambda)$ that satisfy the Sklyanin reflection equation. For each type of boundary they construct the double‑row monodromy matrix $U(\lambda)=T(\lambda)K^{-}(\lambda)T^{-1}(-\lambda)K^{+}(\lambda)$, where $T(\lambda)$ is the bulk Lax‑pair transfer matrix. Expanding the logarithm of the trace of $U(\lambda)$ yields an infinite hierarchy of local conserved charges $I_n$. The authors compute explicitly the first few charges. In the SP case the boundary contributions vanish, so the conserved quantities coincide with those of the periodic NLS. In the SNP case additional boundary terms appear in the momentum and energy, reflecting the fact that a soliton reflected from the boundary can change its internal quantum numbers.

The classical equations of motion are derived from the third conserved charge (the Hamiltonian) and are shown to incorporate the chosen boundary conditions in a fully consistent way. The analysis demonstrates that the integrable structure survives the presence of the boundaries and that the two families of boundary conditions lead to qualitatively different dynamics.

To address the quantum problem the authors discretize the model on a lattice, replacing the continuous field by $gl_N$ spin variables at each site. The bulk interaction is encoded in an $R$‑matrix that satisfies the Yang–Baxter equation, while the same $K^{\pm}(\lambda)$ matrices provide integrable reflections at the ends of the chain. Using the algebraic Bethe Ansatz they construct the eigenstates of the double‑row transfer matrix. The first non‑trivial local integral of motion (the quantum Hamiltonian) is obtained explicitly, and the corresponding Bethe Ansatz equations are derived:

\